第三类Dirichlet边界下四阶抛物方程的紧差分格式  

Compact difference schemes for the fourth-order parabolic equations with the third Dirichlet boundary

作　　者：黄钰 高广花[1] HUANG Yu;GAO Guang-hua(School of Science,Nanjing University of Posts and Telecommunications,Nanjing 210023,Jiangsu,China)

机构地区：[1]南京邮电大学理学院,江苏南京210023

出　　处：《山东大学学报（理学版）》2023年第4期16-28,共13页Journal of Shandong University(Natural Science)

基　　金：江苏省自然科学基金面上资助项目(BK20191375);南京邮电大学自然科学孵化基金资助项目(NY220037)。

摘　　要：应用加权平均和高次Hermite插值等技术,提出逼近四阶导数的几个有用的数值微分公式,并对其截断误差进行分析。在此基础上建立求解第三类Dirichlet边界条件下四阶抛物方程初边值问题的三个高阶紧差分格式,应用Fourier分析方法证明格式的无条件稳定性,并对其进行数值验证。这三个差分格式的差异主要体现在空间导数临近边界处的离散方式不同,所得格式全局精度均达到了时间二阶、空间四阶。Based on some techniques involving the weighted average and high order Hermite interpolation,several useful differentiation formulae for approximating the fourth-order derivatives are derived along with the truncation error analyses.Then three high order compact difference schemes are proposed to solve the initial-boundary value problem of the fourth-order parabolic equations with the third Dirichlet boundary conditions.The unconditional stability is proved by the Fourier analysis method.Numerical experiments are carried out.The major difference of the proposed three schemes lies in the different numerical treatment of spatial derivatives near the boundary.The global accuracy of all presented schemes can attain the order of two in time and four in space.

关 键 词：四阶抛物方程 第三类Dirichlet边界 高精度 紧差分格式 稳定性 

分 类 号：O241.82[理学—计算数学]