[论文解读] Adaptive finite element methods for partial differential equations

Rolf Rannacher|ArXiv.org|May 1, 2003

Advanced Numerical Methods in Computational Mathematics参考文献 5被引用 6

一句话总结

本文提出了一种针对偏微分方程有限元解的后验误差估计方法——双重加权残差（DWR）方法，实现了目标导向的自适应网格加密。通过求解一个对偶问题以加权残差，该方法构建了能够最小化特定输出误差的经济网格——在流体动力学中对阻力计算、阻力最小化及稳定性分析均表现出高效率与高精度。

ABSTRACT

The numerical simulation of complex physical processes requires the use of economical discrete models. This lecture presents a general paradigm of deriving a posteriori error estimates for the Galerkin finite element approximation of nonlinear problems. Employing duality techniques as used in optimal control theory the error in the target quantities is estimated in terms of weighted `primal' and `dual' residuals. On the basis of the resulting local error indicators economical meshes can be constructed which are tailored to the particular goal of the computation. The performance of this {\it Dual Weighted Residual Method} is illustrated for a model situation in computational fluid mechanics: the computation of the drag of a body in a viscous flow, the drag minimization by boundary control and the investigation of the optimal solution's stability.

研究动机与目标

• 开发一种适用于特定输出泛函的伽辽金有限元方法后验误差估计的通用框架。
• 通过最小化目标量在预定精度下的计算成本，实现高效的目标导向网格自适应。
• 将该方法应用于复杂流体动力学问题，包括阻力计算、边界控制下的阻力最小化以及最优流场的稳定性分析。
• 通过自适应加密展示DWR方法在极少自由度下实现高精度的能力。

提出的方法

• 利用最优控制理论中的对偶技术推导后验误差表示，将输出泛函的误差表示为加权原问题残差与对偶问题残差之和。
• 通过求解线性化对偶问题引入双重加权残差（DWR）估计量，其中对偶解为局部误差指标提供权重。
• 将局部误差指标构造为原问题残差（来自离散解）与对偶权重（来自对偶解）的乘积，从而实现局部网格加密。
• 采用分片高阶插值近似对偶解误差，并在单元间平衡误差贡献，以实现最优网格自适应。
• 利用误差平衡策略，优先加密局部误差贡献最大的单元，确保网格上误差分布均匀。
• 在必要时对输出泛函进行正则化，以确保误差估计器的稳定性和收敛性。

实验结果

研究问题

• RQ1 如何为非线性伽辽金有限元近似推导出针对特定输出泛函而非整个解的后验误差估计？
• RQ2 基于双重加权残差的自适应网格加密能否在更少自由度下实现比均匀或启发式加密策略更高的精度？
• RQ3 在边界控制下的阻力计算与阻力最小化中，DWR方法在实际流体动力学问题中的表现如何？
• RQ4 通过边界控制获得的优化流场解的稳定性如何？DWR方法能否支持对此类解的分析？

主要发现

• 在阻力计算中，DWR方法仅用2,244个自由度即实现了阻力系数1%的误差，效率优于均匀加密。
• 在雷诺数Re = 40的阻力最小化问题中，自适应加密在4,264个自由度下达到阻力系数3.16723，而均匀加密需164,928个自由度才能达到3.11800。
• 与均匀加密相比，DWR自适应网格在阻力最小化问题中减少了75%的自由度，同时保持了相近的精度。
• 阻力最小化的最优解位于不稳定边缘，临界特征值的实部趋近于零，表明可能存在霍普夫分岔。
• 流线可视化与特征值分析证实了优化流场中不稳定性的发生，临界特征函数在分岔点附近表现出振荡行为。
• DWR方法实现了对偶解的精确高效计算，从而生成可靠的误差指标，有效指导了所有测试案例中的网格自适应。