Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Fracture Propagation, Fluid Flow, and Geomechanics of Water-Based Hydraulic Fracturing in Shale Gas Systems and Electromagnetic Geophysical Monitoring of Fluid Migration

We investigate fracture propagation induced by hydraulic fracturing with water injection, using numerical simulation. For rigorous, full 3D modeling, we employ a numerical method that can model failure resulting from tensile and shear stresses, dynamic nonlinear permeability, leak-off in all directions, and thermo-poro-mechanical effects with the double porosity approach. Our numerical results indicate that fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics, whereas water saturation is determined by fluid flow. At early times, the water saturation front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with injected water. However, at late times, advance of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, which is filled with reservoir gas. We also find considerable leak-off of water to the reservoir. The inconsistency between the fracture volume and the volume of injected water cannot properly calculate the fracture length, when it is estimated based on the simple assumption that the fracture is fully saturated with injected water. As an example of flow-geomechanical responses, we identify pressure fluctuation under constant water injection, because hydraulic fracturing is itself a set of many failure processes, in which pressure consistently drops when failure occurs, but fluctuation decreases as the fracture length grows. We also study application of electromagnetic (EM) geophysical methods, because these methods are highly sensitive to changes in porosity and pore-fluid properties due to water injection into gas reservoirs. Employing a 3D finite-element EM geophysical simulator, we evaluate the sensitivity of the crosswell EM method for monitoring fluid movements in shaly reservoirs. For this sensitivity evaluation, reservoir models are generated through the coupled flow-geomechanical simulator and are transformed via a rock-physics model into electrical conductivity models. It is shown that anomalous conductivity distribution in the resulting models is closely related to injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for distribution of injected fluids (i.e., migration pathways) during hydraulic fracturing operations

02 PETROLEUM↗
DOE OSTI2022

Sensitivity Analysis of Fluid–Fluid Interfacial Area, Phase Saturation and Phase Connectivity on Relative Permeability Estimation Using Machine Learning Algorithms

Recent studies have shown that relative permeability can be modeled as a state function which is independent of flow direction and dependent upon phase saturation (S), phase connectivity (X), and fluid–fluid interfacial area (A). This study evaluates the impact of each of the three state parameters (S, X, and A) in the estimation of relative permeability. The relative importance of the three state parameters in four separate quadrants of S-X-A space was evaluated using a machine learning algorithm (out-of-bag predictor importance method). The results show that relative permeability is sensitive to all the three parameters, S, X, and A, with varying magnitudes in each of the four quadrants at a constant value of wettability. We observe that the wetting-phase relative permeability is most sensitive to saturation, while the non-wetting phase is most sensitive to phase connectivity. Although the least important, fluid–fluid interfacial area is still important to make the relative permeability a more exact state function.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Coupling fast fluid dynamics and multizone airflow models in Modelica Buildings library to simulate the dynamics of HVAC systems

Multizone models are widely used in building airflow and energy performance simulations due to their fast computing speed. However, multizone models assume that the air in a room is well mixed, consequently limiting their application. In specific rooms where this assumption fails, the use of computational fluid dynamics (CFD) models may be an alternative option. Previous research has mainly focused on coupling CFD models and multizone models to study airflow in large spaces. While significant, most of these analyses did not consider the coupled simulation of the building airflow with the building’s Heating, Ventilation, and Air-Conditioning (HVAC) systems. This paper tries to fill the gap by integrating the models for HVAC systems with coupled multizone and CFD simulations for airflows, using the Modelica simulation platform. To improve the computational efficiency, we incorporated a simplified CFD model named fast fluid dynamics (FFD). We first introduce the data synchronization strategy and implementation in Modelica. Then, we verify the implementation using two case studies involving an isothermal and a non-isothermal flow by comparing model simulations to experiment data. Afterward, we study another three cases that are deemed more realistic. Furthermore, this is done by attaching a variable air volume (VAV) terminal box and a VAV system to previous flows to assess the capability of the models in studying the dynamic control of HVAC systems. Finally, we discuss further research needs on the coupled simulation using the models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Fluid Flow Investigations within a 37 Element CANDU Fuel Bundle Supported by Magnetic Resonance Velocimetry and Computational Fluid Dynamics

The current work presents experimental and computational investigations of fluid flow through a 37 element CANDU nuclear fuel bundle. Experiments based on Magnetic Resonance Velocimetry (MRV) permit three-dimensional, three-component fluid velocity measurements to be made within the bundle with sub-millimeter resolution that are non-intrusive, do not require tracer particles or optical access of the flow field. Computational fluid dynamic (CFD) simulations of the foregoing experiments were performed with the hydra-th code using implicit large eddy simulation, which were in good agreement with experimental measurements of the fluid velocity. Greater understanding has been gained in the evolution of geometry-induced inter-subchannel mixing, the local effects of obstructed debris on the local flow field, and various turbulent effects, such as recirculation, swirl and separation. These capabilities are not available with conventional experimental techniques or thermal-hydraulic codes. Finally, the overall goal of this work is to continue developing experimental and computational capabilities for further investigations that reliably support nuclear reactor performance and safety.

42 ENGINEERING↗
DOE OSTI2021

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

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

A fast non-Fourier method for Landau-fluid operators

An efficient and versatile non-Fourier method for the computation of Landau-fluid (LF) closure operators [Hammett and Perkins, Phys. Rev. Lett. 64, 3019 (1990)] is presented, based on an approximation by a sum of modified-Helmholtz-equation solves (SMHS) in configuration space. This method can yield fast-Fourier-like scaling of the computational time requirements and also provides a very compact data representation of these operators, even for plasmas with large spatial nonuniformity. As a result, the method can give significant savings compared with direct application of “delocalization kernels” [e.g., Schurtz et al., Phys. Plasmas 7, 4238 (2000)], both in terms of computational cost and memory requirements. The method is of interest for the implementation of Landau-fluid models in situations where the spatial nonuniformity, particular geometry, or boundary conditions render a Fourier implementation difficult or impossible. Systematic procedures have been developed to optimize the resulting operators for accuracy and computational cost. The four-moment Landau-fluid model of Hammett and Perkins has been implemented in the BOUT++ code using the SMHS method for LF closure. Excellent agreement has been obtained for the one-dimensional plasma density response function between driven initial-value calculations using this BOUT++ implementation and matrix eigenvalue calculations using both Fourier and SMHS non-Fourier implementations of the LF closures. The SMHS method also forms the basis for the implementation, which has been carried out in the BOUT++ code, of the parallel and toroidal drift-resonance LF closures. The method is a key enabling tool for the extension of gyro-Landau-fluid models [e.g., Beer and Hammett, Phys. Plasmas 3, 4046 (1996)] to codes that treat regions with strong profile variation, such as the tokamak edge and scrapeoff-layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A comparative study of coarse-graining methods for polymeric fluids: Mori-Zwanzig vs. iterative Boltzmann inversion vs. stochastic parametric optimization

We construct effective coarse-grained (CG) models for polymeric fluids by employing two coarse-graining strategies. The first one is a forward-coarse-graining procedure by the Mori-Zwanzig (MZ) projection while the other one applies a reverse-coarse-graining procedure, such as the iterative Boltzmann inversion (IBI) and the stochastic parametric optimization (SPO). More specifically, we perform molecular dynamics (MD) simulations of star polymer melts to provide the atomistic fields to be coarse-grained. Each molecule of a star polymer with internal degrees of freedom is coarsened into a single CG particle and the effective interactions between CG particles can be either evaluated directly from microscopic dynamics based on the MZ formalism, or obtained by the reverse methods, i.e., IBI and SPO. The forward procedure has no free parameters to tune and recovers the MD system faithfully. For the reverse procedure, we find that the parameters in CG models cannot be selected arbitrarily. If the free parameters are properly defined, the reverse CG procedure also yields an accurate effective potential. Moreover, we explain how an aggressive coarse-graining procedure introduces the many-body effect, which makes the pairwise potential invalid for the same system at densities away from the training point. From this work, general guidelines for coarse-graining of polymeric fluids can be drawn.

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

Quasistatic fluid-fluid displacement in porous media: Invasion-percolation through a wetting transition

We study the influence of wettability on the morphology of fluid-fluid displacement through analog porous media in the limit of vanishing flow rates. We introduce an invasion-percolation model that considers cooperative pore filling and corner-flow mechanisms and captures interface motion at the pore scale for all quasistatic flow regimes between strong drainage and strong imbibition. We validate the method against recent experimental observations of wetting transition in microfluidic cells patterned with circular posts and we use it to explore the sensitivity of fluid invasion to wettability heterogeneity, post spacing, and post height. Furthermore, our model therefore extends the Cieplak-Robbins description of quasistatic fluid invasion by reproducing the wetting transition in strong imbibition, a feature that requires incorporating three-dimensional effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Discrete models of fluids: spatial averaging, closure and model reduction

We consider semidiscrete ODE models of single-phase fluids and two-fluid mixtures. In the presence of multiple fine-scale heterogeneities, the size of these ODE systems can be very large. Spatial averaging is then a useful tool for reducing computational complexity of the problem. The averages satisfy exact balance equations of mass, momentum, and energy. These equations do not form a satisfactory continuum model because evaluation of stress and heat flux requires solving the underlying ODEs. To produce continuum equations that can be simulated without resolving microscale dynamics, we recently proposed a closure method based on the use of regularized deconvolution. Here we continue the investigation of deconvolution closure with the long term objective of developing consistent computational upscaling for multiphase particle methods. The structure of the fine-scale particle solvers is reminiscent of molecular dynamics. For this reason we use nonlinear averaging introduced for atomistic systems by Noll, Hardy, and Murdoch-Bedeaux. We also consider a simpler linear averaging originally developed in large eddy simulation of turbulence. We present several simple but representative examples of spatially averaged ODEs, where the closure error can be analyzed. Based on this analysis we suggest a general strategy for reducing the relative error of approximate closure. For problems with periodic highly oscillatory material parameters we propose a spectral boosting technique that augments the standard deconvolution and helps to correctly account for dispersion effects. We also conduct several numerical experiments, one of which is a complete mesoscale simulation of a stratified two-fluid flow in a channel. In this simulation, the operation count per coarse time step scales sublinearly with the number of particles.

Multiscale modeling, discrete models, fluid flow↗
DOE OSTI2025

Characterizing the pore-scale impacts of enhanced-solubilization solutions on organic-liquid morphology and fluid-fluid interfacial area with microtomography

Synchrotron X-ray microtomography was used to directly quantify the impact of enhanced-solubilization solutions on mobilization of a non-wetting fluid in a natural porous medium. Polyoxyethylene sorbitan monooleate (tween 80), hydroxypropyl-β-cyclodextrin (HPCD), sodium dodecyl sulfate (SDS), and ethanol were used as the solubilization reagents. Tetrachloroethene (PCE) was used as the entrapped organic immiscible liquid. Microtomography images were collected prior to and after successive floods with each reagent solution. The total trapping number was varied by more than three orders of magnitude by changing the velocity and viscosity of the aqueous phase as well as the interfacial tension between fluid phases. The results were compared to those obtained from equivalent experiments conducted with water flooding. Inspection of images indicated that mobilization caused a change in PCE morphology and distribution for all reagents with the exception of HPCD, for which a lower trapping number was determined compared to the other experiments. Mobilization of organic liquid was observed at total trapping numbers that were consistent with prior work. The enhanced-solubilization solutions had minimal impact on the configuration or magnitude of fluid-fluid interfacial area.

Interfacial area↗
DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗
DOE OSTI2016

Critical shear rate and torque stability condition for a particle resting on a surface in a fluid flow

In this work, we advance a quantitative description of the critical shear rate $\dot{\gamma}_{c}$ needed to dislodge a spherical particle resting on a surface with a model asperity in laminar and turbulent fluid flows. We have built a cone-plane experimental apparatus which enables measurement of $\dot{\gamma}_{c}$ over a wide range of particle Reynolds number Re p from 10–3 to 1.5 × 10 3 . The condition to dislodge the particle is found to be consistent with the torque balance condition after including the torque component due to drag about the particle center. The data for Re p < 0.5 is in good agreement with analytical calculations of the drag and lift coefficients in the Re p → 0 limit. For higher Rep, where analytical results are unavailable, the hydrodynamic coefficients are found to approach a constant for Re p > 1000. We show that a linear combination of the hydrodynamic coefficients found in the viscous and inertial limits can describe the observed $\dot{\gamma}_{c}$ as a function of the particle and fluid properties.

02 PETROLEUM↗
DOE OSTI2022

Fe Oxidation and Species Distribution at the Rock–Fluid Interface of Marcellus Shale Reacted with Hydraulic Fracturing Fluid

Hydraulic fracturing of shale reservoirs resulted in significant opportunity for increased oil and gas production in the United States. Rock-fluid interactions can cause mineral dissolution and precipitation reactions that lead to permeability changes in the shale matrix, which ultimately may affect transport pathways and hydrocarbon production. Understanding the distribution of secondary precipitates, such as barite and Fe(III) (hydro)oxides, and cation leaching at the rock-fluid interface is an important step to further investigate how these geochemical processes can change permeability and transport pathways. In this study, thin sections of the fracture-matrix interface were made from reacted Marcellus shale cores. The thin sections were characterized using synchrotron X-ray fluorescence imaging and synchrotron X-ray absorption spectroscopy. Fe species with different oxidation states were identified in the maps, together with barite and Ca distribution. The results show that ferrihydrite, as newly formed Fe(III)-bearing precipitates, aligned well with the border of the Ca (e.g., calcite) leaching region in the reaction front. Some Fe-containing clay also dissolved, but the dissolution region for the clay was not as deep as the calcite. Further, the reaction front is about three times deeper in the direction parallel to the shale bedding than that perpendicular to the bedding. The Ca leaching region can be an index for reaction front detection for Marcellus shale. Reactive transport modeling was conducted and the predicted Ca leaching boarder align well with ferrihydrite precipitation, consistent with the experimental observation. The carbonate mineral dissolution can be crucial to promote fluid access into the shale matrix. Together with our previous study on the shale reactive surface, this follow-up study showed similar Ca leaching region and Fe(III) precipitates distribution in the reaction front regardless of barite precipitation on the surface, indicating that the barite coatings on the surface may not pose a significant impact on reactive transport at the shale-fluid interface.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2015

Two-dimensional hybrid Monte Carlo–fluid modelling of dc glow discharges: Comparison with fluid models, reliability, and accuracy

Two-dimensional hybrid Monte Carlo–fluid numerical code is developed and applied to model the dc glow discharge. The model is based on the separation of electrons into two parts: the low energetic (slow) and high energetic (fast) electron groups. Ions and slow electrons are described within the fluid model using the drift-diffusion approximation for particle fluxes. Fast electrons, represented by suitable number of super particles emitted from the cathode, are responsible for ionization processes in the discharge volume, which are simulated by the Monte Carlo collision method. Electrostatic field is obtained from the solution of Poisson equation. The test calculations were carried out for an argon plasma. Main properties of the glow discharge are considered. Current-voltage curves, electric field reversal phenomenon, and the vortex current formation are developed and discussed. The results are compared to those obtained from the simple and extended fluid models. Contrary to reports in the literature, the analysis does not reveal significant advantages of existing hybrid methods over the extended fluid model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The island coalescence problem: Scaling of reconnection in extended fluid models including higher-order moments

As modeling of collisionless magnetic reconnection in most space plasmas with realistic parameters is beyond the capability of today's simulations, due to the separation between global and kinetic length scales, it is important to establish scaling relations in model problems so as to extrapolate to realistic scales. Recently, large scale particle-in-cell simulations of island coalescence have shown that the time averaged reconnection rate decreases with system size, while fluid systems at such large scales in the Hall regime have not been studied. Here, we perform the complementary resistive magnetohydrodynamic (MHD), Hall MHD, and two fluid simulations using a ten-moment model with the same geometry. In contrast to the standard Harris sheet reconnection problem, Hall MHD is insufficient to capture the physics of the reconnection region. Additionally, motivated by the results of a recent set of hybrid simulations which show the importance of ion kinetics in this geometry, we evaluate the efficacy of the ten-moment model in reproducing such results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hidden fluid mechanics: Learning velocity and pressure fields from flow visualizations

For centuries, flow visualization has been the art of making fluid motion visible in physical and biological systems. Although such flow patterns can be, in principle, described by the Navier-Stokes equations, extracting the velocity and pressure fields directly from the images is challenging. We addressed this problem by developing hidden fluid mechanics (HFM), a physics-informed deep-learning framework capable of encoding the Navier-Stokes equations into the neural networks while being agnostic to the geometry or the initial and boundary conditions. We demonstrate HFM for several physical and biomedical problems by extracting quantitative information for which direct measurements may not be possible.

Raissi, Maziar↗