Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 451 records · Page 25

Space Station flexible dynamics under plume impingement

Assembly of the Space Station requires numerous construction flights by the Space Shuttle. A particularly challenging problem is that of control of each intermediate station configuration when the shuttle orbiter is approaching it to deliver the next component. The necessary braking maneuvers cause orbiter thruster plumes to impinge on the station, especially its solar arrays. This in turn causes both overall attitude errors and excitation of flexible-body vibration modes. These plume loads are predicted to lead to CMG saturation during the approach of the orbiter to the SC-5 station configuration, necessitating the use of the station RCS jets for desaturation. They are also expected to lead to significant excitation of solar array vibrations. It is therefore of great practical importance to investigate the effects of plume loads on the flexible dynamics of station configuration SC-5 as accurately as possible. However, this system possesses a great many flexible modes (89 below 5 rad/s), making analysis time-consuming and complicated. Model reduction techniques can be used to overcome this problem, reducing the system model to one which retains only the significant dynamics, i.e. those which are strongly excited by the control inputs or plume disturbance forces and which strongly couple with the measured outputs. The particular technique to be used in this study is the subsystem balancing approach which was previously developed by the present investigator. This method is very efficient computationally. Furthermore, it gives accurate results even for the difficult case where the structure has many closed-spaced natural frequencies, when standard modal truncation can give misleading results. Station configuration SC-5 is a good example of such a structure.

Williams, Trevor↗

Canonical and noncanonical Hamiltonian operator inference

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

97 MATHEMATICS AND COMPUTING↗

A New Approach to Overcoming Spatial Aliasing in Structural Damage Detection

Aircraft, reusable launch vehicles,unmanned aircraft, and other advanced structures are being built using lightweight composite materials/metals with design safety factors as low as 1.25. These advanced structures operate in uncertain and severe environments and are susceptible to damage such as delamination, fiber/matrix damage, hydrothermal strain in composite materials, and fatigue and cracking in metals. To ensure human safety and load-bearing integrity these structures must be inspected to detect and locate often invisible damage and faults before they become catastrophic. Conventional methods of non-destructive evaluation sometimes miss significant damage and are time consuming and expensive to perform. In contrast, vibrometry or vibration signature techniques are a global method of structural integrity monitoring that potentially can efficiently detect damage on large structures, including damage that is away from sensor locations, and in the interior of structures. However, a barrier problem in damage detection using vibration measurements is the need to measure the vibration response at a large number of points on the structure. Typically, model reduction or expansion procedures such as Guyan reduction or dynamic expansion are attempted to overcome the problem of insufficient measurements or spatial aliasing. These approaches depend on using information from the healthy model, and thus put error into the reduction used to represent the damaged structure. In this paper, a Frequency Response Function technique is used to detect damage to a fixed-free beam. The technique uses measured frequency response functions from the healthy structure as reference data, and then monitors vibration measurements during the life of the structure to detect damage. In an analytical simulation using a finite-element model of a beam, damage was located using only sparse measurements because the technique uses both rotation and translation measurements from the damaged structure.

Schulz, Mark J.↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

Gust Load Alleviation Control and Gust Estimation for a High Aspect Ratio Wing Wind Tunnel Model

This paper presents the gust load alleviation (GLA) study of the aspect ratio 13.5 Common Research Model (CRM) wind-tunnel model. This study details the design of the GLA controller in preparation for wind-tunnel testing in the Transonic Dynamics Tunnel at NASA Langley Research Center. An aeroservolastic (ASE) model was first reduced using a model reduction method that takes advantage of the sinusoidal steady-state response. Then, the reduced model was used to design an extended-state Kalman filter which estimates the states and the sinusoidal gust input. The GLA control was then derived using the optimal control solution to a multi-objective cost function. The results of the GLA controller indicate a 75.6% reduction in wing-root strain while maintaining robust stability margins. The final paper will include performance under sub-optimal conditions, e.g., sensor noise, uncertainty in gust frequency, etc..

Christopher Forte↗

Gust Load Alleviation Control and Gust Estimation for a High Aspect Ratio Wing Wind Tunnel Model

This paper presents the gust load alleviation (GLA) study of the aspect ratio 13.5 Common Research Model (CRM) wind-tunnel model. This study details the design of the GLA controller in preparation for wind-tunnel testing in the Transonic Dynamics Tunnel at NASA Langley Research Center. An aeroservolastic (ASE) model is first reduced using a model reduction method that takes advantage of the sinusoidal steady-state response. Then, the reduced model is used to design an extended-state Kalman filter which estimates the states and the sinusoidal gust input. The GLA control is then derived using the optimal control solution to a multi-objective cost function. The results of the GLA controller indicate a 69.07% reduction in wing-root strain without sensor noise and 68.45% reduction with sensor noise, while maintaining robust stability margins. An adaptive GLA controller is developed for uncertain gust frequency and shows a 71.03% reduction in wing root strain compared to the non-adaptive control reduction of just 39.19%.

Christopher Forte↗

Nonlinear proper orthogonal decomposition for convection-dominated flows

Autoencoder techniques find increasingly common use in reduced order modeling as a means to create a latent space. This reduced order representation offers a modular data-driven modeling approach for nonlinear dynamical systems when integrated with a time series predictive model. In this Letter, we put forth a nonlinear proper orthogonal decomposition (POD) framework, which is an end-to-end Galerkin-free model combining autoencoders with long short-term memory networks for dynamics. By eliminating the projection error due to the truncation of Galerkin models, a key enabler of the proposed nonintrusive approach is the kinematic construction of a nonlinear mapping between the full-rank expansion of the POD coefficients and the latent space where the dynamics evolve. We test our framework for model reduction of a convection-dominated system, which is generally challenging for reduced order models. Our approach not only improves the accuracy, but also significantly reduces the computational cost of training and testing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Control of the SCOLE configuration using distributed parameter models

A continuum model for the SCOLE configuration has been derived using transfer matrices. Controller designs for distributed parameter systems have been analyzed. Pole-assignment controller design is considered easy to implement but stability is not guaranteed. An explicit transfer function of dynamic controllers has been obtained and no model reduction is required before the controller is realized. One specific LQG controller for continuum models had been derived, but other optimal controllers for more general performances need to be studied.

Hsiao, Min-Hung↗

Improved physiologically based pharmacokinetic model for oral exposures to chromium in mice, rats, and humans to address temporal variation and sensitive populations

A physiologically based pharmacokinetic (PBPK) model for hexavalent chromium [Cr(VI)] in mice, rats, and humans developed previously (Kirman et al., 2012, 2013), was updated to reflect an improved understanding of the toxicokinetics of the gastrointestinal tract following oral exposures. Improvements were made to: (1) the reduction model, which describes the pH-dependent reduction of Cr(VI) to Cr(III) in the gastrointestinal tract under both fasted and fed states; (2) drinking water pattern simulations, to better describe dosimetry in rodents under the conditions of the NTP cancer bioassay; and (3) parameterize the model to characterize potentially sensitive human populations. Important species differences, sources of non-linear toxicokinetics, and human variation are identified and discussed within the context of human health risk assessment. - Highlights: • An improved version of the PBPK model for Cr(VI) toxicokinetics was developed. • The model incorporates data collected to fill important data gaps. • Model predictions for specific age groups and sensitive subpopulations are provided. • Implications to human health risk assessment are discussed.

60 APPLIED LIFE SCIENCES↗

Computer-Aided Construction of Chemical Kinetic Models

The combustion chemistry of even simple fuels can be extremely complex, involving hundreds or thousands of kinetically significant species. The most reasonable way to deal with this complexity is to use a computer not only to numerically solve the kinetic model, but also to construct the kinetic model in the first place. Because these large models contain so many numerical parameters (e.g. rate coefficients, thermochemistry) one never has sufficient data to uniquely determine them all experimentally. Instead one must work in “predictive” mode, using theoretical rather than experimental values for many of the numbers in the model, and as appropriate refining the most sensitive numbers through experiments. Predictive chemical kinetics is exactly what is needed for computer-aided design of combustion systems based on proposed alternative fuels, particularly for early assessment of the value and viability of proposed new fuels before those fuels are commercially available. This project was aimed at making accurate predictive chemical kinetics practical; this is a challenging goal which requires a range of science advances. The project spanned a wide range from quantum chemical calculations on individual molecules and elementary-step reactions, through the development of improved rate/thermo calculation procedures, the creation of algorithms and software for constructing and solving kinetic simulations, the invention of methods for model-reduction while maintaining error control, and finally comparisons with experiment. Many of the parameters in the models were derived from quantum chemistry calculations, and the models were compared with experimental data measured in our lab or in collaboration with others.

97 MATHEMATICS AND COMPUTING↗

High Penetration Solar PV Deployment Sunshine State Solar Grid Initiative (SUNGRIN)

The report provides results from the Sunshine State Solar Grid Initiative (SUNGRIN) high penetration solar PV deployment project led by Florida State University’s (FSU) Center for Advanced Power Systems (CAPS). FSU CAPS and industry and university partners have completed a five-year effort aimed at enabling effective integration of high penetration levels of grid-connected solar PV generation. SUNGRIN has made significant contributions in the development of simulation-assisted techniques, tools, insight and understanding associated with solar PV effects on electric power system (EPS) operation and the evaluation of mitigation options for maintaining reliable operation. An important element of the project was the partnership and participation of six major Florida utilities and the Florida Reliability Coordinating Council (FRCC). Utilities provided details and data associated with actual distribution circuits having high-penetration PV to use as case studies. The project also conducted foundational work supporting future investigations of effects at the transmission / bulk power system level. In the final phase of the project, four open-use models with built-in case studies were developed and released, along with synthetic solar PV data sets, and tools and techniques for model reduction and in-depth parametric studies of solar PV impact on distribution circuits. Along with models and data, at least 70 supporting MATLAB functions have been developed and made available, with complete documentation.

14 SOLAR ENERGY↗

Modeling and control of beam-like structures

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for beam-like continua. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

Modal cost analysis for simple continua

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations, it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for simple continua such as beam-like structures. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

Modeling and simulations of the DSS 13 antenna control system

A model of the antenna control system for the azimuth and elevation axes of the Deep Space Station (DSS) 13 antenna is developed. This model is used for simulation of elevation and azimuth dynamics, cross-coupled dynamics, and radio-frequency pointing error due to both input commands and wind disturbances. This model also serves as a tool for the antenna controller design. A modal state-space model of the antenna structure was obtained from its finite-element model with a free-rotating tipping structure and alidade. Model reduction techniques applied separately for the antenna structure, elevation, and azimuth drives, and rate-loop model reduce the system order to one-third of that of the original, while preserving its significant dynamic properties. Extensive simulation results illustrate properties of the model.

Gawronski, W.↗

Multi-Physics Surrogate Modeling with Dimensionality Reduction

This summary provides a new algorithm to construct surrogate models for multi-physics, nonlinear, high dimensional models which typically suffer from the curse of dimensionality. To deal with the high computational cost and difficult error estimation, we employ recent advances in reduced order modeling (ROM) techniques to do an error-preserving reduction for the number of influential model parameters, followed by a conventional surrogate construction process. To implement ROM for a multiphysics model, a range finding algorithm (RFA) is used to identify a linear transformation in which the number of influential parameters can be minimized and the errors resulting from the non-influential parameters can be upper-bounded. The proposed algorithm is applied to the MAMMOTH coupled code system developed at the Idaho National Laboratory. In particular, we investigate the reducibility of Rattlesnake and BISON codes and their multi-physics coupling. The most important result of this work is that the reduction obtained with coupled physics is more significant than that with individual physics models. This summary explores the application of ROM techniques on multi-physics coupling codes system. Concise methodology is provided and selected numerical results are provided to show that coupling system provides additional reduction at interfaces. Future work could be invested to determine the optimum surrogate model format. (authors)

97 MATHEMATICS AND COMPUTING↗

Modeling Fire Severity in Black Spruce Stands in the Alaskan Boreal Forest Using Spectral and Non-Spectral Geospatial Data

Biomass burning in the Alaskan interior is already a major disturbance and source of carbon emissions, and is likely to increase in response to the warming and drying predicted for the future climate. In addition to quantifying changes to the spatial and temporal patterns of burned areas, observing variations in severity is the key to studying the impact of changes to the fire regime on carbon cycling, energy budgets, and post-fire succession. Remote sensing indices of fire severity have not consistently been well-correlated with in situ observations of important severity characteristics in Alaskan black spruce stands, including depth of burning of the surface organic layer. The incorporation of ancillary data such as in situ observations and GIS layers with spectral data from Landsat TM/ETM+ greatly improved efforts to map the reduction of the organic layer in burned black spruce stands. Using a regression tree approach, the R2 of the organic layer depth reduction models was 0.60 and 0.55 (pb0.01) for relative and absolute depth reduction, respectively. All of the independent variables used by the regression tree to estimate burn depth can be obtained independently of field observations. Implementation of a gradient boosting algorithm improved the R2 to 0.80 and 0.79 (pb0.01) for absolute and relative organic layer depth reduction, respectively. Independent variables used in the regression tree model of burn depth included topographic position, remote sensing indices related to soil and vegetation characteristics, timing of the fire event, and meteorological data. Post-fire organic layer depth characteristics are determined for a large (N200,000 ha) fire to identify areas that are potentially vulnerable to a shift in post-fire succession. This application showed that 12% of this fire event experienced fire severe enough to support a change in post-fire succession. We conclude that non-parametric models and ancillary data are useful in the modeling of the surface organic layer fire depth. Because quantitative differences in post-fire surface characteristics do not directly influence spectral properties, these modeling techniques provide better information than the use of remote sensing data alone.

Barrett, K.↗

Heating-Rate-Coupled Model for Hydrogen Reduction of JSC-1A

A previously developed and validated model for hydrogen reduction of JSC-1A for a constant reaction-bed temperature is extended to account for reaction during the bed heat-up period. A quasisteady approximation is used wherein an expression is derived for a single average temperature of reaction during the heat-up process by employing an Arrhenius expression for regolith conversion. Subsequently, the regolith conversion during the heat-up period is obtained by using this representative temperature. Accounting for the reaction during heat-up provides a better estimate of the reaction time needed at the desired regolith-bed operating temperature. Implications for the efficiency of the process, as measured by the energy required per unit mass of oxygen produced, are also indicated.

Hegde, U.↗