Search NASA⌕ Search

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.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-equilibrium molecular dynamics simulation to evaluate the effect of confinement on fluid flow in silica nanopores

The flow of fluids in nano-confinement has applications in separations, water purification, medical systems and in the recovery of fluids in petroleum systems. It is believed that fluids do not obey continuum laws when flowing in nanopores (or in confinement). This paper attempts to shed light on various nano confinement effects such as fluid-wall interactions, pore size and molecular geometry, using molecular dynamic simulation. Here, water, hexane, and methanol flow behaviors were simulated for pore diameters ranging from 1 to 8 nm. In addition to the density analysis of the confined fluids, the fluid flows in saturated nanopores were simulated by the sectional flow method. Water and methanol molecules were completely stabilized in the 1 nm pore. Hexane molecules could not enter the 1 nm pore, due to geometric considerations. For the 2 – 8 nm pores, all fluids showed reduced flow rate compared to the Hagen–Poiseuille flow, due to an interfacial molecular layer stabilized at the pore surface. Flow reduction observed in these studies is contrary to significant flow enhancements observed for fluid flow in carbon nanotubes of similar dimensions. Water flow showed almost constant stick length for all the pore sizes. Methanol flow had the largest stick length. Hexane flow was reduced because of overcrowding of molecules at the pore surface. The effect of confinement diminished with an increase in pore diameter.

42 ENGINEERING↗
DOE OSTI2023

Effect of initial water saturation on the performance of fracturing fluids with and without polyallylamine under simulated EGS conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO 2 - reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO 2 , can be used as a potentially less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO 2 fracturing fluids outperformed other fluids such as water, CO 2 , and CO 2 /water in generating large fractures with less fluid consumed. The objective of this work is to investigate the effect of initial water saturation on the performance of StimuFrac by conducting hydraulic fracturing tests with ½ foot cubic rock samples held under representative EGS stress/temperature conditions and by using cyclic injection strategies (under constant injection rate). The resulting fracture hydraulic conductivities, breakdown pressures, and volumes of fluids required are compared. Methods/Procedures/Process: To simulate geothermal reservoir conditions, in all tests, the rock sample was held under triaxial confinement and at 200 °C, and different volumes of water were initially injected into the rock sample before any fracturing processes were initiated. For the single-cycle PAA (or water) alternating CO 2 (PAG or WAG) injection fracturing experiments, one complete cycle consisted of two steps: (1) injecting a PAA slug (or water slug) followed by (2) injecting CO 2 to initiate and propagate the fracture. For experiments involving multiple injection cycles, the CO 2 injection pressure is increased until it peaks and begins to decline (indicating fracture initiation at this moment), and then continued being injected for another 30 seconds to propagate the fracture. Then, these two-step cycles [injection of PAA (or water) followed by CO 2 injection (up to 2-4 mL/min)] are repeated. Applications/Significance/Novelty: The results of this study suggest that water saturation significantly affects the fracturing fluid transmission into the rock pore space, thus affecting the fracture initiation and propagation. In this study, fracturing tests via a single injection cycle or multiple injection cycles were performed. Splitting the rock samples in half after testing reveals that fracture propagation is significantly limited under high water saturation conditions (three-day initial water injection) compared to stimulation experiments performed in hot dry rock. The fractures propagate less than 1/3 of the distance from the wellbore to the outer rock surface, and in some cases, no fracture is generated. This may be caused by leak-off dominating the fracturing process and the fluid injection rate is insufficient to overcome leak-off, even under high injection rate conditions. Additionally, CO 2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate sufficiently high-pressure gradients away from the near-wellbore region. Under low water saturation conditions (dry rock or after 1-day initial water injection), PAA/CO 2 consistently generated significantly larger fractures compared with the other fluids. CO 2 generated large fractures only in the hot dry rock and only when using high injection rates, though data variability is high.

58 GEOSCIENCES↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part I

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In conclusion, these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II

Here, in this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes—wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries—and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid–structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier–Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau–Rayleigh instability in water ripple free surface flow. In these last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Shale-fluid interactions during forced lmbibition and flow-back

One way to minimize the need for disposal of flow-back waters produced during the hydraulic fracturing of unconventional reservoirs is, potentially, to recycle these fluids by recovering and re-injecting them during hydraulic fracturing of nearby new wells. These flow-back waters contain, however, large concentrations of dissolved salts. During hydraulic fracturing of unconventional reservoirs, a well may undergo extended shut-in periods, following the initial hydraulic-fracture stimulation. The concern here, therefore, is that the salts dissolved in the re-injected flow-back waters may precipitate in the subsurface, thus promoting scale formation; or alternately, because of their high salinity, these same fluids may cause the dissolution of various minerals present in the subsurface. Large divalent cations (e.g., Ca²⁺, Ba²⁺ and Sr²⁺) and anions (e.g., CO3²⁻ and SO4²) may chemically interact with and cause the swelling of clay fractions, thus blocking pores in certain pore size ranges. All these phenomena (i.e., scale formation, clay swelling, minerals dissolution) can influence the petrophysical properties of the formation, including its permeability, porosity, contact angle, etc. One way to potentially overcome some of these technical challenges, is to pre-treat these flow-back waters, prior to re-injection (e.g., via chemical precipitation, followed by membrane filtration) in order to remove unwanted salts. However, such pre-treatment is costly, and the degree of salt removal needed to prevent formation damage is important in determining the economic viability of flow-back water recycle via re-injection. A key focus of this work, therefore, is to study the impact that the injection of flow-back waters, that have undergone varying degrees of pre-treatment and clean-up, will have on the flow-through characteristics of a tight-gas reservoir, particularly of the Marcellus Shale formation. To do so, we have carried-out forced-imbibition experiments with shale samples from various depths in the Marcellus formation, whose properties span high to low ranges of TOC, clay content and matrix permeabilities. We have studied the forced-imbibition characteristics, at an injection pressure of 24,233 kPa (3500 psig), of various flow-back fluids into shale cores saturated with methane at a pore pressure of 17,338 kPa (2500 psig) and a confining stress of 27,680 kPa (4000 psig). Following the forced imbibition experiments, we analyzed the cores for potential changes in permeability and porosity. Quantification of the ion content of the injected and recovered fluids was also performed using Ion Chromatography (IC). The IC analysis of the recovered fluids indicates substantial interactions and geochemical reactions between the invading fluids and the shale cores. The injection of flow-back fluids shows significant impact on the flow-through characteristics, with untreated flow-back fluids, which have the highest cation and anion content, also exhibiting the highest observable reductions in permeability (>90%).

defects, mechanical behavior, carbon sequestration↗
DOE OSTI2022

Heat pulse testing at monitoring wells to estimate subsurface fluid velocities in geological CO 2 storage

Monitoring the injected CO 2 during geological CO 2 storage (GCS) is essential to assure containment and identify CO 2 leakage. Here in this work, a new approach is introduced to estimate the evolution of the downhole fluid velocity at a monitoring well and identify CO 2 arrival time using in-well heat pulse/tracer test. The proposed technique involves using a downhole heater to generate a series of heat pulses and measuring their corresponding temperature response. The surface temperature of the downhole heater is controlled by the supplied electrical power and the heat loss by convection to the surroundings. Convective heat transfer is well described using Newton's law of cooling in which the temperature difference between the heater and the surrounding fluids drives the heat transfer, for which the convection heat transfer coefficient (h) controls the magnitude of heat loss. Among various factors that control h, it depends on the type of the flowing fluid and its velocity. Through analyzing the measured temperature at different heat pulses, the changes in h - due to mobilization of the in-situ brine or CO 2 arrival - can be estimated. Consequently, the velocity of the flowing fluid across the heater can be obtained. Since heat transfer by convection is sensitive to the type of the surrounding fluid, intrusion of CO 2 can be detected from the relatively higher surface temperature obtained at CO 2 arrival. Churchill and Bernstein (1977)'s correlation is adopted to estimate the change of fluid velocity in terms of the change in h. To demonstrate the validity of the proposed technique, the results are applied and validated against those of COMSOL Multiphysics simulation tool for single-phase brine (before CO 2 arrival) and single-phase CO 2 (after CO 2 arrival). The observed temperature heating is sensitive to the flowing fluid velocity and fluid type. The temperature signal observed at CO 2 arrival is large and easily detectable using temperature monitoring tool which provides reliable indication for tracking CO 2 arrival at monitoring wells compared with passive temperature monitoring. The results obtained using the proposed technique agree very well with the numerical results obtained from the simulation tool with a maximum estimation error of 7 percent.

02 PETROLEUM↗
DOE OSTI2015

Surface-texture evolution of different chemical-vapor-deposited zinc sulfide flats polished with various magnetorheological fluids

The macro-structure of chemical-vapor-deposited (CVD) zinc sulfide (ZnS) substrates is characterizedby cone-like structures that start growing at the early stages of deposition. As deposition progresses,these cones grow larger and reach centimeter size in height and millimeter size in width. It is challengingto polish out these features from the top layer, particularly for the magnetorheological finishing (MRF)process. A conventional MR fluid tends to leave submillimeter surface artifacts on the finished surface,which is a direct result of the cone-like structure. Here we describe the MRF process of polishing four CVD ZnS substrates, manufactured by four differentvendors, with conventional MR fluid at pH 10 and zirconia-coated-CI (carbonyl iron) MR fluids at pH 4, 5,and 6. We report on the surface–texture evolution of the substrates as they were MRF polished with thedifferent fluids. We show that performances of the zirconia-coated-CI MR fluid at pH 4 are significantlyhigher than that of the same fluid at pH levels of 5 and 6 and moderately higher than that of a conventionalMR fluid at pH 10. An improvement in surface–texture variability from part to part was also observedwith the pH 4 MR fluid.

CVD ZnS↗
DOE OSTI2016

Reduction of hexavalent chromium by fasted and fed human gastric fluid. I. Chemical reduction and mitigation of mutagenicity

Evaluation of the reducing capacity of human gastric fluid from healthy individuals, under fasted and fed conditions, is critical for assessing the cancer hazard posed by ingested hexavalent chromium [Cr(VI)] and for developing quantitative physiologically-based pharmacokinetic models used in risk assessment. In the present study, the patterns of Cr(VI) reduction were evaluated in 16 paired pre- and post-meal gastric fluid samples collected from 8 healthy volunteers. Human gastric fluid was effective both in reducing Cr(VI), as measured by using the s-diphenylcarbazide colorimetric method, and in attenuating mutagenicity in the Ames test. The mean (± SE) Cr(VI)-reducing ability of post-meal samples (20.4 ± 2.6 μg Cr(VI)/mL gastric fluid) was significantly higher than that of pre-meal samples (10.2 ± 2.3 μg Cr(VI)/mL gastric fluid). When using the mutagenicity assay, the decrease of mutagenicity produced by pre-meal and post-meal samples corresponded to reduction of 13.3 ± 1.9 and 25.6 ± 2.8 μg Cr(VI)/mL gastric fluid, respectively. These data are comparable to parallel results conducted by using speciated isotope dilution mass spectrometry. Cr(VI) reduction was rapid, with > 70% of total reduction occurring within 1 min and 98% of reduction is achieved within 30 min with post-meal gastric fluid at pH 2.0. pH dependence was observed with decreasing Cr(VI) reducing capacity at higher pH. Attenuation of the mutagenic response is consistent with the lack of DNA damage observed in the gastrointestinal tract of rodents following administration of ≤ 180 ppm Cr(VI) for up to 90 days in drinking water. Quantifying Cr(VI) reduction kinetics in the human gastrointestinal tract is necessary for assessing the potential hazards posed by Cr(VI) in drinking water. - Highlights: • Cr(VI) reduction capacity was greater in post-meal than paired pre-meal samples. • Cr(VI) reduction was rapid, pH dependent, and due to heat stable components. • Gastric fluid attenuates the mutagenicity of Cr(VI) and reduces risk of GI cancer. • The results are useful to quantitatively evaluate the fate of ingested Cr(VI).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Designing CO 2 -Responsive Multifunctional Nanoscale Fluids with Tunable Hydrogel Behavior for Subsurface Energy Recovery

The ability to direct subsurface fluid flow or achieve a high level of control on permeability in subsurface environments necessitates the development of tunable novel subsurface fluids. These multifunctional fluids should have the potential to form hydrogels to divert flow or enhance fracture networks at elevated pressures, ability to transport proppants, undergo reversible transformations from gel-like to fluid-like in response to chemical- or pressure-based perturbations, and have improved CO 2 carrying capacity for enhancing the miscibility and flowability of oil and gas. In this paper, we discuss the development of multifunctional nanofluids constructed from silica (SiO 2 ) nanoparticles and poly(allylamine) (PAA), amine bearing polymer chains with high affinity for CO 2 . A 2-fold increase in CO 2 absorption in SiO 2 -PAA nanofluids compared to the pure polymer was noted. Upon CO 2 absorption, the weakly interacting polymeric chains around the nanoparticles formed relatively compact hydrogels. Time-resolved ultrasmall-angle X-ray scattering/small-angle X-ray scattering measurements showed the transition from swollen branched polymers to Gaussian coils with an increased exposure to CO 2 . Further, CO 2 -induced hydrogel formation in aqueous fluids bearing 1 wt % SiO 2 -PAA nanofluids occurred at room temperature, unlike in fluids bearing 1 wt % PAA. These observations point to the feasibility of forming hydrogels at lower temperatures and pressures using novel nanofluids as opposed to using the pure polymer. As a result, the ability to tune the structures and morphologies of these fluids expands the potential applications of these fluids in a wide range of subsurface environments.

03 NATURAL GAS↗
DOE OSTI2018

Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation, and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR (Incompressible Adaptive Mesh Refinement) code, which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

42 ENGINEERING↗
DOE OSTI2022

Deep learning-based quasi-continuum theory for structure of confined fluids

Predicting the structural properties of water and simple fluids confined in nanometer scale pores and channels is essential in, for example, energy storage and biomolecular systems. Classical continuum theories fail to accurately capture the interfacial structure of fluids. In this work, we develop a deep learning-based quasi-continuum theory (DL-QT) to predict the concentration and potential profiles of a Lennard-Jones (LJ) fluid and water confined in a nanochannel. The deep learning model is built based on a convolutional encoder–decoder network (CED) and is applied for high-dimensional surrogate modeling to relate the fluid properties to the fluid–fluid potential. The CED model is then combined with the interatomic potential-based continuum theory to determine the concentration profiles of a confined LJ fluid and confined water. Further, we show that the DL-QT model exhibits robust predictive performance for a confined LJ fluid under various thermodynamic states and for water confined in a nanochannel of different widths. The DL-QT model seamlessly connects molecular physics at the nanoscale with continuum theory by using a deep learning model.

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

Wetting transition and fluid trapping in a microfluidic fracture

Immiscible fluid–fluid displacement in confined geometries is a fundamental process occurring in many natural phenomena and technological applications, from geological CO 2 sequestration to microfluidics. Due to the interactions between the fluids and the solid walls, fluid invasion undergoes a wetting transition from complete displacement at low displacement rates to leaving a film of the defending fluid on the confining surfaces at high displacement rates. While most real surfaces are rough, fundamental questions remain about the type of fluid–fluid displacement that can emerge in a confined, rough geometry. Here, we study immiscible displacement in a microfluidic device with a precisely controlled structured surface as an analogue for a rough fracture. We analyze the influence of the degree of surface roughness on the wetting transition and the formation of thin films of the defending liquid. We show experimentally, and rationalize theoretically, that roughness affects both the stability and dewetting dynamics of thin films, leading to distinct late-time morphologies of the undisplaced (trapped) fluid. Finally, we discuss the implications of our observations for geologic and technological applications.

42 ENGINEERING↗
DOE OSTI2023

Induced microseismicity and tremor signatures illuminate different slip behaviours in a natural shale fault reactivated by a fluid pressure stimulation (Mont Terri)

SUMMARY Fault slip induced by fluid perturbation in shale formations may only lead to a sparse seismicity. However, fault slip may strongly impact the integrity of shale formations that serve as caprocks for geological reservoirs holding buoyant fluids such as CO2, natural gas or hydrogen. A better understanding of the fluid reactivation processes of fault and the seismic triggering process is therefore critical for reservoir monitoring and fault stability. Here we analyse the seismic responses of a shale fault exposed to fluid pressurization during an in situ field-scale injection experiment at ∼340 m depth in the Mont Terri underground research laboratory (Switzerland). Two main types of seismic signals are observed as the fault was activated and started to slowly slip. After an aseismic phase, we observed tremor signatures and an increase in noise amplitude, which were directly associated with the slowly propagating fault slip in response to fluid injection. These signatures were later followed by micro-earthquakes that seem to occur further away from the fluid-pressurized area. We interpret these micro-earthquakes to be triggered by stress perturbations from the main slip growth. These two classes of seismic responses therefore highlight two different processes. Tremors seem to be a more direct observation for the fluid-induced slip propagation than micro-earthquakes. Even hidden in the noise, they precede earthquake failures, thus providing a useful tool for monitoring fluid leakage activated by slow deformation on low permeable shale faults, with applications for sealing integrity of caprocks.

De Barros, Louis (ORCID:0000000255419162)↗
DOE OSTI2022

Modelling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

ABSTRACT Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modelling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron–ion or electron–positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann–Vlasov equation. This new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes, not present in previous formulations of two-fluid magnetohydrodynamics. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. We keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The out-of-equilibrium corrections take the form of a collisional 14-moment closure previously described in the context of viscous single fluids. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Evaluation of CO 2 -Fluid-Rock Interaction in Enhanced Geothermal Systems: Field-Scale Geochemical Simulations

Recent studies suggest that using supercritical CO 2 (scCO 2 ) instead of water as a heat transmission fluid in Enhanced Geothermal Systems (EGS) may improve energy extraction. While CO 2 -fluid-rock interactions at “typical” temperatures and pressures of subsurface reservoirs are fairly well known, such understanding for the elevated conditions of EGS is relatively unresolved. Geochemical impacts of CO 2 as a working fluid (“CO 2 -EGS”) compared to those for water as a working fluid (H 2 O-EGS) are needed. The primary objectives of this study are (1) constraining geochemical processes associated with CO 2 -fluid-rock interactions under the high pressures and temperatures of a typical CO 2 -EGS site and (2) comparing geochemical impacts of CO 2 -EGS to geochemical impacts of H 2 O-EGS. The St. John’s Dome CO 2 -EGS research site in Arizona was adopted as a case study. A 3D model of the site was developed. Net heat extraction and mass flow production rates for CO 2 -EGS were larger compared to H 2 O-EGS, suggesting that using scCO 2 as a working fluid may enhance EGS heat extraction. More aqueous CO 2 accumulates within upper- and lower-lying layers than in the injection/production layers, reducing pH values and leading to increased dissolution and precipitation of minerals in those upper and lower layers. Dissolution of oligoclase for water as a working fluid shows smaller magnitude in rates and different distributions in profile than those for scCO 2 as a working fluid. It indicates that geochemical processes of scCO 2 -rock interaction have significant effects on mineral dissolution and precipitation in magnitudes and distributions.

58 GEOSCIENCES↗
DOE OSTI2020

High-temperature and high-pressure NMR investigations of low viscous fluids confined in mesoporous systems

Abstract In this contribution, the relaxation and diffusional behaviors of low viscous fluids, water and methanol confined into mesoporous silica and controlled size pore glass were investigated. The engineered porous systems are relevant to geologically important subsurface energy materials. The engineered porous proxies were characterized by Brunauer–Emmett–Teller (BET) surface analyzer, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (EM) to determine surface area, pore-wall protonation and morphology of these materials, respectively. The confined behavior of the low viscous fluids was studied by varying pore diameter, fluid-to-solid ratio, temperature, and pressure, and then compared to bulk liquid state. Both relaxation and diffusion behaviors for the confined fluids showed increasing deviation from pure bulk fluids as the fluid-to-solid ratio was decreased, and surface-to-volume ratio ( S / V ) was varied. Variable pressure deuteron NMR relaxation of confined D 2 O and confined methanol, deuterated at the hydroxyl or methyl positions, were performed to exploit the sensitivity of the deuteron quadrupole moment to molecular rotation. The methanol results demonstrated greater pressure dependence than those for water only in bulk. The deviations from bulk liquid behavior arise from different reasons such as confinement and the interactions between confined fluid and the nano-pore wall. The results of the present report give insight into the behavior of low viscosity fluid in nano-confined geometries under different state conditions.

Chemistry↗
DOE OSTI2020

Partitioning of Complex Fluids at Mineral Surfaces

This report summarizes the results obtained during the LDRD project entitled "Partitioning of Complex Fluids at Mineral Interfaces." This research addressed fundamental aspects of such interfaces, which are relevant to energy-water applications in the subsurface, including fossil energy extraction and carbon sequestration. This project directly addresses the problem of selectivity of complex fluid components at mineral-fluid interfaces, where complex fluids are defined as a mixture of hydrophobic and hydrophilic components: e.g., water, aqueous ions, polar/nonpolar organic compounds. Specifically, this project investigates how adsorption selectivity varies with surface properties and fluid composition. Both experimental and molecular modeling techniques were used to better understand trends in surface wettability on mineral surfaces. The experimental techniques spanned the macroscale (contact angle measurements) to the nanoscale (cryogenic electronic microscopy and vibrational spectroscopy). We focused on an anionic surfactant and a well-characterized mineral phase representative of clay phases present in oil- and gas-producing shale deposits. Collectively, the results consistently demonstrate that the presence of surfactant in the aqueous fluid significantly affects the mineral-fluid interfacial structure. Experimental and molecular modeling results reveal details of the surfactant structure at the interface, and how this structure varies with surfactant coverage and fluid composition.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Fluid Induced Elemental and Mineralogy Alterations of Caney Shale

ABSTRACT: Studies were conducted on cores and rock cuttings recovered from two wells drilled in Caney shale in the Southern Oklahoma. Mineralogical compositions of these samples were obtained by X-Ray Diffraction (XRD) measurements using two methods of sample preparation: conventional powder (by crushing and grinding of rock fragments) and spray-drying of micronized rock powder method. Microstructure and microanalysis were acquired by Scanning Electron Microscopy/Energy Dispersive Spectroscopy (SEM/EDS). Batch experiments were then conducted using selected rock samples and fracturing fluid whilst deionized water was used as control sample for fluid/rock interaction. These were conducted at 95°C and ambient pressure for various time periods to assess the geochemical reactivity of the Caney shale. The rock samples were subsequently recovered from fluid, dried and analysis repeated using identical XRD and SEM/EDS measurements, following the same conditions. Preliminary experimental results show rock mineralogical compositions are predominantly quartz with relatively high clay content. The clay minerals are mainly illite. Carbonate content is less than 10% by weight for most samples, with other non-clay minerals pyrite, albite, dolomite, and some apatite. In general, sprayed samples showed more intense peaks and less full width at half maxima (FWHM) relative to conventional rock-powder samples. Post-experimental mineralogical changes were observed in samples with clay portions breaking down to amorphous entities. As expected, dissolution of carbonates (calcite and dolomite) was observed and newly detected dissolution of biogenic micro quartz resulted in formation of amorphous silica. 1.0 INTRODUCTION Formation of sedimentary rocks involves processes where rock sediments are emplaced with fluids starting interactions that persist over long periods. During the lithification process, sediments and fluids react and attain equilibrium state. This equilibrium is disturbed during drilling, completion, hydraulic fracturing and other operations that involve the introduction of engineered fluids into the subsurface. To return to equilibrium, dissolution and precipitation of minerals may occur. These reactions have implications on the petrophysical properties of the reservoir. Whiles dissolution and precipitation of minerals can enhance porosity and permeability, the reverse can also be true. For example, when acidic fluids are introduced into formations with high carbonate contents, carbonate minerals dissolve into solution thus increasing the porosity and permeability of the formation. In contrast, when alkaline fluids with pH above 11 are introduced into an illite-rich shale formation, the rock-fluid interaction leads to de-flocculation of illite platelets causing substantial increase in clay fines migration, consequently significant reduction in permeability (Salles et al., 2008).

Awejori, G.↗