An official website of the United States government
This study aims to bridge length scales in immiscible multiphase flow simulation by connecting two published governing equations at the pore-scale and continuum-scale through a novel validation framework. We employ Niessner and Hassnaizadeh's [“A model for two-phase flow in porous media including fluid-fluid interfacial area,” Water Resour. Res. 44(8), W08439 (2008)] continuum-scale model for multiphase flow in porous media, combined with the geometric equation of state of McClure et al. [“Modeling geometric state for fluids in porous media: Evolution of the Euler characteristic,” Transp. Porous Med. 133(2), 229–250 (2020)]. Pore-scale fluid configurations simulated with the lattice-Boltzmann method are used to validate the continuum-scale results. We propose a mapping from the continuum-scale to pore-scale utilizing a generalized additive model to predict non-wetting phase Euler characteristics during imbibition, effectively bridging the continuum-to-pore length scale gap. Continuum-scale simulated measures of specific interfacial area, saturation, and capillary pressure are directly compared to up-scaled pore-scale simulation results. This research develops a numerical framework capable of capturing multiscale flow equations establishing a connection between pore-scale and continuum-scale simulations.
more »
« less
Fully Implicit Dynamic Pore‐Network Modeling of Two‐Phase Flow and Phase Change in Porous Media
Abstract Dynamic pore‐network model (PNM) has been widely used to model pore‐scale two‐phase flow. Numerical algorithms commonly used for dynamic PNM including IMPES (implicit pressure explicit saturation) and IMP‐SIMS (implicit pressure semi‐implicit saturation) can be numerically unstable or inaccurate for challenging flow regimes such as low capillary number (Ca) flow and unfavorable displacements. We perform comprehensive analyses of IMPES and IMP‐SIMS for a wide range of flow regimes under drainage conditions and develop a novel fully implicit (FI) algorithm to address their limitations. Our simulations show the following: (1) While IMPES was reported to be numerically unstable for lowCaflow, using a smoothed local pore‐body capillary pressure curve appears to produce stable simulations. (2) Due to an approximation for the capillary driving force, IMP‐SIMS can deviate from quasi‐static solutions at equilibrium states especially in heterogeneous networks. (3) Both IMPES and IMP‐SIMS introduce mass conservation errors. The errors are small for networks with cubic pore bodies (less than 1.4% for IMPES and 1.2% for IMP‐SIMS). They become much greater for networks with square‐tube pore bodies (up to 45% for IMPES and 46% for IMP‐SIMS). Conversely, the new FI algorithm is numerically stable and mass conservative regardless of the flow regimes and pore geometries. It also precisely recovers the quasi‐static solutions at equilibrium states. The FI framework has been extended to include compressible two‐phase flow, multicomponent transport, and phase change dynamics. Example simulations of two‐phase displacements accounting for phase change show that evaporation and condensation can suppress fingering patterns generated during invasion.
more »
« less
- Award ID(s):
- 2023351
- PAR ID:
- 10452824
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 56
- Issue:
- 11
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this work, we review and describe our computational framework for solving multiphysics phase-field fracture problems in porous media. Therein, the following five coupled nonlinear physical models are addressed: displacements (geo-mechanics), a phase-field variable to indicate the fracture position, a pressure equation (to describe flow), a proppant concentration equation, and/or a saturation equation for two-phase fracture flow, and finally a finite element crack width problem. The overall coupled problem is solved with a staggered solution approach, known in subsurface modeling as the fixed-stress iteration. A main focus is on physics-based discretizations. Galerkin finite elements are employed for the displacement-phase-field system and the crack width problem. Enriched Galerkin formulations are used for the pressure equation. Further enrichments using entropy-vanishing viscosity are employed for the proppant and/or saturation equations. A robust and efficient quasi-monolithic semi-smooth Newton solver, local mesh adaptivity, and parallel implementations allow for competitive timings in terms of the computational cost. Our framework can treat two- and three-dimensional realistic field and laboratory examples. The resulting program is an in-house code named IPACS (Integrated Phase-field Advanced Crack Propagation Simulator) and is based on the finite element library deal.II. Representative numerical examples are included in this document.more » « less
-
Abstract Reconstruction of the magnetic field, electric current, and plasma pressure is provided using a new data mining (DM) method with weighted nearest neighbors (NN) for strong storms with the storm activity indexSym‐H < −300 nT, the Bastille Day event (July 2000), and the 20 November 2003 superstorm. It is shown that the new method significantly reduces the statistical bias of the original NN algorithm toward weaker storms. In the DM approach the magnetic field is reconstructed using a small NN subset of the large historical database, with the subset numberKNN ≫ 1being still much larger than any simultaneous multiprobe observation number. This allows one to fit with observations a very flexible magnetic field model using basis function expansions for equatorial and field‐aligned currents, and at the same time, to keep the model sensitive to storm variability. This also allows one to calculate the plasma pressure by integrating the quasi‐static force balance equation with the isotropic plasma approximation. For strong storms of particular importance becomes the resolution of the eastward current, which prevents the divergence of the pressure integral. It is shown that in spite of the strong reduction of the dominant NN number in the new weighted NN algorithm to capture strong storm features, it is still possible to resolve the eastward current and to retrieve plasma pressure distributions. It is found that the pressure peak for strong storms may be as close as≈2.1REto Earth and its value may exceed 300 nPa.more » « less
-
Abstract Convergence of a finite element method with mass-lumping and flux upwinding is formulated for solving the immiscible two-phase flow problem in porous media. The method approximates directly the wetting phase pressure and saturation, which are the primary unknowns. Well-posedness is obtained in [J. Numer. Math., 29(2), 2021]. Theoretical convergence is proved via a compactness argument. The numerical phase saturation converges strongly to a weak solution in L 2 in space and in time whereas the numerical phase pressures converge strongly to weak solutions in L 2 in space almost everywhere in time. The proof is not straightforward because of the degeneracy of the phase mobilities and the unboundedness of the derivative of the capillary pressure.more » « less
-
Abstract Inorganic aragonite occurs in a wide spectrum of depositional environments and its precipitation is controlled by complex physio‐chemical factors. This study investigates diagenetic conditions that led to aragonite cement precipitation in Cenozoic glaciomarine deposits of McMurdo Sound, Antarctica. A total of 42 sandstones that host intergranular cement were collected from theCIROS‐1 core, located proximal to the terminus of Ferrar Glacier. Standard petrography, Raman spectroscopy and electron microprobe analysis reveal a prominent aragonite cement phase that occurs as a pore‐filling blocky fabric throughout the core. Oxygen isotope compositions (δ18O = −30·0 to −8·6‰ Vienna Pee‐Dee Belemnite) and clumped isotope temperatures (TΔ47 = 13·1 to 31·5°C) determined from the aragonite cements provide precise constraints on isotopic compositions (δ18Ow) of the parent fluid, which mostly range from −10·8 to −7·2‰ Vienna Standard Mean Ocean Water. The fluidδ18Owvalues are consistent with those of pore water, previously identified as cryogenic brine in the nearbyAND‐2A core. Petrographic and geochemical data suggest that aragonite cement in theCIROS‐1 core precipitated from a similar brine. The brine likely formed and infiltrated sediments in flooded glacial valleys along the western margin of McMurdo Sound during the middle Miocene Climatic Transition, and subsequently flowed basinward in the subsurface. Consequently, the brine forms as a longstanding subsurface fluid that has saturated Cenozoic sediments below southern McMurdo Sound since at least the middle Miocene. Aragonite cementation in theCIROS‐1 core is interpreted to reflect its proximal position to sites of brine formation and greater likelihood of experiencing brines with sustained high carbonate saturation states and Mg/Ca ratios. This unusual occurrence expands the range of known natural occurrences of aragonite cement. Given the potential for cryogenic brine formation in glaciomarine settings, blocky aragonite, as the end member of the spectrum of aragonite cement morphology, may be more widespread in glaciomarine sediments than currently thought.more » « less
