Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order Modeling”

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 775 records · Page 43

GrainNN: A neighbor-aware long short-term memory network for predicting microstructure evolution during polycrystalline grain formation

High fidelity simulations of grain formation in alloys are an indispensable tool for process-to-mechanical-properties characterization. Such simulations, however, can be computationally expensive as they require fine spatial and temporal discretizations. Their cost becomes an obstacle to parametric studies and ensemble runs and ultimately makes downstream tasks like optimal control and uncertainty quantification challenging. To enable such downstream tasks, we introduce GrainNN, an efficient and accurate reduced-order model for epitaxial grain growth in additive manufacturing conditions. GrainNN is a sequence-to-sequence long-short-term-memory (LSTM) deep neural network that evolves the dynamics of manually crafted features. Its innovations are (1) an attention mechanism with grain-microstructure-specific transformer architecture; and (2) an overlapping combination of several clones of the network to generalize to grain configurations that are different from those used for training. This design enables GrainNN to predict grain formation for unseen physical parameters, grain number, domain size and geometry. Furthermore, GrainNN not only reconstructs the quantities of interest but also can be pointwise accurate. In our numerical experiments, we use a polycrystalline phase field method to both generate the training data and assess GrainNN. For multiparametric, ensemble simulations with many grains, GrainNN can be orders of magnitude faster than phase field simulations, while delivering 5%–15% pointwise error. Additionally, this speedup includes the cost of the phase field simulations for generating training data.

36 MATERIALS SCIENCE↗

Statistical Correlation of Heliostat Pointing Deviation With Wind

This work was carried out as part of the Heliostat Consortium (HelioCon) Field Deployment subtask with the aim to develop a reduced order model framework for correlating wind speed and pointing deviation of a heliostat facet. There are only sparse field measurements of heliostat pointing deviations and accompanying wind conditions published in the literature. Heliostat test standards, such as IEC 62862-4-3, propose a suite of tests including laser pointing repeatability at wind speeds below 4 m/s, and provide technical requirements for heliostat slope and tracking deviations in coarse average wind speed bins of 4 m/s, 6 m/s, and 8 m/s. In addressing the gap of the variation of heliostat pointing deviation with wind speed, field measurements of laser pointing on a grid target and wind conditions were analyzed in this study at the Third-Party Metrology Platform at the National Laboratory of the Rockies (NLR) Flatirons Campus. Horizontal pointing deviations were found to follow a logarithmic relationship with peak wind speed, whereas vertical pointing deviations follow an exponential relationship with peak wind speed. Both horizontal and vertical pointing deviations also follow a second order polynomial relationship, as expected from the proportionality of elastic loads and deformations with the square of wind speed. The results indicate that heliostat facet pointing deviations in the vertical direction increase at a faster rate than in the horizontal direction with increasing wind speed over the tested range, however these are dependent on the heliostat structural design. Next steps are recommended for additional field measurements to confirm a linear relationship of pointing deviation with applied moment on a heliostat facet, and to distinguish between gravity-induced and wind-induced pointing deviations at different elevation angles. The derived correlations in the preliminary analysis in this report serve as a case study for heliostat developers and plant operators to estimate the wind-induced pointing deviations and their variation with peak gust wind speed. Next steps in future work would recommend higher resolution and longer duration datasets for different elevation angles and wind directions to reduce uncertainties and variance of collected laser beam spot data and their correlations with bin-averaged wind speed.

17 WIND ENERGY↗

AeroFusion: Data Fusion and Uncertainty Quantification for Entry Vehicles

AeroFusion is a NASA Langley initiative to incorporate advances in data science into the aerodynamic modeling process to improve efficiency. The effort can largely be categorized in three components: reduced-order modeling techniques, surrogate modeling techniques, and uncertainty quantification. By combining various methods from these categories, AeroFusion aims to reduce the development cost of aerodynamic models, both in terms of time and money.

Steven Snyder↗

AeroFusion: Data Fusion and Uncertainty Quantification for Entry Vehicles

AeroFusion is a NASA Langley initiative to incorporate advances in data science into the aerodynamic modeling process to improve efficiency. The effort can largely be categorized in three components: reduced-order modeling techniques, surrogate modeling techniques, and uncertainty quantification. By combining various methods from these categories, AeroFusion aims to reduce the development cost of aerodynamic models, both in terms of time and money.

Steven Snyder↗

Krylov model reduction algorithm for undamped structural dynamics systems

Krylov vectors furnish an efficient basis for eigenvalue analysis and model reduction of structural dynamics systems. The reduced-order model obtained by the present Krylov model-reduction algorithm for an undamped structural-dynamics system is found to match low-frequency moments. The transformed system equation in Krylov coordinates reflects the structure of a tandem system.

Craig, Roy R., Jr.↗

Grizzly: A Multi-scale and Multi-Physics Tool to Model Aging of Nuclear Power Plant Components

Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electrochemical state and internal variables estimation using a reduced-order physics-based model of a lithium-ion cell and an extended Kalman filter

This paper addresses the problem of estimating the present value of electrochemical internal variables in a lithium-ion cell in real time, using readily available measurements of cell voltage, current, and temperature. The variables that can be estimated include any desired set of reaction flux and solid and electrolyte potentials and concentrations at any set of one-dimensional spatial locations, in addition to more standard quantities such as state of charge. The method uses an extended Kalman filter along with a one-dimensional physics-based reduced-order model of cell dynamics. Simulations show excellent and robust predictions having dependable error bounds for most internal variables. (C) 2014 Elsevier B.V. All rights reserved.

Stetzel, KD↗

On Computationally Efficient Radiative Transfer Calculations for Three-Dimensional Entry Problems

The current work presents an efficient simulation framework for rigorously modeling radiative fields emanating from non-equilibrium planetary entry flows in complex three-dimensional domains. Key to this endeavor is adoption of finite-volume discretization in lieu of brute-force ray tracing. This change in conjunction with mesh sweeping and Lebedev-type quadrature for angular integration allows spatial-angular resolution of radiative transfer to be performed in a computationally tractable manner. Additionally, a new methodology has been established to build standalone reduced-order spectral databases for non-equilibrium radiative properties that can be applied to a broad range of hypersonic planetary entry problems with minimal loss in accuracy. The efficacy of the new framework has been demonstrated on the atomic nitrogen radiative system. The resulting reduced-order model requires three orders-of-magnitude fewer spectral evaluations and results in a 95% decrease in memory footprint. A comparative study for representative forebody and afterbody lines-of-sight from Stardust, FIRE II, and meteor entries into the Earth atmosphere indicates that both total intensity variation and detailed spectra can be retrieved with as few as 625 reduced-order groups (contrasting with the 100,000 frequencies in the original full set model). Similarly, three-dimensional predictions of radiative heating experienced by the Orion forebody are in excellent agreement with legacy radiation solvers while requiring only a sliver (roughly 0.5%) of computing wall time.

3D radiation↗

On Computationally Efficient Radiative Transfer Calculations for Three-dimensional Entry Problems

The current work presents an efficient simulation framework for rigorously modeling radiative fields emanating from non-equilibrium planetary entry flows in complex three-dimensional domains. Key to this endeavor is adoption of finite-volume discretization in lieu of brute-force ray tracing. This change in conjunction with mesh sweeping and Lebedev-type quadrature for angular integration allows spatial-angular resolution of radiative transfer to be performed in a computationally tractable manner. Additionally, a new methodology has been established to build standalone reduced-order spectral databases for non-equilibrium radiative properties that can be applied to a broad range of hypersonic planetary entry problems with minimal loss in accuracy. The efficacy of the new framework has been demonstrated on the atomic nitrogen radiative system. The resulting reduced-order model requires three orders-of-magnitude fewer spectral evaluations and results in a 95% decrease in memory footprint. A comparative study for representative forebody and afterbody lines-of-sight from Stardust, FIRE II, and meteor entries into the Earth atmosphere indicates that both total intensity variation and detailed spectra can be retrieved with as few as 625 reduced-order groups (contrasting with the 100,000 frequencies in the original full set model). Similarly, three-dimensional predictions of radiative heating experienced by the Orion forebody are in excellent agreement with legacy radiation solvers while requiring only a sliver (roughly 0.5%) of computing wall time.

3D radiation↗

A Hybrid Approach for Model Order Reduction of Barotropic Quasi-Geostrophic Turbulence

We put forth a robust reduced-order modeling approach for near real-time prediction of mesoscale flows. In our hybrid-modeling framework, we combine physics-based projection methods with neural network closures to account for truncated modes. We introduce a weighting parameter between the Galerkin projection and extreme learning machine models and explore its effectiveness, accuracy and generalizability. To illustrate the success of the proposed modeling paradigm, we predict both the mean flow pattern and the time series response of a single-layer quasi-geostrophic ocean model, which is a simplified prototype for wind-driven general circulation models. We demonstrate that our approach yields significant improvements over both the standard Galerkin projection and fully non-intrusive neural network methods with a negligible computational overhead.

42 ENGINEERING↗

A reduced basis approach for PDEs on parametrized geometries based on the shifted boundary finite element method and application to a Stokes flow

Here we propose a model order reduction technique integrating the Shifted Boundary Method (SBM) with a POD-Galerkin strategy. This approach allows to deal with complex parametrized domains in an efficient and straightforward way. The impact of the proposed approach is threefold. First, problems involving parametrizations of complex geometrical shapes and/or large domain deformations can be efficiently solved at full-order by means of the SBM. This unfitted boundary method permits to avoid remeshing and the tedious handling of cut cells by introducing an approximate surrogate boundary. Second, the computational effort is reduced by the development of a Reduced Order Model (ROM) technique based on a POD-Galerkin approach. Third, the SBM provides a smooth mapping from the true to the surrogate domain, and for this reason, the stability and performance of the reduced order basis are enhanced. This feature is the net result of the combination of the proposed ROM approach and the SBM. Similarly, the combination of the SBM with a projection-based ROM gives the great advantage of an easy and fast to implement algorithm considering geometrical parametrization with large deformations. The transformation of each geometry to a reference geometry (morphing) is in fact not required. These combined advantages will allow the solution of PDE problems more efficiently. We illustrate the performance of this approach on a number of two-dimensional Stokes flow problems.

42 ENGINEERING↗

Shock-Stationary Application of Pseudoshock Models During High-Amplitude Combustion-Driven Unsteadiness

The isolator pseudo-shock provides necessary compression within a dual-mode scramjet engine and buffers the engine system against unstart. Quasi-1D flux-conserved models are the state-of-the-art reduced-order model for optimization and online control of dual-mode scramjet engines. The stability and efficacy of this modeling approach is evaluated against data from a combustor-driven direct-connect experiment. The experiment exhibited strong combustor-driven unsteadiness that produced upstream propagating weak shocks into the isolator, interacting with the pseudo-shock. While this configuration resulted in unsteadiness that is atypical of standard operation, the experiment provided an opportunity to evaluate the modeling techniques in highly transient states. Such transients could occur during maneuvering or result from unexpected combustor events. A flexible quasi-1D formulation, the Fievet flux-conserved model is fit using Bayesian inference in the laboratory and shock-stationary reference frames. Model performance is analyzed using the Bayesian posteriors and model evaluations over the measured shock train speed range. It is concluded that to produce consistent isolator pressure profile estimates in this unsteady environment, the model must be implemented in a shock-stationary reference frame. Implementing this conclusion in model-based engine controllers may reduce needed unstart safety margins and increase maximum performance.

Bayesian↗

mphys-surrogate-model

This repository contains python scripts for building and studying reduced-order-modeling representations of droplet coalescence for eventual use in atmospheric models. The included data are generated from high-fidelity superdroplet methods and are utilized by machine learning pipelines to build data-driven models of droplet size distributions that evolve under coalescence. This repository further includes scripts to determine prediction (uncertainty) intervals on the data-driven model products based on conformal prediction.

Katona, JonasE [Lawrence Livermore National Labora↗

On the stability and compressive nonlinearity of a physiologically based model of the cochlea

Hearing relies on a series of coupled electrical, acoustical (fluidic) and mechanical interactions inside the cochlea that enable sound processing. A positive feedback mechanism within the cochlea, called the cochlear amplifier, provides amplitude and frequency selectivity in the mammalian auditory system. The cochlear amplifier and stability are studied using a nonlinear, micromechanical model of the Organ of Corti (OoC) coupled to the electrical potentials in the cochlear ducts. It is observed that the mechano-electrical transduction (MET) sensitivity and somatic motility of the outer hair cell (OHC), control the cochlear stability. Increasing MET sensitivity beyond a critical value, while electromechanical coupling coefficient is within a specific range, causes instability. We show that instability in this model is generated through a supercritical Hopf bifurcation. A reduced order model of the system is approximated and it is shown that the tectorial membrane (TM) transverse mode effect on the dynamics is significant while the radial mode can be simplified from the equations. The cochlear amplifier in this model exhibits good agreement with the experimental data. A comprehensive 3-dimensional model based on the cross sectional model is simulated and the results are compared. It is indicated that the global model qualitatively inherits some characteristics of the local model, but the longitudinal coupling along the cochlea shifts the stability boundary (i.e., Hopf bifurcation point) and enhances stability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Physics-Informed Machine Learning with Application to Solid Oxide Fuel Cell System Modeling and Optimization

Traditional data-driven machine learning (ML) techniques were combined with the physics-based SOFC-MP model toward improving SOFC system-level performance prediction. Four different physics-informed ML methods were proposed to demonstrate better prediction accuracy compared with traditional ML methods. First, it was preferred that the existing Kriging-based surrogate model not be entirely abandoned, so errors between Kriging regression prediction and true solution are learned using ML. The goal was to find the pattern of the error distribution as a function of model input so discrepancies between the Kriging prediction and true solution could be reduced by adding the additional ML-trained error term to the existing prediction. Second, deep neural networks (DNN) coupled with the mass balance model (MBM) have been developed to significantly decrease the reduced order model (ROM) prediction error with fewer training data compared to the traditional Kriging-based ROM and traditional DNN regression approaches. Third, DNN regression coupled with a neural networks (NN) classifier and MBM was developed. It provides superior performance on classifying the physically operational conditions for natural gas fuel cell (NGFC) than traditional classification approaches. Finally, two kinds of DNN transferring were studied and tested. One was transferring DNN for low-fidelity model data to a high-fidelity model. The second one was transferring DNN from state-of-the-art NGFC to advanced NGFC. These physic-informed ML frameworks can serve as a solid foundation for estimating U.S. Department of Energy-targeted fuel cell performance and exhibit great potential to be employed to other engineering applications.

30 DIRECT ENERGY CONVERSION↗

Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations

This paper discusses numerical aspects of computing maximum likelihood estimates for linear dynamical systems in state-vector form. Different gradient-based nonlinear programming methods are discussed in a unified framework and their applicability to maximum likelihood estimation is examined. The problems due to singular Hessian or singular information matrix that are common in practice are discussed in detail and methods for their solution are proposed. New results on the calculation of state sensitivity functions via reduced order models are given. Several methods for speeding convergence and reducing computation time are also discussed.

Gupta, N. K.↗

Discrete-time model reduction in limited frequency ranges

A mathematical formulation for model reduction of discrete time systems such that the reduced order model represents the system in a particular frequency range is discussed. The algorithm transforms the full order system into balanced coordinates using frequency weighted discrete controllability and observability grammians. In this form a criterion is derived to guide truncation of states based on their contribution to the frequency range of interest. Minimization of the criterion is accomplished without need for numerical optimization. Balancing requires the computation of discrete frequency weighted grammians. Close form solutions for the computation of frequency weighted grammians are developed. Numerical examples are discussed to demonstrate the algorithm.

Horta, Lucas G.↗

General modeling methods

The benefits of structural dynamics modeling methods in aerospace structures are reviewed. Four major issues in structural dynamics modeling are discussed which encompass most of its subdisciplines: reduced order modeling, constraints in problems with large motions, computational strategies, and fundamental methods. Directions for future research in these areas are addressed.

Pinson, Larry D.↗