Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic framework”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interactive multiscale modeling to bridge atomic properties and electrochemical performance in Li-CO 2 battery design

Li-CO 2 batteries are promising energy storage systems due to their high theoretical energy density and CO 2 fixation capability, relying on reversible Li 2 CO 3 /C formation during discharge/charge cycles. Here, we present a multiscale modeling framework integrating Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties. The considered Li-CO 2 battery consists of a lithium metal anode, an ionic liquid electrolyte, and a carbon cloth cathode with Sb 0.67 Bi 1.33 Te 3 catalyst. DFT and AIMD determined the electrical conductivities of Sb 0.67 Bi 1.33 Te 3 and Li 2 CO 3 using the Kubo–Greenwood formalism and studied the CO 2 reduction mechanism on the cathode catalyst. MD simulations calculated the CO 2 diffusion coefficient, Li + transference number, ionic conductivity, and Li + solvation structure. The FEA model, parameterized with atomistic simulation data, reproduced the available experimental voltage–capacity profile at 1 mA/cm 2 and revealed spatio-temporal variations in Li 2 CO 3 /C deposition, porosity, and CO 2 concentration dependence on discharge rates in the cathode. Accordingly, Li 2 CO 3 can form large and thin film deposits, leading to dispersed and local porosity changes at 0.1 mA/cm 2 and 1 mA/cm 2 , respectively. The capacity decreases exponentially from 81,570 mAh/g at 0.1 mA/cm 2 to 6200 mAh/g at 1 mA/cm 2 , due to pore clogging from excessive discharge product deposition that limits CO 2 transport to the cathode interior. Therefore, the performance of Li-CO 2 batteries can be improved by enhancing CO 2 transport, regulating Li 2 CO 3 deposition, and optimizing cathode architecture.

Battery performance↗
DOE OSTI2025

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

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

Modeling the impact of smoke from prescribed fire on road visibility

Prescribed fires are planned to achieve conservation and fuel reduction objectives while minimizing smoke ground concentration to limit health impacts and road visibility impairment. Prescribed burns cannot indeed be conducted if those hazards are not within predefined limits. This paper proposes a new framework to evaluate road visibility that overcomes the limitation of the state of the art model, VSMOKE. The framework leverages the fast-running framework QUIC-Fire/QUIC-SMOKE to capture fire and smoke dynamics, and the timing and duration of hazardous conditions on the road network close to the burn unit (within 50–100 km). The paper presents parametric study using a real burn plot at Fort Stewart (GA, USA), under hypothetical wind conditions to understand the interplay between buoyancy and smoke dilution. Here, results showed that faster winds caused fire escape while slower winds did not achieve a complete burn. Furthermore, faster winds featured brief road visibility reduction below braking distance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying structural errors in cloud condensation nuclei activity from reduced representation of aerosol size distributions

Aerosol effects on clouds and radiation are the dominant contribution to uncertainty in radiative forcing relative to the pre-industrial atmosphere. While previous studies have assessed the impact of parametric uncertainty on modeled forcing, structural errors from the numerical representation of particle distributions have not been well quantified. Here we present a framework for quantifying error in aerosol size distributions and cloud condensation nuclei activity, which we apply to the widely used 4-mode version of the Modal Aerosol Module (MAM4). Box model predictions from the MAM4 are evaluated against the Particle Monte Carlo Model for Simulating Aerosol Interactions and Chemistry (PartMC-MOSAIC), a benchmark model that tracks the evolution of individual particles. We show that size distributions simulated by MAM4 diverge from those simulated by PartMC-MOSAIC after only a few hours of aging by condensation and coagulation in polluted conditions, which leads to large errors in modeled cloud condensation nuclei concentrations. We find that differences between MAM4 and PartMC-MOSAIC are largest under polluted conditions, where the size distribution evolves rapidly though aging by condensation of semi-volatile substances and coagulation among particles. These findings suggest that structural error in modeled aerosol properties contributes to the large inter-model variability in aerosol radiative forcing.

Fierce, Laura M.↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Temporal and spatial separations between spin glass and short-range order

Broken-symmetry-induced order parameters account for many phenomena in physics. For spin glasses, this framework dictates the theoretical construction, whereas experiments have only established dynamical behaviors but not the physical entity. Experimental techniques have limitations when the spin glass is probed as an isolated state. Here, we create an evolution from a long-range order using well-controlled non-magnetic substitution on a sublattice. Neutron magnetic diffuse scattering of pico-second timescale reveals that the dynamics of short- and long-range order formation are not affected by disorder, but their spatial ranges are. Across all specimens, the inflection point of spin correlation length’s temperature dependence fully matches with the peak in heat capacity, while spin glass can freeze at millisecond timescale either above or below this characteristic temperature of spin-order formation. Our results identify the spin glass as single, uncorrelated spins at domain walls between spin clusters, and an uncorrelated coexistence of the two.

Dronova, Margarita G [Okinawa Institute of Science↗
DOE OSTI2020

Simulations of the Biodegradation of Citrate Based Polymers for Artificial Scaffolds Using Accelerated Reactive Molecular Dynamics

Here, we explore the reactivity and mechanical properties of poly(1,6-hexanediol-co-Citric acid) via ReaxFF molecular dynamics simulations. We implement an accelerated scheme within the ReaxFF framework to study the hydrolysis reaction of the polymer which is provided with sufficient amount of energy known as the restrain energy after a suitable pre-transition state configuration is obtained to overcome the activation energy barrier and the desired product is obtained. The validity of the ReaxFF force field is established by comparing the ReaxFF energy barriers of ester and ether hydrolysis with benchmark DFT values in literature. We perform chemical and mechanical degradation of polymer chain bundles at 300 K. We find that ester hydrolyzes faster than ether due to lower activation energy barrier of the reaction. The selectivity of the bond-boost scheme has been demonstrated by lowering the boost parameters of the accelerated simulation which almost stops the ether hydrolysis. Mechanical degradation of pre-hydrolyzed and intermittent hydrolyzed polymer bundles is performed along the longitudinal direction at two different strain rates. We find that the tensile modulus of the polymers are increases with increase in strain rates which shows that polymers show a strain dependent behavior. The tensile modulus of the Polyester-ether is higher than Polyester yet reaches yield-stress faster than Polyester. This makes Polyester more ductile than Polyester-ether.

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

Azo(xy) vs Aniline Selectivity in Catalytic Nitroarene Reduction by Intermetallics: Experiments and Simulations

Intermetallic nanoparticles are promising catalysts in hydrogenation and fuel cell technologies. Much is known about the ability of intermetallic nanoparticles to selectively reduce nitro vs alkene, alcohol, or halide functional groups; less is known about their selectivity toward aniline vs azo or azoxy condensation products that result from the reduction of a nitro group alone. Because azo(xy)arenes bear promise as dyes, chemical stabilizers, and building blocks to functional materials but can be difficult to isolate, developing high surface area nanoparticle catalysts that display azo(xy) selectivity is desirable. To address this question, we studied a family of nanocrystalline group 10 metal (Pd, Pt)- and group 14 metal (Ge, Sn, Pb)-containing intermetallics-Pd 2 Ge, Pd 2 Sn, Pd 3 Sn 2 , Pd 3 Pb, and PtSn-in the catalytic reduction of nitroarenes. In contrast to monometallic Au, Pt, and Pd nanoparticles and "random" PdxSn1 – x nanoalloys, which are selective for aniline, nanoparticles of atomically precise intermetallic Pd 2 Ge, Pd 2 Sn, Pd 3 Sn 2 , and PtSn prefer an indirect condensation pathway and have a high selectivity for the azo(xy) products. The only exception is Pd 3 Pb, the most active among the intermetallic nanoparticles studied here, which is instead selective for aniline. Employing a novel application of molecular dynamics-based on machine learned potentials within a DeePMD framework-to heterogeneous catalysis, we are able to identify key reaction species on the different types of catalysts employed, furthering our understanding of the unique selectivity of these materials. By demonstrating how intermetallic nanoparticles can be as active yet more selective than other more traditional catalysts, this work provides new physical insights and opens new opportunities in the use of these materials in other important chemical transformations and applications.

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

Large Ensemble Particle Filter for Spatial Climate Reconstructions Using a Linear Inverse Model

Proxy records that document the last 2000 years of climate provide evidence for the wide range of the natural climate variability from inter-annual to secular timescales not captured by the short window of recent direct observations. Assessing climate models ability to reproduce such natural variations is crucial to understand climate sensitivity and impacts of future climate change. Paleoclimate data assimilation (PDA) offers a powerful way to extend the short instrumental period by optimally combining the physics described by General Circulation Climate Models (GCMs) with information from available proxy records while taking into account their uncertainties. Here we present a new PDA approach based on a sequential importance resampling (SIR) Particle filter (PF) that uses Linear Inverse Modeling (LIM) as an emulator of several CMIP-class GCMs. We examine in a perfect-model framework the skill of the various LIMs to forecast the dynamics of the surface temperatures and provide spatial field reconstructions over the last millennium in a SIR PF. Our results show that the LIMs allow for skillful ensemble forecasts at 1-year lead-time based on GCMs dynamical knowledge with best prediction in the tropics and the North Atlantic. The PDA further provides a set of physically consistent spatial fields allowing robust uncertainty quantification related to climate models biases and proxy spatial sampling. Our results indicate that the LIM yields dynamical memory improving climate variability reconstructions and support the use of the LIM as a GCM-emulator in real reconstruction to propagate large ensembles of particles at low cost in SIR PF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Local and General Patterns of Terrestrial Water‐Carbon Coupling

Abstract Terrestrial carbon uptake and water availability have coupled feedbacks; specifically water uptake for plant growth and soil drying via transpiration. While we might expect this coupling over time at arid sites, climatic water availability also widely covaries geographically with biomass variables that control photosynthetic rates. Using eddy covariance data globally, we find convex, positively‐covarying relations between carbon uptake and a turbulent flux metric controlled by land surface moisture (r = 0.73 monthly across sites) at the site level. We estimate a general, empirical relationship based on site‐wise water‐carbon dynamics. Most sites, and the general relationship, show strong power‐law dependence, implicating the role of sub‐seasonal land‐cover dynamics. We also find that long‐term mean carbon/water states follow a similar convex relationship to the site‐specific temporal dynamics. We discuss opportunities and caveats for space‐for‐time frameworks of carbon/water feedback processes globally.

Geology↗
DOE OSTI2026

Spring Dust in Colorado Plateau: Transport Pathways and Interannual Variability Derived From Two Decades of MERRA‐2 Reanalysis

Spring dust storms in the U.S. Southwest significantly impact environmental and human systems, yet their climatological patterns and driving mechanisms remain poorly understood. Using two decades of MERRA-2 reanalysis data and self-organizing map clustering, we identified four distinct dust transport pathways from surrounding and remote deserts for the top 10% of spring dust events impacting the Colorado Plateau: (a) intense southwesterly near-surface transport where eastward upward transport is blocked by topography, (b) substantial remote transport plus weaker surface emissions, (c) moderate southwesterly transport with stronger upper-level transport, and (d) rare northwesterly pattern linked to the Great Basin sources, which accounts for 38%, 14%, 37%, and 11% of total events, respectively. The frequency and intensity of these events strongly correlate with large-scale climate variability, with increased dust events during La Niña. Our results provide a robust framework for improving seasonal dust forecasts and understanding future dust dynamics under climate change.

aerosol transport↗
DOE OSTI2023

Unravelling crystal growth of nanoparticles

Crystal growth from nanoscale constituents is a ubiquitous phenomenon in biology, geology and materials science. Numerous studies have focused on understanding the onset of nucleation and on producing high-quality crystals by empirically sampling constituents with different attributes and varying the growth conditions. However, the kinetics of post-nucleation growth processes, an important determinant of crystal morphology and properties, have remained underexplored due to experimental challenges associated with real-space imaging at the nanoscale. Here we report the imaging of the crystal growth of nanoparticles of different shapes using liquid-phase transmission electron microscopy, resolving both lateral and perpendicular growth of crystal layers by tracking individual nanoparticles. We observe that these nanoscale systems exhibit layer-by-layer growth, typical of atomic crystallization, as well as rough growth prevalent in colloidal systems. Surprisingly, the lateral and perpendicular growth modes can be independently controlled, resulting in two mixed crystallization modes that, until now, have received only scant attention. Combining analytical considerations with molecular dynamics and kinetic Monte Carlo simulations, we develop a comprehensive framework for our observations, which are fundamentally determined by the size and shape of the building blocks. Furthermore, these insights unify the understanding of crystal growth across four orders of magnitude in particle size and suggest novel pathways to crystal engineering.

Imaging techniques↗
DOE OSTI2025

Perfect spinfluid: A divergence-type approach

We present a new formulation of nondissipative relativistic spin hydrodynamics that incorporates spin degrees of freedom into the divergence-type theory framework. Due to the divergence-type structure, it is straightforward to enforce nonlinear causality and symmetric hyperbolicity of the equations of motion, ensuring local well-posedness of the initial-value problem and stability of the theory. Furthermore, in a specific realization based on spin kinetic theory, we prove that the equations of motion remain nonlinearly causal and symmetric-hyperbolic to all orders in the spin potential, provided a specific thermodynamic constraint is satisfied. Here, this framework can be applied for numerical simulations to study the dynamics of spin-polarized fluids, such as the quark-gluon plasma in heavy-ion collisions.

Chirality↗
DOE OSTI2020

Direct time-domain determination of electron-phonon coupling strengths in chromium

In this work, we report the results of an ultrafast, direct structural measurement of optically pumped phonons in a Cr thin film using ultrashort x-ray pulses from a free-electron laser. In addition to measuring and confirming the known long-wavelength dispersion relation of Cr along a particular acoustic branch, we are able to determine the relative phase of the phonons as they are generated. The Cr sample exhibits two generation mechanisms for the phonons: the releasing of a preexisting charge density wave at higher frequencies, and the creation of an acoustic strain pulse via laser heating that dominates at lower frequencies. Furthermore, for the latter mechanism, we are able to measure the frequency dependence of the time required to generate the phonons. To explain the observed magnitude and slope of the delays, we perform first-principles simulations in the framework of density functional perturbation theory and ab initio molecular dynamics to fit anharmonic phonon models. In conclusion, these results show that the wave-vector dependence of the electron-phonon coupling is the driving mechanism behind the delay times: Phase-space limitation leads to higher times near the zone center. The absolute magnitudes of the delay times measured are found to be much shorter than the equilibrium electron-phonon coupling times we compute, indicating that the coupling strength is greatly enhanced when the electronic system is out of equilibrium with the lattice, as has been seen in bismuth and other systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Barrier-Certificated Reinforcement Learning Approach for Enhancing Power System Transient Stability

Increasing integration of renewable resources brings more flexibility and poses new challenges to modern power systems, leading to highly nonlinear and complex dynamics. Here, this paper aims to provide a general solution framework to traditional control problems, such as frequency control and voltage control, which attempt to maintain the stability of either synchronous generators-governed or inverter-governed systems when subjected to a disturbance and simultaneously guarantee operational constraints, providing a complete complement to existing works on control design. Building on reinforcement learning (RL) and control barrier functions, the framework includes two subsystems, i.e., a model-free controller and a barrier-certification system, which discover RL-based control actions and sequentially filter them using a barrier certificate to satisfy operational constraints. Calculating a barrier function is generally challenging for a complex power system. This is addressed by representing the barrier function using neural networks (NNs) and data-based approaches. An adaptive method is introduced to certify the neural barrier function that perseveres barrier conditions, which is more compatible with online implementation. The proposed framework synthesizes a stabilizing controller that satisfies predefined safety regions. The effectiveness of the proposed framework is demonstrated via several comparative case studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

densegranFoam

densegranFoam is a dense granular material flow solver developed using the OpenFoam Computational Fluid Dynamics (CFD) tool. It is developed using Euler-Euler framework. The transport model is implemented using the shear rate dependant constitutive model. The primary focus of the solver is to simulate shear rate dependency on the dense granular flow bulk behavior.

Ahsan, Syed↗
DOE OSTI2025

LLaMP v0.1.0

Reducing hallucination of Large Language Models (LLMs) is imperative for use in the sciences, where reliability and reproducibility are crucial. However, LLMs inherently lack long-term memory, making it a nontrivial, ad hoc, and inevitably biased task to fine-tune them on domain-specific literature and data. LLaMP is a multimodal retrieval-augmented generation (RAG) framework of hierarchical reasoning and acting (ReAct) agents that can dynamically and recursively interact with Materials Project to ground large language models on high-fidelity materials informatics.

Riebesell, Janosh [Lawrence Berkeley National Labo↗