×

能量扩散动力学方程的精确谱数值格式。 (英语) Zbl 1349.76652号

小结:我们研究了在连续动力学计算中使用与非经典权重函数正交的多项式族离散速度变量的优点。我们考虑一个一维偏微分方程模型,该方程描述了福克-普朗克碰撞引起的速度空间中的能量扩散。这种相对简单的情况允许我们将投影动力学的结果与昂贵但高精度的光谱变换方法进行比较。它还允许我们及时准确地进行积分,并完全关注速度变量离散化的有效性。我们表明,对于固定数量的模式或网格点,非经典多项式可以比经典Hermite多项式或等离子体物理中动力学方程的有限差分求解器精确许多数量级。我们详细分析了这两类多项式在行为和精度上的差异。对于非经典多项式,当解释为三维径向函数时,如果初始条件在原点处不光滑,则精确解会离开多项式子空间一段时间,但返回(直到舍入精度)到投影动力学在该时间演化到的同一点。相比之下,使用经典多项式,当精确解返回到子空间时,它与投影动力学解有很大不同。我们还探讨了投影演化算子的本征函数与全演化算子的(非正规)本征函数之间的联系,以及截断计算域的影响。

MSC公司:

76米25 其他数值方法(流体力学)(MSC2010)
65M70型 偏微分方程初值和初边值问题的谱、配置及相关方法
76层25 湍流输送、混合
76X05型 电磁场中的电离气体流动;浆流
82立方31 随机方法(福克-普朗克、朗之万等)应用于含时统计力学问题
82C80码 时间相关统计力学的数值方法(MSC2010)
82D10号 等离子体的统计力学
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 哈泽尔廷,R。;Meiss,J.,《等离子体约束,物理学前沿》(1992),Addison-Wesley:Addison-Whesley Redwood City
[2] 哈泽尔廷,R。;Waelbroeck,F.,《等离子体物理学框架》(1998),英仙座:英仙座雷丁,马萨诸塞州
[3] Schekochihin,A.A。;南卡罗来纳州考利。;多兰,W。;哈米特,G.W。;Howes,G.G。;普朗克,G.G。;Quataert,E。;Tatsuno,T.,《回旋动力学湍流:通过相空间耗散的非线性途径》,《等离子体物理学》。控制。融合,50124024(2008)
[4] 阿贝尔,I。;巴恩斯,M。;科利,S。;多兰,W。;Schekochihin,A.,用于回转运动模拟的线性化模型Fokker-Planck碰撞算子。一、理论、物理学。等离子体,15122509(2008)
[5] 坎迪,J。;荷兰,C。;华尔兹·R。;法希,M。;Belli,E.,《使用直接回转动力学和新古典模拟进行托卡马克剖面预测》,《物理学》。等离子体,16060704(2009)
[6] 巴恩斯,M。;阿贝尔,I。;多兰,W。;Görler,T。;哈米特,G。;Jenko,F.,传输码与回转湍流码的直接多尺度耦合,Phys。等离子体,17056109(2010)
[7] 海兰德,P。;Sigmar,D.,磁化等离子体中的碰撞输运(2002),剑桥大学出版社:剑桥大学出版社·Zbl 1044.76001号
[8] 布拉塔诺夫,V。;Jenko,F。;哈奇,D。;Brunner,S.,离散系统中线性朗道阻尼的方面,物理学。Plasmas,2022108(2013)
[9] Landreman,M。;Ernst,D.,用于连续动力学计算和Fokker-Planck碰撞的新速度空间离散化,J.Comput。物理。,243, 130-150 (2013) ·Zbl 1349.82094号
[10] Shizgal,B.,用于求解Boltzmann方程和相关问题的高斯求积程序,J.Compute。物理。,41, 309-328 (1981) ·Zbl 0459.65096号
[11] Ball,J.S.,半程广义埃尔米特多项式和相关的高斯象限,SIAM J.Numer。分析。,40, 2311-2317 (2003) ·Zbl 1058.41021号
[12] Ghiroldi,G.P。;Gibelli,L.,基于伪谱速度空间离散化的Boltzmann方程直接方法,J.Compute。物理。,258, 568-584 (2014) ·Zbl 1349.76207号
[13] 巴恩斯,M。;阿贝尔,I。;多兰,W。;Ernst,D。;哈米特,G。;里奇,P。;罗杰斯,B。;舍科钦,A。;Tatsuno,T.,用于回转运动模拟的线性化模型Fokker-Planck碰撞算子。二、。数值实现和测试,物理。等离子体,16072107(2009)
[14] Wilkening,J。;Cerfon,A.,奇异Sturm-Liouville问题的谱变换方法及其在等离子体物理能量扩散中的应用,SIAM J.Appl。数学。,75, 2, 350-392 (2015) ·Zbl 1325.65144号
[15] 卡萨姆,A.-K。;Trefethen,L.N.,刚性偏微分方程的四阶时间步进,SIAM J.Sci。计算。,26, 4, 1214-1233 (2005) ·Zbl 1077.65105号
[16] 肯尼迪,C.A。;Carpenter,M.H.,对流-扩散-反应方程的可加Runge-Kutta格式,应用。数字。数学。,44, 1-2, 139-181 (2003) ·Zbl 1013.65103号
[18] 哈奇,D.R。;Terry,P.W。;Jenko,F。;Merz,F。;Nevins,W.M.,通过阻尼本征模的回转湍流饱和,物理学。修订稿。,106, 115003 (2011)
[19] Gautschi,W.,《高斯-克里斯托菲尔求积公式的构造》,数学。计算。,22, 251-270 (1968) ·兹比尔0187.10603
[20] Gautschi,W.,《关于生成正交多项式》,SIAM J.Sci。统计计算。,3, 3, 289-317 (1982) ·Zbl 0482.65011号
[21] 科丁顿,E.A。;莱文森,N.,《常微分方程理论》(1984),克里格出版社:克里格出版社,佛罗里达州马拉巴尔·Zbl 0042.32602号
[22] Stakgold,I.,格林函数和边值问题(1998),威利:威利纽约·Zbl 0897.35001号
[23] Hajmirzaahmad,M。;Krall,A.M.,奇异二阶算子:最大和最小算子,以及介于两者之间的自伴算子,SIAM Rev.,34,4,614-634(1992)·兹伯利0807.34016
[24] 海尔,E。;诺塞特,S.P。;Wanner,G.,《求解常微分方程I:非刚性问题》(2000),Springer:Springer-Blin
[25] Wilkening,J.,计算Jordan链和反演解析矩阵函数的算法,线性代数应用。,427, 6-25 (2007) ·兹比尔1132.47010
[26] Demmel,J.W.,《应用数值线性代数》(1997),SIAM:SIAM Philadelphia·Zbl 0879.65017号
[27] 考克斯·A·J。;Higham,N.J.,加权最小二乘问题的Householder QR因式分解稳定性,(第17届邓迪双年度会议论文集。第17届丹迪双年度大会论文集,《Pitman Res.Notes Math.Ser.》,第380卷(1998),Addison-Wesley Longman:Addison-Wesley Lengman Harlow,英国埃塞克斯),57-73·Zbl 0903.65036号
[28] Fulton,C.,《具有两个奇异端点的二阶Sturm-Liouville问题的Titchmarsh-Weyl(m)-函数》,数学。纳克里斯。,281, 10, 1418-1475 (2008) ·Zbl 1165.34011号
[29] Kotschenreuther,M。;改写,G。;Tang,W.M.,动力学环形等离子体不稳定性的初始值和特征值代码的比较,计算机。物理。社区。,88, 128-140 (1995) ·Zbl 0923.76198号
[30] 坎迪,J。;Waltz,R.E.,《欧拉回转动力学-麦克斯韦解算器》,J.Compute。物理。,186545-581(2003年)·Zbl 1072.82554号
[31] 巴恩斯,M。;多兰,W。;Tatsuno,T.,《解决连续回转动力学中的速度空间动力学》,物理学。等离子体,17032106(2010)
[32] Numata,R。;Howes,G.G。;Tatsuno,T。;巴恩斯,M。;Dorland,W.,AstroGK:天体物理回转动力学代码,J.Compute。物理。,229, 9347-9372 (2010) ·Zbl 1204.85021号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。