Multiscale Smoothed Particle Hydrodynamics Model Development for Simulating Preferential Flow Dynamics in Fractured Porous Media

Shigorina, Elena; Rüdiger, Florian; Tartakovsky, Alexandre M.; Sauter, Martin; Kordilla, Jannes

doi:10.1029/2020wr027323

Cited by 9 publications

(10 citation statements)

References 49 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In this scheme, the solid matrix is represented as an elasto-plastic material, while the pore fluid is modeled as an incompressible fluid. SPH simulations of fluid flow through fractured porous media were also recently reported by Shigorina et al [50] and by Bui and Nguyen [51].…”

Section: Introductionsupporting

confidence: 58%

Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

Alvarado-Rodríguez

Díaz–Damacillo

Plaza

et al. 2023

Water

View full text Add to dashboard Cite

A popular equation that is often employed to represent the relationship between the pressure loss and the fluid flow in fluidized or packed granular beds is the Ergun equation, which is an extension of Darcy’s law. In this paper, the method of Smoothed Particle Hydrodynamics (SPH) is used to numerically study the flow field across a rectangular channel partially filled with a porous layer both at the Representative Elementary Volume (REV) scale using the Ergun equation and at the pore scale. Since the flow field can be estimated at the REV scale with a much lower cost compared to the pore scale, it is important to evaluate how accurately the pore-scale results can be reproduced at the REV scale. The comparison between both scales is made in terms of the velocity profiles at the outlet of the rectangular channel and the pressure losses across the clear and porous zones for three different arrays of solid grains at the pore scale. The results show that minimum differences in the flow structure and velocity profiles between the REV and the pore scale always occur at intermediate values of the porosity (ϕ=0.44 and 0.55). As the porosity increases, the differences between the REV and the pore scale also increase. The details of the pressure losses are affected by the geometry of the porous medium. In particular, we find that the pressure profiles at the REV scale match those at the pore scale almost independently of the porosity only when the grains are uniformly distributed in a non-staggered square array.

show abstract

Section: Introductionsupporting

confidence: 58%

Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

Alvarado-Rodríguez

Díaz–Damacillo

Plaza

et al. 2023

Water

View full text Add to dashboard Cite

show abstract

“…In more detailed pore scale process descriptions, the elasto-plastic deformations of the soil solid phase due to fluid pressures and interfacial tensions can be accounted for by the DEM (Yuan & Chareyere, 2022) and the smoothed particle hydrodynamics (SPH) as used in fractured rock hydrology (Shigorina et al, 2021). The DEM-based soil modeling was applied to analyze dynamics effects of tillage tools on soil loosening (Tamás et al, 2013).…”

Section: Dynamic Coupling Of Mechanic Hydraulic and Biotic Processesmentioning

confidence: 99%

3–4D soil model as challenge for future soil research: Quantitative soil modeling based on the solid phase

Gerke

Vogel

Weber

et al. 2022

J. Plant Nutr. Soil Sci.

View full text Add to dashboard Cite

A 3-4D soil model represents a logical step forward from one-dimensional soil columns (1D), two-dimensional soil maps (2D), and three-dimensional soil volumes (3D) toward dynamic soil models (4D), with time as the fourth dimension. The challenge is to develop modeling tools that account for the states of soil properties, including the spatial structure of solids and pores, as well as their dynamics, including soil mass and solute transfers in landscapes. Our envisioned 3-4D soil model approach aims at improving the capability to predict fundamental soil functions (e.g., plant growth, storage, matter fluxes) that provide ecosystem services in the socioeconomic context.This study provides a structured overview on current soil models, challenges, open questions, and urgent research needs for developing a 3-4D soil model. A 3-4D soil model should provide an inventory of spatially distributed and temporally variable soil properties. As basis for this, we propose a mass balance model for the solid phase, which needs to be supplemented by a model describing its structure. This should eventually provide adequate 3D parameter sets for the numerical modeling of soil functions (e.g., flow and transport). The target resolution is decameters in the horizontal plane and centimeters to decimeters in the vertical direction to represent characteristic soil properties and soil horizons. The actual state of soils and their properties can be estimated from spatial data that represent the soil forming factors, with the use of machine learning tools. Improved modeling of the dynamics of soil bulk density, biological processes, and the pore structure are required to relate the solid mass balance to matter fluxes.A 3-4D soil model can be built from several types of modeling approaches. We distinguish between (1) process models that simulate mass balances, fluxes and soil structure dynamics, (2) statistical pedometric models using machine learning and geostatistics to estimate the soil inventory within landscapes, and (3) pedotransfer functions to link observable attributes to specific model parameters required to simulate soil functions including water and matter fluxes. This should provide the prerequisites to predict the spatial distribution of soil functions and their changes in response to external forcing. This endeavor can draw upon many already established models and techniques, yet combining them into a newly created 3-4D soil model is a truly an ambitious, but promisingThis is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

show abstract

“…Compared to current experiment techniques, the advantage of the numerical simulation method is that the phase interface can be precisely tracked, and the 3D dynamic evolution of velocity, phase, and concentration field distribution can be well observed. At pore‐scale, the approaches for modeling multiphase flow mainly include the pore network model (An et al., 2022; Blunt, 2001; Y. F. Yang et al., 2019), the lattice Boltzmann method (LBM) (An et al., 2021; Xie et al., 2023; Zhang et al., 2019), smoothed particle hydrodynamics method (Shigorina et al., 2021), and computational fluid dynamics (CFD) method, which solves the Navier‐Stokes equation coupling with interface tracking methods, such as the volume of fluid (VOF) method (Y. F. Yang et al., 2021, 2022), the phase field method (Zhang et al., 2021), and the level set method (Jones & Detwiler, 2019). Among them, LBM and VOF methods have been used extensively to simulate multiphase flow and reactive transport in complex porous media.…”

Section: Introductionmentioning

confidence: 99%

Pore‐Scale Modeling of Coupled CO₂ Flow and Dissolution in 3D Porous Media for Geological Carbon Storage

Yang,

Wang,

Wang

et al. 2023

Water Resources Research

View full text Add to dashboard Cite

Dissolution trapping is one of the crucial trapping mechanisms for geological carbon storage in deep saline aquifers. The injected supercritical CO2 (scCO2) flow and dissolution processes are coupled and interact with each other. Therefore, we performed direct numerical simulations in three‐dimensional micro‐CT images of sandstones using the Volume of Fluid and Continuous Species Transfer (VOF‐CST) method. We investigated the coupled scCO2 flow and dissolution processes at pore scale under different rock structures, capillary numbers and rock wettability conditions. The dynamic evolution of the scCO2/brine phase distribution and scCO2 concentration distribution occurring during the injection period were presented and analyzed. Complicated coupling mechanisms between scCO2‐brine two‐phase flow and interphase mass transfer are also revealed. Our results showed that the scCO2 dissolution was highly dependent on the local distribution of scCO2 clusters. The rock with relatively high porosity and permeability would have more capacity for scCO2 injection resulting in a faster and greater dissolution of scCO2 in brine. The effect of capillary number on the scCO2 dissolution process was related to the range of capillary number. The effective upscaled (macro‐scale) mass transfer coefficient (kA) during scCO2 dissolution was evaluated and the power‐law relationship between kA and Péclet number was obtained. Rock wettability was found to be another factor controlling the scCO2 dissolution processes by affecting the scCO2‐brine interfacial area. Our pore‐scale study provides a deep understanding of the scCO2 dissolution trapping mechanisms, which is important to enhance the prediction of sequestration risk and improve sequestration efficiency.

show abstract

Multiscale Smoothed Particle Hydrodynamics Model Development for Simulating Preferential Flow Dynamics in Fractured Porous Media

Cited by 9 publications

References 49 publications

Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

3–4D soil model as challenge for future soil research: Quantitative soil modeling based on the solid phase

Pore‐Scale Modeling of Coupled CO₂ Flow and Dissolution in 3D Porous Media for Geological Carbon Storage

Contact Info

Product

Resources

About

Multiscale Smoothed Particle Hydrodynamics Model Development for Simulating Preferential Flow Dynamics in Fractured Porous Media

Cited by 9 publications

References 49 publications

Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

Smoothed Particle Hydrodynamics Simulations of Porous Medium Flow Using Ergun’s Fixed-Bed Equation

3–4D soil model as challenge for future soil research: Quantitative soil modeling based on the solid phase

Pore‐Scale Modeling of Coupled CO2 Flow and Dissolution in 3D Porous Media for Geological Carbon Storage

Contact Info

Product

Resources

About

Pore‐Scale Modeling of Coupled CO₂ Flow and Dissolution in 3D Porous Media for Geological Carbon Storage