Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order 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 217 records · Page 12

An application of the IMC software to controller design for the JPL LSCL Experiment Facility

A software package which Integrates Model reduction and Controller design (The IMC software) is applied to design controllers for the JPL Large Spacecraft Control Laboratory Experiment Facility. Modal Cost Analysis is used for the model reduction, and various Output Covariance Constraints are guaranteed by the controller design. The main motivation is to find the controller with the 'best' performance with respect to output variances. Indeed it is shown that by iterating on the reduced order design model, the controller designed does have better performance than that obtained with the first model reduction.

Zhu, Guoming↗

Microbe-iron interactions control lignin decomposition in soil

Lignin decomposition is critically linked to terrestrial carbon (C) cycle due to the enormous C mass of lignin and its importance in controlling overall rates of litter decomposition. Interactions between lignin and iron (Fe) minerals have been increasingly recognized as key mediators of lignin decomposition in experimental studies. However, we still lack a quantitative understanding of how Fe minerals interact with microbes to control lignin decomposition. Here, we leveraged experimental results from an incubation of Fe-rich soil, in which lignin decomposition rates were measured at aerobic conditions after four levels of pre-treated O2 availability, to examine microbe-Fe (MiFe) interactions in lignin decomposition with a MiFe model. We quantified how Fe redox cycling interacted with microbial activities to control lignin decomposition via data-model integration. Our results showed that the MiFe model with time-dependent growth and mortality functions better represented CO2 release from lignin decomposition (R2 ranging from 0.96 to 0.97) than models assuming either first-order or Michaelis-Menten kinetics. Reduction of Fe(III) to Fe(II) after pre-treatments with lower O2 availability stimu-lated the Fenton reaction to break down macro-molecular lignin into small molecules available to microbes. The small molecules of lignin and necromass bounded with oxidized Fe and were protected from decomposition. After 1-year incubation, the model implied that most of C stabilized with Fe minerals was derived from small molecular lignin C. Our quantitative analysis of microbe-Fe interactions sheds new light on lignin decomposition and preservation and helps improve model prediction of soil C persistence under global change.

Liao, Cuijuan↗

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗

Development of a residual acceleration data reduction and dissemination plan

A major obstacle in evaluating the residual acceleration environment in an orbiting space laboratory is the amount of data collected during a given mission: gigabytes of data will be available as SAMS units begin to fly regularly. Investigators taking advantage of the reduced gravity conditions of space should not be overwhelmed by the accelerometer data which describe these conditions. We are therefore developing a data reduction and analysis plan that will allow principal investigators of low-g experiments to create experiment-specific residual acceleration data bases for post-flight analysis. The basic aspects of the plan can also be used to characterize the acceleration environment of earth orbiting laboratories. Our development of the reduction plan is based on the following program of research: the identification of experiment sensitivities by order of magnitude estimates and numerical modelling; evaluation of various signal processing techniques appropriate for the reduction, supplementation, and dissemination of residual acceleration data; and testing and implementation of the plan on existing acceleration data bases. The orientation of the residual acceleration vector with respect to some set of coordinate axes is important for experiments with known directional sensitivity. Orientation information can be obtained from the evaluation of direction cosines. Fourier analysis is commonly used to transform time history data into the frequency domain. Common spectral representations are the amplitude spectrum which gives the average of the components of the time series at each frequency and the power spectral density which indicates the power or energy present in the series per unit frequency interval. The data reduction and analysis scheme developed involves a two tiered structure to: (1) identify experiment characteristics and mission events that can be used to limit the amount of accelerator data an investigator should be interested in; and (2) process the data in a way that will be meaningful to the experiment objectives. A general outline of the plan is given.

Rogers, Melissa J. B.↗

Simplified Modeling of Oxidation of Hydrocarbons

A method of simplified computational modeling of oxidation of hydrocarbons is undergoing development. This is one of several developments needed to enable accurate computational simulation of turbulent, chemically reacting flows. At present, accurate computational simulation of such flows is difficult or impossible in most cases because (1) the numbers of grid points needed for adequate spatial resolution of turbulent flows in realistically complex geometries are beyond the capabilities of typical supercomputers now in use and (2) the combustion of typical hydrocarbons proceeds through decomposition into hundreds of molecular species interacting through thousands of reactions. Hence, the combination of detailed reaction- rate models with the fundamental flow equations yields flow models that are computationally prohibitive. Hence, further, a reduction of at least an order of magnitude in the dimension of reaction kinetics is one of the prerequisites for feasibility of computational simulation of turbulent, chemically reacting flows. In the present method of simplified modeling, all molecular species involved in the oxidation of hydrocarbons are classified as either light or heavy; heavy molecules are those having 3 or more carbon atoms. The light molecules are not subject to meaningful decomposition, and the heavy molecules are considered to decompose into only 13 specified constituent radicals, a few of which are listed in the table. One constructs a reduced-order model, suitable for use in estimating the release of heat and the evolution of temperature in combustion, from a base comprising the 13 constituent radicals plus a total of 26 other species that include the light molecules and related light free radicals. Then rather than following all possible species through their reaction coordinates, one follows only the reduced set of reaction coordinates of the base. The behavior of the base was examined in test computational simulations of the combustion of heptane in a stirred reactor at various initial pressures ranging from 0.1 to 6 MPa. Most of the simulations were performed for stoichiometric mixtures; some were performed for fuel/oxygen mole ratios of 1/2 and 2.

Bellan, Josette↗

A comparison of neural network architectures for data-driven reduced-order modeling

The popularity of deep convolutional autoencoders (CAEs) has engendered new and effective reduced-order models (ROMs) for the simulation of large-scale dynamical systems. Despite this, it is still unknown whether deep CAEs provide superior performance over established linear techniques or other network-based methods in all modeling scenarios. To elucidate this, the effect of autoencoder architecture on its associated ROM is studied through the comparison of deep CAEs against two alternatives: a simple fully connected autoencoder, and a novel graph convolutional autoencoder. Through benchmark experiments, it is shown that the superior autoencoder architecture for a given ROM application is highly dependent on the size of the latent space and the structure of the snapshot data, with the proposed architecture demonstrating benefits on data with irregular connectivity when the latent space is sufficiently large.

42 ENGINEERING↗

A Scalable Space-Time Domain Decomposition Approach for Solving Large Scale Nonlinear Regularized Inverse Ill Posed Problems in 4D Variational Data Assimilation

We address the development of innovative algorithms designed to solve the strong-constraint Four Dimensional Variational Data Assimilation (4DVar DA) problems in large scale applications. We present a space-time decomposition approach which employs the whole domain decomposition, i.e. both along the spacial and temporal direction in the overlapping case, and the partitioning of both the solution and the operator. Starting from the global functional defined on the entire domain, we get to a sort of regularized local functionals on the set of sub domains providing the order reduction of both the predictive and the Data Assimilation models. The algorithm convergence is developed. Performance in terms of reduction of time complexity and algorithmic scalability is discussed on the Shallow Water Equations on the sphere. The number of state variables in the model, the number of observations in an assimilation cycle, as well as numerical parameters as the discretization step in time and in space domain are defined on the basis of discretization grid used by data available at repository Ocean Synthesis/Reanalysis Directory of Hamburg University.

97 MATHEMATICS AND COMPUTING↗

Design of turbofan engine controls using output feedback regulator theory

A multivariable control design procedure based on output feedback regulator (OFR) theory is applied to the F100 turbofan engine. Results for the OFR design are compared to a design based on linear quadratic regulator (LQR) theory. The OFR feedback control is designed in the full order state space and thus eliminates any need for model reduction techniques. Using the performance measure and control structure of the LQR design, an equivalent OFR feedback control is obtained. The flexibility of the OFR as a control design procedure is demonstrated, and differing feedback control structures are evaluated.

Merrill, W. C.↗

Design of turbofan engine controls using output feedback regulator theory

A multivariable control design procedure based on output feedback regulator (OFR) theory is applied to the F100 turbofan engine. Results for the OFR design are compared to a design based on linear quadratic regulator (LQR) theory. This LQR design was obtained as part of the F100 Multivariable Control Synthesis (MVCS) program. In the MVCS program the LQR feedback control was designed in a reduced dimension state space and then applied to the original system. However, the OFR feedback control is designed in the full order state space and thus eliminates any need for model reduction techniques. Using the performance measure and control structure of the MVCS program LQR design, an equivalent OFR feedback control is obtained. The flexibility of the OFR as a control design procedure is demonstrated and differing feedback control structures are evaluated.

Merrill, W. C.↗

Alternatives for jet engine control

The technical progress of researches on alternatives for jet engine control, is reported. The principal new activities involved the initial testing of an input design method for choosing the inputs to a non-linear system to aid the approximation of its tensor parameters, and the beginning of order reduction studies designed to remove unnecessary monomials from tensor models.

Sain, M. K.↗

Effect of Temperature Gradient on Thick Film Selective Emitter Emittance

A temperature gradient across a thick (greater than or equal to .1 mm) film selective emitter will produce a significant reduction in the spectral emittance from the no temperature gradient case. Thick film selective emitters of rare earth doped host materials such as yttrium-aluminum-garnet (YAG) are examples where temperature gradient effects are important. In this paper a model is developed for the spectral emittance assuming a linear temperature gradient across the film. Results of the model indicate that temperature gradients will result in reductions the order of 20% or more in the spectral emittance.

Chubb, Donald L.↗

A Physics-Based Data-Driven Approach for Modeling of Environmental Degradation in Elastomers

Abstract Elastomers are now commonly used in a number of industries, including aerospace, structure, transportation, shipbuilding, and automotive, due to their excellent workability, formability, and flexibility. During their activity, elastomers are subjected to harsh environmental conditions, which decreases their resilience. False predictions made early in their lives can have major financial and environmental implications. Elastomers’ performance and properties, such as strength, durability, and density, are influenced by chemical changes in these materials, known as degradation, which occurs over time. This process can alter the morphology of a polymer matrix as well as cause chain scission and cross-linking, resulting in different behaviors than that of the unaged material. To demonstrate the effect of thermaloxidative aging on the mechanical behavior of elastomers, several experimental and theoretical models have been proposed. In view of the large volume of experimental data available on micro-structural evolution in the course of aging, we propose a physics-based data-driven approach to overcome the shortcomings of both phenomenological and micro-mechanical models. This work presents a novel thermodynamically consistent, multiagent machine-learned model for predicting the constitutive behavior of cross-linked elastomers during environmental aging, such as thermo-oxidative and hydrolytic aging for various states of deformation. Single mechanism degradation changes the polymer matrix over time where it is causing chain scission, reduction of cross-links, and morphology change. To capture the idealized Mullins effect and permanent set due to the effect of single aging mechanisms on nonlinear mechanical responses of elastomers, we propose a data-driven model for simulating inelastic elements in a polymer matrix. By using a sequential order reduction, we were able to reduce the 3D stress-strain tensor mapping problem to a small number of super-constrained 1D mapping problems. To systematically classify such mapping problems into a few categories, an assembly of multiple replicated conditional neural network learning agents (L-agents) is used based on our recent work. Each category is represented by a different type of agent. The effect of deformation history, aging time, and aging temperature is captured by this model. The model is validated using a broad collection of data, ranging from our experimental results to data from the literature. In addition, thermodynamic consistency and frame independence are investigated. The most significant achievements of this model are its precision, simplicity, and prediction of inelasticity under various states of deformation. The model’s accuracy and simplicity make it a good option for commercial and industrial applications. Conveniently, due to the model modular nature, it can be expanded in the future to include viscoelasticity and non-isotropic formation for better precision.

Ghaderi, Aref↗

Ozone reduction in the 1980's - A model simulation of anthropogenic and solar perturbations

An interpretation of global ozone changes deduced from satellite data obtained since 1979 is presented, based on two-dimensional model simulations. The study shows that a depletion in total ozone of the order of 2 percent and a reduction in ozone density near 40 km of 7 to 12 percent over the 1979-1986 period are consistent with the observed increase in trace gas densities (chlorofluorocarbons, methane, nitrous oxide, carbon dioxide) and the simultaneous decrease in solar activity during this period. The model also suggests that ozone variations of solar and anthropogenic origins between 1979 and 1986 were of similar magnitude but that the ozone response to trace gas emissions increases substantially with latitude while the solar signal in ozone is present lower in the atmosphere and is nearly independent of latitude.

Brasseur, Guy↗

Design of an Axisymmetric Afterbody Test Case for CFD Validation

As identified in the CFD Vision 2030 Study commissioned by NASA, validation of advanced RANS models and scale-resolving methods for computing turbulent flow fields must be supported by continuous improvements in fundamental, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, the underpinnings of a new test platform referred to herein as the NASA Axisymmetric Afterbody are presented. The devised body-of-revolution is a modular platform consisting of a forebody section and afterbody section, allowing for a range of flow behaviors to be studied on interchangeable afterbody geometries. A body-of-revolution offers advantages in shape definition and fabrication, in avoiding direct contact with wind tunnel sidewalls, and in tail-sting integration to facilitate access to higher Reynolds number tunnels. The current work is focused on validation of smooth-body turbulent flow separation, for which a six-parameter body has been developed. A priori RANS computations are reported for a risk-reduction test configuration in order to demonstrate critical variation among turbulence model results for a given afterbody, ranging from barely-attached to mild separated flow. RANS studies of the effects of forebody nose (with/without) and wind tunnel boundary (slip/no-slip) on the selected afterbody are presented. Representative modeling issues that can be explored with this configuration are the effect of higher Reynolds number on separation behavior, flow physics of the progression from attached to increasingly-separated afterbody flows, and the effect of embedded longitudinal vortices on turbulence structure.

Disotell, Kevin J.↗

Robust Output Feedback Control Design for Inertia Emulation by Wind Turbine Generators

Wind generation has gained widespread use as a renewable energy source. Most wind turbines and other renewables connected to the grid through converters result in a reduction in the natural inertial response to grid frequency changes. The doubly-fed induction generator (DFIG) can be controlled to compensate for this reduction and, in fact, provide faster response than traditional synchronous machines. This paper proposes to design observer based output feedback linear quadratic regulator (LQR) and H control laws to realize the inertia emulation function and deliver fast frequency support. Furthermore, the aim is to track the reference speed by a diesel synchronous generator (DSG) in order to reach the desired inertia. The control signal is computed based on a reduced order model using the balanced truncation technique. A comparison with selective modal analysis (SMA) and balanced truncation model reduction techniques is presented. Comprehensive results show the effective emulation of synthetic inertia by implementing the control laws on a nonlinear three- phase diesel-wind system. The proposed technique is analyzed for different short circuit ratio (SCR) scenarios.

17 WIND ENERGY↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Study of CORC Conductors With Respect to Individual Tape Properties

CEA has been studying the advantages of a conductor based on an assembly of CORC cables in the high field zone of a hybrid Central Solenoid (CS) magnet for EU-DEMO. To this end, the detailed study of geometrical and electrical parameters of a CORC structure are studied in order to evaluate the cable’s electrical performance as well as to determine a number of important parameters (crossing points, contact surface etc…) as function of the cable structure. The paper first presents these geometrical and performance analyses. It will then try to introduce smeared models in order to reduce the cable’s performance to 1D tape scaling law using effective parameters. That reduction is of importance for use in thermohydraulic models. Finally, the paper will present the case study of a particular CORC conductor that is being procured and is foreseen to be tested in SULTAN in 2025. The paper also will try to give some predictive estimate of the cable performance and identify some of the unknowns related to current redistribution and joint resistance.

CORC↗

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↗