[論文レビュー] 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

ひとこと要約

本稿では、1次元ヘリーショウフロー（線形移動度を伴う第四階の膜方程式）の完全離散ラグランジュ有限体積スキームを提示する。このスキームは質量保存、非負性保存、およびディリクレエネルギーと対数エントロピーの散逸を保つ。スキームは、時間陰性の定式化により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.

研究の動機と目的

• 連続方程式が持つ主要な物理的・構造的性質を保つ1次元ヘリーショウフロー用の完全離散数値スキームの開発。
• メッシュサイズと時間刻みの極限において、離散解が弱解に収束することを保証するが、CFL条件を必要としないこと。
• 離散設定において、ディリクレエネルギーと対数エントロピーの両方をリャプノフ汎関数として同時に保存すること。
• 特にコンパクトな台を持つ解や初期データがゼロに近い場合に、正値性と質量保存を安定に維持できる強固な数値手法の提供。

提案手法

• ヘリーショウフロー方程式の$L^2$＝ワッサーシュタイン距離における勾配流構造に基づき、方程式の変分的定式化を活用する。
• ラグランジュ有限体積離散化を採用し、粒子が質量点を表し、時間的に変分的時間離散化に従って位置が更新される。
• 離散エネルギーとエントロピーが明示的に時間とともに減少することを示し、安定性と収束解析を可能にする。
• 時間陰性の性質により、CFL型の安定性条件を必要とせず、大きな時間刻みでも安定な時間積分が可能である。
• 離散ワッサーシュタイン距離を用いて変分的時間ステッピングを定義し、連続勾配流構造と幾何的整合性を保証する。
• 収束証明は、離散エネルギーとエントロピー散逸の相互作用に依拠し、コンパクト性の議論と$L^\rho$空間における弱*収束を組み合わせる。

実験結果

リサーチクエスチョン

• RQ1 1次元ヘリーショウフロー用に、完全離散的かつ構造保存型のスキームを構築可能か。特に質量保存と非負性を保てるか。
• RQ2 数値スキームにおいて、対数エントロピー散逸の離散版が成立するか。また、正則性の制御に利用可能か。
• RQ3 CFL条件を必要とせず、離散解が弱解に厳密に収束することを証明可能か。
• RQ4 初期データが小さい値またはゼロ値である場合、標準的な有限差分法と比較して、数値的性能はいかが。

主な発見

• 提案されたスキームは、初期データに適切な仮定を置いた場合、メッシュサイズと時間刻みがゼロに近づく極限において、ヘリーショウフロー方程式の弱解に収束することが証明された。
• 離散設定において、ディリクレエネルギーと対数エントロピーが両方ともリャプノフ汎関数として保存され、離散的散逸率は連続的状態と一致する。
• 時間陰性の定式化によりCFL条件の必要がなく、大きな時間刻みでも安定かつ効率的な時間積分が可能である。
• 数値実験により、初期データがゼロまたは近似的にゼロであっても、正値性と質量保存が維持されていることが確認された。標準的な有限差分スキームでは同様の状況で失敗する。
• 特に質量保存と正値性の高精度が求められる状況では、物理的整合性の観点から、標準的な有限差分法を上回る性能を示した。
• 収束速度は$L^p$ノルムにおいて最適であることが観察され、メッシュの細分化に従い期待通りに誤差が減少した。