Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flows”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Velocity Profiling of a Gas–Solid Fluidized Bed Using Electrical Capacitance Volume Tomography

In this study, a method of producing velocity profile maps from electrical capacitance volume tomography (ECVT) measurements by reconstructing displacement from measured changes in capacitance is developed and applied to fluidized bed systems. The mapping of the reconstruction leverages the gradient of the sensitivity distribution of the ECVT sensor to circumvent the need for image cross correlation techniques. Experimental data of both bubbling and slugging fluidized beds are collected in a cold flow model. Adaptation of the technique is discussed in detail, and velocity profiles are obtained for a range of gas flow rates. The produced velocity maps are compared against the established methods of cross correlation and against empirical correlations from the literature and are found to agree well in tracking slug and bubble velocity. The exception is when the tracked object is large relative to the ECVT sensor dimensions, a scenario that can be avoided through proper sensor design. The quantities of average velocity, momentum, and solid and gas volume fraction are derived from the image and velocity profiles. The results demonstrate and extend the power of ECVT as a measurement tool for the study and monitoring of gas–solid fluidized beds by providing a computationally cheaper alternative to 3-D cross correlation for deriving velocity profiles.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

Uniaxial Creep as a Control on Mercury Intrusion Capillary Pressure in Consolidating Rock Salt

The nature of geologic disposal of nuclear waste in salt formations requires validated and verified two-phase flow models of transport of brine and gas through intact, damaged, and consolidating crushed salt. Such models exist in other realms of subsurface engineering for other lithologic classes (oil and gas, carbon sequestration etc. for clastics and carbonates) but have never been experimentally validated and parameterized for salt repository scenarios or performance assessment. Models for waste release scenarios in salt back-fill require phenomenological expressions for capillary pressure and relative permeability that are expected to change with degree of consolidation, and require experimental measurement to parameterize and validate. This report describes a preliminary assessment of the influence of consolidation (i.e. volume strain or porosity) on capillary entry pressure in two phase systems using mercury injection capillary pressure (MICP). This is to both determine the potential usefulness of the mercury intrusion porosimetry method, but also to enable a better experimental design for these tests. Salt consolidation experiments are performed using novel titanium oedometers, or uniaxial compression cells often used in soil mechanics, using sieved run-of-mine salt from the Waste Isolation Pilot Plant (WIPP) as starting material. Twelve tests are performed with various starting amounts of brine pore saturation, with axial stresses up to 6.2 MPa (~900 psi) and temperatures to 90°C. This corresponds to UFD Work Package 15SN08180211 milestone “FY:15 Transport Properties of Run-of-Mine Salt Backfill – Unconsolidated to Consolidated”. Samples exposed to uniaxial compression undergo time-dependent consolidation, or creep, to various degrees. Creep volume strain-time relations obey simple log-time behavior through the range of porosities (~50 to 2% as measured); creep strain rate increases with temperature and applied stress as expected. Mercury porosimetry is used to determine characteristic capillary pressure curves from a series of consolidation tests and show characteristic saturation-capillary pressure curves that follow the common van Genuchten (1978, 1980) formulation at low stresses. Higher capillary pressure data are suspect due to the large potential for sample damage, including fluid inclusion decrepitation and pore collapse. Data are supportive of use of the Leverett “J” function (Leverett, 1941) to use for scaling characteristic curves at different degrees of consolidation, but better permeability determinations are needed to support this hypothesis. Recommendations for further and refined testing are made with the goal of developing a self- consistent set of constitutive laws for granular salt consolidation and multiphase (brine-air) flow.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Development of a conduction-based model for analyzing frozen startup of alkali-metal heat pipes

One key area of interest in heat pipe modeling/simulation is to analyze the startup behavior of the liquid-metal heat pipes (LMHPs) from a frozen state. This so-called ‘frozen startup’ process involves a complex set of nonlinear mass and heat transport phenomena, including phase transitions from solid to liquid and vapor, multiphase interactions, microporous wick flow, and compressible vapor dynamics. The complexity of these processes makes it challenging to simulate LMHP’s frozen startup using conventional numerical methods or commercial computational fluid dynamics (CFD) software. This paper presents a simplified conduction-based modeling approach that can provide practical insights into the entire LMHP frozen startup process, while alleviating the challenges of modeling its complex physics. The theoretical foundation and physical assumptions of the proposed model are based solely on heat-conduction equation, allowing for a more tractable simulation without sacrificing essential physical accuracy. The proposed model was implemented in a commercial CFD software, and its prediction was compared with the experimental data obtained from sodium heat-pipe startup experiments. The comparison highlights the proposed model's ability to capture the transient thermal behavior of LMHP during frozen startup. This study not only validates the conduction-based frozen startup modeling method but also shows its potential as a practical and efficient tool for understanding the startup performance of the LMHP systems.

Microreactor↗
DOE OSTI2021

Critical analysis of velocimetry methods for particulate flows from synthetic data

Particle tracking methods that extract high-fidelity particle velocity data from high speed video of particle laden flows is a common experimental technique applied to chemical processes. These measurements are used to better understand the motion of particles and fluids in complex systems and create data against which computational models are validated. However, the methods, codes, and experimental setups all have limitations. It is imperative that practitioners verify the methods and their implementation as well as understand the limitations of experimental setups. This work focuses on quantifying the visible depth of field in a high particle concentration fluidized bed. Following a precedent set by the particle imaging velocimetry community, a particle velocity field is manufactured using a computational fluid dynamics and discrete element method simulation. Photo realistic high-speed videos are rendered based on the simulated data using the three-dimensional creation software Blender. Particle velocities are extracted from the synthetic high-speed videos using three variants of Particle Tracking Velocimetry and Optical Flow Velocimetry methodologies. Here, the tracked results are then compared to the known solution, quantifying the error associated with the assumed visible depth. The results indicate that at depth of one particle diameter, all three particle tracking codes give accurate measurements, largely within 5%. However, the error increases when the full bed video measurements are compared to the known solution at one particle diameter, i.e., mimicking a validation study. Finally, for some statistics the constant depth assumption only increases the error slightly, for others significantly.

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

Coupled thermal–hydrological–mechanical modeling of CO 2 -enhanced coalbed methane recovery

CO 2 -enhanced coalbed methane recovery, also known as CO 2 -ECBM, is a potential win-win approach for enhanced methane production while simultaneously sequestering injected anthropogenic CO 2 to decrease CO 2 emissions into the atmosphere. Here, CO 2 -ECBM is simulated using a coupled thermal–hydrological–mechanical (THM) numerical model that considers multiphase (gas and water) flow and solubility, multicomponent (CO 2 and CH 4 ) diffusion and adsorption, heat transfer and coal deformation. The coupled model is based on the TOUGH-FLAC simulator, which is applied here for the first time to model CO 2 -ECBM. The capacity of the simulator for modeling methane production is verified by a code-to-code comparison with the general-purpose finite-element solver COMSOL. Then, the TOUGH-FLAC simulator is applied in an isothermal simulation to study the variations in permeability evolution during a CO 2 -ECBM operation while considering four different stress-dependent permeability models that have been implemented into the simulator. Finally, the TOUGH-FLAC simulator is applied in non-isothermal simulations to model THM responses during a CO 2 -ECBM operation.Our simulations show that the permeability evolution, mechanical stress, and deformation are all affected by changes in pressure, temperature and adsorption swelling, with adsorption swelling having the largest effect. The calculated stress changes do not induce any mechanical failure in the coal seam, except near the injection well in one case of a very unfavorable stress field.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Insights from high-fidelity modeling of industrial rotary bell atomization

The global automotive industry sprayed over 2.6 billion liters of paint in 2018, much of which through electrostatic rotary bell atomization, a highly complex process involving the fluid mechanics of rapidly rotating thin films tearing apart into micrometer-thin filaments and droplets. Coating operations account for 65% of the energy usage in a typical automotive assembly plant, representing 10,000s of gigawatt-hours each year in the United States alone. Optimization of these processes would allow for improved robustness, reduced material waste, increased throughput, and significantly reduced energy usage. Here, we introduce a high-fidelity mathematical and algorithmic framework to analyze rotary bell atomization dynamics at industrially relevant conditions. Our approach couples laboratory experiment with the development of robust non-Newtonian fluid models; devises high-order accurate numerical methods to compute the coupled bell, paint, and gas dynamics; and efficiently exploits high-performance supercomputing architectures. These advances have yielded insight into key dynamics, including i) parametric trends in film, sheeting, and filament characteristics as a function of fluid rheology, delivery rates, and bell speed; ii) the impact of nonuniform film thicknesses on atomization performance; and iii) an understanding of spray composition via primary and secondary atomization. These findings result in coating design principles that are poised to improve energy- and cost-efficiency in a wide array of industrial and manufacturing settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

HYDRATE v1.5 OPTION OF TOUGH+ v1.5

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis and Pruess, 2014], a successor to the TOUGH2 [Pruess et al., 1999, 2012] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy's law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

Moridis, GeorgeJ↗
DOE OSTI2020

PFLOTRAN

PFLOTRAN is a GNU LGPL licensed code that simulates single phase water and multiphase (air, co2, water) flow, energy transport, geochemical reaction and biological reaction in the earth's subsurface environment. The problem may be discretized on structured and unstructured grids. The code is designed to run in parallel on supercomputers and leverages the open source PETSc library for solvers and data structures and open source HDF5 for binary file I/O. See www.pflotran.org and/or doc-dev.pflotran.org for further information. SAND2020-12161 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Hammond, Glenn↗
DOE OSTI2014

User's Manual for the Hydrate v1.5 Option of TOUGH+ v1.5: A Code for the Simulation of System Behavior in Hydrate-Bearing Geologic Media

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy’s law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

User's Manual for the RealGasBrine v1.0 Option of TOUGH+ v1.5: A Code for the Simulation of System Behavior in Gas-Bearing Geologic Media

REALGASBRINE v1.0 is a numerical code that for the simulation of the behavior of gas-bearing porous and/fractured geologic media. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. REALGASBRINE v1.0 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. REALGASBRINE v1.0 describes the non-isothermal two- (for pure water) or three-phase (for brine) flow of an aqueous phase and a real gas mixture in a gas-bearing medium, with a particular focus in ultra-tight (such as tight-sand and shale gas) systems. Up to 12 individual real gases can be tracked, and salt can precipitate as solid halite. The capabilities of the code include coupled flow and thermal effects, real gas behavior, Darcy and non-Darcy flow, several isotherm options of gas sorption onto the grains of the porous media, complex fracture descriptions, gas solubility into water, and geomechanical effects on flow properties. REALGASBRINE v1.0 allows the study of flow and transport of fluids and heat over a wide range of time frames and spatial scales not only in gas reservoirs, but also in any problem involving the flow of gases in geologic media, including the geologic storage of greenhouse gas mixtures, the behavior of geothermal reservoirs with multi-component condensable (H 2 O and CO 2 ) and non-condensable gas mixtures, the transport of water and released H 2 in nuclear waste storage applications, etc.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Uncertainty Quantification Analysis of Both Experimental and CFD Simulation Data of a Bench-scale Fluidized Bed Gasifier

Adequate assessment of the uncertainties in modeling and simulation is becoming an integral part of the simulation based engineering design. The goal of this study is to demonstrate the application of non-intrusive Bayesian uncertainty quantification (UQ) methodology in multiphase (gas-solid) flows with experimental and simulation data, as part of our research efforts to determine the most suited approach for UQ of a bench scale fluidized bed gasifier. UQ analysis was first performed on the available experimental data. Global sensitivity analysis performed as part of the UQ analysis shows that among the three operating factors, steam to oxygen ratio has the most influence on syngas composition in the bench-scale gasifier experiments. An analysis for forward propagation of uncertainties was performed and results show that an increase in steam to oxygen ratio leads to an increase in H2 mole fraction and a decrease in CO mole fraction. These findings are in agreement with the ANOVA analysis performed in the reference experimental study. Another contribution in addition to the UQ analysis is the optimization-based approach to guide to identify next best set of additional experimental samples, should the possibility arise for additional experiments. Hence, the surrogate models constructed as part of the UQ analysis is employed to improve the information gain and make incremental recommendation, should the possibility to add more experiments arise. In the second step, series of simulations were carried out with the open-source computational fluid dynamics software MFiX to reproduce the experimental conditions, where three operating factors, i.e., coal flow rate, coal particle diameter, and steam-to-oxygen ratio, were systematically varied to understand their effect on the syngas composition. Bayesian UQ analysis was performed on the numerical results. As part of Bayesian UQ analysis, a global sensitivity analysis was performed based on the simulation results, which shows that the predicted syngas composition is strongly affected not only by the steam-to-oxygen ratio (which was observed in experiments as well) but also by variation in the coal flow rate and particle diameter (which was not observed in experiments). The carbon monoxide mole fraction is underpredicted at lower steam-to-oxygen ratios and overpredicted at higher steam-to-oxygen ratios. The opposite trend is observed for the carbon dioxide mole fraction. These discrepancies are attributed to either excessive segregation of the phases that leads to the fuel-rich or -lean regions or alternatively the selection of reaction models, where different reaction models and kinetics can lead to different syngas compositions throughout the gasifier. To improve quality of numerical models used, the effect that uncertainties in reaction models for gasification, char oxidation, carbon monoxide oxidation, and water gas shift will have on the syngas composition at different grid resolution, along with bed temperature were investigated. The global sensitivity analysis showed that among various reaction models employed for water gas shift, gasification, char oxidation, the choice of reaction model for water gas shift has the greatest influence on syngas composition, with gasification reaction model being second. Syngas composition also shows a small sensitivity to temperature of the bed. The hydrodynamic behavior of the bed did not change beyond grid spacing of 18 times the particle diameter. However, the syngas concentration continued to be affected by the grid resolution as low as 9 times the particle diameter. This is due to a better resolution of the phasic interface between the gases and solid that leads to stronger heterogeneous reactions. This report is a compilation of three manuscripts published in peer-reviewed journals for the series of studies mentioned above.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2018

Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation

Tight oil formations such as those in the Bakken petroleum system are known to hold hundreds of billions of barrels of oil in place; however, the primary recovery factor for these plays is typically less than 10%. Tight oil formations, including the Bakken Formation, therefore, may be attractive candidates for enhanced oil recovery (EOR) using CO 2 . Multiphase fluid behavior and flow in fluid-rich shales can vary substantially depending on the size of pore throats, and properties such as fluid viscosity and density are much different in nanoscale pores than in macroscale pores. Thus it is critical to understand the nature and distribution of nano-, micro-, and macroscale pores and fracture networks. To address these issues, the Energy & Environmental Research Center (EERC) has been conducting a research program entitled “Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation.” The objectives of the project are 1) the use of advanced characterization methods to better understand and quantify the petrophysical and geomechanical factors that control CO 2 and oil mobility within tight oil formation samples, 2) the determination of CO 2 permeation and oil extraction rates in tight reservoir rocks and organic-rich shales of the Bakken, and 3) the integration of the laboratory-based CO 2 permeation and oil extraction data and the characterization data into geologic models and dynamic simulations to develop predictions of CO 2 storage resource and EOR in the Bakken tight oil formation. A combination of standard and advanced petrophysical characterization techniques were applied to characterize samples of Bakken Formation tight reservoir rock and shales from multiple wells. Techniques included advanced computer tomography (CT) imaging, scanning electron microscopy (SEM) techniques, whole-core and micro x-ray CT imaging, field emission (FE) SEM, and focused ion beam (FIB) SEM. Selected samples were also analyzed for geomechanical properties. X-ray CT imaging yielded information on the occurrence of fractures, bedding planes, fossils, and bioturbation in core, as well as data on bulk density and photoelectric factor logs, which were used to interpret porosity, organic content, and mineralogy. FESEM was used for characterization of nano- and microscale features, including nanoscale pore visualization and micropore and pore throat mineralogy. FIBSEM yielded micro- to nanoscale visualization of fracture networks, porosity and pore-size distribution, connected versus isolated porosity, and distribution of organics. Results from the characterization activities provide insight on nanoscale fracture properties, pore throat mineralogy and connectivity, rock matrix characteristics, mineralogy, and organic content. Laboratory experiments demonstrated that CO 2 can permeate the tight matrix of Bakken shale and nonshale reservoir samples and mobilize oil from those samples. Geologic models were created at scales ranging from the core plug to the reservoir, and dynamic simulations were conducted. The data from the characterization and laboratory-based activities were integrated into modeling research activities to determine the fundamental mechanisms controlling fluid transport in the Bakken, which support EOR scheme design and estimation of CO 2 storage potential in tight oil formations. Simulation results suggest a CO 2 storage resource estimate range of 169 million to 1.5 billion tonnes for the Bakken in North Dakota, possibly resulting in 1.8 billion to 16 billion barrels of incremental oil.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

A CONDUCTION-BASED HEAT PIPE MODEL FOR ANALYZING THE ENTIRE PROCESS OF LIQUID-METAL HEAT PIPE STARTUP

The Heat Pipe-cooled Microreactor (HPM) is one of the micro nuclear reactor designs under active study at the U.S. Idaho National Laboratory. Among the major concerns of HPM research is to understand the startup behavior of the heat pipe-cooled system associated with the startup of the liquid-metal heat pipes initialing from frozen state. The startup of liquid-metal heat pipes typically involves a number of nonlinear mass and heat transport processes including the phase change from solid to liquid and vapor. Hence, it is still a huge challenge to simulate the liquid-metal heat pipe startup using conventional CFD methods and software. The major difficulties of numerical CFD modeling come from the phase-change process, multiphase interaction, microporous wick flow, and compressible gas dynamics that occur during startup of the liquid-metal heat pipes. This paper proposes a simplified conduction-based method to provide practical insights into the entire startup process of the liquid-metal heat pipes while mitigating the challenges of addressing all the complex physics. We discuss the theoretical basis and modeling assumptions to analyze the liquid-metal heat pipe startup from frozen state based solely on heat-conduction equations. Then, the proposed model is implemented into the commercial CFD software to verify the model performance. The model prediction results are discussed via the comparison with the experimental data obtained from sodium heat-pipe startup experiments.

42 ENGINEERING↗
DOE OSTI2018

Multiphase Large Eddy Simulation of a Pulsed Sieve-Plate Extraction Column - 18125

Presented here is research resulting from a framework of single-phase and multiphase modelling trials. These are intended to accurately simulate the internal flow phenomena arising in pulsed sieve-plate extraction columns representative of those found in industry. Three-dimensional unsteady flow calculations using large eddy simulation coupled with dynamic sub-grid scale modelling have been used to capture the transient effects of sinusoidal flow conditions using the open source CFD code OpenFOAM{sup R}. Multiphase interface interactions are modelled using the volume of fluid method with appropriate heavy and light phase fluidic properties capturing surface tension, density and viscous effects. The system studied consists of a 150 mm diameter column with two perforated plates with a fractional free area of 25 %, operational under 1 Hz pulsing frequency at a 50 mm amplitude and plate spacing of 300 mm. The materials simulated include water (single-phase) and 3 M nitric acid and dodecane / tributyl phosphate mixture at 30 vol% (multiphase). Results are compared with appropriate correlations from literature. The dispersive mixing efficiency of the system was also evaluated via a mixing index value which allows quantification of the ratio of rotational and irrotational flow components. This study has found that the predictive methods employed are effective at capturing the different scales of turbulence present within pulsed sieve-plate extraction columns, and its influence on their operation. As such, the work highlights the ability of large eddy simulation-based modelling to predict the complex flows in such column geometries. Furthermore, the techniques employed provide valuable insight into the hydrodynamic flow characteristics and mixing mechanisms arising within pulsed sieve-plate extraction columns. (authors)

42 ENGINEERING↗
DOE OSTI2015

On the use of Darcy's law and invasion-percolation approaches for modeling large-scale geologic carbon sequestration

Most large-scale flow and transport simulations for geologic carbon sequestration (GCS) applications are carried out using simulators that solve flow equations arising from Darcy's law. Recently, the computational advantages of invasion-percolation (IP) modeling approaches have been presented. We show that both the Darcy's-law- and the gravity-capillary balance solved by IP approaches can be derived from the same multiphase continuum momentum equation. More specifically, Darcy's law arises from assuming creeping flow with no viscous momentum transfer to stationary solid grains, while it is assumed in the IP approach that gravity and capillarity are the dominant driving forces in a quasi-static two-phase (or more) system. There is a long history of use of Darcy's law for large-scale GCS simulation. However, simulations based on Darcy's law commonly include significant numerical dispersion as users employ large grid blocks to keep run times practical. In contrast, the computational simplicity of IP approaches allows large-scale models to honor fine-scale hydrostratigraphic details of the storage formation which makes these IP models suitable for analyzing the impact of small-scale heterogeneities on flow. However, the lack of time-dependence in the IP models is a significant disadvantage, while the ability of Darcy's law to simulate a range of flows from single-phase- and pressure-gradient-driven flows to buoyant multiphase gravity-capillary flow is a significant advantage. We believe on balance that Darcy's law simulations should be the preferred approach to large-scale GCS simulations.

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

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

The Plant Water Management Experiments: Soil

A simple means of watering plants in the low-g environment aboard orbiting spacecraft is not obvious. Since the beginning of spaceflight, numerous approaches have been pursued to water plants that seek to maximize plant viability and system reliability, while minimizing crew time and system complexity. We are not there yet. The Plant Water Management (PWM) Soil experiments seek to apply recent advances in low-g capillary fluidics phenomena to the challenges faced by plant growth operations aboard spacecraft. The primary challenge is to establish earth-like flows minimizing low-g specific adaptations required of the plants. This is difficult due to the ever-present fluid physics challenges of poorly-wetting multiphase inertial-visco-capillary flows in geometrically complex conduits and containers. In this paper, we present recent flight results for the PWM Soil experiments where arcillite ‘soil reservoirs’ are arranged in a non-wetting host soil that serves as an O2-breathing wetting barrier. In this way, a largely terrestrial water-soil environment is mimicked where, as liquid is evapo-transpired through the growing plant foliage, the effective water table passively ‘falls’ reducing viscous lengths and increasing water uptake for the plant. We present data from 6 days of 24-7 experiments on the ISS testing 3 different plant root models. We also present and correlate a capillary flow model which captures the primary features of the flow. Our summary is valued for the assessment of current and future low-g plant watering systems employing soil media.

microgravity↗
DOE OSTI2017

The Fascinating and Complex Dynamics of Geyser Eruptions

Geysers episodically erupt liquid and vapor. Despite two centuries of scientific study, basic questions persist-why do geysers exist? What determines eruption intervals, durations, and heights? What initiates eruptions? Through monitoring eruption intervals, analyzing geophysical data, taking measurements within geyser conduits, performing numerical simulations, and constructing laboratory models, some of these questions have been addressed. Geysers are uncommon because they require a combination of abundant water recharge, magmatism, and rhyolite flows to supply heat and silica, and large fractures and cavities overlain by low-permeability materials to trap rising multiphase and multicomponent fluids. Eruptions are driven by the conversion of thermal to kinetic energy during decompression. Larger and deeper cavities permit larger eruptions and promote regularity by isolating water from weather variations. Furthermore, the ejection velocity may be limited by the speed of sound of the liquid + vapor mixture.

58 GEOSCIENCES↗