Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order 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 145 records · Page 8

Component-wise reduced-order model design optimization such as for lattice design optimization

Systems and methods for optimizing a lattice structure design are disclosed herein. In some embodiments, a method for optimizing a lattice structure design can include (i) modeling the lattice structure with a component-wise reduced-order model (CWROM) and (ii) optimizing the CWROM based on a selected criterion using a topology optimization algorithm for lattice design. The selected criterion can include a boundary condition and a load applied to the lattice structure. By modeling the lattice structure as a CWROM, the optimization process can be very fast while still permitting the accurate computation of physical quantities of the lattice structure.

Choi, Youngsoo↗

Evaluation of a Reduced-Order Model for IBR Fault Response Representation via OEM Blackbox Models: Preprint

Driven by the need to capture the electromagnetic transients of transmission lines, inverter switching behavior, and detailed control systems, electromagnetic transient (EMT) studies have become increasingly important in industry, such as IBR interconnection study and fault study. However, original equipment manufacturer (OEM) inverter models typically include extensive parameters and proprietary settings that are unavailable to protection engineers. This paper introduces a data-driven, reduced-order model (ROM) developed as a PSCAD library component for use in EMT-based fault studies. The ROM replicates key OEM model behaviors without requiring detailed knowledge of control design or parameterization. The accompanying Python automation scripts streamline data generation, parameter fitting, and validation. The ROM's performance is demonstrated through comparison with both IEEE 2800-compliant and non-compliant OEM models in a real-world power system. Relay responses show nearly identical results, while simulation runtime is reduced by an average of 32.8\%, highlighting the ROM's practicality for protection engineers.

14 SOLAR ENERGY↗

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↗

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↗

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↗

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↗

Bayesian discovery of optimal reduced order models from mechanistic and experimental data: A case study of Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles rely on a silicon carbide (SiC) layer as the primary structural material and barrier to metallic fission products (FPs) release. Accurate prediction of palladium (Pd) transport and penetration is therefore critical for qualifying TRISO fuels for advanced reactors. The empirical correlation for Pd penetration in BISON is derived from historical particle-fuel data, but cannot explain the large scatter in the experimental data that arises from varying experimental conditions. To aid fuel qualification, we previously developed a mechanistic reduced order model (ROM) using BISON that resolves these dependencies. Here, in this work we build on that mechanistic ROM and perform validation and quantify its uncertainty using Bayesian uncertainty quantification (UQ). calibration against a suite of in-pile and out-of-pile experiments spanning particle compositions, geometries, and operating conditions, and we benchmark it against the empirical correlation. Bayesian UQ identifies influential parameters, calibrates them to data, and yields predictive intervals. Results show that while the empirical correlation can be tuned to fit a single experiment type, it transfers poorly; the mechanistic ROM sustains accuracy with credible uncertainty across disparate conditions. This demonstrates a practical path—via Bayesian UQ applied to mechanistic ROMs—to leverage single-effect experiments for inferring in-reactor behavior and supporting TRISO fuel qualification.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Generation of a Reduced-Order Model of an Unmanned Combat Air Vehicle Using Indicial Response Functions

Within the NATO STO AVT-251 Task Group, a generic Unmanned Combat Air Vehicle (UCAV) planform was redesigned based on requirements derived from parts of the flight envelope of the defined mission. Because of the lambda-shape planform and associated flow phenomena including shocks and vortices, the aerodynamic design process relies heavily on high-fidelity predictions from computational fluid dynamics (CFD). These simulations enable accurate flight dynamics predictions, allowing for the identification of any potential performance issues early in the design process. Potentially, a complete flight dynamics model can be derived from CFD calculations provided that the underlying reduced order model (ROM) captures the underlying physics contained within the CFD results. This work investigates the creation of a nonlinear ROM using indicial response functions. The response functions are obtained using a grid motion approach that separates the effects of angle of attack and pitch rate. This approach is demonstrated using three different CFD codes from various organizations: ENSOLV at Royal Netherlands Aerospace Centre, USM3D at NASA Langley Research Center, and Kestrel at the United States Air Force Academy. ROM predictions were generated for a manoeuver resembling a low-speed pull-up. The predictions are found to be sensitive to the quality of the steady-state solutions from which step-response calculations were started, as well as the convergence of the step responses themselves. Nevertheless, the indicial response method is shown to provide aerodynamic predictions with acceptable accuracy including transient effects, and as such is a promising method for computationally efficient flight dynamics predictions.

Michel P.C. van Rooij↗

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.↗

The application of the reduced order model Kalman filter to motion estimation of degraded image sequences

Motion estimation is a field of great interest because of its many applications in areas such as robotics and image coding. The optic flow method is one such scheme which, although fairly accurate, is prone to error in the presence of noise. This thesis describes the use of the reduced order model Kalman filter (ROMKF) in reducing errors in displacement estimation due to degradation of the sequence. The implementation of filtering and motion estimation algorithms on the SUN workstation is also discussed. Results from preliminary testing were used to determine the degrees of freedom available for the ROMKF in the SUN software. The tests indicated that increasing the state to the left leads to slight improvement over the minimum state case. Therefore, the software uses the minimum model, with the option of adding states to the left only. The ROMKF was then used in conjunction with a hierarchical pel recursive motion estimation algorithm. Applying the ROMKF to the degraded displacements themselves generally yielded slight improvements in cases with noise degradation and noise plus blur. Filtering the images of the degraded sequence prior to motion estimation was less effective in these cases. Both methods performed badly in the case of blur alone, resulting in increased displacement errors. This is thought to be due in part to filter artifacts. Some improvements were obtained by varying the filter parameters when filtering the displacements directly. This result suggests that further study in varying filter parameters may lead to better results. The results of this thesis indicate that the ROMKF can play a part in reducing motion estimation errors from degraded sequences. However, more work needs to be done before the use of the ROMKF can be a practical solution.

Simpson, Elizabeth C.↗

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↗

Generating An Advanced Cross-section Library For HTGR Pebble Bed Depletion Calculations Using Reduced-Order Model Generation Techniques

For code development, Advanced Reactor Technologies - Gas Cooled Reactors Program (ART-GCR) rely on a collaboration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, but the cross sections generation and the methodology definition is part of this program area goals. Based on previous studies in FY23, the size of microscopic cross section libraries increases rapidly with the number of tabulations, requiring significant amount of memory and drastically slowing down the Griffin calculations when evaluating cross sections via the multivariate linear interpolation approach. Rising to these challenges, this work investigates constructing Reduced-order Models (ROMs) for the multi-group microscopic cross sections to accelerate the cross section evaluation in Griffin. A database of multigroup cross sections is first collected considering all possible parameters that a designer could change for optimization. Down-selection of the ROM techniques afterward shows Deep Neural Network (DNN) as the best candidate when jointly consider memory efficiency, predictive accuracy, computational cost, scalability, flexibility and ease of implementation of the algorithms in comparison to the multidimensional interpolation. This work develops a specific interface that enables the cross section predictions using pre-trained DNN models into Griffin leveraging the existing ROM capabilities. DNNs have been trained for all isotopes for use in Griffin. Preliminary Griffin testing shows that DNNs exhibit exceptional predictive accuracy and the use of DNNs provides orders of magnitude improvement in memory efficiency compared to conventional interpolation techniques. With such ROM techniques, it holds great promise to further increase the fidelity of the Pebble Bed Reactor (PBR) simulation by increasing the number of tabulations/state variables during cross section evaluation, while maintaining the computational cost affordable in Griffin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗