Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID FLOW”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microfracture-pore structure characterization and water-rock interaction in three lithofacies of the Lower Eagle Ford Formation

The microfracture-pore structures and water-rock interactions of eight samples of three typical lithofacies (three limestone, three wackestone, and two mudstone) from the Lower Eagle Ford Formation have been investigated using an integrated methodology. The methods used were X-ray diffraction (XRD), total organic carbon (TOC) content, pyrolysis, thin-section petrography, scanning electron microscopy (SEM), mercury intrusion porosimetry (MIP), spontaneous water imbibition, and water vapor & nitrogen physisorption. Petrophysical properties such as organic richness, maturity, mineral composition, types and patterns of microfractures and pores, porosity, pore-throat size distribution, pore size distribution, pore connectivity, water adsorption and desorption behavior were determined from these methods, and the results compared among the three lithofacies. For example, the three limestone EAS samples and three wackestone DM samples show a low porosity range from 2.02–3.35%, whereas the two mudstone VRU samples exhibit a much large average porosity of 11.7% by MIP results. Nitrogen physisorption results show that pore size distributions have a good similarity in N2 adsorption volume for the samples with similar lithology. Microfractures are morphologically divided into seven types, with their definitions & pictorial examples presented and possible formation mechanisms proposed. In addition, the water vapor physisorption behavior and their controlling factors are discussed. Here, the water adsorption capacity shows a positive relationship (R 2 = 0.90) with clay content and a negative relationship (R 2 = 0.99) with carbonate content. In general, this integrated study of microfracture-pore structure characterization and water-rock interaction provides some insights of water vapor adsorption behavior in different lithofacies and whether microfractures could be an important pathway for fluid flow.

58 GEOSCIENCES↗
DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2021

Significance and complications of drilling fluid rheology in geothermal drilling: A review

The harsh downhole conditions of high pressure and high temperature (HPHT) encountered in geothermal wells make the drilling operation challenging. Drilling in such environments requires a special drilling mud formulation with high thermal stability and good rheological properties to fulfill the drilling fluid functions. Therefore, great efforts should be put into selecting the suitable drilling fluid, optimize and monitor the drilling fluid properties throughout drilling operations, and predicting its performance under downhole conditions. Rheological properties significantly impact many drilling parameters such as hole cleaning, fluid and wellbore stability, wellbore hydraulics, torque and drag, and other drilling issues. This paper discusses water-based drilling fluids' flow behavior under HPHT conditions and highlights the significance of fluid rheology in geothermal drilling. The common challenges and complications related to fluid rheology encountered in geothermal drilling are addressed in this paper, such as hole cleaning, wellbore hydraulics, and drilling fluid stability. This article also reviews the recent advances in drilling mud systems, rheology enhancement, and rheological properties measurements at surface and subsurface conditions. Furthermore, the rheology models of drilling fluid at elevated temperatures are reviewed to fully understand their flow behavior and establish a method for drilling engineers to optimize fluid formulations for geothermal drilling.

58 GEOSCIENCES↗
DOE OSTI2022

Optimization of flow in additively manufactured porous columns with graded permeability

Chemical engineering systems often involve a functional porous medium, such as in catalyzed reactive flows, fluid purifiers, and chromatographic separations. Ideally, the flow rates throughout the porous medium are uniform, and all portions of the medium contribute efficiently to its function. The permeability is a property of a porous medium that depends on pore geometry and relates flow rate to pressure drop. Additive manufacturing techniques raise the possibilities that permeability can be arbitrarily specified in three dimensions, and that a broader range of permeabilities can be achieved than by traditional manufacturing methods. Using numerical optimization methods, we show that designs with spatially varying permeability can achieve greater flow uniformity than designs with uniform permeability. We consider geometries involving hemispherical regions that distribute flow, as in many glass chromatography columns. By several measures, significant improvements in flow uniformity can be obtained by modifying permeability only near the inlet and outlet.

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

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗
DOE OSTI2024

Dataset for scientific paper "Simulated plant‑mediated oxygen input has strong impacts on fine‑scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands", a modeling study based on field observation at the tidal salt marshes of the Parker River Estuary, Massachusetts, United States

This dataset is the raw and processed data for the paper "Simulated plant ‑ mediated oxygen input has strong impacts on fine ‑ scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands". This study investigated how plant-mediated oxygen input affects subsurface biogeochemical reactions of organic carbon degradation and the resulting methane emissions of coastal wetlands by model simulation. We used the subsurface geochemical simulator PFLOTRAN for the modeling, which produced the simulated changes in porewater chemical substances and methane emissions over 10 days under different scenarios of plant-mediated oxygen input.Specifically, this dataset contains: 1) the input files for PFLOTRAN of all simulation runs conducted in this study. Those files are with an extension of ".in", containing information of the biogeochemical reaction network (stoichiometry, reaction rate, Monod constants, etc), fluid flow rate and oxygen concentration in the fluid which together simulated the plant-mediated oxygen input, the configuration of artificial reactions that simulated the methane fluxes, etc. The PFLOTRAN input files are text files, which can be opened by NotePad, but running these input files will require proper installation of PFLOTRAN (instruction: https://documentation.pflotran.org/user_guide/how_to/installation/installation.html). 2) the raw and processed model output from PFLOTRAN of all simulation runs, and 3) the python scripts used to process the raw model output, including random allocation of root cells, converting raw data into organized formats, calculating the methane fluxes based on the model output, data visualization, etc. The raw and processed model output from PFLOTRAN are in .spydata format, which can be viewed with Python. and 3) the python scripts for data processing and analysis are programming scripts, which can be opened with Python.This modeling work, in particular the model parameterization of root density and initial conditions of porewater concentrations of biogeochemical substances, was based on field measurements at the salt marsh of the Upper Parker River Estuary, Massachusetts, United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Strain shift measured from stress-controlled oscillatory shear: Evidence for a continuous yielding transition and new techniques to determine recovery rheology measures

Understanding the yielding of complex fluids is an important rheological challenge that affects our ability to engineer and process materials for a wide variety of applications. Common theoretical understandings of yield stress fluids follow the Oldroyd–Prager formalism in which the material behavior below the yield stress is treated as solidlike, and above the yield stress as liquidlike, with an instantaneous transition between the two states. This formalism was built on a quasi-static approach to the yield stress, while most applications, ranging from material processing to end user applications, involve a transient approach to yielding over a finite timescale. Using stress-controlled oscillatory shear experiments, we show that yield stress fluids flow below their yield stresses. This is quantified through measuring the strain shift, which is the value about which the strain oscillates during a stress-controlled test and is a function of only the unrecoverable strain. Measurements of the strain shift are, therefore, measurements of flow having taken place. These experimental results are compared to the Herschel–Bulkley form of the Saramito model, which utilizes the Oldroyd–Prager formalism, and the recently published Kamani–Donley–Rogers (KDR) model, in which one constitutive equation represents the entire range of material responses. Scaling relationships are derived, which allow us to show why yield stress fluids will flow across all stresses, above and below their yield stress. Finally, derivations are presented that show strain shift can be used to determine average metrics previously attainable only through recovery rheology, and these are experimentally verified.

Mechanics↗
DOE OSTI2021

Rotational and reflectional equivariant convolutional neural network for data-limited applications: Multiphase flow demonstration

This article deals with approximating steady-state particle-resolved fluid flow around a fixed particle of interest under the influence of randomly distributed stationary particles in a dispersed multiphase setup using convolutional neural network (CNN). The considered problem involves rotational symmetry about the mean velocity (streamwise) direction. Thus, this work enforces this symmetry using SE(3)-equivariant, special Euclidean group of dimension 3, CNN architecture, which is translation and three-dimensional rotation equivariant. This study mainly explores the generalization capabilities and benefits of a SE(3)-equivariant network. Accurate synthetic flow fields for Reynolds number and particle volume fraction combinations spanning over a range of [86.22, 172.96] and [0.11, 0.45], respectively, are produced with careful application of symmetry-aware data-driven approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Pre-Sinkhole Seismicity at the Napoleonville Salt Dome: Implications for Local Seismic Monitoring of Underground Caverns

The formation of a sinkhole at the Napoleonville salt dome located near Bayou Corne in Assumption Parish, Louisiana, was preceded by an increase in seismicity, as observed 10 km away by a continuously recording broadband seismic sensor. The seismicity rate in the days prior to the sinkhole formation contained hundreds of long-period events, some of which could be detected using a power detector (short-term to long-term energy ratio) tuned to the duration and frequency content of the seismicity. A characteristic of some of these events is a resonance near 0.4 Hz, which could be due to either an acoustic resonance within the brine-filled cavity or possibly due to unsteady fluid flow through a crack. In this work, a more exhaustive catalog was obtained using some of the events detected by the power detector as templates in a correlation detector. Furthermore, the full pre-sinkhole sequence of detections could offer lessons in remote monitoring of underground caverns.

58 GEOSCIENCES↗
DOE OSTI2020

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Topology optimization of additively manufactured fluidic components free of internal support structures

This report integrates projection-based approaches for implementing overhang constraints with fluid topology optimization to design additively manufactured fluidic components that do not require internal support structures. Internal support structures are challenging, potentially impossible to remove in complex fluid networks, yet would degrade fluidic performance if left within the part. This presents a major challenge in coupling topology optimization with additive manufacturing in the design of fluid components. The proposed approach leverages past work in overhang constraints for solid mechanics, including projection formulations and adjoint sensitivity analysis, and considers incompressible Navier–Stokes equations for laminar fluid flow. The approach is demonstrated on 2D and 3D fluid topology optimization problems to minimize pressure drop and constrain mass flow rates in pipes and manifolds. Resulting designs are crisp, logical, satisfy performance constraints, and satisfy overhang constraints, eliminating the need for internal support structures.

42 ENGINEERING↗
DOE OSTI2025

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗
DOE OSTI2023

Integrated modeling of cracking during deep penetration laser welding of nickel alloys

During deep penetration laser welding of nickel alloys, interactions between composition and processing conditions can lead to the formation of defects. Inconel 740H, for example, has demonstrated a susceptibility to horizontal fusion zone cracking at locations between 70% and 80% of the weld depth during laser welding at powers above 5 kW. Coupling three-dimensional heat transfer, fluid flow, and stress modeling tools allowed that both the strain rate and stress normal to the solidification direction to be calculated. In the Inconel 740H welds, cracks formed at locations where the strain rate and stress simultaneously reached critical levels. No cracking was observed in the Inconel 690 welds, since the strain rate and stress did not simultaneously reach these critical levels.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Review of heat transfer enhancement techniques for single phase flows

The thermal energy exchange between a flowing fluid and its confining channel is a ubiquitous process in modern society. To enhance the fluid-to-wall or wall-to-fluid heat transfer, several techniques have been developed to maximize the contact area between the fluid and the inner wall and/or disrupt the flow to enhance circulation or induce turbulence. Deployment of channels having features capable of enhancing heat transfer enables the reduction of heat exchanger size while maintaining performance. Reduction in equipment size is critical due to the ability to minimize the required volume of costly working fluids and to mitigate potential safety concerns associated with total system fluid volume. Here, a comprehensive review of single-phase heat transfer enhancement techniques is presented. The article provides a thorough comparison by analyzing the heat transfer rate, pressure drop, and other operational aspects. Single-phase heat transfer enhancement methods are divided into active and passive techniques. Active methods such as electrohydrodynamic (EHD), magnetohydrodynamics (MHD), or mechanical motion require external power to create enhancement. Passive methods such as dimples, fins, or tape inserts do not require external input and rely only on surface modification. Although active methods are more expensive and difficult to implement compared to passive techniques, it enables active control of heat transfer augmentation. Finally, this review develops and summarizes key learning data for design optimization enabled by additive manufacturing and machine learning algorithms, helping to inform these next-generation heat exchanger design methodologies for a plethora of modern applications such as electrification of vehicles, computing, and classical industries.

42 ENGINEERING↗
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2023

Combining synchrotron X-ray diffraction, mechanistic modeling and machine learning for in situ subsurface temperature quantification during laser melting

Laser melting, such as that encountered during additive manufacturing, produces extreme gradients of temperature in both space and time, which in turn influence microstructural development in the material. Qualification and model validation of the process itself and the resulting material necessitate the ability to characterize these temperature fields. However, well established means to directly probe the material temperature below the surface of an alloy while it is being processed are limited. To address this gap in characterization capabilities, a novel means is presented to extract subsurface temperature-distribution metrics, with uncertainty, from in situ synchrotron X-ray diffraction measurements to provide quantitative temperature evolution data during laser melting. Temperature-distribution metrics are determined using Gaussian process regression supervised machine-learning surrogate models trained with a combination of mechanistic modeling (heat transfer and fluid flow) and X-ray diffraction simulation. The trained surrogate model uncertainties are found to range from 5 to 15% depending on the metric and current temperature. The surrogate models are then applied to experimental data to extract temperature metrics from an Inconel 625 nickel superalloy wall specimen during laser melting. The maximum temperatures of the solid phase in the diffraction volume through melting and cooling are found to reach the solidus temperature as expected, with the mean and minimum temperatures found to be several hundred degrees less. The extracted temperature metrics near melting are determined to be more accurate because of the lower relative levels of mechanical elastic strains. However, uncertainties for temperature metrics during cooling are increased due to the effects of thermomechanical stress.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tempering kinetics during multilayer laser additive manufacturing of a ferritic steel

Grade 91 steel forms martensite during additive manufacturing and the extent of tempering of martensite significantly affects the mechanical properties of parts. Currently, there is a lack of quantitative understanding of the tempering kinetics for Grade 91 steel, and as a result, the effects of repeated thermal cycles on properties for different processing conditions cannot be determined. Here we evaluate the tempering kinetics by determining the constant terms in the Johnson Mehl Avrami kinetic equation from the tempering data available in the literature and the thermal cycles computed using a rigorously-tested heat and fluid flow model of multi-layer additive manufacturing. The raw tempering data are cleaned using a neural network to enhance accuracy. The lower layers experience repeating cycles of heating and cooling when the upper layers are added. As a result, the hardness is reduced owing to the tempering of martensite. In contrast, martensite formed in the upper layers is not tempered to the same extent and the hardness remains high. Therefore, the hardness of the part increases with the distance from the substrate. Variations in the heat input at different laser powers and scanning speeds significantly affect the extent of tempering. Finally, since the method used here can provide a quantitative understanding of the tempering of martensite and the spatial variation in hardness, it can be used to tailor the microstructure and hardness of heat treatable printed metallic parts.

36 MATERIALS SCIENCE↗