[论文解读] Convergence of a fully discrete variational scheme for a thin-film equation

Horst Osberger, Daniel Matthes|arXiv (Cornell University)|Sep 4, 2015

Lattice Boltzmann Simulation Studies被引用 4

一句话总结

本文提出了一种用于一维Hele-Shaw流（具有线性迁移率的四阶薄膜方程）的全离散拉格朗日有限体积格式，该格式保持质量守恒、非负性，并耗散Dirichlet能量与对数熵。该格式在离散到连续的极限下被证明收敛于弱解，其理论依据在于能量与熵耗散之间的相互作用，且由于其时间隐式格式，无需满足CFL条件。

ABSTRACT

This paper is concerned with a rigorous convergence analysis of a fully discrete Lagrangian scheme for the Hele-Shaw flow, which is the fourth order thin-film equation with linear mobility in one space dimension. The discretization is based on the equation's gradient flow structure in the $L^2$-Wasserstein metric. Apart from its Lagrangian character --- which guarantees positivity and mass conservation --- the main feature of our discretization is that it dissipates both the Dirichlet energy and the logarithmic entropy. The interplay between these two dissipations paves the way to proving convergence of the discrete approximations to a weak solution in the discrete-to-continuous limit. Thanks to the time-implicit character of the scheme, no CFL-type condition is needed. Numerical experiments illustrate the practicability of the scheme.

研究动机与目标

• 开发一种用于一维Hele-Shaw流的全离散数值格式，以保持连续方程的关键物理与结构特性。
• 在网格尺寸与时间步长趋于零的极限下，确保离散解收敛于弱解，且无需满足CFL条件。
• 在离散设置中同时保持Dirichlet能量与对数熵作为李雅普诺夫泛函。
• 提供一种鲁棒的数值方法，以保持非负性与质量守恒，尤其适用于具有紧支集或接近零的初始数据的解。

提出的方法

• 该格式基于Hele-Shaw方程在$L^2$-Wasserstein度量下的梯度流结构，利用方程的变分公式。
• 采用拉格朗日有限体积离散化，其中粒子代表质量点，其位置通过变分时间离散化随时间演化。
• 离散能量与熵随时间明确递减，确保稳定性并支持收敛性分析。
• 由于时间隐式特性，该格式无需CFL型稳定性条件，可采用大时间步长。
• 该方法使用离散Wasserstein度量来定义变分时间推进，确保与连续梯度流结构的几何一致性。
• 收敛性证明依赖于离散能量与熵耗散之间的相互作用，结合紧致性论证与$L^\rho$空间中的弱*收敛。

实验结果

研究问题

• RQ1 能否为一维Hele-Shaw方程构造一种全离散、结构保持的格式，以保持质量与非负性？
• RQ2 在数值格式中，对数熵耗散的离散版本是否成立，且能否用于控制解的正则性？
• RQ3 能否在无需CFL条件的情况下，严格证明离散解收敛于弱解？
• RQ4 与标准有限差分方法相比，该格式在数值表现上如何，特别是在初始数据值较小或为零时？

主要发现

• 在初始数据满足适当假设的条件下，所提出的格式在网格尺寸与时间步长趋于零的极限下，被证明收敛于Hele-Shaw方程的弱解。
• 该格式保持了Dirichlet能量与对数熵作为李雅普诺夫泛函，其离散耗散率与连续情形一致。
• 时间隐式格式消除了对CFL条件的需求，使得即使采用大时间步长，也能实现稳定高效的时推进。
• 数值实验表明，该格式在初始数据为零或接近零时仍能保持非负性与质量守恒，而标准有限差分格式在此类情形下会失效。
• 在质量与非负性要求较高的情况下，该格式在物理一致性方面优于标准有限差分方法。
• 该格式在$L^p$-范数下的收敛率被观测为最优，误差随网格细化按预期减小。