[论文解读] Mathematical modeling of pattern formation in sub- and supperdiffusive reaction-diffusion systems

V. V. Gafiychuk, Bohdan Datsko+1 more|ArXiv.org|Nov 2, 2006

Nonlinear Dynamics and Pattern Formation参考文献 25被引用 7

一句话总结

本文利用分数阶时间导数（Caputo型，0 < α < 2）研究了亚扩散和超扩散反应-扩散系统中的模式形成。结果表明，分数阶指数 α 充当分岔参数：当 α < α₀ 时，系统表现出类似经典反应-扩散系统的稳定、稳态耗散结构；当 α > α₀ 时，出现振荡性或混沌的时空模式，且均匀振荡在较高 α 值下破坏了模式。临界指数 α₀ 标志着振荡不稳定的起始点，该结果通过解析方法确定，并经数值模拟验证。

ABSTRACT

This paper is concerned with analysis of coupled fractional reaction-diffusion equations. It provides analytical comparison for the fractional and regular reaction-diffusion systems. As an example, the reaction-diffusion model with cubic nonlinearity and Brusselator model are considered. The detailed linear stability analysis of the system with Cubic nonlinearity is provided. 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. Computer simulation and analytical methods are used to analyze possible solutions for a linearized system. A computer simulation of the corresponding nonlinear fractional ordinary differential equations is presented. It is shown that the increase of the fractional derivative index leads to the periodic solutions which become stochastic at the index approaching the value of 2. It is established by computer simulation that there exists a set of stable spatio-temporal structures of the one-dimensional system under the Neumann and periodic boundary condition. The characteristic 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.

研究动机与目标

• 分析分数阶时间导数对反应-扩散系统中模式形成的影响。
• 确定在亚扩散和超扩散系统中振荡不稳定性出现的临界分数阶指数 α₀。
• 比较分数阶反应-扩散系统与经典（整数阶）系统动力学的异同。
• 研究当 α 增大时，从稳态耗散结构向时空振荡的转变过程。
• 探讨在Neumann边界条件和周期性边界条件下，时空结构的鲁棒性与多样性。

提出的方法

• 采用 Caputo 时间导数（阶数 α ∈ (0,2)）构建耦合的分数阶反应-扩散方程。
• 通过在均匀平衡态附近线性化系统，并分析所得特征值问题，对分数阶系统进行线性稳定性分析。
• 通过解析方法推导出临界指数 α₀，即振荡不稳定性出现的阈值，其依赖于特征时间之比和分数阶指数。
• 利用数值模拟求解线性化系统，分析振荡解和混沌解的出现。
• 对非线性分数阶ODE和PDE（包括立方非线性项和Brusselator模型）进行计算机模拟，边界条件为Neumann和周期性边界条件。
• 采用谱方法分析空间模态，并追踪不同 α 值下扰动随时间的演化。

实验结果

研究问题

• RQ1 分数阶导数指数 α 如何影响反应-扩散系统中的稳定性与模式形成？
• RQ2 振荡不稳定性出现的临界值 α₀ 是什么？其如何由系统参数决定？
• RQ3 分数阶系统中的时空模式与经典（α=1）反应-扩散系统中的模式有何异同？
• RQ4 当 α 超过 α₀ 后，系统会涌现出何种动力学行为——振荡、混沌或均匀振荡？
• RQ5 为何在高 α 值下耗散结构会瓦解为均匀振荡？这种现象如何依赖于系统的非线性特性？

主要发现

• 临界分数阶指数 α₀ 充当分岔参数，将稳定、阻尼振荡（α < α₀）与不稳定、振荡或混沌动力学（α > α₀）分隔开。
• 当 α < α₀ 时，系统表现出与经典反应-扩散系统形式相同的稳态耗散结构，模式形态无显著偏差。
• 在 α = α₀ 时，系统发生过渡，进入持续振荡状态，且均匀振荡破坏了立方非线性模型中的空间模式。
• 当 α 超过 α₀ 后，系统从周期性振荡演化为准混沌的时空结构，Brusselator模型中尤为明显。
• 在Brusselator模型中，由于耗散结构具有较大振幅和非光滑轮廓，非均匀振荡在高 α 值下仍能持续存在，未坍缩为均匀振荡。
• 选择性分数阶模拟表明，第一个变量的分数阶导数对振荡的出现具有主导影响。