Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model reduction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Model Reduction for Control System Design

An approach and a technique for effectively obtaining reduced order mathematical models of a given large order model for the purposes of synthesis, analysis and implementation of control systems is developed. This approach involves the use of an error criterion which is the H-infinity norm of a frequency weighted error between the full and reduced order models. The weightings are chosen to take into account the purpose for which the reduced order model is intended. A previously unknown error bound in the H-infinity norm for reduced order models obtained from internally balanced realizations was obtained. This motivated further development of the balancing technique to include the frequency dependent weightings. This resulted in the frequency weighted balanced realization and a new model reduction technique. Two approaches to designing reduced order controllers were developed. The first involves reducing the order of a high order controller with an appropriate weighting. The second involves linear quadratic Gaussian synthesis based on a reduced order model obtained with an appropriate weighting.

Enns, D. F.↗

An eigensystem realization algorithm for modal parameter identification and model reduction

A method called the eigensystem realization algorithm is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular-value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system and noise modes. For illustration of the algorithm, an example is shown using experimental data from the Galileo spacecraft.

Juang, J.-N.↗

An Eigensystem Realization Algorithm (ERA) for modal parameter identification and model reduction

A method, called the Eigensystem Realization Algorithm (ERA), is developed for modal parameter identification and model reduction of dynamic systems from test data. A new approach is introduced in conjunction with the singular value decomposition technique to derive the basic formulation of minimum order realization which is an extended version of the Ho-Kalman algorithm. The basic formulation is then transformed into modal space for modal parameter identification. Two accuracy indicators are developed to quantitatively identify the system modes and noise modes. For illustration of the algorithm, examples are shown using simulation data and experimental data for a rectangular grid structure.

Juang, J. N.↗

Dynamic Model Reduction for Large-Scale Power Systems Using Wide-Area Measurements

To perform faster than the real-time dynamic simulation of large-scale power systems, it is necessary to reduce the simulated system size by using equivalents for surrounding areas of the study area, and existing dynamic model reduction approach could provide the needed structure of the reduced area. However, further parameter optimization is required to achieve the desired accuracy. In this paper, a particle swarm optimization (PSO) based approach is used to solve the above problem. Parameters for the individual dynamic elements in the reduced system are calibrated repeatedly until the wide-area measurements of the reduced model and the original model are very similar to each other with satisfactory accuracy. Results indicate that after optimization, the dynamic response of the reduced model matches better with that of the original one than using existing methods. Under both the generator-trip event and the bus-fault event, the reduced model has a higher frequency match and less power mismatch.

42 ENGINEERING↗

Model reduction and control of flexible structures using Krylov subspaces

Krylov vectors and the concept of parameter-matching are combined to develop a model reduction algorithm for a damped structural dynamics system. The reduced-order model obtained matches a certain number of low-frequency moments of the full-order system. The major application of the present method is to the control of flexible structures. It is shown that, in the control of flexible structures, there generally exist three types of control energy spillover, namely, the control spillover, the observation spillover, and dynamic spillover. The formulation based on Krylov subspaces can eliminate the control and the observation spillover, while leaving only the dynamic spillover to be considered. Two examples are used to illustrate the efficacy of the Krylov method.

Craig, Roy R., Jr.↗

A consistent model reduction of measured modal parameters for reduced-order active control

The problem of synthesizing reduced-order linear models of vibrating structures for the design of fixed-order dynamic feedback control is investigated. The present technique builds on a recently developed procedure for constructing an objective set of mass and stiffness matrices from measured modal parameters that are akin to the Craig-Bampton synthesized ones obtained from finite element models. The constructed mass and stiffness matrices are determined directly from the identification of experimental data, however, rather than through correlation or reconciliation of a finite element model. A model truncation criterion is then applied to the identified minimum-order mass and stiffness model to satisfy certain observability/controllability requirements for the reduced model. Numerical examples illustrate the effectiveness of the proposed technique for synthesizing reduced-order controllers from system realizations of experimental data. The dynamic performance of the resulting closed-loop models is assessed using the known full-order structural dynamics and compared with existing model reduction techniques.

Alvin, K. F.↗

A COMPACT model reduction methodology for articulated, multi-flexible bodies structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one need appropriate reduced-order models for all the flexible components involved. The component modes projection and assembly model reduction methodology is one way to construct these reduced-order component models. However, the generated reduced-order system models typically contain high frequency extraneous modes. The presence of high frequency extraneous modes in a reduced-order system model forces us to use a small stepsize in the numerical integration of that model, resulting in an increase in the elapsed time of the integration process. To eliminate these high-frequency modes, we first identify and discard the component modes that 'contributed' to the generations of these modes. A sensitivity analysis is then used to determine how the remaining component modes should be perturbed in order to still closely capture the selected system modes. The effectiveness of the proposed methodology has been successfully verified using a finite-element model of the Galileo spacecraft.

Lee, Allan Y.↗

Robot arm dynamic model reduction for control

General methods are described by which the mathematical complexities of explicit and exact state equations of robot arms can be reduced to a simplified and compact state equation representation without introducing significant errors into the robot arm dynamic model. The model reduction methods are based on homogeneous coordinates and on the Langrangian algorithm for robot arm dynamics, and utilize matrix, vector and numeric analysis techniques. The derivation of differential vector representation of centripetal and Coriolis forces which has not yet been established in the literature is presented.

Bejczy, A. K.↗

Model reduction in the physical coordinate system

In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.

Yae, K. Harold↗

Model reduction methods for population dynamics with fast-switching environments: Reduced master equations, stochastic differential equations, and applications

Here we study stochastic population dynamics coupled to fast external environments and combine expansions in the inverse switching rate of the environment and a Kramers–Moyal expansion in the inverse size of the population. This leads to a series of approximation schemes, capturing both intrinsic and environmental noise. These methods provide a means of efficient simulation and we show how they can be used to obtain analytical results for the fluctuations of population dynamics in switching environments. We place the approximations in relation to existing work on piecewise-deterministic and piecewise-diffusive Markov processes. Finally, we demonstrate the accuracy and efficiency of these model-reduction methods in different research fields, including systems in biology and a model of crack propagation.

59 BASIC BIOLOGICAL SCIENCES↗

Unsymmetric Lanczos model reduction and linear state function observer for flexible structures

This report summarizes part of the research work accomplished during the second year of a two-year grant. The research, entitled 'Application of Lanczos Vectors to Control Design of Flexible Structures' concerns various ways to use Lanczos vectors and Krylov vectors to obtain reduced-order mathematical models for use in the dynamic response analyses and in control design studies. This report presents a one-sided, unsymmetric block Lanczos algorithm for model reduction of structural dynamics systems with unsymmetric damping matrix, and a control design procedure based on the theory of linear state function observers to design low-order controllers for flexible structures.

Su, Tzu-Jeng↗

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗

Thermochemical reduction modeling in a high-temperature moving-bed reactor for energy storage: 1D model

The design of robust and efficient high-temperature thermochemical reactors and determination of operating conditions are critical steps toward enabling high-efficiency long-duration solar energy storage. This work presents a computational model for the thermal reduction of a metal oxide material (Mg-Mn-O) up to 1450 °C and the coupled complex transport phenomena in a novel tubular thermal reactor design that features the capability for a high extent-of-reduction (high energy storage density) and inherent heat recuperation. A one-dimensional model coupling counter-current gas–solid flow, two-phase heat transfer, thermochemical redox reactions, and species transport in a moving-bed reactor is developed. Simplified versions of the model are validated with published results in the literature for packed beds with both inert and reactive particles; the fully coupled model is also validated with experimental measurements of a moving-bed reactor in terms of local temperatures and oxygen release at the exit. Detailed comparisons on the effects of different boundary conditions in the reaction zone (prescribed wall temperature vs. heat flux conditions) and formulations based on a simple uniform flow assumption vs. plug flow using Ergun equation for gas flow are investigated. The results are compared with experimental measurements, and for all cases, the energy flow components in the reactor system and the thermal to chemical conversion efficiency and overall system efficiency are computed. Finally, the predicted high thermal-to-chemical efficiency ~95% and system efficiency ~30% agree with experimental measurements.

25 ENERGY STORAGE↗

Adaptive Model Reduction for Parareal in Time Method for Transient Stability Simulations

Real time or faster than real time simulation can enable system operators to foresee the effect of crucial contingencies on the power system dynamics and take timely actions to prevent system instability. Parareal in time method uses concurrent computations on different segments of the time domain of interest to speed up the dynamic simulations. This paper describes the application of an adaptive nonlinear model reduction method in improving computational speed of the Parareal solver. The proposed method adaptively switches between a hybrid system with selective linearization and a completely linear system based on the size of a disturbance. The functions in the hybrid system are linearized based on the electrical distance between specific generators and the area where disturbances originated. The proposed method is tested on the 327-machine 2383-bus Polish system.

Osipov, Denis↗

Statistical Learning for Nonlinear Model Reduction from Local Simulations of Stochastic and Particle- and Agent-Based Systems

Stochastic physical systems across the sciences that have very high-dimensional state spaces, with a large number of fast degrees of freedom that force direct simulators to proceed by integration steps that are orders of magnitude smaller than events of interests (e.g., particle collisions). Examples range from molecular motion to dynamics of large populations of cells. A grand challenge in the simulation and understanding of such systems is the systematic construction of accurate, interpretable, reduced models, enabling faster simulations, revealing fundamental properties of the dynamics, and predicting phenomena of interest that the original simulator could not reached with sufficient accuracy or within a given computational budget. In this projected we developed novel statistical estimation/machine learning techniques for analyzing and building empirical reduced models for important families of high-dimensional stochastic systems, in particular: - we developed techniques for estimating interaction kernels in interacting particle- and agent-based systems, which are ubiquitous in Physics, Biology and many other sciences, given observed trajectories of the system; - we developed techniques for nonlinear model reduction for high-dimensional stochastic systems that have a small number of unknown, nonlinear slow variables, and a large number of fast modes, that are possibly of large magnitude, given observed short trajectories of the system in the form of bursts of trajectories from different initial conditions; - we developed novel techniques for estimating linear dynamical systems on graphs when both the dynamics and the underlying graph are unknown, and we have a sparse set of space-time observations; - we considered the problem of estimating an unknown nonlinear observation function of a standard process (e.g. Brownian motion), so that we can recognized if an observed dynamics is "just" a nonlinear version of a known dynamics; we also developed benchmarks for learning algorithms aimed at learning and classifying diffusion processes.

97 MATHEMATICS AND COMPUTING↗

Hybrid CMS methods with model reduction for assembly of structures

Future on-orbit structures will be designed and built in several stages, each with specific control requirements. Therefore there must be a methodology which can predict the dynamic characteristics of the assembled structure, based on the dynamic characteristics of the subassemblies and their interfaces. The methodology developed by CSC to address this issue is Hybrid Component Mode Synthesis (HCMS). HCMS distinguishes itself from standard component mode synthesis algorithms in the following features: (1) it does not require the subcomponents to have displacement compatible models, which makes it ideal for analyzing the deployment of heterogeneous flexible multibody systems, (2) it incorporates a second-level model reduction scheme at the interface, which makes it much faster than other algorithms and therefore suitable for control purposes, and (3) it does answer specific questions such as 'how does the global fundamental frequency vary if I change the physical parameters of substructure k by a specified amount?'. Because it is based on an energy principle rather than displacement compatibility, this methodology can also help the designer to define an assembly process. Current and future efforts are devoted to applying the HCMS method to design and analyze docking and berthing procedures in orbital construction.

Farhat, Charbel↗

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗