[论文解读] Fractional relaxation and wave equations for dielectrics characterized by the Havriliak-Negami response function

Renat T. Sibatov, D. V. Uchaikin|arXiv (Cornell University)|Aug 24, 2010

Fractional Differential Equations Solutions参考文献 29被引用 14

一句话总结

本文利用分数阶微积分推导出具有Havriliak-Negami响应的电介质的分数阶松弛方程，构建了蒙特卡洛算法以数值计算分数阶算子的作用，并通过实验充电-放电曲线验证了该模型。该方法能够捕捉记忆效应和非德拜松弛行为，数值结果与纸介质和电解质电容器的实验数据高度吻合。

ABSTRACT

A fractional relaxation equation in dielectrics with response function of the Havriliak-Negami type is derived. An explicit expression for the fractional operator in this equation is obtained and Monte Carlo algorithm for calculation of action of this operator is constructed. Relaxation functions calculated numerically according to this scheme coincide with analytical functions obtained earlier by other authors. The algorithm represents a numerical way of calculation in relaxation problems with arbitrary initial and boundary conditions. A fractional equation for electromagnetic waves in such dielectric media is obtained. Numerical results are in a good agreement with experimental data.

研究动机与目标

• 开发适用于表现出Havriliak-Negami响应特性的材料中电介质松弛的分数阶微分方程框架。
• 推导控制松弛动力学的分数阶算子的显式解析表达式。
• 构建蒙特卡洛算法，以在任意初始和边界条件下数值计算分数阶算子及其逆算子的作用。
• 在具有Havriliak-Negami响应的介质中建模电磁波传播，通过分数阶微积分引入记忆效应。
• 将理论模型与电介质电容器的实验充电-放电曲线进行验证。

提出的方法

• 推导出形式为 $[1 + ( au\textsf{D}_t)^eta]^eta \rho(t) = \rho(0)\tau^{-1}\rho(t)$ 的分数阶松弛方程，其中 $\textsf{D}_t^\beta$ 表示Caputo分数阶导数。
• 利用泛函分析和Bochner-Phillips关系，将分数阶算子表示为稳定概率分布的形式。
• 基于单边Lévy稳定随机变量的平均，构建蒙特卡洛算法以模拟分数阶算子的作用。
• 将相同的数值框架应用于推导具有Havriliak-Negami响应的电介质中电磁场的分数阶波动方程。
• 通过将模拟的松弛曲线与解析解及纸介质和电解质电容器的实验数据对比，验证数值解的正确性。
• 调整分数阶参数 $\beta$（设为1）和 $\theta$（充电时间）以匹配实验行为，观察到记忆再生效应。

实验结果

研究问题

• RQ1 如何推导出描述具有Havriliak-Negami响应函数的电介质松弛行为的分数阶微分方程？
• RQ2 此类电介质的松弛方程中，分数阶算子的显式形式是什么？
• RQ3 能否构建一种数值算法，以在任意初始和边界条件下计算分数阶算子及其逆算子的作用？
• RQ4 松弛行为如何依赖于充电历史预处理？这能否与记忆效应相联系？
• RQ5 分数阶波动方程的数值解与电介质松弛和波传播的实验数据在多大程度上一致？

主要发现

• 蒙特卡洛算法能够精确计算分数阶算子的作用，并复现了其他作者推导出的解析松弛函数。
• 当 $\beta = 1$ 和 $\beta \to 1$ 时，对纸介质电容器的松弛曲线数值模拟与实验数据高度一致，尤其在 $\theta = 300$ s 且 $\beta \to 1$ 时表现更佳。
• 当 $\theta = 3600$ s 时，模型成功捕捉到放电阶段从 $t^{-\beta}$ 到 $t^{-\beta-1}$ 的幂律衰减过渡，证实了理论预测。
• 该模型表现出一种‘记忆再生效应’，即在交叉时间之后，松弛行为依赖于充电历史，而这一现象在德拜模型中未被观察到。
• 所推导的分数阶波动方程与Jonscher普遍定律的早期结果一致，支持了电磁波传播中的记忆效应。
• 理论与实验放电曲线之间的一致性很强，仅在极长时间段出现微小偏差，表明分数阶模型具有高度鲁棒性。