Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced order systems”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Real-Time Distributed Control of Smart Inverters for Network-level Optimization

The limitations of centralized optimization methods in managing electric power distribution systems operations have led to the distributed paradigm of computing and decision-making. Unfortunately, the existing distributed optimization algorithms are limited in their applicability to managing fast varying phenomena such as those resulting from highly variable Distributed Energy Resource (DER) generation patterns. They require a large number of communication rounds (in the order of 10 2 to 10 3 ) among the computing agents to solve one instance of the optimization problem. Related real-time distributed control methods are equally limited in their applications to power distribution systems with fast-changing DER generation; they require hundreds of rounds of communication and thus are slow in tracking the network-level optimal solutions. In this paper, we propose a novel distributed voltage controller that provides a fast-tracking of rapidly varying DER generation profiles while simultaneously converging to network-level optimal solutions within a few communication rounds. The proposed control algorithm leverages the radial topology of the system, which reduces the required communication rounds to reach the network-level optimum solution by order of magnitude. The novelty lies in carefully reducing the electrical network model from the perspective of each distributed controller and enabling appropriate data sharing among upstream and downstream nodes to achieve fast convergence. The simulation results demonstrate the effectiveness of the proposed approach in minimizing the feeder losses while maintaining the node voltage within the pre-specified limits.

voltage control, optimization, reactive power, inv↗
DOE OSTI2021

Amorphous dispersions of flubendazole in hydroxypropyl methylcellulose: Formulation stability assisted by pair distribution function analysis

Here, we use X-ray pair distribution function (PDF) analysis applied to high-energy synchrotron X-ray powder diffraction data to evaluate the amorphous solid dispersions interactions and their aging stability. The obtained systems are based on hydroxypropyl methylcellulose (hypromellose) derivatives and flubendazole (FBZ) drug dispersions prepared using a spray-dryer technique. We carry out stability studies under aging parameters (40 °C/75% relative humidity) to tune the systems’ recrystallization. The results reveal that ion-base interactions between the drug-polymer matrix are responsible for reducing clustering processes yielding slower recrystallization and different ordering in the hypromellose phthalate (HPMCP/FBZ) and hypromellose acetate succinate (HPMC-AS/FBZ) systems and complete drug clustering in hypromellose (HPMC-E3/FBZ). The structural ordering was accessed using differential X-ray PDFs that revealed the region between 3.5 Å and 5.0 Å could be related to FBZ intermolecular interactions and is more ordered for the least stable system (HPMC-E3/FBZ) and less ordered for the most stable system (HPMCP/FBZ). These results show that the ion-base interactions between drug and matrix occur at these intermolecular distances.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Comparing and contrasting In-Vial and full-scale systems for sparging volatile analytes

Here, in-vial sparging was demonstrated as an effective, practical alternative to a full-scale sparging system for supporting the analysis of volatile constituents. Using elemental mercury (Hg 0 ) and toluene as representative purgeable analytes, the mass removal for various sparge configurations was measured and a reduced order model was developed and validated. In the primary experiments, Hg 0 in the sparge gas was trapped on activated carbon or gold, thermally desorbed, and quantified using atomic absorption or atomic fluorescence spectroscopy. Toluene experiments using the same in-vial sparge apparatus and sparge parameters were performed to demonstrate the applicability of the reduced order model to a broad range of compounds. Toluene removal was tracked by measuring the remaining toluene in sparged aliquots using Ultraviolet-visible (UV–Vis) spectroscopy. For the sparging, flow rates varied from 25 to 75 mL/min for periods from 0 to 30 min. Sparge performance, mass removal as a function of time, and sparge gas volume were measured for both in-vial and full-scale systems. A model based on dimensionless Henry's Law coefficient, normalized sparge gas volume, and fractional extent of equilibrium matched the experimental data for both compounds and provides a practical tool for future applications. For the conditions tested in this study, the calibrated model indicated that the sparge gas in the in-vial system reached approximately 33% of its equilibrium value before exiting the water surface, while a full-scale system reached approximately 100%. The tests validated the quality, reproducibility, and predictability of sparging performance for both full scale and in-vial sparge systems. Related factors such as waste generation, worker risk, and labor were also assessed. Full scale sparge systems provide the advantage of lower detection levels due to larger sample volume, while the in-vial sparge systems provide advantages for most other factors; including automatability, reducing secondary wastes, lessening the need to clean and check the sparge apparatus, and lowering labor and costs. The data and associated reduced order model support continued development and deployment of in-vial sparge platforms as a practical option for analysis of purgeable analytes such as volatile organic compounds and volatile metals/organometallics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A reduced order sulfuric acid decomposition model for a nuclear-powered hybrid sulfur cycle

A reduced order model of sulfuric acid decomposition within a bayonet chemical reactor was developed to support the U. S. Department of Energy Integrated Energy System program, and address the lack of knowledge in scaling and integration for joint chemical and nuclear processes. Sulfuric acid decomposition within a bayonet reactor was modeled to provide chemical and thermodynamic data relevant to advanced nuclear reactor-driven integrated energy systems based on desired operational scale and operational conditions. The temperature range required for high-temperature advanced nuclear reactor integrated energy systems, 750-850 °C, was shown to produce reasonable agreement (within a few percent relative error) with past models and experimental data, and yielded good efficiency results for bayonet reactor operations. The results of the reduced order model agreed with previous work from Savannah River National Labs within a maximum of 3.4% error on the decomposition of sulfur trioxide, and on previous Hybrid Sulfur flowsheets from Gorensek and Summers that showed operational temperature, pressure, and composition ranges for efficiency which made the Hybrid Sulfur cycle competitive with water electrolysis. Here, the agreement with previous high-fidelity models provided a framework for future Integrated Energy System grid evaluations with an advanced nuclear reactor and large-scale hydrogen production using a mathematical model to represent chemical operations.

08 HYDROGEN↗
DOE OSTI2020

Bench Scale Testing of Next Generation Hollow Fiber Membrane Modules

This project seeks to develop a next generation membrane material (PI-2) for application in Air Liquide’s (AL) hybrid process. Combining the cold membrane operation with an integrated CO 2 Compression and Purification Unit (CPU) significantly reduces the overall cost of CO 2 capture. The proposed work will accelerate the membrane bundle fabrication, scale-up and performance validation of this new material. The proposed novel membrane development, combined with the results of the previous Cold Membrane Project, will provide the best opportunity for AL to meet the long term development targets of the DOE with respect to timing, cost, and performance for CO 2 capture. Our previous field test program (DE-FE0013163) showed promising results for the hybrid process utilizing the existing commercial AL membrane bundles (PI-1 material). The project work has also shown the potential for significant improvements through initial laboratory tests with the novel material (PI-2). The initial PI-2 results show a step-change in membrane permeance. This will enable further reduction in the cost of CO 2 capture by reducing the number of membrane modules and associated equipment in the system. In order to capture this value, however, the new material needs to be validated by field testing of large bundles, representative of commercial production. Lastly, a comprehensive evaluation of novel hybrid processes and costs needs to be completed to ensure optimal use of this improved material performance. The current project seeks to develop this next generation membrane to TRL 5 (approximately 600 - 1,000 Nm3/h flue gas, 350 - 600 scfm, 0.2 - 0.3 MWe equivalent). Air Liquide manufacturing methodology will be applied to fabricate prototype (4”) membrane modules from PI-2 material. Once the manufacturing methods prove reliable with the new material, and performance is demonstrated in line with expectation, 6” commercial scale bundles will be fabricated. Recognizing that a low cost membrane solution by itself may not meet DOE’s LCOE target for CO 2 capture, we will also perform simulation studies of novel hybrid processes in combination with the next generation membrane technology. This approach has the best potential for exceeding the DOE targets. The best evaluation of the new commercial units will be performance testing under real flue gas at the NCCC. Finally, a rigorous TEA, which incorporates the novel process studies and field performance, will be completed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Scientific Discovery with Physics-Informed System Identification (Abbreviated Report)

My fellowship research focused on making physics-based simulations faster and more useful through machine learning. Many problems in science and engineering are governed by partial differential equations, but high-fidelity simulations are often too expensive to run repeatedly. I worked on improving Latent Space Dynamics Identification (LaSDI), a reduced-order modeling framework that compresses large simulation data sets into a smaller representation and then learns how that representation evolves over time. The motivation was to develop reduced models that remain accurate for more challenging systems, especially when predictions must remain reliable over long time intervals or when the underlying dynamics are more complicated than standard methods can easily handle. I also contributed to related work on Quandary, a high-performance software effort for simulation and control of open quantum systems, before focusing primarily on Latent Space Dynamics Identification methods. The main outcomes of the fellowship were two new algorithms (both of which were published), Rollout-LaSDI and Higher-Order LaSDI, together with supporting work on multi-stage Latent Space Dynamics Identification. Rollout-LaSDI improved long-term prediction by training the model to stay accurate over extended time horizons, and Higher-Order LaSDI broadened the method so it could model systems with higher-order time dynamics. My contributions to multistage Latent Space Dynamics Identification also helped show that its later training stages could be simplified without losing effectiveness, and that this behavior held across different model architectures and training strategies. Taken together, these advances improved the accuracy, flexibility, and practical value of reduced-order modeling tools for computational science.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Nonlinear dynamics of quadratic gravity in spherical symmetry

We present the first numerically stable nonlinear evolution for the leading-order gravitational effective field theory (quadratic gravity) in the spherically-symmetric sector. The formulation relies on (i) the harmonic gauge to cast the evolution system into quasilinear form (ii) the Cartoon method to reduce to spherical symmetry in keeping with the harmonic gauge, and (iii) order reduction to first order (in time) by means of introducing auxiliary variables. The well posedness of the respective initial-value problem is numerically confirmed by evolving randomly perturbed flat-space and black-hole initial data. In conclusion, our study serves as a proof-of-principle for the possibility of stable numerical evolution in the presence of higher derivatives.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multiscale modeling high-order methods and data-driven modeling

Projection-based reduced-order models (ROMs) comprise a promising set of data-driven approaches for accelerating the simulation of high-fidelity numerical simulations. Standard projection-based ROM approaches, however, suffer from several drawbacks when applied to the complex nonlinear dynamical systems commonly encountered in science and engineering. These limitations include a lack of stability, accuracy, and sharp a posteriori error estimators. This work addresses these limitations by leveraging multiscale modeling, least-squares principles, and machine learning to develop novel reduced-order modeling approaches, along with data-driven a posteriori error estimators, for dynamical systems. Theoretical and numerical results demonstrate that the two ROM approaches developed in this work - namely the windowed least-squares method and the Adjoint Petrov - Galerkin method - yield substantial improvements over state-of-the-art approaches. Additionally, numerical results demonstrate the capability of the a posteriori error models developed in this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An effective field theory for non-maximal quantum chaos

In non-maximally quantum chaotic systems, the exponential behavior of out-of-time-ordered correlators (OTOCs) results from summing over exchanges of an infinite tower of higher “spin” operators. We construct an effective field theory (EFT) to capture these exchanges in (0 + 1) dimensions. The EFT generalizes the one for maximally chaotic systems, and reduces to it in the limit of maximal chaos. The theory predicts the general structure of OTOCs both at leading order in the 1/N expansion (N is the number of degrees of freedom), and after resuming over an infinite number of higher order 1/N corrections. These general results agree with those previously explicitly obtained in specific models. We also show that the general structure of the EFT can be extracted from the large q SYK model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Optimal Operations Management of Mobility-on-Demand Systems

The emergence of the sharing economy in urban transportation networks has enabled new fast, convenient and accessible mobility services referred to as Mobilty-on-Demand systems (e.g., Uber, Lyft, DiDi). These platforms have flourished in the last decade around the globe and face many operational challenges in order to be competitive and provide good quality of service. A crucial step in the effective operation of these systems is to reduce customers' waiting time while properly selecting the optimal fleet size and pricing policy. In this paper, we jointly tackle three operational decisions: (i) fleet size, (ii) pricing, and (iii) rebalancing, in order to maximize the platform's profit or its customers' welfare. To accomplish this, we first devise an optimization framework which gives rise to a static policy. Then, we elaborate and propose dynamic policies that are more responsive to perturbations such as unexpected increases in demand. We test this framework in a simulation environment using three case studies and leveraging traffic flow and taxi data from Eastern Massachusetts, New York City, and Chicago. Our results show that solving the problem jointly could increase profits between 1% and up to 50%, depending on the benchmark. Moreover, we observe that the proposed fleet size yield utilization of the vehicles in the fleet is around 75% compared to private vehicle utilization of 5%.

Wollenstein-Betech, Salomón↗
DOE OSTI2025

New Results on Communication- and Memory-Aware Load Balancing Model and Algorithms

While load balancing in distributed-memory computing has been well-studied, we present an innovative approach to this problem: a unified, reduced-order model that combines three key components to describe “work” in a distributed system: computation, communication, and memory. Our model enables an optimizer to explore complex tradeoffs in task placement, such as augmented parallelism, at the expense of data replication increasing memory usage. We propose a fully distributed, heuristic-based load balancing optimization algorithm, and demonstrate that it quickly finds close-to-optimal solutions. We formalize the complex optimization problem as a mixed-integer linear program, and compare it to our strategy. Finally, we show that when applied to an electromagnetics code, our approach obtains up to 2.3x speedups for the imbalanced execution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluation of pressure reconstruction techniques for Model Order Reduction in incompressible convective heat transfer

This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗
DOE OSTI2025

Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)

The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.

Density functional theory↗
DOE OSTI2025

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗
DOE OSTI2024

Probabilistic locked mode predictor in the presence of a resistive wall and finite island saturation in tokamaks

We present a framework for estimating the probability of locking to an error field in a rotating tokamak plasma. This leverages machine learning methods trained on data from a mode-locking model, including an error field, resistive magnetohydrodynamics modeling of the plasma, a resistive wall, and an external vacuum region, leading to a fifth-order ordinary differential equation (ODE) system. It is an extension of the model without a resistive wall introduced by Akçay et al. [Phys. Plasmas 28, 082106 (2021)]. Tearing mode saturation by a finite island width is also modeled. We vary three pairs of control parameters in our studies: the momentum source plus either the error field, the tearing stability index, or the island saturation term. The order parameters are the time-asymptotic values of the five ODE variables. Normalization of them reduces the system to 2D and facilitates the classification into locked (L) or unlocked (U) states, as illustrated by Akçay et al., [Phys. Plasmas 28, 082106 (2021)]. This classification splits the control space into three regions: L̂, with only L states; Û, with only U states; and a hysteresis (hysteretic) region Ĥ, with both L and U states. In regions L̂ and Û, the cubic equation of torque balance yields one real root. Region Ĥ has three roots, allowing bifurcations between the L and U states. The classification of the ODE solutions into L/U is used to estimate the locking probability, conditional on the pair of the control parameters, using a neural network. We also explore estimating the locking probability for a sparse dataset, using a transfer learning method based on a dense model dataset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗