Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID MECHANICS”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evanescent and inertial-like waves in rigidly rotating odd viscous liquids

Three-dimensional non-rotating odd viscous liquids give rise to Taylor columns and support axisymmetric inertial-like waves ( J. Fluid Mech ., vol. 973, 2023, A30). When an odd viscous liquid is subjected to rigid-body rotation however, there arise in addition a plethora of other phenomena that need to be clarified. In this paper, we show that three-dimensional incompressible or two-dimensional compressible odd viscous liquids, rotating rigidly with angular velocity 𝛺, give rise to both oscillatory and evanescent inertial-like waves or a combination thereof (which we call of mixed type) that can be non-axisymmetric. By evanescent, we mean that along the radial direction, typically when moving away from a solid boundary, the velocity field decreases exponentially. These waves precess in a prograde or retrograde manner with respect to the rotating frame. The oscillatory and evanescent waves resemble respectively the body and wall-modes observed in (non-odd) rotating Rayleigh–Bénard convection ( J. Fluid Mech ., vol. 248, 1993, pp. 583–604). We show that the three types of waves (wall, body or mixed) can be classified with respect to pairs of planar wavenumbers 𝜅 which are complex, real or a combination, respectively. Experimentally, by observing the precession rate of the patterns, it would be possible to determine the largely unknown values of the odd viscosity coefficients. This formulation recovers as special cases recent studies of equatorial or topological waves in two-dimensional odd viscous liquids which provided examples of the bulk–interface correspondence at frequencies 𝜔 < 2⁢𝛺. We finally point out that the two- and three-dimensional problems are formally equivalent. Their difference then lies in the way data propagate along characteristic rays in three dimensions, which we demonstrate by classifying the resulting Poincaré–Cartan equations.

Mechanics↗
DOE OSTI2019

Impact of the Rastall parameter on perfect fluid spheres

We examine the effects of the Rastall parameter on the behaviour of spherically symmetric static distributions of perfect fluid matter. It was claimed by Visser [Physics Letters B, 782, 83, (2018)] that the Rastall proposition is completely equivalent to the Einstein theory. While many authors have raised contrary arguments, our intention is to analyse the properties of Rastall gravity through variation of the Rastall parameter in the context of perfect fluids spheres that may be used to model neutron stars or cold fluid planets. This analysis also serves to counter the claim that Rastall gravity is equivalent to the standard Einstein theory. It turns out that the condition of pressure isotropy is exactly the same as for Einstein gravity and hence that any known solution of the Einstein equations may be used to study the effects of the Rastall dynamical quantities. Moreover, by choosing the well studied Tolman metrics, we discover that in the majority of cases there is substantial deviation from the Einstein case when the Rastall parameter vanishes and in cases where the Einstein model displays defective behaviour, certain Rastall models obey the well known elementary requirements for physical plausibility. These empirical findings do not support the idea that Rastall theory is equivalent to Einstein theory as several deviations in physical behaviour are displayed as counter-examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

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

Tuning the quantumness of simple Bose systems: A universal phase diagram

We present a comprehensive theoretical study of the phase diagram of a system of many Bose particles interacting with a two-body central potential of the so-called Lennard-Jones form. First-principles path-integral computations are carried out, providing essentially exact numerical results on the thermodynamic properties. The theoretical model used here provides a realistic and remarkably general framework for describing simple Bose systems ranging from crystals to normal fluids to superfluids and gases. The interplay between particle interactions on the one hand and quantum indistinguishability and delocalization on the other hand is characterized by a single quantumness parameter, which can be tuned to engineer and explore different regimes. Taking advantage of the rare combination of the versatility of the many-body Hamiltonian and the possibility for exact computations, we systematically investigate the phases of the systems as a function of pressure (P) and temperature (T), as well as the quantumness parameter. Here, we show how the topology of the phase diagram evolves from the known case of 4 He, as the system is made more (and less) quantum, and compare our predictions with available results from mean-field theory. Possible realization and observation of the phases and physical regimes predicted here are discussed in various experimental systems, including hypothetical muonic matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Surfactant inhibition mechanisms of carbonate mineral dissolution in shale

Surfactants are common additives to hydraulic fracturing and enhanced oil recovery (EOR) fluids, and are under consideration for amendment to supercritical carbon dioxide for geological carbon sequestration (GCS). The effect of a common anionic surfactant, internal olefin sulfonate (IOS), on mineral dissolution from shale into brine was evaluated. When added to brine at concentrations exceeding the critical micelle concentration (94 mg/L), IOS inhibited carbonate mineral dissolution in an Eagle Ford shale, as well as dissolution of optical quality calcite (the dominant carbonate in the shale). Laser profilometry images provide spatial resolution across > 3 orders of magnitude, and indicate that IOS addition to brine both enhances the formation of new etch pits in calcite, and impedes their further growth. Time-of-flight secondary ion mass spectrometry surface profiles show for the first time that IOS preferentially adsorbs at calcite pit edges versus flat calcite surfaces (i.e., terraces). Surface pressure calculations, sulfur K-edge near edge X-ray absorption fine structure (NEXAFS) spectroscopy results, and density functional theory (DFT) calculations support this observation; the DFT results indicate that the sulfonate head group of the IOS molecule binds strongly to the calcite step site as compared to the terrace site. The S K-edge NEXAFS results indicate that IOS adsorbed more to etched calcite surfaces compared to smooth calcite surfaces. Overall, the results indicate that weak adsorption on flat calcite surfaces (i.e., terraces) disrupts water structure and enhances mass transfer of dissolution, while strong adsorption on calcite pit edges displaces adsorbed water and inhibits further etch pit growth. This work provides the first direct evidence of preferential adsorption of IOS to etched calcite surfaces and links it to macroscopic dissolution kinetics. Finally, this work has implications for surfactant-containing fluids used in hydraulic fracturing, EOR and potentially GCS for subsurface injection into carbonate rich reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Analytical Solutions for Seismic Fluid-Structure Interaction of Head-Supported Cylindrical Tanks

Fluid-structure interaction (FSI) of earthquake shaking is a design consideration in the civil, nuclear, chemical, and mechanical industries. Preliminary sizing of structures and components in these industries will often use analytical solutions because there is insufficient information to warrant the analysis of numerical models. Here, this study extends prior analytical work on seismic FSI of base-supported tanks in the 1980s to accommodate head-supported tanks with geometries similar to those proposed for new-build, liquid (fluid)-filled, advanced reactors. Analytical solutions are derived for flexible, head-supported, cylindrical tanks subjected to small-amplitude, unidirectional, horizontal seismic motion, and address frequencies of lateral motions of the tank, hydrodynamic pressures, frequencies of waves, wave heights, and reactions at the head support. Solutions are presented for a range of tank dimensions and could be used for a) preliminary design of head-supported cylindrical tanks and their supporting structures with idealized seismic input, and b) verification of numerical models to be used for final seismic analysis, design, and probabilistic risk assessment.

42 ENGINEERING↗
DOE OSTI2022

Design optimization of gap distance for the capillary limitation of a heat pipe with annular-type wick structure

In this study, an experimental investigation was conducted on the rising height and contact angle of fluid in an annular wick-type heat pipe. The annular wick-type heat pipe was characterized by a small gap between the wick structure and tube wall, which compensated for the pressure drop along the porous media and created additional capillary force. To describe and model the advantage of this gap, the rising of a wetting liquid in the gap between a vertical solid plate and a mesh (with a small angle between them) was experimentally measured and analyzed. An additional experiment was performed to investigate the effect of curvature on the capillary rise using tubes and meshes of varying radii. Resultantly, we confirmed that the linear combination of the contact angles of the solid plate and mesh could be applied to calculate the rising height from the Laplace–Young equation. Furthermore, the effect of curvature on the rising height of the liquid was negligible. These results were extended to the investigation of finding the optimal gap distance for the annular wick-type heat pipe by referring to previous studies. We observed that a gap distance of 1.27 mm provided the largest permeability (K) over the effective pore radius (r eff ) value for a heat pipe with ethanol, which in turn resulted in the highest capillary limitation. For a sodium heat pipe, a gap distance of 0.84 mm resulted in the highest capillary limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗
DOE OSTI2018

Experimental Investigation of Time-Dependent Deformation of Fluid Disposal Reservoir Caprocks

For approximately the last two decades low-permeability formations have become a prominent research area because these are the target of unconventional oil shale exploration, and also the primary trapping mechanism for CO2 sequestration and other such subsurface storage. The objective of this research is focused on the sealing ability and induced seismic risk of such tight formations. Disposal of fluid into the subsurface alters the stress state in both the target reservoir and the surrounding formations from their initial in situ state. Therefore, caprocks must be able to deform and accommodate the altered stress and strain state while maintaining the integrity of the seal. The deformation induced in caprock formations often results in a decrease in the normal stress on faults and fractures, which may induce shear slip and therefore create the potential for induced seismicity and leakage. Deformations induced by fluid injection are generally not instantaneous, but typically occur gradually as a pressure front advances through the reservoir and caprock formations. As such, an injection program design can influence both the rate and magnitude of deformation and associated changes in stress. Many common caprock formation types, such as shales and evaporates, are prone to undergo time-dependent deformation such as viscoelastic and/or viscoplastic response. Both the magnitude of the time-dependent deformation and the associated characteristic time scale are important considerations for selecting or operating an injection site. This paper describes tests intended to characterize these effects using caprock samples from a variety of current or potential storage sites. Triaxial and load relaxation test data are presented and used to quantify the time-dependent deformation behavior of rocks with different textures and compositions.

Edelman, Eric↗
DOE OSTI2018

Scale matters

The applicability of Navier–Stokes equations is limited to near-equilibrium flows in which the gradients of density, velocity and energy are small. Here I propose an extension of the Chapman–Enskog approximation in which the velocity probability distribution function (PDF) is averaged in the coordinate phase space as well as the velocity phase space. I derive a PDF that depends on the gradients and represents a first-order generalization of local thermodynamic equilibrium. I then integrate this PDF to derive a hydrodynamic model. Finally, I discuss the properties of that model and its relation to the discrete equations of computational fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Chemical Trapping of CO2 by Clay Minerals at Reservoir Conditions: Two Mechanisms Observed by In Situ High Pressure and Temperature Experiments

The reaction of clay minerals in sedimentary rocks with supercritical CO2-rich fluids is poorly understood, but has potential implications for enhanced oil and gas production and geological C-sequestration. Experiments performed in situ at temperature and pressure relevant to reservoir conditions (T = 323 K and Pfluid = 90 bar) show that trioctahedral clay minerals can react with supercritical CO2 to produce carbonate phases by both ion exchange/precipitation reactions and dissolution/reprecipitation on a timescale of hours under certain conditions. The dissolution/reprecipitation reactions were observed in a synthetic, high surface area/high edge site surface area smectite (laponite) exchanged with Ca2+, Cs+, and tetramethyl ammonium (TMA+), and the ion exchange/precipitation mechanism was observed for a Pb-exchanged natural smectite (hectorite). Novel X-ray diffraction and NMR and infrared spectroscopic tools provide in situ observation of the reactions in real time supported by a suite of ex situ analyses. For the laponites, IR data show that HCO3- ion forms at water contents as small as ~5 H2O molecules/exchangeable cation. When the exchangeable cation is Ca2+, the IR data show the formation of carbonate as well, and the NMR results show formation of amorphous calcium carbonate at low water contents in addition to HCO3-. Laponites equilibrated at 100% R.H. at atmospheric conditions and then exposed to scCO2 generate additional, more mobile HCO3- ion and exhibit evidence of clay dissolution leading to a poorly crystalline or amorphous hydrous magnesium carbonate/bicarbonate that forms from Mg2+ released by dissolution of the octahedral sheet. The 100% R.H. sample with exchangeable Ca2+ also forms calcite, vaterite and aragonite precipitates. Comparison of these results with those already published in the literature suggest that a high edge site surface area is crucial to this process occurring on a short timescale. In the Pb-exchanged hectorite exposed to scCO2, once a critical humidity threshold is reached, cerussite (PbCO3) rapidly forms concurrent with exchange of interlayer Pb2+ with H3O+ formed by reaction of CO2 with water on the clay surface. Given that such a reaction is not observed on a similar timescale with more common Ca2+ or Na+ in hectorite and other smectites, the low solubility of cerussite appears to be the thermodynamic driving force for this process, meaning that such reactions are expected to occur more slowly in most natural samples. However, the results suggest that both types of CO2-capturing processes could occur on short and/or geological timescales during CO2 flooding of oil and gas reservoirs, in C-sequestration reservoirs, and as the result of fracking using CO2-based fluids.

Bowers, Geoffrey M.↗
DOE OSTI2023

Spin polarization induced by magnetic field and the relativistic Barnett effect

First, I study the analogy between the magnetization of a material and the spin polarization of particles in a fluid. Using the relativistic version of the Barnett effect, i.e. the magnetization of a material induced by mechanical rotation, the spin polarization induced by thermal vorticity is obtained within a purely classical model, where spin is treated as an intrinsic magnetic moment and rotation is included as a non-inertial effect. I argue that since spin polarization induced by thermal vorticity can be obtained in a classical theory, it can not be dominated by quantum anomalies. Second, the spin polarization induced by magnetic field is obtained for a fluid at local thermal equilibrium using statistical quantum field theory. The obtained formula is valid beyond the weak field approximation and when contributions from the non-homogeneity of the magnetic field are small. The exact form of spin polarization is studied for a free Dirac field at global equilibrium, and, like magnetic susceptibility, it oscillates according to the de Haas - van Alphen effect. Finally, I briefly review how magnetic field contributes to the difference between the spin polarization of Λ and $\barΛ$ observed in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Effect of Electronic Inertia on the Gravito-Electrostatic Sheath Structure Formation

The gravito-electrostatic sheath (GES) model, previously proposed to address the fundamental issues on the surface emission mechanism of outflowing solar plasma on the basis of plasma−wall interaction processes with inertialess electrons on both bounded and unbounded scales, is reformulated in the light of active electron inertial response amid geometrical curvature effects. We accordingly derive the electron population distribution law considering both weak electron inertia and geometrical curvature effects in a new analytic construct coupled with the GES structure equations in a closed form. The analysis shows that the GES characteristics and hence plasma outflow dynamics are noticeably affected because of electron inertia. As a consequence of the electron inertia inclusion in contrast with the previous GES formalism, it is found that the GES width gets reduced (–5%), the sheath boundary gets contracted (–7%), the net current density at the surface gets reduced (–25%), the GES potential enhances (+17%), the transonic horizon decreases (‒35%), self-gravity enhances (+2%), and so forth. The obtained results are in fair accord with the existing model predictions centered around both the earlier GES formalisms and standard fluid-kinetic predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Study on the Microstructure, Mechanical Properties and Corrosion Behavior of Mg-Zn-Ca Alloy Wire for Biomaterial Application

Due to their excellent biocompatibility and biodegradability, magnesium alloy wires have attracted much attention for biomaterial applications including orthopedic K-wires and sutures in wound closure. In this study, Mg-3Zn-0.2Ca alloy wires were prepared by cold drawing combined with proper intermediate annealing process. Microstructures, texture, mechanical properties and corrosion behavior of Mg-3Zn-0.2Ca alloy wire in a simulated body fluid were investigated. The results showed that the secondary phase and average grain size of the Mg-3Zn-0.2Ca alloy were refined in comparison with the as-extruded alloy and a strong (0002)<10-10>//DD basal fiber texture system was formed after multi-pass cold drawing. After the annealing, most of the basal planes were tilted to the drawing direction (DD) by about 35°, presenting the characteristics of random texture, and the texture intensity decreased. The as-annealed wire shows good mechanical properties with the ultimate tensile strength (UTS), yield strength (YS) and elongation of 253 ± 8.5 MPa, 212 ± 11.3 MPa and 9.2 ± 0.9%, respectively. Electrochemical and hydrogen evolution measurements showed that the corrosion resistance of the Mg-3Zn-0.2Ca alloy wire was improved after the annealing. The immersion test indicated that the Mg-3Zn-0.2Ca wire exhibited uniform corrosion behavior during the initial period of immersion, but then exhibited local corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reactivity of CO 2 with Utica, Marcellus, Barnett, and Eagle Ford Shales and Impact on Permeability

We report that in order to reduce greenhouse gas emissions while recovering hydrocarbons from unconventional shale formations, processes that make use of carbon dioxide to enhance oil recovery while storing carbon dioxide (CO 2 ) should be considered. Here, we examine samples from three shale basins across the United States (Utica and Marcellus Shales in the Appalachian Basin, Barnett Shale in the Bend Arch-Ft. Worth Basin, and Eagle Ford in the Western Gulf Basin) to address the following questions: (1) do changes from reaction with CO 2 and fluids at the micrometer and nanometer scale alter flow pathways and, in turn, impact hydrocarbon production, CO 2 storage, and seal integrity and (2) can CO 2 or fluid reactivity be predicted based on physical or chemical properties of shale formations? Experiments were conducted at 40 °C and 10.3 MPa to characterize the interaction between CO 2 and shale using X-ray diffraction (XRD), carbon and sulfur analysis, in situ Fourier transform infrared spectroscopy (FT-IR), feature relocation scanning electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), mercury (Hg) intrusion porosimetry, and Brunauer-Emmett-Teller (BET) surface area and pore size analysis coupled with density functional theory (DFT) methods. Changes in mechanical, physical, and flow properties of shale cores due to CO 2 exposure were addressed using a New England Research Autolab 1500 and Xenon X-ray computed tomography (CT) scanning. Results showed that CO 2 did not promote significant reactivity with the shale if water was not present; only shales with swelling clays or residual interstitial pore water reacted with dry CO 2 to promote reactivity in shale. When water was added as a reactant, CO 2 formed carbonic acid and reacted with the shale to dissolve carbonate pockets, etched and pitted the shale matrix surfaces, and increased the microporosity and decreased nanoporosity. Porosity and permeability increased appreciably in core shale samples after exposure to CO 2 saturated fluid due to dissolution of carbonate. Shale mechanical properties were not altered. Trends were not observed that could tie CO 2 or fluid reactivity to physical or chemical properties of the shale formations at the basin scale from the samples we examined. However, if the shale contained significant amounts of carbonate and water was available to react with the CO 2 , pore sizes were altered in the matrix and permeability and porosity increased.

04 OIL SHALES AND TAR SANDS↗