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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multitiered computational methodology for extracting three-dimensional rotational diffusion coefficients from x-ray photon correlation spectroscopy data without structural information

X-ray photon correlation spectroscopy (XPCS) is a powerful technique for analyzing particle systems by investigating their dynamics in suspensions across a broad range of temporal and spatial scales. This is done by illuminating samples with coherent x-ray beams and calculating the correlation function of the obtained x-ray scattering images. XPCS is uniquely suited for studying Brownian dynamics, consisting of translational and rotational diffusion. While traditional XPCS image analysis techniques can extract translational diffusion components, they are unable to estimate rotational diffusion coefficients. Here, we introduce a methodology that combines the angular-temporal cross-correlation analysis and a algorithmic framework called Multi-Tiered Estimation for Correlation Spectroscopy in 3D for estimating three-dimensional rotational diffusion coefficients from XPCS images of three-dimensional particle systems. We demonstrate our methodology for extracting rotational diffusion coefficients from XPCS data by applying it to simulated noisy x-ray images of systems of crossing nanotubes and proteins that evolve under translational and rotational Brownian motion for different diffusion rates. Furthermore, our results show that our approach determines rotational diffusion coefficients within a few percent error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Autonomous alloy composition optimization using molecular dynamics guided by a large language model

Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.

Autonomy↗
DOE OSTI2025

Temperature dependence of vacancy/self-interstitial recombination volumes in copper

We use molecular dynamics to calculate rate coefficients for recombination of vacancies with self-interstitial atoms (SIAs) in Cu at temperatures from 300 to 700 K. From these results, we calculate vacancy/SIA recombination volumes and find that they decrease from around 290 (where is one atomic volume) at 300 K to 160 at 500 K and above. By counting the number of distinct pathways by which a stable SIA may migrate to a site of spontaneous recombination with a nearby vacancy, we find a lower bound estimate of 168 for the recombination volume. We furthermore rationalize its temperature dependence based on differences between the activation energies for recombination and SIA migration. Furthermore, our work sheds light on the fundamental nature of the recombination volume and provides information that may be incorporated into multiscale models of radiation response in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

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

Anomalously slow hot carrier cooling and insulator-to-metal transition in a photo-doped Mott insulator

Photo-doped Mott insulators can exhibit novel photocarrier transport and relaxation dynamics and non-equilibrium phases. However, time-resolved real-space imaging of these processes are still lacking. Here, we use scanning ultrafast electron microscopy (SUEM) to directly visualize the spatial-temporal evolution of photoexcited species in a spin-orbit assisted Mott insulator α-RuCl 3 . At low optical fluences, we observe extremely long hot photocarrier transport time over one nanosecond, almost an order of magnitude longer than any known values in conventional semiconductors. At higher optical fluences, we observe nonlinear features suggesting a photo-induced insulator-to-metal transition, which is unusual in a large-gap Mott insulator. Our results demonstrate the rich physics in a photo-doped Mott insulator that can be extracted from spatial-temporal imaging and showcase the capability of SUEM to sensitively probe photoexcitations in strongly correlated electron systems.

42 ENGINEERING↗
DOE OSTI2024

Spatiotemporal analysis of lung immune dynamics in lethal Coccidioides posadasii infection

Coccidioidomycosis, or Valley fever, is a lung disease caused by inhalation of Coccidioides fungi, prevalent in the Southwestern United States, Mexico, and parts of Central and South America. Annually, the United States reports 10,000–20,000 cases, although those numbers are expected to increase as climate change expands the fungal geographic range. While 60% of infections are asymptomatic, 40% symptomatic infections are often misdiagnosed due to similarities with bronchitis or pneumonia. A small subset of infection progress to severe illness, necessitating a better understanding of immune responses during lethal infection. Using single-cell RNA sequencing and spatial transcriptomics, we characterized lung responses during Coccidioides infection. We identified monocyte-derived Spp1-expressing macrophages as potential mediators of tissue remodeling and fibrosis, marked by high expression of profibrotic and proinflammatory transcripts. These macrophages showed elevated TGF-β and IL-6 signaling, pathways involved in fibrosis pathogenesis. Additionally, we observed significant neutrophil infiltration and defective lymphocyte responses, indicating severe adaptive immunity dysregulation in lethal, acute infection. These findings enhance our understanding of Coccidioides infection and suggest new therapeutic targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

An Evaluation of Dynamical Downscaling Methods Used to Project Regional Climate Change

In the past decade, dynamical downscaling using “pseudo-global-warming” (PGW) techniques has been applied frequently to project regional climate change. Such techniques generate signals by adding mean global climate model (GCM)-simulated climate change signals in temperature, moisture, and circulation to lateral and surface boundary conditions derived from reanalysis. An alternative to PGW is to downscale GCM data directly. This technique should be advantageous, especially for simulation of extremes, since it incorporates the GCM's full spectrum of changing synoptic-scale dynamics in the regional solution. Here, we test this assumption, by comparing simulations in Europe and Western North America. We find that for warming and changes in temperature extremes, PGW often produces similar results to direct downscaling in both regions. For mean and extreme precipitation changes, PGW generally also performs surprisingly well in many cases. Moisture budget analysis in the Western North America domain reveals why. Large fractions of the downscaled hydroclimate changes arise from mean changes in large-scale thermodynamics and circulation, that is, increases in temperature, moisture, and winds, included in PGW by design. The one component PGW may have difficulty with is the contribution from changes in synoptic-scale variability. When this component is large, PGW performance could be degraded. Global analysis of GCM data shows there are regions where it is large or dominant. Hence, our results provide a road map to identify, through GCM analyses, the circumstances when PGW would not be expected to accurately regionalize GCM climate signals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Life Cycle Analysis of Growing Canola for Biofuel Production in the United States

This study quantifies and compares the life cycle greenhouse gas (GHG) emissions of renewable diesel (RD), sustainable aviation fuel (SAF), and biodiesel (BD) produced from two U.S. canola production systems: 1) emerging intermediate winter canola, typically grown in double- or relay-cropping systems between the growing seasons of main crops, and 2) main canola, mostly spring canola but also including winter canola, which are grown as primary crops occupying the field for a full growing season. Using the Research and Development version of the Greenhouse gases, Regulated Emissions, and Energy use in Technologies (R&D GREET) model and the most up-to-date life cycle inventory data─field trial data for intermediate winter canola (>37,000 acres) and recent national survey data for spring canola─this life cycle analysis (LCA) estimates the direct emissions from canola cultivation and harvest, the conversion of canola into fuels, fuel transportation, and combustion. In addition, we account for market-mediated emissions associated with a scenario of 0.5 billion gallons per year of spring canola-based biofuels, including induced land use change (ILUC), induced other crop (nonfeedstock) production changes, and induced livestock production changes. For intermediate winter canola, these market-mediated effects were not modeled, as ILUC is expected to be negligible due to its integration into existing rotations, and data are currently insufficient to reliably quantify other market-mediated changes. The estimated life cycle direct emissions of RD/SAF derived from intermediate winter canola and main spring canola are about 32 and 33 g of CO2-equivalent per megajoule of fuel (g CO 2 e/MJ), respectively. Corresponding emissions for BD from intermediate winter canola and main spring canola are about 30 and 31 g of CO 2 e/MJ, respectively. Farming is the dominant emissions source for both canola systems, with intermediate winter canola and main spring canola emitting about 19 and 20 g of CO 2 e/MJ, respectively. ILUC and other induced changes increase emissions of main spring canola-derived RD/SAF and BD by about 18 and 17 g of CO 2 e/MJ, respectively. These results indicate that the GHG emissions of biofuels produced from the two canola systems may differ substantially due to the different land use dynamics of the systems.

biodiesel↗
DOE OSTI2022

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Feasibility and Prospects of the High Resolution Rapid Refresh Model For Dynamic Line Rating

The ampacity of transmission lines is defined as the maximum amount of current the conductor can safely carry. It is necessary for transmission line operators to apply ampacity limits due to the thermal properties of the conductor. Dynamic Line Rating (DLR) is a technology and technique that uses the environmental conditions or a set of the conditions to calculate the ampacity of the conductor. The way the DLR is calculated has depended upon some amount of physical technology to implement the solution. The High Resolution Rapid Refresh (HRRR) Model is an atmospheric model that may provide enough resolution to perform dynamic line rating without any hardware at all. The objective is to use this model to perform Dynamic Line Rating. This will solve two of the most pervasive problems in Dynamic Line Rating, how to instrument a transmission line hundreds of miles long and then how to predict near future ampacity sufficient for transmission operators to effectively use, This technique will use the Idaho National Laboratory site as the test location for analysis. The High Resolution Rapid Refresh model is a reliable, freely available, and well-maintained weather forecasting model for the entire United States of America. The model provides enough special resolution that additional local weather monitoring devices at location are not necessary to conduct dynamic line rating. The HRRR ratings are lower on average than the locally measured ratings which allows for an additional margin of safety. The adoption of using purely digital methods for dynamic line rating could allow for wide scale adoption of the technology with very little overhead expenditure of implementation.

17 WIND ENERGY↗
DOE OSTI2021

On the Development of Order and Interfaces during the Growth of Ultrathin La 2 CuO 4 Films by Molecular Beam Epitaxy

While the atomic structure of interfaces in complex oxide heterostructures created by epitaxial growth has been investigated extensively, few studies have been conducted on how interfaces form and restructure at the initial stage of film growth. The dynamic aspects of the growth behavior can strongly influence the final interfacial atomic structure, which may lead to the emergence of interface-specific properties, such as the rise of interfacial superconductivity between certain Ruddlesden-Popper oxide materials. In this report, the structural development of La2CuO 4 thin films grown by molecular beam epitaxy on a LaSrAlO 4 substrate is investigated by X-ray diffraction measurements with rapid scans over a volume of reciprocal space. This method provides far more detailed information on interface formation than traditional fixed-point measurements. The results show that the atomic structure of the interface becomes fully established after just a single unit cell of growth. Interestingly, restructuring continues to occur within the topmost half to one unit cell of the film during the deposition process. However, diffraction intensity oscillations from both reflection high-energy electron and X-ray measurements stabilize only after the growth of two unit cells, indicating that the growth front morphology continues to evolve until the start of the third unit cell. This multimodal investigation provides insights into the atomic processes taking place during layered oxide interface formation, including the dynamical rearrangement of LaO and CuO 2 layers. Such information is not only relevant to the engineering and optimization of functional layer structures but can also be critical for ultrathin films.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impacts of A‐Site Composition on the Cation Dissolution‐Mediated Surface Restructuring of Layered Nickelate Oxide Electrocatalysts During Alkaline Oxygen Evolution Reaction

The oxygen evolution reaction (OER) is a key anodic counter‐reaction for electrochemical production of fuels and chemicals. It is hindered by sluggish four‐electron transfer kinetics requiring highly oxidative operating potentials to achieve commercially relevant rates. NiFeO x H y electrocatalysts are among the most promising for OER in alkaline electrolytes. The Ni(OH) 2 /NiOOH redox couple has been reported as the active phase in Ni‐based electrocatalysts; however, its activity is often hindered by deactivation arising from the formation of OER‐inactive insulating species. Limited strategies exist for mitigating this deactivation. This study aims to address this by interrogating the evolution of OER active sites as a function of precatalyst composition and structural properties using a series of layered, crystalline Ni‐based Ruddlesden–Popper (RP) oxides (A 2 NiO 4+δ ). In situ evolution of the active NiO x H y surface is probed through Ni‐site OER turnover frequency analysis, and electrochemical impedance spectroscopy coupled with scanning transmission electron microscopy. We show that the stability of layered nickelate oxide electrocatalysts is governed by the dynamic competition between cation dissolution and Ni‐site reversibility, which can be tuned through the A‐site composition of RP oxides. These findings yield insights toward engineering OER oxide precatalysts that optimize the stability of in situ‐generated OER active phases.

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

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2020

A Predictive Model for Survival of Escherichia coli O157:H7 and Generic E. coli in Soil Amended with Untreated Animal Manure

Abstract This study aimed at developing a predictive model that captures the influences of a variety of agricultural and environmental variables and is able to predict the concentrations of enteric bacteria in soil amended with untreated Biological Soil Amendments of Animal Origin (BSAAO) under dynamic conditions. We developed and validated a Random Forest model using data from a longitudinal field study conducted in mid‐Atlantic United States investigating the survival of Escherichia coli O157:H7 and generic E. coli in soils amended with untreated dairy manure, horse manure, or poultry litter. Amendment type, days of rain since the previous sampling day, and soil moisture content were identified as the most influential agricultural and environmental variables impacting concentrations of viable E. coli O157:H7 and generic E. coli recovered from amended soils. Our model results also indicated that E. coli O157:H7 and generic E. coli declined at similar rates in amended soils under dynamic field conditions.The Random Forest model accurately predicted changes in viable E. coli concentrations over time under different agricultural and environmental conditions. Our model also accurately characterized the variability of E. coli concentration in amended soil over time by providing upper and lower prediction bound estimates. Cross‐validation results indicated that our model can be potentially generalized to other geographic regions and incorporated into a risk assessment for evaluating the risks associated with application of untreated BSAAO. Our model can be validated for other regions and predictive performance also can be enhanced when data sets from additional geographic regions become available.

Pang, Hao↗
DOE OSTI2025

Improved Representations of Land‐Atmosphere Interactions Over the Continental U.S. Through Dynamic Root Modeling

Previous studies have identified the oversimplified root system representation as a key factor leading to inaccuracies in vegetation‐atmosphere feedbacks. In this study, a dynamic root water uptake scheme in the Noah‐MP land surface model has been coupled to the Weather Research and Forecasting (WRF) model to investigate its impact on the surface hydroclimate variables and land‐atmosphere interactions. To evaluate the impact of the dynamic root, two coupled simulations were conducted, one with the dynamic root water uptake scheme (DynRt) and one with the static root water uptake scheme (StcRt), which is based on the default root representation in Noah‐MP, with slight modifications, primarily in vegetation‐related parameters. Both DynRt and StcRt simulations were conducted with a small ensemble of three members to account for variations in physical parameterizations, initial and boundary forcing and model setup. When compared with reference data sets, the DynRt simulations show improved results than the StcRt simulations, reducing biases in the simulated leaf area index, surface energy fluxes, soil moisture and precipitation. Two different mechanisms through which roots affect land‐atmosphere coupling have been identified. Over the transitional climate zone between the dry and wet climate, the dynamic root scheme affects surface climate and land‐atmosphere coupling mainly through changes in soil moisture through hydraulic redistribution by plant root system. Over the energy‐limited mesic zone, the dynamic root affects regional land‐atmosphere coupling mainly through changes in carbon allocation. This work highlights the importance of dynamic root representation in improving vegetation‐atmosphere simulations by enhancing predictions of water, energy, and carbon fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗