对流扩散方程的一种显式有限体积——有限元方法精选.pdfVIP

维普资讯 2001年 12月 应用数学与计算数学学报 第 1j卷 第 2期 Dec 2001 COMM ON APPL．MATH AND COM PUT l i．j5 No2 对流扩散方程的一种显式有限体积 一 有限元方法 寞 f安徽大学数学 系 合肥， 230039) 摘 要：本文给 出非线性对流扩散问题的一种有限体积和有限元方法相 结合的显式离散 方法，证 明了数值解 的稳定性，并给出了一个 实际算例 。 关键词：对流扩散方程，有限体积方法，有限元方法，显式方法 ，稳定性 1．引 言 对流扩散方程是描述有粘可压缩流 NavieStokes方程的最简单模型，它 出现在流体力学，气动力学，环境保护 ，化学工程等各种科技领域 ，因而 ， 如何更加准确地求解该方程 的数值解就成为重要 问题 求 以散度形式表示 的偏微分方程的弱解大致有两种常用的方法，即有限元方法 (FE)和 80年 代发展起来 的有限体积方法 (FV)有限元方法对应于能量方法，自然地导出 了椭 圆和抛物 问题即扩散 问题 的单元离散方法 有限体积方法 ，从积分守 恒形式出发，采用单元剖分，选择控制元离散 ，通过引入 以Pdenlann问题近 似解为基础 的数值流通量 ，以一种 自然 的方式反映出流体力学中的质量 ， 动量 ，能量等物理量 的守恒规律 ，可 以很好地处理非线性守恒律 问题 ，如 Euler方程等 ．特别是 80年代 以来 ，由于 自适应 网格 ，结构和无结构网格技 术的发展，有限体积方法得到 了更长足的进步，在处理大变形，激波和各种 复杂流体动力学 问题 的能力，以及在方法 的精度和收敛性的理论研究方面 都有 了实质性的进展 把这两种方法相结合，充分发挥各 自的特点，非线 性对流项利用 FV方法逼近，而扩散项利用有限元方法离散，更好地求解具 有初边值条件 的对流扩散方程．对这种用于对流扩散方程的 FV—FM 格式的 讨论，在 f1,2]中证 明了一种半隐式 FV．FM离散格式的收敛性和误差估计结 果， 3『1给出另一种半隐式 FV—FM格式．本文给 出一种显式 FV—FM 格式， 证 明了格式的稳定性 ，并用一个算例说 明了方法 的有效性 ． 本文 2001年 3月 1日收到 安徽省教育厅 自然科学研究项 目f2001kj012) 维普资讯 应用数学与计算数学学报 15卷 2．对流扩散方程的有限体积法和有限元方法 相结合 的显式离散方法 设 n ∈僻 是多边形有界 区域 ， = (1，x2)∈ ，记 Q =n×0【，T)(0 T。。]、本文考虑如下形式的非线性对流扩散方程 的初边值 问题 娄 ～ () 删 )： 0． 0)= ． z∈Q 其 中，。0为给定的常数，^∈C1【．满足 fdo)=0，9∈clio，T】：W (Q))， ∈ ’f【2)这里 P，g2． 为给出方程的弱形式，记 =硎 (n)，用 ”∈V乘方程的两端，并利用 Green公式，得到 一 娄 +n… 记b(u，)一J蚤2^(u)箍d·((nt·，))如Vu’ z·则问题()的弱解为 满足下列条件的函数 ： a)tI∈L (0，：)nL。。(0 ) b) (“．)+b(“、)十。((，川 = (9， ∈ 以(0．