摘要
针对圆域上的二阶抛物问题,提出了基于高阶多项式逼近的一种有效的数值方法。该方法的主要思想是利用极坐标变换及Fourier基函数展开,将原问题分解为一系列解耦的一维二阶抛物问题。然后,对每个一维二阶抛物问题,建立了一种弱形式及其离散格式,并从理论上证明了该格式的稳定性,弱解和逼近解的存在唯一性以及它们之间的误差估计。最后,给出了一些数值算例,数值结果表明了算法的稳定性和收敛性。
For the second-order parabolic problem in a circular domain,we propose in this paper an effective numerical method based on high-order polynomial approximation.The main idea of this method is to use polar coordinate transformation and Fourier basis function expansion to decompose the original problem into a series of decoupled one-dimensional second-order parabolic problems.Then,for each one-dimensional second-order parabolic problem,a weak form and its discrete scheme is established,and theoretically prove the stability of the schemes,the existence and uniqueness of the weak solution and the approximate solution,as well as the error estimate between them.Finally,some numerical examples is presented,and the numerical results show the stability and convergence of our algorithm.
作者
秦鸿
潘珍兰
安静
QIN Hong;PAN Zhenlan;AN Jing(School of Mathematical Science,Guizhou Normal University,Guiyang,Guizhou 550025,China)
出处
《贵州师范大学学报（自然科学版）》
CAS
北大核心
2024年第2期82-90,共9页
Journal of Guizhou Normal University：Natural Sciences
基金
国家自然科学基金项目(12061023)。
关键词
二阶抛物方程
差分谱逼近
稳定性和误差估计
圆域
Second-order parabolic equation
difference spectral approximation
stability and error estimation
circular domain
