[論文レビュー] A BGK approximation to scalar conservation laws with discontinuous flux

Florent Berthelin, Julien Vovelle|arXiv (Cornell University)|Mar 22, 2010

Navier-Stokes equation solutions参考文献 11被引用数 5

ひとこと要約

本稿は、空間的不連続なフラックスを有するスカラー保存則のエントロピー解へのBGK緩和モデルの収束を確立する。キネティック定式化と一般化解のための新規還元定理を用いてBGK方程式を分析することで、緩和パラメータ $\varepsilon \to 0$ の際にBGK近似が唯一のエントロピー解に強い $L^p$ 収束することを証明し、不連続フラックス問題における選択の曖昧さを解消する。

ABSTRACT

We study the BGK approximation to first-order scalar conservation laws with a flux which is discontinuous in the space variable. We show that the Cauchy Problem for the BGK approximation is well-posed and that, as the relaxation parameter tends to 0, it converges to the (entropy) solution of the limit problem.

研究の動機と目的

• 不連続フラックス関数を有するスカラー保存則の弱解における一意性の欠如に対処すること。
• このような方程式におけるBGK近似の適切な定式化（well-posedness）を確立すること。
• BGKモデルの緩和極限が元の保存則のエントロピー解に強く収束することを証明すること。
• 一般化キネティック解が古典的エントロピー解に一致することを示す還元定理を導入し、証明すること。
• 不連続フラックス問題におけるBGKに基づく数値スキームの厳密な基礎を提供すること。

提案手法

• フラックスの空間的不連続性を考慮し、不連続係数 $k(x)$ と速度に依存するフラックス関数 $a(\xi)$ を用いて、BGKモデルを緩和キネティック方程式として定式化する。
• 平衡関数 $\chi_{u^\varepsilon}(\xi)$ を導入し、これは区間 $(0, u^\varepsilon)$ の特性関数である。これにより平衡への緩和をモデル化する。
• 文献 [BV06] にインspiredされたキネティック定式化の手法を用い、非平衡効果を表すキネティック測度 $m^\varepsilon$ を導入する。
• エネルギー推定と $L^\infty$ における弱コンパクト性を用いて、BGK方程式の解の存在と一意性を確立する。
• 任意の一般化キネティック解が古典的エントロピー解に還元されることを示す還元定理（定理7）を構築する。
• 還元定理を用いて、$\varepsilon \to 0$ の際に $f^\varepsilon$ が $\chi_u$ に強い $L^p$ 収束することを導出する。ただし、$L^\infty$ における一様有界性を仮定する。

実験結果

リサーチクエスチョン

• RQ1 BGKモデルは、不連続フラックスを有するスカラー保存則に対して適切な近似を提供するか？
• RQ2 BGK緩和極限は、元の保存則のエントロピー解に強く収束するか？
• RQ3 フラックスの不連続性に起因する複数の弱解の中から一意な解を選択する厳密な方法はあるか？
• RQ4 キネティック測度を用いて定義された一般化キネティック解が、古典的エントロピー解に還元可能であるか？
• RQ5 初期データが平衡でない場合でも、BGKモデルの収束は成立するか？

主な発見

• すべての $\varepsilon > 0$ に対してBGK方程式は適切に定式化されており、解は $L^\infty$ に属し、一様有界性を満たす。
• BGK近似は、すべての $1 \leq p < \infty$ に対して $L^p(Q \times \mathbb{R})$ で収束し、$\varepsilon \to 0$ の際にエントロピー解 $\chi_u$ に強く収束する。
• 収束は、任意の一般化キネティック解が古典的エントロピー解に一致することを示す還元定理（定理7）を介して確立される。
• BGKモデルに付随するキネティック測度 $m^\varepsilon$ は弱*-収束して非負の測度 $m$ に収束し、エントロピー条件と整合的であることが保証される。
• 初期データ $f_0^\varepsilon$ が平衡でない場合でも、$L^1$ で弱収束し、$u_0(x) = \int f_0(x,\xi) d\xi$ を満たす初期状態 $f_0$ に収束する限り、結果は成立する。
• 元の [BV06] の証明におけるわずかなギャップを、キネティック定式化における剰余項を明示的に追跡することで是正した。