[论文解读] 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

一句话总结

本文通过三种不同的试探函数——群不变性、行波剖面和自相似 Ansatz，对一维扩散方程的一般自相似解进行了全面分析。主要贡献是一个具有幂律时间衰减 $ C(x,t) \sim t^{-(\alpha+1)} $ 的精确解族，适用于 $ \alpha \geq 0 $，包括 $ \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 $）和异常扩散）联系起来。
• 在未来工作中将分析扩展至具有时间与空间依赖扩散系数的情况。

提出的方法

• 采用三种试探函数：群不变性、行波剖面和自相似 Ansatz，以生成精确解。
• 使用自相似 Ansatz $ C(x,t) = t^{-\alpha} f(\eta) $，其中 $ \eta = x / \sqrt{D t} $，将偏微分方程约化为常微分方程。
• 利用合流超几何函数求解得到的 $ f(\eta) $ 的常微分方程，并对整数 $ \alpha $ 使用多项式修正。
• 推导出 $ \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 自相似 Ansatz 如何系统地应用于推导具有多项式修正的精确解？
• RQ5 这些解对具有时间或空间依赖扩散系数的系统有何影响？

主要发现

• 当 $ D = 1 $ 时，对于 $ \alpha = 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 $ 下可观察到爆破行为。
• 该方法可产生用合流超几何函数表示的解，并可推广至非线性和时变扩散方程。