×

基于FFT的交错边界泊松方程高阶中心差分格式。 (英语) Zbl 1456.65140号

总结:这项工作涉及到快速高阶算法的发展,用于求解具有矩形区域和均匀网格但涉及交错边界的单变量泊松方程。此处交错边界表示边界位于两个相邻网格节点之间的中间。由于交错网格在求解多变量偏微分方程(PDE)的科学计算中的广泛应用,计划中的开发值得进一步研究,但很少在文献中报道,因为在谱方法、紧凑有限差分和快速傅里叶变换(FFT)方面存在巨大挑战处理交错边界的算法。本文引入了一种系统方法来处理这方面的各种开放问题,这是最近开发的用于非交错边界的增广匹配界面和边界(AMIB)方法的自然推广。AMIB方法通过带零加解的浸没边界问题来表示,它将任意高阶中心差分与FFT反演相结合。在交错边界上,所提出的AMIB方法可以处理Dirichlet、Neumann、Robin或任何边界条件组合。对AMIB方法在二维和三维中的四阶、六阶和八阶收敛阶进行了数值验证。此外,所提出的AMIB方法对一些具有挑战性的问题表现良好,例如边界附近的低正则性解、PDE解不满足边界条件以及两端同时包含交错和非交错边界。

MSC公司:

65号06 含偏微分方程边值问题的有限差分方法
65T50型 离散快速傅立叶变换的数值方法
65N85型 含偏微分方程边值问题的虚拟域方法
35J05型 拉普拉斯算子、亥姆霍兹方程(约化波动方程)、泊松方程
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 阿比德,S。;Zeghmati,B.,泊松方程的多重网格缺陷修正和四阶紧致格式,计算。数学。申请。,73, 1433-1444 (2017) ·Zbl 1372.65321号 ·doi:10.1016/j.camwa.2017.01.016
[2] Averbuch,A。;以色列,M。;Vozovoi,L.,矩形区域中任意阶精度的快速泊松解算器,SIAM J.Sci。计算。,19, 933-952 (1998) ·Zbl 0910.35027号 ·doi:10.1137/S1064827595288589
[3] Boisvert,RF,三维亥姆霍兹方程的四阶精确傅里叶方法,ACM Trans。数学。柔和。(TOMS),13,221-234(1987)·Zbl 0626.65104号 ·数字对象标识代码:10.1145/29380.29863
[4] Braverman,E。;以色列,M。;Averbuch,A。;Vozovoi,L.,任意阶精度的快速3D泊松解算器,J.Compute。物理。,144, 109-136 (1998) ·Zbl 1398.65365号 ·doi:10.1006/jcph.1998.6001
[5] Braverman,E。;以色列,M。;Averbuch,A.,三维亥姆霍兹方程的快速光谱求解器,SIAM J.Sci。计算。,2237-2260(1999年)·兹伯利0939.65127 ·doi:10.1137/S1064827598334241
[6] 布鲁格,A。;尼尔森,J。;Kress,W.,《不可压缩Navier-Stokes方程的紧凑高阶有限差分方法》,J.Sci。计算。,17, 551-560 (2002) ·Zbl 1001.76064号 ·doi:10.1023/A:1015166529060
[7] 冯·H。;Long,G。;赵,S.,解决椭圆界面问题的增广匹配界面和边界(MIB)方法,J.Compute。申请。数学。,361, 426-433 (2019) ·Zbl 1422.65311号 ·doi:10.1016/j.cam.2019.05.004
[8] 冯·H。;Zhao,S.,基于FFT的三维泊松方程的高阶中心差分格式,各种边界条件,J.Compute。物理。,410, 109391 (2020) ·Zbl 1436.65159号 ·doi:10.1016/j.jcp.2020.109391
[9] 冯·H。;赵,S.,用FFT加速求解椭圆界面问题的四阶有限差分方法,J.Compute。物理。,419, 109677 (2020) ·Zbl 07507238号 ·doi:10.1016/j.jcp.2020.109677
[10] Fornberg,B.,《有限差分公式中权重的计算》,SIAM Rev.,40,685-691(1998)·Zbl 0914.65010号 ·doi:10.1137/S0036144596322507
[11] Golub,生长激素;黄,LC;西蒙,H。;Tang,W.,不可压缩Navier-Stokes方程有限差分解的快速泊松解算器,SIAM J.Compute。,19, 1606-1624 (1998) ·Zbl 0957.76053号 ·doi:10.1137/S1064827595285299
[12] 古普塔,MM;Kouatchou,J。;Zhang,J.,二阶和四阶离散化多重网格泊松解的比较,J.Compute。物理。,132, 226-232 (1997) ·Zbl 0881.65120号 ·doi:10.1006/jcph.1996.5466
[13] Ge,Y.,三维泊松方程的多重网格法和不等网格四阶紧致差分离散格式,J.Compute。物理。,229, 6381-6391 (2010) ·兹比尔1197.65169 ·doi:10.1016/j.jp.201004.048
[14] Haidvoge博士。;Zang,T.,通过Chebyshev多项式展开精确求解泊松方程,J.Compute。物理。,30, 167-180 (1979) ·Zbl 0397.65077号 ·doi:10.1016/0021-9991(79)90097-4
[15] 北卡罗来纳州坎帕尼斯;Ekaternaris,JA,《不可压缩Navier-Stokes方程的交错网格高精度方法》,J.Compute。物理。,215, 589-613 (2006) ·Zbl 1173.76375号 ·doi:10.1016/j.jcp.2005.11.014
[16] Lai,M-C,极几何上的简单紧凑四阶泊松解算器,J.Compute。物理。,182, 337-345 (2002) ·Zbl 1016.65093号 ·文件编号:10.1006/jcph.2002.7172
[17] 马,ZH;钱,L。;因果关系,DM;顾,HB;Mingham,CG,《不可压缩流的笛卡尔重影单元多重网格泊松解算器》,国际期刊Numer。方法。工程,85,230-246(2011)·Zbl 1217.76052号 ·数字对象标识代码:10.1002/nme.2967
[18] Nagel,J.R.:使用有限差分法(FDM)求解广义泊松方程。犹他大学盐湖城分校电气与计算机工程系(2011年)
[19] 压力,WH;Teukolsky,SA;韦特林,WT;英国石油公司弗兰纳里,《Fortran中的数字配方:科学计算的艺术》(1992),剑桥:剑桥大学出版社,剑桥·Zbl 0778.65002号
[20] 舒曼,美国。;Sweet,R.,在任意大小的交错网格上求解具有neumann边界条件的泊松方程的直接方法,J.Compute。物理。,20, 171-182 (1976) ·doi:10.1016/0021-9991(76)90062-0
[21] 舒曼,美国。;Sweet,R.,《带交错边界条件的泊松方程直接解的快速傅里叶变换》,J.Compute。物理。,75, 123-137 (1988) ·Zbl 0642.65070号 ·doi:10.1016/0021-9991(88)90102-7
[22] 沈杰。;Tang,T。;Wang,LL,《谱方法:算法、分析和应用》(2011),柏林:施普林格出版社,柏林·兹比尔1227.65117 ·doi:10.1007/978-3-540-71041-7
[23] Sun,X.H.,Zhuang,Y.:具有neumann边界条件的亥姆霍兹方程的高阶直接求解器。技术报告。科学与工程计算机应用研究所(ICASE)(1997年)
[24] Swarztrauber,PN,对称FFT,数学。计算。,47, 323-346 (1986) ·Zbl 0615.42008号 ·doi:10.1090/S0025-5718-1986-0842139-3
[25] Swarztrauber,P。;Sweet,R.,算法541:求解可分离椭圆偏微分方程的高效Fortran子程序,ACM Trans。数学。柔和。(TOMS),5352-364(1979)·Zbl 0419.35043号 ·doi:10.1145/355841.355850
[26] 美国特罗滕贝格。;Oosterlee,CW,Multigrid(2001),剑桥:学术出版社,剑桥·Zbl 0976.65106号
[27] Wang,H。;Zhang,Y。;马,X。;邱,J。;Liang,Y.,带Dirichlet边界条件的泊松方程四阶紧致差分格式的有效实现,计算。数学。申请。,71, 1843-1860 (2016) ·Zbl 1443.65288号 ·doi:10.1016/j.camwa.2016.02.022
[28] 王,Y。;张,J.,结合多重网格方法和二维泊松方程外推技术的六阶紧致格式,J.Compute。物理。,228, 137-146 (2009) ·Zbl 1157.65469号 ·doi:10.1016/j.jcp.2008.09.002
[29] Wesseling,P.,《多重网格方法简介》。《纯粹与应用数学》(1992),霍博肯:威利·Zbl 0760.65092号
[30] Zhang,KKO;肖托班,B。;明科维茨,WJ;Mashayek,A.,《Navier-Stokes流交错网格上的紧凑有限差分法》,国际期刊Numer。《液体方法》,52867-881(2006)·Zbl 1370.76112号 ·文件编号:10.1002/fld.1207
[31] 赵,S。;Wei,GW,通过材料界面麦克斯韦方程的导数匹配的高阶FDTD方法,J.Compute。物理。,200, 60-103 (2004) ·Zbl 1050.78018号 ·doi:10.1016/j.jcp.2004.03.008
[32] 赵,S。;魏,GW;Xiang,Y.,用迭代匹配边界法对自由边光束进行DSC分析,J.Sound Vib。,284, 487-493 (2005) ·Zbl 1237.74198号 ·doi:10.1016/j.jsv.2004.08.037
[33] 赵,S.,关于高阶有限差分方法中的伪解,计算。方法应用。机械。工程师,1965031-5046(2007)·兹比尔1173.74462 ·doi:10.1016/j.cma.2007.06.18文件
[34] Zhao,S.,弯曲理想导电边界波导的四阶有限差分方法,计算。方法应用。机械。工程,1992655-2662(2010)·Zbl 1231.78045号 ·doi:10.1016/j.cma.2010.05.005
[35] 赵,S。;Wei,GW,《在高阶中心有限差分中实现边界条件的匹配界面和边界(MIB)》,国际期刊Numer。方法工程,77,1690-1730(2009)·Zbl 1161.65080号 ·doi:10.1002/nme.2473
[36] 周,YC;赵,S。;Feig,M。;Wei,GW,具有间断系数和奇异源的椭圆方程的高阶匹配界面和边界方法,J.Compute。物理。,213, 1-30 (2006) ·Zbl 1089.65117号 ·doi:10.1016/j.jcp.2005.07.022
[37] 庄,Y。;Sun,X.,奇异泊松方程的高阶快速直接求解器,J.Compute。物理。,20, 79-94 (2001) ·Zbl 0990.65120号 ·doi:10.1006/jcph.2001.6771
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。