Effect of Pore-Wall Roughness and Péclet Number on Conservative Solute Transport in Saturated Porous Media

While modeling solute transport has been an active subject of research in the past few decades, the influence of pore-wall roughness on contaminant migration has not yet been addressed. We therefore conduct particle tracking simulations in three porous domains that have different pore-wall roughness characteristics. Specifically, we consider five surface fractal dimensions ds = 1.0, 1.1, 1.2, 1.4, and 1.6, and four different Péclet numbers Pe = 10, 102, 103, and 105. Overall, arrival time distributions are simulated for 60 scenarios (3 domains x 5 surface fractal dimensions x 4 Péclet numbers) some of which show heavy-tailed patterns indicating non-Fickian transport. To interpret the simulations and quantify the transport behavior, we analyze the resulting arrival time distributions by the continuous time random walk (CTRW) approach. Results show that, on average, as the surface fractal dimension increases from 1.0 to 1.6, the CTRW model parameters 𝛽, an exponent showing the degree of anomalous transport, v, the average solute velocity, and t2, the cut-off time to Fickian transport, remain nearly constant. However, the dispersion coefficient, D, increases and the characteristic transition time, t1, decreases. We found t1 and D are more sensitive to pore-wall roughness compared to the other CTRW parameters. We also found that as the Péclet number increases from 10 to 105, on average, v and D increase, t1 and 𝛽 decrease, and t2 remains nearly constant. The simulations demonstrate that the exponent 𝛽 and the dispersion coefficient are correlated to the average solute velocity.

An edited version of this paper was published by AGU. Copyright (year) American Geophysical Union [Effect of Pore‐Wall Roughness and Péclet Number on Conservative Solute Transport in Saturated Porous Media. Water Resources Research 59, 2 (2023)]

Files

Metadata

Work Title Effect of Pore-Wall Roughness and Péclet Number on Conservative Solute Transport in Saturated Porous Media
Access
Open Access
Creators
  1. Behzad Ghanbarian
  2. Yashar Mehmani
  3. Brian Berkowitz
Keyword
  1. Particle tracking
  2. Péclet number
  3. Roughness
  4. Solute transport
  5. Surface fractal dimension
License In Copyright (Rights Reserved)
Work Type Article
Publisher
  1. Water Resources Research
Publication Date January 22, 2023
Publisher Identifier (DOI)
  1. https://doi.org/10.1029/2022WR033119
Deposited April 25, 2024

Versions

Analytics

Collections

This resource is currently not in any collection.

Work History

Version 1
published

  • Created
  • Added Ghanbarian_Mehmani_Berkowitz__2023__-_roughness_solute__preprint_.pdf
  • Added Creator Behzad Ghanbarian
  • Added Creator Yashar Mehmani
  • Added Creator Brian Berkowitz
  • Published
  • Updated
  • Updated Keyword, Description, Publication Date Show Changes
    Keyword
    • Particle tracking, Péclet number, Roughness, Solute transport, Surface fractal dimension
    Description
    • <p>While modeling solute transport has been an active subject of research in the past few decades, the influence of pore-wall roughness on contaminant migration has not yet been addressed. We therefore conduct particle tracking simulations in three porous domains that have different pore-wall roughness characteristics. Specifically, we consider five surface fractal dimensions d<sub>s</sub> = 1.0, 1.1, 1.2, 1.4, and 1.6, and four different Péclet numbers Pe = 10, 10<sup>2</sup>, 10<sup>3</sup>, and 10<sup>5</sup>. Overall, arrival time distributions are simulated for 60 scenarios (3 domains (Formula presented.) 5 surface fractal dimensions (Formula presented.) 4 Péclet numbers) some of which show heavy-tailed patterns indicating non-Fickian transport. To interpret the simulations and quantify the transport behavior, we analyze the resulting arrival time distributions by the continuous time random walk (CTRW) approach. Results show that, on average, as the surface fractal dimension increases from 1.0 to 1.6, the CTRW model parameters (Formula presented.), an exponent showing the degree of anomalous transport, v, the average solute velocity, and t<sub>2</sub>, the cut-off time to Fickian transport, remain nearly constant. However, the dispersion coefficient, D, increases and the characteristic transition time, t<sub>1</sub>, decreases. We found t<sub>1</sub> and D are more sensitive to pore-wall roughness compared to the other CTRW parameters. We also found that as the Péclet number increases from 10 to 10<sup>5</sup>, on average, v and D increase, t<sub>1</sub> and (Formula presented.) decrease, and t<sub>2</sub> remains nearly constant. The simulations demonstrate that the exponent (Formula presented.) and the dispersion coefficient are correlated to the average solute velocity.</p>
    • <p>While modeling solute transport has been an active subject of research in the past few decades, the influence of pore-wall roughness on contaminant migration has not yet been addressed. We therefore conduct particle tracking simulations in three porous domains that have different pore-wall roughness characteristics. Specifically, we consider five surface fractal dimensions d<sub>s</sub> = 1.0, 1.1, 1.2, 1.4, and 1.6, and four different Péclet numbers Pe = 10, 10<sup>2</sup>, 10<sup>3</sup>, and 10<sup>5</sup>. Overall, arrival time distributions are simulated for 60 scenarios (3 domains x 5 surface fractal dimensions x 4 Péclet numbers) some of which show heavy-tailed patterns indicating non-Fickian transport. To interpret the simulations and quantify the transport behavior, we analyze the resulting arrival time distributions by the continuous time random walk (CTRW) approach. Results show that, on average, as the surface fractal dimension increases from 1.0 to 1.6, the CTRW model parameters 𝛽, an exponent showing the degree of anomalous transport, v, the average solute velocity, and t<sub>2</sub>, the cut-off time to Fickian transport, remain nearly constant. However, the dispersion coefficient, D, increases and the characteristic transition time, t<sub>1</sub>, decreases. We found t<sub>1</sub> and D are more sensitive to pore-wall roughness compared to the other CTRW parameters. We also found that as the Péclet number increases from 10 to 10<sup>5</sup>, on average, v and D increase, t<sub>1</sub> and 𝛽 decrease, and t<sub>2</sub> remains nearly constant. The simulations demonstrate that the exponent 𝛽 and the dispersion coefficient are correlated to the average solute velocity.</p>
    Publication Date
    • 2023-02-01
    • 2023-01-22