Search NASA⌕ Search

SEARCH · Search NASA

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Co-optimization of CO 2 -EOR and Storage Processes under Geological Uncertainty

This paper presents an integrated numerical framework to co-optimize EOR and CO 2 storage performance in the Farnsworth field unit (FWU), Ochiltree County, Texas. The framework includes a field-scale compositional reservoir flow model, an uncertainty quantification model and a neural network optimization process. The reservoir flow model has been constructed based on the field geophysical, geological, and engineering data. A laboratory fluid analysis was tuned to an equation of state and subsequently used to predict the thermodynamic minimum miscible pressure (MMP). A history match of primary and secondary recovery processes was conducted to estimate the reservoir and multiphase flow parameters as the baseline case for analyzing the effect of recycling produced gas, infill drilling and water alternating gas (WAG) cycles on oil recovery and CO 2 storage. A multi-objective optimization model was defined for maximizing both oil recovery and CO 2 storage. The uncertainty quantification model comprising the Latin Hypercube sampling, Monte Carlo simulation, and sensitivity analysis, was used to study the effects of uncertain variables on the defined objective functions. Uncertain variables such as bottom hole injection pressure, WAG cycle, injection and production group rates, and gas-oil ratio among others were selected. The most significant variables were selected as control variables to be used for the optimization process. A neural network optimization algorithm was utilized to optimize the objective function both with and without geological uncertainty. The vertical permeability anisotropy (Kv/Kh) was selected as one of the uncertain parameters in the optimization process. The simulation results were compared to a scenario baseline case that predicted CO 2 storage of 74%. The results showed an improved approach for optimizing oil recovery and CO 2 storage in the FWU. The optimization process predicted more than 94% of CO 2 storage and most importantly about 28% of incremental oil recovery. The sensitivity analysis reduced the number of control variables to decrease computational time. A risk aversion factor was used to represent results at various confidence levels to assist management in the decision-making process. The defined objective functions were proved to be a robust approach to co-optimize oil recovery and CO 2 storage. The Farnsworth CO 2 project will serve as a benchmark for future CO 2 –EOR or CCUS projects in the Anadarko basin or geologically similar basins throughout the world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Enhanced energy conversion efficiency promoted by cavitation in gasoline direct injection

High-pressure direct fuel injection plays the most crucial role in energy conversion and improving engine combustion efficiency and emission. The optimization of turbulent and multiphase fuel injection has focused on controlling hydrodynamic parameters such as injection pressure. While the thermodynamic influence is often considered in the flash boiling situation, we inquire into how gasoline-type fuel’s hydro- and thermodynamic properties impact the injection dynamics by fuel-temperature induced cavitation. The turbulent and cavitating flows emanating from the direct-injection nozzle is visualized by ultrafast x-ray imaging with an unprecedented spatiotemporal resolution. The ultrafast liquid-fuel dynamics are dominated by injection pressure as well as fuel temperature through cavitation, an important thermodynamic parameter but often difficult to control in engine combustion. With the most direct and quantitative measurement, we discovered that the near-nozzle fuel-jet dynamics can be perfectly scaled by a single dimensionless parameter, cavitation number, particularly sensitive to the fuel temperature, in a wide operation range. This universal scaling shows that cavitation can be harnessed to elevate the pneumatic-hydraulic to kinetic energy conversion efficiency, critical for promoting fuel atomization and engine combustion performance. This enhancement effect will have even more impact on engine combustion using alternative low-emission fuels with higher saturated vapor pressure.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

BMINN: Learning chemical potentials and parameters from voltage data for multi-phase battery modeling

Free-energy landscapes and chemical potentials govern the dynamics of phase transitions, transport, and stability in functional materials, yet they remain experimentally inaccessible under realistic operating conditions. Here we introduce a Bayesian model-integrated neural network (BMINN) that embeds physics-based formulations of non-autonomous partial differential-algebraic equations into probabilistic learning. This approach reconstructs hidden thermodynamics directly from macroscopic current-voltage data, providing quantitative access to metastable states, staging transitions, and energy barriers without synchrotron probes. Demonstrated on lithium-graphite electrodes, BMINN recovers full Gibbs free-energy landscapes with fidelity validated against operando X-ray diffraction. The framework generalizes across dynamical regimes, enabling accurate voltage prediction, internal state estimation, and inference of governing parameters. Beyond batteries, BMINN exemplifies a broadly applicable strategy for learning missing physics in multiphase, non-equilibrium systems, offering a new pathway to uncover hidden thermodynamic functions across condensed matter and materials physics.

25 ENERGY STORAGE↗
DOE OSTI2018

Carbonate stability in the reduced lower mantle

Carbonate minerals are important hosts of carbon in the crust and mantle with a key role in the transport and storage of carbon in Earth’s deep interior over the history of the planet. Whether subducted carbonates efficiently melt and break down due to interactions with reduced phases or are preserved to great depths and ultimately reach the core-mantle boundary remains controversial. In this study, experiments in the laser-heated diamond anvil cell (LHDAC) on layered samples of dolomite (Mg,Ca)CO3 and iron at pressure and temperature conditions reaching those of the deep lower mantle show that carbon-iron redox interactions destabilize the MgCO3 component, producing a mixture of diamond, Fe7C3, and (Mg,Fe)O. However, CaCO3 is preserved, supporting its relative stability in carbonate-rich lithologies under reducing lower mantle conditions. These results constrain the thermodynamic stability of redox-driven breakdown of carbonates and demonstrate progress towards multiphase mantle petrology in the LHDAC at conditions of the lowermost mantle.

58 GEOSCIENCES↗
DOE OSTI2020

HIMU geochemical signature originating from the transition zone

Plume volcanism may sample mantle sources deeper than mid-ocean ridge and arc volcanism. Ocean island basalts (OIBs) are commonly related to plume volcanism, and their diverse isotopic and elemental compositions can be described using a limited number of mantle endmembers. However, the origins and depths of these mantle endmembers are highly debated. Here we show that the HIMU (high μ, μ = 238 U/ 204 Pb) endmember may reside in the transition zone. Specifically, we report the geochemical signature of a high-pressure multiphase diamond inclusion, entrapped at 420–440 km depth and 1450 ±50 K, which matches exactly the geochemical patterns of the HIMU-rich OIBs. Since the HIMU component is variably sampled by almost all OIBs, our finding implies that the transition zone causes a major overprint of the geochemical features of mantle plumes. Furthermore, some mantle plumes, like those feeding Bermuda, St Helena, Tubuai and Mangaia, appear to be dominated by this source. Furthermore, our finding highlights the importance of the transition zone in highly incompatible element budget of the mantle.

58 GEOSCIENCES↗
DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗
DOE OSTI2017

Internal and near nozzle measurements of Engine Combustion Network “Spray G” gasoline direct injectors

Gasoline direct injection (GDI) sprays are complex multiphase flows. When compared to multi-hole diesel sprays, the plumes are closely spaced, and the sprays are more likely to interact. The effects of multi-jet interaction on entrainment and spray targeting can be influenced by small variations in the mass fluxes from the holes, which in turn depend on transients in the needle movement and small-scale details of the internal geometry. In this paper, we present a comprehensive overview of a multi-institutional effort to experimentally characterize the internal geometry and near-nozzle flow of the Engine Combustion Network (ECN) Spray G gasoline injector. In order to develop a complete picture of the near-nozzle flow, a standardized setup was shared between facilities. A wide range of techniques were employed, including both X-ray and visible-light diagnostics. The novel aspects of this work include both new experimental measurements, and a comparison of the results across different techniques and facilities. The breadth and depth of the data reveal phenomena which were not apparent from analysis of the individual data sets. We show that plume-to-plume variations in the mass fluxes from the holes can cause large-scale asymmetries in the entrainment field and spray structure. Both internal flow transients and small-scale geometric features can have an effect on the external flow. The sharp turning angle of the flow into the holes also causes an inward vectoring of the plumes relative to the hole drill angle, which increases with time due to entrainment of gas into a low-pressure region between the plumes. In conclusion, these factors increase the likelihood of spray collapse with longer injection durations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2017

Oil/water/rock wettability: Influencing factors and implications for low salinity water flooding in carbonate reservoirs

Wettability of the oil/brine/rock system is an essential petro-physical parameter which governs subsurface multiphase flow behaviour and the distribution of fluids, thus directly affecting oil recovery. Recent studies [1–3] show that manipulation of injected brine composition can enhance oil recovery by shifting wettability from oil-wet to water-wet. However, what factor(s) control system wettability has not been completely elucidated due to incomplete understanding of the geochemical system. To isolate and identify the key factors at play we used in this paper SO 4 2 —free solutions to examine the effect of salinity (formation brine/FB, 10 times diluted formation brine/10 dFB, and 100 times diluted formation brine/100 dFB) on the contact angle of oil droplets at the surface of calcite. We then compared contact angle results with predictions of surface complexation by low salinity water using PHREEQC software. We demonstrate that the conventional dilution approach likely triggers an oil-wet system at low pH, which may explain why the low salinity water EOR-effect is not always observed by injecting low salinity water in carbonated reservoirs. pH plays a fundamental role in the surface chemistry of oil/brine interfaces, and wettability. Our contact angle results show that formation brine triggered a strong water-wet system (35°) at pH 2.55, yet 100 times diluted formation brine led to a strongly oil-wet system (contact angle = 175°) at pH 5.68. Surface complexation modelling correctly predicted the wettability trend with salinity; the bond product sum ([>CaOH 2 + ][–COO - ] + [>CO 3 - ][–NH + ] + [>CO 3 - ][–COOCa + ]) increased with decreasing salinity. Finally, at pH < 6 dilution likely makes the calcite surface oil-wet, particularly for crude oils with high base number. Yet, dilution probably causes water wetness at pH > 7 for crude oils with high acid number.

02 PETROLEUM↗
DOE OSTI2019

Highly resolved Eulerian simulations of fuel spray transients in single and multi-hole injectors: Nozzle flow and near-exit dynamics

In high pressure fuel injectors, needle opening and closing transients cause complex off-design fluid dynamics behaviors that profoundly impact the spray and mixture formation processes. These dynamics are completely different from what is known to occur in steady state conditions. In this study, diesel spray transients have been investigated in single-hole and 3-hole nozzles, encompassing internal and external nozzle flow and including needle motion, performing highly resolved (2.5 μm) computational fluid dynamics (CFD) simulations. We focused on end-of-injection (EOI) and start-of-injection (SOI) processes, in order to provide insights in to the physics. The liquid fuel, vapor and gas species are modeled with a single-fluid multiphase mixture approach, with diffuse interface, and with large eddy simulations (LES) of the turbulence. Occurrence of phase change due to cavitation is accounted for, and the spray dispersion is described with a turbulent dispersion model. Detailed needle motion data and orifice internal surface are available from xray synchrotron source measurements carried out at Argonne National Laboratory, and shared through the Engine Combustion Network (ECN) community. Simulations are compared against x-ray phase contrast imaging and radiography of the internal and near-exit flow, in addition to optical microscopy data of the near-exit sprays. Simulation results are found to agree well with available experimental data, and are able to realistically capture local and global features. The simulations allow to gain insight into the physics of gas ingestion and dribbles at EOI, for different hole diameters, operating conditions and number of holes. At SOI, timing of liquid appearance out of the injector and spray tip penetration are adequately predicted, by using the EOI flow field as in-nozzle initialization, and by prescribing the measured tip needle displacement with an informed effective valve opening point inferred from the x-ray observations. Lastly, the variation of spreading angle over time is also discussed in detail for the multi-hole case, including hole-to-hole variations. Due to real geometry features and asymmetric needle motion with eccentric components, it is found that the three holes exhibit swirling flows of increasing intensity as the lift decreases, causing the near cone angle to open and spread, in a quasi-hollow cone structure. These features are not observed in axial single-hole injectors because of their relative simplicity and intrinsic symmetry.

42 ENGINEERING↗
DOE OSTI2019

Model calibration and optimization of a post-combustion CO 2 WAG pilot in a mature oil field

A carbon dioxide (CO 2 ) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO 2 flood in the West Ranch Field as part of the Petra Nova project, the world's largest postcombustion CO 2 capture and utilization initiative. With a fluvial formation geology, this is a challenging and novel application of CO 2 enhanced oil recovery. The aim of this study is two-fold. First, to build a predictive dynamic model that incorporates the data acquired during the pilot operation. Second, to utilize this model to provide guidance for optimal design for a field-wide CO 2 EOR operation. Our study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. Next, local permeability field is fine-tuned using a streamline-based method to match CO 2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO 2 recovery as a fraction of the injected CO 2 . Sensitivity studies were conducted using the calibrated pilot model to evaluate the effects of various operational parameters. Finally, we performed a Design of Experiment based multi-objective optimization of the CO 2 WAG processes using the influential operational parameters identified from the sensitivity analysis.

02 PETROLEUM↗
DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗
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 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↗