对流扩散方程的半拉格朗日有限差分法.pdfVIP

鲁东大学学报 (自然科学版) LudongUniversityJournal(NaturalScienceEdition) 2014，30(2)：114—120 对流扩散方程的半拉格朗 日有限差分法 邵光茹 ，龙晓瀚 (1．大连理工大学 数学科学学院，辽宁 大连 116023；2．鲁东大学 数学与统计科学学院，山东 烟台 264039) 摘要：针对二维对流扩散方程提出了半拉格朗日Crank—Nicolson有限差分格式．化对流扩散方程成拉格朗 日形 式 ，然后进行时间和空间离散．时间的离散采用Crank—Nicolson格式，空间离散采用有限差分格式．数值实验表 明，新建格式是稳定和收敛的． 关键词：对流扩散方程；半拉格朗日法；Crank-Nicolson格式；有限差分法 中图分类号：0241 文献标志码：A 文章编号：1673-8020(2014)02—0114-07 对流扩散方程的数值方法是微分方程数值解法和计算流体力学的研究重点 J．对流扩散方程的 半拉格朗日(semi—Lagrangian，SL)有限差分方法与其它类似的方法，如改进的特征方法、欧拉一拉格朗 日 法 等，通称为特征方法．sL方法与Crank—Nicolson(CN)格式结合来离散对流扩散方程的时间变量， 得到对流扩散方程的SL-CN有限差分法．SL—CN方法与有限元和谱方法结合分别构成SL—CN有限元法 和SL—CN谱方法，这方面已有许多研究成果 J．相对于有限元或谱方法，SL—CN有限差分方法中扩散 项的处理较为复杂．虽然可以预料SL．CN有限差分方法有类似于有限元和谱方法的稳定性和收敛性 ， 但对此还未见数值验证或理论证明．本文针对两个对流扩散问题分别建立SL-CN有限差分格式，并用 算例验证格式的稳定性和收敛性． 1 对流扩散方程的SL方法 1．1 拉格朗 日形式 在处理强对流扩散问题时，欧拉法和拉格朗日法各 自具有局限性．一方面，虽然欧拉法使离散区域 的网格保持不变，但要保证数值方法的稳定性，要求时间步长与空间步长满足很强的约束条件．另一方 面，虽然拉格朗El格式有较好的稳定性，但网格会随着流体的流动而不断发生变化，流体的流动无规则， 使得先前有规则的空间网格变形，甚至发生严重扭曲，这给数值计算带来了很大困难．sL方法介于欧拉 法和拉格朗日法之间，既克服了欧拉法的不稳定，又保持了空间网格不变．因此，sL方法在对流扩散方 程数值求解中具有明显优势． 考虑如下的对流扩散问题： + 。V 一 。△ = n力×( )， (，0)： ，in ， () =0，on02／X(0，T)． 其中，力cR为有界区域，为一个常数，：力×(0，)一R为未知函数，是一个常数；：力×(0，)一 R f~lf：O×(0，) R为已知函数，其中 为流体的速度．假设流体是不可压缩的，3／,满足如下条件： 收稿 日期：2013—11—29；修回日期：2014．01．03 基金项 目：山东省自然科学基金的资助项 目(ZR2010AL021) 作者简介：邵光茹(1989一)，女 ，山东滕州人。硕士研究生，主要从事微分方程理论与数值计算方面的研究。E—mail：2502359000@ qq·COm o 通讯作者：龙晓瀚(1965一)，男，山东烟台人。副教授，硕士研究生导师，主要从事微分方程数值解法研究。E-mail：1ong669@163．corn。 第2期 邵光茹，等 ：对流扩散方程的半拉格朗 日有限差分法 115 u=O,on 3~ , (2) ： 。，in 对于给定的 (，丁)∈ ×(o，)，经过 (，)的特征线对应一个向量函数 (，丁；￡)，满足常微分方程初 值问题 ：