[论文解读] Variational structure of the optimal artificial diffusion method for the advection-diffusion equation

K. B. Nakshatrala, Albert J. Valocchi|arXiv (Cornell University)|May 29, 2009

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

一句话总结

本文通过从变分原理推导出稳定弱格式，为对流-扩散问题中的最优人工扩散方法建立了变分基础。该方法在常系数和常数源项的一维问题中产生精确的节点解，证明了最优人工扩散方法等价于基于指数权函数的加权内积的变分格式。

ABSTRACT

In this research note we provide a variational basis for the optimal artificial diffusion method, which has been a cornerstone in developing many stabilized methods. The optimal artificial diffusion method produces exact nodal solutions when applied to one-dimensional problems with constant coefficients and forcing function. We first present a variational principle for a multi-dimensional advective-diffusive system, and then derive a new stable weak formulation. When applied to one-dimensional problems with constant coefficients and forcing function, this resulting weak formulation will be equivalent to the optimal artificial diffusion method. We present representative numerical results to corroborate our theoretical findings.

研究动机与目标

• 为最优人工扩散方法提供变分基础，该方法此前缺乏变分表述。
• 解决经典伽辽金方法在对流主导问题中出现的不稳定性和振荡问题。
• 利用变分原理推导出对流-扩散方程的稳定弱格式。
• 证明所推导的格式在一维情况下与最优人工扩散方法产生相同的差分方程。
• 强调变分原理在非自伴算子（如对流-扩散算子）中的理论重要性。

提出的方法

• 将对流-扩散方程表述为非自伴算子，并利用带权函数 α(x) = exp(−vx/k) 的加权内积推导出变分原理。
• 应用Vainberg定理，将加权残差表述与标量泛函联系起来，从而实现变分格式。
• 通过最小化涉及加权内积的泛函，推导出稳定弱格式，确保对流主导问题的稳定性。
• 使用基于三节点带宽的线性有限元形状函数，推导出离散差分方程。
• 引入依赖于佩克莱特数的权函数 α(x)，确保一维情况下的稳定性和节点精确性。
• 将所得离散系统与经典最优人工扩散方法进行比较，确认在一维情况下具有等价性。

实验结果

研究问题

• RQ1 能否从变分原理推导出最优人工扩散方法？
• RQ2 是否存在一种对流-扩散方程的稳定弱格式，能在常系数的一维问题中产生精确的节点解？
• RQ3 能否利用Vainberg定理为对流-扩散算子构造一个标量泛函？
• RQ4 所推导的变分格式是否在一维情况下产生与最优人工扩散方法相同的离散系统？
• RQ5 在非自伴问题中，使用非均匀权函数的变分格式具有何种理论意义？

主要发现

• 从变分原理导出的稳定弱格式，对于常系数和常数源项的一维对流-扩散问题，与最优人工扩散方法产生相同的差分方程。
• 所得离散系统中的系数与最优人工扩散方法完全一致：γ−1 = −v/(2h)(1 + coth(Ph/e))，γ0 = v/h coth(Ph/e)，γ+1 = v/(2h)(1 − coth(Ph/e))。
• 数值结果证实，新格式在一维情况下对所有佩克莱特数均实现节点精确解，验证了理论等价性。
• 通过推导出的泛函和权函数 α(x) = exp(−vx/k)，证明了最优人工扩散方法具有坚实的变分基础。
• 该方法在理论上具有鲁棒性，并为开发新型稳定格式提供了基础，尽管在大规模二维和三维问题中可能计算成本较高。
• 本研究强调了Tonti关于非自伴系统变分原理的理论工作，在稳定有限元方法背景下的相关性。