求解Fisher型方程的混合Jacobi球面调和谱方法.pdfVIP

2009年 12月 应用数学与计算数学学报 第 23卷第 2期 Dec．，2009 COMM ．ON APPL．MATH．AND C0MPUT Vo1．23 No．2 求解 Fisher型方程的混合 Jacobi一球面调和谱方法 刘智翔 黄 伟 摘要本文提出了两种数值求解单位球内Fisher型方程的混合 Jacobi一球面调和谱格 式，并分别给出了格式的收敛性及相关的数值结果． 关键词 混合 Jacobi．球面调和谱方法，Fisher型方程，球形区域 M ixed Jacobi—SphericalHarmonicSpectralM ethod fortheFisher—likeEquation LiuZhixiang HuangW ei Abstract Inthispaper，twomixedJacobi—sphericalharmonicspectralschemesrae proposed for the Fisher—likeequation in a bal1． Theconvergenceofthe schemesand numericalresultsra egivenrespectively． KeywordsmixedJacobi—sphericalhramonicspectralmethod，Fisher—likeequation， sphericalgeometry 1 引 言 对诸如气象科学和海洋科学等领域 中某些 问题的研究往往归结为球形区域上偏 微分方程的数值求解 [13,16】，基于灵活方便的考虑通常会选择有限差分法或有限元素 法．但 由于这类方法的数值精度受到格式本身的限制，在 问题的精确解 比较光滑的情 况下，谱方法更有高精度的优势 【-4，． 在上世纪90年代，郭本瑜等 5[-6,9]发展了以球面调和函数为基函数的逼近方法， 并成功应用于球面上的涡度方程和流体低马赫数流动．此后，郭本瑜等 [8,11-12】建立 了 一 套 比较完整的Jacobi逼近理论，提出了适合于一类奇异问题的Jacobi谱方法，为用 谱方法解决球 内问题提供 了可能．最近，郭本瑜和黄伟 [10,14]发表了有关混合 Jacobi一 球面调和正交逼近 的基本结果 ，并在此基础上解决 了单位球 内用混合谱方法数值求 解 Navier．stokes方程的问题．本文将混合 Jacobi球面调和谱方法应用到单位球 内的 Fisher型方程上，取得 了很好的效果． 收稿 日期： 2009年 9月 10 日． 基金项 目：国家 自然科学基金资助项 目，上海市重点学科建设资助项目(J50101)及上海大学 创新基金资助项 目(A．10—0101—07—408) 1．上海大学数学系，上海 200444；DepartmentofMathematics，ShanghaiUniversity，Shanghai200444， China 18 应用数学与计算数学学报 我们用r， 和 分别表示半径，经度和纬度，并令Q：IxS，其中I= r10≤r1) 以及S={(，)}0≤ ≤27r，一号≤ ≤号)．再令T0，根据 [15]，在单位球内Fisher型 方程可用如下的初边值问题描述 ： f 一AUr(,，8，t)一r(，)、，，t)(1一(7l，)、，，t))：9r(，入，，t)， 2 在 ×(0，]内， (1．1) Iu(r，，，t)：0， 在0fl×[0，上， IU(r，，，0)=Yo(r，，)， 在 上， 其中9(r，，，￡)和 (7’，，)都是已知函数． 2 混合 Jacobi．球面调和逼近 在本节中，我们简要回