基于Taylor展开的分步有限元格式.PDF

第 29 卷第 6 期 河 海 大 学 学 报 Vol 29 No. 6 2001 年 11 月 JOURNAL OF HOHAI UNIVERSITY Nov. 2001 基于 Taylor 展开的分步有限元格式 江春波, 杜丽惠, 张庆海 ( ) 清华大学水利水电工程系, 北京 　100084 摘要 :描述了一种基于 Taylor 展开的分步有限元数值格式, 该法在时间上进行分步计算, 在空间上 采用标准的 Galerkin 有限元格式. 对该数值格式的稳定性分析表明, 该法在时间和空间上均具有三 阶精度, 数值稳定性好, 在库朗数 0～1 的范围内均收敛. 相比于 TaylorGalerkin 法, 本有限元法不包 含高阶微分项, 适用于非线性多维问题及具有复杂边界形状的流动. 该法具有计算简便、精度高、数 值稳定性好等优点. 关键词 :分步有限元法; 稳定性分析; 对流扩散方程 ( ) 中图分类号: TV131. 4 　　　文献标识码 :A 　　　文章编号 :1000 1980 2001 06 0099 04 ( ) 有限元法适用于求解具有复杂边界形状的流动. NavierStokes NS 方程和对流扩散方程由于非线性对 流项的存在, 而使空间离散格式成为非线性非对称的形式, 而常规的有限元法实质上是一种中心格式, 当对 流项较强时, 如果有限元网格不适当, 在数值解中就会出现振荡失真. 为消除这种振荡失真, 可以采用加密有 限元网格的途径, 但这势必会增加计算机的内存和计算量, 往往是不经济的. 这时人们就计算格式本身提出 了各种解决方法, 如迎风有限元法. ( 迎风有限元法的基本思想是在加权函数中采用考虑上游效应的不对称权函数, 其中 SUPG Streamline ) [ 1] [2] Upwind/ PetrovGalerkin 法 精度较高; TaylorGalerkin 有限元法 将对流扩散方程分解为对流方程和扩散方 程两部分, 将对流部分在时间展开至二阶, 再结合标准的 Galerkin 有限元空间离散. 这种格式稳定性较好, 且 提高了计算精度, 相比于 SUPG法, 不必修改自由参数. 但 TaylorGalerkin 有限元法在时间上只有二阶精度, 为 了取得稳定性的数值解, 计算时间步长必须较小, 即库朗数 Cr 必须小于 1/ 3 . 另一种再如LaxWendroff 有限 元格式[ 3] , 其出发点是在对流项中以空间导数代替 Taylor 展开式中的时间导数; 而三阶精度的 TaylorGalerkin 有限元法具有更大的稳定域[2] , 但由于高阶微分项的存在而使得数值离散和边界处理比较困难, 特别是在求 解非线性多维问题时, 这种困难尤其显著. 目前还没有见到应用三阶的TaylorGalerkin 法求解 NavierStokes 流 动的文献[4 ,