×

一种利用自适应性降低降阶模型装配成本的超简化方法。 (英语) Zbl 1506.65179号

摘要:在非线性或含时系统的一些降阶建模方案的在线阶段的每个迭代或时间步,必须组装大型线性系统,然后将其投影到小维降阶基础上。投影的小型线性系统求解成本低,但装配和投影成为主要的计算成本。在本文中,我们介绍了一种新的超还原策略,称为简化装配(RA),它可以大幅降低这些成本。RA包含一个三角剖分自适应算法,该算法使用局部误差指示器构造一个特别适合于降阶基的简化装配三角剖分。重要的是,这种减少的装配三角剖分比原始的单元更少,从而降低装配和投影成本。我们证明了RA对Galerkin-POD型降阶模型(RAPOD)的有效性。对于用半隐式时间步法求解的非线性反应扩散问题,我们的性能比基线Galerkin-POD方法提高了五倍,对于用连续Newton-Raphson算法求解的三维超弹性问题,性能提高了七倍。这些示例在DOLFIN有限元求解器中使用线性代数的PETSc和SLEPc实现。作为补充材料,提供了生成本文结果的完整代码和数据文件。

MSC公司:

65M99型 偏微分方程、初值和含时初边值问题的数值方法
65平方米 含偏微分方程初值和初边值问题离散方程的数值解
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Morin,P.等人。;诺切托,R.H。;Siebert,K.G.,自适应有限元方法的收敛性,SIAM Rev.,44,4,631-658(2002)·Zbl 1016.65074号
[2] Niroomandi,S。;González,D。;阿尔法罗,I。;Bordeu,F。;Leygue,A。;库托,E。;Chinesta,F.,《生物软组织的实时模拟:PGD方法》,《国际数值杂志》。方法生物识别。工程师,29,5586-600(2013)
[3] Constantine,P.G.,《主动子空间:参数研究中降维的新思路》(2015),SIAM·Zbl 1431.65001号
[4] Maday,Y。;Rönquist,E.,《一种简化基础单元法》,《科学杂志》。计算。,17, 1, 447-459 (2002) ·Zbl 1014.65119号
[5] Prud'Homme,C。;罗瓦斯,D。;Veroy,L。;机械,L。;Maday,Y。;Patera,A。;Turinici,G.,参数化偏微分方程的可靠实时解:降基输出界方法,J.Fluids Eng.,124,1,70-80(2002)
[6] 君士坦丁,P。;Wang,Q.,参数化非线性动力学系统的残差最小化模型插值,SIAM J.Sci。计算。,34、4(2012),A21118-A2144·Zbl 1259.65110号
[7] Schenone,E。;韦斯,S。;Prud'Homme,C.,简化基方法的高性能计算。自然对流应用,ESAIM:Proc。,2013年12月43日,255-273·Zbl 1329.76189号
[8] Pearson,K.,《关于最接近点系的直线和平面》,Phil.Mag.,2,6,559-572(1901)
[9] Hotelling,H.,《将复杂统计变量分析为主成分》,J.Educ。心理医生。,23, 417-441 (1933)
[10] Lumley,J.L.,非均匀湍流的结构,大气。湍流。无线电波传播。,166-178 (1967)
[11] Sirovich,L.,湍流和相干结构动力学。第一部分:相干结构,夸特。申请。数学。,45,561-571(1987),网址https://www.jstor.org/stable/i40145482 ·Zbl 0676.76047号
[12] Sirovich,L。;Sirovich,C.H.,复杂现象的低维描述,Contemp。数学。,99, 277-305 (1989) ·Zbl 0681.76061号
[13] Boulakia,M。;Schenone,E。;Gerbeau,J.,心脏电生理学的降阶建模。参数识别的应用,国际期刊数字。方法生物识别。工程,28,6-7,727-744(2012)
[14] Radermacher,A。;Reese,S.,基于POD的模型简化和应用于非线性弹性的经验插值,国际。J.数字。方法工程(2015)
[15] Chinesta,F。;阿马尔。;Leygue,A。;Keunings,R.,《适当广义分解及其在计算流变学中的应用概述》,J.Non-Newton。流体力学。,166, 11, 578-592 (2011) ·兹比尔1359.76219
[16] Barrault,M。;Maday,Y。;Nguyen,N。;Patera,A.,《一种“经验插值”方法:应用于偏微分方程的有效降基离散化》,C.R.Math。阿卡德。科学。巴黎,339,667-672(2004)·Zbl 1061.65118号
[17] 格雷普,M。;Maday,Y。;Nguyen,N。;Patera,A.,《非仿射非线性偏微分方程的高效降基处理》,ESAIM Math。模型。数字。分析。,41, 3, 575-605 (2007) ·Zbl 1142.65078号
[18] Chaturantabut,S。;Sorensen,D.,通过离散经验插值进行非线性模型简化,SIAM J.Sci。计算。,32, 5, 2737-2764 (2010) ·Zbl 1217.65169号
[19] Chaturantabut,S。;Sorensen,D.,POD-DEIM非线性模型约简的状态空间误差估计,SIAM J.Numer。分析。,50, 1, 46-63 (2012) ·Zbl 1237.93035号
[20] Astrid,P。;Weiland,S。;Willcox,K。;Backx,T.,通过适当的正交分解描述的模型中的缺失点估计,IEEE Trans。自动化。控制,53,10,2237-2251(2008)·Zbl 1367.93110号
[21] Hernández,J.A。;Caicedo,医学硕士。;Ferrer,A.,《通过经验容积法对非线性有限元模型进行尺寸超缩》,计算。方法应用。机械。工程,313687-722(2017)·兹比尔1439.74419
[22] Ryckelynck,D.,涉及内部变量的力学模型的超简化,国际。J.数字。方法工程,77,1,75-89(2009)·Zbl 1195.74299号
[23] 埃弗森,R。;Sirovich,L.,Karhunen-Loève gappy data程序,J.Opt。Soc.Amer公司。A、 1657年8月12日(1995年)
[24] 卡尔伯格,K。;Bou-Mosleh,C。;Farhat,C.,《通过最小二乘Petrov-Galerkin投影和压缩张量近似实现高效非线性模型简化》,国际。J.数字。方法工程,86,2,155-181(2011)·Zbl 1235.74351号
[25] Taddei,T.,参数化泛函双范数计算的离线/在线程序:经验求积和经验测试空间,高级计算。数学。,2429-2462年5月45日(2019年)·Zbl 1435.65210号
[26] 亚诺,M。;Patera,A.T.,参数化非线性偏微分方程降基处理的LP经验求积程序,计算。方法应用。机械。工程,344,1104-1123(2019)·Zbl 1440.65238号
[27] Choi,Y。;库姆斯,D。;Anderson,R.,SNS:用于时间相关模型降阶的基于解的非线性子空间方法,SIAM J.Sci。计算。,42、2、A1116-A1146(2020)·Zbl 1442.37111号
[28] Reyes,R。;Codina,R.,流量问题降阶模型中的元素边界项:区域分解和自适应粗网格超缩,计算。方法应用。机械。工程师,368,第113159条,第(2020)页·Zbl 1506.65162号
[29] Reyes,R。;Codina,R.,基于投影的流动问题降阶模型:变分多尺度方法,计算。方法应用。机械。工程,363,第112844条pp.(2020)·Zbl 1436.76034号
[30] Fritzen,F。;Haasdonk,B。;Ryckelynck,D。;Schöps,S.,非线性热问题的超还原和离散经验插值方法的算法比较,数学。计算。申请。,23, 1, 8 (2018) ·Zbl 1390.80017号
[31] Tadmor,E。;飞利浦股份有限公司。;Ortiz,M.,固体变形的混合原子论和连续模型,Langmuir,12,4529-4534(1996)
[32] 克纳普,J。;Ortiz,M.,《准连续介质方法分析》,J.Mech。物理学。固体,49,1899-1923(2001)·Zbl 1002.74008号
[33] Gunzburger,M。;Zhang,Y.,拟包含方法的求积型近似,多尺度模型。模拟。,8, 571-590 (2010) ·Zbl 1188.70046号
[34] 杨琼。;Biyikli,E。;To,A.,多分辨率分子力学:静力学,计算。方法应用。机械。工程,258,26-38(2013)·Zbl 1296.74006号
[35] 蜜蜂。;Peerlings,R。;Geers,M.,《准连续介质方法中的中心求和》,J.Mech。物理学。固体,70242-261(2014)·Zbl 1328.74004号
[36] 洛杉矶Beex。;科尔弗里登,P。;Rabczuk,T。;Bordas,S.P.A.,《承受平面内和平面外变形的板状梁格基于准连续性的多尺度方法》,计算。方法应用。机械。工程,279,附录C,348-378(2014)
[37] Amelang,J。;文丘里尼,G。;Kochmann,D.,基于完全非局部能量的准连续介质方法的求和规则,J.Mech。物理学。固体,82378-413(2015)
[38] 阿尔内斯,M。;布莱希塔,J。;Hake,J。;Johansson,A。;Kehlet,B。;Logg,A。;Richardson,C。;Ring,J。;罗杰斯,M.E。;Wells,G.N.,《围栏项目1.5版》,Arch。数字。软质。,3, 100 (2015)
[39] 阿尔泽塔,G。;阿恩特,D。;班杰斯,W。;博德杜,V。;品牌,B。;达维多夫,D。;Gassmoeller,R。;Heister,T。;赫尔泰,L。;科尔曼,K。;Kronbichler,M。;梅尔,M。;佩尔特,J.-P。;Turcksin,B。;Wells,D.,The library,版本9.0,J.Numer。数学。,26, 4, 173-183 (2018) ·Zbl 1410.65363号
[40] Burstede,C。;Wilcox,L.C.公司。;Ghattas,O.,P4est:八叉树森林上并行自适应网格细化的可缩放算法,SIAM J.Sci。计算。,33, 3, 1103-1133 (2011) ·Zbl 1230.65106号
[41] Logg,A。;Wells,G.N.,DOLFIN:自动化有限元计算,ACM Trans。数学。软件,37,2,20:1-20:28(2010)·兹比尔1364.65254
[42] B.、S。;Abhyankar,S。;M.F.亚当斯。;Brown,J。;布鲁纳,P。;Buschelman,K。;达尔星。;埃伊霍特,V。;格罗普,W.D。;考希克,D。;Knepley,M.G。;McInnes,L.C.公司。;鲁普,K。;B.F.史密斯。;扎皮尼,S。;张,H。;Zhang,H.,PETSc用户手册技术。代表ANL-95/11-3.7版(2016年),阿贡国家实验室,URLhttp://www.mcs.anl.gov/petsc
[43] 埃尔南德斯,V。;罗曼,J.E。;Vidal,V.,SLEPc:用于解决特征值问题的可扩展且灵活的工具包,ACM Trans。数学。软质。,31, 3, 351-362 (2005) ·Zbl 1136.65315号
[44] Gonnet,P.,自适应正交误差估计综述,ACM计算。调查。,44, 4, 22:1-22:36 (2012) ·Zbl 1293.65037号
[45] Dörfler,W.,泊松方程的收敛自适应算法,SIAM J.Numer。分析。,331106-1124(2006年)·Zbl 0854.65090号
[46] Hale,J.S。;李,L。;C.N.理查森。;Wells,G.N.,便携式、生产性和性能科学计算容器,计算。科学。工程师,19,6,40-50(2017)
[47] Wriggers,P.,非线性有限元方法(2008),Springer Science&Business Media·Zbl 1153.74001号
[48] Alns,医学硕士。;Logg,A。;Ølgard,K.B。;罗杰斯,M.E。;Wells,G.N.,《统一形式语言:偏微分方程弱公式的领域特定语言》,ACM Trans。数学。软件,40,2,9:1-9:37(2014)·Zbl 1308.65175号
[49] Gräßle,C.公司。;Hinze,M.,利用任意有限元离散化对演化方程进行POD降阶建模,高级计算。数学。,44, 6, 1941-1978 (2018) ·Zbl 1404.65264号
[50] Hale,J.S。;Schenone,E。;巴罗利,D。;洛杉矶Beex。;Bordas,S.P.A.,简化装配本征正交分解方法的实施(2018)
[51] 瓦雷特,S。;Bouvry,P。;Cartiaux,H。;Georgatos,F.,《学术HPC集群的管理:UL经验》,(2014年国际高性能计算与仿真大会(HPCS 2014)(2014),IEEE:IEEE博洛尼亚,意大利),959-967
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。