×

非线性哈密顿微分方程的高效不变量守恒显式Runge-Kutta格式。 (英语) Zbl 07506164号

摘要:利用不变能量二次化技术,提出了非线性哈密顿常微分方程和偏微分方程的不变守恒显式Runge-Kutta格式的统一框架。首先,将非线性哈密顿微分方程转化为一个允许二次能量不变量的等价形式。对于非线性常微分方程,然后使用一类松弛Runge-Kutta格式离散化重新公式。对于非线性偏微分方程,首先采用傅里叶伪谱离散化方法在空间方向上离散格式,该方法保留了半离散二次守恒定律。然后利用显式松弛Runge-Kutta方法在时间方向上积分半离散系统。该方法可以保持不同的二次守恒定律对机器精度的影响,提高了长时间计算的数值稳定性。此外,该方法在进行时间步长缩放计算时保持了与标准Runge-Kutta格式相同的收敛速度,并且在不进行步长缩放的情况下,收敛速度最多减少了一个。通过对几个常微分方程和偏微分方程的数值实验,说明了所提算法在长时间内的优势,并验证了理论分析。

理学硕士:

65-XX岁 数值分析
35-XX年 偏微分方程

软件:

RRK_rr(R)
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 海尔,E。;卢比奇,C。;Wanner,G.,《几何-数值积分:常微分方程的结构保持算法》,第31卷(2006),Springer·Zbl 1094.65125号
[2] Feng,K。;秦,M.,哈密顿系统的辛差分格式,(哈密顿体系的辛几何算法(2010),施普林格),187-211
[3] 布里奇斯,T.J。;Reich,S.,《多辛积分器:保持辛选择性的哈密顿偏微分方程的数值格式》,Phys。莱特。A、 284、4-5、184-193(2001)·兹比尔0984.37104
[4] Furihata,D。;Matsuo,T.,《离散变分导数方法:偏微分方程的保结构数值方法》(2010),Chapman和Hall/CRC
[5] 秦,M。;王毅,《偏微分方程的保结构方法》(2011),浙江科技出版社
[6] Cooper,G.,轨道问题Runge-Kutta方法的稳定性,IMA J.Numer。分析。,7, 1, 1-13 (1987) ·Zbl 0624.65057号
[7] Sanz-Serna,J.,哈密顿系统的Runge-Kutta格式,BIT-Numer。数学。,28, 4, 877-883 (1988) ·Zbl 0655.70013号
[8] 秦,M。;Zhang,M.,哈密顿系统的辛Runge-Kutta格式,JCM增刊,205-215(1992)·Zbl 0793.65054号
[9] 海尔,E。;Wanner,G.,《求解常微分方程II》(1996),施普林格出版社:施普林格-柏林-海德堡·Zbl 0859.65067号
[10] Morten Dahlby,B.O.,推导偏微分方程保积分数值方法的一般框架,SIAM J.Sci。计算。,33, 5, 2318-2340 (2011) ·Zbl 1246.65240号
[11] Celledoni,E。;格林·V。;麦克拉克伦,R.I。;迈凯轮,D。;奥尼尔,D。;奥雷恩,B。;Quispel,G.,《保存能量》。使用“平均向量场”方法的数值偏微分方程中的耗散,J.Compute。物理。,231, 20, 6770-6789 (2012) ·Zbl 1284.65184号
[12] 蔡伟(Cai,W.)。;李,H。;Wang,Y.,分区平均向量场方法,J.Compute。物理。,370, 25-42 (2018) ·Zbl 1398.65331号
[13] 布鲁尼亚诺。;伊韦纳罗,F。;Trigante,D.,哈密顿边值方法(能量保持离散线积分方法),J.Numer。分析。Ind.申请。数学。,5, 1-2, 17-37 (2010) ·Zbl 1432.65182号
[14] Miyatake,Y.,哈密顿系统的能量保持指数填充连续级Runge-Kutta方法,BIT-Numer。数学。,54, 3, 777-799 (2014) ·兹比尔1304.65263
[15] 王,B。;Wu,X.,保守或耗散系统的指数配置方法,J.Compute。申请。数学。,360, 99-116 (2019) ·Zbl 1418.65093号
[16] Celledoni,E。;麦克拉克伦,R.I。;奥雷恩,B。;Quispel,G.R.W.,《卡汉方法的几何性质》,J.Phys。A、 数学。理论。,46,2,第025201条pp.(2012)·Zbl 1278.65192号
[17] 卡尔沃,M。;拉伯塔,M。;蒙蒂亚诺,J。;Rández,L.,保留Lyapunov函数的投影方法,BIT Numer。数学。,50, 2, 223-241 (2010) ·Zbl 1194.65090号
[18] 海尔,E。;Lubich,C.,梯度系统的能量递减整合,IMA J.Numer。分析。,34, 2, 452-461 (2013) ·Zbl 1321.65115号
[19] Dekker,K.,刚性非线性微分方程Runge-Kutta方法的稳定性,CWI专著,第2卷(1984)·Zbl 0571.65057号
[20] Del Buono,N。;Mastroserio,C.,基于一类四阶段四阶Runge-Kutta方法的保持二次定律的显式方法,J.Compute。申请。数学。,140, 1-2, 231-243 (2002) ·Zbl 1007.65048号
[21] 卡尔沃,M。;Hernández-Abreu,D。;蒙蒂亚诺,J.I。;Rández,L.,《关于用显式Runge-Kutta方法保存不变量》,SIAM J.Sci。计算。,28, 3, 868-885 (2006) ·Zbl 1118.65085号
[22] Ketcheson,D.I.,松弛龙格-库塔方法:内生规范的守恒和稳定性,arXiv预印本·Zbl 1427.65115号
[23] 拉诺查,H。;Sayyari,M。;达尔星。;帕萨尼,M。;Ketcheson,D.I.,《松弛Runge-Kutta方法:可压缩Euler和Navier-Stokes方程的全离散显式熵稳定格式》,SIAM J.Sci。计算。,42、2、A612-A638(2020)·Zbl 1432.76207号
[24] 杨,X。;Ju,L.,二元流体表面活性剂相场模型的线性和无条件能量稳定格式,计算。方法应用。机械。工程,3181005-1029(2017)·Zbl 1439.76029号
[25] 杨,X。;Ju,L.,相场弹性弯曲能量模型的无条件能量稳定性高效线性格式,计算。方法应用。机械。工程,315691-712(2017)·Zbl 1439.74165号
[26] 杨,X。;赵,J。;Wang,Q.,基于不变能量求积方法的分子束外延生长模型的数值近似,J.Compute。物理。,333, 104-127 (2017) ·Zbl 1375.82121号
[27] 杨,X。;Yu,H.,使用不变能量求积方法的相场移动接触线模型的有效二阶无条件稳定格式,SIAM J.Sci。计算。,40、3、B889-B914(2018)·Zbl 1395.65095号
[28] 江,C。;蔡伟(Cai,W.)。;Wang,Y.,基于不变能量求积方法的Sine-Gordon方程的线性隐式局部能量保持格式,J.Sci。计算。,80, 3, 1629-1655 (2019) ·Zbl 1428.65028号
[29] 李,H。;Hong,Q.,洛伦兹力系的一种高效能量保持算法,应用。数学。计算。,358, 161-168 (2019) ·兹比尔1429.65297
[30] 江,C。;Wang,Y。;Gong,Y.,Camassa-Holm方程的任意高阶能量守恒格式,arXiv预印本·Zbl 1439.37079号
[31] 张,H。;钱,X。;Song,S.,非线性哈密顿常微分方程的新型高阶能量守恒对角隐式Runge-Kutta格式,应用。数学。莱特。,102,第106091条pp.(2020)·Zbl 1524.65949号
[32] 陈,C。;Yang,X.,带熔体对流的枝晶凝固相场模型的高效数值格式,J.Compute。物理。,388, 41-62 (2019) ·Zbl 1452.65195号
[33] 杨,X。;Zhang,G.-D.,求解具有一般非线性势的Cahn-Hilliard和Allen-Cahn方程的不变能量积分(IEQ)格式的收敛性分析,J.Sci。计算。,82, 3, 1-28 (2020) ·Zbl 1434.65222号
[34] Cheng,Q。;沈杰。;Yang,X.,使用SAV方法的外延薄膜生长模型的高效和精确的数值方案,J.Sci。计算。,78, 3, 1467-1487 (2019) ·Zbl 1444.76020号
[35] 龚,Y。;赵,J。;Wang,Q.,使用标量辅助变量方法的梯度流模型的任意高阶无条件能量稳定格式,arXiv预印本
[36] 刘,Z。;Li,X.,相场晶体方程的高效修正稳定不变能量求积方法,数值。算法,1-26(2019)
[37] 刘,Z。;Li,X.,基于标量辅助变量的分步求解方案和梯度流的不变能量求积方法,arXiv预印本·兹比尔1481.65195
[38] Cheng,Q.,梯度流的广义标量辅助变量方法(G-SAV),arXiv预印本
[39] Yang,Z。;Dong,S.,基于广义辅助变量的耗散系统离散能量稳定方案路线图,保证正,J.Compute。物理。,404,第109121条pp.(2020)·Zbl 1453.65276号
[40] 格里菲斯,D.F。;Higham,D.J.,《常微分方程的数值方法:初值问题》(2010),Springer Science&Business Media·Zbl 1209.65070号
[41] 卡尔沃,M。;Iserles,A。;Zanna,A.,等谱流的数值解,数学。计算。,66, 220, 1461-1486 (1997) ·Zbl 0907.65067号
[42] Bridges,T.J.,《多符号结构与波传播》,《剑桥哲学学会数学学报》,第121卷,第147-190页(1997),剑桥大学出版社·Zbl 0892.35123号
[43] 赫塞文,J.S。;哥特利布,S。;Gottlieb,D.,《时间相关问题的谱方法》,第21卷(2007),剑桥大学出版社·Zbl 1111.65093号
[44] 张,H。;钱,X。;严,J。;Song,S.,正则对数薛定谔方程的新型高阶质量和能量守恒Runge-Kutta积分器,提交出版,2019
[45] 吉,B。;Liao,H.-l。;龚,Y。;Zhang,L.,时间分数分子束外延生长模型的自适应二阶Crank-Nicolson时间步进方案,arXiv预印本·Zbl 1464.35229号
[46] 江,C。;Wang,Y。;Gong,Y.,一般哈密顿偏微分方程的显式高阶能量保持方法,arXiv预印本·Zbl 1456.65090号
[47] 布鲁格纳诺,L。;伊韦纳罗,F。;Trigante,D.,一种具有能量守恒特性的两步四阶方法,计算。物理学。社区。,183, 9, 1860-1868 (2012) ·Zbl 1305.65238号
[48] 李,H。;Wang,Y。;秦,M.,六阶平均向量场法,J.Compute。数学。,34, 5, 479-498 (2016) ·Zbl 1374.65206号
[49] 朱,H。;陈,Y。;Song,S。;Hu,H.,二维薛定谔方程的辛和多符号小波配置方法,应用。数字。数学。,61, 3, 308-321 (2011) ·Zbl 1207.65120号
[50] 宋,M。;钱,X。;张,H。;Song,S.,非线性薛定谔方程和Korteweg-de-Vries方程的哈密顿边值方法,Adv.Appl。数学。机械。,9, 4, 868-886 (2017) ·Zbl 1488.65518号
[51] 舒,C.-W。;Osher,S.,《本质上非振荡冲击捕获方案的有效实现》,J.Compute。物理。,77, 2, 439-471 (1988) ·Zbl 0653.65072号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。