[論文レビュー] Hodograph Method and Numerical Solution of the Two Hyperbolic Quasilinear Equations System. Part II. Zonal Electrophoresis Equations

E. V. Shiryaeva, M. Yu. Zhukov|arXiv (Cornell University)|Oct 17, 2014

Nonlinear Waves and Solitons参考文献 9被引用数 5

ひとこと要約

本稿では、ゾーナル電気泳動に由来する2階非線形双曲型PDEのコーシー問題を解くための、新しいホドグラフに基づく手法を提示する。リーマン–グリーン関数の分解を用いて、この系をODE問題に還元する。この手法により、数値的粘性や差分スキームを用いずに、近似を伴わず、多値的プロファイルを含む高精度な数値解が得られる。周期的、ガウス型、パルス状初期条件のシミュレーションにより、高い精度で妥当性が検証された。

ABSTRACT

The paper presents the solutions for the zonal electrophoresis equations are obtained by analytical and numerical methods. The method proposed by the authors is used. This method allows to reduce the Cauchy problem for two hyperbolic quasilinear PDE's to the Cauchy problem for ODE's. In some respect, this method is analogous to the method of characteristics for two hyperbolic equations. The method is effectively applicable in all cases when the explicit expression for the Riemann-Green function of some linear second order PDE, resulting from the use of the hodograph method for the original equations, is known. One of the method advantages is the possibility of constructing a multi-valued solutions. Compared with the previous authors paper, in which, in particular, the shallow water equations are studied, here we investigate the case when the Riemann-Green function can be represent as the sum of the terms each of them is a product of two multipliers depended on different variables. The numerical results for zonal electrophoresis equations are presented. For computing the different initial data (periodic, wave packet, the Gaussian distribution) are used.

研究の動機と目的

• ゾーナル電気泳動方程式、すなわち2つの双曲型非線形PDE系を解くための解析的かつ数値的手法を開発すること。
• 従来のホドグラフ法で生じる陰的で超越的な方程式を解く困難さを克服し、問題をODE系に還元すること。
• 衝撃形成や波の破壊の状況において物理的に重要な多値的解の構築を可能にすること。
• 人工的粘性や有限差分離散化を含まない、数値的に安定した近似を伴わない解法フレームワークを提供すること。

提案手法

• ホドグラフ法を用いて、リーマン–グリーン関数の分解を経て、元の双曲型PDE系をODE系に変換する。
• リーマン–グリーン関数が、異なる変数に依存する関数の積の有限和として表現可能であると仮定し、解析的取り扱いを可能にする。
• PDEのコーシー問題は、元のデータから導かれる初期条件を有するホドグラフ変数（a, b）におけるODE系の解に還元される。
• ODEの数値積分には、標準的手法（例：ルンゲ＝クッタ法）を用い、固定ステップサイズ δτ = 0.05 を使用する。
• 等時線 t(a,b) および x(a,b) 上で解を再構築し、特定の時刻における濃度プロファイルの可視化を可能にする。
• この手法は数値的拡散を回避し、リーマンソルバーや有限体積法、人工的粘性を必要としない。

実験結果

リサーチクエスチョン

• RQ1 ゾーナル電気泳動方程式（双曲的かつ非線形）に、ホドグラフ法が効果的に適用可能であり、解析的・数値的解法を可能にするか。
• RQ2 リーマン–グリーン関数の特定構造（可分積の和）をどのように活用して、解法プロセスを簡略化できるか。
• RQ3 周期的、ガウス型、パルス状分布の初期条件に対して、解の挙動はどのように変化するか。
• RQ4 人工的近似や人工的粘性を用いずに、多値的解や衝撃形成を正確に捉えられるか。
• RQ5 元の超越方程式を直接解くのと比較して、削減されたODE系の数値的解法は、効率性と精度の点で優れているか。

主な発見

• 2階非線形双曲型PDEのコーシー問題が、ODE系にうまく還元され、元のPDEを近似せずに数値的に解けることが確認された。
• ガウス型摂動を伴う初期データに対しては、1つの濃度摂動が他方の時間発展に影響を与えることが示され、系の結合性が裏付けられた。
• γ = +1 の場合、時間の経過とともに初期ガウス摂動が2つの明確なピークに分裂し、非線形波干渉の兆候が観察された。
• γ = -1 の場合、t ≈ 19.942 で衝撃波が形成され、多値的濃度プロファイルが出現した。これは、本手法による不連続性の捉え込み能力を示している。
• シミュレーションにおけるグリッド拡散の欠如から、人工的粘性を用いない高精度な解が得られていることが確認された。
• t ≈ 6.358 以降の数値結果には明確な波の破壊と多値的挙動が観察され、本手法の複雑な領域における頑健性が検証された。