Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unit dynamics”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Status Report on the INL IES Plug-and-Play Framework

This report discusses the advancements and status of the flexible plug-and-play framework development currently ongoing that aims to integrate Modelica/Dymola with the Risk Analysis and Virtual ENvironment (RAVEN) software in terms of both Functional Mock-Up Interface (FMI)/Functional Mock-Up Unit (FMU) construction and usage, which aims to ease the sharing and simulation of complex dynamic models. This report discusses the FMI/FMU development advancements, overall focusing on the deployment of methodologies for RAVEN to export and use FMI/FMU of both Python-based models and advanced AI-constructed algorithms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Perspectives on Next Steps in Arms Control and Non-Proliferation

Much has happened in arms control over the last year as this chapter has been in development. The prospects for arms control as a viable national security tool have, if anything, grown darker. Russia has falsely equated tearing down the remaining arms control edifice with creating potential leverage vis-à-vis its actions in Ukraine. China has failed to update its talking points on being a responsible arms control actor with its massive nuclear force expansion. The United States has stated clearly that it is open for business on arms control and risk reduction, but it unfortunately has no customers in Russia and China. It remains unclear what specifically the United States is selling in terms of what it wants in a potential agreement. Competitive dynamics and an atmosphere of great power competition are driving all parties further away from the table.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Reviewing Flexibility in Industrial Electrification: U.S. Green Ammonia and Steel Industries [Slides]

The renewable energy transition in the power sector involves a paradigm shift for flexibility. Supply flexibility faces new constraints due to the increased share of variable renewable resources. Increased demand flexibility can allow less use of peaking power plants and delay need for additional capacity and transmission. Industrial customers are larger on average than residential and commercial consumers and have typically provided the largest share of demand response in the United States. We consider industrial demand, studying characteristics of flexible industrial loads. We examine the dynamics of change occurring around industrial load flexibility by focusing on two case studies: green ammonia and steel production via electric arc furnaces. Electric arc furnace steel production is an important component of current demand response programs, whereas green ammonia and green fuels offer new paradigms for flexibility. We analyze the structure and functions of the technological innovation systems of load flexibility in those two industries via interviews with twenty-two stakeholders. We conclude that the technological innovation systems are not well-functioning for flexibility in EAF based steelmaking, but are in the green ammonia space. Explicit connections between scope two greenhouse gas emissions reporting and flexibility are lacking, and industry stakeholders do not appear to make a connection between decarbonization and load flexibility.

25 ENERGY STORAGE↗
DOE OSTI2026

Modeling and Power-Hardware-in-the-Loop Validation of Synchronous Machine Governor

This paper introduces the development of a high-fidelity gas turbine governor model using a programmable logic controller for power-hardware-in-the-loop (PHIL) validation. The governor model is integrated with the National Renewable Energy Laboratory's (NREL) PHIL test bed, featuring a 2-MVA synchronous machine and a 2.5-MW variable-speed drive, to emulate NG-driven HRSGs and CTs under various operational scenarios. The primary objective of this research is to study the grid-connected and islanding operations of conventional generation sources, with representative startup sequences including turbine purge, ignition, speed ramp-up, synchronization, and breaker closure. Preliminary results of the generator governor model on NREL PHIL platform, particularly using the 2.5-MW dynamometer system, offered significant insights into the modeling techniques, hardware integration, scaling, and real-world simulation dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

ENSO Impact on Winter Precipitation in the Southeast United States through a Synoptic Climate Approach

The ENSO impact on winter precipitation in the Southeast United States was analyzed from the perspective of daily weather types (WTs). We calculated the dynamic contribution associated with the change in frequency of the WTs and the thermodynamic contribution due to changes in the spatial patterns of the environmental fields of the WTs. Six WTs were obtained using a k-means clustering analysis of 850 hPa winds in reanalysis data from November to February of 1948–2022. All the WTs can only persist for a few days. The most frequent winter weather type is WT1 (shallow trough in Eastern U.S.), which can persist or likely transfer to WT4 (Mississippi River Valley ridge). WT1 becomes less frequent in El Niño years, while the frequency of WT4 does not change much. WTs 2–6 correspond to a loop of eastward propagating waves with troughs and ridges in the mid-latitude westerlies. Three WTs with a deep trough in the Southeast U.S., which are WT2 (east coast trough), WT3 (off east coast trough) and WT6 (plains trough), become more frequent in El Niño years. The more frequent deep troughs (WTs 2, 3 and 6) and less frequent shallow trough (WT1) result in above-normal precipitation in the coastal Southeast U.S. in the winter of El Niño years. WT5 (off coast Carolina High), with maximum precipitation extending from Mississippi Valley to the Great Lakes, becomes less frequent in El Niño years, which corresponds to the below-normal precipitation from the Great Lakes to Upper Mississippi and Ohio River Valley in El Niño years, and vice versa in La Niña years. The relative contribution of the thermodynamic and dynamic contribution is location dependent. On the east coast, the two contributions are similar in magnitude.

Qian, Jian-Hua↗
DOE OSTI2026

Transient Catalytic Reaction Analysis Through Signal Defragmentation

The Temporal Analysis of Products (TAP) pulse response technique provides valuable insights into catalytic function and reaction kinetics. However, complex fragmentation patterns in the TAP mass spectrometry signals can complicate precise quantification, particularly when analyzing transient gas flux data typical of TAP experiments. This work demonstrates a standard defragmentation method that deconvolves transient TAP signals while maintaining the temporal resolution of the experiment. First, the integrals of calibration gas fluxes are used to determine the fingerprint fragmentation pattern and construct a fragmentation matrix. This matrix is then used to defragment experimental flux data at each recorded time point via a non-negative least squares regression. The effectiveness of this method is demonstrated using virtual data and control experiments with a TAP reactor system. The defragmentation is then applied to the more complex propane dehydrogenation reaction on a chromia/alumina catalyst, which can contain up to ten significant gas species in the reactor outlet. Initial propane pulsing reveals an induction period during which propane is fully oxidized to CO2, followed by partial reduction to CO. Afterwards, there is a transition in chemistries towards coking and propylene production. Our example illustrates a practical method for the accurate determination of the time-dependent reactant/product concentrations and rates for a thorough analysis of the propane dehydrogenation kinetics. This approach can be broadly applied to any transient mass spectrometry experiment for a better understanding of catalyst-reaction dynamics.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Utilization of an Advanced Sensor network to determine fuel heating value and Real-Time net unit heat rate during transient operation

Coal-fired utility boilers are being increasingly used as variable electricity generation to resolve the imbalance in the energy market from the expansion of intermittent renewable energy. The frequent transient operation required to meet residual energy demand has created a challenge for coal-fired units to operate efficiently. This work utilizes an advanced sensor network (ASN) to calculate net unit heat rate (NUHR) of a coal-fired boiler in real time through combustion calculations and statistical correlations to provide the tools for optimizing dynamic operation. Real-time heating values that were necessary to determine fuel input energy to calculate accurate NUHR were found using both fundamental and data-driven methods. Real-time NUHR shows distinct shifts that reflect changes in process conditions that will improve the ability to optimize transient operation. Data-driven heating value correlations had 24% lower root mean square error (RMSE) than the fundamental combustion calculation approach when compared to daily retrospective proximate analysis. Furthermore, the data-driven method RMSE improved by 7% with the inclusion of ASN data. Future work is to validate by comparing unit performance with and without the inclusion of NUHR as a control parameter for the dynamic neural network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Understanding the Influence of Chain Architecture on the Transport Quantities of Polymer Electrolytes with Covalently Bonded Anions

Here, we use a combination of experiments and coarse-grained molecular dynamics simulations to elucidate the structure–property relationships in polymer electrolytes obtained by the copolymerization of poly(vinyl ethylene carbonate─lithium styrene bis(trifluoromethanesulfonyl)imide) or p(VEC-LiSTFSI). Experiments show that the conductivity reduces with increasing anion (i.e., STFSI) fraction on the chain, and the cation transference number (t + ) is found to be dependent on the anion fraction. Furthermore, a significant fraction of unpolymerized VEC monomers are observed. Since it is inherently difficult to experimentally control the chain architecture and the amount of unpolymerized VEC in these systems, we perform coarse-grained molecular dynamics simulations on model polymer systems with different chain architectures to mimic the plausible experimental systems. Specifically, we look at the differences in transference numbers arising from (i) a random copolymer of VEC and STFSI monomers; (ii) a blend of VEC-STFSI copolymer with VEC monomers; and (iii) a ternary blend of the VEC homopolymer, STFSI homopolymers, and VEC monomers. The ternary blend model demonstrates the closest resemblance with the experimental transference numbers and diffusivities. The lithium diffusivity obtained from the coarse-grained models with VEC monomers (plasticizers) is about 1.5 times that of the model without VEC monomers, showing that the plasticizing effect of VEC monomers is modest. We rationalize the experimental observations based on aggregate and cluster analyses obtained from molecular simulations. This work reveals that polymer electrolyte chain architecture and plasticizers can critically influence the transport properties, and these parameters should be considered when designing single ion conducting polymeric electrolytes.

cluster distribution↗
DOE OSTI2025

Nodal hybridization in a two-dimensional heavy-fermion material

Metals with partially filled core atomic shells can form quasiparticles at a low temperature arising from the hybridization of the core level and conduction electrons. The thermodynamic and spectroscopic properties of these metals can be understood as those of a simple metal, but with a significant mass enhancement over the free electron mass—commonly referred to as heavy fermions. In most heavy-fermion materials, the hybridization is approximately isotropic in position and momentum space. However, a combination of low dimensionality and symmetry properties of the core-level wavefunctions can give rise to highly anisotropic electronic interactions with the conduction electrons. Here, in this study, we demonstrate anisotropic hybridization that vanishes along specific directions in momentum space—referred to as nodes—in a lanthanide-based two-dimensional van der Waals heavy-fermion compound, CeSiI. Quasiparticle interference measurements reveal a set of discrete hotspots with high spectral intensity on the Fermi surface. Theoretical modelling and comparison with the quasiparticle interference pattern of the non-heavy-fermion isostructural analogue LaSiI suggest that these features arise from an unconventional electron interaction involving hybridization nodes unique to CeSiI. As a result, the effective mass of the quasiparticles varies by orders of magnitude depending on their direction in momentum space.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes↗
DOE OSTI2024

Observational evidence for groundwater influence on crop yields in the United States

As climate change shifts crop exposure to dry and wet extremes, a better understanding of factors governing crop response is needed. Recent studies identified shallow groundwater—groundwater within or near the crop rooting zone—as influential, yet existing evidence is largely based on theoretical crop model simulations, indirect or static groundwater data, or small-scale field studies. Here, we use observational satellite yield data and dynamic water table simulations from 1999 to 2018 to provide field-scale evidence for shallow groundwater effects on maize yields across the United States Corn Belt. We identify three lines of evidence supporting groundwater influence: 1) crop model simulations better match observed yields after improvements in groundwater representation; 2) machine learning analysis of observed yields and modeled groundwater levels reveals a subsidy zone between 1.1 and 2.5 m depths, with yield penalties at shallower depths and no effect at deeper depths; and 3) locations with groundwater typically in the subsidy zone display higher yield stability across time. We estimate an average 3.4% yield increase when groundwater levels are at optimum depth, and this effect roughly doubles in dry conditions. Groundwater yield subsidies occur ~35% of years on average across locations, with 75% of the region benefitting in at least 10% of years. Overall, we estimate that groundwater-yield interactions had a net monetary contribution of approximately $10 billion from 1999 to 2018. This study provides empirical evidence for region-wide groundwater yield impacts and further underlines the need for better quantification of groundwater levels and their dynamic responses to short- and long-term weather conditions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗
DOE OSTI2026

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗
DOE OSTI2026

Linking Community‐Climate Disequilibrium to Ecosystem Function

Turnover in species composition often lags behind the pace of climate change, resulting in mismatches between climate and communities. However, the impact of these community‐climate disequilibria on ecosystem functions is rarely considered, and current methods for measuring disequilibria assume that species ranges were, until recently, in equilibrium with climate. Here, in this work, we develop a simple theoretical model to address both of these problems by linking community‐climate disequilibrium with ecosystem functioning. We show how disequilibrium can impair functioning in the near‐term even when climate change is expected to enhance functioning in the long‐term. Responses are most likely to change over time in communities where turnover is slow, the impact of disequilibrium counteracts the direct effects of climate on ecosystem function, and pre‐existing disequilibrium is large. These findings emphasise the importance of precise and unbiased estimates of community‐climate disequilibria for improving ecological forecasts. By fitting our model to time series of both climate and ecosystem function from a metacommunity simulation, we show the potential for community‐climate disequilibrium to be inferred without direct knowledge about species' distributions or climatic tolerances. We end by outlining a research agenda to apply dynamic disequilibrium concepts and test novel hypotheses across diverse ecosystems.

climate change↗
DOE OSTI2026

Bounded-Confidence Models of Multidimensional Opinions with Topic-Weighted Discordance

People’s opinions on a wide range of topics often evolve over time through their interactions with others. Models of opinion dynamics primarily focus on one-dimensional opinions, which represent opinions on one topic. However, opinions on various topics are rarely isolated; instead, they can be interdependent and correlated. In a bounded-confidence model (BCM) of opinion dynamics, agents are receptive to each other only if their opinions are sufficiently similar. Here, we extend classical agent-based BCMs—namely, the Hegselmann–Krause BCM, which has synchronous interactions, and the Deffuant–Weisbuch BCM, which has asynchronous interactions—to a multidimensional setting, in which the opinions are multidimensional vectors representing opinions of different topics and opinions on different topics are interdependent. To measure opinion differences between agents, we introduce topic-weighted discordance functions that account for opinion differences in all topics. We define regions of receptiveness for our models, and we use them to characterize the steady-state opinion clusters and provide an analytical approach to compute these regions. In addition, we numerically simulate our models on various networks with initial opinions drawn from a variety of distributions. When initial opinions are correlated across different topics, our topic-weighted BCMs yield significantly different results in both transient and steady states compared to baseline models, where the dynamics of each opinion topic are independent.

Mathematics and Computing↗
DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗