Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: Computational analysis

Eric Morales-Casique, Shlomo P Neuman, Alberto Guadagnini

Research output: Contribution to journalArticle

32 Citations (Scopus)

Abstract

In a companion paper we presented exact (though implicit and not closed) nonlocal conditional first and second moment equations for nonreactive advective-dispersive transport under both steady state and transient flow regimes in bounded, randomly heterogeneous porous domains. To allow solving our nonlocal equations we developed recursive moment equations in Laplace space for the special case of steady state flow to second order in σY (a measure of the standard deviation of natural log hydraulic conductivity Y, which is generally nonhomogeneous) and proposed a higher-order iterative closure scheme. We also showed that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that generally varies in space-time. The purpose of this paper is to explore the behavior and assess the accuracy and computational efficiency of our moment solutions in comparison to conditional and unconditional Monte Carlo simulations. To do so, we present a high-accuracy computational algorithm for our iterative nonlocal and recursive localized moment equations and corresponding computational results in two spatial dimensions conditional on measurements of Y. Our algorithm solves the moment equations by finite elements in Laplace-transformed space and inverts the solution numerically back into the time domain. Conditional results obtained with our iterative algorithm compare well with Monte Carlo simulations for σY2 = 0.3 and Peclet number Pe = 100 defined in terms of the integral scale of Y, and for Pe = 10 in the unconditional case. As σY2, Pe and time increase the quality of our iterative moment solution deteriorates. We show that this is due to our disregarding velocity moments of order higher than two and propose that including such moments should render our iterative solution workable over a wider range of these parameters. Second-order recursive nonlocal and space-localized results are considerably less accurate than those obtained with our iterative nonlocal algorithm. Even though our moment solution does not require computing (space-time localized) macrodispersion coefficients, we nevertheless do so to examine the influence of boundaries and conditioning on their behavior. Our results support an earlier observation by the second author [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: theoretical framework. Adv Water Resour, in press., Neuman SP. Universal scaling of hydraulic conductivities and dispersivities in geologic media. Water Resour Res 1990;26(8):1749-58], based on world-wide tracer test results, that the rate at which apparent longitudinal dispersivity increases with scale diminishes with conditioning. In preliminary runs conducted on a relatively small grid without optimizing our algorithms and without parallelization, the moment solutions required considerably less computer time than did the Monte Carlo simulations.

Original languageEnglish (US)
Pages (from-to)1399-1418
Number of pages20
JournalAdvances in Water Resources
Volume29
Issue number9
DOIs
StatePublished - Sep 2006

Fingerprint

solute transport
porous medium
dispersivity
conditioning
hydraulic conductivity
advection
simulation
transient flow
analysis
tracer
water

Keywords

  • Advection
  • Dispersion
  • Moments
  • Perturbation solution
  • Random media
  • Stochastic
  • Transport

ASJC Scopus subject areas

  • Earth-Surface Processes

Cite this

Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media : Computational analysis. / Morales-Casique, Eric; Neuman, Shlomo P; Guadagnini, Alberto.

In: Advances in Water Resources, Vol. 29, No. 9, 09.2006, p. 1399-1418.

Research output: Contribution to journalArticle

@article{e03a9e1e63f84b3992ebeffd61864e86,
title = "Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: Computational analysis",
abstract = "In a companion paper we presented exact (though implicit and not closed) nonlocal conditional first and second moment equations for nonreactive advective-dispersive transport under both steady state and transient flow regimes in bounded, randomly heterogeneous porous domains. To allow solving our nonlocal equations we developed recursive moment equations in Laplace space for the special case of steady state flow to second order in σY (a measure of the standard deviation of natural log hydraulic conductivity Y, which is generally nonhomogeneous) and proposed a higher-order iterative closure scheme. We also showed that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that generally varies in space-time. The purpose of this paper is to explore the behavior and assess the accuracy and computational efficiency of our moment solutions in comparison to conditional and unconditional Monte Carlo simulations. To do so, we present a high-accuracy computational algorithm for our iterative nonlocal and recursive localized moment equations and corresponding computational results in two spatial dimensions conditional on measurements of Y. Our algorithm solves the moment equations by finite elements in Laplace-transformed space and inverts the solution numerically back into the time domain. Conditional results obtained with our iterative algorithm compare well with Monte Carlo simulations for σY2 = 0.3 and Peclet number Pe = 100 defined in terms of the integral scale of Y, and for Pe = 10 in the unconditional case. As σY2, Pe and time increase the quality of our iterative moment solution deteriorates. We show that this is due to our disregarding velocity moments of order higher than two and propose that including such moments should render our iterative solution workable over a wider range of these parameters. Second-order recursive nonlocal and space-localized results are considerably less accurate than those obtained with our iterative nonlocal algorithm. Even though our moment solution does not require computing (space-time localized) macrodispersion coefficients, we nevertheless do so to examine the influence of boundaries and conditioning on their behavior. Our results support an earlier observation by the second author [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: theoretical framework. Adv Water Resour, in press., Neuman SP. Universal scaling of hydraulic conductivities and dispersivities in geologic media. Water Resour Res 1990;26(8):1749-58], based on world-wide tracer test results, that the rate at which apparent longitudinal dispersivity increases with scale diminishes with conditioning. In preliminary runs conducted on a relatively small grid without optimizing our algorithms and without parallelization, the moment solutions required considerably less computer time than did the Monte Carlo simulations.",
keywords = "Advection, Dispersion, Moments, Perturbation solution, Random media, Stochastic, Transport",
author = "Eric Morales-Casique and Neuman, {Shlomo P} and Alberto Guadagnini",
year = "2006",
month = "9",
doi = "10.1016/j.advwatres.2005.10.014",
language = "English (US)",
volume = "29",
pages = "1399--1418",
journal = "Advances in Water Resources",
issn = "0309-1708",
publisher = "Elsevier Limited",
number = "9",

}

TY - JOUR

T1 - Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media

T2 - Computational analysis

AU - Morales-Casique, Eric

AU - Neuman, Shlomo P

AU - Guadagnini, Alberto

PY - 2006/9

Y1 - 2006/9

N2 - In a companion paper we presented exact (though implicit and not closed) nonlocal conditional first and second moment equations for nonreactive advective-dispersive transport under both steady state and transient flow regimes in bounded, randomly heterogeneous porous domains. To allow solving our nonlocal equations we developed recursive moment equations in Laplace space for the special case of steady state flow to second order in σY (a measure of the standard deviation of natural log hydraulic conductivity Y, which is generally nonhomogeneous) and proposed a higher-order iterative closure scheme. We also showed that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that generally varies in space-time. The purpose of this paper is to explore the behavior and assess the accuracy and computational efficiency of our moment solutions in comparison to conditional and unconditional Monte Carlo simulations. To do so, we present a high-accuracy computational algorithm for our iterative nonlocal and recursive localized moment equations and corresponding computational results in two spatial dimensions conditional on measurements of Y. Our algorithm solves the moment equations by finite elements in Laplace-transformed space and inverts the solution numerically back into the time domain. Conditional results obtained with our iterative algorithm compare well with Monte Carlo simulations for σY2 = 0.3 and Peclet number Pe = 100 defined in terms of the integral scale of Y, and for Pe = 10 in the unconditional case. As σY2, Pe and time increase the quality of our iterative moment solution deteriorates. We show that this is due to our disregarding velocity moments of order higher than two and propose that including such moments should render our iterative solution workable over a wider range of these parameters. Second-order recursive nonlocal and space-localized results are considerably less accurate than those obtained with our iterative nonlocal algorithm. Even though our moment solution does not require computing (space-time localized) macrodispersion coefficients, we nevertheless do so to examine the influence of boundaries and conditioning on their behavior. Our results support an earlier observation by the second author [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: theoretical framework. Adv Water Resour, in press., Neuman SP. Universal scaling of hydraulic conductivities and dispersivities in geologic media. Water Resour Res 1990;26(8):1749-58], based on world-wide tracer test results, that the rate at which apparent longitudinal dispersivity increases with scale diminishes with conditioning. In preliminary runs conducted on a relatively small grid without optimizing our algorithms and without parallelization, the moment solutions required considerably less computer time than did the Monte Carlo simulations.

AB - In a companion paper we presented exact (though implicit and not closed) nonlocal conditional first and second moment equations for nonreactive advective-dispersive transport under both steady state and transient flow regimes in bounded, randomly heterogeneous porous domains. To allow solving our nonlocal equations we developed recursive moment equations in Laplace space for the special case of steady state flow to second order in σY (a measure of the standard deviation of natural log hydraulic conductivity Y, which is generally nonhomogeneous) and proposed a higher-order iterative closure scheme. We also showed that, under a limited set of conditions, the mean transport equation can be localized to yield a familiar-looking advection-dispersion equation with a conditional macrodispersion tensor that generally varies in space-time. The purpose of this paper is to explore the behavior and assess the accuracy and computational efficiency of our moment solutions in comparison to conditional and unconditional Monte Carlo simulations. To do so, we present a high-accuracy computational algorithm for our iterative nonlocal and recursive localized moment equations and corresponding computational results in two spatial dimensions conditional on measurements of Y. Our algorithm solves the moment equations by finite elements in Laplace-transformed space and inverts the solution numerically back into the time domain. Conditional results obtained with our iterative algorithm compare well with Monte Carlo simulations for σY2 = 0.3 and Peclet number Pe = 100 defined in terms of the integral scale of Y, and for Pe = 10 in the unconditional case. As σY2, Pe and time increase the quality of our iterative moment solution deteriorates. We show that this is due to our disregarding velocity moments of order higher than two and propose that including such moments should render our iterative solution workable over a wider range of these parameters. Second-order recursive nonlocal and space-localized results are considerably less accurate than those obtained with our iterative nonlocal algorithm. Even though our moment solution does not require computing (space-time localized) macrodispersion coefficients, we nevertheless do so to examine the influence of boundaries and conditioning on their behavior. Our results support an earlier observation by the second author [Morales Casique E, Neuman SP, Guadagnini A. Nonlocal and localized analyses of nonreactive solute transport in bounded randomly heterogeneous porous media: theoretical framework. Adv Water Resour, in press., Neuman SP. Universal scaling of hydraulic conductivities and dispersivities in geologic media. Water Resour Res 1990;26(8):1749-58], based on world-wide tracer test results, that the rate at which apparent longitudinal dispersivity increases with scale diminishes with conditioning. In preliminary runs conducted on a relatively small grid without optimizing our algorithms and without parallelization, the moment solutions required considerably less computer time than did the Monte Carlo simulations.

KW - Advection

KW - Dispersion

KW - Moments

KW - Perturbation solution

KW - Random media

KW - Stochastic

KW - Transport

UR - http://www.scopus.com/inward/record.url?scp=33746612848&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33746612848&partnerID=8YFLogxK

U2 - 10.1016/j.advwatres.2005.10.014

DO - 10.1016/j.advwatres.2005.10.014

M3 - Article

AN - SCOPUS:33746612848

VL - 29

SP - 1399

EP - 1418

JO - Advances in Water Resources

JF - Advances in Water Resources

SN - 0309-1708

IS - 9

ER -