Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Regional inertia dynamics of U.S. interconnections: An event-based measurement approach

Power grid inertia plays a vital role in frequency stability following large disturbances, yet its distribution across the U.S. grid is highly uneven. While interconnection-wide inertia benchmarks are useful, they can mask regional variability driven by resource mix, network coupling, and geographic separation. This paper extends event-driven inertia estimation to the regional scale using field measurements from the Frequency Monitoring Network (FNET/GridEye). Starting from balancing authority and independent system operator footprints, candidate regions are refined using a composite coherency score that combines frequency-trajectory shape similarity, timing spread, and lead/lag behavior to ensure dynamic consistency. A filtered sliding difference method (FSDM) is then used to construct regional frequency trajectories, detect disturbance onset, and compute robust regional rate-of-change of frequency (RoCoF). Regional, local, and interconnection inertia are estimated by combining RoCoF with event power imbalance, and additional indicators (regional-to-system inertia ratio and inertial-support arrival time) quantify regional-to-interconnection coupling and relative regional contributions. The method is demonstrated on eleven regions across the Eastern Interconnection (EI) and the Western Electricity Coordinating Council (WECC), with the Electric Reliability Council of Texas (ERCOT) used for validation. In ERCOT, estimates compared against energy management system (EMS) values achieve a mean absolute percentage error of 17.94%. WECC exhibits consistently shorter inertial-support arrival times (0.15–0.3 s) than EI (0.7–1.1 s), highlighting contrasting coupling and disturbance-propagation behavior. Overall, the results reveal pronounced spatial heterogeneity in inertia and coupling, underscoring the value of regional monitoring for both operational decision-making and long-term system planning.

Disturbance events↗
DOE OSTI2020

Sonic horizons and causality in phase transition dynamics

A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes nonadiabatic near the critical point. In the simplest approximation, that stage can be regarded as an “impulse” where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such a “freeze-out” might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the sonic horizon: Domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent z = 1 and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with z > 1 where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with increasing transition time. In the extreme case, the dynamical exponent z can diverge such as in the Griffiths region of the random Ising chain where the localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with z < 1 , the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the nonadiabatic stage it lags behind the adiabatic evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Unexpected Scaling of Peak Acceleration for a Yielding Mass-Spring System Subjected to a Triangular Base Acceleration Pulse

The use of bounding scenarios is a common practice that greatly simplifies the design and qualification of structures. However, this approach implicitly assumes that the quantities of interest increase monotonically with the input to the structure, which is not necessarily true for nonlinear structures. This paper surveys the literature for observations of nonmonotonic behavior of nonlinear systems and finds such observations in both the earthquake engineering and applied mechanics literature. Numerical simulations of a single degree-of-freedom mass-spring system with an elastic–plastic spring subjected to a triangular base acceleration pulse are then presented, and it is shown that the relative acceleration of this system scales nonmonotonically with the input magnitude in some cases. The equation of motion for this system is solved symbolically and an approximate expression for the relative acceleration is developed, which qualitatively agrees with the nonmonotonic behavior seen in the numerical results. The nonmonotonicity is investigated and found to be a result of dynamics excited by the discontinuous derivative of the base acceleration pulse, the magnitude of which scales nonmonotonically with the input magnitude due to the fact that the first yield of the spring occurs earlier as the input magnitude is increased. The relevance of this finding within the context of defining bounding scenarios is discussed, and it is recommended that modeling be used to perform a survey of the full range of possible inputs prior to defining bounding scenarios.

42 ENGINEERING↗
DOE OSTI2025

Data-Driven Compositional Optimization in Misspecified Regimes

With a manifold growth in the scale and intricacy of systems, the challenges of parametric misspecification become pronounced. These concerns are further exacerbated in compositional settings, which emerge in problems complicated by modeling risk and robustness. In “Data-Driven Compositional Optimization in Misspecified Regimes,” the authors consider the resolution of compositional stochastic optimization problems, plagued by parametric misspecification. In considering settings where such misspecification may be resolved via a parallel learning process, the authors develop schemes that can contend with diverse forms of risk, dynamics, and nonconvexity. They provide asymptotic and rate guarantees for unaccelerated and accelerated schemes for convex, strongly convex, and nonconvex problems in a two-level regime with extensions to the multilevel setting. Surprisingly, the nonasymptotic rate guarantees show no degradation from the rate statements obtained in a correctly specified regime and the schemes achieve optimal (or near-optimal) sample complexities for general T-level strongly convex and nonconvex compositional problems.

Business & Economics↗
DOE OSTI2023

Advancing uncertainty characterization for understanding projected water scarcity in multi-sector, multi-actor river basins across scales.

Invited Talk at AGU Fall Meeting 2023 Modeling how human institutions and infrastructure interact with the water cycle is essential to better understand the vulnerability, and resilience of water resources systems from the local to the global scale. This is especially true when investigating multi-sector, multi-actor responses to the effects of long-term changes and short-term shocks. It is also well recognized that uncertainty present throughout the modeling cycle (e.g., in data, functional relations, and model coupling approaches) limits our ability to trace the interactive dynamics of human-water systems, as well as to quantify their implications for management and planning. Systems with large numbers of diverse stakeholders further compound this challenge, as uncertain drivers and complex dynamics might have very disparate effects on water users. The work presented is conducted through the Integrated Multisector Multiscale Modeling (IM3) Science Focus Area, which explores how human and natural systems co-evolve in response to change. Through this multi-year effort, we have developed exploratory modeling methods to better characterize how the uncertain human and natural drivers of water scarcity yield consequential vulnerabilities in institutionally complex, multi-sectoral systems. This talk will specifically present on a series of uncertainty characterization experiments performed in the Upper Colorado River Basin, a sub-basin of the Colorado with thousands of water users. These complementary and systematic experiments are aimed at understanding: How are various uncertain stressors (e.g., climate change, demand growth) affecting the diverse water users of this basin in terms of water shortage? What are the key drivers of this shortage for each user? Methods and results from this work are used to address additional questions on the ability of adaptation to modulate the effects of these uncertain drivers, and on their compounding effects across spatial scales and through sectoral interactions.

Hadjimichael, Antonia↗
DOE OSTI2025

Experimental analysis of gas dynamics in the reactor cavity section of a High-Temperature Gas-Cooled Reactor during an accident scenario

Here, this study investigates the behavior of a High-Temperature Gas-Cooled Reactor (HTGR) during a break of the helium pressure boundary event. Understanding the gas dynamics that occur during a break could help develop safety systems that reduce the probability of air entering the reactor core. To achieve this objective, experimental studies were conducted in a scaled-down HTGR experimental facility to provide information on the air-helium gas mixture within the confinement building and the impact of the ventilation system on the helium/air gas concentration following the depressurization event. The study evaluated different configurations, including active ventilation time scales, break sizes, and locations. The analysis of the active ventilation time revealed that a ventilation time of 22 seconds produced positive results, and smaller break sizes resulted in longer depressurization times, which improved the ventilation process. The orientation and elevation of the breaks had little effect on the oxygen concentration in the cavity but did impact the gas velocities.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Magnetic Energy Release, Plasma Dynamics, and Particle Acceleration in Relativistic Turbulent Magnetic Reconnection

In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be the primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue in studies of relativistic reconnection is the consequence of three-dimensional dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of three-dimensional, fully kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron–electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long-wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. Additionally, we find that plasma dynamics in RTMR is vastly different from its 2D counterpart in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted by the secondary kink instability. This turbulent evolution leads to superdiffusive behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and the timescale for energy release can be very fast and do not depend strongly on the turbulence amplitude. Additionally, the main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the nonideal electric field acceleration plays a subdominant role. We also discuss possible observational implications of three-dimensional RTMR in high-energy astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2021

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Budget-Constrained Sizing of Renewable and Energy Storage Systems for Farm-Scale Ammonia Production Within the Food-Energy-Water Nexus

In the transition toward sustainable agriculture, farms have emerged as eco-friendly pioneers, harnessing clean hybrid wind and solar systems to improve farm performance. A concern in this paradigm is the effective sizing of renewable energy systems to ensure optimal energy use within budget considerations. This research focuses on optimizing renewable energy sizing in small-scale ammonia production to meet specific farm demands and enhance local resilience, emphasizing the interplay between environmental and economic factors. These findings promise increased energy efficiency and sustainability in this innovative agricultural sector. Additionally, our approach considers small-scale ammonia plant needs and the dynamic relationships between ammonia, water, and farm demands. Simulations demonstrate substantial cost savings in farm electricity consumption. Specifically, scenarios with renewable energy integration in the farm can reduce at least 13% electricity cost compared to a grid-dependent system in the 15-year simulation.

ammonia↗
DOE OSTI2020

Fast Magnetic Reconnection with Turbulence in High Lundquist Number Limit

Here, we use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number S limit (above ~104), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate ~0.01–0.13 is weakly dependent on S. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a superdiffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
DOE OSTI2021

Hamiltonian reduction of Vlasov–Maxwell to a dark slow manifold

In this work, we show that non-relativistic scaling of the collisionless Vlasov–Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov–Maxwell phase space. Vlasov–Maxwell dynamics restricted to the slow manifold recovers the Vlasov–Poisson and Vlasov–Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov–Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov–Maxwell solutions that do not excite light waves, and are therefore ‘dark’. We provide a heuristic lower bound for the time interval over which Vlasov–Maxwell solutions initialized optimally near the slow manifold remain dark. We also show how the dynamics on the slow manifold naturally inherits a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov–Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evaluating ecosystem water use efficiency and recovery dynamics during flash droughts: insights from observations and model simulations

Flash droughts (FD), rapidly emerging in a warming future, disrupt ecosystems, agriculture, and water security. Ecosystem water use efficiency (WUE), the ratio of gross primary production (GPP) to actual evapotranspiration (AET), balances carbon assimilation and water loss. FD rapidly disrupts this balance, making WUE critical for assessing plant stress and recovery. Here, this study investigates the dynamics of landscape-scale WUE, and the components of GPP and AET under FD utilizing both observed data from the Missouri Ozark AmeriFlux site (US-MOz) and version 2 of the U.S. Department of Energy’s Earth, Energy, Exascale System Model (E3SM) Land Model (ELMv2). Observations and simulations reveal GPP as dominant for WUE during earlier FD events (2005, 2007, 2012), shifting to AET in recent events (2014, 2018). This agreement indicates that the ELM can capture the shifting dynamics of GPP and AET in regulating WUE under FD conditions. However, the ELM systematically underestimates both GPP and AET and does so in a manner that does not preserve their ratio. As a result, WUE is also underestimated, suggesting that GPP is more strongly underestimated than AET. Furthermore, the ELM also underestimates the speed of GPP recovery, producing an artificially prolonged GPP recovery time following FD events. Observed environmental drivers such as vapor pressure deficit (VPD), soil moisture (SM), and predawn leaf water potential (PLWP) effectively predict WUE, but ELM primarily highlights SM, underestimating VPD’s role. This study demonstrates that relying solely on soil moisture fails to capture the rapid hydraulic recovery observed in PLWP, underscoring the necessity of integrating plant hydraulics into land surface models to improve flash drought predictability.

Evapotranspiration↗
DOE OSTI2021

Accelerated annealing of colloidal crystal monolayers by means of cyclically applied electric fields

Abstract External fields are commonly applied to accelerate colloidal crystallization; however, accelerated self-assembly kinetics can negatively impact the quality of crystal structures. We show that cyclically applied electric fields can produce high quality colloidal crystals by annealing local disorder. We find that the optimal off-duration for maximum annealing is approximately one-half of the characteristic melting half lifetime of the crystalline phase. Local six-fold bond orientational order grows more rapidly than global scattering peaks, indicating that local restructuring leads global annealing. Molecular dynamics simulations of cyclically activated systems show that the ratio of optimal off-duration for maximum annealing and crystal melting time is insensitive to particle interaction details. This research provides a quantitative relationship describing how the cyclic application of fields produces high quality colloidal crystals by cycling at the fundamental time scale for local defect rearrangements; such understanding of dynamics and kinetics can be applied for reconfigurable colloidal assembly.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic Model of Back-to-Back Converter for System Level Phasor Simulation

The power system is expected to evolve rapidly with increasing deployment of power electronic interface and conditioning systems, microgrids, and hybrid AC/DC grids. Among power electronic systems, back-to-back (BTB) converters can be a powerful interface to integrate microgrids and networked microgrids. However, in system level dynamic simulations for synchronous machine dominated system with BTB converters, a detailed electromagnetic sinusoidal models is not warranted as these simulations are focused on electromechanical transients. This paper focus on developing a simple phasor model for back- to-back converters that can be easily integrated into powerflow solvers to facilitate large-scale power system simulations. The model is implemented using C++ language and integrated into GridLAB-D software as a potential new capability. The GridLAB-D phasor domain model is validated against the electromagnetic transient (EMT) simulation of the detailed switching model. Simulation results show that the phasor model successfully captures the dominant dynamics of the converters.

Mahmood, Hisham↗