[论文解读] Numerical Methods for Solving Convection-Diffusion Problems

Alexander Churbanov, Petr N. Vabishchevich|arXiv (Cornell University)|Aug 28, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 36被引用 79

一句话总结

本文提出了一类无条件稳定、单调的有限差分与有限元格式，用于对流-扩散方程，采用对流输运的守恒型、非守恒型及反对称形式。该研究在希尔伯特空间与巴拿赫空间中建立了稳定性与收敛性，提出新颖的局部一维（LOD）与加法平均分裂格式，确保在适当扩散与速度场条件下，任意时间步长下均保持单调性与最优误差界。

ABSTRACT

Convection-diffusion equations provide the basis for describing heat and mass transfer phenomena as well as processes of continuum mechanics. To handle flows in porous media, the fundamental issue is to model correctly the convective transport of individual phases. Moreover, for compressible media, the pressure equation itself is just a time-dependent convection-diffusion equation. For different problems, a convection-diffusion equation may be be written in various forms. The most popular formulation of convective transport employs the divergent (conservative) form. In some cases, the nondivergent (characteristic) form seems to be preferable. The so-called skew-symmetric form of convective transport operators that is the half-sum of the operators in the divergent and nondivergent forms is of great interest in some applications. Here we discuss the basic classes of discretization in space: finite difference schemes on rectangular grids, approximations on general polyhedra (the finite volume method), and finite element procedures. The key properties of discrete operators are studied for convective and diffusive transport. We emphasize the problems of constructing approximations for convection and diffusion operators that satisfy the maximum principle at the discrete level --- they are called monotone approximations. Two- and three-level schemes are investigated for transient problems. Unconditionally stable explicit-implicit schemes are developed for convection-diffusion problems. Stability conditions are obtained both in finite-dimensional Hilbert spaces and in Banach spaces depending on the form in which the convection-diffusion equation is written.

研究动机与目标

• 为流体动力学与传热/传质问题中出现的对流-扩散方程开发鲁棒的、无条件稳定的数值方法。
• 通过对方流与扩散算子的单调逼近，确保离散最大值原理的满足。
• 构建具有最优稳定性与收敛性质的两层与三层时间格式。
• 在有限维希尔伯特空间与巴拿赫空间中，针对对流-扩散方程的不同形式，分析稳定性与单调性。
• 设计高效、可并行化的分裂格式——局部一维与加法平均分裂格式，以利于实际实现。

提出的方法

• 在矩形网格上使用有限差分格式，在一般多面体网格上使用有限体积方法，对空间离散化采用兼容的有限元方法。
• 采用对流输运算子的守恒型（散度型）、非守恒型（特征型）及反对称形式，以模拟不同的物理行为。
• 应用加权时间推进格式（如Crank-Nicolson型），引入参数σ，以实现二阶精度与无条件稳定性。
• 通过演化算子的乘法分裂引入局部一维（LOD）格式，确保在时间步长约束下保持单调性。
• 提出基于加法分裂的加法平均LOD格式，实现各分量的独立计算，提升并行性。
• 通过矩阵分析推导稳定性与单调性条件，要求迭代矩阵为M-矩阵且其范数不超过1。

实验结果

研究问题

• RQ1 能否为对流-扩散方程在守恒与非守恒形式下构造无条件稳定且单调的数值格式？
• RQ2 对流输运算子的不同形式（散度型、非散度型、反对称型）如何影响数值解的稳定性和精度？
• RQ3 在有限维希尔伯特空间与巴拿赫空间中，时间分裂格式的单调性与稳定性的充要条件是什么？
• RQ4 局部一维或加法平均格式能否在任意时间步长下保持单调性与最优误差界？
• RQ5 与标准方法相比，所提格式在可并行性与计算效率方面表现如何？

主要发现

• 采用乘法分裂的局部一维格式在σ = 1时，对任意τ > 0均保持单调性，并满足L∞或L1范数下的先验估计；当σ < 1时，在时间步长约束下亦保持单调性。
• 加法平均LOD格式在类似条件下确保单调性与相同的先验估计，且由于分量可独立计算，显著提升了并行性。
• 通过验证迭代矩阵S₁与S₂为M-矩阵且其范数≤1，证明了稳定性和单调性，从而确保满足离散最大值原理。
• 当σ = 1时，格式在时间方向达到一阶精度且无条件稳定；当σ < 1时，稳定性需满足涉及扩散与速度项加权和最大值的时间步长限制。
• 该格式适用于对流-扩散方程的多种形式，包括热传导方程、浓度输运方程以及类似纳维-斯托克斯的动量方程。
• 理论分析通过L∞与L1范数下的先验估计得到验证，证实了在问题数据一般假设下的收敛性与鲁棒性。