Debojyoti Ghosh；米哈伊尔·多夫（Mikhail A.Dorf）。；米洛·多尔。；杰弗里·A·F·希丁格。 半隐式时间积分碰撞磁化等离子体的动力学模拟。 （英语） Zbl 1407.65146号 科学杂志。计算。 77，第2期，819-849（2018）. 概述：具有不同碰撞性的等离子体在许多应用中都会发生，例如托卡马克边缘区域，在那里，流动的特点是密度和温度发生显著变化。虽然弱碰撞高温等离子体需要动力学模型，但较冷区域的高碰撞性使得方程由于不同的时间尺度而在数值上变得僵硬。本文针对这种情况，我们提出了一种隐式-显式算法，其中碰撞项在时间上隐式积分，而平流项在时间中显式积分，从而允许与平流时间尺度相当的时间步长。与使用显式时间积分器的算法相比，这种划分产生了更高效的算法，其中时间步长受到刚性碰撞时间尺度的限制。我们在COGENT（一种高阶有限体积回转动力学代码）中实现了半隐式可加性Runge-Kutta方法，并测试了这些半隐式方法在高碰撞等离子体测试用例中的准确性、收敛性和计算成本。 引用于4文件 MSC公司： 6500万08 含偏微分方程初值和初边值问题的有限体积法 65升06 常微分方程的多步、Runge-Kutta和外推方法 第35季度83 弗拉索夫方程 84年第35季度 福克-普朗克方程 82天75 核反应堆理论；中子输运 82D10号 等离子体统计力学 82立方31 随机方法（福克-普朗克、朗之万等）应用于含时统计力学问题 关键词：IMEX时间集成；等离子体物理学；回转运动模拟；Vlasov-Fokker-Planck方程 软件：ELMFIRE公司；GYSELA公司；FPPAC公司；硝基苯砜；炒作；COGENT公司 PDF格式BibTeX公司 XML格式引用 \textit｛D.Ghosh｝等人，《科学杂志》。计算。77，No.2，819--849（2018；Zbl 1407.65146） 全文： DOI程序 arXiv公司 参考文献： [1] 阿舍尔，UM；Ruuth，SJ；Spiteri，RJ，含时偏微分方程的隐式显式Runge-Kutta方法，应用。数字。数学。，25, 151-167, (1997) ·Zbl 0896.65061号 ·doi:10.1016/S0168-9274（97）00056-1 [2] 银行，JW；Brunner，S。；Berger，RL；Tran，TM，Vlasov电子-离子碰撞对电子等离子体波阻尼影响的模拟，Phys。等离子，23，032108，（2016）·数字对象标识代码：10.1063/1.4943194 [3] EA Belli；Candy，J.，《新古典输运模拟中的完全线性化福克-普朗克碰撞》，等离子体物理学。控制。融合，54015015，（2012）·doi:10.1088/0741-3335/54/015015 [4] Berezin，Y。；库迪克，V。；Pekker，M.，Fokker-Planck方程的保守有限差分格式，不违反熵增定律，J.Compute。物理。，69, 163-174, (1987) ·Zbl 0644.76091号 ·doi:10.1016/0021-9991（87）90160-4 [5] Braginskii，SI，血浆中的转运过程，血浆物理学评论。，1, 205, (1965) [6] Brunner，S.，Tran，T.，Hittinger，J.：欧拉-弗拉索夫模拟中非线性Landau碰撞算子的数值实现。第一部分：Rosenbluth势的计算。技术代表LLNL-SR-459135，加利福尼亚州利弗莫尔劳伦斯·利弗莫雷国家实验室（2010） [7] 布埃特，C。；Cordier，S.，各向同性Fokker-Planck-Landau方程的守恒和熵衰减数值格式，J.Compute。物理。，145, 228-245, (1998) ·Zbl 0917.76045号 ·doi:10.1006/jcph.1998.6015 [8] 布埃特，C。；科迪尔，S。；德贡，P。；Lemou，M.，Fokker-Planck-Landau方程数值、保守和熵近似的快速算法，J.Compute。物理。，133, 310-322, (1997) ·Zbl 0880.65112号 ·文件编号：10.1006/jcph.1997.5669 [9] Butcher，J.：常微分方程的数值方法。威利，纽约（2003年）·Zbl 1040.65057号 ·doi:10.1002/0470868279 [10] 坎迪，J。；Waltz，RE，超级计算机模拟匹配的DIII-D托卡马克中的异常输运定标，Phys。修订稿。，91, 045001, (2003) ·doi:10.1103/PhysRevLett.91.045001 [11] 卡萨诺娃，M。；Larroche，O。；Matte，JP，等离子体中碰撞冲击波的动力学模拟，物理。修订稿。，67, 2143-2146, (1991) ·doi:10.1103/PhysRevLett.67.2143 [12] Chacón，L。；巴恩斯特区；诺尔，DA；Miley，GH，《隐式能量守恒2D Fokker-Planck算法：I.差分格式》，J.Compute。物理。，157, 618-653, (2000) ·Zbl 0961.76057号 ·doi:10.1006/jcph.1999.6394 [13] Chacón，L。；巴恩斯特区；诺尔，DA；Miley，GH，隐式能量守恒2D Fokker-Planck算法：II。无Jacobian牛顿-克利洛夫解算器，J.Compute。物理。，157, 654-682, (2000) ·Zbl 0961.76058号 ·doi:10.1006/jcph.1999.6395 [14] Chang，CS；Ku，S.，静态托卡马克边缘等离子体中的自发旋转源，物理学。等离子体，15062510，（2008）·数字对象标识代码：10.1063/1.2937116 [15] Chang，J。；Cooper，G.，福克-普朗克方程的实用差分格式，J.Compute。物理。，6, 1-16, (1970) ·Zbl 0221.65153号 ·doi:10.1016/0021-9991（70）90001-X [16] 陈，G。；Chacón，L.，《能量和电荷守恒，非线性隐式，电磁1D-3V Vlasov-Darwin粒子-细胞算法》，计算。物理。社区。，185, 2391-2402, (2014) ·Zbl 1360.78014号 ·doi:10.1016/j.cpc.2014.05.010 [17] 陈，G。；Chacón，L.，《多维能量和电荷守恒非线性隐式电磁Vlasov-Darwin粒子-细胞算法》，计算。物理。社区。，197, 73-87, (2015) ·Zbl 1352.65405号 ·doi:10.1016/j.cpc.2015.08.008 [18] 陈，G。；Chacón，L。；Barnes，D.，一种能量和电荷守恒的隐式静电粒子在细胞中的算法，J.Comput。物理。，2307018-7036，（2011年）·Zbl 1237.78006号 ·doi:10.1016/j.jp.2011.05.031 [19] Cheng，C。；Knorr，G.，《组态空间中Vlasov方程的积分》，J.Compute。物理。，22, 330-351, (1976) ·doi:10.1016/0021-9991（76）90053-X [20] 科恩，RH；Xu，XQ，边缘等离子体动力学模拟进展，等离子体物理。，48, 212-223, (2008) ·doi:10.1002/ctp.200810038 [21] 科尔拉，P。；多尔，M。；Hittinger，J。；Martin，D.，映射坐标中的高阶有限体积方法，J.Compute。物理。，230, 2952-2976, (2011) ·Zbl 1218.65119号 ·doi:10.1016/j.jcp.2010.12.044 [22] 北卡罗来纳州克鲁塞尔斯。；Respaud，T。；Sonnendrücker，E.，Vlasov方程数值解的正向半拉格朗日方法，计算。物理。社区。，180, 1730-1745, (2009) ·Zbl 1197.82012年 ·doi:10.1016/j.cpc.2009.04.024 [23] Dennis，J.，Schnabel，R.：无约束优化和非线性方程的数值方法。费城工业和应用数学学会（1996年）。https://doi.org/10.1137/1.9781611971200 ·兹伯利0847.65038 ·doi:10.1137/1.9781611971200 [24] 马萨诸塞州多夫；科恩，RH；康普顿，JC；多尔，M。；罗格利恩，TD；安格斯，J。；Krasheninikov，S。；科尔拉，P。；马丁·D·。；McCorquodale，P.，《COGENT边缘动力学代码的进展：碰撞算子选项》，《控制等离子体物理》。，52, 518-522, (2012) ·doi:10.1002/ctpp.201210042 [25] 马萨诸塞州多夫；科恩，RH；多尔，M。；Hittinger，J。；TD Rognlien，《COGENT边缘动力学代码的进展：实现福克-普朗克碰撞算子》，Contrib.Plasma Phys。，54, 517-523, (2014) ·doi:10.1002/cpp.201410023 [26] 多尔夫，马萨诸塞州；科恩，RH；多尔，M。；Rognlien，T。；Hittinger，J。；康普顿，J。；科尔拉，P。；马丁·D·。；McCorquodale，P.，用连续回转动力学代码COGENT模拟新古典输运，物理学。等离子，20，012513，（2013）·doi:10.1063/1.4776712 [27] Dorr，M.R.，Colella，P.，Dorf，M.A.，Ghosh，D.，Hittinger，J.A.F.，Schwartz，P.O.：边缘等离子体几何中旋回动力学Vlasov模型的高阶离散化。提交日期（2017年）。arXiv:1712.01978年 [28] 杜兰，DR；Blossey，PN，快波-慢波问题的隐式-显式多步方法，Mon。《天气评论》，140，1307-1325，（2012）·doi:10.1175/MWR-D11-00088.1 [29] Epperlein，E.，Fokker-Planck方程的隐式和保守差分格式，J.Compute。物理。，112, 291-297, (1994) ·Zbl 0806.76050号 ·doi:10.1006/jcph.1994.1101 [30] 埃珀林，EM；GJ Rickard；贝尔，AR，激光等离子体中的二维非局域电子输运，物理学。修订稿。，61, 2453-2456, (1988) ·doi:10.1103/PhysRevLett.61.2453 [31] Falgout，R.D.，Yang，U.M.：《炒作：高性能预处理器库》，第632-641页。施普林格，柏林（2002）。https://doi.org/10.1007/3-540-47789-6_66 ·Zbl 1056.65046号 ·doi:10.1007/3-540-47789-666 [32] 费尔贝特，F。；Pareschi，L.，非齐次情况下Fokker-Planck-Landau方程精确解的数值方法，J.Compute。物理。，179, 1-26, (2002) ·Zbl 1003.82011年 ·doi:10.1006/jcph.2002.7010 [33] 费尔贝特，F。；Sonnendrücker，E。；Bertrand，P.，Vlasov方程的保守数值格式，J.Compute。物理。，172, 166-187, (2001) ·Zbl 0998.65138号 ·doi:10.1006/jcph.2001.6818 [34] Ghosh，D。；Constantinescu，EM，《基于特征通量分配的大气流半隐式时间积分》，SIAM J.Sci。计算。，38，a1848-a1875，（2016）·Zbl 1419.76478号 ·doi:10.1137/15M1044369 [35] 吉拉尔多，FX；凯利，JF；Constantinescu，E.，三维非流体静力统一大气模型（NUMA）的隐式-显式公式，SIAM J.Sci。计算。，35，b1162-b1194，（2013）·Zbl 1280.86008号 ·doi:10.137/120876034 [36] 吉拉尔多，FX；Restelli，M。；Läuter，M.，《Navier-Stokes方程的半隐式公式：在非静力大气建模中的应用》，SIAM J.Sci。计算。，32, 3394-3425, (2010) ·Zbl 1237.76153号 ·doi:10.1137/090775889 [37] 格兰吉拉德，V。；布鲁内蒂，M。；Bertrand，P。；北贝塞。；加贝，X。；Ghendrih，P。；Manfredi，G。；萨拉津，Y。；O.绍特。；Sonnendrücker，E。；瓦茨拉维克，J。；Villard，L.，离子湍流模拟的漂移动力学半拉格朗日4D代码，J.Compute。物理。，217, 395-423, (2006) ·Zbl 1160.76385号 ·doi:10.1016/j.jcp.2006.01.023 [38] Hahm，TS，核心输运屏障湍流的非线性回转动力学方程，物理学。等离子体，34658-4664（1996）·数字对象标识代码：10.1063/1.872034 [39] 海基宁，J。；Kiviniemi，T。；Kurki Suonio，T。；Peeters，A。；Sipilä，S.，《新古典等离子体的粒子模拟》，J.Compute。物理。，173, 527-548, (2001) ·Zbl 1051.82522号 ·doi:10.1006/jcph.2001.6891 [40] 日本海基宁；Henriksson，S。；Janhunen，S。；Kiviniemi，TP；Ogando，F.，边缘等离子体中存在宽轨道和强剖面变化时粒子和热传输的回转动力学模拟，Contrib.plasma Phys。，46, 490-495, (2006) ·doi:10.1002/ctpp.200610035 [41] Huba，J.D.：NRL血浆处方。华盛顿特区海军研究实验室技术代表（2016年） [42] Y.小村。；艾达，M。；卡诺，T。；Aiba，N。；Tokuda，S.，保守型全球回转动力学环形全f五维Vlasov模拟，计算。物理。社区。，179, 391-403, (2008) ·Zbl 1197.82110号 ·doi:10.1016/j.cpc.2008.04.005 [43] Y.小村。；艾达，M。；Tokuda，S.，保守旋翼动力学Vlasov模拟，Commun。非线性科学。数字。模拟。，13227-233（2008年）·兹比尔1130.35368 ·doi:10.1016/j.cnsns.2007.05.015 [44] Jackson，J.D.：《经典电动力学》，第三版。威利，纽约（1999）·Zbl 0920.00012号 [45] 雅各布斯，G。；Hesthaven，J.，具有双曲散度清理的高阶粒子-细胞方法的隐式-显式时间积分，计算。物理。社区。，180, 1760-1767, (2009) ·Zbl 1197.76081号 ·doi:10.1016/j.cpc.2009.05.020 [46] James，R.，孤立源分布泊松方程的解，J.Compute。物理。，25, 71-93, (1977) ·Zbl 0373.65051号 ·doi:10.1016/0021-9991（77）90013-4 [47] 卡迪奥格鲁，SY；诺尔，DA；罗伊，RB；Rauenzahn，RM，辐射流体力学的二阶自洽IMEX方法，J.Compute。物理。，229, 8313-8332, (2010) ·Zbl 1381.76262号 ·doi:10.1016/j.jcp.2010.07.019 [48] 加利福尼亚州肯尼迪；Carpenter，MH，对流-扩散-反应方程的可加Runge-Kutta格式，应用。数字。数学。，44, 139-181, (2003) ·Zbl 1013.65103号 ·doi:10.1016/S0168-9274（02）00138-1 [49] Kho，TH，带电粒子系统的弛豫，物理学。版本A，32，666-669，（1985）·doi:10.1103/PhysRevA.32.666 [50] Kingham，R。；Bell，A.，用隐式Vlasov-Fokker-Planck代码模拟二维磁场中的非局域电子输运，J.Compute。物理。，194, 1-34, (2004) ·Zbl 1136.76400号 ·doi:10.1016/j.jcp.2003.08.017 [51] Knoll，D。；Keyes，D.，《无Jacobian牛顿-克利洛夫方法：方法和应用调查》，J.Compute。物理。，193, 357-397, (2004) ·Zbl 1036.65045号 ·doi:10.1016/j.jcp.2003.08.010 [52] Kumar，H。；Mishra，S.，双流体等离子体方程的熵稳定数值格式，J.Sci。计算。，52, 401-425, (2012) ·Zbl 1311.76150号 ·doi:10.1007/s10915-011-9554-7 [53] 朗道，LD，库仑相互作用情况下的动力学方程，Zh。埃克斯佩尔。i特奥雷特。Fiz.公司。，7, 203-209, (1937) ·Zbl 0016.28602号 [54] Landreman，M。；Ernst，DR，《连续体动力学计算和福克-普朗克碰撞的新速度空间离散化》，J.Compute。物理。，243, 130-150, (2013) ·Zbl 1349.82094号 ·doi:10.1016/j.jcp.2013.02.041 [55] Larroche，O.，等离子体碰撞实验的动力学模拟，Phys。流体B，52816-2840，（1993）·数字对象标识代码：10.1063/1.860670 [56] Larroche，O.，《ICF靶内爆中燃料离子传输的动力学模拟》，《欧洲物理学》。J.D At.Mol.Opt.公司。血浆物理学。，27, 131-146, (2003) ·doi:10.1140/epjd/e2003-00251-1 [57] Larsen，E。；利弗莫尔，C。；Pomraning，G。；Sanderson，J.，《一维Fokker-Planck算子的离散化方法》，J.Compute。物理。，61, 359-390, (1985) ·Zbl 0581.65091号 ·doi:10.1016/0021-9991（85）90070-1 [58] Lee，WW，粒子模拟中的回旋动力学方法，物理。流体，26556-562，（1983）·Zbl 0576.76120号 ·数字对象标识代码：10.1063/1.864140 [59] Lemou，M.，Fokker-Planck-Landau算子的多极展开，数值。数学。，78, 597-618, (1998) ·Zbl 0898.65090号 ·doi:10.1007/s00211005327 [60] 勒莫，M。；Mieussens，L.，Fokker-Planck-Landau方程的快速隐式格式，C.R.数学。，338, 809-814, (2004) ·Zbl 1081.82022号 ·doi:10.1016/j.crma.2004.03.010 [61] 勒莫，M。；Mieussens，L.，Fokker-Planck-Landau方程的隐式格式，SIAM J.Sci。计算。，27, 809-830, (2005) ·兹比尔1096.82015 ·doi:10.1137/040609422 [62] Maeyama，S。；石泽，A。；渡边，TH；Nakajima，N。；Tsuji-Iio，S。；Tsutsui，H.，回转动力学Vlasov模拟的半拉格朗日和加性半隐式Runge-Kutta格式的混合方法，计算。物理。社区。，183, 1986-1992, (2012) ·doi:10.1016/j.cpc.2012.04.028 [63] McCorquodale，P。；多尔，M。；Hittinger，J。；Colella，P.，映射多块网格上双曲守恒律的高阶有限体积方法，J.Compute。物理。，288, 181-195, (2015) ·Zbl 1354.65178号 ·doi:10.1016/j.jp.2015.01.006 [64] McCoy，M。；米林，A。；Killeen，J.，FPPAC：二维多物种非线性Fokker-Planck软件包，计算。物理。社区。，24, 37-61, (1981) ·doi:10.1016/0010-4655（81）90105-3 [65] 穆索，V。；Knoll，D.，碰撞等离子体的完全隐式动力学解，J.Compute。物理。，136, 308-323, (1997) ·Zbl 0896.76057号 ·doi:10.1006/jcph.1997.5736 [66] Pareschi，L。；Russo，G.，隐式显式Runge-Kutta格式及其在带松弛双曲方程组中的应用，J.Sci。计算。，25, 129-155, (2005) ·Zbl 1203.65111号 ·doi:10.1007/BF02728986 [67] Pataki，A。；Greengard，L.，圆柱坐标系中的快速椭圆解算器和库仑碰撞算子，J.Compute。物理。，230, 7840-7852, (2011) ·Zbl 1432.31002号 ·doi:10.1016/j.jcp.2011.07.005 [68] 佩妮斯，M。；Walker，HF，NITSOL:非线性系统的牛顿迭代求解器，SIAM J.Sci。计算。，19, 302-318, (1998) ·Zbl 0916.65049号 ·doi:10.1137/S1064827596303843 [69] 波特，GD；爱尔兰共和国。；Boedo，J。；Rognlien，TD，DIII-D刮除层和边缘等离子体中模拟和测量等离子体轮廓的详细比较，Phys。等离子体，73663-3680，（2000）·doi:10.1063/1.1286509 [70] 邱，JM；Christlieb，A.，Vlasov方程的保守高阶半拉格朗日WENO方法，J.Comput。物理。，2291130-1149（2010年）·兹比尔1188.82069 ·doi:10.1016/j.jcp.2009.10.016 [71] 邱，JM；Shu，CW，保正半拉格朗日间断Galerkin公式：理论分析和在Vlasov-Poisson系统中的应用，J.Compute。物理。，230, 8386-8409, (2011) ·Zbl 1273.65147号 ·doi:10.1016/j.jcp.2011.07.018 [72] 明尼苏达州罗森布鲁斯；麦克唐纳，WM；贾德，DL，福克-普朗克反平方力方程，物理学。修订版，107，1-6，（1957）·Zbl 0077.44802号 ·doi:10.1103/PhysRev.107.1 [73] 马萨诸塞州罗斯马尼思；Seal，DC，用于Vlasov-Poisson方程的保正高阶半拉格朗日间断Galerkin格式，J.Compute。物理。，230, 6203-6232, (2011) ·Zbl 1419.76506号 ·doi:10.1016/j.jcp.2011.04.018 [74] Saad，Y.：稀疏线性系统的迭代方法，第2版。费城工业与应用数学学会（2003年）·Zbl 1002.65042号 ·数字对象标识代码：10.1137/1.9780898718003 [75] 萨阿德，Y。；Schultz，MH，GMRES:求解非对称线性系统的广义最小残差算法，SIAM J.Sci。统计计算。，7, 856-869, (1986) ·Zbl 0599.65018号 ·doi:10.1137/0907058 [76] Salmi，S。；Toivanen，J。；Sydow，L.，《带跳跃的随机波动率模型下期权定价的IMEX方案》，SIAM J.Sci。计算。，36，b817-b834，（2014）·Zbl 1308.91193号 ·doi:10.1137/130924905 [77] Scott，B.，边缘剪切层的回旋动力学研究，等离子体物理学。控制。融合，48，a387，（2006）·doi:10.1088/0741-3335/48/5A/S39 [78] 西塔塔诺。；Chacón，L。；Simakov，A.，适用于任意不同质量和温度状态的自适应保守0D-2V多物种Rosenbluth-Fokker-Planck解算器，J.Compute。物理。，318, 391-420, (2016) ·Zbl 1349.76393号 ·doi:10.1016/j.jcp.2016.03.071 [79] 西塔塔诺。；Chacón，L。；西马科夫，A。；Molvig，K.，多维多物种Rosenbluth-Fokker-Planck方程的质量、动量和能量守恒、完全隐式、可扩展算法，J.Comput。物理。，297357-380（2015）·Zbl 1349.65384号 ·doi:10.1016/j.jcp.2015.05.025 [80] 托马斯。；Tzoufras，M。；罗宾逊，A。；Kingham，R。；里奇斯，C。；夏洛克，M。；Bell，A.，《惯性约束聚变等离子体的Vlasov-Fokker-Planck数值模拟综述》，J.Compute。物理。，231, 1051-1079, (2012) ·Zbl 1385.76015号 ·doi:10.1016/j.jcp.2011.09.028 [81] AGR托马斯；金姆，RJ；Ridgers，CP，非线性各向异性不稳定性对等离子体中激光散斑的快速自磁化，新物理学杂志。，11, 033001, (2009) ·doi:10.1088/1367-2630/11/3/033001 [82] 沃尔克，R。；Knoth，O.，应用于大气化学传输建模的隐式显式Runge-Kutta方法，环境。模型。软质。，15, 711-719, (2000) ·doi:10.1016/S1364-8152（00）00034-7 [83] Xu，X。；Bodi，K。；科恩，R。；Krasheninikov，S。；Rognlien，T.，单零托卡马克几何中等离子体输运的TEMPEST模拟，Nucl。融合，50064003，（2010）·doi:10.1088/0029-5515/50/6/064003 [84] Xu，X。；熊，Z。；多尔，M。；Hittinger，J。；Bodi，K。；坎迪，J。；科恩，B。；科恩，R。；科尔拉，P。；科尔贝尔，G。；Krasheninikov，S。；内文斯，W。；秦，H。；Rognlien，T。；斯奈德，P。；Umansky，M.，Edge回转动力学理论和连续模拟，Nucl。融合，47，809，（2007）·doi:10.1088/0029-5515/47/8/1011 [85] XQ Xu；熊，Z。；高，Z。；内文斯，WM；McKee，GR，TEMPEST对边缘等离子体基座测地线声学模式无碰撞阻尼的模拟，Phys。修订稿。，100, 215,001, (2008) ·doi:10.1103/PhysRevLett.100.215001 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。