Laplace-transform finite element solution of nonlocal and localized stochastic moment equations of transport

Eric Morales-Casique, Shlomo P Neuman

Research output: Contribution to journalArticle

5 Scopus citations


Morales-Casique et al. (Adv. Water Res., 29 (2006), pp. 1238-1255) developed exact first and second nonlocal moment equations for advective-dispersive transport in finite, randomly heterogeneous geologic media. The velocity and concentration in these equations are generally nonstationary due to trends in heterogeneity, conditioning on site data and the influence of forcing terms. Morales-Casique et al. (Adv. Water Res., 29 (2006), pp. 1399-1418) solved the Laplace transformed versions of these equations recursively to second order in the standard deviation σy of (natural) log hydraulic conductivity, and iteratively to higher-order, by finite elements followed by numerical inversion of the Laplace transform. They did the same for a space-localized version of the mean transport equation. Here we recount briefly their theory and algorithms; compare the numerical performance of the Laplace-transform finite element scheme with that of a high-accuracy ULTIMATE-QUICKEST algorithm coupled with an alternating split operator approach; and review some computational results due to Morales-Casique et al. (Adv. Water Res., 29 (2006), pp. 1399-1418) to shed light on the accuracy and computational efficiency of their recursive and iterative solutions in comparison to conditional Monte Carlo simulations in two spatial dimensions.

Original languageEnglish (US)
Pages (from-to)131-161
Number of pages31
JournalCommunications in Computational Physics
Issue number1
Publication statusPublished - Jul 2009



  • Conditioning
  • Finite elements
  • Geologic media
  • Laplace transform
  • Localization
  • Moment equations
  • Nonlocality
  • Porous media
  • Random heterogeneity
  • Stochastic transport
  • Ultimate-quickest

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)

Cite this