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

Interpretable artificial intelligence and exascale molecular dynamics simulations to reveal kinetics: Applications to Alzheimer's disease

The rapid increase in computing power, especially with the integration of graphics processing units, has dramatically increased the capabilities of molecular dynamics simulations. To date, these capabilities extend from running very long simulations (tens to hundreds of microseconds) to thousands of short simulations. However, the expansive data generated in these simulations must be made interpretable not only by the investigator who performs them but also by others as well. Here, we demonstrate how integrating learning techniques, such as artificial intelligence, machine learning, and neural networks, into analysis pipelines can reveal the kinetics of Alzheimer's disease (AD) protein aggregation. Finally, we review select AD targets, describe current simulation methods, and introduce learning concepts and their application in AD, highlighting limitations and potential solutions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Polymers in Deep Eutectic Solvents

The project investigated the behavior of polymers in ionic liquids such as deep eutectic solvents using an array of techniques. These included the development of atomistic force fields, coarse graining these force fields to the united atom level, large scale molecular dynamics (MD) simulations using new thermostat algorithms, and machine learning (ML) methods for the phase behavior. The project demonstrated the feasibility and accuracy of first principles force fields for ionic liquids, deep-eutectic solvents, and urea-water mixtures. A novel hierarchical coarse graining method was then used to develop accurate and efficient united-atom models, using which microsecond simulations were performed for polymers in ionic liquids. These simulations were in quantitative agreement with experiment, thus resolving previous controversies. Methods were also developed to obtain the potential of mean force between complex ions in solution. Finally, supervised ML methods were developed for the phase behavior of polymers in ionic liquids.

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

Coarsening dynamics of Ising-nematic order in a frustrated Heisenberg antiferromagnet

We study the phase ordering dynamics of the classical antiferromagnetic 𝐽 1 −𝐽 2 (nearest-neighbor and next-nearest-neighbor couplings) Heisenberg model on the square lattice in the strong frustration regime (𝐽 2 /𝐽 1 > 1/2). While thermal fluctuations preclude any long-range magnetic order at finite temperatures, the system exhibits a long-range spin-driven nematic phase at low temperatures. The transition into the nematic phase is further shown to belong to the two-dimensional Ising universality class based on the critical exponents near the phase transition. Our large-scale stochastic Landau-Lifshitz-Gilbert simulations find a two-stage phase ordering when the system is quenched from a high-temperature paramagnetic state into the nematic phase. In the early stage, collinear alignments of spins lead to a locally saturated Ising-nematic order. Once domains of well-defined Ising order are developed, the late-stage relaxation is dominated by curvature-driven domain coarsening, as described by the Allen-Cahn equation. The characteristic size of Ising-nematic domains scales as the square root of time, similar to the kinetic Ising model described by the time-dependent Ginzburg-Landau theory. Our results confirm that the late-stage ordering kinetics of the spin-driven nematic, which is a vestigial order of the frustrated Heisenberg model, belongs to the dynamical universality class of a nonconserved Ising order. Interestingly, the system shows no violation of the superuniversality hypothesis under weak bond disorder. The dynamic scaling invariance is preserved in the presence of weak bond disorder. Here, we also discuss possible applications of our results to materials for which vestigial Ising-nematic order is realized.

Antiferromagnets↗
DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

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

Introducing Kynema, an Open-Source Performance-Portable Flexible-Multibody-Dynamics Solver

In this talk we introduce Kynema, an open-source general flexible-multibody-dynamics solver that is well suited for simulating wind turbine structural dynamics. Kynema uses a Lie-group time integrator for constrained systems and runs on both CPUs and GPUs. Timing results for simulations are presented for the IEA 15-MW turbine with and without aerodynamic forces.

17 WIND ENERGY↗
DOE OSTI2025

Assessment of ESM Readiness Level for Exascale HPC

Advancement of Earth System Models (ESMs) is becoming increasingly challenging due to a confluence of factors including increasing model complexity – to more fully represent the earth system, increasing spatial resolution - to achieve higher accuracy by resolving fine-scale dynamical to physical, biological, and chemical processes and their interaction, increasing ensemble size - to more accurately represent predictive uncertainty, and increased computing requirements – to enable more accurate and timely weather predictions and climate projections for societal benefit. The belief by many that computing will take care of itself is no longer valid given the disruptive changes in HPC that are driving up the cost of computing, increasing the difficulty of using emerging HPC effectively, and exposing limits in parallelism, portability and scalability of the ESM applications themselves.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fallowed agricultural lands dominate anthropogenic dust sources in California

Air pollution remains a major problem in many parts of California, significantly impacting public health and regional climate. However, the contribution of anthropogenic dust from agricultural sources, among major pollutants in California’s semi-arid Central Valley, remains largely unclear. Here, we used the Cropland Data Layer from the U.S. Department of Agriculture to identify crop types and land use/cover and leveraged satellite-derived estimates of major dust events between 2008 and 2022 over California. We identified fallowed land—an unplanted agricultural land parcel—as a key anthropogenic dust source in California. Specifically, we find that the Central Valley accounts for about 77% of total fallowed land areas in California, where they are associated with about 88% of major anthropogenic dust events. We also find that the geographic coverage of these fallowed lands expanded between 2008 and 2022 with associated increasing anthropogenic dust activities. Additionally, these anthropogenic dust activities are sensitive to the drought severity over the fallowed lands, with potential cumulative effects on downstream dust burden during prolonged multi-year drought conditions. Overall, our results have important implications for public health, including increased risk for Valley fever and for regional climates, such as increases in extreme precipitation and snowmelt over the Sierra Nevada.

Climate sciences↗
NASA NTRS1995

Analysis methods for wind turbine control and electrical system dynamics

The integration of new energy technologies into electric power systems requires methods which recognize the full range of dynamic events in both the new generating unit and the power system. Since new energy technologies are initially perceived as small contributors to large systems, little attention is generally paid to system integration, i.e. dynamic events in the power system are ignored. As a result, most new energy sources are only capable of base-load operation, i.e. they have no load following or cycling capability. Wind turbines are no exception. Greater awareness of this implicit (and often unnecessary) limitation is needed. Analysis methods are recommended which include very low penetration (infinite bus) as well as very high penetration (stand-alone) scenarios.

Hinrichsen, E. N.↗
DOE OSTI2023

DP-Based Optimization of BESS to Substitute RICE Reserves for Improved Economic Benefits

In regions with a significant proportion of renewable energy sources, reciprocating internal combustion engines (RICE or IC) are frequently employed as grid reserve solutions. This study suggests a scheduling technique for battery energy storage systems (BESS) that utilizes dynamic programming (DP) to replace IC units, adding versatility to the total reserve capacity and resulting in higher BESS earnings compared to traditional scheduling approaches. The capacity sweeping approach to optimize BESS’s revenue sometimes leads to operation which is counter intuitive to conventional battery sizing approaches.

25 ENERGY STORAGE↗
NASA NTRSDate not supplied

Autonomous Nanosatellite Swarming (ANS) Using Radio Frequency and Optical Navigation

The Autonomous Nanosatellite Swarming (ANS) subsystem is a dynamics, guidance, navigation, and control software unit based on commercial-off-the-shelf avionics, which is distributed and deployed on a cluster of cooperative nanosatellites to enable future exploration missions to near-Earth objects (NEOs). ANS allows multiple nanosatellites to act together as a single large spacecraft through radio-frequency (RF) cross-links, optical navigation, and low-thrust control. In particular, ANS-equipped nanosatellites can utilize images of a NEO taken from multiple viewpoints to achieve stereo-vision capabilities and characterize the target body’s shape, gravity field, and dynamical properties. All this is enabled by the fusion of major advances to the state of the art in the astrodynamics of relative motion, RF and vision-based navigation, as well as low-thrust control.

Simone D'Amico↗
DOE OSTI2020

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

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

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗
DOE OSTI2020

Three-Dimensional Dynamics of a Magnetic Hopfion Driven by Spin Transfer Torque

Magnetic hopfion is a three-dimensional (3D) topological soliton with novel spin structure that would enable exotic dynamics. Herein, we observe and research the current-driven 3D dynamics of a magnetic hopfion with a unit Hopf index in a frustrated magnet. Attributed to the spin Berry phase and symmetry of the hopfion, the phase space entangles multiple collective coordinates, thus the hopfion exhibits rich dynamics including longitudinal motion along the current direction, transverse motion perpendicular to the current direction, rotational motion, and dilation. Furthermore, the characteristics of hopfion dynamics is determined by the ratio between the nonadiabatic spin transfer torque parameter and the damping parameter. Such peculiar 3D dynamics of magnetic hopfion could shed light on understanding the universal physics of hopfions in different systems and boost the prosperous development of 3D spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Quantifying atmospheric carbon removal at pulp and paper mills: a life cycle assessment across system boundaries

The pulp and paper industry is a promising yet underexplored platform for large-scale carbon dioxide removal (CDR) due to its use of biogenic feedstocks and production of concentrated CO 2 emissions from point sources. This study presents the first comprehensive life cycle assessment (LCA) of retrofitting an amine-based carbon capture and storage (CCS) system into a representative virgin kraft pulp and paper mill in the Southeastern U.S. We evaluate carbon removal across five system configurations, applying both static and dynamic LCA methods under multiple functional units: CO 2 captured, biomass input, and paper output. Results show that CCS retrofits can convert a conventional mill from a net emitter into a net carbon sink, with total removal efficiencies from 17% to 92% (metric tonnes of CO 2 removed per metric tonne of CO 2 available for removal under selected boundary conditions). When carbon removal is normalized to the quantity of biogenic CO 2 captured—a narrow, gate-to-gate system boundary that considers only CCS facility emissions—removal efficiencies reached as high as 92%. The use of such narrow boundaries aligns with precedents in traditional LCA methodology, where gate-to-gate assessments are commonly applied to isolate process-level performance and allocate emissions accordingly, providing a consistent basis for comparison across technologies. Under broader cradle-to-grave boundaries—which begin tracking carbon at the point of its physical removal from the atmosphere via photosynthesis in the forest, and extend to include upstream forest operations, mill-wide emissions, and downstream product decomposition—efficiencies declined, ranging from 17% to 46% under static assumptions and dropping to 12% when accounting for dynamic biogenic carbon fluxes over time. These results underscore how system boundary definitions influence reported outcomes, while also illustrating the complementary roles of narrow and broad perspectives for different decision-making contexts.

09 BIOMASS FUELS↗
DOE OSTI2022

Subduction, Underplating, and Return Flow Recorded in the Cycladic Blueschist Unit Exposed on Syros, Greece

Exhumed high-pressure/low-temperature (HP/LT) metamorphic rocks provide insights into deep (~20–70 km) subduction interface dynamics. On Syros Island (Cyclades, Greece), the Cycladic Blueschist Unit preserves blueschist-to-eclogite facies oceanic- and continental-affinity rocks that record the structural and thermal evolution linked to Eocene subduction. Despite decades of research, the metamorphic and deformation history (P-T-D) and timing of subduction and exhumation are matters of ongoing discussion. We suggest that Syros comprises three coherent tectonic slices and that each slice underwent subduction, underplating, and syn-subduction return flow along similar P-T trajectories, but at progressively younger times. Subduction and exhumation are distinguished by lineations and ductile fold axis orientations, and are kinematically consistent with previous studies that document top-to-the-S-SW shear (prograde-to-peak subduction), top-to-the-NE shear (blueschist facies exhumation), and then E-W coaxial stretching (greenschist facies exhumation). Amphibole zonations record cooling during decompression, indicating return flow above a cold slab. Multi-mineral Rb-Sr isochrons and compiled metamorphic geochronology show that the three slices record distinct stages of peak subduction (53–52, ~50, and 45 Ma) that young with structural depth. Retrograde blueschist and greenschist facies fabrics span ~50–40 and ~43–20 Ma, respectively, and also young with structural depth. Synthesized data sets support a revised tectonic framework for Syros, involving subduction of structurally distinct coherent slices and simultaneous return flow of previously accreted tectonic slices in the subduction channel shear zone. Distributed, ductile, dominantly coaxial return flow in an Eocene-Oligocene subduction channel proceeded at rates of ~1.5–5 mm/yr and accommodated ~80% of the total exhumation of this HP/LT complex.

58 GEOSCIENCES↗
NASA NTRS2014

Benefits Analysis of Multi-Center Dynamic Weather Routes

Dynamic weather routes are flight plan corrections that can provide airborne flights more than user-specified minutes of flying-time savings, compared to their current flight plan. These routes are computed from the aircraft's current location to a flight plan fix downstream (within a predefined limit region), while avoiding forecasted convective weather regions. The Dynamic Weather Routes automation has been continuously running with live air traffic data for a field evaluation at the American Airlines Integrated Operations Center in Fort Worth, TX since July 31, 2012, where flights within the Fort Worth Air Route Traffic Control Center are evaluated for time savings. This paper extends the methodology to all Centers in United States and presents benefits analysis of Dynamic Weather Routes automation, if it was implemented in multiple airspace Centers individually and concurrently. The current computation of dynamic weather routes requires a limit rectangle so that a downstream capture fix can be selected, preventing very large route changes spanning several Centers. In this paper, first, a method of computing a limit polygon (as opposed to a rectangle used for Fort Worth Center) is described for each of the 20 Centers in the National Airspace System. The Future ATM Concepts Evaluation Tool, a nationwide simulation and analysis tool, is used for this purpose. After a comparison of results with the Center-based Dynamic Weather Routes automation in Fort Worth Center, results are presented for 11 Centers in the contiguous United States. These Centers are generally most impacted by convective weather. A breakdown of individual Center and airline savings is presented and the results indicate an overall average savings of about 10 minutes of flying time are obtained per flight.

NAS based Benefits Analysis↗
DOE OSTI2025

Role of Diurnal Cycle of Insolation on the MJO Propagation in the Maritime Continent

The diurnal cycle of convection in the Maritime Continent (MC) has been hypothesized to act as a barrier to the eastward propagation of the Madden‐Julian oscillation (MJO). To test this hypothesis, we use a regional model with realistic MJO to simulate an event from the boreal spring of 2013 that weakened and stalled over the MC. Two simulations are conducted: one that includes the diurnal cycle of insolation (CTL), and another without it (NO_DC). The MJO in the simulations was identified and tracked using a large‐scale precipitation tracking method that distinguishes propagation and non‐propagation unlike the usual Real‐time Multivariate MJO method. In the NO_DC simulation, the absence of diurnal heating reduces land precipitation, allowing more continuous eastward MJO propagation. An analysis of moist static energy budget reveals that MJO maintenance in NO_DC is due to increased longwave heating and reduced advection, whereas the persistent MJO propagation in NO_DC is due to increased advection and reduced longwave heating and surface latent heat flux. These processes, however, may vary across different parts of the MC, emphasizing the complexity of MJO propagation across the MC.

Zhou, Xin [National Center for Atmospheric Researc↗
DOE OSTI2026

A scalable variational method for estimating the latent infection-rate field of an outbreak

In this paper, we explore whether the infection-rate of a disease can serve as a robust monitoring variable in epidemiological surveillance algorithms. The infection-rate is dependent on population mixing patterns that do not vary erratically day-to-day; in contrast, daily case-counts used in contemporary surveillance algorithms are corrupted by reporting errors. The technical challenge lies in estimating the latent infection-rate from case-counts. Here we devise a Bayesian method to estimate the infection-rate across multiple adjoining areal units, and then use it, via an anomaly detector, to discern a change in epidemiological dynamics. We extend an existing model for estimating the infection-rate in an areal unit by incorporating a Markov random field model, so that we may estimate infection-rates across multiple areal units, while preserving spatial correlations observed in the epidemiological dynamics. To carry out the high-dimensional Bayesian inverse problem, we develop an implementation of mean-field variational inference specific to the infection model and integrate it with the random field model to incorporate correlations across counties. The method is tested on estimating the COVID-19 infection-rates across all 33 counties in New Mexico using data from the summer of 2020, and then employing them to detect the arrival of the Fall 2020 COVID-19 wave. We perform the detection using a temporal algorithm that is applied county-by-county. We also show how the infection-rate field can be used to cluster counties with similar epidemiological dynamics.

60 APPLIED LIFE SCIENCES↗