[论文解读] Numerical computations of split Bregman method for fourth order total variation flow

Yoshikazu Giga, Yuki Ueda|arXiv (Cornell University)|Jun 11, 2019

Fluid Dynamics and Turbulent Flows参考文献 51被引用 10

一句话总结

本文提出了一种新颖的分裂Bregman框架，用于通过后向欧拉时间离散化和分段常数空间离散化求解四阶总变差流和Spohn的四阶模型。该方法为Spohn模型引入了新的软阈值算子，并在周期性边界条件下展示了有效的数值解法，实现了2D模拟中准确的平坦晶面和表面松弛效果，参数缩放满足O(h⁻⁴)和O(h⁻²)。

ABSTRACT

The split Bregman framework for Osher-Sol\'e-Vese (OSV) model and fourth order total variation flow are studied. We discretize the problem by piecewise constant function and compute $ abla(-\Delta_{\mathrm{av}})^{-1}$ approximately and exactly. Furthermore, we provide a new shrinkage operator for Spohn's fourth order model. Numerical experiments are demonstrated for fourth order problems under periodic boundary condition.

研究动机与目标

• 开发一种高效的数值格式，用于求解四阶总变差流和Spohn的四阶模型，这些模型具有高度奇异性且计算困难。
• 将已在二阶问题中得到充分验证的分裂Bregman框架扩展至四阶区域，采用后向欧拉时间推进方法。
• 为Spohn模型设计一种新型软阈值算子，以实现晶体表面松弛的稳定且精确的数值求解。
• 提供一种鲁棒的数值方法，能够保持关键物理特征（如平坦晶面和不连续性），与理论预期一致。

提出的方法

• 在1D或2D环面上使用分段常数有限元方法对问题进行离散化，以处理四阶PDE中的奇异性。
• 对次微分形式应用后向欧拉方法，将PDE转化为包含总变差和H⁻¹保真项的最小化问题。
• 利用分裂Bregman框架解耦L1和L2部分，实现高效的交替最小化。
• 实现两种方案：一种使用∇(−Δₐᵥ)⁻¹的近似，另一种使用通过二阶B样条实现的精确逆拉普拉斯算子，适用于分段常数函数。
• 通过在分裂Bregman框架下求解欧拉-拉格朗日方程，为Spohn模型推导出新的软阈值算子，同时结合ℓ¹和ℓᵖ范数。
• 采用一致的参数缩放：在2D中λ = O(h⁻⁴)，µ = O(h⁻²)，确保数值实验中的稳定性和收敛性。

实验结果

研究问题

• RQ1 分裂Bregman框架能否成功扩展至四阶总变差流和Spohn模型，这些模型比二阶对应物更具奇异性？
• RQ2 在四阶PDE背景下，如何高效且准确地计算分段常数函数的逆拉普拉斯算子？
• RQ3 为处理Spohn四阶模型中ℓ¹与ℓᵖ结构的组合，需要哪些新型软阈值算子？它们如何提升数值稳定性？
• RQ4 数值解是否能再现理论预测的物理特征，如平坦晶面和不连续性？

主要发现

• 所提出的分裂Bregman格式在周期性边界条件下成功计算了四阶总变差流，生成的解具有平坦晶面和不连续性，与理论预期一致。
• 各向同性和各向异性四阶总变差流的数值结果表现出稳定的演化过程，具有锐利的过渡区，且各向异性情形保持了分段常数结构，支持理论猜想。
• 为Spohn模型设计的新软阈值算子有效实现了晶体表面松弛的模拟，生成了晶面形成和光滑演化过程，与先前研究结果一致。
• 在2D中采用λ = 1.25h⁻⁴和µ = 5h⁻²时，该格式表现出良好的性能，证明了其在复杂四阶PDE上的鲁棒性和收敛性。
• 使用二阶B样条实现精确逆拉普拉斯算子的方法，相较于近似方法显著提升了精度，尤其在捕捉非局部效应方面表现更优。