D. del-Castillo-Negrete,B. A. Carreras,V. E. Lynch
Abstract:The use of reaction-diffusion models rests on the key assumption that the underlying diffusive process is Gaussian. However, a growing number of studies have pointed out the prevalence of anomalous diffusion, and there is a need to understand the dynamics of reactive systems in the presence of this type of non-Gaussian diffusion. Here we present a study of front dynamics in reaction-diffusion systems where anomalous diffusion is due to the presence of asymmetric Levy flights. Our approach consists of replacing the Laplacian diffusion operator by a fractional diffusion operator, whose fundamental solutions are Levy $\alpha$-stable distributions. Numerical simulation of the fractional Fisher-Kolmogorov equation, and analytical arguments show that anomalous diffusion leads to the exponential acceleration of fronts and a universal power law decay, $x^{-\alpha}$, of the tail, where $\alpha$, the index of the Levy distribution, is the order of the fractional derivative.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: the behavior of front dynamics in reaction - diffusion systems in the presence of non - Gaussian diffusion (especially anomalous diffusion caused by Levy flights). Specifically, the author focuses on how the front in the reaction - diffusion system propagates when the diffusion process no longer follows the traditional Gaussian assumption, especially how this asymmetric Levy flight affects the dynamic characteristics of right - moving and left - moving fronts.
### Background and Problem Description of the Paper
Traditional reaction - diffusion models rely on the assumption of Gaussian diffusion, that is, the random motion of particles can be described by the Laplacian operator \(\partial_x^2\). However, an increasing number of studies have shown that in some physical, chemical, and biological systems, the diffusion process may be anomalous, that is, the growth rate of the mean - square displacement \(\langle [x - \langle x \rangle]^2 \rangle \sim t^\gamma\) is different from that of Gaussian diffusion (\(\gamma = 1\)). In particular, when \(\gamma > 1\) it is called super - diffusion, and when \(\gamma < 1\) it is called sub - diffusion.
In this case, using the standard Gaussian diffusion model to describe the front dynamics of these systems may not be accurate enough. Therefore, understanding the front behavior of reaction - diffusion systems under non - Gaussian diffusion conditions has become an important research topic.
### Research Methods and Models
To study this problem, the author introduced the fractional - order diffusion operator \(aD_x^\alpha\) to replace the traditional Laplacian operator \(\partial_x^2\), where \(\alpha\) is the order of the fractional derivative and \(1 \leq \alpha < 2\). The solution of this fractional - order diffusion operator is the Levy - \(\alpha\) stable distribution, rather than the Gaussian distribution. In particular, the author considered asymmetric Levy flights, that is, the case where Levy flights are exhibited in the positive direction and Gaussian diffusion is exhibited in the negative direction.
The specific model equation is:
\[
\frac{\partial \varphi}{\partial t} = \chi aD_x^\alpha \varphi+F(\varphi)
\]
where \(F(\varphi)=\gamma \varphi(1 - \varphi)\) is the Fisher - Kolmogorov - type reaction term, which describes bistable reaction kinetics.
### Main Findings
Through numerical simulation and theoretical analysis, the author found the following:
1. **Right - moving Front**: The right - moving front exhibits exponential acceleration, and the tail shows algebraic decay \(\varphi \sim x^{-\alpha}\).
2. **Left - moving Front**: The left - moving front exhibits exponential decay and propagates at a constant speed.
3. **Universal Law**: For the right - moving front, its Lagrangian velocity \(v_L\) grows exponentially with time \(v_L \sim e^{\gamma t / \alpha}\), while for the left - moving front, its speed depends on the steepness of the initial conditions.
### Conclusions
This study reveals the complex behavior of front dynamics in reaction - diffusion systems under non - Gaussian diffusion conditions. These results not only help to understand the transport phenomena in plasma physics but may also be applied to processes involving long - range diffusion in biology and other fields.
### Application Prospects
The results of this study are of great significance for understanding the turbulence - shear flow interaction in plasma physics and the dynamics of the L - H transition front. In addition, the fractional - order Fisher - Kolmogorov equation has a certain similarity with the integro - differential equations used to describe long - range diffusion and spatial patterns in biology, and may be helpful in describing strongly non - local processes.