一维常系数对流扩散方程的样条子域精细积分法.docx

一维常系数对流扩散方程的样条子域精细积分法林丽烽1,2,刘利斌2,刘桂利2(1.福建农林大学计算机与信息学院,福建福州350002;2.广西民族大学数学与计算机科学学院,广西南宁530006)摘要:基于子域精细积分的理论,提出求解对流扩散方程初边值问题中含参数α0(ατ)样条子域精细积分(SSPI)的方法,其中τ是时间步长;并分析了该方法的稳定性.数值试验结果表明,与三次样条配置法相比,SSPI方法的精度更高,应用也更广.关键词:对流扩散方程;样条子域精细积分(SSPI);三次样条配置;稳定性文章编号:167125470(2008)0420444205中图分类号:O241.8文献标识码:ASplinesubdomainpreciseintegrationschemeforconvection2diffusionequationwithconstantcoefficientLINLi2feng1,2,LIULi2bin2,LIUGui2li2(1.CollegeofComputerandInformation,FujianAgricultureandForestryUniversity,Fuzhou,Fujian350002,China;2.CollegeofMathematicsandComputerScience,GuangxiUniversityforNationalities,Nanning,Guangxi530006,China)Abstract:Bysubdomainpreciseintegrationmethod,thesplinesubdomainpreciseintegration(SSPI)schemecontainingparametersuchasα0(ατ)forinitial2boundaryvalueproblemofconvection2diffusionequationanditsstabilityanalysiswerepresented(τwastimestep).ThenumericalexampleshowedthatSSPImethodhadhigheraccuracyandbetterpracticabilitythanthemethodofcubicsplinecollocation.Keywords:convection2diffusionequation;splinesubdomainpreciseintegration(SSPI);cubicsplinecollocation;stability对流扩散方程可描述众多物理现象,如大气、海洋、湖泊、江河及地下水污染物的迁移扩散过程等,它具有N2S方程的基本特征,故对其进行数值研究具有重要意义.国内外许多学者都对此进行了大量研究,提出了一些高精度的计算格式.Rubinetal[1]使用三次样条近似求解流体力学问题,还将差分法中的交替方向隐式(ADI)法推广到样条中,形成样条交替方向隐式(SADI)法.王璞和Kahawita将求解矩阵系统的过程转化为求解仅包含节点的函数值、一阶导数值、二阶导数值的三类等价的标量方程组[2-5].根据样条近似和配置法的特点,王璞[6]又提出了按坐标方向分步计算Navier2Stokes方程的方法,有效地克服了SA2DI法中每一步都需确定附加边界条件的缺点.王璞又应用三次样条配置法,给出了流体力学方程中的主要的样条函数关系和解算步骤[7].随后,林建国等[8]利用三次样条差分方法对正则化长波方程进行了模拟,取得了较好的结果,进而又给出了含源项非定常非线性对流扩散方程的三次样条四阶差分格式[9].子域精细积分法也得到许多学者的关注[10-13].钟万勰等[10-11]利用子域精细积分方法求解偏微分方程.金承日和刘明珠基于子域精细积分法的理论,针对对流扩散方程和对流扩散方程初值问题,提出了交替分组显式迭代(AGEI)方法[12-13].本文基于子域精细积分理论,提出了一种新的求解对流扩散方程的样条子域精细积分(SSPI)方法.收稿日期:2007-09-18修回日期:2007-12-25基金项目:广西民族大学研究生教育创新基金资助项目(GXUN2CHX0756).作者简介:林丽烽(1980-),女,助教.研究方向:偏微分方程数值解.第4期林丽烽等:一维常系数对流扩散方程的样条子域精细积分法1求解对流扩散方程的SSPI方法1.1SSPI方法的构造考虑如下对流扩散方程的初边值问题:29u=a9u+ε9u9x29t9x(1)u(0,t)=α1(t),u(1,t)u(x,0)=ψ(x)=α2(t)t0x∈[0,1]其中,u(x,t)是x、t的二元未知函数,a,ε是常数,且ε0.分别取时间步长为τ;空间步长为h,h=1