skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 10 until 2:00 AM ET on Friday, October 11 due to maintenance. We apologize for the inconvenience.


Title: A macroscopic two-length-scale model for natural convection in porous media driven by a species-concentration gradient
The modelling of natural convection in porous media is receiving increased interest due to its significance in environmental and engineering problems. State-of-the-art simulations are based on the classic macroscopic Darcy–Oberbeck–Boussinesq (DOB) equations, which are widely accepted to capture the underlying physics of convection in porous media provided the Darcy number, $Da$ , is small. In this paper we analyse and extend the recent pore-resolved direct numerical simulations (DNS) of Gasow et al. ( J. Fluid Mech , vol. 891, 2020, p. A25) and show that the macroscopic diffusion, which is neglected in DOB, is of the same order (with respect to $Da$ ) as the buoyancy force and the Darcy drag. Consequently, the macroscopic diffusion must be modelled even if the value of $Da$ is small. We propose a ‘two-length-scale diffusion’ model, in which the effect of the pore scale on the momentum transport is approximated with a macroscopic diffusion term. This term is determined by both the macroscopic length scale and the pore scale. It includes a transport coefficient that solely depends on the pore-scale geometry. Simulations of our model render a more accurate Sherwood number, root mean square (r.m.s.) of the mass concentration and r.m.s. of the velocity than simulations that employ the DOB equations. In particular, we find that the Sherwood number $Sh$ increases with decreasing porosity and with increasing Schmidt number $(Sc)$ . In addition, for high values of $Ra$ and high porosities, $Sh$ scales nonlinearly. These trends agree with the DNS, but are not captured in the DOB simulations.  more » « less
Award ID(s):
2042834
NSF-PAR ID:
10299405
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
926
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Natural convection in porous media is a fundamental process for the long-term storage of CO 2 in deep saline aquifers. Typically, details of mass transfer in porous media are inferred from the numerical solution of the volume-averaged Darcy–Oberbeck–Boussinesq (DOB) equations, even though these equations do not account for the microscopic properties of a porous medium. According to the DOB equations, natural convection in a porous medium is uniquely determined by the Rayleigh number. However, in contrast with experiments, DOB simulations yield a linear scaling of the Sherwood number with the Rayleigh number ( $Ra$ ) for high values of $Ra$ ( $Ra\gg 1300$ ). Here, we perform direct numerical simulations (DNS), fully resolving the flow field within the pores. We show that the boundary layer thickness is determined by the pore size instead of the Rayleigh number, as previously assumed. The mega- and proto-plume sizes increase with the pore size. Our DNS results exhibit a nonlinear scaling of the Sherwood number at high porosity, and for the same Rayleigh number, higher Sherwood numbers are predicted by DNS at lower porosities. It can be concluded that the scaling of the Sherwood number depends on the porosity and the pore-scale parameters, which is consistent with experimental studies. 
    more » « less
  2. The study of thermal convection in porous media is of both fundamental and practical interest. Typically, numerical studies have relied on the volume-averaged Darcy–Oberbeck–Boussinesq (DOB) equations, where convection dynamics are assumed to be controlled solely by the Rayleigh number ( Ra ). Nusselt numbers ( Nu ) from these models predict Nu – Ra scaling exponents of 0.9–0.95. However, experiments and direct numerical simulations (DNS) have suggested scaling exponents as low as 0.319. Recent findings for solutal convection between DNS and DOB models have demonstrated that the ‘pore-scale parameters’ not captured by the DOB equations greatly influence convection. Thermal convection also has the additional complication of different thermal transport properties (e.g. solid-to-fluid thermal conductivity ratio k s / k f and heat capacity ratio σ ) in different phases. Thus, in this work we compare results for thermal convection from the DNS and DOB equations. On the effects of pore size, DNS results show that Nu increases as pore size decreases. Mega-plumes are also found to be more frequent and smaller for reduced pore sizes. On the effects of conjugate heat transfer, two groups of cases (Group 1 with varying k s / k f at σ  = 1 and Group 2 with varying σ at k s / k f  = 1) are examined to compare the Nu – Ra relations at different porosity ( ϕ ) and k s / k f and σ values. Furthermore, we report that the boundary layer thickness is determined by the pore size in DNS results, while by both the Rayleigh number and the effective heat capacity ratio, $\bar{\phi } = \phi + (1 - \phi )\sigma$ , in the DOB model. 
    more » « less
  3. This work employs single-mode equations to study convection and double-diffusive convection in a porous medium where the Darcy law provides large-scale damping. We first consider thermal convection with salinity as a passive scalar. The single-mode solutions resembling steady convection rolls reproduce the qualitative behavior of root-mean-square and mean temperature profiles of time-dependent states at high Rayleigh numbers from direct numerical simulations (DNS). We also show that the single-mode solutions are consistent with the heat-exchanger model that describes well the mean temperature gradient in the interior. The Nusselt number predicted from the single-mode solutions exhibits a scaling law with Rayleigh number close to that followed by exact 2D steady convection rolls, although large aspect ratio DNS results indicate a faster increase. However, the single-mode solutions at a high wavenumber predict Nusselt numbers close to the DNS results in narrow domains. We also employ the single-mode equations to analyze the influence of active salinity, introducing a salinity contribution to the buoyancy, but with a smaller diffusivity than the temperature. The single-mode solutions are able to capture the stabilizing effect of an imposed salinity gradient and describe the standing and traveling wave behaviors observed in DNS. The Sherwood numbers obtained from single-mode solutions show a scaling law with the Lewis number that is close to the DNS computations with passive or active salinity. This work demonstrates that single-mode solutions can be successfully applied to this system whenever periodic or no-flux boundary conditions apply in the horizontal. 
    more » « less
  4. null (Ed.)
    In this paper, we consider an important problem for modeling complex coupled phenomena in porous media at multiple scales. In particular, we consider flow and transport in the void space between the pores when the pore space is altered by new solid obstructions formed by microbial growth or reactive transport, and we are mostly interested in pore-coating and pore-filling type obstructions, observed in applications to biofilm in porous media and hydrate crystal formation, respectively. We consider the impact of these obstructions on the macroscopic properties of the porous medium, such as porosity, permeability and tortuosity, for which we build an experimental probability distribution with reduced models, which involves three steps: (1) generation of independent realizations of obstructions, followed by, (2) flow and transport simulations at pore-scale, and (3) upscaling. For the first step, we consider three approaches: (1A) direct numerical simulations (DNS) of the PDE model of the actual physical process called BN which forms the obstructions, and two non-DNS methods, which we call (1B) CLPS and (1C) LP. LP is a lattice Ising-type model, and CLPS is a constrained version of an Allen–Cahn model for phase separation with a localization term. Both LP and CLPS are model approximations of BN, and they seek local minima of some nonconvex energy functional, which provide plausible realizations of the obstructed geometry and are tuned heuristically to deliver either pore-coating or pore-filling obstructions. Our methods work with rock-void geometries obtained by imaging, but bypass the need for imaging in real-time, are fairly inexpensive, and can be tailored to other applications. The reduced models LP and CLPS are less computationally expensive than DNS, and can be tuned to the desired fidelity of the probability distributions of upscaled quantities. 
    more » « less
  5. Reactive transport modeling of subsurface environments plays an important role in addressing critical problems of geochemical processes, such as dissolution and precipitation of minerals. Current transport models for porous media span various scales, ranging from pore-scale to continuum-scale. In this study, we established an upscaling method connecting pore-scale and continuum-scale models by employing a deep learning methodology of Convolutional Neural Networks (CNNs). We applied Darcy-Brinkmann-Stokes (DBS) method to simulate the fluid flow and reactive transport in pore-scale models, which would act as constituents of a continuum-scale model. The datasets of spatial pore distribution of subvolume samples were used as the input for the deep learning model, and the continuum (Darcy)-scale parameters such as permeability, effective surface area, and effective diffusion coefficient were figured out as outputs (i.e., labels). By applying the trained models of the subvolumes in the entire sample volume, we generated the initial field of porosity, permeability, effective diffusion coefficient, and effective surface area for continuum-scale simulation of a mineral dissolution problem. We took an acid dissolution case as an example to utilize the outcomes of trained deep learning models as input data in the continuum-scale simulation. This work presents a comprehensive upscaling workflow, as bridging the findings of microscale simulations to the continuum-scale simulations of a reactive transport problem. 
    more » « less