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
-
This paper presents a second-order, implicit numerical model for one-dimensional, large strain thaw consolidation of ice-rich, fine-grained permafrost. The phase composition of permafrost at sub-freezing temperatures is determined using an unfrozen water content model that accounts for both capillary and adsorptive unfrozen water. The model incorporates secondary compression to improve the accuracy of long-term thaw consolidation simulations. The algorithm incorporates conduction, advection, and phase change in heat transfer and simultaneous occurrence of primary consolidation and secondary compression. Benchmarking and verification of the model show good agreement with existing numerical models. The proposed model is validated against experimental observations. The model indicates that adsorbed unfrozen water dominates over a wide range of sub-freezing temperatures, while capillary unfrozen water freezes at temperatures just below the freezing point. Numerical simulations suggest that ignoring secondary compression can lead to underestimation of excess pore pressure and settlement during both thaw and post-thaw consolidation. Void ratio and average degree of consolidation are overestimated when secondary compression is not considered. The effect of secondary compression on excess pore pressure and void ratio during thawing becomes more pronounced in thicker, field-scale permafrost layers. Results from this study highlight the importance of considering adsorptive and capillary unfrozen water to determine permafrost composition and incorporating secondary compression in thaw consolidation modeling and thaw settlement estimation for long-term civil infrastructure planning in cold regions. The proposed model provides a comprehensive framework for simulating thaw consolidation processes in permafrost regions.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
