×

刚性方程组的成对显式Runge-Kutta格式。 (英语) Zbl 1452.65128号

摘要:在本文中,我们介绍了一类适用于求解刚性方程组的显式Runge-Kutta方法,称为配对显式Runge-Kutta-(P-ERK)格式。P-ERK方法允许在域的刚性部分使用具有大量导数求值和较大绝对稳定性区域的Runge-Kutta格式,而在非刚性部分使用导数求值相对较少的格式,以降低计算成本。重要的是,可以根据局部刚度要求选择具有不同导数评估次数的不同P-ERK方案,并彼此无缝配对。然后,我们使用带有任意方案组合的欧拉方程验证了P-ERK方案获得了其设计精度。然后,我们证明P-ERK方案可以在使用Navier-Stokes方程(包括SD7003翼型上的层流和湍流)进行模拟时实现约5的加速因子。这些结果表明,当使用显式方法时,P-ERK格式可以显著加速刚性方程组的求解,并且与传统的Runge-Kutta方法和可用的参考数据相比,它们保持了准确性。

MSC公司:

65升06 常微分方程的多步、Runge-Kutta和外推方法
65升04 刚性方程的数值方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] Stetter,H.J.,常微分方程离散化方法分析(1973),Springer·Zbl 0276.65001号
[2] 海尔,E。;诺塞特,S.P。;Wanner,G.,《求解常微分方程I:非刚性问题》,计算数学中的Springer级数(1993),Springer-Verlag:Springer-Verlag Berlin-Heidelberg·Zbl 0789.65048号
[3] 海尔,E。;Wanner,H.,《求解常微分方程II:刚性和微分代数问题》(1996),Springer·Zbl 0859.65067号
[4] 阿舍尔,U.M。;Ruuth,S.J。;Spiteri,R.J.,含时偏微分方程的隐式显式Runge-Kutta方法,应用数值数学,25,2-3,151-167(1997年11月)·Zbl 0896.65061号
[5] Shoeybi,M。;Svärd先生。;哈姆,F.E。;Moin,P.,非结构化网格上可压缩流DNS和LES的自适应隐式显式格式,计算物理杂志,229,17,5944-5965(2010年8月)·Zbl 1425.76108号
[6] Hundsdorfer,W。;Jaffré,J.,空间不连续有限元的隐式显式时间步进,应用数值数学,45,2-3,231-254(2003)·Zbl 1029.65115号
[7] Kanevsky,A。;Carpenter,M.H。;Gottlieb,D。;Hesthaven,J.S.,隐式显式高阶Runge-Kutta方法在间断Galerkin格式中的应用,计算物理杂志,225,2,1753-1781(2007年8月)·Zbl 1123.65097号
[8] 佛梅尔,公元前。;Nadarajah,S.,高阶非结构化方法的自适应IMEX方案,计算物理杂志,280,261-286(2015)·Zbl 1349.76545号
[9] 阿兰帕利,V。;Hixon,R。;纳拉萨米,M。;Sawyer,S.D.,计算气动声学的高精度大步长显式Runge-Kutta(HALE-RK)格式,计算物理杂志,228,10,3837-3850(2009年6月)·兹比尔1171.76039
[10] 图洛奇,T。;Desmet,W.,应用于波传播问题的间断Galerkin空间离散的最优Runge-Kutta格式,计算物理杂志,231,42067-2091(2012年2月)·Zbl 1242.65190号
[11] 迪尔,R。;Busch,K。;Niegemann,J.,《Maxwell方程非连续Galerkin时域模拟低存储Runge-Kutta格式的比较》,《计算与理论纳米科学杂志》,7,8,1572-1580(2010年8月)
[12] 尼格曼,J。;迪尔,R。;Busch,K.,具有优化稳定区域的高效低存储Runge-Kutta格式,计算物理杂志,231,2364-372(2012年1月)·Zbl 1243.65113号
[13] 贝尔纳迪尼,M。;Pirozzoli,S.,波传播现象的Runge-Kutta格式优化的一般策略,计算物理杂志,228,11,4182-4199(2009年6月)·Zbl 1273.76294号
[14] 帕萨尼,M。;Ketcheson,D。;Deconick,W.,《用于波传播问题的谱差分法的优化显式Runge-Kutta格式》,SIAM科学计算杂志,35,2,A957-A986(2013年1月)·Zbl 1266.65157号
[15] Kubatko,E.J。;Yeager,B.A。;Ketcheson,D.I.,非连续Galerkin方法的最优保强稳定性Runge-Kutta时间离散,科学计算杂志,60,2,313-344(2014年8月)·Zbl 1304.65219号
[16] van Der Houwen,P.J。;Sommeijer,B.P.,《关于大m值显式m阶段Runge-Kutta方法的内部稳定性》,Zeitschrift für Angewandte Mathematik und Mechanik,60,10,479-485(1980)·Zbl 0455.65052号
[17] Medovikov,A.A.,抛物方程的高阶显式方法,BIT数值数学,38,2,372(1998年6月)·Zbl 0909.65060号
[18] 阿卜杜勒。;Medovikov,A.A.,基于正交多项式的二阶Chebyshev方法,数值数学,90,1,1-18(2001年11月)·Zbl 0997.65094号
[19] Abdulle,A.,具有递推关系的四阶Chebyshev方法,SIAM J.Sci。计算。,2041-2054年6月23日(2001年6月)·Zbl 1009.65048号
[20] 阿卜杜勒。;Vilmart,G.,PIROCK:瑞士刀分割隐式-显式正交Runge-Kutta-Chebyshev积分器,用于有或无噪声的刚性扩散-平流-反应问题,计算物理杂志,242869-888(2013)·Zbl 1297.65110号
[21] Martín-Vaquero,J。;Janssen,B.,刚性问题的二阶稳定显式Runge-Kutta方法,计算机物理通信,180,10,1802-1810(2009)·兹比尔1197.65006
[22] O’Sullivan,S.,一类高阶Runge-Kutta-Chebyshev稳定性多项式,计算。物理。,300、665-678(2015年11月)·Zbl 1349.65225号
[23] O'Sullivan,S.,因式分解Runge-Kutta-Chebyshev方法,《物理杂志:会议系列》,837,第012020页(2017年5月)
[24] S.O'Sullivan,《对流扩散问题的Runge-Kutta-Gegenbauer方法》,2017年。;S.O'Sullivan,《对流扩散问题的Runge-Kutta-Gegenbauer方法》,2017年。
[25] Verwer,J.G。;Sommeijer,B.P。;Hundsdorfer,W.,《对流-扩散-反应问题的RKC时间步长》,J.Compute。物理。,201、1、61-79(2004年11月)·Zbl 1059.65085号
[26] Ketcheson,D。;Ahmadia,A.,初值问题数值积分的最优稳定性多项式,应用数学和计算科学通讯,7,2,247-271(2012)·Zbl 1259.65114号
[27] Constantinescu,E.M。;Sandu,A.,双曲守恒律的多速率时间步长方法,科学计算杂志,33,3,239-278(2007)·Zbl 1127.76033号
[28] 塞尼,B。;兰布雷奇斯,J。;康布伦,R。;Legat,V。;Remacle,J.F.,用于加速显式不连续Galerkin计算并应用于地球物理流的多速率时间步进,流体数值方法国际期刊,71,1,41-64(2012)·Zbl 1431.65178号
[29] 施莱格尔,M。;诺思,O。;阿诺德,M。;Wolke,R.,对流方程的多速率Runge-Kutta格式,计算与应用数学杂志,226,2,345-357(2009)·Zbl 1167.65428号
[30] Huynh,H.T.,高阶格式的通量重建方法,包括间断Galerkin方法,(第18届AIAA计算流体动力学会议,第AIAA 2007-4079卷(2007),美国航空航天研究所)
[31] Huynh,H.T.,高阶方案的重建方法,包括用于扩散的间断Galerkin,(第47届AIAA航空航天科学会议,包括新视野论坛和航空航天博览会(2009年1月),美国航空航天研究所)
[32] Matlab 2016a,The MathWorks,马萨诸塞州纳蒂克,美国,2016年。;Matlab 2016a,MathWorks,马萨诸塞州纳蒂克,美国。
[33] M.Grant,S.Boyd CVX,用于严格凸编程的Matlab软件,2.1版。http://cvxr.com/cvx; M.Grant,S.Boyd CVX,用于严格凸编程的Matlab软件,2.1版。http://cvxr.com/cvx
[34] 格兰特,M。;Boyd,S.,非光滑凸程序的图形实现,(Blondel,V.;Boyd·兹比尔1205.90223
[35] 里德,W.H。;Hill,T.R.,《中子输运方程的三角网格法》(1973),洛斯阿拉莫斯科学实验室:洛斯阿拉莫斯科学实验室NM,美国,技术报告LA-UR-73-479,CONF-730414-2
[36] Cockburn,B。;Shu,C.,TVB Runge-Kutta局部投影非连续Galerkin守恒有限元方法,II:一般框架,计算数学,52,186,411-435(1989)·Zbl 0662.65083号
[37] 刘伟。;Vinokur,M。;王振杰,非结构网格的谱差分方法,I:基本公式,计算物理杂志,216,2,780-801(2006)·Zbl 1097.65089号
[38] 王振杰。;刘,Y。;May,G。;Jameson,A.,非结构化网格的谱差分方法,II:欧拉方程的扩展,科学计算杂志,32,1,45-71(2007)·Zbl 1151.76543号
[39] 文森特,体育。;卡斯通圭,P。;Jameson,A.,von Neumann高阶通量重建方案分析的见解,计算物理杂志,230,22,8134-8154(2011年9月)·Zbl 1343.65117号
[40] 佛梅尔,公元前。;Vincent,P.E.,关于隐式大涡模拟的能量稳定通量重建方案的性质,计算物理杂志,327368-388(2016)·Zbl 1373.76066号
[41] 佛梅尔,公元前。;Vincent,P.E.,《关于全离散通量重建方案的行为》,《应用力学与工程中的计算机方法》,3151053-1079(2017)·Zbl 1439.65121号
[42] 文森特,体育。;卡斯通圭,P。;Jameson,A.,一类新的高阶能量稳定通量重建方案,《科学计算杂志》,47,1,50-72(2011年4月)·Zbl 1433.76094号
[43] 王振杰。;Fidkowski,K。;Abgrall,R。;Caraeni,D。;Bassi,F。;A.卡里。;Deconick,H。;哈特曼,R。;Hillewaert,K。;Huynh,H.T。;北卡罗尔。;梅,G。;佩尔森,P。;van Leer,B。;Visbal,M.,《高阶CFD方法:现状与展望》,《国际流体数值方法杂志》,72,8,811-845(2013)·Zbl 1455.76007号
[44] F.D.威瑟顿。;Farrington,A.M。;PyFR,P.E.Vincent,《使用通量重建方法解决流结构上的平流-扩散型问题的开源框架》,《计算机物理通信》,185,11,3028-3040(2014)·Zbl 1348.65005号
[45] Uranga,A。;佩尔松,P。;Drela,M。;Peraire,J.,使用间断Galerkin方法对低雷诺数下湍流过渡的隐式大涡模拟,《国际工程数值方法杂志》,87,1-5,232-261(2011)·Zbl 1242.76085号
[46] 洛佩兹,M.R。;Sheshadri,A。;布尔,J.R。;Economon,T.D。;罗梅罗,J。;Watkins,J.E。;威廉姆斯,D.M。;帕拉西奥斯,F。;A.詹姆逊。;Manosalvas,D.E.,《HiFiLES的验证与验证:多GPU平台上的高阶LES非结构化解算器》,(第32届AIAA应用空气动力学会议,编号AIAA 2014-3168(2014),美国航空航天研究所)
[47] Bassi,F。;Rebay,S.,可压缩Navier-Stokes方程数值解的高精度间断有限元方法,计算物理杂志,131,2,267-279(1997年3月)·Zbl 0871.76040号
[48] 佛梅尔,公元前。;F.D.威瑟顿。;Vincent,P.E.,《关于GPU加速高阶方法在非恒定流模拟中的应用:与行业标准工具的比较》,《计算物理杂志》,334497-521(2017)
[49] 佛蒙特哥伦比亚省。;Nadarajah,S。;Tucker,P.G.,通过重构方案使用高阶校正程序进行隐式大涡模拟,《流体数值方法国际期刊》,82,5,231-260(2015)
[50] 贝克,A.D。;Bolemann,T。;弗拉德·D。;弗兰克·H。;Gassner,G.J。;辛登朗,F。;Munz,C.D.,用于过渡和湍流模拟的高阶不连续伽辽金谱元方法,国际流体数值方法杂志,76,8,522-548(2014)
[51] Garmann,D.J。;维斯巴尔,M.R。;Orkwis,P.D.,低雷诺数翼型应用隐式和亚脊模型大涡模拟技术的比较研究,国际流体数值方法杂志,71,12,1546-1565(2013)·Zbl 1430.76304号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。