[論文レビュー] Analytic solutions for the Burgers equation with source terms

Gino I. Montecinos|arXiv (Cornell University)|Mar 29, 2015

Nonlinear Waves and Solitons参考文献 3被引用数 11

ひとこと要約

本稿では、剛性および非線形形を含むソース項を有するバーガース方程式を解くための特徴線に基づく解析的手法を提示する。偏微分方程式（PDE）を、特徴線に沿った進化を支配する等価ODE（特徴線に沿った進化を支配するODE）と、空間時間内での曲線を定義する特徴線ODE（特徴線を定義するODE）の2つの可解な常微分方程式（ODE）に変換することで、ソース項の逆数に原始関数が存在する場合に正確な解が得られる。主な貢献は、数値スキームの検証に役立てるための解析的解を体系的に生成する手順を提供することであり、特に特徴線曲線を求める際に数値的根の探索（バイセクション法）が必要な場合でも有効である。

ABSTRACT

Analytic solutions for Burgers equations with source terms, possibly stiff, represent an important element to assess numerical schemes. Here we present a procedure, based on the characteristic technique to obtain analytic solutions for these equations with smooth initial conditions.

研究の動機と目的

• ソース項を有するバーガース方程式の正確な解析的解を体系的に導出するための手法を開発すること。
• 非線形性および剛性を捉えた正確な解を提供することで、数値スキームの厳密な検証を可能にすること。
• 特徴線に基づくODE系が可解となる条件を確立すること。特に、ソース項およびその逆数に原始関数が存在する条件に焦点を当てる。
• 初期条件に依存する解の暗黙的依存性に対処するため、特徴線曲線の定義から生じる代数方程式を解くためにバイセクション法を導入すること。

提案手法

• 本手法は特徴線法を用いてPDEを2つの常微分方程式（ODE）の系に簡約する。
• PDEに沿った特徴線に沿った進化を支配する等価ODEは、$ \frac{d}{dt}\hat{q}(t) = s(\hat{q}(t)) $ で与えられ、初期条件 $ \hat{q}(0) = h_0(y) $ を満たす。
• 等価ODEが可解であるためには、$ s(q)^{-1} $ が原始関数 $ G(q) $ を持つ必要がある。これにより、解は $ \mathcal{E}(t, h_0(y)) = G^{-1}(t + G(h_0(y))) $ と表現可能である。
• 特徴線ODEは $ \frac{d}{dt}x(t) = \mathcal{E}(t, h_0(y)) $ であり、$ \mathcal{E} $ を積分することで解かれ、$ x = y + \mathcal{F}(t, h_0(y)) $ が得られる。ここで $ \mathcal{F} $ は $ \mathcal{E} $ の原始関数である。
• 完全な解は、$ q(x,t) = \mathcal{E}(t, h_0(y)) $ で定義され、$ y $ は $ x = y + \mathcal{F}(t, h_0(y)) $ を満たす。非明示的ケースでは数値的根の探索が必要となる。
• 非明示的 $ y $ の場合、特徴線曲線の定義から生じる代数方程式に対してバイセクション法が適用される。

実験結果

リサーチクエスチョン

• RQ1 剛性または非線形ソース項を有するバーガース方程式に対して、特徴線に基づくアプローチを用いて解析的解を導出できるか？
• RQ2 ソース項にどのような条件下で等価ODEが可解となるか？
• RQ3 特徴線曲線方程式が暗黙的である場合、初期条件を用いて解をどのように表現できるか？
• RQ4 特徴線曲線定式化において生じる暗黙的代数方程式を解くために有効な数値法は何か？

主な発見

• 等価ODEは、ソース項 $ s(q)^{-1} $ が原始関数 $ G(q) $ を持つ場合に可解であり、解は $ \mathcal{E}(t, h_0(y)) = G^{-1}(t + G(h_0(y))) $ と表現可能である。
• 特徴線ODEは、解 $ \mathcal{E}(t, h_0(y)) $ が原始関数 $ \mathcal{F}(t, h_0(y)) $ を持つ場合に可解であり、曲線は $ x = y + \mathcal{F}(t, h_0(y)) $ と与えられる。
• 線形ソース項 $ s(q) = \beta q $ の場合、解は $ \hat{q}(t) = \frac{h_0(y)}{1 - \beta t h_0(y)} $ であり、特徴線曲線は $ x = y - \frac{\ln(1 - \beta t h_0(y))}{\beta} $ である。
• 2次ソース項 $ s(q) = \beta q^2 $ の場合、解は $ \hat{q}(t) = \frac{h_0(y)}{1 - \beta t h_0(y)} $ であり、特徴線曲線は明示的に可解である。
• 非線形ソース項 $ s(q) = e^{\beta q} $ の場合、解は $ \hat{q}(t) = -\frac{1}{\beta} \ln(e^{-\beta h_0(y)} - \beta t) $ であり、特徴線曲線は複雑な原始関数を用いて表現される。
• 特徴線曲線方程式が $ y $ について明示的に解けない場合、バイセクション法を用いて $ y(x,t) $ を数値的に決定し、任意の $ (x,t) $ における解の評価が可能となる。