# American Institute of Mathematical Sciences

• Previous Article
Persistent two-dimensional strange attractors for a two-parameter family of Expanding Baker Maps
• DCDS-B Home
• This Issue
• Next Article
Convergence rate and stability of the split-step theta method for stochastic differential equations with piecewise continuous arguments
February  2019, 24(2): 671-694. doi: 10.3934/dcdsb.2018202

## Invasion fronts on graphs: The Fisher-KPP equation on homogeneous trees and Erdős-Réyni graphs

 1 Franklin W. Olin College of Engineering, Needham, MA 02492, USA 2 Department of Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA

Received  June 2017 Revised  January 2018 Published  June 2018

We study the dynamics of the Fisher-KPP equation on the infinite homogeneous tree and Erdős-Réyni random graphs. We assume initial data that is zero everywhere except at a single node. For the case of the homogeneous tree, the solution will either form a traveling front or converge pointwise to zero. This dichotomy is determined by the linear spreading speed and we compute critical values of the diffusion parameter for which the spreading speed is zero and maximal and prove that the system is linearly determined. We also study the growth of the total population in the network and identify the exponential growth rate as a function of the diffusion coefficient, α. Finally, we make predictions for the Fisher-KPP equation on Erdős-Rényi random graphs based upon the results on the homogeneous tree. When α is small we observe via numerical simulations that mean arrival times are linearly related to distance from the initial node and the speed of invasion is well approximated by the linear spreading speed on the tree. Furthermore, we observe that exponential growth rates of the total population on the random network can be bounded by growth rates on the homogeneous tree and provide an explanation for the sub-linear exponential growth rates that occur for small diffusion.

Citation: Aaron Hoffman, Matt Holzer. Invasion fronts on graphs: The Fisher-KPP equation on homogeneous trees and Erdős-Réyni graphs. Discrete & Continuous Dynamical Systems - B, 2019, 24 (2) : 671-694. doi: 10.3934/dcdsb.2018202
##### References:

show all references

##### References:
The linear spreading speed for (5), calculated numerically as a function of $\alpha$ for $k = 3$ (red), $k = 4$ (black) and $k = 5$ (blue). Note the critical values $\alpha_2(k)$ for which the spreading speed is zero and $\alpha_1(k)$ where the speed is maximal. Also note that as $\alpha\to 0$, these spreading speeds appear to approach a common curve
Critical rates of diffusion for period trees with period $m = 2$. On the left, we plot $\alpha_1$ as a function of $k_1$ with $k_2$ fixed to preserve the mean degree. On the right, we plot $\alpha_2$ as a function of $k_1$. Note that in both case the periodic heterogeneity increases the critical diffusion rates
Numerical simulations of (2) with $k = 3$ and for $\alpha = 0.2$ (left), $\alpha = 0.8$ (middle) and $\alpha = 2.2$ (right). The blue curves are $u_n(t)$ while the red curves depict the normalized population at each level, i.e. $w_n(t)/\max_n(w_n(t))$. Note that $0.2 < \alpha_1(3) < 0.8 < \alpha_2(3) < 2.2$. For $\alpha = 0.2$, we observe that the maximal population is concentrated at the front interface. For $\alpha = 0.8$, the maximal population is concentrated ahead of the front interface. Finally, for $\alpha = 2.2$ the local population at any fixed node converges to zero, but the total population grows and eventually is concentrated at the final level of the tree
On the left, we compare predictions for the exponential growth rate of the maximum of $w_n(t)$ as a function of $\alpha$ (blue line) against the exponential growth rates of $M(t)$ observed in direct numerical simulations (asterisks) for $k = 5$. On the right, we compare numerically observed spreading speeds for $w_n(t)$ (asterisks) versus linear spreading speeds determined numerically from the pinched double root criterion applied to $\tilde{d}_s(\gamma,\lambda)$ (blue line). Here we have taken $k = 5$
Arrival times for an Erdős-Réyni graph with $N = 60,000$ and expected degree $k_{ER} = 2$. Various values of $\alpha$ are considered. In green is the best fit linear approximation for the mean arrival times for nodes with distance between $3$ and $12$ from the initial location
Speed associated to the mean arrival times in numerical simulations on an Erdős-Réyni graph with $N = 60,000$ and expected degree $k_{ER} = 2$ are shown in asterisks. The blue curve is the spreading speed predicted by the analysis in Section 2 for the homogeneous tree with $k = 2.54$, found by numerically computing roots of (6). This value is chosen since it is one less than the mean degree of the network over those nodes with distance between $3$ and $12$ from the original location
Growth rate of the total population for Erdős-Rényi graph. On the left, $N = 500,000$ and $\alpha = 0.1,0.35,0.6,0.85$. Larger values of $\alpha$ correspond to faster growth rates. On the right is the case of $N = 60,000$ with the same values of $\alpha$
Numerically calculated exponential growth rate for the Erdős-Rényi graph. On the left, $N = 500,000$ and observed growth rates are plotted as circles. The asterisks are the corresponding growth rates in the homogeneous tree with depth $13$. The lower curve is degree $k = 3$ while the larger curve is degree $k = 4$. On the right are the same computations, but for the Erdős-Rényi graph with $N = 60,000$ and for homogeneous trees with $k = 2$ and $k = 3$
 [1] Gregoire Nadin. How does the spreading speed associated with the Fisher-KPP equation depend on random stationary diffusion and reaction terms?. Discrete & Continuous Dynamical Systems - B, 2015, 20 (6) : 1785-1803. doi: 10.3934/dcdsb.2015.20.1785 [2] Lina Wang, Xueli Bai, Yang Cao. Exponential stability of the traveling fronts for a viscous Fisher-KPP equation. Discrete & Continuous Dynamical Systems - B, 2014, 19 (3) : 801-815. doi: 10.3934/dcdsb.2014.19.801 [3] Matt Holzer. A proof of anomalous invasion speeds in a system of coupled Fisher-KPP equations. Discrete & Continuous Dynamical Systems - A, 2016, 36 (4) : 2069-2084. doi: 10.3934/dcds.2016.36.2069 [4] Matthieu Alfaro, Arnaud Ducrot. Sharp interface limit of the Fisher-KPP equation. Communications on Pure & Applied Analysis, 2012, 11 (1) : 1-18. doi: 10.3934/cpaa.2012.11.1 [5] Wenxian Shen, Zhongwei Shen. Transition fronts in nonlocal Fisher-KPP equations in time heterogeneous media. Communications on Pure & Applied Analysis, 2016, 15 (4) : 1193-1213. doi: 10.3934/cpaa.2016.15.1193 [6] Hiroshi Matsuzawa. A free boundary problem for the Fisher-KPP equation with a given moving boundary. Communications on Pure & Applied Analysis, 2018, 17 (5) : 1821-1852. doi: 10.3934/cpaa.2018087 [7] Christian Kuehn, Pasha Tkachov. Pattern formation in the doubly-nonlocal Fisher-KPP equation. Discrete & Continuous Dynamical Systems - A, 2019, 39 (4) : 2077-2100. doi: 10.3934/dcds.2019087 [8] Matthieu Alfaro, Arnaud Ducrot. Sharp interface limit of the Fisher-KPP equation when initial data have slow exponential decay. Discrete & Continuous Dynamical Systems - B, 2011, 16 (1) : 15-29. doi: 10.3934/dcdsb.2011.16.15 [9] Benjamin Contri. Fisher-KPP equations and applications to a model in medical sciences. Networks & Heterogeneous Media, 2018, 13 (1) : 119-153. doi: 10.3934/nhm.2018006 [10] François Hamel, James Nolen, Jean-Michel Roquejoffre, Lenya Ryzhik. A short proof of the logarithmic Bramson correction in Fisher-KPP equations. Networks & Heterogeneous Media, 2013, 8 (1) : 275-289. doi: 10.3934/nhm.2013.8.275 [11] Margarita Arias, Juan Campos, Cristina Marcelli. Fastness and continuous dependence in front propagation in Fisher-KPP equations. Discrete & Continuous Dynamical Systems - B, 2009, 11 (1) : 11-30. doi: 10.3934/dcdsb.2009.11.11 [12] James Nolen, Jack Xin. KPP fronts in a one-dimensional random drift. Discrete & Continuous Dynamical Systems - B, 2009, 11 (2) : 421-442. doi: 10.3934/dcdsb.2009.11.421 [13] Feng Cao, Wenxian Shen. Spreading speeds and transition fronts of lattice KPP equations in time heterogeneous media. Discrete & Continuous Dynamical Systems - A, 2017, 37 (9) : 4697-4727. doi: 10.3934/dcds.2017202 [14] Patrick Martinez, Jean-Michel Roquejoffre. The rate of attraction of super-critical waves in a Fisher-KPP type model with shear flow. Communications on Pure & Applied Analysis, 2012, 11 (6) : 2445-2472. doi: 10.3934/cpaa.2012.11.2445 [15] Aijun Zhang. Traveling wave solutions with mixed dispersal for spatially periodic Fisher-KPP equations. Conference Publications, 2013, 2013 (special) : 815-824. doi: 10.3934/proc.2013.2013.815 [16] Karel Hasik, Sergei Trofimchuk. Slowly oscillating wavefronts of the KPP-Fisher delayed equation. Discrete & Continuous Dynamical Systems - A, 2014, 34 (9) : 3511-3533. doi: 10.3934/dcds.2014.34.3511 [17] Tzong-Yow Lee and Fred Torcaso. Wave propagation in a lattice KPP equation in random media. Electronic Research Announcements, 1997, 3: 121-125. [18] Michiel Bertsch, Danielle Hilhorst, Hirofumi Izuhara, Masayasu Mimura, Tohru Wakasa. A nonlinear parabolic-hyperbolic system for contact inhibition and a degenerate parabolic fisher kpp equation. Discrete & Continuous Dynamical Systems - A, 2019, 0 (0) : 1-26. doi: 10.3934/dcds.2019226 [19] Mei Li, Zhigui Lin. The spreading fronts in a mutualistic model with advection. Discrete & Continuous Dynamical Systems - B, 2015, 20 (7) : 2089-2105. doi: 10.3934/dcdsb.2015.20.2089 [20] Yanni Zeng, Kun Zhao. On the logarithmic Keller-Segel-Fisher/KPP system. Discrete & Continuous Dynamical Systems - A, 2019, 39 (9) : 5365-5402. doi: 10.3934/dcds.2019220

2018 Impact Factor: 1.008