Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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 199 records · Page 11

Operator inference for non-intrusive model reduction with quadratic manifolds

This paper proposes a novel approach for learning a data-driven quadratic manifold from high-dimensional data, then employing this quadratic manifold to derive efficient physics-based reduced-order models. The key ingredient of the approach is a polynomial mapping between high-dimensional states and a low-dimensional embedding. This mapping consists of two parts: a representation in a linear subspace (computed in this work using the proper orthogonal decomposition) and a quadratic component. The approach can be viewed as a form of data-driven closure modeling, since the quadratic component introduces directions into the approximation that lie in the orthogonal complement of the linear subspace, but without introducing any additional degrees of freedom to the low-dimensional representation. Combining the quadratic manifold approximation with the operator inference method for projection-based model reduction leads to a scalable non-intrusive approach for learning reduced-order models of dynamical systems. Further, applying the new approach to transport-dominated systems of partial differential equations illustrates the gains in efficiency that can be achieved over approximation in a linear subspace.

42 ENGINEERING↗

Dynamic substructure analysis of the International Ultraviolet Explorer (IUE) spacecraft

The results are presented of a steady state vibration analysis of the IUE spacecraft simulating its sinusoidal vibration test. The model of the spacecraft, including solar arrays and the scientific instrument, consisted of three separate substructure models one of which was used to represent the two identical solar arrays. Substructuring techniques were used since the large overall size of the problem precluded solving it utilizing a single model. The models used for each substructure are discussed, including reduction to an acceptable size for the combined dynamic analysis. The DMAP alters needed for performing the modal analysis and the subsequent modal frequency response and substructure data recovery are included. Comparison of the results with data obtained during vibration tests of the spacecraft are also included.

Case, W. R.↗

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Langley 14- by 22-foot subsonic tunnel test engineer's data acquisition and reduction manual

The Langley 14- by 22-Foot Subsonic Tunnel is used to test a large variety of aircraft and nonaircraft models. To support these investigations, a data acquisition system has been developed that has both static and dynamic capabilities. The static data acquisition and reduction system is described; the hardware and software of this system are explained. The theory and equations used to reduce the data obtained in the wind tunnel are presented; the computer code is not included.

Quinto, P. Frank↗

Simulation of aircraft crash and its validation

A joint FAA/NASA program is discussed which is aimed at developing a reliable technology for the design of crashworthy light aircraft. This program encompasses the development of analytical methods, the definition of a survivable crash envelope, and the design of improved seat and restraint systems. A facility for full-scale crash-simulation testing is described along with the test method and results of five full-scale crash tests of twin-engine light aircraft. The major goals of the analytical portion of the program are outlined, including the development and validation of the analytical technique using simplified structural specimens that approximate aircraft components, as well as the mathematical modeling of the complete airframe and its subsequent dynamic analysis by substructuring and matrix reduction techniques.

Alfaro-Bou, E.↗

Estimability And Regulability Of Linear Systems

Report presents two new properties of systems characterized by linear state space models (e.g., dynamical systems and associated control systems): estimability and regulability. Provides criteria for reductions of errors and cost functions in control systems. Useful in design of optimal controllers and estimators.

Baram, Yoram↗

A family of dynamic models for large-eddy simulation

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Carati, D.↗

Enforcing constraints for time series prediction in supervised, unsupervised and reinforcement learning

We assume that we are given a time series of data from a dynamical system and our task is to learn the flow map of the dynamical system. We present a collection of results on how to enforce constraints coming from the dynamical system in order to accelerate the training of deep neural networks to represent the flow map of the system as well as increase their predictive ability. In particular, we provide ways to enforce constraints during training for all three major modes of learning, namely supervised, unsupervised and reinforcement learning. In general, the dynamic constraints need to include terms which are analogous to memory terms in model reduction formalisms. Such memory terms act as a restoring force which corrects the errors committed by the learned flow map during prediction. For supervised learning, the constraints are added to the objective function. For the case of unsupervised learning, in particular generative adversarial networks, the constraints are introduced by augmenting the input of the discriminator. Finally, for the case of reinforcement learning and in particular actor-critic methods, the constraints are added to the reward function. In addition, for the reinforcement learning case, we present a novel approach based on homotopy of the action-value function in order to stabilize and accelerate training. We use numerical results for the Lorenz system to illustrate the various constructions.

Stinis, Panagiotis↗

Fractal Model for Drag Reduction on Multiscale Nonwetting Rough Surfaces

Rough surfaces in contact with a flow of fluid exhibit alternating no-slip and free shear boundary conditions at the solid–liquid and air–liquid interfaces, respectively, thereby potentially offering drag reduction benefits. The balance between the dynamic pressure in the flow and the restoring capillary pressure in the interasperity spaces determines the stability of the Cassie state of wettability and is a function of the relative extent of no-slip and free shear regions per unit surface area. In the present study, using a fractal representation of rough surface topography, an analytical model is developed to quantify the stability of the Cassie state of wettability as well as drag reduction and the friction factor for laminar flow in a rectangular channel between nonwetting multiscale rough surfaces. A systematic study is conducted to quantify the effects of fractal parameters of the surfaces and the flow Reynolds number on drag reduction and the friction factor. The studies are used to develop friction factor curves extending the classical Moody diagram to hydrophobic and superhydrophobic surfaces. On the basis of the studies, regime maps are derived for estimating the extent of drag reduction offered by hydrophobic and superhydrophobic surfaces, revealing that superhydrophobic surfaces do not always offer the best drag reduction performance. The application of the fractal model to practical topographies of nonwetting surfaces of copper, aluminum, and zinc oxide fabricated via electrodeposition and etching is also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics-Informed Active Learning With Simultaneous Weak-Form Latent Space Dynamics Identification

The parametric greedy latent space dynamics identification (gLaSDI) framework has demonstrated promising potential for accurate and efficient modeling of high-dimensional nonlinear physical systems. However, it remains challenging to handle noisy data. Here, to enhance robustness against noise, we incorporate the weak-form estimation of nonlinear dynamics (WENDy) into gLaSDI. In the proposed weak-form gLaSDI (WgLaSDI) framework, an autoencoder and WENDy are trained simultaneously to discover intrinsic nonlinear latent-space dynamics of high-dimensional data. Compared with the standard sparse identification of nonlinear dynamics (SINDy) employed in gLaSDI, WENDy enables variance reduction and robust latent space discovery, therefore leading to more accurate and efficient reduced-order modeling. Furthermore, the greedy physics-informed active learning in WgLaSDI enables adaptive sampling of optimal training data on the fly for enhanced modeling accuracy. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including viscous and inviscid Burgers' equations, time-dependent radial advection, and the Vlasov equation for plasma physics. With data that contains 5%–10% Gaussian white noise, WgLaSDI outperforms gLaSDI by orders of magnitude, achieving 1%–7% relative errors. Compared with the high-fidelity models, WgLaSDI achieves 121 to 1779x speed-up.

97 MATHEMATICS AND COMPUTING↗

Redox poise in R. rubrum phototrophic growth drives large-scale changes in macromolecular pathways

During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H 2 ), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H 2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.

59 BASIC BIOLOGICAL SCIENCES↗

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗

A cavitation and dynamic void growth model for a general class of strain-softening amorphous materials

In this work, we derive a homogenized framework for cavitation and dynamic void growth in a general class of strain-softening materials with particular emphasis on amorphous materials. Analytic solutions are derived by considering an assemblage of spherical shells subject to dynamic hydrostatic tensile pressure. The isochoric response of the matrix material is captured by a free volume theory based amorphous viscoplasticity model. The framework accounts for the coupled rate-dependent, strain-softening behavior typically exhibited in amorphous materials. Validation of the homogenized theory is carried out against numerous direct numerical simulations with excellent agreement. Furthermore, we find that a reduction in cavitation strength and significant acceleration of dynamic void growth rates is induced as a consequence of strain-softening in the amorphous matrix material. Lastly, the model is utilized to interpret some non-intuitive experimental observations of amorphous materials subject to shock compression and subsequent spall failure.

36 MATERIALS SCIENCE↗

Control and dynamics of a flexible spacecraft during stationkeeping maneuvers

A case study of a spacecraft having flexible solar arrays is presented. A stationkeeping attitude control mode using both earth and rate gyro reference signals and a flexible vehicle dynamics modeling and implementation is discussed. The control system is designed to achieve both pointing accuracy and structural mode stability during stationkeeping maneuvers. Reduction of structural mode interactions over the entire mode duration is presented. The control mode using a discrete time observer structure is described to show the convergence of the spacecraft attitude transients during Delta-V thrusting maneuvers without preloading thrusting bias to the onboard control processor. The simulation performance using the three axis, body stabilized nonlinear dynamics is provided. The details of a five body dynamics model are discussed. The spacecraft is modeled as a central rigid body having cantilevered flexible antennas, a pair of flexible articulated solar arrays, and to gimballed momentum wheels. The vehicle is free to undergo unrestricted rotations and translations relative to inertial space. A direct implementation of the equations of motion is compared to an indirect implementation that uses a symbolic manipulation software to generate rigid body equations.

Liu, D.↗

The Adjoint Petrov–Galerkin method for non-linear model reduction

Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.

42 ENGINEERING↗