Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Reduced order models for assessing CO 2 impacts in shallow unconfined aquifers

Risk assessment studies of potential CO 2 sequestration projects consider many factors, including the possibility of brine and/or CO 2 leakage from the storage reservoir. Detailed multiphase reactive transport simulations have been developed to predict the impact of such leaks on shallow groundwater quality; however, these simulations are computationally expensive and thus difficult to directly embed in a probabilistic risk assessment analysis. Here we present a process for developing computationally fast reduced-order models which emulate key features of the more detailed reactive transport simulations. A large ensemble of simulations that take into account uncertainty in aquifer characteristics and CO 2 /brine leakage scenarios were performed. Twelve simulation outputs of interest were used to develop response surfaces (RSs) using a MARS (multivariate adaptive regression splines) algorithm (Milborrow, 2015). A key part of this study is to compare different measures of ROM accuracy. We then show that for some computed outputs, MARS performs very well in matching the simulation data. The capability of the RS to predict simulation outputs for parameter combinations not used in RS development was tested using cross-validation. Again, for some outputs, these results were quite good. For other outputs, however, the method performs relatively poorly. Performance was best for predicting the volume of depressed-pH-plumes, and was relatively poor for predicting organic and trace metal plume volumes. We believe several factors, including the non-linearity of the problem, complexity of the geochemistry, and granularity in the simulation results, contribute to this varied performance. The reduced order models were developed principally to be used in probabilistic performance analysis where a large range of scenarios are considered and ensemble performance is calculated. We demonstrate that they effectively predict the ensemble behavior. But, the performance of the RSs is much less accurate when used to predict time-varying outputs from a single simulation. If an analysis requires only a small number of scenarios to be investigated, computationally expensive physics-based simulations would likely provide more reliable results. Finally, if the aggregate behavior of a large number of realizations is the focus, as will be the case in probabilistic quantitative risk assessment, the methodology presented here is relatively robust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Fines migration during CO 2 injection: Experimental results interpreted using surface forces

The South West Hub project is one of the Australian Flagship Carbon Capture and Storage projects located in the south-west of Western Australia. To evaluate the injectivity potential during the forthcoming full-scale CO 2 injection, we conducted three core-flooding experiments using reservoir core plugs from the well Harvey-1. We aimed to investigate in this paper whether the injection of CO 2 leads to fines migration and permeability reduction due to the relatively high kaolinite content (up to 13%) in the injection interval of the target formation (i.e. the Wonnerup Member of the Lesueur Formation). We imaged the core samples before flooding to verify the presence of kaolinite at the pore-scale using scanning electron microscopy (SEM). We also examined the pore network of the core plugs before and after the core-flooding experiments using Nuclear Magnetic Resonance (NMR). Moreover, to gain a better understanding of any kaolinite fines migration, we delineated surface force using two models based on Derjaguin-Landau-Verwey-Overbeek (denoted by DLVO) theory coupled hydrodynamic force: (1) sphere/flat model representing interaction between kaolinite/quartz, and (2) flat/flat model representing interaction between kaolinite/kaolinite. Our core-flooding experimental results showed that CO 2 /brine injection triggered moderate to significant reduction in the permeability of the core samples with a negligible porosity change. NMR measurements supported the core-flooding results, suggesting that the relatively large pores disappeared in favour of a higher proportion of the medium to small pores after flooding. The DLVO calculations showed that some kaolinite particles probably lifted off and detached from neighbouring kaolinite particles rather than quartz grains. Moreover, the modelling results showed that the kaolinite fines migration would not occur under normal reservoir multiphase flow conditions. This is not because of the low hydrodynamic force. It is rather because the geometries of the particles dominate their interplay. Finally and overall, both of the experimental and analytical modelling results point to the fines migration to be the most likely cause of the permeability impairment observed during core-flooding experiments.

58 GEOSCIENCES↗
DOE OSTI2019

The effect of CO 2 -induced dissolution on flow properties in Indiana Limestone: An in situ synchrotron X-ray micro-tomography study

The injection of CO 2 -rich fluids in carbonate rocks results in an evolution of the pore space, with consequent changes in the hydraulic properties of the reservoir; how these properties evolve, particularly for parameters relevant to multiphase flow e.g. P c (s), remains a topic of active research despite several decades of study. We have carried out an in situ synchrotron X-ray microtomography experiment to monitor pore structure evolution during dissolution of an Indiana Limestone core; the experiment involved flowing CO 2 -saturated water through the core for 36 h and resulted in 10 volumes corresponding to different temporal stages of the dissolution process. The injection parameters corresponded to the flow velocities expected near the well-bore region of a shallow aqueous CO 2 injection; fast flow rates with high reactant availability. Analysis of the tomographic data shows flow-enhanced dissolution i.e. channeling, and provides a time-resolved map of pore space alteration. Using the resulting 4D pore space volume, we modeled the evolution of capillary-pressure curves; this exercise demonstrates how pore structure evolution could impact the invasion and remobilization of non-wetting fluids, dramatically decreasing the entry pressure and the P C in some parts of the sample. The modeling of permeability, using a Stokes solver approach, quantified the relationship of porosity vs. permeability; we found that a modest increase in porosity, especially when the channeling system is more developed, greatly affects permeability. Furthermore, these results demonstrate how movement of CO 2 saturated brine near injected plumes might alter drainage dynamics near the plume boundary, thus leading to mobilization across subtle capillary barriers.

58 GEOSCIENCES↗
DOE OSTI2019

Impact of CO2 injection on wettability of coal at elevated pressure and temperature

Injection of carbon dioxide (CO 2 ) into coal seams has been demonstrated as an effective technology for enhanced methane recovery and CO 2 storage. However, the impacts of the geochemical reactions between CO 2 and coal on the wettability of coal pore surfaces, which influences immiscible multiphase displacement, are not yet well understood. We studied wettability alterations of coal surfaces resulting from reactions with gas, liquid and supercritical (sc)CO 2 under varied pressure (1–141 bar) and temperatures (˜25 - 60 °C) through measuring static and dynamic contact angles with anthracite coal plates. We found that reactions with gas CO 2 only slightly changed the wettability of coal surfaces from water-wet to intermediate-wet with static contact angles from ˜60° to 70°-90°. However, reactions with liquid and scCO 2 altered the coal surfaces to strongly CO 2 -wet, with the contact angles up to 115-180°. We also found that both static and dynamic contact angles increase significantly with increasing pressure. Temperature affects the contact angles reversely especially under supercritical pressure conditions. These relationships of contact angles with pressure and temperature may be explained by the CO 2 density dependence on pressure and temperature.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Testing the efficacy of active pressure management as a tool to mitigate induced seismicity

Manipulation of subsurface stresses through fluid injection/extraction may induce seismicity through modulation of overpressure along faults. Globally, most industrial operations do not induce significant earthquakes. However, the possibility of inducing moderate to large earthquakes is recognized as a significant risk for basin-scale carbon dioxide (CO 2 ) storage operations (Niass et al., 2017). One possible method to reduce the occurrence of induced seismicity is to maintain only small (or zero) fluid overpressure changes in the subsurface through extraction of formation brine. For active pressure management (APM), operators may employ different procedures including extraction of various brine volumes and different configurations of injection/extraction well locations. However, the impact each these operational procedures on induced seismicity must be understood in order to develop management strategies. Here, we describe a numerical modeling workflow to evaluate the efficacy of APM techniques, that when applied to a specific site may be used to inform operational management strategies. We couple the multiphase reservoir model, NUFT, with the earthquake simulator, RSQSim and simulate seismicity induced by CO 2 injection into a hypothetical brine saturated reservoir. Reservoir overpressure is modified by co-extraction of brine at various fractions of the injected CO 2 volume from five potential wells. Our goal is to understand how the well-pair configuration and the extracted brine volume impact the total number and the maximum event magnitude. Additionally, the performance of each scenario is also evaluated in terms of the volume of CO 2 that can be stored before an earthquake above a tolerance threshold is induced.

58 GEOSCIENCES↗
DOE OSTI2015

Lattice Boltzmann modeling of boiling heat transfer: The boiling curve and the effects of wettability

A hybrid thermal lattice Boltzmann (LB) model is presented to simulate thermal multiphase flows with phase change based on an improved pseudopotential LB approach (Li et al., 2013). The present model does not suffer from the spurious term caused by the forcing-term effect, which was encountered in some previous thermal LB models for liquid–vapor phase change. Using the model, the liquid–vapor boiling process is simulated. The boiling curve together with the three boiling stages (nucleate boiling, transition boiling, and film boiling) is numerically reproduced in the LB community for the first time. The numerical results show that the basic features and the fundamental characteristics of boiling heat transfer are well captured, such as the severe fluctuation of transient heat flux in the transition boiling and the feature that the maximum heat transfer coefficient lies at a lower wall superheat than that of the maximum heat flux. Moreover, the effects of the heating surface wettability on boiling heat transfer are investigated. It is found that an increase in contact angle promotes the onset of boiling but reduces the critical heat flux, and makes the boiling process enter into the film boiling regime at a lower wall superheat, which is consistent with the findings from experimental studies.

42 ENGINEERING↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2026

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗
DOE OSTI2020

Scaling of two-phase water-steam relative permeability and thermal fluxes in porous media

Two-phase water-steam flow conditions are frequently encountered in many engineering applications, including geothermal reservoirs. Although routine calculations are based on the multiphase Darcy’s law, the role of the topology of the flowing phases at the pore-scale is usually neglected in the estimation of relative permeabilities. Instead, the latter are frequently computed using empirical models like the Corey correlation. In this work, we first apply the model for relative permeabilities based on pore-scale flow regimes developed by Picchi and Battiato (2019), Relative permeability scaling from pore-scale flow regimes, Water Resour. Res. 55, 3215–3233, to scenarios typical of geothermal reservoirs and then extend it by deriving the scaling laws for the transmissibilities and the thermal properties as a function of temperature. First, we discuss the scaling behavior of normalized relative permeabilities in terms of viscosity ratio and capillary number of water-steam systems and, then, we provide a validation of the model against experimental data available in the literature. The model captures well the data trends collected in real 3D porous media. These results suggest that water-steam relative permeabilities follow the same scaling behavior of gas-liquid systems where the non-wetting phase is much less viscous than the wetting phase. Finally, we investigate the impact that relative permeabilities have on heat transfer rates at two-phase flow conditions and the scaling of mass and energy transmissibility and thermal properties of the mixture. An estimation of the exergy carried by the two-phase water-steam mixture is also included.

42 ENGINEERING↗
DOE OSTI2021

Vortex Flow and Cavitation in Liquid Injection: A Comparison between High-Fidelity CFD Simulations and Experimental Visualizations on Transparent Nozzle Replicas

Experimental instantaneous shadowgraph visualizations on transparent glass nozzle replicas of high-pressure fuel injectors have been used to validate a novel in-house high-fidelity LES-VOF multiphase solver, to study the evolution of vortex flow and fuel cavitation. Both experiments and simulations capture the formation of an unsteady vapor structure inside the nozzle volume, which is referred to as 'string-cavitation'; strings are found at the core of the recirculation zones. Furthermore, high-fidelity simulations provide a very detailed insight into the vortex generation in the injector nozzle; strings appear within the time scales that are relevant for fast injection events (on the order of 0.1 milliseconds) and, for the problem under consideration, their generation seems mostly related to the flow pattern in the sac. It is also shown that vortexes interact, merge till they disrupt and favor the temporary inception of shear cavitation.

3-phase LES-VOF solver↗
DOE OSTI2022

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗
DOE OSTI2023

A numerical study of an atomizing jet in a resonant acoustic field

By applying state-of-the-art, high-fidelity methods for compressible multiphase flows, this work features a parametric study of a turbulent liquid jet atomizing under the influence of standing acoustic waves with five different amplitude values. We perform these simulations in a semi-periodic domain to isolate interactions between acoustic and atomization processes and to emulate a pure liquid fuel spray placed at the pressure node of a resonant acoustic field. Characterization of the flow includes measurements of the mean dimensions, droplet number, and surface area of the liquid distribution, alongside qualitative depictions of the interface evolution, breakup events, and hydrodynamic fields. As acoustic forcing is provided to the flow, the added energy provokes atomization when shear would otherwise be insufficient for breakup. The magnitude of the forcing determines the rate of droplet production, which accelerates as the evolving liquid surface modifies the induced motion of the gas phase. At the nodal plane, the acoustic field produces an organizing effect on the liquid structures.

42 ENGINEERING↗
DOE OSTI2026

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗
DOE OSTI2017

A crystal plasticity model incorporating the effects of precipitates in superalloys: Application to tensile, compressive, and cyclic deformation of Inconel 718

An elasto-plastic polycrystal plasticity model is developed and applied to an Inconel 718 (IN718) superalloy that was produced by additive manufacturing (AM). The model takes into account the contributions of solid solution, precipitates shearing, and grain size and shape effects into the initial slip resistance. Non-Schmid effects and backstress are also included in the crystal plasticity model for activating slip. The hardening law for the critical resolved shear stress is based on the evolution of dislocation density. In using the same set of material and physical parameters, the model is compared against a suite of compression, tension, and large-strain cyclic mechanical test data applied in different AM build directions. We demonstrate that the model is capable of predicting the particularities of both monotonic and cyclic deformation to large strains of the alloy, including decreasing hardening rate during monotonic loading, the non-linear unloading upon the load reversal, the Bauschinger effect, the hardening rate change during loading in the reverse direction as well as plastic anisotropy and the concomitant microstructure evolution. It is anticipated that the general model developed here can be applied to other multiphase alloys containing precipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Development of a coupled thermo-hydro-mechanical model in discontinuous media for carbon sequestration

Geomechanical alteration of porous media is generally ignored for most shallow subsurface applications, whereas CO 2 injection, migration, and trapping in deep saline aquifers will be controlled by coupled multifluid flow, energy transfer, and geomechanical processes. The accurate assessment of the risks associated with potential leakage of injected CO 2 and the design of effective injection systems requires that we represent these coupled processes within numerical simulators. The objectives of this study were to develop a coupled thermal-hydro-mechanical model into a single software, and to examine the coupling of thermal, hydrological, and geomechanical processes for simulation of CO 2 injection into the subsurface for carbon sequestration. A numerical model is developed to couple nonisothermal multiphase hydrological and geomechanical processes for prediction of multiple interconnected processes for carbon sequestration in deep saline aquifers. The geomechanics model was based on Rigid Body-Spring Model (RBSM), one of the discrete methods to model discontinuous rock system. Poisson’s effect that was often ignored by RBSM was considered in the model. The simulation of large-scale and long-term coupled processes in carbon capture and storage projects requires large memory and computational performance. Global Array Toolkit was used to build the model to permit the high performance simulations of the coupled processes. The model was used to simulate a case study with several scenarios to demonstrate the impacts of considering coupled processes and Poisson’s effect for the prediction of CO 2 sequestration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A micromechanical hydro-mechanical-damage coupled model for layered rocks considering multi-scale structures

An anisotropic Hydro-Mechanical-Damage coupled model was established in the framework of thermodynamics for saturated layered rocks, which contain arbitrarily-distributed cracks of smaller scales and orientated bedding planes of much larger scales. Here, the proposed model considers the anisotropic damage growth and frictional sliding of cracks and shear sliding and dilatancy of bedding planes, as well as the interaction between cracks and bedding planes. Numerical implementation of the nonlinear coupled model was developed based on TOUGHREACT, which is a well-established simulator for multiphase fluid flow and reactive transport analysis. The proposed model was first validated with laboratory water injection test results on rock blocks, and then used to study the effects of injection rates, in situ stresses and bedding planes on fluid injection processes at the field scale, demonstrating the importance to consider the variation of multi-scale structures for better understanding the coupled Hydro-Mechanical responses in layered rocks.

58 GEOSCIENCES↗
DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗