### Abstract

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 language | English (US) |
---|---|

Pages (from-to) | 131-161 |

Number of pages | 31 |

Journal | Communications in Computational Physics |

Volume | 6 |

Issue number | 1 |

Publication status | Published - Jul 2009 |

### Fingerprint

### Keywords

- 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

*Communications in Computational Physics*,

*6*(1), 131-161.