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

Unraveling the Dynamics of Nucleosome Arrays

The organization of genomic DNA into chromatin is a fundamental determinant of genome stability, regulation, and cellular function. Nucleosomes, the basic repeating units of chromatin, assemble into higher-order structures whose organization and heterogeneity remain difficult to characterize using conventional ensemble-averaged techniques. A key need in the field is the development of experimental approaches capable of directly visualizing nucleosome assemblies and their structural variability at the single-molecule level. This LDRD Lab-Wide project focused on establishing and evaluating atomic force microscopy (AFM)–based approaches for the characterization of nucleosome assemblies. The work emphasized experimental workflows for preparing, imaging, and assessing multi-nucleosome systems, rather than isolated single nucleosomes. Through method development and exploratory measurements, the project demonstrated the feasibility of applying scanning probe microscopy to investigate chromatin-relevant assemblies and provided preliminary insight into the strengths and limitations of this approach for future quantitative studies. Results and lessons learned from this effort were disseminated to the broader scientific community through multiple national conference presentations, helping to position LLNL for continued work in chromatin and genome organization research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Sequential Bayesian Parameter Estimation of Stochastic Dynamic Load Models

In this paper we focus on the parameter estimation of dynamic load models with stochastic terms-in particular, load models where protection settings are uncertain, such as in aggregated air conditioning units. We show how the uncertainty in the aggregated protection characteristics can be formulated as a stochastic differential equation with process noise. We cast the parameter inversion within a Bayesian parameter estimation framework, and we present methods to include process noise. We demonstrate the benefits of considering stochasticity in the parameter estimation and the risks of ignoring it.

Bayesian Statistics↗
DOE OSTI2025

Comprehensive assessment of deep reinforcement learning approaches for economic dispatch in nuclear-driven microgrids

As the electrical grid integrates more variable renewable energy sources such as wind and solar, the demand for distributed and flexible systems to address this increased variability becomes critical. Nuclear-driven microgrids provide a promising solution by offering stable generation to complement intermittent renewables, ensuring grid reliability and operating efficiency. This paper proposes a recurrent deep reinforcement learning framework for optimal economic dispatch in a nuclear-powered microgrid integrating renewable energy sources, small modular reactors, battery storage systems, and balance-of-plant dynamics. A three-agent control architecture is developed, where demand and renewable energy agents act as forecasters, and a reinforcement learning-based dispatch agent performs real-time energy allocation. A nonlinear programming formulation is first used to generate an optimal baseline for benchmarking. The proposed dispatch controller, based on Proximal Policy Optimization enhanced with Long Short-Term Memory networks, exploits temporal correlations in system dynamics by taking advantage of the time series used as inputs to improve policy robustness under uncertainty. Comparative analysis against established deep reinforcement learning methods, including Proximal Policy Optimization with a feedforward architecture, Soft Actor-Critic, and Twin Delayed Deep Deterministic Policy Gradient, demonstrates superior performance. Numerical results indicate that the proposed controller achieves a 0.39% cost reduction relative to the nonlinear programming benchmark and outperforms other learning-based methods by generating additional revenue of up to 0.35%. All reinforcement learning controllers compute dispatch actions in less than 0.3 s, resulting in a computational speedup of more than three orders of magnitude over the nonlinear programming baseline. The findings of this paper highlight their applicability for real-time operation and control in nuclear-integrated microgrids under volatile operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Atomistic origin of the entropy of melting from inelastic neutron scattering and machine learned molecular dynamics

The latent heat, L, is central to melting, but its atomic origin remains elusive. It is proportional to the entropy of fusion, ΔS fus = L/T m (T m is the melting temperature), which depends on changes of atom configurations, atom vibrations, and thermal electron excitations. Here, we combine inelastic neutron scattering and machine-learned molecular dynamics to separate ΔS fus into these components for Ge, Si, Bi, Sn, Pb, and Li. When the vibrational entropy of melting, ΔS vib , is zero, ΔS fus ≃ 1.2 kB per atom. This result provides a baseline for ΔSconfig and nearly coincides with “Richard’s Rule” of melting. The ΔSfus deviates from this value for most elements, however, and we show that this deviation originates with extra ΔS vib and extra ΔS config . These two components are correlated for positive and negative deviations from Richard’s rule – the extra ΔS config is consistently ~ 80% of ΔS vib . Our results, interpreted with potential energy landscape theory, imply a correlation between the change in the number of basins and the change in the inverse of their curvature for the melting of pure elements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Butterfly valve performance factors using the multiphysics object oriented simulation environment

Butterfly valves are typically used in nuclear reactors to control incompressible fluid flow with high inlet velocities. Performance factors for butterfly valves include the pressure drop across the valve and the loss coefficient from which hydrodynamic torque and flow coefficients can be computed. This work explores a computational fluid dynamics approach for butterfly valve performance factors using the open-source Multiphysics Object Oriented Simulation Environment (MOOSE) framework. While MOOSE is often used in the nuclear energy modeling and simulation community for simulations ranging from fuel characterization to heat pipe simulation, this work employs the MOOSE open-source Navier–Stokes solver capability for simulating butterfly valve performance factors and compares those to experimentally measured results from the Advanced Test Reactor at Idaho National Laboratory at Reynolds numbers in the order of 10 6 for the partially opened configuration. The MOOSE framework results are compared against experimentally measured butterfly valve performance factors across five valve opening angles using meshes with order 10 4 – 10 5 elements. This validation serves to enable MOOSE-based multiphysics simulations incorporating the open-source Navier–Stokes module.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Decarbonization Scenarios in the United States: Comparing Biofuels Growth in Two Models - GCAM and BSM

Scenarios for deep decarbonization rely on biomass for biofuels, biopower, and bioproducts, often including negative emissions via carbon capture and storage or utilization. Despite the prominence of biomass in many deep decarbonization pathways, critical questions remain about biomass allocation, effects of transportation electrification, the pace of growth, and implications for agriculture and land use. We address these questions through a unique comparison of carbon pricing effects on the growth of biomass utilization and its effects on land use in the United States by comparing results from a multisectoral integrated assessment model, the Global Change Analysis Model [GCAM], with results from a biomass-to-biofuels system dynamics model, the Biomass Scenario Model [BSM]. We contribute to model comparison efforts by analyzing the biomass deployment needed for a scenario consistent with a "Middle of the Road" Shared Socioeconomic Pathway [SSP2] and a representative concentration pathway of 2.6 W/m2. The GCAM scenarios solve for global equilibrium conditions that are consistent with this pathway, including demands for biomass across all economic sectors and representing bioenergy with carbon capture and storage as a technology option. The BSM scenarios assess those biomass and biofuel results for the United States and identify challenges associated with that pace and amount of expansion. In the scenario analysis, we harmonize key factors such as carbon price trajectory, domestic ethanol fuel demand, ethanol blending, and arable land availability, and vary them in both models. In GCAM, we vary the carbon price, transportation electrification, ethanol blending constraints, and arable land availability inputs and the value of the carbon in land; in BSM, in addition to directly inputting certain GCAM results, we vary the maximum rate of biorefinery construction, flexibility of feedstock types across conversion processes, and policy incentives such as tax credits and renewable identification number payments. The selected carbon price trajectory results in a rapid increase in biofuel production in the United States, reaching about 9.4 EJ/year in 2060 in the highest scenario analyzed in GCAM. Results differ between the two models in timing and ultimate quantity of biomass and biofuel production. GCAM biofuel quantities generally exceed BSM amounts because CCS is applied to biofuel pathways in GCAM, and because of differences in capacity expansion and related dynamics of land allocation, biomass production, and price dynamics. These dynamics include rapid biorefinery capacity expansion in high demand cases. To satisfy this biomass demand, GCAM rapidly equilibrates land allocation, but the BSM limits the rate at which this re-allocation can occur. A further contrast with the equilibrium approach in GCAM is that the BSM represents a delay between planting and harvesting woody biomass resources. As a result of these model contrasts, feedstock costs in BSM increase more than in GCAM, and the absence of CCS in the BSM also reduces the relative economic attractiveness of biofuels production. The bottlenecks, lags, and price increases also lead to potential for volatility in feedstock price and land allocation to biomass in the BSM. GCAM has more biomass production than BSM in all scenarios, partly because of the broader, economy-wide coverage of GCAM, in contrast to BSM's exclusive focus on biofuels. In both models, trends like those of biofuels production were observed for biomass production: minimal growth without a carbon price and policy incentives, and increases with a carbon price, particularly with carbon capture and storage, because the inputs assume that biopower and biofuels decrease greenhouse gas emissions. In high policy scenarios, biomass demand is high, and the consequent high biomass prices due to the land re-allocation bottleneck in the BSM limit biofuel production even if the biorefinery capacity is expanded. However, because biomass prices do not increase as much in the low policy scenario, growth is slower and the land-reallocation bottleneck no longer dominates, such that the effect of increased capacity can be seen. Across both the models, a change in assumptions from less to more land availability increases biofuel production in both GCAM and BSM, as the upward pressure on feedstock price and volatility are both reduced.

biofuels↗
DOE OSTI2024

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2026

Out-of-time ordered correlation functions for the localized 𝑓 electrons in the Falicov-Kimball model

We provide an exact evaluation of the out-of-time correlation (OTOC) functions for the localized 𝑓-particle states in the Falicov-Kimball model within dynamical mean-field theory. Different regimes of quantum chaos and quantum scrambling are distinguished by the winding numbers of the block Toeplitz matrices used in the calculation. The similarities of these fermionic OTOCs and their logarithmic derivatives for time evolution with the OTOCs for quantum spin models with disorder are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Mechanism of the Non‐Kasha Fluorescence in Pyrene

The high-energy shoulder in the gas-phase fluorescence emission spectrum of pyrene is a well-known example of non-Kasha emission. Here, we comparatively assess two approaches, vibronic perturbation theory and nonadiabatic dynamics, in their ability to predict and explain the gas-phase fluorescence spectrum of pyrene. While both methods qualitatively capture the non-Kasha emission, they differ in their computational requirements, accuracy, and physical interpretation. Vibronic perturbation theory and nonadiabatic dynamics are complementary and can be combined in a two-step approach to non-Kasha fluorescence.

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

Filling the Gaps: A Bayesian Mixture Model for Imputing Missing Soil Water Content Data

ABSTRACT Soil water content (SWC) data are central to evaluating how soil moisture varies over time and space and influences critical plant and ecosystem functions, especially in water‐limited drylands. However, sensors that record SWC at high frequencies often malfunction, leading to incomplete timeseries and limiting our understanding of dryland ecosystem dynamics. We developed an analytical approach to impute missing SWC data, which we tested at six eddy flux tower sites along an elevation gradient in the southwestern United States. We impute missing data as a mixture of linearly interpolated SWC between the observed endpoints of a missing data gap and SWC simulated by an ecosystem water balance model (SOILWAT2). Within a Bayesian framework, we allowed the relative utility (mixture weight) of each component (linearly interpolated vs. SOILWAT2) to vary by depth, site and gap characteristics. We explored “fixed” weights versus “dynamic” weights that vary as a function of cumulative precipitation, average temperature, and time since the start of the gap. Both models estimated missing SWC data well ( R 2 = 0.70–0.88 vs. 0.75–0.91 for fixed vs. dynamic weights, respectively), but the utility of linearly interpolated versus SOILWAT2 values depended on site and depth. SOILWAT2 was more useful for more arid sites, shallower depths, longer and warmer gaps and gaps that received greater precipitation. Overall, the mixture model reliably gap‐fills SWC, while lending insight into processes governing SWC dynamics. This approach to impute missing data could be adapted to accommodate more than two mixture components and other types of environmental timeseries.

Ogle, Kiona [School of Informatics, Computing, and↗
DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗
DOE OSTI2026

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗
DOE OSTI2023

Recurrent Convolutional Deep Neural Networks for Modeling Time-Resolved Wildfire Spread Behavior

The increasing incidence and severity of wildfires underscores the necessity of accurately predicting their behavior. While high-fidelity models derived from first principles offer physical accuracy, they are too computationally expensive for use in real-time fire response. Low-fidelity models sacrifice some physical accuracy and generalizability via the integration of empirical measurements, but enable real-time simulations for operational use in fire response. Machine learning techniques have demonstrated the ability to bridge these objectives by learning first-principles physics while achieving computational speedups. While deep learning approaches have demonstrated the ability to predict wildfire propagation over large time periods, time-resolved fire-spread predictions are needed for active fire management. Here, in this work, we evaluate the ability of deep learning approaches in accurately modeling the time-resolved dynamics of wildfires. We use an autoregressive process in which a convolutional recurrent deep learning model makes predictions that propagate a wildfire over 15 min increments. We apply the model to four simulated datasets of increasing complexity, containing both field fires with homogeneous fuel distribution as well as real-world topologies sampled from the California region of the United States. We show that even after 100 autoregressive predictions representing more than 24 h of simulated fire spread, the resulting models generate stable and realistic propagation dynamics, achieving a Jaccard score between 0.89 and 0.94 when predicting the resulting fire scar. The inference time of the deep learning models are examined and compared, and directions for future work are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗
DOE OSTI2025

Theory of x-ray photon correlation spectroscopy for multiscale flows

Complex multiscale flows associated with instabilities and turbulence are commonly induced under high-energy density (HED) conditions, but accurate measurement of their transport properties has been challenging. X-ray photon correlation spectroscopy (XPCS) with coherent x-ray sources can, in principle, probe material dynamics to infer transport properties using time autocorrelation of density fluctuations. Here we develop a theoretical framework for utilizing XPCS to study material diffusivity in multiscale flows. We extend single-scale shear flow theories to broadband flows using a multiscale analysis that captures shear and diffusion dynamics. Our theory is validated with simulated XPCS for Brownian particles advected in multiscale flows. We demonstrate the versatility of the method over several orders of magnitude in timescale using sequential-pulse XPCS, single-pulse x-ray speckle visibility spectroscopy (XSVS), and double-pulse XSVS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

An Overview of the Circular Economy Lifecycle Assessment and VIsualization (CELAVI) Framework

A circular economy emphasizes the efficient use of all resources (e.g., materials, land, water). Despite anticipated overall benefits to society, the transition to a circular economy is likely to create regional differences in impacts. Current tools are unable to fully evaluate these potential externalities, which will be important for informing research prioritization and regional decision making. This presentation will review and contrast nine current impact assessment methods and describe their granularity, scope, data requirements, and capabilities within the context of identifying regional and sectoral transfers of impacts that could result from transitioning to a circular economy for energy systems. The review of current approaches highlights the need for a new hybrid circularity assessment framework that leverages multiple methods. We describe one such approach: the Circular Economy Lifecycle Assessment and VIsualization (CELAVI) framework. CELAVI uses system dynamics to model material flows for multiple circular economy pathways, network theory to track the spatial and sectoral flow of functional units across a graph, and discrete event simulation to to step through time and evaluate lifecycle assessment data at each time step. The framework is designed to be flexible and scalable enough to accommodate multiple energy materials and multiple energy technologies. CELAVI's hybridization of multiple existing methods yields additional capabilities, which may help answer questions about how flows between the technosphere and ecosphere may evolve if the circularity of energy systems is modified.

circular economy↗