Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mobility simulation”

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

Nonsolvating Fluoroaromatic Cosolvent Enabled Long-Term Cycling of High-Voltage Lithium-Ion Batteries with Organosulfur Electrolytes

Here, the structure–activity relationships of nonsolvating cosolvents for organosulfur-based electrolyte systems were revealed. The performance of nonsolvating dilutant fluorobenzene (FB) was compared to various fluorinated ether dilutants in high-voltage electrolytes containing a concentration of 1.2 M LiPF 6 dissolved in fluoroethylene carbonate (FEC), ethyl methyl sulfone (EMS), and the dilutant. In a high-voltage and high-loading LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) full cell configuration, the organosulfur-based electrolyte containing FB dilutant enabled superior electrochemical performance compared to the electrolytes using other nonsolvating fluorinated ether formulations. Moreover, the FB-containing electrolyte exhibited the highest ionic conductivity and lowest viscosity among all organosulfur-based electrolytes containing nonsolvating dilutant. These improvements are attributed to the enhanced physical properties of electrolyte and lithium-ion mobility. Furthermore, by employing first-principles simulations, the observed suppression of side reactions at high voltage is linked to FB’s lower reactivity toward singlet dioxygen, which is likely produced at the NMC interface. Overall, FB is considered an excellent diluent that does not impede cell operation by mass decomposition at the cathode.

25 ENERGY STORAGE↗
DOE OSTI2025

Skyrmion-skyrmionium phase separation and laning transitions via spin-orbit torque currents

Many driven binary systems can exhibit laning transitions when the two species have different mobilities, such as colloidal particles with opposite charges in electric fields. Another example is pedestrian or active-matter systems, where particles moving in opposite directions form a phase-separated state that enhances the overall mobility. In this work, we use atomistic simulations to demonstrate that mixtures of skyrmions and skyrmioniums also exhibit pattern formation and laning transitions. Skyrmions move more slowly and at a finite skyrmion Hall angle compared to skyrmioniums, which move faster and without a skyrmion Hall effect. At low drives, the system forms a partially jammed phase where the skyrmioniums are dragged by the surrounding skyrmions, resulting in a finite angle of motion for the skyrmioniums. At higher drives, the system transitions into a laned state, but unlike colloidal systems, the lanes in the skyrmion-skyrmionium mixture are tilted relative to the driving direction due to the intrinsic skyrmion Hall angle. In the laned state, the skyrmionium angle of motion is reversed when it aligns with the tilted lane structure. At even higher drives, the skyrmioniums collapse into skyrmions. Below a critical skyrmion density, both textures can move independently with few collisions, but above this density, the laning state disappears entirely, and the system transitions to a skyrmion-only state. We map out the velocity and Hall responses of the different textures and identify three distinct phases: partially jammed, laned, and skyrmion-only moving crystal states. In conclusion, we compare our results to recent observations of tilted laning phases in pedestrian flows, where chiral symmetry breaking in the particle interactions leads to similar behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Interactions Between Climate Policy and Technology-influenced Travel Behavior: Mitigating Induced Demand from CACC

Advances in vehicle technology have influenced the development of automated vehicle systems, where vehicles that do not require human intervention are already deployed in the roadway networks. While these advances are proved to increase roadway safety and highway capacity, more research is needed to understand the long-term and regional-level impacts on mobility, land use, energy consumption, and emissions. This study proposes a multi-model approach to analyze the effect of vehicle automation and deep decarbonization policies over a period from 2020 to 2040 in Austin, Texas. We use the Global Change Analysis Model (GCAM) to develop internally the scenarios that are then passed to the SMART Mobility modeling workflow, a large-scale simulation framework combining the POLARIS activity-based travel demand model and mesoscopic traffic simulator with the Autonomie vehicle energy consumption model and the UrbanSim land use simulator. Results suggest that the introduction of vehicles with advanced automation could increase fuel consumption when no decarbonization policies are implemented. Also, advances in vehicle technology research and development could lead to a decline in energy use in the long-term. Energy pricing and vehicle electrification incentives could help reduce the impact of vehicle automation. Finally, our analysis indicates the relevance of introducing land use processes in longterm vehicle automation studies.

land use↗
DOE OSTI2025

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

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

Perspective: Are glasses really frozen?

A glass is obtained by cooling a liquid fast enough to avoid crystallization. It is a solid with elastic constants and density similar to those of a crystal with the same composition. It is customarily assumed that the structure of a glass is frozen, and atoms are largely locked into a random, but specific, structure; while a small number of loosely bound atoms may be mobile, most atoms are frozen at room temperature and below. However, recent simulation studies on metallic glasses present a rather different picture. A very significant fraction (~20%) of atoms are not frozen and locally mobile, even down to T = 0 partly because of quantum effect. Here, we describe these observations and discuss their implications on the nature of the glass structure and the origin of the glass transition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomistic Insights into the Reactive Diffusion of CO 2 in Guanidine-Based Facilitated Transport Membranes

The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.

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

Structural origin of disorder-induced ion conduction in NaFePO 4 cathode materials

Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.

Christensen, Rasmus↗
DOE OSTI2024

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗
DOE OSTI2024

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

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

What is the Role of Relative Humidity on Conductivity in Polymer Electrolytes?

Ion transport properties in polymer electrolytes have been widely studied under rigorously dry and highly water-swollen conditions. However, the transition between these extremes, i.e., the low hydration regime, is still poorly understood at the molecular level and relevant to applications ranging from battery electrolytes to electrophoretic separations. Here, in this study, we apply atomistic molecular dynamics simulations to probe the role of hydration on ion mobilities and conductivities in LiTFSI-doped polyethers at low water content (less than 10% water by volume), which corresponds to 0−80% relative humidity conditions. With increasing water content, Li + ions exhibit two distinct regimes in their mobilities. At low water content, the majority of the Li + ions are weakly hydrated (two or fewer water molecules per Li + ion), maintaining strong interactions with their polymer solvation cage, and therefore exhibit only a weak increase in mobility relative to dry conditions. As water content increases, Li + ions form complete hydration shells, residing within water-rich domains that promote faster mobilities. Further, Li + mobilities increase more significantly in polymer electrolytes composed of more hydrophilic polymers and higher salt concentrations due to the formation of larger water-rich domains. In contrast, TFSI − has comparatively weak interactions with both the polymer and water, exhibiting a monotonic increase in mobility as a function of water content. Overall, these results help clarify the molecular mechanisms underlying ionic mobilities and conductivities in polymer electrolytes at low water content.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Cellulose acetate membranes exhibit exceptional monovalent to divalent cation selectivities

Salt transport properties of cellulose acetate membranes are reported for a series of chloride salts with mono- valent and divalent cations (LiCl, NaCl, MgCl 2 , CaCl 2 ). Measurements include salt permeability and sorption, with diffusivity values calculated from the permeability and sorption results. We report an exceptionally high LiCl/MgCl 2 selectivity of 750:1. Salts with similar valence (LiCl and NaCl; MgCl 2 and CaCl 2 ) have similar transport properties. The high monovalent/divalent selectivity arises from differences in both sorption and diffusion, with a LiCl/MgCl 2 solubility selectivity of about 11 and a diffusivity selectivity of about 70. Atomistic molecular dynamics simulations show that ions tend to reside in isolated clusters of water. Increasing ion charge strengthens ion–water interactions relative to ion–polymer interactions, explaining the reduced sorption of divalent ions. Diffusion of ions through the membrane occurs via hop-like motion between water clusters. Lithium diffuses faster than magnesium due to weaker ion–water coordination for lithium, which allows for greater mobility within water clusters and more frequent hopping. Altogether, our atomistic simulations suggest that the high LiCl/MgCl 2 selectivity is linked to cellulose acetate’s high water/salt selectivity and is a consequence of low water content and relatively uniform water distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2025

Spatial Impacts of Electric Vehicle Charging on Power Grid Stability: A Downtown Atlanta Case Study

The rapid increase in electric vehicle (EV) charging demand poses a potential risk to power grid stability, particularly as the spatial distribution of this demand remains underexplored. Existing research often focuses on technical optimization models while overlooking the geographic and human dynamics that affect energy consumption. This study addresses this gap by incorporating mobility data to estimate both building energy use and EV charging demand while also considering geographic factors for a better understanding of grid load. Using agent-based simulations and the Open-Source Distribution System Simulator, the study evaluates the effect of various EV penetration scenarios on grid voltage and unbalance. The results show that, although voltage remains within acceptable limits at lower EV penetration rates, significant voltage drop and unbalance occur as EV penetration exceeds 40%, particularly in residential areas with high charging demand. This study offers a framework for integrating spatial analysis and mobility data in power network simulations, providing insights for future EV infrastructure planning.

Pan, Melrose [ORNL] (ORCID:000000031627448X)↗
DOE OSTI2025

The CO2 mineralization potential of New Mexico basalts and associated critical element behavior and mobility.

This presentation reports on laboratory experiments assessing the carbon dioxide (CO₂) mineralization potential of New Mexico basalts, including trachybasalt and basalt lithologies. Using simulated groundwater and controlled conditions, the study evaluates the kinetics of CO₂ sequestration, mineral saturation indices, and secondary mineral formation. It also investigates the mobility of critical and trace elements such as Ni, Zn, Cu, and Li during basalt dissolution. Results indicate that New Mexico basalts can effectively sequester CO₂ through carbonate formation, with grain size and mineral composition influencing reaction rates and elemental behavior. These findings support the viability of geologic carbon storage in mafic volcanic terrains and inform strategies for critical element recovery.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Role of effective mass and long-range interactions in the band-gap renormalization of photoexcited semiconductors

Understanding how to control changes in the electronic structure and related dynamical renormalizations by external driving fields is the key for understanding ultrafast spectroscopy and applications in electronics. Here, we focus on the band gap's modulation by external electric fields and uncover the effect of band dispersion on the gap renormalization. We employ the Green's function formalism using the real-time Dyson expansion to account for dynamical correlations induced by photodoping. The many-body formalism captures the dynamics of systems with long-range interactions, carrier mobility, and variable electron and hole effective mass. We also demonstrate that mean-field simulations based on the Hartree-Fock Hamiltonian, which lacks dynamical correlations, yields a qualitatively incorrect picture of band-gap renormalization. We find the trend that increasing effective mass, thus decreasing mobility, leads to as much as a 6% enhancement in band-gap renormalization. Further, the renormalization is strongly dependent on the degree of photodoping. As the screening induced by free electrons and holes effectively reduces any long-range and interband interactions for highly excited systems, we show that there is a specific turnover point with a minimal band gap. Here, we further demonstrate that the optical gap renormalization follows the same trend though its magnitude is altered by the Moss-Burstein effect.

Approximation methods for many-body systems↗
DOE OSTI2026

NLR Core Modeling & Decision Support Capabilities: FASTSim, RouteE, T3CO & OpenPATH

This project is part of the program area to develop and improve core capabilities for the Energy-Efficient Mobility Systems (EEMS) program that enable research, development and deployment of advanced mobility solutions and enhance the EEMS Program's ability to address system-level transportation challenges. Advancements to the Future Automotive Systems Technology Simulator (FASTSim), Route Energy Prediction Model (RouteE), Transportation Technology Total Cost of Ownership (T3CO) and Open Platform for Agile Trip Heuristics (OpenPATH) core capabilities under this project supports the overall EEMS Program goals to effectively evaluate energy and mobility impacts of future transportation technologies and services, and to identify the most promising pathways to reduce transportation costs and environmental harms, and to improve mobility access. This presentation was prepared for the 2026 Annual Merit Review of this project.

33 ADVANCED PROPULSION SYSTEMS↗