[论文解读] Nonlinear oscillations and stability domains in fractional reaction-diffusion systems

V. V. Gafiychuk, Bohdan Datsko+1 more|ArXiv.org|Feb 7, 2007

Nonlinear Dynamics and Pattern Formation参考文献 24被引用 13

一句话总结

本文研究具有三次非线性的分数阶反应-扩散系统，表明分数阶导数指数 α 与特征时间比 τ₁/τ₂ 共同决定了振荡不稳定性的发生。通过结合线性稳定性分析与数值模拟，本文表明当 α 超过临界阈值 α₀ 时，系统会涌现出均匀或非均匀振荡，从而破坏稳态耗散结构的稳定性。

ABSTRACT

We study a fractional reaction-diffusion system with two types of variables: activator and inhibitor. The interactions between components are modeled by cubical nonlinearity. Linearization of the system around the homogeneous state provides information about the stability of the solutions which is quite different from linear stability analysis of the regular system with integer derivatives. It is shown that by combining the fractional derivatives index with the ratio of characteristic times, it is possible to find the marginal value of the index where the oscillatory instability arises. The increase of the value of fractional derivative index leads to the time periodic solutions. The domains of existing of periodic solutions for different parameters of the problem are obtained. A computer simulation of the corresponding nonlinear fractional ordinary differential equations is presented. It is established by computer simulation that there exists a set of stable spatio-temporal tructures of the one-dimensional system under the Neumann and periodic boundary condition. The features of these solutions consist in the transformation of the steady state dissipative structures to homogeneous oscillations or space temporary structures at a certain value of fractional index and the ratio of characteristic times of system.

研究动机与目标

• 分析具有三次非线性的分数阶反应-扩散系统中的稳定性及非线性振荡的出现。
• 确定分数阶导数指数 α 与特征时间比 τ₁/τ₂ 如何影响从稳态耗散结构到振荡行为的转变。
• 通过线性和非线性分析，研究均匀振荡与时空振荡出现的条件。
• 在 Neumann 与周期性边界条件下，建立稳定时空结构的存在性及其性质。

提出的方法

• 采用阶数满足 0 < α < 2 的 Caputo 分数阶导数来描述反应-扩散系统的时间演化。
• 利用 Fréchet 导数与系统矩阵的特征值分解，对均匀平衡态进行线性稳定性分析。
• 通过相似变换将线性化系统转化为对角形式，以特征值 λ₁、λ₂ 分析稳定性。
• 通过改变 α 与 τ₁/τ₂ 的数值模拟对应的分数阶常微分方程（FODEs），研究非线性动力学行为。
• 在一维区域上施加 Neumann 与周期性边界条件，并以小的空间非均匀扰动作为初始条件。
• 追踪解的演化过程，以识别稳态结构、脉动非均匀振荡与均匀振荡之间的转变。

实验结果

研究问题

• RQ1 分数阶导数指数 α 如何影响具有三次非线性的反应-扩散系统中振荡不稳定的 onset？
• RQ2 激活剂与抑制剂特征时间之比 τ₁/τ₂ 在决定稳定性及所出现振荡类型方面起什么作用？
• RQ3 在何种条件下，稳态耗散结构会转变为均匀或非均匀振荡模式？
• RQ4 由 α、τ₁/τ₂ 与均匀解 ̄n₁ 定义的临界稳定性曲面如何确定稳定与不稳定区域之间的边界？

主要发现

• 存在一个临界值 α₀，当 α > α₀ 时，系统从稳态耗散结构转变为振荡动力学，振荡在临界稳定性曲面上出现。
• 当分数阶导数指数 α 与比值 τ₁/τ₂ 共同越过某一阈值时，振荡不稳定性产生，导致周期性时间解。
• 当 α < 1 时，亚扩散行为导致形成稳态耗散结构，随着 τ₁/τ₂ 减小，这些结构可演化为脉动的非均匀振荡。
• 当 α > 1 时，超扩散行为加速了振荡的出现，使稳态结构在较低 τ₁/τ₂ 值下失稳，形成均匀振荡。
• 均匀振荡的吸引域比耗散结构更广，其出现导致空间模式的破坏。
• 第二变量 n₂ 的振幅几乎保持不变，仅在其平衡值附近有微小偏离，表明其在动态转变中作用微弱。