利用多级蒙特卡洛法量化亨利问题中的不确定性

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":null,"pages":null},"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\":null,\"pages\":null},\"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}
引用次数: 0

摘要

我们研究了著名的多级蒙特卡洛(MLMC)方法在密度驱动流动问题中的适用性,特别是沿海含水层盐碱化问题。作为测试案例,我们求解了不确定的亨利盐水入侵问题。未知的孔隙度、渗透率和补给参数通过随机场来模拟。经典的确定性亨利问题是非线性和随时间变化的,很容易耗费几个小时的计算时间。在不确定的情况下,需要求解确定性问题的多次变现,总计算成本会急剧增加。标准方法,如蒙特卡罗方法或基于代理的方法,是一个不错的选择,但它们在同一网格上计算所有随机变现,通常网格非常细。它们也无法平衡随机和离散误差。这些事实促使我们应用 MLMC 方法。我们证明,通过在多级空间和时间网格上求解亨利问题,MLMC 方法降低了总体计算和存储成本。为了进一步降低计算成本,我们在物理空间和随机空间进行了并行化处理。为了解决每个确定性场景,我们以黑盒方式运行并行多网格求解器 ug4。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
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.
求助全文
通过发布文献求助,成功后即可免费获取论文全文。 去求助
来源期刊
自引率
0.00%
发文量
0
期刊最新文献
Censored Data Forecasting: Applying Tobit Exponential Smoothing with Time Aggregation How to survive the Squid Games using probability theory Cross-sectional personal network analysis of adult smoking in rural areas Modeling information spread across networks with communities using a multitype branching process framework Asymptotic confidence intervals for the difference and the ratio of the weighted kappa coefficients of two diagnostic tests subject to a paired design
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1