SEARCH · Engineering Papers
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.
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.
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
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
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.
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 .
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.
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)
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.
Numerous studies have shown that advanced control of a wave energy converter's (WEC's) power take off (PTO) can provide significant increases (on the order of 200-300%) in WEC energy absorption. Transitioning these control approaches from simplified paper studies to application in full-scale devices remains an open and extremely challenging problem will be central to creating economically competitive WECs and delivering clean renewable energy to the US electrical grid. The Advanced WEC Dynamics and Controls project is targeted on assisting WEC developers to apply novel control systems for their devices and thus achieving major increases in performance and economic viability. The success of any control strategy is based directly upon the availability of a reduced-order model with the ability to accurately capture the dynamics of the system with sufficient accuracy. A model-scale WEC was designed and fabricated for use in studies to advance the state-of-the-art in WEC controls. This test, which is the first in a series of planned tests, focused on system identification (system ID) and model validation.
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).
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.
During the first quarter of FY 2015, the following technical progress has been made toward project milestones: 1) Autodesk delivered a new research version of ASMI to PNNL. This version includes the improved 3D fiber orientation solver, and the reduced order model (ROM) for fiber length distribution using the proper orthogonal decomposition (POD) implemented in the mid-plane, dual-domain and 3D solvers. 2) Autodesk coordinated a conference paper with PNNL reporting ASMI mid-plane fiber orientation predictions compared with the measured data for two PlastiComp plaques. This paper was accepted for presentation at the 2015 Society for Plastics Engineers (SPE) ANTEC conference. 3) The University of Illinois (Prof. Tucker) assisted team members from Purdue with fiber orientation measurement techniques, including interpretation of off-axis cross sections. 4) The University of Illinois assisted Autodesk team members with software implementation of the POD approach for fiber length modeling, and with fiber orientation modeling. 5) The University of Illinois co-authored in the SPE ANTEC paper, participated with the team in discussions of plaque data and model results, and participated in the definition of go/no-go experiments and data. 6) Purdue University (Purdue) conducted fiber orientation measurements for 3 PlastiComp plaques: fast-fill 30wt% LCF/PP center-gated, fast-fill 50wt% LCF/PA66 edge-gated and fast-fill 50wt% LCF/PA66 center-gated plaques, and delivered the fiber orientation data for these plaques at the selected locations (named A, B, and C) to PNNL. However, the data for the fast-fill 50wt% LCF/PA66 edge-gated plaque exhibited unusual variations and could not be used for the model validation. Purdue will re-measure fiber orientation for this plaque. 7) Based on discussions with the University of Illinois Purdue explained the ambiguity in the measurements of the fiber orientation components. 8) PNNL discussed with team members to establish a go/no-go decision plan for the project and submitted the established plan to DOE. 9) PNNL performed ASMI mid-plane analyses for the fast-fill center-gated 30wt% LCF/PP and 50wt% LCF/PA66 plaques and compared the predicted fiber orientations with the measured data provided by Purdue at Locations A, B, and C on these plaques. 10) Based on discussions with the University of Illinois and Autodesk, PNNL proposed a procedure to adjust fiber orientation data for Location A of the center-gated plaques so that the data can be expressed and interpreted in the flow/cross-flow direction coordinate system. 11) PNNL tested the new ASMI version received from Autodesk, examined and discussed 3D fiber orientation predictions for PlastiComp plaques. 12) PlastiComp, Inc. (PlastiComp), Toyota Research Institute North America (Toyota) and Magna Exteriors and Interiors Corp. (Magna) participated in discussions with team members on the go/no-go plan and the issues related to fiber length measurements. Toyota continued the discussion with Magna on tool modification for molding the complex part in order to achieve the target fiber length in the part.
The goal of this effort was to apply four potential control analysis/design approaches to the design of distributed grid control systems to address the impact of latency and communications uncertainty with high penetrations of photovoltaic (PV) generation. The four techniques considered were: optimal fixed structure control; Nyquist stability criterion; vector Lyapunov analysis; and Hamiltonian design methods. A reduced order model of the Western Electricity Coordinating Council (WECC) developed for the Matlab Power Systems Toolbox (PST) was employed for the study, as well as representative smaller systems (e.g., a two-area, three-area, and four-area power system). Excellent results were obtained with the optimal fixed structure approach, and the methodology we developed was published in a journal article. This approach is promising because it offers a method for designing optimal control systems with the feedback signals available from Phasor Measurement Unit (PMU) data as opposed to full state feedback or the design of an observer. The Nyquist approach inherently handles time delay and incorporates performance guarantees (e.g., gain and phase margin). We developed a technique that works for moderate sized systems, but the approach does not scale well to extremely large system because of computational complexity. The vector Lyapunov approach was applied to a two area model to demonstrate the utility for modeling communications uncertainty. Application to large power systems requires a method to automatically expand/contract the state space and partition the system so that communications uncertainty can be considered. The Hamiltonian Surface Shaping and Power Flow Control (HSSPFC) design methodology was selected to investigate grid systems for energy storage requirements to support high penetration of variable or stochastic generation (such as wind and PV) and loads. This method was applied to several small system models.
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.
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.
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.
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.
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.