Manuscript Published in UNCECOMP 2023, 5th International Conference on Uncertainty Quantification in Computational Science and Engineering
MULTILEVEL AND MULTIGRID METHODS FOR SOLVING HENRY PROBLEM WITH UNCERTAIN COEFFICIENTS: A joint work of Alexander Litvinenko, Dmitry Logashenko, Raul Tempone, Ekaterina Vasilyeva and Gabriel Wittum has been published in UNCECOMP 2023 proceedings which were organized from 12-14 June. The conference focused on analyzing and designing processes with uncertainty quantification, emphasizing multiscale analysis and design of complex systems.
We are solving a problem of salinization of coastal aquifers. As a test case example, we consider the Henry saltwater intrusion problem. Since porosity, permeability and recharge are unknown or only known at a few points, we model them using random fields. The Henry problem describes a two-phase flow and is nonlinear and time-dependent. The solution to be found is the expectation of the salt mass fraction, which is uncertain and time-dependent. To estimate this expectation we use the well known multilevel Monte Carlo (MLMC) method. The MLMC method takes just a few samples on computationaly expensive (fine) meshes and more samples on cheap (coarse) meshes. Then, by building a telescoping sum, the MLMC method estimates the expected value at a much lower cost than the classical Monte Carlo method. The deterministic solver used here is the well-known parallel and scalable ug4 solver.