×

一种完全隐式、保守、非线性、电磁混合粒子-离子/流体-电子算法。 (英语) Zbl 1416.76189号

摘要:含有动力学离子和流体电子的准中性混合模型是一种很有希望的方法,可以用于桥接空间和实验室等离子体中许多问题的固有多尺度性质。在这里,针对具有多离子物种的全三维电磁问题,导出了一种新的基于粒子-细胞的隐式混合模型方案,该方案具有全局质量、动量和能量守恒的特点。该方案包括离子的亚循环和轨道平均,以细胞为中心的有限差分和隐式中点时间推进。为了减少离散粒子噪声,该方案允许任意阶形状函数用于粒子-网格插值和保守二项式平滑的应用。通过对大量测试问题的验证,证明了新方案实现的正确性、独特的守恒性和良好的稳定性。特别是,在几种常用的非保守方案高度不稳定的情况下,没有迹象表明漂移通过均匀空间网格的冷离子群的有限网格不稳定性不稳定增长。

MSC公司:

76M20码 有限差分方法在流体力学问题中的应用
76米25 其他数值方法(流体力学)(MSC2010)
76×05 电磁场中的电离气体流动;浆流
6500万06 含偏微分方程初值和初边值问题的有限差分方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 温斯克,D。;尹,L。;Omidi,N。;Karimabadi,H。;Quest,K.,《混合仿真代码:过去、现在和未来——教程》,(Büchner,J.;Dum,C.;Scholer,M.,《空间等离子体模拟》,《空间等离子模拟》,物理讲义,第615卷(2003),柏林斯普林格出版社),136-165
[2] Lipatov,A.S.,《混合多尺度模拟技术:天体物理和实验室等离子体应用简介》,科学计算(2002年),柏林-海德堡施普林格出版社·兹比尔1015.76001
[3] Karimabadi,H。;罗特斯提恩,V。;Vu,H.X。;Omelchenko,Y.A。;斯库德,J。;多顿,W。;丁莫克,A。;Nykyri,K。;Wan,M。;Sibeck,D。;Tatineni,M。;马久姆达尔,A。;Loring,B。;Geveci,B.,《无碰撞等离子体中激波、湍流和磁重联之间的联系》,Phys。Plasmas,21,6,第062308条pp.(2014)
[4] Stanier,A。;多顿,W。;Chacón,L。;Karimabadi,H。;Ng,J。;黄,Y.-M。;哈基姆,A。;Bhattacharjee,A.,离子动力学物理学在磁通绳相互作用中的作用,Phys。修订稿。,115、17、第175004条pp.(2015)
[5] Sturdevant,B.J。;陈,Y。;Parker,S.E.,环形离子温度梯度不稳定性的低频全动力学模拟,Phys。Plasmas,24,8,第081207条,pp.(2017)
[6] Birdsall,C.K。;Langdon,A.B.,《通过计算机模拟的等离子体物理》(1991年),McGraw-Hill:McGraw-Hill纽约
[7] 霍克尼,R.W。;Eastwood,J.W.,《使用粒子的计算机模拟》(1981年),McGraw-Hill:McGraw-Hill纽约·Zbl 0662.76002号
[8] Karimabadi,H。;Krauss-Varban,D。;胡巴,J.D。;Vu,H.X.,《磁场重联机制和相关三维不对称:混合、无霍尔混合和霍尔-MHD模拟》,J.Geophys。研究空间物理学。,109,第A09205条pp.(2004)
[9] Harned,D.S.,宏观等离子体现象的准中性混合模拟,J.Compute。物理。,47552-462(1982年)·Zbl 0493.76122号
[10] 温斯克,D。;Quest,K.B.,《电磁离子束不稳定性——一维和二维模拟的比较》,J.Geophys。研究,91,8789-8797(1986)
[11] Kunz,M.W。;斯通,J.M。;Bai,X.N.,Pegasus:天体物理等离子体动力学的新混合动力学粒子-细胞代码,J.Compute。物理。,259, 154-174 (2014) ·Zbl 1349.82149号
[12] 温斯克,D。;Quest,K.B.,垂直超临界无碰撞冲击下的磁场和密度波动,J.Geophys。研究,93,9681-9693(1988)
[13] Matthews,A.P.,2D多光谱混合等离子体模拟的当前前进方法和循环跳跃,计算机杂志。物理。,112, 1, 102-116 (1994) ·Zbl 0801.76057号
[14] Fujimoto,M.,磁层速度剪切层中的不稳定性(1990),东京大学:东京大学,博士论文
[15] 托马斯,V.A。;温斯克,D。;Omidi,N.,重塑超临界准平行激波:1。一维和二维模拟,J.Geophys。研究空间物理学。,95,A1118809-18819(1990)
[16] Kunz,M.W。;Lesur,G.,《霍尔主导磁旋转湍流中的磁自组织》,孟买。不是。R.阿斯顿。Soc.,434,2295-2312(2013)
[17] Swift,D.W.,使用混合代码模拟地球磁层,Geophys。Res.Lett.公司。,22311-314(1995年)
[18] Hewett,D.W.,在零惯性-电子混合等离子体模拟代码中求解电子场方程的全局方法,J.Compute。物理。,38, 378-395 (1980) ·Zbl 0462.76107号
[19] 天野之弥。;东岛,K。;白川方明,K.,《无碰撞等离子体混合模拟中处理低密度区域的稳健方法》,J.Compute。物理。,275, 197-212 (2014) ·Zbl 1349.82068号
[20] Cheng,J。;帕克,S.E。;陈,Y。;Uzdensky,D.A.,混合模拟的二阶半隐式(δ)f方法,J.Compute。物理。,245, 364-375 (2013) ·Zbl 1349.76595号
[21] Sturdevant,B.J。;帕克,S.E。;陈,Y。;Hause,B.B.,《洛伦兹离子亚循环和轨道平均的隐式颗粒-细胞方法》,J.Compute。物理。,316519-533(2016)·Zbl 1349.82153号
[22] 陈,G。;Chacón,L。;Barnes,D.C.,《能量和电荷守恒的隐式静电粒子-细胞算法》,J.Compute。物理。,230, 7018-7036 (2011) ·Zbl 1237.78006号
[23] Markidis,S.等人。;Lapenta,G.,《节能颗粒-细胞法》,J.Compute。物理。,230, 7037-7052 (2011) ·Zbl 1231.82067号
[24] Lewis,H.R.,Vlasov等离子体的能量守恒数值近似,J.Compute。物理。,6, 136-141 (1970) ·Zbl 0205.57801号
[25] 陈,G。;Chacón,L.,《多维能量和电荷守恒非线性隐式电磁Vlasov-Darwin粒子-细胞算法》,计算。物理。社区。,19773-87(2015年)·Zbl 1352.65405号
[26] Chacón,L。;Chen,G.,多维中的曲线、完全隐式、保守电磁PIC算法,J.Compute。物理。,316, 578-597 (2016) ·Zbl 1349.82139号
[27] Squire,J。;秦,H。;Tang,W.M.,Vlasov-Maxwell系统与变分粒子-细胞方案的几何积分,Phys。Plasmas,第19、8章,第084501页(2012年)
[28] Evstatiev,E.G。;Shadwick,B.A.,动力学等离子体模拟粒子算法的变分公式,J.Compute。物理。,245, 376-398 (2013) ·Zbl 1349.82078号
[29] 北卡罗来纳州克鲁塞尔斯。;Einkemmer,L。;Faou,E.,Vlasov-Maxwell方程的哈密顿分裂,J.Comput。物理。,283, 224-240 (2015) ·Zbl 1351.35223号
[30] 秦,H。;刘,J。;肖,J。;张,R。;何毅。;Wang,Y。;孙,Y。;Burby,J.W。;Ellison,L。;Zhou,Y.,Vlasov-Maxwell方程长期大规模模拟的正则辛粒子-细胞方法,Nucl。Fusion,56,1,第014001条pp.(2016)
[31] 克劳斯,M。;科尔曼,K。;莫里森,P.J。;Sonnendrucker,E.,《Gempic:几何电磁粒子-细胞方法》,J.Plasma Phys。,第83、4条,第905830401页(2017年)
[32] Okuda,H.,空间网格导致的等离子体模拟中的非物理噪声和不稳定性,J.Compute。物理。,10, 475-486 (1972)
[33] Brackbill,J.U.,《粒子-细胞等离子体模拟中的能量和动量守恒》,J.Compute。物理。,317405-427(2016)·Zbl 1349.76916号
[34] 弗拉德,G。;布里古利奥,S。;Fogaccia,G。;Martino,B.D.,无网格有限尺寸粒子等离子体模拟,计算机。物理。社区。,134, 58-77 (2001) ·Zbl 0998.76071号
[35] 黄,C.-K。;曾勇。;Wang,Y。;医学博士迈耶斯。;Yi,S。;奥尔布赖特,B.J.,《静电粒子-细胞算法的有限网格不稳定性和光谱保真度》,计算。物理。社区。,207, 123-135 (2016) ·Zbl 1375.76146号
[36] Rambo,P.W.,准中性混合模拟中的有限网格不稳定性,J.Compute。物理。,118152-158(1995年)·Zbl 0824.76056号
[37] 穆尼奥斯,P.A。;Jain,N。;Kilian,P。;Büchner,J.,一种新的无需近似即可实现惯性电子流体方程的混合代码(CHIEF),ArXiv e-prints·Zbl 07694308号
[38] Le,A。;多顿,W。;Karimabadi,H。;Egedal,J.,《磁重联与动力学离子和流体电子压力各向异性的混合模拟》,Phys。Plasmas,23,3,文章032114 pp.(2016)
[39] Ohia,O。;埃杰达尔,J。;Lukin,V.S。;多顿,W。;Le,A.,捕获磁重联动力学结构的各向异性流体闭合演示,Phys。修订稿。,109,11,第115004条第(2012)页
[40] Ng,J。;黄,Y.-M。;哈基姆,A。;巴塔查吉,A。;Stanier,A。;多顿,W。;王,L。;Germaschewski,K.,《岛屿合并问题:包括高阶矩在内的扩展流体模型中重联的缩放》,Phys。Plasmas,22,11,第112104条pp.(2015)
[41] Chacón,L.,曲线几何中三维隐式扩展磁流体力学的非交错、保守、div(B)=0的有限体积格式,计算。物理。社区。,163, 143-171 (2004) ·Zbl 1196.76040号
[42] 维拉塞纳,J。;Buneman,O.,局部电磁场解算器的严格电荷守恒,计算机。物理。社区。,69, 306-316 (1992)
[43] Mardahl,P.J。;Verbonceur,J.P.,电磁PIC码中的电荷守恒;Boris/DADI和Langdon Marder方法的光谱比较,Comput。物理。社区。,106, 219-229 (1997) ·Zbl 0939.78009号
[44] Esirkepov,T.,《具有任意形状因子的颗粒-细胞模拟的精确电荷守恒方案》,计算。物理。社区。,135, 2, 144-153 (2001) ·Zbl 0981.78014号
[45] Knoll,D.A。;Keyes,D.E.,《无雅可比牛顿-克利洛夫方法:方法和应用调查》,J.Compute。物理。,193, 357-397 (2004) ·Zbl 1036.65045号
[46] Chacón,L.,三维粘阻磁流体力学的最优、并行、全隐式Newton-Krylov解算器,Phys。Plasmas,第15、5页,第056103条,pp.(2008)
[47] Saad,Y.,《一种灵活的内外预处理GMRES算法》,SIAM J.Sci。计算。,14, 2, 461-469 (1993) ·Zbl 0780.65022号
[48] 科恩,B.I。;弗雷斯·R·P。;Thomas,V.,《轨道平均隐式粒子代码》,J.Compute。物理。,45, 3, 345-366 (1982) ·Zbl 0478.65072号
[49] Vu,H.等人。;Brackbill,J.,《磁场中带电粒子运动的精确数值解》,J.Compute。物理。,116, 2, 384-387 (1995) ·Zbl 0824.65130号
[50] Hirt,C.W.,有限差分方程的启发式稳定性理论,J.Compute。物理。,2, 339-355 (1968) ·Zbl 0187.12101号
[51] Gaskell,P.H。;Lau,A.K.C.,《曲率补偿对流传输:SMART,一种新的有界性保护传输算法》,国际数值杂志。液体方法,8617-641(1988)·Zbl 0668.76118号
[52] Gary,S.P.,空间等离子体微不稳定性理论(1993)
[53] 托尔德,D。;Cookmeyer,J。;Astfalk,P。;Jenko,F.,等离子体物理混合动离子/流电子模型的线性色散关系,新物理学杂志。,18,7,文章075001第(2016)页
[54] 价格,成本加运费。;斯威夫特,D.W。;Lee,L.-C.,地球磁鞘非振荡镜波的数值模拟,J.Geophys。研究空间物理学。,91、A1、101-112(1986)
[55] 加里,S.P。;Fuselier,S.A。;Anderson,B.J.,磁鞘中的离子各向异性不稳定性,J.Geophys。Res.,98,1481-1488(1993)
[56] 休伯特,D。;Perche,C。;C.C.哈维。;拉科姆,C。;Russell,C.T.,准垂直激波下游镜像波的观测,地球物理学。Res.Lett.公司。,1619-162(1989年)
[57] Birn,J。;Drake,J.F。;Shay,医学硕士。;罗杰斯,B.N。;丹顿,R.E。;黑塞,M。;库兹涅佐娃,M。;马振伟。;巴塔查吉,A。;奥托,A。;Pritchett,P.L.,《地球空间环境建模(GEM)磁重联挑战》,J.Geophys。决议,106,3715-3720(2001)
[58] 帕克,S.E。;Lee,W.W.,《回转运动模拟的完全非线性特征方法》,Phys。流体B,5,77-86(1993)
[59] 王,B。;Miller,G.H。;Colella,P.,《动力学等离子体的自适应相空间重映射粒子-细胞方法》,SIAM J.Sci。计算。,33, 6, 3509-3537 (2011) ·兹比尔1232.76046
[60] 佐科,A。;Chacón,L。;Simakov,A.N.,《物理学》。Plasmas,第16、11章,第110703页(2009年)
[61] Stanier,A。;Simakov,A.N。;Chacón,L。;Daughton,W.,《带强导向场的流体与动力磁重联》,Phys。Plasmas,22,10,第101203条pp.(2015)
[62] Chacón,L。;Knoll,D.A.,《二维高(β)霍尔MHD隐式非线性解算器》,J.Compute。物理。,188, 573-592 (2003) ·Zbl 1127.76375号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。