×

高维参数偏微分方程解的自适应非侵入重建。 (英语) Zbl 1530.65014号

摘要:随机参数偏微分方程的数值方法可以从自适应精化方案中受益匪浅,特别是当函数近似像随机Galerkin和随机配置方法那样计算时。本文研究了基于残差估计的自适应Galerkin有限元方法的非侵入性推广。它结合了随机最小二乘法的非侵入性和随机Galerkin方法的后验误差分析。该方法使用变分蒙特卡罗方法,以高效的层次张量格式获得Galerkin投影的准最优低阶近似。我们推导了一种由可靠的误差估计器控制的自适应细化算法。与随机Galerkin方法相反,该方法很容易应用于广泛的问题,能够完全自动调整所有离散化参数。仿射和(无界)对数正态系数域的基准示例说明了非侵入自适应算法的性能,显示了单层策略的预期收敛速度。

MSC公司:

65立方米 随机微分和积分方程的数值解
65号35 偏微分方程边值问题的谱、配置及相关方法
15A69号 多线性代数,张量演算
65奈拉 涉及偏微分方程的边值问题的误差界
65年20月 数值算法的复杂性和性能
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Alns,M.、Blechta,J.、Hake,J.,Johansson,A.、Kehlet,B.、Logg,A.、Richardson,C.、Ring,J.和Rognes,M.E.以及Wells,G.N.,FEniCS项目版本1.5,Arch。数字。软质。,3(2015),第9-23页,doi:10.11588/ANS.2015.100.20553。
[2] Bachmayr,M.、Cohen,A.、DeVore,R.和Migliorati,G.,参数椭圆偏微分方程的稀疏多项式近似。第二部分:对数正态系数,ESAIM数学。模型。数字。分析。,51(2017),第341-363页·Zbl 1366.41005号
[3] Bachmayr,M.、Schneider,R.和Uschmajew,A.,《求解高维偏微分方程的张量网络和层次张量》,Found。计算。数学。,16(2016),第1423-1472页,doi:10.1007/s10208-016-9317-9·Zbl 1357.65153号
[4] Bachmayr,M.和Voulis,I.,基于随机域多级展开的自适应随机Galerkin方法:收敛性和最优性,ESAIM:数学。模型。数字。分析。,56,(2022),第1955-1992页,doi:10.1051/m2an/2022062·Zbl 1501.35147号
[5] Beck,J.,Liu,Y.,von Schwerin,E.和Tempone,R.,面向目标的自适应有限元多级蒙特卡罗收敛率,计算。方法应用。机械。工程师,402,(2022),第115582页,doi:10.1016/j.cma.2022.115582,https://arxiv.org/abs/2206.10314。 ·Zbl 1507.65008号
[6] Bespalov,A.、Powell,C.E.和Silvester,D.,参数算子方程的能量范数后验误差估计,SIAM J.Sci。计算。,36(2014),第A339-A363页,doi:10.1137/130916849·Zbl 1294.35199号
[7] Bespalov,A.、Praetorius,D.、Rocchi,L.和Ruggeri,M.,自适应随机Galerkin FEM的收敛性,SIAM J.Numer。分析。,57(2019),第2359-2382页,doi:10.1137/18m1229560·Zbl 1425.65146号
[8] Bespalov,A.、Praetorius,D.和Ruggeri,M.,自适应多级随机Galerkin FEM的收敛性和速率优化,IMA J.Numer。分析。,42(2021),第2190-2213页,doi:10.1093/imanum/drab036·Zbl 1502.65188号
[9] Bespalov,A.、Praetorius,D.和Ruggeri,M.,自适应多层随机Galerkin FEM的两层后验误差估计,SIAM/ASA J.Uncer。数量。,9,(2021),第1184-1216页,doi:10.1137/20m1342586,https://arxiv.org/abs/2006.02255。 ·Zbl 1516.65125号
[10] Bespalov,A.和Rocchi,L.,随机数据椭圆偏微分方程的高效自适应算法,SIAM/ASA J.不确定性。数量。,6(2018),第243-272页,doi:10.1137/17m1139928·Zbl 1398.65289号
[11] Bespalov,A.和Silvester,D.,参数算子方程的高效自适应随机Galerkin方法,SIAM J.Sci。计算。,38(2016),第A2118-A2140页,doi:10.1137/15m1027048·Zbl 1416.65435号
[12] Bryant,C.M.、Prudhomme,S.和Wildey,T.,具有参数不确定性的PDE响应面近似的误差分解和自适应性,SIAM/ASA J.不确定性。数量。,3(2015),第1020-1045页,doi:10.1137/140962632·Zbl 1327.65011号
[13] Carstensen,C.、Eigel,M.、Hoppe,R.H.和Löbhard,C.,统一后验有限元误差控制综述,数值。数学。理论方法应用。,5(2012年),第509-558页·Zbl 1289.65249号
[14] Cohen,A.和DeVore,R.,高维参数偏微分方程的近似,Acta Numer。,24(2015),第1-159页·Zbl 1320.65016号
[15] Cohen,A.、DeVore,R.和Schwab,C.,一类椭圆sPDE的最佳N项Galerkin近似的收敛速度,发现。计算。数学。,10(2010),第615-646页,doi:10.1007/s10208-010-9072-2·Zbl 1206.60064号
[16] Cohen,A.、DeVore,R.和Schwab,C.,参数和随机椭圆偏微分方程的解析正则性和多项式近似,Ana。申请。,9(2011),第11-47页,doi:10.1142/s0219530511001728·Zbl 1219.35379号
[17] Cohen,A.和Migliorati,G.,最优加权最小二乘法,SMAI J.计算。数学。,3(2017),第181-203页,doi:10.5802/smai-jcm.24·Zbl 1416.62177号
[18] Crowder,A.J.,Powell,C.E.和Bespalov,A.,使用隐式后验误差估计的有效自适应多级随机Galerkin近似,SIAM J.Sci。计算。,41(2019),第A1681-A1705页,doi:10.137/18m1194420·Zbl 1420.65120号
[19] Dolgov,S.、Anaya-Izquierdo,K.、Fox,C.和Scheichl,R.,《张量列分解中多元概率分布的近似和抽样》,统计计算。,30(2019),第603-625页,doi:10.1007/s11222-019-09910-z·Zbl 1436.62192号
[20] Dolgov,S.、Khoromskij,B.N.、Litvinenko,A.和Matthies,H.G.,以张量列数据格式计算响应面,预印本,arXiv:1406.28162014,https://arxiv.org/abs/1406.2816。 ·Zbl 1329.65271号
[21] Dolgov,S.、Khoromskij,B.N.、Litvinenko,A.和Matthies,H.G.,随机系数的多项式混沌展开和张量列格式的随机偏微分方程的解,SIAM/ASA J.不确定性。数量。,3(2015),第1109-1135页,doi:10.1137/140972536·Zbl 1329.65271号
[22] Dolgov,S.和Scheichl,R.,参数偏微分方程的交替最小二乘TT-交叉混合算法,SIAM/ASA J.不确定性。数量。,7(2019),第260-291页,doi:10.1137/17m113881·Zbl 1418.65170号
[23] Eigel,M.、Farchmin,N.、Heidenreich,S.和Trunschke,P.,《用张量网络高效逼近高维指数》,国际。J.Uncer。数量。,13,(2023),第25-51页,doi:10.1615/int.j.不确定性量化.2022039164,https://arxiv.org/abs/2105.09064。
[24] Eigel,M.、Gittelson,C.、Schwab,C.和Zander,E.,一种具有准最优空间网格的收敛自适应随机Galerkin有限元方法,技术报告2014-01,应用数学研讨会,瑞士苏黎世理工大学,2014,https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2014/2014-01.pdf。 ·Zbl 1335.65006号
[25] Eigel,M.、Gittelson,C.J.、Schwab,C.和Zander,E.,自适应随机Galerkin FEM,计算。方法应用。机械。工程,270(2014),第247-269页,doi:10.1016/j.cma.2013.11.015·Zbl 1296.65157号
[26] Eigel,M.、Gruhlke,R.和Marschall,M.,集中密度的低秩张量重建及其在贝叶斯反演中的应用,统计计算。,32,(2022),doi:10.1007/s11222-022-10087-1,https://arxiv.org/abs/2008.04264。 ·Zbl 1484.62002年
[27] Eigel,M.、Marschall,M.和Multerer,M.,随机扰动域的自适应随机Galerkin张量列离散化,SIAM/ASA J.Uncer。数量。,8,(2020),第1189-1214页,doi:10.1137/19m1246080,https://arxiv.org/abs/1902.07753。 ·Zbl 1447.65137号
[28] Eigel,M.、Marschall,M.,Pfeffer,M.和Schneider,R.,分层张量表示中对数正态系数的自适应随机Galerkin FEM,Numer。数学。,145(2020),第655-692页,doi:10.1007/s00211-020-01123-1·Zbl 1439.65160号
[29] Eigel,M.、Marschall,M.和Schneider,R.,自适应分层张量表示的无采样贝叶斯反演,反演问题,34(2018),035010,doi:10.1088/1361-6420/aaa998·Zbl 1404.65261号
[30] Eigel,M.、Merdon,C.和Neumann,J.,一种自适应多级蒙特卡罗方法,具有不确定数据感兴趣量的随机边界,SIAM/ASA J.不确定性。数量。,4(2016),第1219-1245页,doi:10.1137/15m1016448·Zbl 1398.35306号
[31] Eigel,M.、Neumann,J.、Schneider,R.和Wolf,S.,高维随机偏微分方程的非侵入张量重建,计算。方法应用。数学。,19(2019),第39-53页,doi:10.1515/cmam-2018-0028·Zbl 1419.35259号
[32] Eigel,M.、Pfeffer,M.和Schneider,R.,具有层次张量表示的自适应随机伽辽金有限元法,Numer。数学。,136(2016),第765-803页,doi:10.1007/s00211-016-0850-x·Zbl 1397.65262号
[33] Eigel,M.、Schneider,R.和Trunschke,P.,经验非线性最小二乘法的收敛界,ESAIM数学。模型。数字。分析。,56(2022),第79-104页,doi:10.1051/m2an/2021070·Zbl 1482.62071号
[34] Eigel,M.、Schneider,R.、Trunschke,P.和Wolf,S.,《机器学习和高维偏微分方程的变分蒙特卡罗桥接概念》,高级计算。数学。,45(2019),第2503-2532页,doi:10.1007/s10444-019-09723-8·Zbl 1435.68258号
[35] Ernst,O.G.和Sprungk,B.,《随机数据下椭圆偏微分方程的随机配置:对数正态情况》,摘自《稀疏网格和应用》,慕尼黑2012年,施普林格,查姆,2014年,第29-53页·Zbl 1316.65008号
[36] Espig,M.、Hackbusch,W.、Litvinenko,A.、Matthies,H.G.和Wähnert,P.,张量格式随机Galerkin矩阵的有效低阶近似,计算。数学。申请。,67(2014),第818-829页·Zbl 1350.65005号
[37] Galvis,J.和Sarkis,M.,在无均匀椭圆的情况下逼近无穷维随机Darcy方程,SIAM J.Numer。分析。,47(2009),第3624-3651页,doi:10.1137/080717924·Zbl 1205.60121号
[38] Ghanem,R.G.和Spanos,P.,《随机有限元:谱方法》,Springer,纽约,1991年·Zbl 0722.73080号
[39] Gittelson,C.J.,对数正态各向同性扩散问题的随机Galerkin离散化,数学。模型方法应用。科学。,20(2010),第237-263页,doi:10.1142/s021820251004210·Zbl 1339.65216号
[40] Hackbusch,W.和Schneider,R.,张量空间和层次张量表示,in,Springer,Cham,2014,第237-261页,doi:10.1007/978-3319-08159-5_12·Zbl 1317.65102号
[41] Hoang,V.H.和Schwab,C.,对数正态高斯随机输入椭圆偏微分方程的N项Wiener混沌近似率,数学。模型方法应用。科学。,24(2014),第797-826页,doi:10.1142/s0218202513500681·Zbl 1294.65010号
[42] Holtz,S.、Rohwedder,T.和Schneider,R.,关于固定TT-rank张量的流形,Numer。数学。,120(2011),第701-731页,doi:10.1007/s00211-011-0419-7·兹比尔1242.15022
[43] Huber,B.和Wolf,S.,Xerus——通用Tensor图书馆,2014-2021年,https://libxerus.org/。
[44] Maêtre,O.Le和Knio,O.M.,《不确定性量化的光谱方法:在计算流体动力学中的应用》,施普林格,多德雷赫特,荷兰,2010年·Zbl 1193.76003号
[45] Mugler,A.,Verallgemeinentes Polynomilels Chaos zur Lösung Stationärer Diffusionsproblem mit Zufälligen Koeffizienten,博士论文,英国科技大学科特布斯-森芬伯格分校,2013年。
[46] Mugler,A.和Starkloff,H.-J.,关于随机椭圆偏微分方程随机Galerkin方法的收敛性,ESAIM Math。模型。数字。分析。,47(2013),第1237-1263页·兹比尔1297.65010
[47] Nobile,F.、Tempone,R.和Webster,C.G.,《随机输入数据偏微分方程的各向异性稀疏网格随机配置方法》,SIAM J.Numer。分析。,46(2008),第2411-2442页,doi:10.1137/070680540·Zbl 1176.65007号
[48] Nobile,F.,Tempone,R.和Webster,C.G.,《随机输入数据偏微分方程的稀疏网格随机配置方法》,SIAM J.Numer。分析。,46(2008),第2309-2345页,doi:10.1137/060663660·Zbl 1176.65137号
[49] Nouy,A.,《第4章:高维近似和模型降阶的低秩方法》,载于《模型降阶和近似》,SIAM,费城,2017年,第171-226页,doi:10.1137/1.9781611974829.ch4。
[50] Oseledets,I.和Tyrtyshnikov,E.,多维数组的TT-交叉近似,线性代数应用。,432(2010),第70-88页,doi:10.1016/j.laa.209.07.24·兹比尔1183.65040
[51] Oseledets,I.V.,张量-应变分解,SIAM J.Sci。计算。,33(2011),第2295-2317页,doi:10.1137/090752286·Zbl 1232.15018号
[52] Oseledets,I.V.和Tyrtyshnikov,E.E.,《打破维度诅咒,或如何在多个维度中使用SVD》,SIAM J.Sci。计算。,31(2009),第3744-3759页,doi:10.1137/090748330·Zbl 1200.65028号
[53] Prudhomme,S.和Bryant,C.M.,响应面近似的自适应代理建模及贝叶斯推理应用,高级模型。模拟。工程科学。,2(2015),22,doi:10.1186/s40323-015-0045-5。
[54] Schwab,C.和Gittelson,C.J.,高维参数和随机PDE的稀疏张量离散化,Acta Numer。,20(2011),第291-467页,doi:10.1017/S0962492911000055·Zbl 1269.65010号
[55] 史蒂文森,R.,由二等分创建的局部精化单形分区的完成,数学。公司。,77(2008),第227-241页,doi:10.1090/s0025-5718-07-01959-x·Zbl 1131.65095号
[56] Trunschke,P.,非线性最小二乘的收敛界限及其在张量恢复中的应用,预印本,arXiv:2108.052372021,https://arxiv.org/abs/1208.05237。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。