Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗
DOE OSTI2019

Thicknesses of Chemically Altered Zones in Shale Matrices Resulting from Interactions with Hydraulic Fracturing Fluid

Hydraulic fracturing of unconventional shale reservoirs increases the fracture network surface area to access hydrocarbons from the low permeability rock matrix. Porosity and permeability of the matrix, through which hydrocarbons migrate to fractures, are important for determining production efficiency and can be altered by chemical interactions between shale and hydraulic fracturing fluids (HFFs). Here, we present results from an experimental study that characterizes the thickness of the alteration zone in the shale matrix after shale–HFF interactions. Experiments were conducted with whole cores submerged in HFF both with and without added barium and sulfate to promote barite scale formation. After 3 weeks of reaction at 77 bar and 80°C, the cores were characterized using X-ray microtomography, synchrotron X-ray fluorescence microprobe imaging, and synchrotron X-ray absorption spectroscopy. Our results show that the thickness of the altered zone depends on shale mineralogical composition and varies for different chemical reactions. Reactions between the low-carbonate Marcellus shale and HFF, pyrite (FeS 2 ) oxidation manifests as both a thick zone of sulfur oxidation (>0.5 cm) and a thinner zone of iron oxidation (100–150 $μ$m). Carbonate dissolution extended 100–200 μm into the matrix, with the resulting observable secondary porosity localized at the shale–fluid interface where mineral grains were removed by either dissolution or mechanical erosion. In solutions oversaturated with respect to barite, barite precipitates were observed in the reaction fluid and at the shal–HFF interface. In contrast, the carbonate dissolution zone in the high-carbonate Eagle Ford was only 30–40 μm thick, within which a uniform texture of increased porosity was observed. Pyrite oxidation in the Eagle Ford was evident from an iron oxidation zone (150–200 μm thick), while sulfur oxidation was minor and hard to observe. Barite precipitation extended 1–2 mm into the matrix when the initial HFF was oversaturated with respect to barite, filling shale microcracks down to the submicrometer length scale. Our findings provide a scientific basis to predict the extent of chemical alteration in shale reservoirs during hydraulic fracturing and its impacts on hydrocarbon production.

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

Chemical Speciation and Stability of Uranium in Unconventional Shales: Impact of Hydraulic Fracture Fluid

Uranium and other radionuclides are prominent in many unconventional oil/gas shales and is a potential contaminant in flowback/produced waters due to the large volumes/types of chemicals injected into the subsurface during stimulation. In this work, to understand the stability of U before and after stimulation, a geochemical study of U speciation was carried out on three shales (Marcellus, Green River, and Barnett). To do so, two types of samples for each shale were subjected to sequential chemical extractions: unreacted and shale-reacted with a synthetic hydraulic fracture fluid. Furthermore, a significant proportion of the total U (20–57%) was released from these three shales after reaction with fracture fluid, indicating that U is readily leachable. The total U released exceeds labile water-soluble and exchangeable fractions in unreacted samples, indicating that fluids leach more recalcitrant phases in the shale. Radiographic analysis of unreacted Marcellus shale thin sections shows U associated with detrital quartz and the clay matrix in the shale. Detrital zircon and TiO2 identified by an electron microprobe could account for the hot spots. This study shows that significant proportions of U in three shales are mobile upon stimulation. In addition, the extent of mobilization of U depends on the U species in these rocks.

58 GEOSCIENCES↗
DOE OSTI2017

Boundary Conditions for Diffusion-Mediated Processes within Linear Nanopores: Exact Treatment of Coupling to an Equilibrated External Fluid

In this paper, we consider a variety of diffusion-mediated processes occurring within linear nanopores, but which involve coupling to an equilibrated external fluid through adsorption and desorption. By determining adsorption and desorption rates through a set of tailored simulations, and by exploiting a spatial Markov property of the models, we develop a formulation for performing efficient pore-only simulations of these processes. Coupling to the external fluid is described exactly through appropriate nontrivial boundary conditions at the pore openings. This formalism is applied to analyze the following: (i) tracer counter permeation (TCP) where different labeled particles adsorb into opposite ends of the pore and establish a nonequilibrium steady state; (ii) tracer exchange (TE) with exchange of differently labeled particles within and outside the pore; (iii) catalytic conversion reactions where a reactant in the external fluid adsorbs into the pore and converts to a product which may desorb. The TCP analysis also generates a position-dependent generalized tracer diffusion coefficient, the form of which controls behavior in the TE and catalytic conversion processes. We focus on the regime of single-file diffusion within the pore which produces the strongest correlations and largest deviations from mean-field type behavior. Finally, behavior is quantified precisely via kinetic Monte Carlo simulations but is also captured with appropriate analytic treatments.

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

Continuum Model of Gas Uptake for Inhomogeneous Fluids

We describe a continuum model of gas uptake for inhomogeneous fluids (CMGIF) and use it to predict fluid adsorption in porous materials directly from gas-substrate interaction energies determined by first principles calculations or accurate effective force fields. The method uses a perturbation approach to correct bulk fluid interactions for local inhomogeneities caused by gas substrate interactions, and predicts local pressure and density of the adsorbed gas. The accuracy and limitations of the model are tested by comparison with the results of Grand Canonical Monte Carlo simulations of hydrogen uptake in metal-organic frameworks (MOFs). We show that the approach provides accurate predictions at room temperature and at low temperatures for less strongly interacting materials. As a result, the speed of the CMGIF method makes it a promising candidate for high-throughput materials discovery in connection with existing databases of nano-porous materials.

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

Flow Reduction in Pore Networks of Packed Silica Nanoparticles: Insights from Mesoscopic Fluid Models

A modified many-body dissipative particle dynamics (mDPD) model is rigorously calibrated to achieve realistic fluid–fluid/solid interphase properties and applied for mesoscale flow simulations to elucidate the transport mechanisms of heptane liquid and water, respectively, through pore networks formed by packed silica nanoparticles with a uniform diameter of 30 nm. Two million CPU core hours were used to complete the simulation studies. Results show reduction of permeability by 54–64% in heptane flow and by 88–91% in water flow, respectively, compared to the Kozeny–Carman equation. In these nanopores, a large portion of the fluids are in the near-wall regions and thus not mobile due to the confinement effect, resulting in reduced hydraulic conductivity. Moreover, intense oscillations in the calculated flow velocities also indicate the confinement effect that contests the external driven force to flow. Here, the generic form of Darcy’s law is considered valid for flow through homogeneous nanopore networks, while permeability depends collectively on pore size and surface wettability. This fluid-permeability dependency is unique to flow in nanopores. In addition, potential dependence of permeability on pore connectivity is observed when the porosity remains the same in different core specimens.

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

Effects of fluid shear stress on polyelectrolyte multilayers by neutron scattering studies

The structure of layer-by-layer (LbL) deposited nanofilm coatings consists of alternating polyethylenimine (PEI) and polystyrenesulfonate (PSS) films deposited on a single crystal quartz substrate. LbL-deposited nanofilms were investigated by neutron reflectomery (NR) in contact with water in the static and fluid shear stress conditions. The fluid shear stress was applied through a laminar flow of the liquid parallel to the quartz/polymer interface in a custom-built solid–liquid interface cell. The scattering length density profiles obtained from NR results of these polyelectrolyte multilayers (PEM), measured under different shear conditions, showed proportional decrease of volume fraction of water hydrating the polymers. For the highest shear rate applied (ca. 6800 s –1 ) the water volume fraction decreased by approximately 7%. The decrease of the volume fraction of water was homogeneous through the thickness of the film. Since there were not any significant changes in the total polymer thickness, it resulted in negative osmotic pressures in the film. The PEM films were compared with the behavior of thin films of thermoresponsive poly(N-isopropylacrylamide) (pNIPAM) deposited via spin-coating. The PEM and pNIPAM differ in their interactions with water molecules, and they showed opposite behaviors under the fluid shear stress. In both cases the polymer hydration was reversible upon the restoration of static conditions. Furthermore, a theoretical explanation is given to explain this difference in the effect of shear on hydration of polymeric thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Stability of Nd Hydroxyl Complexes at Near Neutral to Alkaline pH from 25 to 75 °C: Implications for Monazite Solubility in Hydrothermal Aqueous Fluids

Speciation of rare earth elements (REE) in aqueous fluids controls their mobilization during fluid-rock interaction. Thermodynamic modeling provides important insights into the factors controlling REE mobility and fractionation in aqueous fluids, in which aqueous REE hydroxyl complexes are particularly important across a broad pH range. However, the stability of these REE hydroxyl complexes remains poorly determined due to limited experimental data above ambient temperature. In this work, UV–vis spectrophotometry experiments were conducted to determine the cumulative hydrolysis constants for neodymium (Nd) from 25 to 75 °C in near-neutral to alkaline solutions. The color indicator m -cresol purple was used to determine in situ pH. Alkaline NaOH-bearing aqueous solutions were doped with varying initial Nd concentrations (0 to ∼0.155 mmol/kg), which resulted in the release of protons (H + ) and pH decrease during hydrolysis of Nd 3+ according to Nd 3+ + n H 2 O = Nd(OH) n 3– n + H + . The average OH – ligand number coordinated to Nd 3+ was found to increase from 1.0 to 1.6 at 25 °C and from 1.2 to 2.8 at 75 °C over a pH range from 6.3 to 9.0. The measured speciation shows an increased predominance of Nd(OH) 3 0 over the Nd(OH) 2 + and Nd(OH) 2+ species with increased temperature and pH. The increasing cumulative formation constants (log*β° n , n = 1 to 3) retrieved from 25 to 75 °C differ by 0.1 to 1.9 logarithmic units in comparison to existing thermodynamic databases. These updated thermodynamic data have important implications for geochemical modeling of the speciation of Nd hydroxyl complexes and the solubility of monazite as a function of pH.

58 GEOSCIENCES↗
DOE OSTI2019

Insights into a Greener Stimuli-Responsive Fracturing Fluid for Geothermal Applications

StimuFracTM, is a non-toxic stimuli-responsive fracturing fluid consisting on a CO2-reactive polymer which has shown at the lab-scale to consistently fracture rock cores at significantly lower net pressures in a range of representative geothermal pressure/temperature conditions. In this work, we report on the two main mechanisms responsible for fracturing rock at lower net pressures with this novel fracturing fluid. The first mechanism is the additional work / overpressure generated by aqueous polymer solution during its CO2-triggered volume expansion. The second mechanism is associated to a reduction in pore invasion pressure as a result of the lower interfacial tension of StimuFrac. This reduction in pore invasion pressure results in a reduction in fracturing net pressure as previously reported. The results presented in this work illustrate the potential of StimuFrac as a greener and less-energy intensive fracturing fluid for geothermal and fossil energy production.

Fernandez, Carlos A.↗
DOE OSTI2013

Characterizing Microbial Community and Geochemical Dynamics at Hydrothermal Vents Using Osmotically Driven Continuous Fluid Samplers

Microbes play a key role in mediating all aquatic biogeochemical cycles, and ongoing efforts are aimed at better understanding the relationships between microbial phylogenetic and physiological diversity, and habitat physical and chemical characteristics. Establishing such relationships is facilitated by sampling and studying microbiology and geochemistry at the appropriate spatial and temporal scales, to access information on the past and current environmental state that contributes to observed microbial abundances and activities. A modest number of sampling systems exist to date, few of which can be used in remote, harsh environments such as hydrothermal vents, where the ephemeral nature of venting underscores the necessity for higher resolution sampling. We have developed a robust, continuous fluid sampling system for co-registered microbial and biogeochemical analyses. The osmosis-powered bio-osmosampling system (BOSS) use no electricity, collects fluids with daily resolution or better, can be deployed in harsh, inaccessible environments and can sample fluids continuously for up to five years. Here we present a series of tests to examine DNA, RNA and protein stability over time, as well as material compatability, via lab experiments. We also conducted two field deployments at deep-sea hydrothermal vents to assess changes in microbial diversity and protein expression as a function of the physico-chemical environment. Our data reveal significant changes in microbial community composition co-occurring with relatively modest changes in the geochemistry. These data additionally provide new insights into the distribution of an enigmatic sulfur oxidizing symbiont in its free-living state. Data from the second deployment reveal differences in the representation of peptides over time, underscoring the utility of the BOSS in meta-proteomic studies. In concert, these data demonstrate the efficacy of this approach, and illustrate the value of using this method to study microbial and geochemical phenomena.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2019

Drift Instabilities in Thin Current Sheets Using a Two–Fluid Model With Pressure Tensor Effects

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. We use a two–fluid 10–moment model, which includes the pressure tensor and has been used to study reconnection, to study the drift kink and lower hybrid drift instabilities. Using a nonlocal linear eigenmode analysis, we find that for the kink mode, the 10–moment model shows good agreement with kinetic calculations with the same closure model used in reconnection simulations, while the electromagnetic and electrostatic lower hybrid instabilities require modeling the effects of the ion resonance using a Landau fluid closure. As a result, comparisons with kinetic simulations and the implications of the results for global magnetospheric simulations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The Seismic Response to Injected Carbon Dioxide: Comparing Observations to Estimates Based Upon Fluid Flow Modeling

Time-lapse seismic amplitude differences and travel time shifts, obtained while monitoring enhanced oil recovery at Cranfield, Mississippi, reveal coherent changes that are associated with the injection of carbon dioxide. Rock physics modeling highlights the importance of the oil, brine, and gas content of pore fluids prior to the injection of carbon dioxide. For example, compressional velocity changes due to the injection of carbon dioxide can drop from 300 m/s to less than 100 m/s as the percentage of oil increases from 1% to 50%. Predictions based upon a new technique for modeling wave propagation in a poroelastic medium containing an arbitrary number of fluids, coupled with multicomponent numerical reservoir modeling at Cranfield, reproduce the general pattern of observed seismic amplitude changes and travel time shifts. In particular, time-lapse amplitude changes suggest a significant and widespread lowering of compressional velocities due to the injection of CO 2 into an aquifer bounding the oil rim of the reservoir. It appears that the large-scale variations in preexisting pore fluid content have a major influence on seismic velocity changes, even in the highly heterogeneous reservoir at Cranfield.

58 GEOSCIENCES↗
DOE OSTI2019

Nonhysteretic Capillary Pressure in Two‐Fluid Porous Medium Systems: Definition, Evaluation, Validation, and Dynamics

Abstract A closure relation for capillary pressure plays an important role in the formulation of both traditional and evolving models of two‐fluid‐phase flow in porous medium systems. We review the traditional approaches to define capillary pressure, to describe it mathematically, to determine parameters for this relation, and to constrain the domain of applicability of this relation. In contrast to the traditional approach, we provide a rigorous, multiscale definition of capillary pressure, define the state domain of interest in practice, summarize computational and experimental approaches to investigate the system state, and apply the methods for two‐fluid states in a model ink bottle system, the classical Finney pack of spheres, and a synthetic sphere pack system. The results of these applications show that a state equation exists that describes capillary pressure without hysteresis. This state equation parameterizes a function that describes the nonwetting phase volume fraction in terms of the capillary pressure, the interfacial area, and the specific Euler characteristic of the nonwetting phase. Furthermore, this state equation applies over the complete range of conditions encountered in practice, and it applies under both equilibrium and dynamic conditions. This state equation involving capillary pressure forms an important foundation for the development of the next generation of macroscale two‐fluid‐phase flow models in porous medium systems.

Miller, C. T.↗
DOE OSTI2020

Microbial Community Composition in Deep‐Subsurface Reservoir Fluids Reveals Natural Interwell Connectivity

Abstract The identification of natural fractures and the wells they connect is crucial for the development of geological reservoirs because it may have an important impact on reservoir model construction and hydraulic fracture propagation. In this study we investigated the use of a novel data source, the microbial community composition in the reservoir formation fluids, for identification of interwell connectivity caused by natural fractures. We verified this concept at a newly developed mesoscale enhanced geothermal system (EGS) field testbed located 4,850 ft (1,478 m) beneath ground surface at the Sanford Underground Research Facility in Lead, SD. Fluids produced at or near the EGS testbed were sampled and subjected to high‐throughput 16S rRNA gene amplicon sequencing to analyze the microbial community profile therein. Despite the typically substantial heterogeneity across the community profiles of samples spatially distributed (10 m to 1.9 km apart) throughout the site, samples from two wells at the EGS testbed showed highly similar microbial community composition, suggesting the two wells intersected the same natural fracture. This evidence of natural connectivity between the two wells at the EGS testbed was later corroborated by core log analysis and sewer camera surveys into the boreholes. Besides the field case described in this study, microbial community composition as a reservoir diagnostic tool would be applicable in a much broader context such as unconventional hydrocarbon exploitation, groundwater reservoir characterization, and environmental remediation, adding valuable “hard” data capable of pinpointing the origins of fluids unambiguously.

Zhang, Yuran↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2025

Eocene-Oligocene Metamorphism, Fluid Flow and Deformation in the Ruby Mountains-East Humboldt Range Metamorphic Core Complex

Metamorphic core complexes are ubiquitous in collapsed orogens globally and play a primary role in crustal exhumation. Here, we investigate the metamorphic history of the Ruby Mountains-East Humboldt Range metamorphic core complex, Nevada, using petrochronology to understand how magmatism, metamorphism and deformation interact to modulate crustal rheology, and the timing of exhumation within the Sevier orogenic belt. Field study, microstructural analysis, thermobarometry and laser ablation split-stream monazite, titanite and allanite petrochronology of the mylonitic footwall were integrated to elucidate the pressure-temperature-time (P-T-t) evolution. Major-, trace-element and quartz-in-garnet thermobarometry show peak metamorphism occurred at 5.5–6 kbar and 600–650°C across the study area. Monazite and titanite U-(Th)-Pb petrochronology constrain this episode of metamorphism to 88–81 Ma with the structurally deepest samples yielding dates down to ca. 71 Ma. The preservation of peak metamorphic mineral assemblages indicates these rocks remained at near-peak temperature conditions into the Cenozoic. Cenozoic metamorphism occurred during a punctuated episode contemporaneous with the Great Basin ignimbrite flare up. All samples show a distinct 39–37 Ma population of monazite, titanite and allanite dates that correlate with the emplacement of quartz diorite and gabbro intrusions throughout the footwall. Accessory phases contain zoning textures in backscattered electron images that are consistent with fluid-mediated dissolution-reprecipitation reactions. The U-(Th)-Pb analyses from these zones yield 39–27 Ma dates interpreted as the timing of fluid flow. The results of this study elucidate the integrated importance of magmatism, melting and fluid flow in driving the coupled evolution of Late Cretaceous metamorphism within the Sevier hinterland and the subsequent Cenozoic metamorphic core complex development that exhumed the middle crust.

58 GEOSCIENCES↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
DOE OSTI2022

Transient evolution of permeability and friction in a slowly slipping fault activated by fluid pressurization

The mechanisms of permeability and friction evolution in a natural fault are investigated in situ. During three fluid injection experiments at different places in a fault zone, we measured simultaneously the fluid pressure, fault displacements and seismic activity. Changes in fault permeability and friction are then estimated concurrently. Results show that fault permeability increases up to 1.58 order of magnitude as a result of reducing effective normal stress and cumulative dilatant slip, and 19-to-60.8% of the enhancement occurs without seismic emissions. When modeling the fault displacement, we found that a rate-and-state friction and a permeability dependent on both slip and slip velocity together reasonably fit the fault-parallel and fault-normal displacements. This leads to the conclusion that the transient evolution of fault permeability and friction caused by a pressure perturbation exerts a potentially dominant control on fault stability during fluid flow.

58 GEOSCIENCES↗