解非线性对流分散方程的部分迎风有限元法与离散极值原理.pdf

解 部分迎风 原理 王庆晟 (南开大学数学学院，天津，300071夕、， 摘要 本文对二维非线性对流扩散方程建立了～个线性化修正的部分 迎风有限元格式以及一个线性化修正的部 校正格式，证明其分别满足离散极值原理， 式的误差估计． ONPARTIALUPWINDFINITEELEMENTMETHODAND DISCRETEMAXIMUMPRINCIPLEFORNONLINEAR CONVECTl0N—DIFFUSION EQUATIONS Qingsheng Wang SchoolofMathematical Sciences，Nankai abstract Inthis nonlinearconvection．diflusion paper，two．dimensional equations arestudiedAhnearized anda finiteelementscheme pattialupwind hnearized finiteelement axe scheme partialupwind predictor．corrector constructed isshownthatthenumericalsolutionsof respectivelv．It theseschemes discretemaximum The preserve principleconvergence arealso anderrorestimate forbothschemes given o l引言 在流体力学计算中，迎风有限元方法是一类适用于对流占优情形的有限垂 方法．该方法易于实现，适用于多维问题，尤其是它们保留了原问题的两个尊 要性质：极值原理和质量守恒原理，显然，无论从物理或力学方面，还是从数值 分析方面，这两个性质都十分重要．然而，上述迎风格式为了保持解的稳定性， 常常带有过多的人工粘性，从而抹平了解的变化，降低了解的精度．文…一f41 分别以不同方式给出了一个部分迎风有限元格式，降低了人工粘性，进而提高 了数值解的精度． 上面提到的迎风格式几乎都是针对线性对流扩散问题的，没有涉及常见的 非线性方程．文f6】对一类二维非线性Burgers方程给出了部分迎风格式，并给 出了误差估计和收敛性证明．本文针对更一般的二维非线性对流扩散方程，给 出了一类线性化修正的部分迎风格式，证明其满足离散极值原理，并给出了格 式的误差估计；在此基础上，针对上述格式中时间方向上的误差精度较低，又 构造了一类线性化的预测．校正格式，证明其满足离散极值原理，并证明了只 要适当的选取迎风参数，其在时间方向上能达到二阶精度