Dmitry Logashenko, Alexander Litvinenko, Raul Tempone, Ekaterina Vasilyeva, Gabriel Wittum
{"title":"利用多级蒙特卡洛法量化亨利问题中的不确定性","authors":"Dmitry Logashenko, Alexander Litvinenko, Raul Tempone, Ekaterina Vasilyeva, Gabriel Wittum","doi":"arxiv-2403.17018","DOIUrl":null,"url":null,"abstract":"We investigate the applicability of the well-known multilevel Monte Carlo\n(MLMC) method to the class of density-driven flow problems, in particular the\nproblem of salinisation of coastal aquifers. As a test case, we solve the\nuncertain Henry saltwater intrusion problem. Unknown porosity, permeability and\nrecharge parameters are modelled by using random fields. The classical\ndeterministic Henry problem is non-linear and time-dependent, and can easily\ntake several hours of computing time. Uncertain settings require the solution\nof multiple realisations of the deterministic problem, and the total\ncomputational cost increases drastically. Instead of computing of hundreds\nrandom realisations, typically the mean value and the variance are computed.\nThe standard methods such as the Monte Carlo or surrogate-based methods is a\ngood choice, but they compute all stochastic realisations on the same, often,\nvery fine mesh. They also do not balance the stochastic and discretisation\nerrors. These facts motivated us to apply the MLMC method. We demonstrate that\nby solving the Henry problem on multi-level spatial and temporal meshes, the\nMLMC method reduces the overall computational and storage costs. To reduce the\ncomputing cost further, parallelization is performed in both physical and\nstochastic spaces. To solve each deterministic scenario, we run the parallel\nmultigrid solver ug4 in a black-box fashion.","PeriodicalId":501323,"journal":{"name":"arXiv - STAT - Other Statistics","volume":"70 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-02-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method\",\"authors\":\"Dmitry Logashenko, Alexander Litvinenko, Raul Tempone, Ekaterina Vasilyeva, Gabriel Wittum\",\"doi\":\"arxiv-2403.17018\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"We investigate the applicability of the well-known multilevel Monte Carlo\\n(MLMC) method to the class of density-driven flow problems, in particular the\\nproblem of salinisation of coastal aquifers. As a test case, we solve the\\nuncertain Henry saltwater intrusion problem. Unknown porosity, permeability and\\nrecharge parameters are modelled by using random fields. The classical\\ndeterministic Henry problem is non-linear and time-dependent, and can easily\\ntake several hours of computing time. Uncertain settings require the solution\\nof multiple realisations of the deterministic problem, and the total\\ncomputational cost increases drastically. Instead of computing of hundreds\\nrandom realisations, typically the mean value and the variance are computed.\\nThe standard methods such as the Monte Carlo or surrogate-based methods is a\\ngood choice, but they compute all stochastic realisations on the same, often,\\nvery fine mesh. They also do not balance the stochastic and discretisation\\nerrors. These facts motivated us to apply the MLMC method. We demonstrate that\\nby solving the Henry problem on multi-level spatial and temporal meshes, the\\nMLMC method reduces the overall computational and storage costs. To reduce the\\ncomputing cost further, parallelization is performed in both physical and\\nstochastic spaces. To solve each deterministic scenario, we run the parallel\\nmultigrid solver ug4 in a black-box fashion.\",\"PeriodicalId\":501323,\"journal\":{\"name\":\"arXiv - STAT - Other Statistics\",\"volume\":\"70 1\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2024-02-22\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - STAT - Other Statistics\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2403.17018\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Other Statistics","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2403.17018","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Uncertainty quantification in the Henry problem using the multilevel Monte Carlo method
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.