[論文レビュー] Numerical modelling of convection-diffusion-reaction problems with free boundary in 1D

Gabriela Kacurova|ArXiv.org|Sep 2, 2009

Fluid Dynamics and Thin Films参考文献 7被引用数 3

ひとこと要約

本稿では、自由境界を伴う1次元の対流拡散反応問題に対して、局所的な移動波解の形状から導かれたODEを用いて界面の時間発展をモデル化する、新しい数値解法を提示する。問題を固定領域に変換し、界面付近での適応的空間離散化を施すことにより、鋭い波面を高精度に捉えることが可能となり、多孔質媒体や反応拡散方程式の既存の解析解および数値研究と比較することで検証された。

ABSTRACT

We discuss a numerical method for convection-diffusion-reaction problems with a free boundary in 1D. The method is based on the numerical modelling of the interface evolution, the transformation to a fixed domain problem and the approximation by an ODE system. The interface evolution is modelled by means of the local shape of the corresponding travelling wave solution. The method can be applied to many free boundary problems with a finite speed of the interface. The presented method can also approximate some problems with an infinite speed of the interface for damped travelling wave type solutions. In the numerical experiments we compare our numerical solution with the analytical ones for some problems.

研究の動機と目的

• 自由境界を伴う1次元の対流拡散反応問題に対して、特に有限速度で移動する界面を有する場合に適した効率的かつ高精度な数値解法の開発。
• 局所的な移動波解から導かれたODEを用いて、界面 $ s(t) $ の時間発展をモデル化し、高価な適応メッシュ生成を回避する。
• 変数変換 $ y = x/s(t) $ を用いて移動境界問題を固定領域問題に変換し、安定的かつ高精度な空間離散化を可能にする。
• 非線形拡散、多孔質媒体、反応拡散方程式に本手法を適用し、非退化および減衰する移動波解を有する場合を含む。
• 多孔質媒体方程式や乱流流れモデルなどのベンチマーク問題に対して、解析解および既存の数値結果と比較して手法を検証する。

提案手法

• 多孔質媒体型方程式を想定し、局所的な移動波解の形状から導かれたODE $ \dot{s}(t) = -\frac{n}{n-1} \partial_x u^{n-1} \big|_{x \nearrow s(t)} $ を用いて界面発展 $ s(t) $ をモデル化する。
• 移動領域 $ (0, s(t)) $ を変数変換 $ y = x/s(t) $ により固定領域 $ (0, 1) $ に変換し、標準的な空間離散化を可能にする。
• 特に $ y = 1 $（すなわち $ x = s(t) $）付近に高い解像度を持つ特別な空間離散化を適用し、鋭い波面を高精度に捉える。
• 変換後のPDEを時間方向に離散化し、界面ODEを含むODE系に変換し、剛性のあるODE系として扱い、安定な数値解法で解けるようにする。
• 非退化ケース（例：線形分配等温線）を近似するために、非線形性を $ p \nearrow 1 $ に近づけることで摂動を加え（例：$ \Psi(s^p) $）、界面の形成を回避しつつも精度を維持する。
• 数値解および界面を、解析解（例：Barenblatt-Pattle解）および既に発表済みの数値ベンチマークと比較することで、結果の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 局所的な移動波解の形状から導かれたODEを用いることで、1次元の対流拡散反応問題における界面発展を高精度にモデル化できるか？
• RQ2 変数変換 $ y = x/s(t) $ を用いて移動境界問題を固定領域に変換することで、界面付近での適応的空間解像度を有する安定かつ高精度な数値解法が可能になるか？
• RQ3 本手法は、非退化（無限速度、減衰波を伴う）および非線形性に依存する界面を有するケース（例：分配等温線）を、両方とも高精度に近似できるか？
• RQ4 多孔質媒体方程式や乱流流れモデルなどのベンチマーク問題において、既存の手法と比較して精度および効率に優れているか？
• RQ5 非線形性を $ p \nearrow 1 $ に近づけることで非退化問題を近似可能であり、その近似はどの程度安定か？

主な発見

• 多孔質媒体系（$ n=6 $）において、数値解と解析解の界面がほぼ完全に一致しており、鋭い波面を高精度に捉えていることが裏付けられた。
• 多孔質媒体方程式において、数値的界面 $ s(t) $ が解析解 $ s(t) = [\frac{2n(n+1)}{n-1}(t+1)]^{1/(n+1)} $ とほとんど誤差なく一致し、界面ODEモデルの妥当性が確認された。
• 非退化拡散（$ p=0.95 $）のケースでは、線形分配を安定かつ高精度に近似でき、[10]および[24]の結果とよく一致した。
• 乱流流れモデル $ \partial_t u = \partial_x^2 u^{3/2} - (u^{3/2} - u^{1/2}) $ において、数値的界面と解析的界面は視覚的に区別できず、時間経過に伴っても誤差は低く保たれた。
• 泡の排水モデル $ \partial_t u = \partial_x^2 u^{3/2} + \partial_x(u^2) $ および粘性液体の前進モデル $ \partial_t u = \partial_x^2 u^4 + \partial_x(u^3) $ においても、期待される物理的挙動（対称的拡散、波面の進行）を正確に再現した。
• 文献[24]のオペレータ分割法の結果と比較したところ、$ D=0.1 $ および $ D=0.001 $ の両ケースで極めて良好な一致が得られ、異なる拡散係数に対しても本手法の頑健性が確認された。