[论文解读] An exponential fitting scheme for general convection-diffusion equations on tetrahedral meshes

Raytcho Lazarov, Ludmil Zikatanov|arXiv (Cornell University)|Nov 5, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 20被引用 21

一句话总结

该论文提出了一种针对四面体网格上对流-扩散-反应方程的新型指数拟合有限元格式，通过将通量结构直接嵌入变分形式，避免了迎风化处理。该方法在解的正则性要求最低的情况下保证了稳定性与一阶收敛性，且在依赖网格的条件下建立了单调性，为具有各向异性扩散的奇异摄动问题提供了鲁棒的替代方案。

ABSTRACT

This paper contains construction and analysis a finite element approximation for convection dominated diffusion problems with full coefficient matrix on general simplicial partitions in $R^d$, $d=2,3$. This construction is quite close to the scheme of Xu and Zikatanov (Math. Comp. 1999) where a diagonal coefficient matrix has been considered. The scheme is of the class of exponentially fitted methods that does not use upwind or checking the flow direction. It is stable for sufficiently small discretization step-size assuming that the boundary value problem for the convection-diffusion equation is uniquely solvable. Further, it is shown that, under certain conditions on the mesh the scheme is monotone. Convergence of first order is derived under minimal smoothness of the solution.

研究动机与目标

• 为非结构化四面体网格上的对流主导型对流-扩散问题开发一种稳定且单调的有限元方法。
• 解决标准伽辽金FEM的局限性，如边界层和内部层中出现的非物理振荡及单调性丧失。
• 构建一种能内在尊重对流-扩散结构的格式，无需依赖迎风化或流场方向检测。
• 在对解和问题数据的正则性假设最弱的条件下，确保一阶收敛性。
• 将指数拟合方法推广至具有完整、可能为各向异性的扩散矩阵以及一般单形剖分的问题。

提出的方法

• 该方法使用一种修改后的双线性形式，将通量 $ J(u) = D\nabla u + \mathbf{b}u $ 直接嵌入变分形式。
• 通过在单元边界上对扩散项和对流项使用调和平均，构造离散问题，以保持一致性和稳定性。
• 该格式源于一种包含基于通量梯度的扰动项的变分形式，从而确保指数拟合行为。
• 通过基于通量重构的方式将对流项的方向性行为嵌入到离散双线性形式中，从而避免迎风化处理。
• 分析依赖于离散的inf-sup条件，以及涉及通量和解的 $ W^{1,p} $-半范数的估计，其中 $ p > d $。
• 引入一种基于缩放函数 $ \alpha(x) $ 的变量变换，以改善离散系统的条件数，尤其适用于病态的扩散矩阵。

实验结果

研究问题

• RQ1 能否为四面体网格上的通用对流-扩散问题构造一种避免迎风化处理的有限元格式，同时保持稳定性和单调性？
• RQ2 何种网格和问题数据条件可确保离散inf-sup条件，从而保证格式的唯一可解性？
• RQ3 与标准伽辽金FEM相比，指数拟合机制在对流主导问题中如何提升稳定性和收敛性？
• RQ4 该方法是否可在对解和数据的正则性假设最弱的条件下实现一阶收敛？
• RQ5 该格式如何适应处理各向异性和病态的扩散矩阵，如多孔介质流中出现的情况？

主要发现

• 在足够小的网格尺寸 $ h $ 下，只要连续问题唯一可解，所提出的格式即保持稳定。
• 在涉及 $ \|\mathbf{b}\|_{1,\infty,T} $ 的依赖网格的条件下，建立了离散格式的单调性，确保了离散最大值原理。
• 在假设 $ u \in W^{1,p}(T) $，$ J(u) \in (W^{1,p}(T))^d $，且 $ \gamma \in C(\overline{T}) $ 的条件下，证明了能量范数下的一阶收敛性，其中 $ p > d $。
• 收敛估计依赖于通量 $ J(u) $ 和解 $ u $ 的 $ W^{1,p} $-半范数，反映了边界层和内部层的影响。
• 当扩散矩阵 $ D(x) $ 病态时，该方法依然保持鲁棒性，尤其在结合合适的缩放函数 $ \alpha(x) $ 以改善条件数时效果更佳。
• 该格式可自然推广至具有矩阵型反应项和非线性对流的系统，因此适用于多组分传输模型。