×

一种精确、高效、稳定的三维磁流体力学非分裂交错网格格式。 (英文) Zbl 1349.76494号

摘要:本文扩展了二维磁流体力学(MHD)的非分裂交错网格格式(USM)[作者和A.E.迪恩同上,第228号,第4期,第952-975页(2009年;Zbl 1330.76093号)]全3D MHD方案。该方案是一种有限体积Godunov方法,由约束传输(CT)方法和高效准确的单步定向无分裂多维数据重建进化算法组成,扩展了Colella的原始二维角点传输迎风(CTU)方法[P.科尔拉同上,87,第1号,171-200(1990年;Zbl 0694.65041号)]. 我们提出了两种3D数据重建进化算法:(1)简化CTU方案和(2)完整CTU方案。简化的3D CTU方案是Collela的2D CTU方法的简单3D扩展的变体,被视为2D USM方案的直接扩展。完整的3D CTU方案是我们的主要3D解算器,其中包括用于稳定性的所有多维交叉导数项。后一种方法在逻辑上类似于3D未拼接CTU方法J.萨尔茨曼[同上,115,第1号,153-168(1994年;Zbl 0813.65111号)]. 我们算法的主要创新点有两个。首先,我们将简化的CTU方案推广到能够在CFL数接近于一的情况下运行的完整CTU方案。这两种方法都利用二维USM算法中开发的横向更新技术来计算横向通量没有解决中间黎曼问题,从而通过减少黎曼解的总数提供具有成本效益的3D方法。所提出的算法简单有效,尤其是当包含维持平面内磁场动力学的多维MHD项时。其次,我们介绍了一种新的CT方案,它利用适当的逆风信息来获取电场的平均值。我们的三维USM方案可以很容易地与各种重建方法(例如,一阶Godunov、二阶MUSCL-Hancock、三阶PPM和五阶WENO)以及多种基于一维的黎曼解算器(例如,局部Lax-Friedrichs、HLLE、HLLC、HLLD和Roe)相结合。3D USM-MHD解算器在芝加哥大学Flash中心的官方Flash版本中提供。

MSC公司:

76M20码 有限差分方法在流体力学问题中的应用
2006年6月65日 含偏微分方程初值和初边值问题的有限差分方法
76周05 磁流体力学和电流体力学
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Abgrall,R.,用Roe型方法逼近可压缩流体力学的多维Riemann问题,C.R.Acad。科学。巴黎Sér I数学。,319, 499 (1994) ·Zbl 0813.76074号
[2] Almgren,A.S。;贝尔,J.B。;伦德曼,C.A。;Zingale,M.,Ia型超新星的低马赫数建模,I.流体动力学,天体物理学。J.,637922-936(2006)
[3] Almgren,A.S。;贝尔,J.B。;伦德曼,C.A。;Zingale,M.,CASTRO:一个新的可压缩天体物理解算器I.流体动力学和自重,天体物理学。J.,715,1221-1238(2010)
[4] Balsara,D.S。;Spicer,D.S.,在磁流体动力学模拟中使用高阶Godunov磁通确保螺线管磁场的交错网格算法,J.Compute。物理。,149, 270-292 (1999) ·Zbl 0936.76051号
[5] Balsara,D.S。;Shu,C.-W.,精度越来越高的保单调加权非振荡格式,J.Compute。物理。,160405-452(2000年)·Zbl 0961.65078号
[6] Balsara,D.S.,《无发散重建磁流体动力学的二阶精确方案》,天体物理学。《补充期刊》,151149-184(2004)
[7] Balsara,D.S。;阿尔特曼,C。;蒙兹,C.D。;Dumbser,M.,RKDG方案中基于子小区的问题区域指标和一种新型的混合型RKDG+HWENO方案,J.Comput。物理。,226, 586-620 (2007) ·Zbl 1124.65072号
[8] Balsara,D.S。;臀部,T。;Dumbser,M。;Munz,C.D.,《流体力学和无发散磁流体力学的高效、高精度ADER-WENO格式》,J.Compute。物理。,228, 2480-2516 (2009) ·Zbl 1275.76169号
[9] Balsara,D.S.,《多维HLLE-Riemann解算器:欧拉和磁流体力学流动的应用》,J.Compute。物理。,229, 1970-1993 (2010) ·Zbl 1303.76140号
[10] Balsara,D.S.,《守恒定律的二维HLLC-Riemann解算器:应用于欧拉和磁流体力学流动》,J.Compute。物理。,231, 7476-7503 (2012) ·Zbl 1284.76261号
[11] Balsara,D.S。;梅耶,C。;Dumbser,M。;杜,H。;Xu,Z.,结构网格上Euler和磁流体动力学流动的ADER方案的有效实现——与Runge-Kutta方法的速度比较,J.Compute。物理。,235, 934-969 (2013) ·Zbl 1291.76237号
[13] 贝尔·J。;Dawson,C.N。;Shubin,G.R.,《多维标量守恒定律的一种非分裂的高阶Godunov方法》,J.Compute。物理。,74, 1-24 (1988) ·Zbl 0684.65088号
[14] 布里奥,M。;Zakharian,A.R。;Webb,G.M.,欧拉气体动力学方程的二维黎曼解算器,J.Compute。物理。,167, 177-195 (2001) ·兹比尔1043.76042
[15] Colella,P.,双曲守恒定律的多维迎风方法,J.Compute。物理。,87, 171-200 (1990) ·Zbl 0694.65041号
[16] 科尔拉,P。;Woodward,P.,气体动力学模拟的分段抛物型方法(PPM),J.Comput。物理。,54, 174-201 (1984) ·Zbl 0531.76082号
[17] Crockett,R.K。;科尔拉,P。;Fisher,R.T。;Klein,R.I。;McKee,C.F.,《理想MHD的非分裂、以细胞为中心的Godunov方法》,J.Compute。物理。,203, 422-448 (2005) ·Zbl 1143.76599号
[18] Dai,W。;Woodward,P.,多维磁流体力学的简单有限差分格式,计算机J。物理。,142, 331-369 (1998) ·Zbl 0932.76048号
[19] A.杜比。;Antypas,K。;Ganapathy,M.K。;里德·L.B。;莱利,K.M。;Sheeler,D。;西格尔,A。;Weide,K.,基于可扩展组件的FLASH架构,大规模并行多物理仿真代码,并行计算。,35, 10-11, 512-522 (2009)
[20] Dumbser,M。;Käser,M.,线性双曲方程组非结构网格上的任意高阶非振荡有限体积格式,J.Compute。物理。,221, 693-723 (2007) ·Zbl 1110.65077号
[21] Dumbser,M。;Käser,M。;Titarev,V.A。;Toro,E.F.,非线性双曲方程组非结构网格上的无四次非振荡有限体积格式,J.Compute。物理。,226, 204-243 (2007) ·Zbl 1124.65074号
[22] Dumbser,M。;巴尔萨拉,D。;托罗,E.F。;Munz,C.D.,在非结构化网格上构造一步有限体积和间断Galerkin格式的统一框架,J.Compute。物理。,227, 8209-8253 (2008) ·Zbl 1147.65075号
[23] Dumbser,M。;埃诺,C。;Toro,E.F.,刚性双曲平衡定律的高精度有限体积格式,J.Compute。物理。,227, 3971-4001 (2008) ·Zbl 1142.65070号
[24] 埃文斯,C.R。;Hawley,J.F.,《磁流体动力学流动模拟-约束传输方法》,天体物理学。J.,332659-677(1988)
[25] Fey,M.,多维上绕1。求解欧拉方程的传输方法,J.Compute。物理。,143, 159-180 (1998) ·Zbl 0932.76050号
[26] Fey,M.,《多维上绕2》。将欧拉方程分解为平流方程,J.Compute。物理。,143, 181-199 (1998) ·Zbl 0932.76051号
[28] 加德纳,T。;Stone,J.,通过约束传输实现理想MHD的非分裂Godunov方法,J.Compute。物理。,205, 509-539 (2005) ·Zbl 1087.76536号
[29] 加德纳,T。;Stone,J.,通过三维约束传输实现理想磁流体动力学的非分裂Godunov方法,J.Compute。物理。,227, 4123-4141 (2008) ·Zbl 1317.76057号
[30] 吉尔金,H。;劳伦斯,J。;Rosier,C.,线性双曲型系统的多维黎曼问题,数字注释。流体力学。,43, 284 (1993) ·Zbl 0921.35090号
[31] Hanawa,T。;米卡米,H。;Matsmoo,T.,《基于MHD方程特性改善激波不规则性》,J.Compute。物理。,227, 7952-7976 (2008) ·Zbl 1268.76038号
[32] Harten,A。;Osher,S。;Engquist,B。;Charkaravarthy,S.,一致高阶精确本质非振荡格式的一些结果,应用。数字。数学。,2, 347-377 (1987) ·Zbl 0627.65101号
[33] Harten,A。;Engquist,B。;Osher,S。;Charkaravarthy,S.,《一致高阶基本无振荡格式III》,J.Compute。物理。,71, 231-303 (1987) ·Zbl 0652.65067号
[34] Helzel,C。;Rossmanith,J.A。;Taetz,B.,三维理想磁流体动力学方程的非聚集约束输运方法,J.Compute。物理。,230, 3803-3829 (2011) ·Zbl 1369.76061号
[35] 胡,C。;Shu,C.W.,三角形网格上的加权本质非振荡格式,J.Compute。物理。,150, 97-127 (1999) ·Zbl 0926.65090号
[36] 江,G。;Shu,C.W.,加权ENO方案的高效实现,J.Compute。物理。,126, 202-228 (1996) ·Zbl 0877.65065号
[37] Kraichnan,R.H.,二维湍流中的惯性范围,物理学。流体,101417-1423(1967)
[38] Lee,D。;Deane,A.E.,多维磁流体力学的非分裂交错网格格式,J.Compute。物理。,228, 952-975 (2009) ·Zbl 1330.76093号
[42] 刘晓东。;Osher,S。;Chan,T.,加权基本非振荡格式,J.Compute。物理。,115, 200-212 (1994) ·Zbl 0811.65076号
[43] Londrilo,P。;Del Zanna,L.,《关于理想磁流体力学Godunov型格式中的无发散条件:迎风约束输运方法》,J.Compute。物理。,195, 14-48 (2004) ·Zbl 1087.76074号
[44] 米农,A。;Tzeferacos,P.,《以细胞为中心的MHD的二阶非分裂Godunov方案——CTU-GLM方案》,J.Compute。物理。,229, 2117-2138 (2010) ·Zbl 1303.76142号
[45] 米尼亚蒂,F。;Martin,D.F.,CHARM中自适应网格细化的约束传输磁流体力学,天体物理学。J.Suppl.,195,5-25(2011)
[46] 三好,T。;Kusano,K.,理想磁流体力学的多状态HLL近似黎曼解算器,J.Compute。物理。,208, 315-344 (2005) ·Zbl 1114.76378号
[47] Orszag,S.A.公司。;唐春明,二维磁流体动力学湍流的小尺度结构,流体力学杂志。,90, 129-143 (1979)
[48] 潘多尔菲,M。;D’Ambrosio,D.,迎风方法中的数值不稳定性:红玉现象的分析和治疗,J.Compute。物理。,166, 271-301 (2001) ·Zbl 0990.76051号
[50] Quirk,J.,《对伟大的Riemann解算器辩论的贡献》,国际J.Numer出版社。液体方法,18,555-574(1994)·Zbl 0794.76061号
[51] Roe,P.L.,近似黎曼解算器,参数向量和差分格式,J.Compute。物理。,43, 357-372 (1981) ·Zbl 0474.65066号
[52] Rumsey,C.L。;van Leer,B。;Roe,P.L.,一个多维通量函数及其在欧拉方程和纳维-斯托克斯方程中的应用,计算杂志。物理。,105, 306-323 (1993) ·Zbl 0767.76039号
[53] Saltzman,J.,《双曲守恒定律的非分裂三维迎风方法》,J.Compute。物理。,115153-168(1994年)·Zbl 0813.65111号
[54] 桑德斯,R。;莫拉诺,E。;Druguet,M.,迎风格式的多维耗散:稳定性和气体动力学应用,J.计算。物理。,145, 511-537 (1998) ·Zbl 0924.76076号
[55] 舒,C.-W。;Osher,S.J.,《本质上非振荡激波捕获方案的有效实现》,J.Compute。物理。,77, 439-471 (1988) ·兹伯利0653.65072
[56] 舒,C.-W。;Osher,S.J.,《本质上非振荡激波捕获方案的有效实现II》,J.Compute。物理。,83,32-78(1989年)·Zbl 0674.65061号
[57] 斯通,J。;加德纳,T。;Teuben,P。;霍利,J。;Simon,J.,《雅典娜:天体物理MHD的新代码,营养不良》。《补充期刊》,178137-177(2008)
[58] Suresh,A。;Huynh,H.T.,带Runge-Kutta时间步进的精确单调性保持格式,J.Compute。物理。,136, 83-99 (1997) ·Zbl 0886.65099号
[59] Titarev,V.A。;Toro,E.F.,ADER:任意高阶Godunov方法,科学杂志。计算。,17, 1-4, 609-618 (2002) ·兹比尔1024.76028
[60] Titarev,V.A。;Toro,E.F.,三维非线性双曲方程组的ADER格式,J.Compute。物理。,204, 715-736 (2005) ·Zbl 1060.65641号
[61] Toro,E.F.,Riemann解算器和流体动力学数值方法,实用简介(2009),Springer·Zbl 1227.76006号
[62] 托罗,E.F。;Titarev,V.A.,《守恒定律和ADER方法系统的导数黎曼解算器》,J.Compute。物理。,212, 1, 150-165 (2006) ·Zbl 1087.65590号
[63] Tóth,G.,冲击捕获磁流体动力学代码中的约束,J.Compute。物理。,161, 605-656 (2000) ·Zbl 0980.76051号
[64] van Leer,B.,《可压缩流的迎风和高分辨率方法:从供体细胞到残余分布方案》,Commun。计算。物理。,1, 192-206 (2006) ·Zbl 1114.76049号
[65] Wendroff,B.,气体动力学的二维HLLE-Riemann解算器和相关Godunov型差分格式,计算。数学。申请。,38, 175-185 (1999) ·Zbl 0984.76064号
[66] Zachary,A。;Malagoli,A。;Colella,P.,多维理想磁流体力学的高阶Godunov方法,J.Compute。物理。,2, 263-284 (1994) ·Zbl 0797.76063号
[67] Zhang,Y.-T。;Shu,C.-W.,三角网格上Hamilton-Jacobi方程的高阶WENO格式,SIAM J.Sci。计算。,24, 1005-1030 (2003) ·Zbl 1034.65051号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。