[論文レビュー] General self-similar solutions of diffusion equation and related constructions

Imre Ferenc Barna, László Mátyás|arXiv (Cornell University)|Apr 19, 2021

Fractional Differential Equations Solutions参考文献 25被引用数 23

ひとこと要約

本稿では、群不変性、移動波形プロファイル、自己相似アンザッツの3つの異なる試行関数を用いて、1次元拡散方程式の一般化された自己相似解について包括的な分析を提示する。主な貢献は、$ \alpha \geq 0 $ に対して有効なべき乗則時間減衰 $ C(x,t) \sim t^{-(\alpha+1)} $ を示す正確な解の族であり、$ \alpha = 1, 2, \dots $ の明示的形と、$ \alpha = 0 $ の極限ケース（流体蒸発に関連）を含む。本研究は、拡散系の長時間漸近挙動に関する新たな知見を提供する。

ABSTRACT

Transport phenomena plays an important role in science and technology. In the wide variety of applications both advection and diffusion may appear. Regarding diffusion, for long times, different type of decay rates are possible for different non-equilibrium systems. After summarizing the existing solutions of the regular diffusion equation, we present not so well known solution derived from three different trial functions, as a key point we present a family of solutions for the case of infinite horizon. By this we tried to make a step toward understanding the different long time decays for different diffusive systems.

研究の動機と目的

• 標準的なガウス型とは異なる一般化された自己相似解を導出し、その解析を行うこと。
• 非ガウス的でべき乗則時間減衰を示す拡散系の長時間挙動を調査すること。
• 導出された解を、流体蒸発（$ \alpha = 0 $）や異常拡散といった物理現象と結び付けること。
• 今後、時間的・空間的依存性を示す拡散係数を含む場合の解析を拡張すること。

提案手法

• 群不変性、移動波形プロファイル、自己相似アンザッツの3つの試行関数を用いて正確な解を生成する。
• 自己相似アンザッツ $ C(x,t) = t^{-\alpha} f(\eta) $（$ \eta = x / \sqrt{D t} $）を用い、偏微分方程式を常微分方程式に還元する。
• 整数 $ \alpha $ に対して、退化超幾何関数を用い、多項式補正を加えて $ f(\eta) $ の結果する常微分方程式を解く。
• $ \alpha = 1, 2, \dots $ に対して明示的解を導出し、$ f(\eta) \propto \eta e^{-\eta^2/(4D)} \left( \kappa_0 + \kappa_1 \frac{\eta^2}{4D} + \cdots \right) $ の形をとる。
• 漸近的挙動を分析：$ t $ が大きいとき $ C(x,t) \sim t^{-(\alpha+1)} $ となり、$ t \to t_0 - t $ の代入により爆発解を特定する。
• 解と物理系（流体蒸発、異常輸送など）との関連を確立する。

実験結果

リサーチクエスチョン

• RQ1 標準的なガウス型とは異なる拡散方程式の一般化された自己相似解は何か？
• RQ2 異なる自己相似指数 $ \alpha $ は、濃度分布の長時間減衰挙動にどのように影響を与えるか？
• RQ3 $ \alpha = 0 $ および $ \alpha < 0 $ の極限ケースは、どのような物理現象に対応するか？
• RQ4 多項式補正を含む自己相似アンザッツを体系的に応用して正確な解を導出する方法は何か？
• RQ5 これらの解は、時間的または空間的依存性を示す拡散係数を有する系にどのような意味を持つのか？

主な発見

• $ \alpha = 1 $ の場合、$ D = 1 $ のとき解は $ C(x,t) = \frac{1}{t^2} \eta e^{-\eta^2/(4D)} \kappa_0 \left(1 - \frac{1}{6} \eta^2 \right) $ となり、特定のべき乗則減衰を示す。
• 整数 $ \alpha = n > 2 $ の場合、解は $ C(x,t) = \frac{1}{t^n} \eta e^{-\eta^2/(4D)} \sum_{k=0}^{n-1} \kappa_k \left( \frac{\eta^2}{4D} \right)^k $ の形を取り、ガウス型分布に対する多項式補正が現れる。
• 有限な $ x $ と $ \alpha > 0 $ の場合の漸近的減衰は $ C(x,t) \sim t^{-(\alpha+1)} $ であり、べき乗則緩和が確認される。
• $ \alpha = 0 $ のケースでは、固定された $ x $ において時間に依存しない解が得られ、流体蒸発に関連する極限ケースと解釈される。
• $ \alpha < 0 $ の解は時間が経つにつれて発散するが、$ \alpha \geq 0 $ の解は減衰する。時間反転 $ t \to t_0 - t $ の下で爆発挙動が観察される。
• 本手法により、退化超幾何関数で表せる解が得られ、非線形および時間依存性を示す拡散方程式へも拡張可能である。