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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Developing a system for in vivo imaging of maize roots containing iodinated contrast media in soil using synchrotron XCT and XRF

We sought to develop a novel experimental system which enabled application of iodinated contrast media to in vivo plant roots intact in soil and was compatible with time-resolved synchrotron X-ray computed tomography imaging. The system was developed to overcome issues of low contrast to noise within X-ray computed tomography images of plant roots and soil environments, the latter of which can complicate image processing and result in the loss of anatomical information. To demonstrate the efficacy of the system we employ the novel use of both synchrotron X-ray computed tomography and synchrotron X-ray fluorescence mapping to capture the translocation of the contrast media through root vasculature into the leaves. With the application of contrast media we identify fluid flow in root vasculature and visualise anatomical features, which are otherwise often only observable in ex vivo microscopy, including: the xylem, metaxylem, pith, fibres in aerenchyma and leaf venation. We are also able to observe interactions between aerenchyma cross sectional area and solute transport in the root vasculature with depth. Our novel system was capable of successfully delivering sufficient contrast media into root and leaf tissues such that anatomical features could be visualised and internal fluid transport observed. We propose that our system could be used in future to study internal plant transport mechanisms and parameterise models for fluid flow in plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗
DOE OSTI2023

Measuring higher-order moments of neutron-time-of-flight data for cryogenic inertial confinement fusion implosions on OMEGA

Ion temperatures serve as an important diagnostic for inertial confinement fusion (ICF) implosions. In direct-drive ICF experiments on OMEGA, neutron-time-of-flight (nTOF) data are used to infer the ion temperature of the fusing plasma produced in the implosion experiment. The analysis of the nTOF data requires an assumption about the shape of the underlying source signal. Since the source nTOF signal is a near-replica of the neutron energy spectrum, an ideal Gaussian shape, corresponding to the neutron energy spectrum of a uniform temperature plasma, is routinely employed. However, spatial and temporal variations of the ion temperature in the plasma give rise to higher-order moments, which were first described by Munro [Nucl. Fusion 56, 036001 (2016)]. In this work, we show a simpler alternative analysis to derive moments of the neutron energy spectrum for a plasma with variations in ion temperature. We also present a revised analysis of measured nTOF signals that uses a model with an additional degree of freedom to take into account the effect of ion temperature variations on the shape of the spectrum. Compared to presently used nTOF analysis, the revised analysis yields on average ≈2× more accurate fits to the data and up to 15% higher ion temperatures for cryogenic experiments. As a result, we quantify the ion temperature inflation caused by radially symmetric fluid flows, which are present even in a symmetric implosion, and which serve as a lower bound on the ion temperature inflation in real implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.

42 ENGINEERING↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

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

Optimal renormalization of multiscale systems

While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.

Price, Jacob↗
DOE OSTI2021

Acoustic velocity and permeability of acidized and propped fractures in shale

We report that from geochemical reactions to proppant emplacement, hydraulic fracturing induces various chemomechanical fracture alterations in shale reservoirs. Hydraulic fracturing through the injection of a vast amount and variety of fluids and proppants has substantial impacts on fluid flow and hydrocarbon production. There is a strong need to improve our understanding on how fracture alterations affect flow pathways within the stimulated rock volume and develop monitoring tools. We have conducted time-lapse rock-physics experiments on clay-rich (carbonate-poor) Marcellus shales to characterize the acoustic velocity and permeability responses to fracture acidizing and propping. Acoustic P- and S-wave velocities and fracture permeability were measured before and after laboratory-induced fracture alterations along with microstructural imaging through X-ray computed tomography and scanning electron microscopy. Our experiments indicate that the S-wave velocity is an important geophysical observable, particularly the S-wave polarized perpendicular to fractures because it is sensitive to fracture stiffness. The acidizing and propping of a fracture decrease its elastic stiffness. This effect is stronger for acidizing, so it is possible that proppant monitoring will be masked by chemical alteration except when propping is highly efficient (i.e., most fractures are propped). However, fracture permeability is undermined by the softening of fracture surfaces due to acidizing, while being greatly enhanced by propping. These contrasting effects on fluid flow in combination with similar seismic attributes indicate the importance of experiments to improve existing rock-physics models, which must include changes to the rock frame. Such improvements are necessary for a correct interpretation of seismic velocity monitoring of flow pathways in stimulated shales.

58 GEOSCIENCES↗
DOE OSTI2022

Multiple experimental studies of pore structure and mineral grain sizes of the Woodford shale in southern Oklahoma, USA

Pore structure study is an important part of unconventional shale reservoir characterization, since the pore system provides the primary petroleum storage space and fluid flow pathways. Previous studies have suggested that the pore structure is related to the total organic carbon (TOC) content, mineral compositions, and the maturity of the organic matter (OM). However, few studies have focused on the mineral grains, the primary grains being deposited but before cementation, which are the building blocks of shale. Eight Woodford Shale outcrop samples from southern Oklahoma were chosen to study the effects of mineral grain size on the pore structure characterization, using multiple and complementary experimental approaches, including laser diffraction, mineralogy, TOC, pyrolysis, liquid immersion porosimetry, mercury intrusion porosimetry, gas physisorption, (ultra) small angle X-ray scattering, scanning electron microscopy, and spontaneous imbibition. The results from different experiments of eight samples show that the Woodford Shale has the mean mineral grain diameters at 3–6 μm, a wide range of porosity at 3–40% and pore diameters at 50–1,000 nm, and various pore connectivity. Grain size variation was probably caused by the sea-level fluctuation during its deposition, which affect the porosity, pore size distribution, and pore connectivity. With decreasing mineral grain sizes, the porosity tends to increase while the pore connectivity worsens. The results also indicate that OM and carbonates in this low-maturity Woodford Shale could block the pores and decrease the porosity. Coupling with the grain size analyses, the control of depositional environment on grain sizes and subsequent effects on pore structure is identified. The pore structure characteristics over a wide pore-diameter range provided by multiple experiments could improve the understanding of storage space and fluid flow in the Woodford Shale to further increase its petroleum production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Comparative Study of Deep Learning Models for Fracture and Pore Space Segmentation in Synthetic Fractured Digital Rocks

This study focuses on the comparative study of deep learning (DL) models for pore space and discrete fracture networks (DFNs) segmentation in synthetic fractured digital rocks, specifically targeting low-permeability rock formations, such as shale and tight sandstones. Accurate characterization of pore space and DFNs is critical for subsequent property analysis and fluid flow modeling. Four DL models, SegNet, U-Net, U-Net-wide, and nested U-Net (i.e., U-Net++), were trained, validated, and tested using synthetic datasets, including input and label image pairs with varying properties. The model performance was assessed regarding pixel-wise metrics, including the F1 score and pixel-wise difference maps. In addition, the physics-based metrics were considered for further analysis, including sample porosity and absolute permeability. Particularly, We first simulated the permeability of porous media containing only pore space and then simulated the permeability of porous media with DFNs added. The difference between these two values is used to quantify the connectivity of segmented DFNs, which is an important parameter for low-permeability rocks. The pixel-wise metrics showed that the nested U-Net model outperformed the rest of the DL models in pore space and DFNs segmentation, with the SegNet model exhibiting the second-best performance. Particularly, nested U-Net enhanced segmentation accuracy for challenging boundary pixels affected by partial volume effects. The U-Net-wide model achieved improved accuracy compared to the U-Net model, which indicated the influence of parameter numbers. Similarly, nested U-Net has the closest match to the ground truth of physics-based metrics, including the porosity of pore space and DFNs, and the permeability difference quantifying the connectivity of DFNs. The findings highlight the effectiveness of DL models, especially the U-Net++ model with nested architecture and redesigned skip connections, in accurately segmenting pore spaces and DFNs, which are crucial for pore-scale fluid flow and transport simulation in low-permeability rocks.

Wang, Hongsheng↗
DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2023

Utah FORGE: Fault Reactivation Through Fluid Injection Induced Seismicity Laboratory Experiments

Included are results from shear reactivation experiments on laboratory faults pre-loaded close to failure and reactivated by the injection of fluid into the fault. The sample comprises a single-inclined-fracture (SIF) transecting a cylindrical sample of Westerly granite. All experiments are conducted at ambient temperature and follow a similar protocol: (i) application of confining stresses (3MPa) on the fault fully saturated with DI water, (ii) shear-mobilization through the increase of axial loading at a constant displacement rate until a post-peak steady-state condition is reached, (iii) reduction of axial loading and related shear stress to a prescribed fraction of the peak steady-state frictional strength (typically 60% to 90%, representing intermediate to high magnitudes) and (iv), fault reactivation triggered by a stepwise increase of pore pressure on the fault in 0.1 MPa increments held constant for 1-5 minutes. Mechanical data from three ISCO pumps connected to a Temco pressure vessel measure axial, confining, and fault-related parameters, including fluid pressure (kPa), fluid flow rate (mL/min), and axial displacement (mm). See included code for initial data analysis and visualization for select experiments. Resource names represent experiment numbers found in the "Read Me" file, which describes each experimental setup and parameters.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Sorption heat pump and control method

A method operates an absorption heat pump system, specifically the flow of hydronic cooling fluid through the condenser during system start-ups, or when the cooling fluid temperature is low. To minimize the time for an absorption heat pump to reach full cooling or heating capacity, it is desirable for the high side pressure to increase as fast as possible, and the low side pressure to decrease as fast as possible. Since the high side pressure is a function of the temperature of the refrigerant exiting the condenser, if the condenser cooling fluid temperature is low, the corresponding high side pressure will be low, which may not permit adequate working fluid flow rates from the high pressure side of the system to the low pressure side.

Garrabrant, Michael A.↗