English

Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method

Computational Engineering, Finance, and Science 2024-03-27 v1 Numerical Analysis Mathematical Physics math.MP Numerical Analysis Other Statistics

Abstract

We investigate the applicability of the well-known multilevel Monte Carlo (MLMC) method to the class of density-driven flow problems, in particular the problem of salinisation of coastal aquifers. As a test case, we solve the uncertain Henry saltwater intrusion problem. Unknown porosity, permeability and recharge parameters are modelled by using random fields. The classical deterministic Henry problem is non-linear and time-dependent, and can easily take several hours of computing time. Uncertain settings require the solution of multiple realisations of the deterministic problem, and the total computational cost increases drastically. Instead of computing of hundreds random realisations, typically the mean value and the variance are computed. The standard methods such as the Monte Carlo or surrogate-based methods is a good choice, but they compute all stochastic realisations on the same, often, very fine mesh. They also do not balance the stochastic and discretisation errors. These facts motivated us to apply the MLMC method. We demonstrate that by solving the Henry problem on multi-level spatial and temporal meshes, the MLMC method reduces the overall computational and storage costs. To reduce the computing cost further, parallelization is performed in both physical and stochastic spaces. To solve each deterministic scenario, we run the parallel multigrid solver ug4 in a black-box fashion.

Cite

@article{arxiv.2403.17018,
  title  = {Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method},
  author = {Dmitry Logashenko and Alexander Litvinenko and Raul Tempone and Ekaterina Vasilyeva and Gabriel Wittum},
  journal= {arXiv preprint arXiv:2403.17018},
  year   = {2024}
}

Comments

20 pages, 11 figures, 3 tables. arXiv admin note: substantial text overlap with arXiv:2302.07804

R2 v1 2026-06-28T15:33:06.875Z