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.

SourceResearch connectionFrom yearThrough year

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

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Modeling of wave processes in blocky media with porous and fluid-saturated interlayers

The wave processes in blocky media are analyzed by applying different mathematical models, wherein the elastic blocks interact with each other via pliant interlayers with the complex mechanical properties. Four versions of constitutive equations are considered. In the first version, an elastic interaction between the blocks is simulated within the framework of linear elasticity theory, and the model of elastic–plastic interlayers is constructed to take into account the appearance of irreversible deformation of interlayers at short time intervals. In the second one, the effects of viscoelastic shear in the interblock interlayers are taken into the consideration using the Poynting–Thomson rheological scheme. In the third option, the model of an elastic porous material is used in the interlayers, where the pores collapse if an abrupt compressive stress is applied. In the fourth case, the model of a fluid-saturated material with open pores is examined based on Biot's equations. The collapse of pores is modeled by the generalized rheological approach, wherein the mechanical properties of a material are simulated using four rheological elements. Three of them are the traditional elastic, viscous and plastic elements, the fourth element is the so-called rigid contact, which is used to describe the behavior of materials with the different resistance to tension and compression. It was shown that the thermodynamically consistent model is provided, which means that the energy balance equation is fulfilled for an entire blocky structure, where the kinetic and potential energy of the system is the sum of the kinetic and potential energies of the blocks and interlayers. Under numerical implementation of the interlayers models, the dissipationless finite difference Ivanov's method was used. The splitting method by spatial variables in the combination with the Godunov gap decay scheme was applied in the blocks. As a result, robust and stable computational algorithms are built and tested. Using MPI technology, the parallel software was designed for the modeling of wave processes in 2D setting. The numerical results are presented, discussed and future studies are outlined. - Highlights: • We model the propagation of waves in a blocky medium with compliant interlayers. • Models of porous and fluid saturated interlayers are proposed using Biot's approach. • A parallel computational algorithm is developed for supercomputers of cluster type. • Problems of instant rotation of central and boundary blocks are solved numerically. • Wave motions are compared in the cases of elastic and fluid-saturated interlayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Characterization of the mud displacement in an enlarged wellbore: An integrated rock-fluid model

Cement-mud displacement plays a crucial role in the sealability of cement sheaths. Irregular geometric features of a wellbore due to washout can have a negative impact on mud and cement mobilization. An unstable interface between two fluids always leads to mud channeling, interfluid mixing, and cement contamination, degrading the cement quality. Many factors, such as mechanical and rheological properties of fluids, annulus geometry, flow pattern, and flow rate, significantly influence the displacement efficiency. This study investigates the characterization of the mud displacement in an irregular horizontal well using a 3D computational fluid dynamics (CFD) model. Mud is displaced in an enlarged wellbore by geopolymer and neat class G cement. The wellbore geometry is developed based on the caliper log data from an unconventional shale well in the Tuscaloosa Marine Shale (TMS) lithology. The effects of pump rate, density difference, and mud contamination are evaluated by numerical simulations. The results present those residual muds mainly exist in the upper annulus of the enlarged section. Geopolymer has a better sealing performance and can resist more water-based mud (WBM) contaminations than neat class G cement. The scenario with a low mud-cement density difference and high cement injection rate results in a high cement volume fraction, mitigating the gas migration.

58 GEOSCIENCES↗

DOE OSTI2021

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2017

Local and global Λ polarization in a vortical fluid

We compute the fermion spin distribution in the vortical fluid created in off-central high energy heavy-ion collisions. We employ the event-by-event (3+1)D viscous hydrodynamic model. The spin polarization density is proportional to the local fluid vorticity in quantum kinetic theory. As a result of strong collectivity, the spatial distribution of the local vorticity on the freeze-out hyper-surface strongly correlates to the rapidity and azimuthal angle distribution of fermion spins. We investigate the sensitivity of the local polarization to the initial fluid velocity in the hydrodynamic model and compute the global polarization of Λ hyperons by the AMPT model. The energy dependence of the global polarization agrees with the STAR data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Device-scale computational fluid dynamics modeling of carbon dioxide absorption using encapsulated sorbents

Micro-encapsulated carbon dioxide sorbent (MECS) is considered to be a promising material for carbon capture system because of high capture rate, low fabrication cost, and efficient energy consumption. To accelerate the development of MECS technology for industrial deployment, device-scale computational fluid dynamics (CFD) simulations are performed to study hydrodynamics and CO2 absorption behavior in a conceptual fluidized bed reactor. The two-fluid model solver, built in the Multiphase Flow with Interphase eXchanges (MFIX), is adopted to solve the multiphase flow hydrodynamics by representing both fluid and solids as continuous phases. Filtered models are implemented to capture the sub-grid effects that normally cannot be directly considered because MECS particles are too small to be explicitly resolved at such a large scale. Using established CFD tools, this study aims to characterize MECS behavior and predict its carbon capture performance in a conceptual device-scale absorber. The MECS particle size and gas flow rate have significant impact on CO2 capture efficiency. The CO2 capture efficiency decreases with increasing gas flow rate. Smaller particle size lead to better gas-solid homogeneity, as a results CO2 absorption efficiency increases.

Wang, Chao (ORCID:0000000279406614)↗

DOE OSTI2017

From 2D to 3D in fluid turbulence: unexpected critical transitions

How do the laws of physics change with changes in spatial dimension? Maybe not at all in some cases, but in important cases, the changes are dramatic. Fluid turbulence – the fluctuating, intermittent and many-degree-of-freedom state of a highly forced fluid – determines the transport of heat, mass and momentum and is ubiquitous in nature, where turbulence is found on spatial scales from microns to millions of kilometres (turbulence in stars) and beyond (galactic events such as supernovae). When the turbulent degrees of freedom are suppressed in one spatial dimension, the resulting turbulent state in two dimensions (2D) is remarkably changed compared with the turbulence in three dimensions (3D) – energy flows to small scales in 3D but towards large scales in 2D. Although this result has been known since the 1960s due to the pioneering work of Kraichnan, Batchelor and Leith, how one transitions between 3D and 2D turbulence has remained remarkably unexplored. For real physical systems, this is a highly significant question with important implications about transport in geophysical systems that determine weather on short time scales and climate on longer scales. Is the transition from 3D to 2D smooth or are there sharp transitions that signal a threshold of the dominance of one type of turbulence over another? Finally, recent results by Benavides & Alexakis (J. Fluid Mech., vol. 822 (2017), pp. 364–385) suggest that the latter may be the case – a surprising and provocative discovery.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗