Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , , , |
|---|---|
| Format: | Preprint |
| Publié: |
2024
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866911812843208704 |
|---|---|
| author | Logashenko, Dmitry Litvinenko, Alexander Tempone, Raul Vasilyeva, Ekaterina Wittum, Gabriel |
| author_facet | Logashenko, Dmitry Litvinenko, Alexander Tempone, Raul Vasilyeva, Ekaterina Wittum, Gabriel |
| contents | 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2403_17018 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method Logashenko, Dmitry Litvinenko, Alexander Tempone, Raul Vasilyeva, Ekaterina Wittum, Gabriel Computational Engineering, Finance, and Science Numerical Analysis Mathematical Physics Other Statistics 65M55, 65M06, 65M22, 65M50 G.1.3; G.1.8; G.1.10; I.6.5 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. |
| title | Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method |
| topic | Computational Engineering, Finance, and Science Numerical Analysis Mathematical Physics Other Statistics 65M55, 65M06, 65M22, 65M50 G.1.3; G.1.8; G.1.10; I.6.5 |
| url | https://arxiv.org/abs/2403.17018 |