摘要:本文是通过优化的KOND算法[1]对二维热传导方程进行数值离散化, 从而得到二维热传导方程的数值解。该算法所用的网格点数少,步长均匀一致。利用网格点之间的对称性,减少了计算机资源的存储量,提高了计算结果的精度与计算过程的速度。该方法与其它方法比较,具有计算精度高、速度快等优点,是求解偏微分方程数值解的比较优秀的方法,也是当前对地震波探究中常采用的方法。
关键词:KOND算法;二维热传导方程;近似解;泰勒展开
本文采用Y.KONDOH的优化近似解析离散化方法“Kernal Optimum Nearly- Analytical discretization Algorithm”(简称KOND算法[8、9、10])求解二维的热传导方程.KOND算法发表于1994年英国大不列颠杂志《计算数学与运用》.经文献查阅,在地震波勘探技术-正演与反演研究中常采用该方法[11], 从杨顶辉教授的研究成果[12]中可看出,KOND方法技巧独特,图像清晰,频散小,数值精度高;在数据处理过程中,减小了数据量存储,提高了计算速度,是方程数值解中较优秀的算法之一.本文就是选择KOND算法求解二维热传导方程的近似解从研究过程了解数学应用的思想,应用的方法技巧,提高数学价值与意义的认识.
在经典的教科书中,热传导系数设为常数a,但在实际问题中,人们往往要确定某时刻,某位置的传导性质,需要设定a为变系数,一般应该设定为a=a(t,x,y),本文仅设为a(x),以方便论述.
由于数学专业的特殊性,可能有很多公式在网页简介里显示不了,在原文中是有的。