[论文解读] Finite element formulation of general boundary conditions for incompressible flows

Roland Becker, Daniela Capatina+2 more|arXiv (Cornell University)|May 6, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 27被引用 7

一句话总结

本论文提出了一种基于Nitsche方法的有限元格式，用于处理不可压缩流中的一般边界条件，该方法在统一的弱形式中同时适用于粘性（Navier-Stokes）和无粘性（Euler）流动 regimes。通过采用通量雅可比矩阵的平衡谱分解，并引入一致的动能稳定项，该方法在所有粘度水平下（包括双曲极限）均保证了稳定性与一致性。数值实验表明，该方法在鲁棒性和尺度不变性方面表现优异，相较于强约束方法，在捕捉复杂流动特征时无振荡现象，表现更优。

ABSTRACT

We study the finite element formulation of general boundary conditions for incompressible flow problems. Distinguishing between the contributions from the inviscid and viscid parts of the equations, we use Nitsche's method to develop a discrete weighted weak formulation valid for all values of the viscosity parameter, including the limit case of the Euler equations. In order to control the discrete kinetic energy, additional consistent terms are introduced. We treat the limit case as a (degenerate) system of hyperbolic equations, using a balanced spectral decomposition of the flux Jacobian matrix, in analogy with compressible flows. Then, following the theory of Friedrich's systems, the natural characteristic boundary condition is generalized to the considered physical boundary conditions. Several numerical experiments, including standard benchmarks for viscous flows as well as inviscid flows are presented.

研究动机与目标

• 开发一种适用于不可压缩流中一般边界条件的统一有限元格式，确保在所有粘度水平（包括无粘性Euler极限）下均保持稳定与一致。
• 解决在高Péclet数和双曲型流动区域中，标准方法失效时离散动能的控制问题。
• 通过通量雅可比矩阵的平衡谱分解，将可压缩流理论中的特征边界条件推广至不可压缩流。
• 确保数值模拟中的鲁棒性与尺度不变性，尤其在复杂流动（如射流冲击与后向台阶流）中表现优异。
• 对比Nitsche方法的弱约束与强约束方法，证明弱约束在边界层与内角区域具有更优的稳定性与更少的振荡。

提出的方法

• 该格式将无粘性（Euler）与粘性（Stokes）贡献分开处理，采用Nitsche方法对所有粘度值（包括µ = 0）实现边界条件的弱约束。
• 通过通量雅可比矩阵的平衡谱分解定义边界矩阵B，确保在退化双曲极限下与内部方程保持量纲一致性。
• 引入额外的一致稳定项以控制离散动能，防止数值解中出现非物理解的能增长。
• 通过将Stokes流的弱格式与基于雅可比矩阵特征值绝对值的一致稳定项相结合，将该方法推广至Navier-Stokes方程。
• 采用连续有限元进行实现，仅关注空间离散化，并与采用SUPG稳定化的强约束变体进行对比。
• 通过在缩放域中求解射流冲击问题并相应变换边界数据，测试尺度不变性，验证该方法在几何与物理尺度变换下保持解的不变性（精度达机器精度）。

实验结果

研究问题

• RQ1 如何使单一有限元格式在粘性与无粘性不可压缩流（包括Euler极限）下均保持鲁棒性与一致性？
• RQ2 动能稳定项在高Péclet数或双曲特征的不可压缩流中，对维持离散能量界起何作用？
• RQ3 通量雅可比矩阵的平衡谱分解如何实现特征边界条件向不可压缩流的自然推广？
• RQ4 与强约束相比，Nitsche方法的弱约束在捕捉涡流与边界层不稳定性等复杂流动特征方面有何数值优势？
• RQ5 所提方法在离散解中在多大程度上保持了尺度不变性，特别是在射流冲击等基准问题中？

主要发现

• 基于Nitsche方法的弱格式成功处理了所有粘度水平，包括无粘性Euler极限，且未损失一致性或稳定性。
• 引入一致的动能稳定项有效抑制了非物理解能量增长，尤其在高Péclet数与双曲区域表现显著。
• 在后向台阶与射流冲击问题的数值实验中，弱约束格式相较于强约束格式产生更少的虚假振荡，尤其在内角与壁面附近表现更优。
• 该方法实现了尺度不变性：当域与边界数据按适当方式缩放时，离散解在机器精度范围内保持不变，而强约束方法则不具备此特性。
• 射流冲击问题中的压力等值线与速度幅值分布图表明，弱格式将振荡局域化于边界附近，而强格式则导致虚假振荡向域内传播。
• 通量雅可比矩阵的谱分解使特征边界条件向不可压缩流的物理解释性推广成为可能，与Friedrichs理论对一阶系统的一致性相吻合。