求解非定常N-S方程的稳定化分数步长法研究-论文.pdfVIP

第卷 第期 工 程 数 学 学 报 年月 CHINESE JOURNAL OF ENGINEERING MATHEMATICS doi: 10.3969/j.issn.1005-3085.2017.02.005 文章编号: 1005-3085(2017)02-0155-16 求解非定常 方程的稳定化分数步长法研究∗ 李 倩, 贾慧勇 贾宏恩† (太原理工大学数学学院，太原 030024) 摘 要: 本文研究了求解非定常Navier-Stokes 方程的稳定化分数步长法．首先，通过一阶精度 的算子分裂，将非线性项和不可压缩条件分裂到两个不同的子问题中，并对非线性项 采用Oseen 迭代．格式分为两步：第一步求解一个线性椭圆问题；第二步求解一个广 义的Stokes 问题．这两个子问题关于速度都满足齐次Dilichlet 边界条件．同时，在格 式的第二步添加了局部稳定化项，使用等阶序对来加强数值解的稳定性．通过能量估 计方法，对速度与压力做了收敛性分析和误差估计．最后，数值实验验证了方法的有 效性． 关键词: 投影方法；稳定化有限元方法；Navier-Stokes 方程；误差分析 分类号: (2000) 65N30 中图分类号: O241.82 文献标识码: A 引引引言言言 对粘性不可压缩流体的流动进行有限元分析一直是一个热点问题．为了获得稳定收 敛的有限元逼近，通常的速度/压力空间序对需要满足LBB 条件，如Taylor-Hood 元， 但其所形成的总刚度矩阵的维数较大，求解时间较长．低阶协调元(P1-P0, Q1-P0, P1- P1 和Q1-Q1) 序对有着简单的、正则的数据结构，可以节约计算时间，降低计算成本，在 工程实践中有着广泛的应用．然而由于LBB 条件不能被满足，低阶协调元是不稳定的． [1-8] 为了克服这一缺陷，大量的稳定化方法被提出 ．但是，这些方法或者依赖于稳定化参 数的选取，或者涉及到压力的高阶导数，边界积分和网格的相互嵌套等．另一方面，对于 时间相关的微分方程，分数步长方法可以将数值求解所遇到的困难分裂到不同的子问题 上，从而达到降低求解难度的目的[9-13] ． 在文献[14,15] 中，采用分数步长方法对于求解非定常N-S 方程进行了研究，但其对 于三线性项的处理均采用隐式处理．在本文中，我们用分数步长法来求解非定常的N-S 方 程，将非线性项和不可压缩条件分裂到两个不同的子问题中，对于非线性项的处理采用显 式处理．两个子问题关于速度都满足齐次Dilichlet 边界条件，使得中间速度和最终速度在 2 1 (Ω) 均收敛到一个连续解．同时在格式的第二步增加了一个稳定化项， 空间L (Ω) 和H0 使得计算变得稳定、高效．最后，通过数值算例验证了方法的有效性． 收稿日期 作者简介 李倩年月生，女，硕士生 研究方向：偏微分方程数值解与有限元方法 ∗基金项目 国家自然科学基金；山西省自然科学基金 †通讯作者 贾宏恩 万方数据 156 工 程 数 学 学 报 第卷 模模模型型型与与与方方方法法法 Navier-Stokes 方程的基本理论 考虑的模型如下   - νΔ + ( · ∇) + ∇p = , (x, t) ∈ Ω × (0, T], 