[论文解读] Relaxation Equations: Fractional Models

Ester C. F. A. Rosa, E. Capelas de Oliveira|arXiv (Cornell University)|Oct 6, 2015

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

一句话总结

本文提出了一种基于 Riemann-Liouville 分数阶导数（阶数为 $\gamma$，$0 < \gamma \leq 1$）的广义分数阶动力学模型，用于描述介电松弛行为。通过拉普拉斯变换求解所得方程，得到以 Mittag-Leffler 函数表示的解。其主要贡献在于：当 $\gamma = 1$ 时，经典 Debye、Cole-Cole、Cole-Davidson 和 Havriliak-Negami 松弛函数均成为该模型的特例，从而在分数阶框架下统一了这些模型，增强了对异常松弛行为的灵活性描述。

ABSTRACT

The relaxation functions introduced empirically by Debye, Cole-Cole, Cole-Davidson and Havriliak-Negami are, each of them, solutions to their respective kinetic equations. In this work, we propose a generalization of such equations by introducing a fractional differential operator written in terms of the Riemann-Liouville fractional derivative of order $γ$, $0 &lt; γ\leq 1$. In order to solve the generalized equations, the Laplace transform methodology is introduced and the corresponding solutions are then presented, in terms of Mittag-Leffler functions. In the case in which the derivative's order is $γ=1$, the traditional relaxation functions are recovered. Finally, we presente some 2D graphs of these function.

研究动机与目标

• 通过在动力学方程中引入阶数为 $\\gamma$（$0 < \\gamma \leq 1$）的分数阶导数，推广经典介电松弛模型。
• 利用拉普拉斯变换方法，推导这些分数阶动力学方程的解析解。
• 将解表示为 Mittag-Leffler 函数，以推广经典松弛函数。
• 证明当 $\gamma = 1$ 时，经典松弛模型（Debye、Cole-Cole、Cole-Davidson、Havriliak-Negami）可作为特例被恢复。
• 生成不同参数取值下分数阶松弛函数的二维图形表示。

提出的方法

• 使用 Riemann-Liouville 分数阶导数（阶数为 $\gamma$）构建分数阶动力学方程，以建模松弛过程。
• 基于对应于每个经典松弛模型（Debye、Cole-Cole、Cole-Davidson、Havriliak-Negami）的积分核构造记忆函数。
• 对分数阶动力学方程应用拉普拉斯变换，推导变换后松弛函数的代数表达式。
• 利用已知的 Mittag-Leffler 函数的拉普拉斯变换恒等式，特别是 $\mathscr{L}[t^{\\gamma-1}E_{\\alpha,\\gamma}^{\\beta}(at^{\\alpha})](s) = \frac{s^{\\alpha\beta - \\gamma}}{(s^{\\alpha} - a)^{\\beta}}$，实现解的逆变换。
• 通过设定 $\gamma = 1$ 及特定参数值（如 $\beta = 1$，$\alpha = 1$）恢复经典模型。
• 生成分数阶松弛函数的二维图像，以可视化其在 $\gamma$、$\alpha$、$\beta$ 和 $\sigma$ 参数变化下的行为。

实验结果

研究问题

• RQ1 分数阶微分方程（使用 Riemann-Liouville 导数，阶数 $\gamma \in (0,1]$）能否推广经典介电松弛模型？
• RQ2 当 $\gamma = 1$ 时，这些分数阶方程的解是否能还原为已知的经典松弛函数（Debye、Cole-Cole、Cole-Davidson、Havriliak-Negami）？
• RQ3 当 $\gamma < 1$ 时，基于 Mittag-Leffler 函数的解表现出何种行为？其单调性如何？是否具有物理解释性？
• RQ4 参数 $\alpha$、$\beta$、$\gamma$ 和 $\sigma$ 在分数阶框架中对松弛动力学的塑造作用是什么？
• RQ5 能否通过拉普拉斯变换关系 $\tilde{\\varepsilon}(s) = 1 - s^{\\gamma}\tilde{\\varphi}(s)$ 一致地从分数阶松弛函数推导出复介电常数？

主要发现

• 使用 Riemann-Liouville 导数（阶数为 $\gamma$）的分数阶动力学方程可推广经典松弛模型，其解以 Mittag-Leffler 函数表示。
• 当 $\gamma = 1$ 时，分数阶 Debye 模型恢复经典 Debye 松弛函数 $\varphi(t) = e^{-t/\sigma}$。
• 分数阶 Cole-Cole 模型给出 $\varphi(t)_{CCF} = \frac{t^{\gamma-1}}{\Gamma(\gamma)} - \sigma^{-\alpha}t^{\alpha + \gamma - 1}E_{\alpha,\alpha + \gamma}^{\alpha}\left(-\frac{t}{\sigma}\right)$（当 $\alpha = 0.5$ 时），呈现拉伸指数衰减特征。
• 分数阶 Cole-Davidson 模型表示为 $\varphi(t)_{CDF} = \frac{t^{\gamma-1}}{\Gamma(\gamma)} - \sigma^{-\beta}t^{\beta + \gamma - 1}E_{1,\beta + \gamma}^{\beta}\left(-\frac{t}{\sigma}\right)$，在 $\gamma = 1$ 时恢复经典 C-D 函数。
• 分数阶 Havriliak-Negami 模型表达为 $\varphi(t)_{HNF} = \frac{t^{\gamma-1}}{\Gamma(\gamma)} - \sigma^{-\alpha\beta}t^{\alpha\beta + \gamma - 1}E_{\alpha,\alpha\beta + \gamma}^{\beta}\left[-\left(\frac{t}{\sigma}\right)^\alpha\right]$，当 $\gamma = 1$ 时退化为经典 H-N 模型。
• 图形分析表明，分数阶松弛函数在 $t > 0$ 时保持完全单调性，且表现出依赖于 $\gamma$、$\alpha$ 和 $\beta$ 的非指数性、幂律类衰减特性。