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 163 records · Page 9

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↗

Scalar QED as a toy model for higher-order effects in classical gravitational scattering

Quantum Electrodynamics (QED) serves as a useful toy model for classical observables in gravitational two-body systems with reduced complexity due to the linearity of QED. We investigate scattering observables in scalar QED at the sixth order in the charges (two-loop order) in a classical regime analogous to the post-Minkowskian expansion in General Relativity. We employ modern scattering amplitude tools and extract classical observables by both eikonal methods and the formalism of Kosower, Maybee, and O’Connell (KMOC). In addition, we provide a simplified approach to extracting the radial action beyond the conservative sector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling Methods for Capturing System Interactions of Combined Technologies: A Study of PV + Battery

The costs of solar photovoltaics (PV) have been dropping in recent years, leading to increasing installations of solar PV systems and growing interest in how the technology will impact the electric grid at higher penetrations. Additionally, as battery costs decline it becomes important to understand the benefits and limitations of battery storage for the grid as well as potential benefits of co-locating battery storage and PV, particularly within the realm of future system planning. However, it is non-trivial to capture these potential benefits within a linearized capacity expansion model (CEM). This paper presents methodological developments to more fully represent the value and limitations of coupled PV and battery systems (PV + Battery) in CEMs using the Resource Planning Model (RPM), which co-optimizes capacity investments, transmission investments, and reduced-order dispatch in the Western Interconnection of North America through 2045. We use the model to simulate the evolution of the generation and transmission system under two core scenarios - a baseline scenario and a high renewable penetration scenario - coupled with sensitivities assuming low and midline PV and battery cost projections. When incorporating PV + battery, we find that it is important for CEMs to capture the ability of the coupled technology to provide firm capacity and reduce expected curtailment, compared to independent systems. These interactions can have even more dramatic impacts at higher solar penetrations.

14 SOLAR ENERGY↗

Energetically consistent model reduction for metriplectic systems

The metriplectic formalism is useful for describing complete dynamical systems which conserve energy and produce entropy. This creates challenges for model reduction, as the elimination of high-frequency information will generally not preserve the metriplectic structure which governs long-term stability of the system. Based on proper orthogonal decomposition, a provably convergent metriplectic reduced-order model is formulated which is guaranteed to maintain the algebraic structure necessary for energy conservation and entropy formation. Further, numerical results on benchmark problems show that the proposed method is remarkably stable, leading to improved accuracy over long time scales at a moderate increase in cost over naive methods.

42 ENGINEERING↗

Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations

This work presents and compares efficient implementations of high-order discontinuous Galerkin methods: a modal matrix-free discontinuous Galerkin (DG) method, a hybridizable discontinuous Galerkin (HDG) method, and a primal formulation of HDG, applied to the implicit solution of unsteady compressible flows. The matrix-free implementation allows for a reduction of the memory footprint of the solver when dealing with implicit time-accurate discretizations. HDG reduces the number of globally-coupled degrees of freedom relative to DG, at high order, by statically condensing element-interior degrees of freedom from the system in favor of face unknowns. The primal formulation further reduces the element-interior degrees of freedom by eliminating the gradient as a separate unknown. This paper introduces a p-multigrid preconditioner implementation for these discretizations and presents results for various flow problems. Benefits of the p-multigrid strategy relative to simpler, less expensive, preconditioners are observed for stiff systems, such as those arising from low-Mach number flows at high-order approximation. The p-multigrid preconditioner also shows excellent scalability for parallel computations. Finally, additional savings in both speed and memory occur with a matrix-free/reduced version of the preconditioner.

97 MATHEMATICS AND COMPUTING↗

A comparison of neural network architectures for data-driven reduced-order modeling

The popularity of deep convolutional autoencoders (CAEs) has engendered new and effective reduced-order models (ROMs) for the simulation of large-scale dynamical systems. Despite this, it is still unknown whether deep CAEs provide superior performance over established linear techniques or other network-based methods in all modeling scenarios. To elucidate this, the effect of autoencoder architecture on its associated ROM is studied through the comparison of deep CAEs against two alternatives: a simple fully connected autoencoder, and a novel graph convolutional autoencoder. Through benchmark experiments, it is shown that the superior autoencoder architecture for a given ROM application is highly dependent on the size of the latent space and the structure of the snapshot data, with the proposed architecture demonstrating benefits on data with irregular connectivity when the latent space is sufficiently large.

42 ENGINEERING↗

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗

Magnetic pair distribution function and half polarized neutron powder diffraction at the HB-2A powder diffractometer

Local magnetic order and anisotropy are often central for understanding fundamental behavior and emergent functional properties in quantum materials and beyond. Advances in neutron powder diffraction experiments and analysis tools now allow for quantitative determination. Here, we demonstrate this here with complementary total neutron scattering and polarized neutron measurements on the HB-2A neutron powder diffractometer at the High Flux Isotope Reactor (HFIR). In recent years, magnetic pair distribution function (mPDF) analysis has emerged as a powerful technique for probing local magnetic spin ordering of magnetic materials. This method can be broadly applied to any magnetic material but is particularly effective for studying systems with short-range magnetic order, such as materials with reduced dimensionality, geometrically frustrated magnets, thermoelectrics, multiferroics, and correlated paramagnets. Magnetic anisotropy often underpins the short-range order adopted. Half-polarized neutron powder diffraction (pNPD) can be used to determine the local susceptibility tensor on the magnetic sites to quantify the magnetic anisotropy. Combining the techniques of mPDF and pNPD can therefore provide valuable insights into local magnetic behavior. A series of measurements optimized for these techniques are presented as exemplar cases focused on frustrated materials where short-range order dominates, these include measurements to ultra-low temperature (<100 mK) not typically accessible for such experiments.

Half polarized neutron scattering↗

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗

Effect of random antiferromagnetic exchange on the spin waves in a three-dimensional Heisenberg ferromagnet

Neutron scattering is used to study spin waves in the three-dimensional Heisenberg ferromagnet YTiO 3 , with spin-spin exchange disorder introduced via La substitution at the Y site. No significant changes are observed in the spin-wave dispersion up to a La concentration of 20%. However, a strong broadening of the spectrum is found, indicative of shortened spin-wave lifetimes. Density-functional theory calculations predict minimal changes in exchange constants as a result of average structural changes due to La substitution, in agreement with the data. The absence of significant changes in the spin-wave dispersion, the considerable lifetime effect, and the reduced ordered magnetic moment previously observed in the La-substituted system are qualitatively captured by an isotropic, nearest-neighbor, three-dimensional Furthermore, Heisenberg ferromagnet model with random antiferromagnetic exchange. We therefore establish Y 1-$x$ La $x$ TiO 3 as a model system for the study of antiferromagnetic spin-exchange disorder in a three-dimensional Heisenberg ferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bias-pulsed atomic layer etching of 4H-silicon carbide producing subangstrom surface roughness

A new approach to atomic layer etching (ALE) has been demonstrated, and its application to 4H-SiC is reported here. By pulsing only the DC bias for an Ar/Cl2 inductively coupled plasma-reactive ion etching system, the etch cycle duration is reduced by more than an order of magnitude relative to conventional ALE processes. Gas flows are not changed throughout the ALE process. With this process protocol, we achieved an etch rate of 2.48±0.09 Å/cycle with 6 s cycles, an RMS surface roughness (Rq) of 0.83±0.08 Å, and an ALE synergy value of S = 99%. The parameters explored within this ALE process demonstrate effective subangstrom smoothening of 4H-SiC surfaces and is well-suited for a variety of classical and quantum device nanofabrication.

42 ENGINEERING↗

Specification Guidance for School Lighting Systems [Slides]

Lighting is a key component to creating a functional and high-quality learning environment for students and staff. LED lighting is a common lighting technology being specified today and many options are available for replacement or integration with existing fluorescent lighting systems. This document intends to help school districts retrofit or replace existing lighting systems in order to improve the classroom and school environment and reduce energy consumption and associated carbon emissions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Early-Time Exponential Instabilities in Nonchaotic Quantum Systems

The majority of classical dynamical systems are chaotic and exhibit the butterfly effect: a minute change in initial conditions has exponentially large effects later on. But this phenomenon is difficult to reconcile with quantum mechanics. One of the main goals in the field of quantum chaos is to establish a correspondence between the dynamics of classical chaotic systems and their quantum counterparts. In isolated systems in the absence of decoherence, there is such a correspondence in dynamics, but it usually persists only over a short time window, after which quantum interference washes out classical chaos. In this work, we demonstrate that quantum mechanics can also play the opposite role and generate exponential instabilities in classically nonchaotic systems within this early-time window. Our calculations employ the out-of-time-ordered correlator (OTOC)—a diagnostic that reduces to the Lyapunov exponent in the classical limit but is well defined for general quantum systems. We show that certain classically nonchaotic models, such as polygonal billiards, demonstrate a Lyapunov-like exponential growth of the OTOC at early times with Planck’s-constant-dependent rates. This behavior is sharply contrasted with the slow early-time growth of the analog of the OTOC in the systems’ classical counterparts. These results suggest that classical-to-quantum correspondence in dynamics is violated in the OTOC even before quantum interference develops.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterization and Stimulation of Selenium Reducing Microbial Consortia from Real Flue Gas Desulfurization Wastewater Biological Treatment System

Selenium is released into the environment via anthropogenic activities such as coal combustion. The U.S. Environmental Protection Agency’s Effluent Limitation Guideline (ELG) Rule limits the release of selenium species into U.S. waterways. To meet the ELG limit for selenium, coal-fired power plants must treat the wastewater produced by the flue gas desulfurization (FGD) systems. Selenium removal is of interest because small changes of the concentrations of selenium in the environment can cause detrimental ecological consequences. Biological treatment is used to remediate this wastewater, but little is known about the individual species in the microbial consortium that reduce toxic water-soluble selenium oxyanions to insoluble elemental selenium. To optimize current biological treatment capabilities of this water, the selenium-reducing capabilities must be leveraged. In order to do so, the selenium-reducing bacteria (SeRB) in biological treatment systems must be characterized. This study will be the first to characterize the microbial consortia from a real FGD WW biological treatment system and provide insights on how to better achieve selenium removal via biological treatment. This work will provide insights into the microorganisms predominantly responsible for the reduction of selenium oxyanions in FGD wastewater biological treatment systems, which can be used to better optimize more efficient biological treatment strategies of this industrial waste.

geomicrobiology↗

Optimal renormalization of multiscale systems

While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.

Price, Jacob↗