Sample selection based on sensitivity analysis in parameterized model order reduction

Many applications concerning physical and technical processes employ dynamical systems for simulation purposes. The increasing demand for a more accurate and detailed description of realistic phenomena leads to high dimensional dynamical systems and hence, simulation often yields an increased computational effort. An approximation, e.g. with model order reduction techniques, of these large-scale systems becomes therefore crucial for a cost efficient simulation. This paper focuses on a model order reduction method for linear time in-variant (LTI) systems based on modal approximation via dominant poles. There, the original large-scale LTI system is projected onto the left and right eigenspaces corresponding to a specific subset of the eigenvalues of the system matrices, namely the dominant poles of the system's transfer function. Since these dominant poles can lie anywhere in the spectrum, specialized eigenvalue algorithms that can compute eigentriplets of large and sparse matrices are required. The Jacobi-Davidson method has proven to be a suitable and competitive candidate for the solution of various eigenvalue problems and hence, we discuss how it can be incorporated into this modal truncation approach. Generalizations of the reduction technique and the application of the algorithms to second-order systems are also investigated. The computed reduced order models obtained with this modal approximation can be combined with the ones constructed with Krylov subspace or balanced truncation based model order reduction methods to get even higher accuracies.

Microgrids (MGs) concepts are a systematic way to integrate the supply and demand side along with the storage devices. These mostly use renewable energy resources (RESs) viz. solar PV and wind turbines, interfaced to the low voltage distribution network via an inverter. As several components are involved, MGs has a system of systems (SoSs) architecture and its modeling leads to form a high order system. Hence, for the controllers' design purposes, these models had to be brought down to a reduced lower order system. This paper performs a comparative study of various model order reduction (MOR) techniques viz. aggregation, balanced truncation (BT), singular perturbation (SP) and Hankel norm approximation (HNA) on a two-source autonomous MG system. The standard performance evaluation methods used for a multi-input multi-output (MIMO) system i.e. the step and impulse responses, eigenvalue analysis and Bode plots were considered to determine the effectiveness of these MOR method.

Many new promising model order reduction (MOR) methods and algorithms were developed during the last decade. Industry and academic research institutions intend to test, validate, compare, and use these new promising MOR techniques with their own models. Therefore, an MOR toolbox bridging the gap between theoretical, algorithmic, and numerical developments to an end-user-oriented program, usable by non-experts, was developed called ‘Model Order Reduction of Elastic Multibody Systems’ (Morembs). A C++ implementation as well as a Matlab implementation including an intuitive graphical user interface is available. Import from various FE programs is possible, and the reduced elastic bodies can be exported to a variety of programs to simulate the compact models. In the course of the various projects, many improvements on the algorithmic side were added. As we learned over the years, there is not one ‘optimal’ MOR method. ‘Optimal’ MOR depends on circumstances, like boundary conditions, excitation spectra, further model usage. The toolbox is now used, e.g., in solid mechanics, biomechanics, vehicle dynamics, control of flexible structures, or crash simulations. In all these use cases, the toolbox allows the user to facilitate their well-known modeling and simulation environment. Only the critical MOR process during preprocessing is performed with Morembs, which helps to compare the various MOR techniques to find the most suited one.

During lamination, residual thermo-mechanical stresses are induced in the encapsulated solar cells composing photovoltaic (PV) modules. Depending on the material and geometrical configuration of the layers of the laminate, this residual stress field can be beneficial since it may lead to a compressive stress state in Silicon and therefore crack closure effects in the presence of cracks, with a recovery of electrical conductivity in cracked solar cells. It is therefore important to investigate the distribution of thermo-mechanical stresses within the PV laminate with a view to optimizing the coupling between the electrical response and elastic deformation in the operation of PV modules. A promising approach proposed in the present thesis regards the prediction of residual stresses in composite laminates by using a shear-lag theory to model the epoxy-vinil-acetate polymeric layers, accounting for their thermo-visco-elastic response. Moreover, it will be shown that thermomechanical formulations for stress analysis of a PV laminate lead to a system of higher order ordinary differential equations or partial differential equations in which the exact solutions may be impossible to be determined in closed form and hence numerical schemes become desirable. However, the computational cost associated with the implementation of the numerical scheme may be significantly expensive. Therefore, a method to reduce the computational complexity is expected to be very important. To this aim, Model Order Reduction (MOR) techniques are applied hierarchically, first to the thermal system of a PV module in service, and then extended to coupled thermo-mechanical problems. A combination of proper orthogonal decomposition (POD) and discrete empirical interpolation method (DEIM) with a modified formulation is proposed for the first-order thermal equations of photovoltaic system during service and a new coupled second-order Krylov based formulation is developed for model order reduction of the coupled thermo-mechanical model of the photovoltaic module. The results of these reduction schemes show a huge computational gain in the reduced system solutions and a high accuracy of the reduced system outputs

Numerical analysis of wide band large antenna array plays a crucial role in the prediction of design performance prior to construction. Since these arrays mostly are electrically large and include complex geometries, accurate simulation is a great challenge. However, recent developments of domain decomposition method (DDM) [1][2] provides us an efficient way to obtain field solution at each frequency point. However, for many real-life phased array applications, the spectral responses over a wide frequency band are usually needed. For frequency domain CEM approaches, the model order reduction (MORe) techniques are often employed for fast frequency sweep. In this paper, we study a few MORe techniques and apply them for large finite antenna array. Through our numerical studies and plots of pole distributions, we found that the number of expansion points increases significantly with the size of array. By using the current existing MORe method, the costs for modeling electrically large electromagnetic problems are still prohibitive.

The reduction of the large order models of complex systems into a sufficiently low dimensional order in such a way that important dynamic characteristics of the real system are preserved, is a nontrivial task. The motivations for such a model order reduction are either to reduce computations for analysis and practical control design or to simplify the control system structure. In the present article, major techniques in model order reduction are briefly discussed. The characteristic similarities and differences are highlighted, their advantages and disadvantages are underlined, and further research requirements and trends are commented upon. Introduction and Perspective One of the central issues in the active control of complex systems such as large flexible space structures (LFSS) is the derivation of a correct mathematical model for both the controlled and the uncontrolled dynamical systems. Theoretically, there are infinitely many elastic modes or degrees of freedom (DOF) in the distributed parameter (DP) models of LFSS, usually with very low natural damping. Moreover, the flexible modes contribute to the actual deformation of the structure.' The truly infinite dimensional character of LFSS models has to be approximated by some fidelity finite (but usually very large) dimensional model. The normal approach taken by engineers to achieve this end is via modal models with a large number of modes that provide a reasonable representation of the spacecraft dynamic characteristics. However, a difficult problem still remains to be the development of a model of lowenough dimensional order that it can be utilized by the onboard controller, yet high enough dimensional order that it preserves the dynamic characteristics of the real system represented and controlled. The motivations for such a reduction are either to reduce computations for analysis and practical control design or to simplify the control system structure. Discretization procedures will not be discussed herein and it will be assumed that a large finite dimensional model has been generated somehow and for implementation and other practical considerations, the model needs to be further reduced. *This work was-supported, in part, by the Air Force Wright Aeronautical, Air Force Systems Command, Flight Dynamics Lab, Structures and Dynamics Division, while the author was employed with HR-Textron, Valencia, California. There are various approaches to model order reduction; some are or pseudo-optimal and others are ad-hoc methods based on practical considerations and engineering experience and judgment. These techniques are known by different names such as reduction, condensation, economization, aggregation, optimal projection, and other combinations thereof. It has often been pointed out that such techniques generally constitute application of Rayleigh-Ritz/Galerkin optimization, and matrix transformation methods to the eigenvalue/eigenvector formulation for structural dynamic problems. The reduction or the condensation methods are based on transformations of the coordinates in the equations of motion that essentially maintain the invariance of the quadratic forms of the potential and kinetic energies. An important feature of these methods is that the reduced order model often loses the basic characteristics of the original system. Considerable progress was made by Likins, Ohkami, and wong5 in this respect. However, they failed to develop a general enough criterion that could reduce the system model in an sense without significantly affecting the eigenvalues of the original model. The optimization or the mathematical reduction procedure, on the other hand, is based on the reduction of the eigenvalue problem to a smaller size based on some optimality criterion (usually quadratic). In these techniques, the reductions are carried out without actually resorting to approximations, i .e. , without truncatiag any coordinates or states of the original system. Balanced model reduction of linear timeinvariant dynamical systems is essentially based on the controllability and observability relations of the states of the system. Subsystem models are obtained by deleting those states that contribute the least to the controllability and observability (or the impulse response) of the original system7 and thus are only in this sense. A more recent approach to model reduction is proposed in Skelton's work,* where each state of the system model is assigned a relative to a given basis, via a quadratic criterion, and the states with the least cost are deleted in a systematic manner. The resulting reduced model is a function of the state-space basis, and thus there is no guarantee for optimality for all choices. The latest development in model order reduction techniques is the work by Hyland and erns stein.^ Herein first order necessary conditions for reduced order modeling of linear time-invariant systems are derived via a pair of modified Lyapunov equations coupled by a nonorthogonal projection. This approach reveals the possibility of multiple extrema forsome of the abovementioned methods. Copyright @ 1988 by the American Institute of Aeronautics and Astronautics. Inc. All Rights Reserved.

Numerical simulations for flow and transport in subsurface porous media often prove computationally prohibitive due to property data availability at multiple spatial scales that can vary by orders of magnitude. A number of model order reduction approaches are available in the existing literature that alleviate this issue by approximating the solution at a coarse scale preserving fine scale features. We attempt to present a comparison between two such model order reduction techniques, namely (1) adaptive numerical homogenization and (2) generalized multiscale basis functions. We rely upon a non-linear, multi-phase, black-oil model formulation, commonly encountered in the oil and gas industry, as the basis for comparing the aforementioned two approaches. An expanded mixed finite element formulation is used to separate the spatial scales between non-linear, flow, and transport problems. A numerical benchmark is setup using fine scale property information from the 10th SPE comparative project dataset for the purpose of comparing accuracies of these two schemes. An adaptive criterion is employed by both the schemes for local enrichment that allows us to preserve solution accuracy compared to the fine scale benchmark problem. The numerical results indicate that both schemes are able to adequately capture the fine scale features of the model problem at hand.

In this work, a model order reduction (MOR) technique for heat-transfer system-level modeling is presented. A heat transfer solver using the finite difference method (FDM), which is appropriate for MEMS thermal analysis, is implemented. The numerical models generated by the FDM solver then were successfully reduced into low-order compact models by an Arnoldi-based model-order-reduction technique. The compact macromodels can be easily inserted into system-level or circuit simulators. We have demonstrated that the macromodel results for a thermal actuator are in good agreement with the experimental results. We also use the same technique on periodic surface heat flux problems to observe the relationship between frequency and minimum required order of the reduced model under specific error tolerance.

A new mixed method for relative error model order reduction is proposed. In the proposed method the frequency domain balanced stochastic truncation method is improved by applying the generalized singular perturbation method to the frequency domain balanced system in the reduction procedure. The frequency domain balanced stochastic truncation method, which was proposed in [15] and [17] by the author, is based on two recently developed methods, namely frequency domain balanced truncation within a desired frequency bound and inner-outer factorization techniques. The proposed method in this paper is a carry over of the frequency-domain balanced stochastic truncation and is of interest for practical model order reduction because in this context it shows to keep the accuracy of the approximation as high as possible without sacrificing the computational efficiency and important system properties. It is shown that some important properties of the frequency domain stochastic balanced reduction technique are extended to the proposed reduction method by using the concept and properties of the reciprocal systems. Numerical results show the accuracy, simplicity and flexibility enhancement of the method.

<div class="section abstract"><div class="htmlview paragraph">Battery temperature is known to have a critical influence on overall battery pack performance, from electrochemical behavior, charge acceptance, power availability, trip efficiency, safety, reliability and life-cycle costs. Temperature monitoring is critical to ensure safe and reliable battery pack operation.</div><div class="htmlview paragraph">Monitoring of cell temperatures in battery packs allows for control and estimation algorithms that can ensure homogenous pack temperature distribution, prevent excessive pack temperature rise and even infer cell core temperature, potentially allowing to both predict and mitigate onset of thermal runaway.</div><div class="htmlview paragraph">The increasing need for improved accuracy requires inclusion of more detail in the modelling stage, leading inevitably to ever larger-scale, ever more complex dynamical systems. Simulations in such large-scale settings lead in turn to unmanageably large demands on computational resources, which is the main motivation for Model Order Reduction.</div><div class="htmlview paragraph">Model Order Reduction is focused on reducing the complexity of large-scale dynamical systems, while preserving their input-output behavior. Reduced models mimic the complex behavior of large-scale dynamical systems, and can be efficiently used for design, optimization and sensitivity analysis. The resulting reduced order model can be used to replace the original system as a component in a larger simulation or it might be used to develop a simpler and faster controller suitable for real time applications.</div><div class="htmlview paragraph">Three model order reduction techniques are compared. The first method, Balanced Truncation, assumes access to the high-fidelity model structure is available whilst the other two, Nyquist and neural network approximation, assume only input/output behavioral knowledge is available. The reduced order models are compared, in simulation, to a high-fidelity model of a battery pack in terms of simulation speed-up, model accuracy and property preservation. The applicability of each technique for real-time applications, particularly in the context of xEV (grid connected electric vehicle) thermal management, together with pros and cons of each method is discussed.</div></div>

This paper presents three types of linear model order reduction (MOR) technique, namely singular value decomposition (SVD)-based, Krylov-based, and modal truncationbased type applied to large-scale wind farm models. The first type includes a Balanced Truncation (BT) and Alternating Direction Implicit (ADI)-based BT method, while the second type encompasses a Rational Krylov (RK), and Iterative Rational Krylov Algorithm (IRKA) method. In the third type, a Subspace Accelerated MIMO Dominant Pole Algorithm (SAMDP) method is used. The effectiveness of these methods are tested on practicalsized wind farms with 90, 120 and 210 doubly-fed induction generators (DFIGs). Merits and demerits of each method are discussed in detail. The reduced order model (ROM) of wind farm is validated against the full order model (FOM) in term of frequency domain indices and waveform agreement at the point of common coupling (PCC).

Sensitivity analysis of random linear differential–algebraic equations using system norms

In steam turbines (STs), evaluation technology of thermal deformations during transient operations is increasingly important, because of a demand for improvement of operability such as shortening startup time. However, it is still difficult to predict temperatures of metal parts with sufficient accuracy due to complexity of heat transfer phenomena. The present study puts its focus on developing a method to identify metal temperatures of STs during start-up operations using particle filter (PF), a kind of sequential Bayesian filter utilizing ensemble approximation, and model order reduction (MOR) technique. In this method, a self-organizing state-space model is used to calculate time evaluations of metal temperatures, and heat transfer coefficient (HTC) parameters and steam bulk temperature parameters. And an optimal estimate of these is obtained through iterations of a short-time prediction and a correction of variables using measured temperatures. In this state-space model, state variables including HTC and bulk temperature parameters are treated as random variables, and the time evolution of each variable is modeled as follows; (1) Metal temperature is modeled using a reduced order model (ROM) in order to reduce computational time of PF. The ROM is constructed from a finite element model for unsteady heat transfer analysis using MOR technique. (2) HTC is modeled as a random walk model where its changes during one time step is randomly determined using Gaussian noise. Here, the magnitude of this noise is adjusted automatically. (3) Bulk temperature is modeled using prediction formulas with Gaussian noise added. As a validation problem, a cold start-up operation of a ST unit of a gas turbine combined cycle (GTCC) is considered. The proposed method is applied to the problem, and boundary conditions (BCs) are identified using measured temperatures at 68 measurement points. Then, the heat transfer analysis based on finite element analysis (FEA) is performed using the identified BCs. As a result of this FEA, it is confirmed that the calculated metal temperature tends to agree better with measurements compared with that of initial FEA, and that the errors of the calculated temperature at measurement points reduce by 41% on average compared with initial FEA. From this result, it is concluded that this proposed method is effective for metal temperature identification during start-up operations.

In the present work, multiple non existing model order reduction (MOR) techniques for unstable second‐order form systems (SOSs) are proposed. For unstable SOSs, continuous‐time algebraic Lyapunov equations get unsolvable that halt the reduction process. To avoid this problem, unstable SOS is first decomposed into stable and unstable portions and balanced truncation is applied to the stable part. The obtained reduced order model (ROM) for the stable portion is augmented with the unstable portion to obtain the overall reduced system. It is observed that the second‐order structure in ROM for the first technique gets lost as well as augmented unstable dynamics degrade the ROM performance. To cater to these constraints, two structure‐preserving second‐order balanced truncation techniques for unstable SOSs are proposed in second part. System is first stabilized using the Bernoulli feedback stabilization procedure and then gramians are computed for the stabilized system. Further, gramians are partitioned into position and velocity portions to achieve structure preservation in ROM and balanced truncation is applied. As second technique retains second‐order structure as well as involves stabilized system dynamics, this technique far closely approximates original system behavior. Proposed techniques are tested on multiple systems and results certify the correct development of proposed techniques that can be utilized for MOR applications of unstable SOSs. © 2021 Institute of Electrical Engineers of Japan. Published by Wiley Periodicals LLC.

Transient Electromagnetic problems constitute an area of significant investigative effort. The principal computational issue in these problems is the solution of large system of differential algebraic equations (DAE's) obtained after Finite Element (FE) discretization. Model Order Reduction (MOR) Techniques provide a mechanism to generate reduced order models from the detailed description of the original FE network. This is achieved by using moment matching techniques, where the reduced order model matches the moments of the original system to approximate the response with a low order transfer function. However, these numerical techniques all conserve the original system moments only with respect to frequency. While this provides a significant CPU cost advantage when performing a single frequency sweep, a new reduced order model is required each time a parameter is varied in the structure under study. This necessitates the use of parametric MOR strategies so as to expedite optimization and design space exploration cycles. In this work, we present a methodology for transient electromagnetic field simulations through parameterized model order reduction (pMOR). The proposed methodology is illustrated for a generic system with promising results and a significant saving in computational effort.