[论文解读] Finite Element Technique for Solving the Stream Function Form of a Linearized Navier-Stokes Equations Using Argyris Element

Faisal Fairag, Noah Al-Mulla|ArXiv.org|Jun 3, 2004

Advanced Numerical Methods in Computational Mathematics参考文献 12被引用 3

一句话总结

本文提出一种使用Argyris单元的有限元方法，求解线性化Navier-Stokes方程的流函数形式，重点通过节点重排最小化矩阵带宽，以提高计算效率。该方法实现了最优收敛性和稳定性，数值结果表明在单位正方形区域上使用Re=1和自适应高斯积分规则时，对一个构造的精确解具有高精度的解。

ABSTRACT

The numerical implementation of finite element discretization method for the stream function formulation of a linearized Navier-Stokes equations is considered. Algorithm 1 is applied using Argyris element. Three global orderings of nodes are selected and registered in order to conclude the best banded structure of matrix and a fluid flow calculation is considered to test a problem which has a known solution. Visualization of global node orderings, matrix sparsity patterns and stream function contours are displayed showing the main features of the flow.

研究动机与目标

• 开发线性化Navier-Stokes方程流函数形式的稳定且高效的有限元离散化方法。
• 评估不同全局节点排序对所得刚度矩阵带宽和稀疏结构的影响。
• 实现并测试用于求解有限元离散化产生的非对称、稀疏线性系统的稳定双共轭梯度法（BiCGSTAB）。
• 通过在单位正方形区域上使用已知解析行为的构造精确解，验证数值格式。
• 分析不同网格细化和积分规则下的计算成本、收敛率及迭代求解器性能。

提出的方法

• 流函数形式消除了压力项，并自动满足不可压缩性，将问题简化为求解四阶双调和型PDE。
• 基于Sobolev空间 $ H^2_0(\Omega) $ 推导线性化问题的弱形式，其中双线性形式 $ a(\psi,\varphi) $ 和 $ b(\psi^*;\psi,\varphi) $ 分别表示拉普拉斯项和对流项。
• 采用Argyris单元作为 $ H^2_0(\Omega) $ 中的协调有限元，确保四阶问题所需的 $ C^1 $ 连续性。
• 测试三种全局节点重排策略，以最小化矩阵带宽并改善刚度矩阵的条件数。
• 对所得稀疏、非对称线性系统应用BiCGSTAB迭代求解器，通过残差和迭代值范数容差监控收敛性。
• 数值积分采用四点和六点高斯积分点，以评估计算元素矩阵时的精度和计算成本。

实验结果

研究问题

• RQ1 节点重排如何影响 $ H^2 $-协调离散化中有限元矩阵的带宽和稀疏模式？
• RQ2 积分点选择（4点与6点）对Argyris有限元解的精度和计算成本有何影响？
• RQ3 对于线性化问题，BiCGSTAB迭代次数和CPU时间如何随网格细化（$ h = 1/3, 1/5, 1/9 $）变化？
• RQ4 Argyris有限元方法能否在构造精确解下对流函数形式实现最优收敛率？
• RQ5 初始猜测 $ \psi^* $ 的选择（通过求解双调和方程获得）如何影响迭代求解器的收敛性？

主要发现

• 节点重排1产生的矩阵带宽最窄，CPU时间最低（28.45秒），在计算效率方面优于其他方法。
• 对于 $ h = 1/9 $，$ L^2 $-范数下的相对误差为 $ 1.423 \times 10^{-4} $，$ H^1 $-范数下的误差为 $ 6.986 \times 10^{-4} $，BiCGSTAB迭代次数为567.5次。
• 算术运算量随网格细化显著增加：从 $ h=1/3 $ 时的约1750万次增加到 $ h=1/9 $ 时的约35亿次，反映出 $ H^2 $-协调单元的 $ O(h^{-4}) $ 缩放特性。
• 使用六点积分相比四点积分略微增加了误差和运算量，但改善了稳定性和收敛行为。
• 流线可视化结果在所有网格尺寸下均显示物理解释一致的流场模式，验证了该方法捕捉涡结构的能力。
• 该方法在所有测试的网格尺寸和积分规则下均表现出鲁棒的收敛性，残差和迭代值范数均低于容差 $ 10^{-5} $。