Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reduced models”

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 253 records · Page 14

ROMNet: Neural-Network-Based Surrogates for Reduced-Order Model Dynamics

SAND2023-05829O ROMNet is a toolbox built on the Tensorflow library for training reduced-order, neural operator surrogate models. The software is used for training reduced order neural operator surrogate models for dynamical systems applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗

On Ohm’s law in reduced plasma fluid models

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced-Order Models for the Aeroelastic Analysis of Ares Launch Vehicles

This document presents the development and application of unsteady aerodynamic, structural dynamic, and aeroelastic reduced-order models (ROMs) for the ascent aeroelastic analysis of the Ares I-X flight test and Ares I crew launch vehicles using the unstructured-grid, aeroelastic FUN3D computational fluid dynamics (CFD) code. The purpose of this work is to perform computationally-efficient aeroelastic response calculations that would be prohibitively expensive via computation of multiple full-order aeroelastic FUN3D solutions. These efficient aeroelastic ROM solutions provide valuable insight regarding the aeroelastic sensitivity of the vehicles to various parameters over a range of dynamic pressures.

Silva, Walter A.↗

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗

PIV Measurements of Chevrons on F400-Series Tactical Aircraft Nozzle Model

Reducing noise of tactical jet aircraft has taken on fresh urgency as core engine technologies allow higher specific-thrust engines and as society become more concerned for the health of its military workforce. Noise reduction on this application has lagged the commercial field as incentives for quieting military aircraft have not been as strong as in their civilian counterparts. And noise reduction strategies employed on civilian engines may not be directly applicable due to the differences in exhaust system architecture and mission. For instance, the noise reduction technology of chevrons, examined in this study, will need to be modified to take into account the special features of tactical aircraft nozzles. In practice, these nozzles have divergent slats that are tied to throttle position, and at take off the jet flow is highly overexpanded as the nozzle is optimized for cruise altitude rather than sea level. In simple oil flow visualization experiments conducted at the onset of the current test program flow barely stays attached at end of nozzle at takeoff conditions. This adds a new twist to the design of chevrons. Upon reaching the nozzle exit the flow shrinks inward radially, meaning that for a chevron to penetrate the flow it must extend much farther away from the baseline nozzle streamline. Another wrinkle is that with a variable divergence angle on the nozzle, the effective penetration will differ with throttle position and altitude. The final note of realism introduced in these experiments was to simulate the manner in which bypass flow is bled into the nozzle wall in real engines to cool the nozzle, which might cause very fat boundary layer at exit. These factors, along with several other issues specific to the application of chevrons to convergent-divergent nozzles have been explored with particle image velocimetry measurements and are presented in this paper.

Bridges, James↗

Analysis of Exhaust Plume Effects on Sonic Boom for a 59-Degree Wing Body Model

Reducing or eliminating the operational restrictions of supersonic aircraft over populated areas has led to extensive research at NASA. Restrictions are due to the disturbance of the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Recent work has been performed to reduce the magnitude of the sonic boom N-wave generated by airplane components with focus on shock waves caused by the exhaust nozzle plume. Previous Computational Fluid Dynamics (CFD) analyses showed how the shock wave formed at the nozzle lip interacted with the nozzle boat-tail expansion wave. The nozzle lip shock moved with increasing nozzle pressure ratio (NPR) and reduced the nozzle boat-tail expansion. Lip shock movement caused a favorable change in the observed pressure signature. These results were applied to a simplified supersonic vehicle geometry with no inlets and no tail, in which the goal was to demonstrate how under-expanded nozzle operation reduced the sonic boom signature by twelve percent. A secondary goal was to demonstrate the use of the Cart3D inviscid code for off-body pressure signatures including the nozzle plume effect.

Castner, Raymond S.↗

Stress-constrained topology optimization of lattice-like structures using component-wise reduced order models

We report lattice-like structures can provide a combination of high stiffness with light weight that is useful in many applications, but a resolved finite element mesh of such structures results in a computationally expensive discretization. This computational expense may be particularly burdensome in many-query applications, such as optimization. We develop a stress-constrained topology optimization method for lattice-like structures that uses component-wise reduced order models as a cheap surrogate, providing accurate computation of stress fields while greatly reducing run time relative to a full order model. We demonstrate the ability of our method to produce large reductions in mass while respecting a constraint on the maximum stress in a pair of test problems. The ROM methodology provides a speedup of about 150x in forward solves compared to full order static condensation and provides a relative error of less than 5% in the relaxed stress.

97 MATHEMATICS AND COMPUTING↗

Reduced-order modeling for efficient cross section library development in high-temperature gas reactor pebble-bed depletion analysis

Accurate modeling of running-in and equilibrium conditions in pebble-bed reactors (PBRs) requires precise microscopic multigroup neutron cross sections. In Griffin, deterministic neutronics calculations rely on multivariate interpolation over large cross section libraries, resulting in significant memory usage and performance bottlenecks. This work, together with a companion paper on Griffin integration, explores reduced-order models (ROMs) to replace interpolation with lightweight surrogates. Several ROM techniques are benchmarked, with deep neural networks (DNNs) demonstrating superior memory efficiency, scalability, and predictive accuracy. A total of 295 DNNs were trained to build a comprehensive isotope library, integrated into Griffin through a custom LibTorch interface for depletion analysis. Initial results demonstrate that DNN-based ROMs drastically reduce memory demands while preserving accuracy, enabling finer tabulations and additional state variables without overhead. In conclusion, the framework also supports online cross section generation and real-time DNN updates through transfer learning, improving fidelity by capturing self-shielding and evolving nuclide compositions during burnup.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation↗

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

Si–Cl 2 –Ar + Atomic Layer Etching Window: A Fundamental Study Using Molecular Dynamics Simulations and a Reduced Order Model

Silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ) is studied by using molecular dynamics (MD) simulations and a reduced order model (ROM). Here, the purpose of this study is to elucidate the properties of the ALE window, a range of ion energies where the amount of Si etched over a series of cycles is nonzero and nearly independent of ion energy. Experimental studies of the Si–Cl 2 –Ar + ALE system report contradictory results related to the ALE window’s ion energy range. Both MD simulations and the ROM show that there is an ALE window present from approximately 15 to 20 eV for normal incidence argon ions. The Si–Cl 2 –Ar + system, therefore, exhibits a narrow ALE window. The amount of Si etched per cycle is less than one atomic layer because of the higher etch yield of Cl atoms relative to atomic Si and silicon chlorides. A modified version of the ROM with an artificially increased Si physical sputtering threshold energy expands the ALE window, illustrating the importance of the difference in chemical and physical sputtering threshold energies in the ALE window energy range. The ROM is also used to examine the dependence of the EPC on the Ar + ion fluence.

energy↗

CSI compensation for reduced-order model based control of a flexible robot manipulator

In controller design for flexible structures, certain system modes are extremely important for the overall performance of the structure. A reduced-order model (ROM) based control focuses on these modes, providing a viable, active control algorithm for large systems. Unfortunately, unmodeled structure dynamics can interact with the ROM controller (CSI) and cause crippling deterioration of system performance, possibly to the point that system stability is lost. A residual model filter (RMF) eliminates one channel of control structure interaction (CSI), while adding only a simple, second-order filter to the control loop. Thus, the ROM controller can be designed independently, based strictly on performance criteria, and residual mode filters can then be selected to compensate for CSI. A flexible robot manipulator is used for preliminary experimentation with the ROM/RMF design methodology. Since the controller was to be implemented both with, and without compensation for CSI, the ROM control gains are carefully chosen such that closed loop stability is never compromised. In this way, RMF effectiveness is easily evaluated in terms of the improvement in system performance resulting from CSI compensation.

Reisenauer, Brian T.↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Advanced Cross Section Library Generation using Reduced Order Models

Deterministic neutronics calculations rely on multigroup neutron cross section libraries, which consist of databases of tabulated values, used to calculate the neutron cross sections through multivariate linear interpolation. However, interpolation of the multidimensional cross section data becomes memory inefficient and time consuming as the number of tabulations increases, significantly slowing down the neutronics calculation, especially in the case of microscopic cross section libraries where every isotope (on the order of hundreds) has its own set of specific reactions and cross sections. In order to address this challenge, this work constructs efficient and robust reduced-order models (ROMs) of the multi-group cross sections to support the Griffin simulation of high-temperature gas-cooled reactors (HTGRs). The first part of the study investigates the linearity of the multigroup cross section data across isotopes, reaction types, and energy groups on pre-generated datasets for the purpose of dimensionality reduction. Secondly, a down-selection of ROM techniques is presented on representative classical machine learning (ML) techniques, including variants of linear regression, kernel-based methods, tree-based algorithms, and artificial neural networks. The selection criteria jointly consider the memory efficiency, predictive accuracy, prediction speed, and scalability in comparison to the multidimensional interpolation. Among all the ML techniques, deep neural networks (DNNs) have proven to be the best selection with sufficient accuracy, high robustness, good memory efficiency, great scalability, and superior flexibility. DNNs have been trained for all isotopes in this work and systematic Griffin testing is ongoing to ensure the feasibility of this ROM technique for predicting cross section and reducing memory requirements without a significant sacrifice in computational performance.

42 - ENGINEERING↗