[Paper Review] Propagation in Fisher-KPP type equations with fractional diffusion in periodic media

Xavier Cabré, Anne-Charline Coulon+1 more|arXiv (Cornell University)|Sep 21, 2012

Mathematical and Theoretical Epidemiology and Ecology ModelsMedicine2 references3 citations

TL;DR

This paper studies Fisher-KPP equations with fractional diffusion in periodic media, proving that level sets of solutions propagate exponentially in time with a speed determined by a periodic principal eigenvalue. Unlike classical results, the propagation speed is independent of spatial direction, a key distinction from the Freidlin-Gärtner formula for standard diffusion.

ABSTRACT

We are interested in the time asymptotic location of the level sets of solutions to Fisher-KPP reaction-diffusion equations with fractional diffusion in periodic media. We show that the speed of propagation is exponential in time, with a precise exponent depending on a periodic principal eigenvalue, and that it does not depend on the space direction. This is in contrast with the Freidlin-Gärtner formula for the standard Laplacian.

Motivation & Objective

• To analyze the asymptotic spreading behavior of solutions to Fisher-KPP equations with fractional Laplacian in periodic media.
• To determine how the propagation speed of level sets depends on the structure of the periodic medium and the fractional order α.
• To establish that the propagation speed is exponential in time and direction-independent, contrasting with the directional dependence in classical Laplacian-based models.
• To construct explicit sub- and supersolutions to rigorously bound the spreading front and derive the precise asymptotic speed.
• To extend known results from the homogeneous case (constant μ) to the heterogeneous, periodic setting with fractional diffusion.

Proposed method

• The authors use a sub- and supersolution method based on explicit ansatz functions involving the principal eigenfunction φ₁ of the operator (−Δ)^α − μ(x)I.
• They construct sub- and supersolutions of the form u̲(x,t) = aφ₁(x) / (|λ₁|⁻¹ + b̲(t)|x|^{d+2α}) and ū(x,t) = ↑φ₁(x) / (|λ₁|⁻¹ + b̄(t)|x|^{d+2α}), where b̲(t) and b̄(t) are time-dependent coefficients.
• The time evolution of the coefficients b̲(t) and b̄(t) is designed to match the exponential growth rate e^{|λ₁|t/(d+2α)} observed in the spreading front.
• The proof relies on verifying that these functions satisfy the PDE inequality (sub/super-solution property) using estimates on the fractional Laplacian and the principal eigenvalue λ₁.
• The maximum principle is applied to compare the constructed sub- and supersolutions with the actual solution u(x,t), yielding bounds on the level sets.
• The analysis is carried out for α < 1/2 in detail, with the general case deferred to a forthcoming paper [4].

Experimental results

Research questions

• RQ1 How does the presence of fractional diffusion (α ∈ (0,1)) affect the spreading speed of solutions in periodic media compared to classical diffusion?
• RQ2 Does the propagation speed in periodic media with fractional diffusion depend on the direction of propagation, as in the classical Freidlin-Gärtner formula?
• RQ3 Can explicit sub- and supersolutions be constructed to capture the exponential-in-time spreading behavior of level sets?
• RQ4 What role does the periodic principal eigenvalue λ₁ of (−Δ)^α − μ(x)I play in determining the asymptotic spreading rate?
• RQ5 How does the decay rate of the initial data influence the long-time behavior and spreading speed?

Key findings

• The level sets of solutions to the Fisher-KPP equation with fractional diffusion in periodic media propagate exponentially in time at a rate e^{|λ₁|t/(d+2α)}, where λ₁ < 0 is the principal eigenvalue.
• The spreading speed is independent of the spatial direction, a significant contrast to the Freidlin-Gärtner formula for standard Laplacian diffusion.
• For any λ ∈ (0, min μ), there exist constants c_λ > 0 and t_λ > 0 such that the set {x : u(x,t) = λ} is contained in the annular region {x : c_λ e^{|λ₁|t/(d+2α)} ≤ |x| ≤ c_λ⁻¹ e^{|λ₁|t/(d+2α)}} for t ≥ t_λ.
• The construction of sub- and supersolutions relies on a carefully chosen time-dependent coefficient b(t) that grows as e^{|λ₁|t/(d+2α)}, matching the spreading rate.
• The principal eigenvalue λ₁ determines the exponent of the exponential spreading, with |λ₁| scaling the growth rate.
• The method is robust and extends to general nonlinearities and more general diffusions, as noted in the concluding remarks for future work.