[论文解读] Stability of an upwind Petrov Galerkin discretization of convection diffusion equations

Snorre H. Christiansen, Tore Gunnar Halvorsen+1 more|arXiv (Cornell University)|Jun 2, 2014

Differential Equations and Numerical Methods参考文献 3被引用 5

一句话总结

该论文针对对流主导区域（小粘性系数）中的对流-扩散方程，建立了迎风Petrov-Galerkin有限元方法的稳定性。通过使用各向异性的 $ H^{1/2} $ 范数和能量型项，证明了在网格尺寸和粘性系数下均保持一致的连续性与 inf-sup 条件，最多存在对数因子，即使在边界层存在的情况下亦成立。

ABSTRACT

We study a numerical method for convection diffusion equations, in the regime of small viscosity. It can be described as an exponentially fitted conforming Petrov-Galerkin method. We identify norms for which we have both continuity and an inf-sup condition, which are uniform in mesh-width and viscosity, up to a logarithm, as long as the viscosity is smaller than the mesh-width or the crosswind diffusion is smaller than the streamline diffusion. The analysis allows for the formation of a boundary layer.

研究动机与目标

• 建立对流主导区域中对流-扩散方程的Petrov-Galerkin方法在粘性系数趋于零时的统一稳定性。
• 分析在出流边界可能出现边界层形成条件下的方法性能。
• 识别出在这些条件下，连续性与 inf-sup 条件均能一致保持的范数，特别是各向异性的 $ H^{1/2} $ 范数与能量范数。
• 将一维稳定性技术推广至含横向扩散的多维问题。
• 为基于有限元外微分计算与指数拟合的有限元方法提供理论基础，尤其适用于对齐网格。

提出的方法

• 该方法采用符合的Galerkin试函数空间，并通过求解局部伴随问题获得迎风测试函数空间。
• 采用粘性系数依赖的测试函数的Petrov-Galerkin弱形式，以模拟对流主导流场中的指数拟合行为。
• 分析基于[3]的抛物模型，其中时间被视为与流动方向对齐的空间变量。
• 关键工具包括希尔伯特变换和一种变体的Besov空间（$ H^{1/2} $ 的变体），该空间可处理不连续函数并捕捉边界层行为。
• 通过将函数投影到分片常数函数上，以及从分片仿射函数到迎风函数的节点插值，推导出连续性与 inf-sup 估计。
• 在多维情形下，通过引入佩克莱特数条件并控制各向异性扩散项，处理横向扩散的影响。

实验结果

研究问题

• RQ1 具有迎风测试函数的Petrov-Galerkin方法是否能在对流主导区域中实现统一稳定性？
• RQ2 何种范数可确保该方法在粘性系数趋于零时，同时满足连续性与 inf-sup 稳定性？
• RQ3 该方法如何处理出流区域中的边界层，特别是当边界层厚度小于网格尺寸时？
• RQ4 一维稳定性框架在多维对流-扩散问题中可被推广到何种程度？
• RQ5 当存在横向扩散时，佩克莱特数在确保稳定性方面起到何种作用？

主要发现

• 该方法在范数 $\|u\|_{\alpha} = \|u\|_{\mathrm{H}^{1/2}_{00}(0,T)\otimes\mathrm{L}^{2}(V)} + \alpha^{1/2}\|\nabla u\|_{\mathrm{L}^{2}(U)}$ 下，连续性与 inf-sup 条件均保持一致，最多存在对数因子依赖于粘性系数。
• 当佩克莱特数足够大或横向扩散相对于流线扩散足够小时，inf-sup 条件在网格尺寸与粘性系数下保持一致。
• 即使存在边界层，稳定性估计依然稳健，这由指数型边界层轮廓在 $ H^{1/2}_{00} $ 范数下的行为所证实。
• 范数 $\|v\|_{\alpha}^{\prime} = \|v\|_{\mathrm{H}^{1/2}(0,T)\otimes\mathrm{L}^{2}(V)} + \alpha^{1/2}\|\nabla v\|_{\mathrm{L}^{2}(U)}$ 满足 $\|v\|_{\alpha}^{\prime} \preccurlyeq |\log(\alpha)|\|u\|_{\alpha}$，对构造的测试函数成立。
• 双线性形式满足 $\int \beta \dot{u} v + \int \gamma u v + \alpha \int \nabla u \cdot \nabla v \succcurlyeq \|u\|_{\alpha}^2$，确认了 inf-sup 条件最多存在对数因子。
• 分析表明，由于典型边界层解中 $ H^{1/2}_{00} $ 范数的 $ |\log(\alpha)|^{1/2} $ 发散，因此无法实现与 $\alpha$ 无关的 inf-sup 条件。