Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical Fluids”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

AM Powder Flowability Capabilities at NASA Marshall Space Flight Center (MSFC)

Additive Manufacturing (AM) powder flowability is critical to metal 3D printing, because the more fluid the powder is, the better powder spreads. AM Powder Flowability Capabilities at NASA MSFC studies density, flow, particle size distribution, and morphology for programs such as Moon-to-Mars Planetary Autonomous Construction Technology (MMPACT), ASTM Proficiency Testing for AM and Powder Metallurgy, and MSFC’s AM team. Examining and measuring powder characteristics is essential for improving flowability of AM powder and ensuring lot-to-lot consistency, which will help to prevent defects in manufactured parts. The Contamination Control Team (CCT) uses several methods to characterize AM powder. Optical particle size distribution and morphology analyses determine particle parameters including size, circularity, convexity, and dimensions. The Carney and Hall Flowmeter Funnel measures time it takes powder to flow through a funnel to compare relative flowability (free-flowing and non-free-flowing). In addition to this, the device utilizes a density cup which determines apparent density. Lastly, the Revolution Powder Analyzer measures dynamic powder flowability and behavior over time via digital imaging. Using the CCT’s lab, the team has produced ASTM Powder Proficiency Testing results with flow rates, apparent density, particle size distribution, and tap density on titanium-based and nickel-based powder and compared data with other companies/labs. Morphology data was also collected for various regolith simulant powder for the MMPACT program. Lastly, the CCT used morphology data to examine the flowability of Inconel 718 powder for Selective Laser Melting (SLM) and Directed Energy Deposition (DED) printers at MSFC. The CCT’s essential work in studying and enabling AM powder characterization has helped and will continue to help study new and refined powders in the AM industry.

Additive Manufacturing↗
DOE OSTI2025

FLEX-FUEL MIXING CONTROLLED COMBUSTION ENABLED BY PRECHAMBER IGNITION

There is an imminent need to displace fossil diesel fuel with cleaner burning, domestically produced, renewable fuels for use in heavy-duty engines. Bioethanol is a prime candidate as it widely adopted in the U.S. as a gasoline additive ranging in volume percentage from 10% (E10) up to 85% (E85). Direct substitution of market available ethanol-gasoline blends for diesel fuel is not plausible as the stark reactivity differences would not constitute the same ignition quality nor achieve auto-ignition at all. This work focuses on the development of prechamber enabled mixing-controlled combustion (PC-MCC) as an advanced combustion strategy to facilitate reliable ignition and diffusion style combustion ethanol-gasoline fuel blends. PC-MCC involves integration of an actively fueled prechamber (PC) into a conventional compression ignition combustion system. When ignited, the PC ejects hot turbulent jets into the main combustion chamber that then interact with the direct injected fuel, prompting immediate ignition. The PC jet flames provide a robust thermal ignition source that allows the engine to operate agnostic of fuel composition, or flex-fuel. Computational fluid dynamics (CFD) modeling was used to assess critical design features of the PC while garnering insights into the ignition strategies that facilitate robust performance. A key finding was the ignition performance benefits of fuel-rich PC operation which yield exothermic jets. Based on the numerical findings, a prototype igniter was tested experimentally on both single and multi-cylinder engine platforms at a variety of operating conditions. The experimental results indicate flex-fuel PC-MCC is well capable of diesel-like combustion processes by demonstrating matched or improved gross thermal efficiencies and load variability within 2%. Fuel grade ethanol (E98) exhibited consistently lower NOx and immeasurable soot across the load space. E98 also demonstrated a significant improvement in thermal efficiency at light loads.

Zeman, Jared↗
NASA NTRS2004

Convective Instability in Ice I: Application to Callisto and Ganymede

Laboratory experiments measuring ice rheology suggest that it deforms under the influence of several nonNewtonian creep mechanisms, where the viscosity depends on both strain rate and temperature [Goldsby & Kohlstedt, 2001]. Whether or not a fluid with a purely temperature-dependent viscosity convects can be determined by comparing the Rayleigh number of the system to the critical Rayleigh number (Racr), which depends on rheological, thermal, and physical parameters of the fluid layer. However, in a nonNewtonian fluid where viscosity depends on the strain rate (i.e. velocity), convection can only occur if a temperature or velocity perturbation is issued to the system to lower the viscosity and permit fluid motions. Therefore, whether convection occurs in an ice I layer depends on initial conditions, in addition to rheological, thermal, and physical properties of the layer. We show new results for a scaling between the critical Rayleigh number and perturbation amplitude for grain boundary sliding rheology. This scaling can be used to determine the conditions required to initiate convection in the ice I shell of a generic icy satellite. We use this scaling to judge the convective instability of Ganymede and Callisto's ice shells in the absence of tidal dissipation.

Barr, A. C.↗
NASA NTRS1983

Studies of flows through N-sequential orifices

Critical mass flux and axial pressure profile data for fluid nitrogen are presented for N = 20, 15, 10, and 7 N-sequential orifice inlet configurations uniformly spaced at 15.5 cm. These data correlate well over a wide range in reduced temperature (0.7 T sub r,o ambient) and reduced pressure (to P sub r = 2) and are in general agreement with studies of one to four inlets and labyrinth seal gas data for 33 tooth equivalent sequential inlets. Experimental and theoretical agreement is good for liquid and gas critical mass flux.

Hendricks, R. C.↗
DOE OSTI2022

Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks

Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Physics-informed machine learning with differentiable programming for heterogeneous underground reservoir pressure management

Abstract Avoiding over-pressurization in subsurface reservoirs is critical for applications like CO $$_2$$ 2 sequestration and wastewater injection. Managing the pressures by controlling injection/extraction are challenging because of complex heterogeneity in the subsurface. The heterogeneity typically requires high-fidelity physics-based models to make predictions on CO $$_2$$ 2 fate. Furthermore, characterizing the heterogeneity accurately is fraught with parametric uncertainty. Accounting for both, heterogeneity and uncertainty, makes this a computationally-intensive problem challenging for current reservoir simulators. To tackle this, we use differentiable programming with a full-physics model and machine learning to determine the fluid extraction rates that prevent over-pressurization at critical reservoir locations. We use DPFEHM framework, which has trustworthy physics based on the standard two-point flux finite volume discretization and is also automatically differentiable like machine learning models. Our physics-informed machine learning framework uses convolutional neural networks to learn an appropriate extraction rate based on the permeability field. We also perform a hyperparameter search to improve the model’s accuracy. Training and testing scenarios are executed to evaluate the feasibility of using physics-informed machine learning to manage reservoir pressures. We constructed and tested a sufficiently accurate simulator that is 400 000 times faster than the underlying physics-based simulator, allowing for near real-time analysis and robust uncertainty quantification.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

Forced Convection Boiling and Critical Heat Flux of Ethanol in Electrically Heated Tube Tests

Electrically heated tube tests were conducted to characterize the critical heat flux (transition from nucleate to film boiling) of subcritical ethanol flowing at conditions relevant to the design of a regeneratively cooled rocket engine thrust chamber. The coolant was SDA-3C alcohol (95% ethyl alcohol, 5% isopropyl alcohol by weight), and tests were conducted over the following ranges of conditions: pressure from 144 to 703 psia, flow velocities from 9.7 to 77 ft/s, coolant subcooling from 33 to 362 F, and critical heat fluxes up to 8.7 BTU/in(exp 2)/sec. For the data taken near 200 psia, critical heat flux was correlated as a function of the product of velocity and fluid subcooling to within +/- 20%. For data taken at higher pressures, an additional pressure term is needed to correlate the critical heat flux. It was also shown that at the higher test pressures and/or flow rates, exceeding the critical heat flux did not result in wall burnout. This result may significantly increase the engine heat flux design envelope for higher pressure conditions.

Meyer, Michael L.↗
DOE OSTI2023

In-situ thermodynamics measurements at metal oxides-solution interfaces using Flow Adsorption Microcalorimetry.

Mineral-fluid interfaces are the principal sites of geochemical processes near Earth’s surface, hosting chemical reactions that play a fundamental role in (bio)geochemical cycles and in the fate and transport of anthropogenic contaminants, and that effectively control the compositions of soil and water environments. These complex interfaces are critical for our energy and environmental future. The mineral-fluid interface has been studied in an unprecedented level of detail with both experimental and computational approaches, separately and in combination. However, conspicuously missing from studies of the mineral-fluid interface are direct measurements of the energies of ion sorption and exchange. The literature on energetics and enthalpies of exchange, adsorption, dissolution, precipitation, and surface protonation reactions, especially those directly supported by experimental data, remains scarce despite their fundamental nature. The overarching goal of this project is to complete a systematic study of the thermodynamics properties of interfacial reactions at four MO surfaces (Rutile (α-TiO2), Quartz (SiO2), boehmite (γ-AlOOH) and goethite (α-FeOOH)) through the application and construction of novel flow adsorption microcalorimetry techniques and instrumentations. These unique and specialized microcalorimeters will operate at various temperatures and solution chemical compositions allowing for in-situ measurements across metal oxides and ligands of various characteristics. The overall research goal will be accomplished by completing the following three specific objectives (O): O1) Determine the energetics of surface protonation and deprotonation, ion exchange and ligand sorption reactions; O2) Investigate the surface charge thermodynamic properties under a range of temperature and solution chemical compositions; and O3) Develop predictive trends of the interplay between MO structure, surface coverage and surface reactivity. In addition to key thermodynamics parameters, calorimetric measurements provide a wealth of mechanistic information about reactions energetics and kinetics, surface charge characteristics, and structure-reactivity or selectivity relationships, all obtained in-situ and in real-time. This report includes science highlights from various projects completed over the performance period of the project. Also listed are dissemination opportunities, people supported on the grant and the impact on available physical resources and the discipline as a whole.

58 GEOSCIENCES↗
NASA NTRS1990

Critical point wetting drop tower experiment

The 100 m Drop Tower at NASA-Marshall was used to provide the step change in acceleration from 1.0 to 0.0005 g. An inter-fluid meniscus oscillates vertically within a cylindrical container when suddenly released from earth's gravity and taken into a microgravity environment. Oscillations damp out from energy dissipative mechanisms such as viscosity and interfacial friction. Damping of the oscillations by the later mechanism is affected by the nature of the interfacial junction between the fluid-fluid interface and the container wall. In earlier stages of the project, the meniscus shape which developed during microgravity conditions was applied to evaluations of wetting phenomena near the critical temperature. Variations in equilibrium contact angle against the container wall were expected to occur under critical wetting conditions. However, it became apparent that the meaningful phenomenon was the damping of interfacial oscillations. This latter concept makes up the bulk of this report. Perfluoromethyl cyclohexane and isopropanol in glass were the materials used for the experiment. The wetting condition of the fluids against the wall changes at the critical wetting transition temperature. This change in wetting causes a change in the damping characteristics of the interfacial excursions during oscillation and no measurable change in contact angle. The effect of contact line friction measured above and below the wetting transition temperature was to increase the period of vertical oscillation for the vapor-liquid interface when below the wetting transition temperature.

Kaukler, William F.↗
DOE OSTI2023

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
NASA NTRS1997

Noise and Dynamical Pattern Selection in Solidification

The overall goal of this project was to understand in more detail how a pattern-forming system can adjust its spacing. "Pattern-forming systems," in this context, are nonequilibrium contina whose state is determined by experimentally adjustable control parameter. Below some critical value of the control system then has available to it a range of linearly stable, spatially periodic steady states, each characterized by a spacing which can lie anywhere within some band of values. These systems like directional solidification, where the solidification front is planar when the ratio of growth velocity to thermal gradient is below its critical value, but takes on a cellular shape above critical. They also include systems without interfaces, such as Benard convection, where it is the fluid velocity field which changes from zero to something spatially periodic as the control parameter is increased through its critical value. The basic question to be addressed was that of how the system chooses one of its myriad possible spacings when the control parameter is above critical, and in particular the role of noise in the selection process. Previous work on explosive crystallization had suggested that one spacing in the range should be preferred, in the sense that weak noise should eventually drive the system to that spacing. That work had also suggested a heuristic argument for identifying the preferred spacing. The project had three main objectives: to understand in more detail how a pattern-forming system can adjust its spacing; to investigate how noise drives a system to its preferred spacing; and to extend the heuristic argument for a preferred spacing in explosive crystallization to other pattern-forming systems.

Kurtze, Douglas A.↗
DOE OSTI2026

Mobility of heavy rare earth elements in magmatic-hydrothermal systems: Experimental determination of DyPO 4 solubility in supercritical fluids

Dysprosium (Dy) is part of the heavy rare earth elements (REE) and a critical component in the production of high-performance super magnets used in wind turbines, electric vehicles, and other green technologies. The HREE phosphate, xenotime, is commonly associated with hydrothermal systems, where the mobility, fractionation, and deposition of REE are controlled by the stability of aqueous complexes and solubility of REE minerals which depend on temperature, pressure, and fluid composition (salinity, pH, and ligand availability). Here, we conducted a series of synthetic DyPO 4 solubility experiments in NaCl-HCl-bearing aqueous solutions (pH of 2 at 25 °C) from 500 to 700 °C and 1.5 kbar using externally heated pressure vessels. DyPO 4 displays a strong prograde solubility with temperature across all studied fluid salinities (0–0.5 m NaCl). The DyPO₄ solubility isotherms at 500 °C, 600 °C, and 700 °C show a complex relationship with salinity where DyPO₄ solubility is highest (39.2–1552 ppm at 500–700 °C) in 0 m NaCl solutions, decreases strongly (21–331 ppm at 500–700 °C) in the 0.06 m NaCl solutions and remains relatively constant in the 0.25 m (42.8–248 ppm) and 0.5 m NaCl (40.2–124 ppm) solutions. Thermodynamic modeling suggests that pH increases as a function of temperature and salinity, and controls solubility due to the stability of chloride complexes (i.e., DyCl +2 and DyCl 2 + ) at 500–600 °C, transitioning to hydroxyl complexes (i.e., Dy(OH) 3 0 ) dominating at 700 °C and relatively alkaline pH. Comparison with NdPO 4 solubility data indicates that the light REE in monazite is significantly more soluble than the heavy REE in xenotime in saline fluids at high temperature (500–700 °C), implying efficient light/heavy REE fractionation in supercritical fluids. The lower solubility of DyPO 4 relative to NdPO 4 suggests that at high temperatures, acidic (pH 2–4), low-salinity (∼5 wt% NaCl) fluids will favor the precipitation of xenotime over monazite as it will cool from 700 °C to below 500 °C. These observations align with paragenetic mineral assemblages in NYF pegmatites (e.g., Baveno, Italy), where xenotime commonly immobilizes heavy REE during the high-temperature hydrothermal stages.

58 GEOSCIENCES↗
NASA NTRS1976

Stability of flow of a thermoviscoelastic fluid between rotating coaxial circular cylinders

The stability problem of thermoviscoelastic fluid flow between rotating coaxial cylinders is investigated using nonlinear thermoviscoelastic constitutive equations due to Eringen and Koh. The velocity field is found to be identical with that of the classical viscous case and the case of the viscoelastic fluid, but the temperature and pressure fields are found to be different. By imposing some physically reasonable mechanical and geometrical restrictions on the flow, and by a suitable mathematical analysis, the problem is reduced to a characteristic value problem. The resulting problem is solved and stability criteria are obtained in terms of critical Taylor numbers. In general, it is found that thermoviscoelastic fluids are more stable than classical viscous fluids and viscoinelastic fluids under similar conditions.

Ghandour, N. N.↗
NASA NTRS2015

Resource Prospector Propulsion System Cold Flow Testing

Resource Prospector (RP) is a NASA mission being led by NASA Ames Research Center with current plans to deliver a scientific payload package aboard a rover to the lunar surface. As part of an early risk reduction activity, Marshall Space Flight Center (MSFC) and Johnson Space Flight Center (JSC) have jointly developed a government-version concept of a lunar lander for the mission. The spacecraft consists of two parts, the lander and the rover which carries the scientific instruments. The lander holds the rover during launch, cruise, and landing on the surface. Following terminal descent and landing the lander portion of the spacecraft become dormant after the rover embarks on the science mission. The lander will be equipped with a propulsion system for lunar descent and landing, as well as trajectory correction and attitude control maneuvers during transit to the moon. Hypergolic propellants monomethyl hydrazine and nitrogen tetroxide will be used to fuel sixteen 70-lbf descent thrusters and twelve 5-lbf attitude control thrusters. A total of four metal-diaphragm tanks, two per propellant, will be used along with a high-pressure composite-overwrapped pressure vessel for the helium pressurant gas. Many of the major propulsion system components are heritage missile hardware obtained by NASA from the Air Force. In parallel with the flight system design activities, a simulated propulsion system based on flight drawings was built for conducting a series of water flow tests to characterize the transient fluid flow of the propulsion system feed lines and to verify the critical operation modes such as system priming, waterhammer, and crucial mission duty cycles. The primary objective of the cold flow testing was to simulate the RP propulsion system fluid flow operation through water flow testing and to obtain data for anchoring analytical models. The models will be used to predict the transient and steady state flow behaviors in the actual flight operations. All design and build efforts, including the analytical modeling, have been performed. The cold flow testing of the propulsion system was set up and conducted at a NASA MSFC test facility. All testing was completed in the summer of 2014, and this paper documents the results of that testing and the associated fluid system modeling efforts.

Williams, Hunter↗
NASA NTRS1990

Numerical simulation of a catastrophe model for coronal mass ejections

Van Tend and Kuperus (1978) proposed a simple catastrophe model for magnetically driving coronal mass ejections, prominence eruptions, and two-ribbon flares. Their model, which is based on simple circuit concepts, suggests that a stable configuration containing a current filament will lose equilibrium when the filament current exceeds a critical value. Here, a two-dimensional numerical simulation is used to test how the Van Tend-Kuperus model works in an ideal MHD fluid. The simulation exhibits the expected loss of mechanical equilibrium near the predicted critical value, but the current filament moves only a short distance upward before coming to rest at a new equilibrium. However, this new equilibrium contains a current sheet which is resistively unstable to magnetic reconnection, and if magnetic reconnection occurs rapidly, the filament can continue to move upward at Alfvenic speeds.

Forbes, T. G.↗
NASA NTRS1997

The Neuroendocrinology of Thirst and Salt Appetite: Visceral Sensory Signals and Mechanisms of Central Integration

This review examines recent advances in the study of the behavioral responses to deficits of body water and body sodium that in humans are accompanied by the sensations of thirst and salt appetite. Thirst and salt appetite are satisfied by ingesting water and salty substances. These behavioral responses to losses of body fluids, together with reflex endocrine and neural responses, are critical for reestablishing homeostasis. Like their endocrine and neural counterparts, these behaviors are under the control of both excitatory and inhibitory influences arising from changes in osmolality, endocrine factors such as angiotensin and aldosterone, and neural signals from low and high pressure baroreceptors. The excitatory and inhibitory influences reaching the brain require the integrative capacity of a neural network which includes the structures of the lamina terminalis, the amygdala, the perifornical area, and the paraventricular nucleus in the forebrain, and the lateral parabrachial nucleus (LPBN), the nucleus tractus solitarius (NTS), and the area postrema in the hindbrain. These regions are discussed in terms of their roles in receiving afferent sensory input and in processing information related to hydromineral balance. Osmoreceptors controlling thirst are located in systemic Viscera and in central structures that lack the blood-brain barrier. Angiotensin and aldosterone act on and through structures of the lamina terminalis and the amygdala to stimulate thirst and sodium appetite under conditions of hypovolemia. The NTS and LPBN receive neural signals from baroreceptors and are responsible for inhibiting the ingestion of fluids under conditions of increased volume and pressure and for stimulating thirst under conditions of bypovolemia and hypotension. The interplay of multiple facilitory influences within the brain may take the form of interactions between descending angiotensinergic systems originating in the forebrain and ascending adrenergic systems emanating from the hindbrain. Oxytocin and serotonin are additional candidate neuro- chemicals with postulated inhibitory central actions and with essential roles in the overall integration of sensory input within the neural network devoted to maintaining hydromineral balance.

Johnson, Alan Kim↗
DOE OSTI2025

The Fluid Dynamics Uncertainty Quantification Challenge Problem: XFOIL vs. MFOIL

Uncertainty quantification (UQ) has become more critical in aerospace engineering due to the growing dependence on computational tools for design optimization and performance analyses of aerospace vehicles. Even though the significance of UQ in assessing the credibility of computational analyses is well recognized, its costs and complexity impede its integration into standard practices, particularly in computational fluid dynamics (CFD) and other fluid analyses. This paper presents a UQ study for low-fidelity computational aerodynamics analyses with XFOIL and mfoil (i.e., the MATLAB version of XFOIL with several implementation modifications); these tools are utilized widely in both research and education. The main contributions of this paper are as follows: 1) improved precision in quantifying the uncertainty of the baseline Monte Carlo results used to benchmark surrogate modeling techniques for UQ, 2) quantification of the effect of the implementation differences between XFOIL and mfoil on solution quantities of interest (QoIs), such as lift and pitching moment coefficients, and 3) development of an open-source UQ library for use with XFOIL and mfoil, which has educational values and helps promote UQ for fluid analyses with aerospace applications. Results and discussions revolve around cases 1-4 of the challenge problem posed by the AIAA Fluid Dynamics Technical Committee’s Uncertainty Quantification Discussion Group (UQDG). In case 3, this work employs CFDverify, an open-source solution verification software, to quantify the discretization error and evaluate the extrapolated QoIs based on the grid convergence index (GCI). This UQ study differentiates itself from previous studies in the rigor of handling baseline Monte Carlo uncertainty and in including mfoil, which is a more accessible alternative to XFOIL. Finally, despite the growing computing power, low-fidelity computational tools remain valuable, such as for aerodynamic shape optimization at Mach numbers below 0.65 and low-to-mid Reynolds numbers.

Lay, Aidan S [University of Tennessee, Knoxville (↗
DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗