SEARCH · Engineering Papers

Results for “Reduced Order Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Analysis of the Effect of Communication Latencies on HVDC-Based Damping Control

A wide-area controller to damp inter-area oscillations in the North American Western Interconnection (WI) by modulating power transfers in a HVDC link is used in this paper to investigate the effects that latencies in its feedback signals have on its performance. This controller uses two feedback measurements to perform its control action. The analysis show that the stabilizing effect of the controller in transient stability and small signal stability is compromised as the feedback measurements experience higher delays. The results show that one of the feedback signals can tolerate more delay than the other. The analysis was performed with Bode plots and time domain simulations on a reduced order model of the WI from which a linear version was obtained.

Wilches-Bernal, Felipe↗

DOE OSTI2022

Provably-Stable Overload Ride-Through Control for Grid-Forming Inverters Using System-Wide Lyapunov Function Analysis

A key challenge associated with a grid-forming (GFM) inverter based resource (IBR) is its behavior during severe grid disturbances: since a GFM inverter regulates voltage in the fast timescale instead of current or power, it may experience a transient overload of current, power and/or energy during a severe grid disturbance. While many promising control strategies for overload ride-through have been proposed over the past two decades, transient stability of the system during and after the transition to an overload ride-through control mode remains difficult to guarantee. In this work, a novel overload ride-through control strategy is proposed for a system of grid-forming inverters that takes both self-protection and system-wide transient stability into account. A proposed system-level supervisory control uses slow communication to pre-emptively assign a set of local ride through control parameters to individual GFM IBR, including a current-limiting virtual reactance, that guarantees that synchronism is still preserved for any set of anticipated grid disturbances. At the core of the supervisory control lies a Lyapunov-function-based routine capable of establishing a strong, albeit conservative, transient stability guarantee for the system. Here, the proposed overload ride-through control strategy is validated via numerical integration of a reduced-order model, as well as through detailed electromagnetic transient (EMT) simulation.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2018

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↗

DOE OSTI2020

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↗

DOE OSTI2017

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↗

DOE OSTI2019

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↗

DOE OSTI2019

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↗

DOE OSTI2019

Pressio v.1.0

Blonigan, Patrick↗

DOE OSTI2021

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↗

DOE OSTI2023

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↗

DOE OSTI2024

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↗

DOE OSTI2022

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↗

DOE OSTI2020

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↗

DOE OSTI2016

Advanced WEC Dynamics and Controls, Test 1

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.

16 TIDAL AND WAVE POWER↗

DOE OSTI2026

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↗

DOE OSTI2026

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↗

DOE OSTI2015

Predictive Engineering Tools for Injection-Molded Long-Carbon-Fiber Thermoplastic Composites - FY 2015 First Quarterly Report

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.

42 ENGINEERING↗

DOE OSTI2015

Final Technical Report: Distributed Controls for High Penetrations of Renewables

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.

14 SOLAR ENERGY↗