SEARCH · Search NASA

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2022

Effect of oxidative breakers on organic matter degradation, contaminant mobility and critical mineral release during shale-fracturing fluid interactions in the Marcellus Shale

Oil and gas production from organic-rich shale formations has become viable through advancements in multistage hydraulic fracturing. However, fluctuations in oil and gas prices, coupled with the sharp decline in gas production after the initial days of fracturing operations, drive operators to devise strategies for enhancing hydrocarbon production. The use of highly reactive fracturing fluids that include strong oxidizing agents, also known as breakers, can potentially increase well productivity. The oxidative breakers are used to reduce the viscosity of gel-based fluids after the proppant is transported into fracture zones of the target formation. These breakers can also degrade the organic matter, enhancing the release of hydrocarbons. However, chemical byproducts generated by the interaction of oxidative breakers with the shale matrix have not been extensively studied. Here, this study investigated the fluid-rock interactions between Marcellus Shale and synthetic hydraulic fracturing fluid (HFF) solutions comprising three different oxidative breakers, i.e., ammonium persulfate, sodium bromate, and sodium hypochlorite, commonly used in the Appalachian Basin, USA. Our results demonstrate that the type of oxidizing breakers used in the HFF controlled the type and amount of volatile organic compounds (VOCs) and semi-volatile organic compounds (SVOCs) produced. In all HFF reacted effluents, we observed the transformation of VOCs and the presence of organic acids in variable amounts. However, effluents from HFF containing sodium bromate and sodium hypochlorite breakers showed the presence of several halogenated organic compounds. Changes in major ions and mineralogy indicate that carbonate dissolution and barite precipitation were ubiquitous in all shale reacted effluents. Our results also demonstrate that the addition of oxidative breakers increased the concentration of several major and trace elements in the effluents. These elements fall under the critical mineral (CM) or critical element category due to their high demand in emerging technologies and susceptibility to supply chain disruption due to variety of factors. However, before oxidative breakers can be used at a larger field scale to enhance the release and recovery of CM and hydrocarbons, a better understanding is required of the potential environmental impacts associated with the generation and transformation of contaminants during breaker-fluid-rock interaction.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2022

Effects of flow on uranium speciation in soils impacted by acidic waste fluids

Radioactive acidic liquid waste is a common byproduct of uranium (U) and plutonium (Pu) enrichment and recycling processes whose accidental and planned release has led to a significant input of U into soils and sediments across the world, including at the U.S. DOE's Hanford site (WA, USA). Because of the particularly hazardous nature of U, it is important to predict its speciation when introduced into soils and sediments by acidic waste fluids. Of fundamental importance are the coupled effects of acid-driven mineral transformation and reactive transport on U speciation. To evaluate the effect of waste-fluid residence time and co-associated dissolved phosphate concentrations on U speciation in impacted soils and sediments, uncontaminated surface materials (from the Hanford Site) were reacted with U-containing synthetic acidic waste fluids (pH 2) amended with dissolved phosphate concentrations in both batch (no flow) and flow-through column systems for 7–365 days. By comparing dissolved U behavior and solid phase speciation as a function of flow regimen, we found that the availability of proton-promoted dissolution products (such as Si) to sequester U into uranyl silicates was dependent on waste fluid-sediment contact time as uranyl silicates were not detected in short contact time flow-through systems but were detected in no-flow, long contact time, reactors. Moreover, the dominance of uranyl phosphate as neoprecipitate U scavenger (principally in the form of meta-ankoleite) in phosphate amended systems confirmed the importance of phosphate amendments for an efficient sequestration of U in the soils and sediments. Overall, our experiments suggest that the formation of uranyl silicates in soils impacted by acidic waste fluids is likely to be limited unless reaction products are allowed to accumulate in soil pores, highlighting the importance of investigating soil U speciation in flow-through, transport-driven systems as opposed to no-flow, batch systems. Finally, this study provides insights into uranium speciation and its potential changes under acidic conditions for better prediction of risks and subsequent development of efficient remediation strategies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Validation of SOLPS-ITER simulations with kinetic, fluid, and hybrid neutral models for JET-ILW low-confinement mode plasmas

For JET L-mode plasmas in low-recycling conditions (electron temperature at the outer strike point, T e,ot ≳ 30 eV ), SOLPS-ITER simulations agree within the error bars for the experimental profiles at the low-field side (LFS) divertor target. The peak Balmer-α (Dα ) emission in the LFS divertor agrees within the error bars of the KS3 filterscope diagnostic, but is approximately 30% lower than the peak value of the KT1 spectrometer. Simulations have been performed with fluid, kinetic, and hybrid models for the neutrals. The large fluid-kinetic discrepancies of more than a factor 2 are successfully corrected by using a hybrid fluid-kinetic approach, for which kinetic atoms are transferred to the fluid population when the local Knudsen number of the atom becomes smaller than a user-defined transition Knudsen number Kn t . The hybrid-kinetic discrepancies are limited to a few % for Kn t ≤ 100 . When increasing the upstream density to high-recycling conditions, at the onset of detachment ( T e,ot ≈ 5 eV ), the simulations predict more than a factor 2 lower peak ion saturation current to the LFS divertor than the experiments. Also the Dα emission is underpredicted with approximately a factor 2. For these high-recycling conditions, the fluid-kinetic discrepancies are limited to maximum 50%, which are again corrected by using the hybrid approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗

DOE OSTI2015

Nanoscopic dynamics of phospholipid in unilamellar vesicles: Effect of gel to fluid phase transition

Dynamics of phospholipids in unilamellar vesicles (ULV) is of interest in biology, medical, and food sciences since these molecules are widely used as biocompatible agents and a mimic of cell membrane systems. We have investigated the nanoscopic dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) phospholipid in ULV as a function of temperature using elastic and quasielastic neutron scattering (QENS). The dependence of the signal on the scattering momentum transfer, which is a critical advantage of neutron scattering techniques, allows the detailed analysis of the lipid motions that cannot be carried out by other means. In agreement with a differential scanning calorimetry measurement, a sharp rise in the elastic scattering intensity below ca. 296 K indicates a phase transition from the high-temperature fluid phase to the low-temperature solid gel phase. The microscopic lipid dynamics exhibits qualitative differences between the solid gel phase (in a measurement at 280 K) and the fluid phase (in a measurement at a physiological temperature of 310 K). The data analysis invariably shows the presence of two distinct motions: the whole lipid molecule motion within a monolayer, or lateral diffusion, and the relatively faster internal motion of the DMPC molecule. The lateral diffusion of the whole lipid molecule is found to be Fickian in character, whereas the internal lipid motions are of localized character, consistent with the structure of the vesicles. The lateral motion slows down by an order of magnitude in the solid gel phase, whereas for the internal motion not only the time scale, but also the character of the motion changes upon the phase transition. In the solid gel phase, the lipids are more ordered and undergo uniaxial rotational motion. However, in the fluid phase, the hydrogen atoms of the lipid tails undergo confined translation diffusion rather than uniaxial rotational diffusion. The localized translational diffusion of the hydrogen atoms of the lipid tails is a manifestation of the flexibility of the chains acquired in the fluid phase. Because of this flexibility, both the local diffusivity and the confinement volume for the hydrogen atoms increase linearly from near the lipid s polar head group to the end of its hydrophobic tail. Our results present a quantitative and detailed picture of the effect of the gel-fluid phase transition on the nanoscopic lipid dynamics in ULV. Lastly, the data analysis approach developed here has a potential for probing the dynamic response of lipids to the presence of additional cell membrane components.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Supercritical Fluid Adsorption to Weakly Attractive Solids: Universal Scaling Laws

In this study, adsorption of light hydrocarbons (C1–C3) confined in the pores of silica aerogel (nominal density 0.2 g/cm 3 ) was studied experimentally and by computer simulation. Five isotherms of ethane (25.0, 32.0, 32.3, 33.0, and 35.0 °C) were obtained over a temperature range including the critical temperature of the bulk fluid (T c = 32.18 °C) to ~220 bar by the vibrating tube densimetry (VTD) method, which yields unique information impossible to obtain by excess sorption measurements. Two isotherms (35 and 50 °C) of methane adsorption in SiO 2 aerogel were measured far above its critical temperature, T c = -82.59 °C, at 0–150 bar using a gravimetric method. The experimental pore fluid properties of methane, ethane, propane, and CO 2 (obtained earlier by VTD for the same SiO 2 aerogel sample) show behavior conforming to the principle of corresponding states over the entire covered reduced temperature ranges from T r = 0.976 to T r = 1.696. The magnitudes of the excess adsorption maxima and the reduced fluid densities where the maxima occur for C1–C3 hydrocarbons and CO 2 are similar functions of the reduced temperature. This behavior indicates that, at least under the conditions covered in this work, the differences in molecular size and shape have only minor impacts on the solid–fluid interactions underlying the adsorption behavior. Lastly, these findings are supported by large-scale GCMC simulations of a lattice gas confined in the pores of a computer representation of the SiO 2 aerogel framework, yielding excess adsorption isotherms in good agreement with experimental measurements, as well as the underlying microstructures of the adsorbed fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Olivine Dissolution in Simulated Lung and Gastric Fluid as an Analog to the Behavior of Lunar Particulate Matter Inside the Human Respiratory and Gastrointestinal Systems

With the Artemis III mission scheduled to land humans on the Moon in 2025, work must be done to understand the hazards lunar dust inhalation would pose to humans. In this study, San Carlos olivine was used as an analog of lunar olivine, a common component of lunar dust. Olivine was dissolved in a flow-through apparatus in both simulated lung fluid and 0.1 M HCl (simulated gastric fluid) over a period of approximately 2 weeks at physiological temperature, 37°C. Effluent samples were collected periodically and analyzed for pH, iron, silicon, and magnesium ion concentrations. The dissolution rate data derived from our measurements allow us to estimate that an inhaled 1.0 μm diameter olivine particle would take approximately 24 years to dissolve in the human lungs and approximately 3 weeks to dissolve in gastric fluid. Results revealed that inhaled olivine particles may generate the toxic chemical, hydroxyl radical, for up to 5–6 days in lung fluid. Olivine dissolved in 0.1 M HCl for 2 weeks transformed to an amorphous silica-rich solid plus the ferric iron oxy-hydroxide ferrihydrite. Olivine dissolved in simulated lung fluid shows no detectable change in composition or crystallinity. Equilibrium thermodynamic models indicate that olivine in the human lungs can precipitate secondary minerals with fibrous crystal structures that have the potential to induce detrimental health effects similar to asbestos exposure. Our work indicates that inhaled lunar dust containing olivine can settle in the human lungs for years and could induce long-term potential health effects like that of silicosis.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2024

Fluid-Driven Mass Transfer During Retrograde Metamorphism and Exhumation of the UHP Western Gneiss Region Terrane, Norway

Dehydration reactions within subducted oceanic crust are important for fluid-mediated element transfer within the subducting plate and potentially to the mantle wedge. The effects of metamorphic reactions and fluid flow on element recycling that occur during retrogression and exhumation of subducted continental crust from mantle depths are poorly understood. We study two metabasite pods with fresh eclogite cores and retrogressed amphibolite-facies rims and surrounding host gneiss within the Western Gneiss Region (WGR), Norway, to better understand element mobility and mass transfer during exhumation of subducted continental crust. Bulk-rock data were collected from samples taken across the pod and into the host gneiss. Phengite breakdown in eclogite and epidote recrystallization in veins and/or gneiss within pod cores contributed large ion lithophile elements and REE to retrogressed eclogite closest to the pod cores. In gneiss hosting the pods, phengite and epidote breakdown provided fluid that mediated elemental transfer and redistribution to the pod rim or tail. Compared to the studied pod in the southern WGR, the pod in the northern WGR underwent higher P-T conditions, partial melting and higher strain rates. This resulted in the infiltration of external fluid farther into the pod interior from the rim and facilitated larger mass gain in trace elements in the amphibolite tail of the pod relative to fresh eclogite in the core. The results show clear evidence for retrogression dehydration reactions driving significant fluid-mediated element redistribution as observed on the outcrop scale during exhumation following ultrahigh-pressure metamorphism, which directly impacts element signatures within the exhuming crust.

58 GEOSCIENCES↗

DOE OSTI2024

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

Recycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO 2 ) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ 34 S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H 2 O-NaCl solutions to calibrate a high-pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO 2 from FMQ-2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO 2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur-rich and high δ 34 S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34 S-rich subducted lithologies.

58 GEOSCIENCES↗

DOE OSTI2015

Diamond formation due to a pH drop during fluid–rock interactions

Diamond formation has typically been attributed to redox reactions during precipitation from fluids or magmas. Either the oxidation of methane or the reduction of carbon dioxide has been suggested, based on simplistic models of deep fluids consisting of mixtures of dissolved neutral gas molecules without consideration of aqueous ions. The role of pH changes associated with water–silicate rock interactions during diamond formation is unknown. Here we show that diamonds could form due to a drop in pH during water–rock interactions. We use a recent theoretical model of deep fluids that includes ions, to show that fluid can react irreversibly with eclogite at 900 °C and 5.0 GPa, generating diamond and secondary minerals due to a decrease in pH at almost constant oxygen fugacity. Overall, our results constitute a new quantitative theory of diamond formation as a consequence of the reaction of deep fluids with the rock types that they encounter during migration. Diamond can form in the deep Earth during water–rock interactions without changes in oxidation state.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Probing flow-induced nanostructure of complex fluids in arbitrary 2D flows using a fluidic four-roll mill (FFoRM)

Abstract Engineering flow processes to direct the microscopic structure of soft materials represents a growing area of materials research. In situ small-angle neutron scattering under flow (flow-SANS) is an attractive probe of fluid microstructure under simulated processing conditions, but current capabilities require many different sample environments to fully interrogate the deformations a fluid experiences in a realistic processing flow. Inspired by recent advances in microfluidics, we present a fluidic four-roll mill (FFoRM) capable of producing tunable 2D flow fields for in situ SANS measurements, that is intended to allow characterization of complex fluid nanostructure under arbitrary complex flows within a single sample environment. Computational fluid dynamics simulations are used to design a FFoRM that produces spatially homogeneous and sufficiently strong deformation fields. Particle tracking velocimetry experiments are then used to characterize the flows produced in the FFoRM for several classes of non-Newtonian fluids. Finally, a putative FFoRM-SANS workflow is demonstrated and validated through the characterization of flow-induced orientation in a semi-dilute cellulose nanocrystal dispersion under a range of 2D deformations. These novel experiments confirm that, for steady state straining flows at moderate strain rates, the nanocrystals orient along the principal strain-rate axis, in agreement with theories for rigid, rod-like Brownian particles in a homogeneous flow.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Gravity-driven controls on fluid and carbonate precipitation distributions in fractures

Abstract Many challenges related to carbon-dioxide ( $$\hbox {CO}_2$$ CO 2 ) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90 $$^\circ$$ ∘ (vertical fracture plane) to 30 $$^\circ$$ ∘ . The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0 $$^\circ$$ ∘ ). However, for fracture inclinations greater than 10 $$^\circ$$ ∘ , the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester $$\hbox {CO}_2$$ CO 2 through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Turbulent mixing of a slightly supercritical van der Waals fluid at low-Mach number

Supercritical fluids near the critical point are characterized by liquid-like densities and gas-like transport properties. These features are purposely exploited in different contexts ranging from natural products extraction/fractionation to aerospace propulsion. Large part of studies concerns this last context, focusing on the dynamics of supercritical fluids at high Mach number where compressibility and thermodynamics strictly interact. Despite the widespread use also at low Mach number, the turbulent mixing properties of slightly supercritical fluids have still not investigated in detail in this regime. This topic is addressed here by dealing with Direct Numerical Simulations of a coaxial jet of a slightly supercritical van der Waals fluid. Since acoustic effects are irrelevant in the low Mach number conditions found in many industrial applications, the numerical model is based on a suitable low-Mach number expansion of the governing equation. According to experimental observations, the weakly supercritical regime is characterized by the formation of finger-like structures – the so-called ligaments – in the shear layers separating the two streams. The mechanism of ligament formation at vanishing Mach number is extracted from the simulations and a detailed statistical characterization is provided. Ligaments always form whenever a high density contrast occurs, independently of real or perfect gas behaviors. The difference between real and perfect gas conditions is found in the ligament small-scale structure. More intense density gradients and thinner interfaces characterize the near critical fluid in comparison with the smoother behavior of the perfect gas. A phenomenological interpretation is here provided on the basis of the real gas thermodynamics properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Adapting SAFT-γ perturbation theory to site-based molecular dynamics simulation. II. Confined fluids and vapor-liquid interfaces

In this work, a new classical density functional theory is developed for group-contribution equations of state (EOS). Details of implementation are demonstrated for the recently-developed SAFT-γ WCA EOS and selective applications are studied for confined fluids and vapor-liquid interfaces. The acronym WCA (Weeks-Chandler-Andersen) refers to the characterization of the reference part of the third-order thermodynamic perturbation theory applied in formulating the EOS. SAFT-γ refers to the particular form of “statistical associating fluid theory” that is applied to the fused-sphere, heteronuclear, united-atom molecular models of interest. For the monomer term, the modified fundamental measure theory is extended to WCA-spheres. A new chain functional is also introduced for fused and soft heteronuclear chains. The attractive interactions are taken into account by considering the structure of the fluid, thus elevating the theory beyond the mean field approximation. The fluctuations of energy are also included via a non-local third-order perturbation theory. The theory includes resolution of the density profiles of individual groups such as CH{sub 2} and CH{sub 3} and satisfies stoichiometric constraints for the density profiles. New molecular simulations are conducted to demonstrate the accuracy of each Helmholtz free energy contribution in reproducing the microstructure of inhomogeneous systems at the united-atom level of coarse graining. At each stage, comparisons are made to assess where the present theory stands relative to the current state of the art for studying inhomogeneous fluids. Overall, it is shown that the characteristic features of real molecular fluids are captured both qualitatively and quantitatively. For example, the average pore density deviates ∼2% from simulation data for attractive pentadecane in a 2-nm slit pore. Another example is the surface tension of ethane/heptane mixture, which deviates ∼1% from simulation data while the theory reproduces the excess accumulation of ethane at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

An EQT-cDFT approach to determine thermodynamic properties of confined fluids

We present a continuum-based approach to predict the structure and thermodynamic properties of confined fluids at multiple length-scales, ranging from a few angstroms to macro-meters. The continuum approach is based on the empirical potential-based quasi-continuum theory (EQT) and classical density functional theory (cDFT). EQT is a simple and fast approach to predict inhomogeneous density and potential profiles of confined fluids. We use EQT potentials to construct a grand potential functional for cDFT. The EQT-cDFT-based grand potential can be used to predict various thermodynamic properties of confined fluids. In this work, we demonstrate the EQT-cDFT approach by simulating Lennard-Jones fluids, namely, methane and argon, confined inside slit-like channels of graphene. We show that the EQT-cDFT can accurately predict the structure and thermodynamic properties, such as density profiles, adsorption, local pressure tensor, surface tension, and solvation force, of confined fluids as compared to the molecular dynamics simulation results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Self-similarity and scaling transitions during rupture of thin free films of Newtonian fluids

Rupture of thin liquid films is crucial in many industrial applications and nature such as foam stability in oil-gas separation units, coating flows, polymer processing, and tear films in the eye. In some of these situations, a liquid film may have two free surfaces (referred to here as a free film or a sheet) as opposed to a film deposited on a solid substrate that has one free surface. The rupture of such a free film or a sheet of a Newtonian fluid is analyzed under the competing influences of inertia, viscous stress, van der Waals pressure, and capillary pressure by solving a system of spatially one-dimensional evolution equations for film thickness and lateral velocity. The dynamics close to the space-time singularity where the film ruptures is asymptotically self-similar and, therefore, the problem is also analyzed by reducing the transient partial differential evolution equations to a corresponding set of ordinary differential equations in similarity space. Here, for sheets with negligible inertia, it is shown that the dominant balance of forces involves solely viscous and van der Waals forces, with capillary force remaining negligible throughout the thinning process in a viscous regime. On the other hand, for a sheet of an inviscid fluid for which the effect of viscosity is negligible, it is shown that the dominant balance of forces is between inertial, capillary, and van der Waals forces as the film evolves towards rupture in an inertial regime. Real fluids, however, have finite viscosity. Hence, for real fluids, it is further shown that the viscous and the inertial regimes are only transitory and can only describe the initial thinning dynamics of highly viscous and slightly viscous sheets, respectively. Moreover, regardless of the fluid’s viscosity, it is shown that for sheets that initially thin in either of these two regimes, their dynamics transition to a late stage or final inertial-viscous regime in which inertial, viscous, and van der Waals forces balance each other while capillary force remains negligible, in accordance with the results of Vaynblat, Lister, and Witelski.

36 MATERIALS SCIENCE↗