Search NASA⌕ Search

SEARCH · Search NASA

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

JANUS: Resilient and Adaptive Data Transmission for Enabling Timely and Efficient Cross-Facility Scientific Workflows

In modern science, the growing complexity of large-scale scientific projects has led to an increasing reliance on cross-facility scientific workflows, where resources and expertise from multiple institutions and geographic locations are leveraged to accelerate scientific discovery. These workflows often require transmitting huge amounts of scientific data through wide-area networks. Although high-speed networks like ESnet and transfer services such as Globus have improved data mobility, several challenges remain. The sheer volume of data can overwhelm network bandwidth, widely used transport protocols such as TCP suffer from inefficiencies due to retransmissions triggered by packet loss, and existing fault-tolerance mechanisms like erasure coding introduce substantial overhead. In this paper, we propose Janus, a resilient and adaptable data transmission approach designed for cross-facility scientific workflows. Unlike traditional TCP-based methods, Janus leverages UDP, integrates erasure coding for fault tolerance, and combines it with error-bounded lossy compression to reduce overhead. This novel design allows users to balance data transmission time and accuracy, optimizing transfer performance based on specific scientific requirements. Additionally, Janus dynamically adjusts erasure coding parameters in response to real-time network conditions, ensuring efficient data transfers even in fluctuating environments. We develop optimization models for determining ideal configurations and implement adaptive data transfer protocols to enhance reliability. Through extensive simulations and real-network experiments, we demonstrate that Janus significantly improves transfer efficiency while maintaining data fidelity.

Esaulov, Vladislav [Georgia State University, Atla↗
DOE OSTI2025

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗
DOE OSTI2025

Real-Time, Adaptive Radiological Anomaly Detection and Isotope Identification Using Non-Negative Matrix Factorization

Spectroscopic anomaly detection and isotope identification algorithms are integral components in nuclear nonproliferation applications such as search operations. The task is especially challenging in the case of mobile detector systems because the observed gamma-ray background changes more than for a static detector system, and a pretrained background model can easily find itself out of domain. The result is that algorithms may exceed their intended false alarm rate or sacrifice detection sensitivity to maintain the desired false alarm rate. Non-negative matrix factorization (NMF) is a powerful tool for spectral anomaly detection and identification, but, like many similar algorithms that rely on data-driven background models, in its conventional implementation, it is unable to update in real time to account for environmental changes that affect the background spectroscopic signature. Here, we have developed a novel NMF-based algorithm that periodically updates its background model to accommodate changing environmental conditions. The adaptive NMF algorithm involves fewer assumptions about its environment, making it more generalizable than existing NMF-based methods while maintaining or exceeding detection performance on simulated and real-world datasets.

Anomaly detection↗
DOE OSTI2025

Sliding dynamics of skyrmion molecular crystals

Using both atomistic and particle-based simulations, we investigate the current-driven dynamics of skyrmions on two-dimensional periodic substrates when there are multiple skyrmions per substrate minimum. At zero drive, the system forms pinned skyrmion molecular crystal states consisting of dimers, trimers, or dimer-trimer mixtures that have both positional and orientational order. On a square substrate lattice, the motion above depinning occurs via a running soliton that travels completely transverse to the applied current. This motion is generated by a torque from the Magnus force, which rotates the n-mer states perpendicular to the applied current. At higher drives, the flow becomes disordered while the Hall angle diminishes and gradually approaches the intrinsic value. In some cases, we also find directional locking where the Hall angle becomes locked to certain symmetry directions of the substrate over a range of currents. The transitions into and out of directionally locked states are accompanied by negative differential mobility in which the net velocity decreases as the drive increases. In conclusion, on a triangular substrate, we find no transverse mobility effects, but still observe directionally locked motion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Studying Aerosol, Clouds, and Air Quality in the Coastal Urban Environment of Southeastern Texas

A multi-agency succession of field campaigns was conducted in southeastern Texas during July 2021 through October 2022 to study the complex interactions of aerosols, clouds and air pollution in the coastal urban environment. As part of the Tracking Aerosol Convection interactions Experiment (TRACER), the TRACER- Air Quality (TAQ) campaign the Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) and the Convective Cloud Urban Boundary Layer Experiment (CUBE), a combination of ground-based supersites and mobile laboratories, shipborne measurements and aircraft-based instrumentation were deployed. These diverse platforms collected high-resolution data to characterize the aerosol microphysics and chemistry, cloud and precipitation micro- and macro-physical properties, environmental thermodynamics and air quality-relevant constituents that are being used in follow-on analysis and modeling activities. We present the overall deployment setups, a summary of the campaign conditions and a sampling of early research results related to: (a) aerosol precursors in the urban environment, (b) influences of local meteorology on air pollution, (c) detailed observations of the sea breeze circulation, (d) retrieved supersaturation in convective updrafts, (e) characterizing the convective updraft lifecycle, (f) variability in lightning characteristics of convective storms and (g) urban influences on surface energy fluxes. The work concludes with discussion of future research activities highlighted by the TRACER model-intercomparison project to explore the representation of aerosol-convective interactions in high-resolution simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Selenium interaction with iron minerals: Quantitative comparison of sorption and coprecipitation impacts on mobility

Given the significance of selenium (Se) as a micronutrient, the radioactive nature of some of its isotopes, and its affinity to iron (Fe) minerals, extensive research has been conducted on the sorption mechanisms between Se and these minerals. Here, in this study, we employ sorption data sourced from the L-SCIE database and coprecipitation data from available literature to achieve the following objectives: i) establish coherence between adsorption and coprecipitation processes, ii) quantitatively evaluate the importance of these processes in nuclear waste repository science, and iii) propose a forward-looking approach for integrating coprecipitation into reactive transport models. Our findings indicate that a correlation between Se adsorption and coprecipitation can be established using the λ formalism. The comparable log(λ Se(IV) /λ Se(VI) ) ratios derived from adsorption and coprecipitation experiments suggest that these processes can be quantitatively compared and evaluated using our numerical approach. Across all iron oxide phases examined, coprecipitation leads to significantly greater immobilization of Se compared to adsorption. Specifically, for hydrous ferric oxide, hematite, and goethite, coprecipitation is predicted to result in 100–1000 times more Se immobilization compared to adsorption, irrespective of the Se oxidation state (Se(IV) or Se(VI)); notably stronger immobilization potential via coprecipitation was observed for magnetite. The modeling approach and quantitative analysis presented herein clearly highlight the importance of including coprecipitation processes when simulating Se (and other elements) transport, particularly under conditions where mineral compositions are transient or evolving with time. Neglecting coprecipitation in models is likely to lead to significant overestimates of migration.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Automated Signal Timing Plan Reconstruction Using High-Resolution Event-Based Controller Data for Digital Twins

Transportation digital twins are essential tools for evaluating emerging technologies such as connected and automated vehicles, adaptive traffic signal control, and mobility optimization strategies. Realistic digital twins require accurate emulation of real-world signal controllers and detailed signal timing plans. However, signal timing plans are often unavailable or difficult to access, forcing researchers and modelers to rely on assumed fixed timings or halt their analysis. To overcome this challenge, we present a method that directly estimates signal timing plan parameters using high-resolution, event-based data from traffic signal controllers. The proposed method extracts key parameters, including cycle length, offset, phase sequence, coordinated phases, phase-specific minimum and maximum green durations, vehicle extensions, and splits under coordination. A rule-based deterministic signal timing reconstruction algorithm based on traffic signal operation rules, such as those outlined in the Signal Timing Manual, is developed and validated. We evaluate this method, which uses high-resolution controller event logs and verified signal timing plans, on 94 signalized intersections in Nashville, Tennessee, demonstrating their ability to generate accurate, simulation-ready signal timing plans for tools such as SUMO and Vissim.

Saroj, Abhilasha [ORNL] (ORCID:0000000191178063)↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2026

Simultaneous Reduction of NOx and Fuel Consumption for Off-Road Powertrains

Simultaneously reducing criteria pollutants and fuel consumption is important for clean air and improving vehicle total cost of ownership. The goal of this effort was focused on a 90% NOx reduction and 10% fuel savings for an off-road 407 kW diesel engine. The baseline was a production Fiat Powertrain 13L engine and aftertreatment system meeting 0.4 g/kW-hr NOx. The baseline system was quantified over the NRTC, RMC, new low load cycle and five field cycles. A next generation engine was built incorporating several fuel-efficient design features, including a higher compression ratio, increased fuel-rail pressure, low-friction piston rings, and a high-efficiency variable-geometry turbocharger. Cylinder deactivation and EGR pump technologies were added to this engine as well. The combination was optimized prior to adding advanced aftertreatment systems, showing the trade-off of engine out NOx and exhaust temperature. Two next-generation catalyst technologies were employed into a LO-SCR plus main SCR system, both with and without an electric heater upstream of the LO-SCR. These catalysts were hydrothermally aged to simulate significant field use. Dual SCR dosing with newly developed controls played a critical role in achieving the proper split between the upstream LO-SCR and the downstream main SCR. Adding a next generation mixer for the downstream SCR proved essential in obtaining the final results. The optimal configuration required adding an electric heater to elevate the exhaust temperature at the LO-SCR for early cycle NOx reduction. The final results showed a 94.8% NOx reduction and 15.7% fuel savings on the composite NRTC.

McCarthy, James [Eaton Corporation]↗
DOE OSTI2025

The effect of ion pairing on speciation and transport in ion exchange membranes at varying hydration levels: A four-state model

Understanding ion pairing in ion exchange membranes (IEMs) is essential for advancing IEM applications in energy and environmental technologies. Here, this study introduces a four-state molecular dynamics model to quantify speciation and transport within Nafion-117, specifically examining the role of ion pairing in monovalent and divalent counterions (NaCl, Na 2 SO 4 , and MgSO 4 ). By analyzing radial distribution functions (RDFs) and molecular snapshots, we distinguish ion pairing modes and classify counterions into four states: condensed counterion, condensed ion pair, free ion pair, and free counterion. A key finding is that while divalent counterions (e. g., Mg 2+ ) maintain stable speciation across hydration levels, monovalent counterions (e.g., Na + ) show notable speciation shifts with hydration. Both monovalent and divalent counterions are not diffusive when condensed onto the polymer (sorbed to membrane functional groups). In contrast, free counterions are diffusive across all hydration levels. To evaluate the overall diffusivity of counterions, four-state fractions and diffusivities are computed, each contributing to counterion transport. The condensed/free ion speciation for multivalent sulfate salts aligns with previous revisions to the Donnan-Manning framework that include ion pairing, thereby validating its relevance to established membrane theories. The four-state model's diffusivity results support several current ion exchange assumptions, including that the condensed counterions are immobile, while uncondensed counterions are mobile. The four-state model offers insights into contact ion pairing within IEMs, highlighting its potential even when undetected in aqueous solution experiments. This work advances the theoretical understanding of counterion speciation in IEMs while identifying model limitations that suggest avenues for refinement, such as distinguishing water-mediated ion pairs between fully hydrated ions.

Ion exchange membranes↗
DOE OSTI2024

Atomistic simulations to reveal HIP-bonding mechanisms of Al6061/Al6061

Molecular dynamics simulations were employed to understand the diffusion bonding process during hot isostatic pressing (HIP) of Al6061/Al6061 alloy. Simulations of the HIP process reveal atomistic phenomena that are difficult or unlikely to be observed experimentally and provide useful insights into the mechanism of diffusion and bonding. Here, the results reveal that at the start of the HIP process, a massive incursion of oxygen atoms occurs from the pre-existing γ-Al 2 O 3 to the 6061 region across the interphase interface. These oxygen atoms interact with the enriched Mg atom layer present at the existing γ-Al 2 O 3 and 6061 matrix to form a secondary complex Mg 2 Al 2 O 5 phase. Diffusion calculations also show that transport of atoms due to the applied pressure is 4–5 orders of magnitude higher than would occur in the absence of HIP conditions. The Mg 2 Al 2 O 5 phase also provides efficient pathways for the rapid transport of Mg atoms. Because of the higher diffusion coefficients observed for Mg within the phase, Mg atoms can move more swiftly compared to their diffusion within other phases such as γ-Al 2 O 3 . This accelerated mobility facilitates the rapid movement of Mg atoms across the interface, leading to changes in the local composition and the potential growth of the Mg 2 Al 2 O 5 phase.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elucidating the Role of Hydrogen at c‐In 2 O 3 /a‐In 2 O 3− x Interface

Ab initio molecular dynamics simulations and hybrid functional electronic structure calculations are employed to determine the formation, the structural and electronic properties, and the dynamics of covalent (In—OH) and ionic (In—H—In) hydrogen defects at crystalline-In 2 O 3 /amorphous-In 2 O 3−x interface. This comprehensive computational study considers (i) various interstitial and substitutional hydrogen site locations within the crystalline, amorphous, and interfacial regions; (ii) several oxygen-to-hydrogen ratios; and (iii) possible defect charged states. The results reveal hydrogen's inability to fully passivate the undercoordinated under-shared in atoms in amorphous highly substoichiometric oxide, giving rise to the formation of deep electron traps even in net-charge neutral structures. These trap defects are found to be sensitive to photoexcitation, in contrast to In—OH with H's electronic states located below the valence band and to In—H—In, where H is found to maintain its charge state upon illumination. Nevertheless, H plays a critical role in photoinduced conductivity and its relaxation by promoting In—O coordination transformations at the interfacial region, including deterioration of the crystalline layer. Finally, the results help identify mobile H species and metastable H defect complexes (such as In—H—H—In, In—OH–H—In, and In—OH–O—In) that are responsible for long relaxation times of the conductivity decay.

H defects↗
DOE OSTI2024

Adsorption of neodymium, dysprosium, and ytterbium to goethite under varying aqueous chemistry conditions [dataset]

The adsorption of rare earth elements (REEs) to iron oxides can regulate the mobility of REEs in the environment and is heavily influenced by water chemistry. This study utilized batch experiments to examine the adsorption of Nd, Dy, and Yb to goethite under varying pH, electrolyte (type and concentration), and concentrations of dissolved inorganic carbon and citrate. REE adsorption was strongly influenced by pH, with an increase from essentially no adsorption at pH 3.0 to nearly complete adsorption at pH 6.5 and higher. Citrate enhanced the adsorption of REEs at low pH (< 5.0), likely by forming goethite-REE-citrate ternary surface complexes. However, citrate inhibited the adsorption of REEs at higher pH (> 5.0) by forming aqueous REE-citrate complexes. Ionic strength had a small influence on REE adsorption, and the presence of dissolved inorganic carbon had no discernible effect. Equilibrium adsorption was interpreted with a triple layer surface complexation model (SCM). The selection of surface complexation reactions was guided by extended X-ray absorption fine structure spectra. A SCM with a single bidentate inner-sphere surface complexation reaction for Nd and two inner-sphere surface complexation reactions (one monodentate and one bidentate reaction) for Dy and Yb effectively simulated adsorption across a broad range of conditions in the absence of citrate. Accounting for the effects of citrate on REE adsorption required the addition of up to two ternary REE-citrate-goethite surface complexes. The SCM can enable predictions of REE transport in subsurface environments that have goethite as an important adsorbent mineral. This predictive capability could contribute to identifying potential REE sources and facilitating efficient extraction of REEs.

58 GEOSCIENCES↗
DOE OSTI2024

Adsorption of Neodymium, Dysprosium, and Ytterbium to Goethite under Varying Aqueous Chemistry Conditions

The adsorption of rare earth elements (REEs) to iron oxides can regulate the mobility of REEs in the environment and is heavily influenced by water chemistry. This study utilized batch experiments to examine the adsorption of Nd, Dy, and Yb to goethite under varying pH, electrolyte (type and concentration), and concentrations of dissolved inorganic carbon and citrate. REE adsorption was strongly influenced by pH, with an increase from essentially no adsorption at pH 3.0 to nearly complete adsorption at pH 6.5 and higher. Citrate enhanced the adsorption of REEs at low pH (<5.0), likely by forming goethite-REE-citrate ternary surface complexes. However, citrate inhibited the adsorption of REEs at higher pH (>5.0) by forming aqueous REE-citrate complexes. Ionic strength had a small influence on REE adsorption, and the presence of dissolved inorganic carbon had no discernible effect. Equilibrium adsorption was interpreted with a triple-layer surface complexation model (SCM). The selection of surface complexation reactions was guided by extended X-ray absorption fine structure spectra. An SCM with a single bidentate inner-sphere surface complexation reaction for Nd and two inner-sphere surface complexation reactions (one monodentate and one bidentate reaction) for Dy and Yb effectively simulated adsorption across a broad range of conditions in the absence of citrate. Accounting for the effects of citrate on REE adsorption required the addition of up to two ternary REE-citrate-goethite surface complexes. The SCM can enable predictions of REE transport in subsurface environments that have goethite as an important adsorbent mineral. Furthermore, this predictive capability could contribute to identifying potential REE sources and facilitating efficient extraction of REEs.

58 GEOSCIENCES↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2026

Impact of Rock Wettability and Mineralization on CO 2 Storage Efficiency in Basalt Reservoirs

Basalt formations have emerged as highly promising targets for geological CO 2 storage due to their abundance of reactive silicate minerals that can rapidly convert dissolved CO 2 into stable carbonate minerals. The efficiency of CO 2 trapping in basalts, however, depends not only on their geochemical reactivity but also on reservoir wettability, which governs how CO 2 partitions among structural, residual, solubility, and mineral trapping mechanisms. However, wettability and mineral kinetics have largely been examined separately, wettability for plume migration and kinetics in terms of mineralization, leaving their coupled impact unquantified. Here, we systematically integrate th ese processes in a reservoir-scale reactive transport model to quantify the impact of three distinct wettability conditions (water-wet, mixed-wet, and CO 2 -wet) on plume migration, mineral precipitation, pore-scale blockage, and long-term immobilization. Results show that under ideal conditions (no mineralization), water-wet systems exhibited plume aspect ratios up to 28% higher than CO 2 -wet systems and promoted at least 60% greater residual trapping. When mineralization was incorporated, mineral trapping accounted for 57%, 53%, and 46% of the total immobilized CO 2 in water-wet, mixed-wet, and CO 2 -wet systems, respectively. Dissolved CO 2 contributed 43-46%, while residual trapping remained below 5%, and mobile CO 2 was substantially reduced due to mineralization, demonstrating that as the system becomes more water-wet, CO 2 is more effectively converted into stable carbonate minerals, although extensive mineral precipitation may lead to pore clogging. A global Sobol-Morris sensitivity analysis further reveals that the kinetic parameters of reactive minerals are the most influential factor, followed by wettability, while salinity has only a minor influence. Furthermore, these findings highlight that wettability is not merely a flow property but a fundamental control on CO 2 mineralization in basalts, with direct implications for site selection for storage projects.

Basalts↗