×

大气扩散模拟中的多级蒙特卡罗和改进的时间步长方法。 (英语) Zbl 1380.65012号

摘要:模拟污染物在大气中的传输和扩散的一种常见方法是使用随机拉格朗日扩散模型。在数学上,这些模型用随机微分方程(SDE)描述湍流输运过程。计算瓶颈是蒙特卡罗算法,它模拟了湍流速度场中大量模型粒子的运动;对于每个粒子,使用数值时间步长方法计算轨迹。在操作应急应用中,选择有效的数值方法尤其重要,例如跟踪核事故中的放射性云或预测火山灰云对国际航空的影响,而准确及时的预测至关重要。在本文中,我们研究了多层蒙特卡罗(MLMC)方法在模拟大气边界层中具有代表性的一维弥散场景中粒子传播的应用。与目前用于大气弥散建模的标准蒙特卡罗(StMC)方法相比,MLMC可以产生渐近更高的计算复杂度和更低的计算成本。为了在非渐近情况下降低该方法的绝对成本,同样重要的是在每个级别上选择尽可能好的数值时间步长方法。为了研究这一点,我们还将许多操作模型中使用的标准辛Euler方法与两种基于SDE分裂方法的改进时间步长算法进行了比较。

MSC公司:

65二氧化碳 蒙特卡罗方法
86A10美元 气象学和大气物理学
60 H10型 随机常微分方程(随机分析方面)
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 德雷克斯勒,R。;阿诺德·D·。;奇诺,M。;Galmarini,S。;霍特,M。;琼斯,A。;Leadbetter,S。;马洛,A。;Maurer,C。;罗尔夫,G。;齐藤,K。;Servanckx,R。;Shimbori,T。;索拉佐,E。;Wotawa,G.,世界气象组织福岛第一核电站事故放射性核素扩散和沉积的模型模拟,J.Environ。放射性。,139, 172-184 (2015)
[2] Dacre,H.F。;Grant,A.L。;霍根,R.J。;贝尔彻,S.E。;汤姆森,D。;Devenish,B。;马兰科,F。;霍特,M。;海伍德,J.M。;Ansmann,A.,《使用激光雷达观测和NAME模拟评估2010年埃亚菲亚德拉火山喷发初始阶段火山灰羽的结构和大小》,J.Geophys。大气研究。,116,D20(2011)
[3] 韦伯斯特,H.N。;汤姆森·D·J。;约翰逊,B.T。;赫德,I.P.C。;Turnbull,K。;马兰科,F。;新泽西州克里斯蒂安森。;多尔西,J。;Minikin,A。;魏泽尔,B。;舒曼,美国。;斯帕克斯,R.S.J。;Loughlin,南卡罗来纳州。;Hort,医学博士。;Leadbetter,S.J。;Devenish,B.J。;曼宁,A.J。;威瑟姆,C.S。;海伍德,J.M。;Golding,B.W.,《2010年埃亚菲亚德拉火山喷发远端火山云中火山灰浓度的操作预测》,J.Geophys。第117号决议,D20(2012年)
[4] Jones,A.,英国气象局大气扩散模型,天气,59311-316(2004)
[5] 琼斯,A。;汤姆森,D。;霍特,M。;Devenish,B.,英国气象局的下一代大气扩散模型,NAME III,(空气污染建模及其应用十七(2007)),580-589
[6] F.B.Smith,M.J.Clark,《切尔诺贝利核电站事故空中碎片的运输和沉积,特别强调对英国的影响》,气象局科学论文第42号。;F.B.Smith,M.J.Clark,《切尔诺贝利核电站事故空中碎片的运输和沉积,特别强调对英国的影响》,气象局科学论文第42号。
[7] 雷丁顿,A.L。;德温特,R.G。;威瑟姆,C.S。;Manning,A.J.,模拟硫酸盐和硝酸盐气溶胶对云、pH和氨排放的敏感性,大气。环境。,43, 20, 3227-3234 (2009)
[8] Heinrich,S.,《多级蒙特卡罗方法》,(Margenov,S.;Wa-si-niewski,J.;Yalamov,P.,《大尺度科学计算》,《计算机科学讲义》,第2179卷(2001),施普林格:施普林格-柏林-海德堡),58-67·Zbl 1031.65005号
[9] Giles,M.B.,多层蒙特卡罗路径模拟,Oper。研究,56,3,607-617(2008)·Zbl 1167.65316号
[10] Giles,M.B.,多层蒙特卡罗方法,数值学报。,24, 259 (2015) ·Zbl 1316.65010号
[11] 贾尔斯,M。;Szpruch,L.,金融应用的多级蒙特卡罗方法,arXiv预印本·Zbl 1277.91193号
[12] 格雷厄姆,I.G。;Scheichl,R。;Ullmann,E.,对数正态扩散和多级蒙特卡罗方法的混合有限元分析,Stoch。部分差异。Equ.、。,分析。计算。,4, 1, 41-75 (2016) ·Zbl 1371.65008号
[13] 穆勒,E.H。;Scheichl,R。;Shardlow,T.,《随机微分方程的改进多级蒙特卡罗及其在朗之万方程中的应用》,Proc。R.Soc.A,数学。物理学。工程科学。,471, 2176 (2015) ·Zbl 1372.65020号
[14] Skorokhod,A.V.,有界区域扩散过程的随机方程,理论概率论。申请。,264-274年6月3日(1961年)·Zbl 0215.53501号
[15] 伯纳多,M。;巴德,C。;Champneys,A.R。;Kowalczyk,P.,《分段-光滑动力系统:理论与应用》,第163卷(2008),施普林格科学与商业媒体·Zbl 1146.37003号
[16] Wilson,J.D。;Flesch,T.K.,《随机飞行色散模型中的流动边界:实施良好混合条件》,J.Appl。美托洛尔。,32, 11, 1695-1707 (1993)
[17] Rodean,H.C.,湍流扩散的随机拉格朗日模型,第26卷(1996),美国气象学会
[18] 贾尔斯,M.B。;Nagapetyan,T。;Ritter,K.,分布函数和密度的多级蒙特卡罗近似,SIAM/ASA J.不确定性。量化。,3, 1, 267-295 (2015) ·Zbl 1322.65014号
[19] Kloeden,P。;Platen,E.,随机微分方程的数值解,(随机建模和应用概率(2011),施普林格:施普林格-柏林-海德堡)
[20] 莱姆库勒,B。;Matthews,C.,《分子动力学:确定性和随机数值方法》,Interdiscip。申请。数学。(2015),施普林格国际出版公司·Zbl 1351.82001号
[21] Maruyama,G.,连续马尔可夫过程和随机方程,Rend。循环。马特·巴勒莫,4148-90(1955)·Zbl 0053.40901号
[22] 北卡罗来纳州Bou-Rabee。;Owhadi,H.,随机背景下变分积分器的长运行精度,SIAM J.Numer。分析。,48, 1, 278-297 (2010) ·Zbl 1215.65012号
[23] Hoel,H。;Von Schwerin,E。;塞佩西,A。;Tempone,R.,自适应多级蒙特卡罗算法的实现与分析,蒙特卡罗方法应用。,20, 1, 1-41 (2014) ·Zbl 1284.65011号
[24] 贾尔斯,M.B。;莱斯特,C。;Whittle,J.,《多级蒙特卡罗计算中的非测试自适应时间步长》,303-314(2016),Springer International Publishing:Springer国际出版公司Cham·Zbl 1356.65017号
[25] E.H.Mueller,G.Katsiolides,T.Shardlow,MLMCLangevin代码,2017年5月,https://doi.org/10.5281/zenodo.580662; E.H.Mueller,G.Katsiolides,T.Shardlow,MLMCLangevin代码,2017年5月,https://doi.org/10.5281/zenodo.580662
[26] Shardlow,T.,随机微分方程的修正方程,BIT Numer。数学。,46, 1, 111-125 (2006) ·Zbl 1091.65005号
[27] Zygalakis,K.C.,《关于随机微分方程修正方程的存在性和应用》,SIAM J.Sci。计算。,33, 1, 102-130 (2011) ·Zbl 1236.60069号
[28] 汤姆森,D.J.,湍流中颗粒轨迹随机模型的选择标准,流体力学杂志。,180, 529-556 (1987) ·Zbl 0644.76065号
[29] Kolmogorov,A.N.,《非常大雷诺数下不可压缩粘性流体湍流的局部结构》,(Dokl.Akad.Nauk SSSR.Dokl.Akad.Nauk SSSR,JSTOR,第30卷(1941)),301-305
[30] Wilson,J。;Yee,E。;Ek,N。;d'Amours,R.,《城市环境中风力运输的拉格朗日模拟》,Q.J.R.Meteorol。Soc.,135,643,1586-1602(2009)
[31] H.Webster,D.Thomson,N.Morrison,NAME的新湍流剖面,气象局湍流和扩散注288。;H.Webster,D.Thomson,N.Morrison,NAME的新湍流剖面,英国气象局湍流和扩散注释288。
[32] Khasminskii,R.,微分方程的随机稳定性,第66卷(2011),Springer Science&Business Media
[33] 罗德,G.J。;鲍威尔,C.E。;Shardlow,T.,《计算随机偏微分方程导论》,第50卷(2014),剑桥大学出版社:剑桥大学出版社,纽约·Zbl 1327.60011号
[34] 科利尔,N。;哈吉·阿利,A.-L。;诺比尔,F。;von Schwerin,E。;Tempone,R.,一种连续多层蒙特卡罗算法,BIT-Numer。数学。,55, 2, 399-432 (2015) ·Zbl 1317.65030号
[35] Hoel,H。;Schwerin,E。;塞佩西,A。;Tempone,R.,自适应多级蒙特卡罗模拟,(多尺度计算的数值分析:2009年班夫国际研究站冬季研讨会论文集(2012),施普林格:施普林格-柏林-海德堡),217-234·Zbl 1245.65002号
[36] 方,W。;Giles,M.B.,非全局Lipschitz漂移SDE的自适应Euler-Maruyama方法:第一部分,有限时间间隔,arXiv预打印
[37] Kuo,F.Y。;Scheichl,R。;施瓦布,C。;斯隆,I.H。;Ullmann,E.,对数正态扩散问题的多层准蒙特卡罗方法,arXiv预印本·Zbl 1368.65005号
[38] Kebaier,A。;Lelong,J.,使用样本平均近似的耦合重要性采样和多级蒙特卡罗(2015年10月13日)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。