Search NASA⌕ Search

SEARCH · Search NASA

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 937 records · Page 52

Optimizing DER Participation in Inertial and Primary-Frequency Response

This paper develops an approach to enable the optimal participation of distributed energy resources (DERs) in inertial and primary-frequency response alongside conventional synchronous generators. Leveraging a reduced-order model description of frequency dynamics, DERs' synthetic inertias and droop coefficients are designed to meet time-domain performance objectives of frequency overshoot and steady-state regulation. Furthermore, an optimization-based method centered around classical economic dispatch is developed to ensure that DERs share the power injections for inertial- and primary-frequency response in proportion to their power ratings. Simulations for a modified New England test-case system composed of ten synchronous generators and six instances of the IEEE 37-node test feeder with frequency-responsive DERs validate the design strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Probabilistic Power Flow Based on a Gaussian Process Emulator

In this letter, a novel Gaussian process emulator is proposed, for the first time, to conduct the probabilistic power-flow calculation. Based on Bayesian inference, a Gaussian process emulator is trained and served as a nonparametric, reduced-order model of the nonlinear power-flow model. This emulator has allowed us to evaluate the time-consuming power-flow solver at the sampled values with a negligible computational cost. The simulations reveal the excellent performance of this method.

42 ENGINEERING↗

Passive Suppression Mechanisms in Laminar Vortex-Induced Vibration of a Sprung Cylinder With a Strongly Nonlinear, Dissipative Oscillator

We study cross-flow vortex-induced vibration (VIV) of a linearly sprung circular cylinder equipped with a dissipative oscillator with cubic stiffness nonlinearity, restrained to move in the direction of travel of the cylinder. The dissipative, essentially nonlinear coupling between the cylinder and the oscillator allows for targeted energy transfer (TET) from the former to the latter, whereby the oscillator acts as a nonlinear energy sink (NES) capable of passively suppressing cylinder oscillations. For fixed values of the Reynolds number (Re = 48, slightly above the fixed-cylinder Hopf bifurcation), cylinder-to-fluid density ratio, and dimensionless cylinder spring constant, spectral-element simulations of the Navier–Stokes equations coupled to the rigid-body motion show that different combinations of NES parameters lead to different long-time attractors of the dynamics. We identify four such attractors which do not coexist at any given point in the parameter space, three of which lead to at least partial VIV suppression. We construct a reduced-order model (ROM) of the fluid–structure interaction (FSI) based on a wake oscillator to analytically study those four mechanisms seen in the high-fidelity simulations and determine their respective regions of existence in the parameter space. Asymptotic analysis of the ROM relies on complexification-averaging (CX-A) and slow–fast partition of the transient dynamics and predicts the existence of complete and partial VIV-suppression mechanisms, relaxation cycles, and Hopf and Shilnikov bifurcations. These outcomes are confirmed by numerical integration of the ROM and comparisons with spectral-element simulations of the full system.

Blanchard, Antoine↗

libROM

libROM is a collection of C++ classes that compute reduced order models and hyperreduced order models for systems of ordinary differential equations. libROM includes parallel, adaptive methods for proper orthogonal decomposition, and parallel, non-adaptive methods for hyperreduction using the discrete empirical interpolation method.

Choi, Youngsoo↗

Pressio

Reduced-order Models (ROMs). ROMs are "compressed" approximations of high-fidelity computational models. Implementing ROM methods in production codes is typically intrusive, i.e. the native codes must be adapted and modified, with individual ROM methods reimplemented in every application code. This clearly constitutes an obstacle to users who want to use ROMs. To address this, we have developed Pressio, which is a C++ header-only ROM library that can be used by an application code in a non-intrusive fashion. The main objectives behind the design of Pressio are: (a) support for general computational models (data structures, PDE types); (b) heterogeneous hardware (CPU, GPU); (c) performance; and (d) extensibility. Using metaprogramming techniques, Pressio detects and leverages the application?s native linear algebra data structures (e.g. vector and matrix) to run the ROM methods. The application only needs to expose the velocity vector and Jacobian matrix via an adapter class. Pressio checks at compile time that the target application meets the correct API. The library is structured in several modular packages including (a) non-linear least-square solvers (e.g. Gauss-Newton with and without line-search); (b) singular value decomposition (for the computation of the ROM generalized coordinates) and (c) algebraic hyper-reduction.SAND2019-5878 M

Blonigan, Patrick↗

Pressio v.1.0

Reduced-order Models (ROMs). ROMs are "compressed" approximations of high-fidelity computational models. Implementing ROM methods in production codes is typically intrusive, i.e. the native codes must be adapted and modified, with individual ROM methods reimplemented in every application code. This clearly constitutes an obstacle to users who want to use ROMs. To address this, we have developed Pressio, which is a C++ header-only ROM library that can be used by an application code in a non-intrusive fashion. The main objectives behind the design of Pressio are: (a) support for general computational models (data structures, PDE types); (b) heterogeneous hardware (CPU, GPU); (c) performance; and (d) extensibility. Using metaprogramming techniques, Pressio detects and leverages the application?s native linear algebra data structures (e.g. vector and matrix) to run the ROM methods. The application only needs to expose the velocity vector and Jacobian matrix via an adapter class. Pressio checks at compile time that the target application meets the correct API. The library is structured in several modular packages including (a) non-linear least-square solvers (e.g. Gauss-Newton with and without line-search); (b) singular value decomposition (for the computation of the ROM generalized coordinates) and (c) algebraic hyper-reduction.SAND2019-5878 M

Blonigan, Patrick↗

tesuract v.1.0

SAND2021-14370 O Tesuract (tensor surrogate approximation and computation) provides Python tools to build machine learning regressors that include polynomial chaos expansions. The software also offers methods for creation regressors for tensor outputs such as spatially varying fields and using a reduced order modeling methodology. The library is built on and compatible with the scikit-learn API, which allows integration with one of the most ubiquitous machine learning Python libraries out there. 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.

Chowdhary, Kamaljit↗

python binder for libROM

pylibROM introduces python binder for libROM through pybind11. Through python interface, users who are familiar with python can take advantage of capabilities available in libROM that is fully parallelized C++ library for reduced order modeling. libROM itself is an open source code developed at LLNL. libROM is available at https://github.com/LLNL/libROM .

Choi, Youngsoo↗

GP-BayesOpInf

SAND2025-01851O GP-BayesOpInf is a software tool that uses algorithms to combine Gaussian process regression, principal component analysis, and linear Bayesian inference to produce a probabilistic reduced-order model for time-dependent systems. Numerical examples include the compressible Euler equations for an ideal gas, a heat diffusion process with a nonlinear reaction term, and a set of ordinary differential equations describing a compartmental model in epidemiology. 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↗

Proposal and application of ROM-Lasso method for sensitivity coefficient evaluation

We propose a novel method for evaluating sensitivity coefficients of neutronics parameters to cross sections, so-called the reduced-order modeling technique ROM-Lasso. In this method, cross sections of interest are randomly sampled, and corresponding perturbed core analyses are performed. Then, the sensitivity coefficient vector of the higher-level model is expanded via the active subspace bases obtained with the lower-level model whose dimensional complexity is smaller than that of the higher-level model, and the expansion coefficients are estimated by the Lasso regression. A unique feature of the ROM-Lasso method allows the use of different bases optimized for each neutronics parameter. We conducted a verification calculation for an accelerator-driven system and demonstrated that the ROM-Lasso method can reproduce the sensitivity coefficients with a much smaller number of forward calculations than the direct method. The proposed method can be used to practically evaluate sensitivity coefficients. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Applying Constrained Bayesian Optimization to the Design of Critical Experiments

Often when planning a criticality experiment, many design configurations are iteratively investigated with a Monte Carlo transport code. The goal is that the experiment will be optimal with respect to some variable, like the fraction of fissions occurring at a certain energy range, while simultaneously being critical. Unfortunately, the Monte Carlo transport simulations are expensive, which can ultimately limit the number of configurations that can be explored. In this work, we present how Gaussian processes (GPs) can be used as a reduced-order model in a constrained Bayesian optimization (CBO) algorithm to design a criticality experiment. The GPs replace the Monte Carlo transport simulations that explore the design space. The CBO algorithm efficiently identifies new points in the design space to run the Monte Carlo transport code while respecting the criticality constraint. It does so in a manner that both improves the accuracy of the GP and finds the approximate global optimum. We demonstrate the performance of CBO with the design of a Thermal Epithermal eXperiment (TEX) for the criticality safety validation of nuclear waste models of the Hanford Tank Farm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-Fidelity Building Emulator

This dataset provides high-fidelity time series data for an emulated commercial office building sited in the Chicago, IL area during a Typical Meteorological Year (TMY). This dataset consists of air-side HVAC measurements and control inputs, and it includes normal operations as well as various implemented faults (with associated ground truth measurements) implemented on selected days. This data could be used to quantify and compare the impacts of different faults, and it could also be used as training or validation data for machine learning algorithms (e.g., reduced-order modelling, fault detection and diagnosis).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Rapid Analyses of Sparse Seismoacoustic Data Reveals the Timing and Size of the Accurate Energetic Systems Explosion

On 10 October 2025 an explosion occurred at a facility operated by Accurate Energetic Systems in Humphreys County, Tennessee. The incident resulted in 16 fatalities and created a debris field over several square kilometers. To address remaining questions about explosion timing and size, we collected about 20 seismic and 19 acoustic records of the blast from sensors up to hundreds of kilometers away. We then deployed 10 distinct physics-based, reduced order models (ROMs) that used validated geological structure and atmospheric conditions from the time of the event, along with observations of body- and surface-wave energy, as well as acoustic overpressure and phase duration. Each ROM predicted either timing, yield estimates, or both. We binned these estimates and their uncertainties according to each ROMs’ assumptions about confinement (aboveground, buried fully coupled, and buried partially coupled) and combined these estimates with other forensic data to conclude that the event occurred as a single, aboveground explosion on 10 December 2025 12:47:50.8 ±0.1 s with a yield equivalent to 11.8 [2.3,16.5] tons of Trinitrotoluene. Our estimates align with the Bureau of Alcohol, Tobacco, Firearms and Explosives inventory reports of 11–13 tons. This multimethod approach demonstrates the use of remotely observed geophysical data to rapidly aid conventional forensic investigations of accidental explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Advanced Computational Methods for Thermal Radiative Heat Transfer

Participating media radiation (PMR) in weapon safety calculations for abnormal thermal environments are too costly to do routinely. This cost may be s ubstantially reduced by applying reduced order modeling (ROM) techniques. The application of ROM to PMR is a new and unique approach for this class of problems. This approach was investigated by the authors and shown to provide significant reductions in the computational expense associated with typical PMR simulations. Once this technology is migrated into production heat transfer analysis codes this capability will enable the routine use of PMR heat transfer in higher - fidelity simulations of weapon resp onse in fire environments.

97 MATHEMATICS AND COMPUTING↗

Local Variability of the Peierls Barrier of Screw Dislocations in Ta-10W.

It is well know that the addition of substitutional elements changes the mechanical behavior of metals, a effect referred to solid solution hardening. For body-centered-cubic (BCC) metals, screw dislocation play a key role in the mechanical properties. Here the detailed modification of the Peierls barrier for screw dislocation motion in Ta with W substitutional atoms is computing using density functional theory (DFT). A reduced order model (ROM) of the influence of W substitution on the Peierls barrier is developed. The mean field change in the Peierls barrier for a Ta10W alloy is determined and shown to be larger than anticipated based on simple elasticity considerations. The ROM could be used in future calculations to determine the local variability of the Peierls barrier and the resultant influence on the motion of screw dislocation in this alloy.

42 ENGINEERING↗

Fluid-Structure Interaction of Deforming Porous Media

Sandia National Laboratories (SNL) and the Procter and Gamble Company (P&G) collaborated to develop a computational tool that predicts the in-use performance of consumer good products based on specific loading, geometries, materials, and boundary conditions. This computational tool is particularly useful in exploring the consumer-product interface and product performance. The benefit to P&G was many-fold, as the numerous project-tasks that were funded under this agreement led to advanced analysis software that allowed P&G to improve the performance and manufacturing of many of their products across several business units. The benefit to SNL was many-fold as well, including improved reduced-order models for thin-region dynamics (used for customers through DSA, NW, and EC PMUs), improved computational mechanics capability for capillary hydrodynamics, porous flow, structural mechanics, and more. Example applications of these improved capabilities at SNL include, but are not limited to, encapsulation, welding, soldering, and brazing models for NW, interface sliding models for DoD SPP programs, and porous-elastic subsidence modeling for energy and climate applications. Moreover, many current code developers on the ASC program’s SIERRA applications gained valuable skills working on this CRADA.

97 MATHEMATICS AND COMPUTING↗

Universal Design Curves for Elastomeric Direct Ink-Writing Architectures

Direct ink writing of silicone elastomers enables printing with precise control of porosity and mechan- ical properties of ordered cellular solids, suitable for shock absorption and stress mitigation applications. With the ability to manipulate structure and feedstock stiffness, the design space becomes challenging to parse to obtain a solution producing a desired mechanical response. Here we derive an analytical design approach for a specific architecture. Results from finite element simulations and quasi-static mechanical tests of two different parallel strand architectures were analyzed to understand the structure-property relationships under compressive strain. Combining effective stiffness-density scaling with least squares optimization of the stress responses yielded general response curves parameterized by resin modulus and strand spacing. An analytical expression of these curves serves as a reduced order model, which, when optimized, provides a rapid design capability for filament-based 3D printed structures. As a demonstration, the optimal design of a face-centered tetragonal architecture is computed that satisfies prescribed minimum and maximum load constraints.

36 MATERIALS SCIENCE↗

Assessment of the Area of Review and Leakage Impact for Site 7 using the NRAP-IAM-CCS Tool, Northern Michigan Basin-CarbonSAFE Phase 1 Pre-Feasibility Study

U.S. Environmental Protection Agency’s (EPA’s) Class VI regulations require owners or operators of carbon storage projects to determine an Area of Review (AoR) representative of project risk to underground sources of drinking water (USDWs). The AoR is an estimate of the project footprint and is used to develop monitoring plans to ensure protection of USDWs. In this syudy, the NRAP-IAM-CS software tool was applied to estimate the AoR and the leakage potential of legacy wells located within the AoR to impact groundwater quality at a carbon storage screening site for the NMB-CS, Phase 1 project. The NRAP-IAM-CS is a science-based toolset developed by the U.S. Department of Energy for quantitative risk assessment of geologic sequestration of carbon dioxide (CO 2 ) (Pawar et al., 2016). The toolset adopts a stochastic approach in which predictions include site uncertainties using storage reservoir, leakage scenario, and shallow groundwater impact reduced order models (ROMs). Risk-based analysis done using the NRAP-IAM-CS yielded an AoR that was comparable to estimates defined by the critical pressure needed move fluid from the reservoir to the overlying USDW through an open wellbore. The risk-based AoR was slightly smaller than that based on the critical pressure (234 km 2 compared to 269 km 2 ), because small fluxes did not impact groundwater quality. Leakage from two legacy wells located within the AoR should not impact groundwater quality over the 30-year injection period. Legacy Well 1 penetrates the simulated CO 2 plume and would require a permeability of 5 x 10 -12 m 2 (~5 Darcy) to impact groundwater quality after about 20 years of injection. Legacy Well 2 falls outside of the CO 2 plume footprint, where reservoir pressures are too small to generate large enough leakage flux to change groundwater quality even with well permeabilities as high as 5 x 10 -11 m 2 (~50 Darcy). This work represents one of the first applications of the NRAP toolset for the screening of potential CO 2 storage sites. The toolset provides a risk-based method of evaluating the AoR and the impact of CO 2 or brine leakage through legacy wells. The following recommendations will strengthen the use of probabilistic assessments for site selection and permitting of Class VI CO 2 injection wells.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗