3 The free boundary problem
3.3 The criteria for spreading and vanishing
Here we first give the comparison principle. The proof is similar to the proof of [7, Lemma 3.5].
Lemma 3.9. Let h¯ ∈ C1([0,∞)), ¯u, ¯v ∈ C(D)×C2,1(D), with D := {(x,t) ∈ R2 : 0 < x <
h¯(t),t>0}. Assume that di >0,ri >0,ai >0,bi >0and(u, ¯¯ v, ¯h)satisfies
¯
ut−d1u¯xx+β1u¯x ≥u¯(r1−a1u¯), 0< x<h¯(t),
¯
vt−d2v¯xx+β2v¯x ≥v¯(r2−a2v¯), 0< x<h¯(t),
¯
u(0,t)≥0, ¯v(0,t)≥0, ¯u(h¯(t),t) =0, ¯v(h¯(t),t) =0, t>0, h¯0(t)≥ −µ[u¯x(h¯(t),t) +ρv¯x(h¯(t),t)], t>0.
(3.12)
Assume that di >0,ri >0,ai >0,b1> 0,b2 <0and(u, ¯¯ v, ¯h)satisfies
¯
ut−d1u¯xx+β1u¯x ≥u¯(r1−a1u¯), 0<x< h¯(t),
¯
vt−d2v¯xx+β2v¯x≥ v¯(r2−a2v¯−b2u¯), 0<x< h¯(t),
¯
u(0,t)≥0, ¯v(0,t)≥0, ¯u(h¯(t),t) =0, ¯v(h¯(t),t) =0, t>0, h¯0(t)≥ −µ[u¯x(h¯(t),t) +ρv¯x(h¯(t),t)], t>0.
(3.13)
Ifh¯(0) ≥ h0, ¯u(x, 0) ≥ u0(x)andv¯(x, 0) ≥ v0(x) on [0,h0], then the solution (u,v,h)of problem (1.5)satisfiesh¯(t)≥h(t)on[0,∞)andu¯ ≥u, ¯v≥v on[0,h(t)]×[0,∞).
Letλ(1i)(l)be the principle eigenvalue of the following problem fori=1, 2 (−diφxx+βiφx =λ(1i)(l)φ, 0< x<l,
φ(0) =φ(l) =0. (3.14)
It is well known thatλ1(i)(l) = 2dβi
i
2
+di πl2
is a strictly decreasing and continuous function inland
liml→0λ(1i)(l) =∞, lim
l→∞λ(1i)(l) = βi
2di 2
.
Theorem 3.10. Suppose that di >0,ri > 0,ai > 0,b1 > 0,b2 ∈ R, 0≤ βi <2√
diri for i=1, 2. If h∞ <∞, then h∞ ≤h∗ =min{Lri,i=1, 2}, where Lri satisfiesλ1(i)(Lri) =ri.
Proof. Due to Theorem3.6, limt→∞ku(·,t),v(·,t)kC1([0,h(t)]) = 0 ifh∞ < ∞. Assume h∞ > h∗ to get a contradiction. Ifh∞ > Lr1, then there exists ε >0 such thath∞ > Lr1−b1ε. For suchε, there existsT01 such thath(T0) =l> Lr1−b1ε andv(x,t)≤εfort ≥T0, 0≤ x≤h(t).
Letz= z(x,t)be the unique solution of
zt−d1zxx+β1zx =z(r1−a1z−b1ε), 0<x <l,t ≥T0, z(0,t) =z(l,t) =0, t ≥T0,
z(x,T0) =u(x,T0), 0≤x ≤l.
Applying the comparison principle,z(x,t)≤ u(x,t) fort ≥ T0, 0 ≤ x ≤ l. Since l > Lr1−b1ε, thenkz(·,t)−Z(·)kC([0,l])→0 ast →∞, whereZ(x)is the unique positive solution of
(−d1Zxx+β1Zx = Z(r1−a1Z−b1ε), 0< x<l, Z(0) =Z(l) =0.
lim inft→∞u(x,t)≥limt→∞z(x,t) =Z(x)>0 in(0,l). This is a contradiction. Ifh∞ >Lr2, we can get a contradiction by using the similar argument.
Lemma 3.11. Suppose that ri >0,ai > 0,bi >0, 0≤ βi <2√
diri,di =1for i=1, 2and h0 <h∗, then there existsµ>0depending on u0and v0such that h∞ < ∞ifµ≤µ.
Proof. Sinceh0 <h∗, then λ(1i)(h0) = 2dβi
i
2
+di 2πh
0
2
> ri (i =1, 2). We can chooseδ,γsmall such that
1 h0(1+δ)
diπ2
h0(1+δ)−δγh0
+ β
2i
4 −γ−ri >0.
Define where M is a positive constant to be determined.
Direct computations yield
Lemma 3.12. Suppose that di >0,ri >0,ai > 0,bi ∈R, 0≤ βi <2√
diri for i= 1, 2and h0 <h∗, there existsµ¯ >0such that h∞ =∞ifµ≥µ.¯
Proof. Due to the boundedness ofuandv, there existsδ∗ such that
u(r1−a1u−b1v)≥ −δ∗u, v(r2−a2v−b2u)≥ −δ∗v.
Consider the following problem
wt−d1wxx+β1wx =−δ∗w, 0< x<r(t),t>0, zt−d2zxx+β2zx =−δ∗z, 0< x<r(t),t>0, w(0,t) =z(0,t) =0, t>0,
w(r(t),t) =z(r(t),t) =0, t>0, r0(t) =−µ[wx(r(t),t) +ρzx(r(t),t)], t>0, w(x, 0) =u0(x),z(x, 0) =v0(x),r(0) =h0, 0< x<h0.
(3.15)
Similar to Lemma 3.1, such problem admits a unique global solution (w,z,r). Applying the comparison principle, it follows that
u(x,t)≥w(x,t), v(x,t)≥ z(x,t), h(t)≥r(t), forx ∈[0,r(t)], t> 0. (3.16) Next, we prove that for all largeµ,r(1)≥h∗. Choose a smooth functionr(t)such that
r(0) =h0/2, r(1) =h∗, r0(t)>0, fort>0.
Consider the following initial-boundary value problem
wt−d1wxx+β1wx =−δ∗w, 0< x<r(t),t>0, zt−d2zxx+β2zx =−δ∗z, 0< x<r(t),t>0, w(0,t) =z(0,t) =0, t>0,
w(r(t),t) =z(r(t),t) =0, t>0,
w(x, 0) =w0(x),z(x, 0) =z0(x), 0< x<h0/2,
(3.17)
here(w0(x),z0(x))satisfies
0<w0(x)≤u0(x) on [0,h0/2], w0(0) =w0(h0/2) =0, w00(h0/2)<0, 0<z0(x)≤v0(x) on[0,h0/2], z0(0) =z0(h0/2) =0, z00(h0/2)<0.
The standard theory for parabolic equations ensures that (3.17) has an unique positive solution (w,z) and wx(r(t),t) < 0,zx(r(t),t) < 0 for all t ∈ [0, 1] due to the Hopf Lemma.
Then there exists a constant ¯µ>0 such that for allµ≥µ,¯
r0(t)≤ −µ[wx(r(t),t) +ρzx(r(t),t)] fort∈[0, 1]. (3.18) Since the choice of initial values and (3.15)–(3.18), we have
r(t)≥r(t), w(x,t)≥w(x,t),z(x,t)≥z(x,t), for x∈[0,r(t)], t ∈[0, 1],
which implies r(1) ≥ r(1) = h∗. In view of (3.16), h∞ > h(1) ≥ h∗. Together with Theo-rem3.10, derives the desired result.
Theorem 3.13. Assume that ri >0,ai >0,bi >0, 0≤ βi <2√
diri,di =1for i=1, 2and h0 <h∗, there exists µ∗ ≥ µ∗ > 0 such that vanishing happens (h∞ < ∞) if 0 < µ ≤ µ∗ or µ = µ∗, and spreading happens (h∞= ∞) ifµ>µ∗.
Proof. DefineΓ := {µ > 0 : h∞ ≤ h∗}. Due to Lemma3.11, Γ 6= ∅. In view of Lemma3.12, µ∗ :=supΓ∈ [µ, ¯µ]. By the definition ofµ∗ and Theorem 3.10, we get thath∞= ∞ifµ>µ∗.
Next, we prove that h∞ < ∞ if µ = µ∗. If not, h∞ = ∞ for µ = µ∗. So there exists T such thath(T)> h∗. Since the solution(u,v,h)depends onµ, we write(uµ,vµ,hµ)instead of (u,v,h). By the continuous dependence of(uµ,vµ,hµ)on µ, for smallε >0,hµ(T)>h∗for all [µ∗−ε,µ∗+ε]. Then supΓ≤ µ∗−ε, which contradicts to the definition ofµ∗. Henceµ∗ ∈ Γ.
DenoteΛ := {ν > 0 : ν ≥ µsuch thath∞ ≤ h∗ for all 0 < µ≤ ν}andµ∗ := supΛ≤ µ∗. Using the similar way to the above, we obtain thatµ∗ ∈Λ. The proof is completed.
4 Acknowledgments
This work is supported by the National Natural Science Foundation of China (11771407), the Innovative Research Team of Science and Technology in Henan Province (17IRTSTHN007), and the National Key R&D Program of China (2017YFB0702504).
References
[1] A. A. Berryman, The orgins and evolution of predator–prey theory, Ecology 73(1992), No. 5, 1530–1535.https://doi.org/10.2307/1940005
[2] C. Bianca, Mathematical modeling for keloid formation triggered by virus: malignant effects and immune system competition,Math. Models Methods. Appl. Sci.21(2011), No. 2, 389–419.https://doi.org/10.1142/S021820251100509X;MR2776673
[3] C. Bianca, M. Pennisi, S. Motta, M.A. Ragusa, Immune system network and cancer vaccine, AIP Conference Proceedings 1389(2011), No. 1, 945–948. https://doi.org/10.
1063/1.3637764
[4] C. Bianca, F. Pappalardo, S. Motta, M. A. Ragusa, Persistence analysis in a Kolmogorov-type model for cancer-immune system competition,AIP Conference Proceed-ings 1558(2013), No. 1, 1797–1800.https://doi.org/10.1063/1.4825874
[5] R. S. Cantrell, C. Cosner, Spatial ecology via reaction–diffusion equations, Wiley Series in Mathematical and Computational Biology, John Wiley & Sons Ltd., Chichester, 2003.
https://doi.org/10.1002/0470871296;MR2191264
[6] J. Dockery, V. Hutson, K. Mischaikow, M. Pernarowski, The evolution of slow dis-persal rates: a reaction diffusion model, J. Math. Biol. 37(1998), No. 1, 61–83. https:
//doi.org/10.1007/s002850050120;MR1636644
[7] Y. Du, Z. Lin, Spreading–vanishing dichotomy in the diffusive logistic model with a free boundary,SIAM J. Math. Anal.42(2010), No. 1, 377–405.https://doi.org/10.1137/
090771089;MR2607347
[8] J. Guo, C. Wu, On a free boundary for a two-species weak competition system, J. Dyn. Differential Equations 24(2012), No. 4, 873–895. https://doi.org/10.1007/
s10884-012-9267-0;MR3000608
[9] X. He, W. Ni, Global dynamics of the Lotka–Volterra competition–diffusion system: dif-fusion and spatial heterogeneity I, Comm. Pure. Appl. Math. 69(2016), No. 5, 981–1014.
https://doi.org/10.1002/cpa.21596;MR3481286
[10] C. Lei, Z. Lin, H. Wang, The free bondary problem describing information diffusion in online social networks, J. Differential Equations 254(2013), No. 3, 1326–1341. https:
//doi.org/10.1016/j.jde.2012.10.021;MR2997373
[11] Y. Kaneko, Y. Yamada, A free boundary problem for a reaction diffusion equation ap-pearing in ecology,Adv. Math. Sci. Appl.21(2011), No. 2, 467–492.MR2953128
[12] Y. Lou, F. Lutscher, Evolution of dispersal in open advective environments, J. Math.
Biol. 69 (2014), No. 6–7, 1319–1342. https://doi.org/10.1007/s00285-013-0730-2;
MR3275198
[13] F. Lutscher, E. McCauley, M. A. Lewis, Spatial patterns and coexistence mechanisms in systems with unidirectional flow, Theor. Popul. Biol. 71(2007), No. 3, 267–277. https:
//doi.org/10.1016/j.tpb.2006.11.006
[14] Y. Kuang, E. Beretta, Global qualitative analysis of a ratio-dependent predator–prey sys-tem, J. Math. Biol.36(1998), No. 4, 389–406. https://doi.org/10.1007/s002850050105;
MR1624192
[15] K. Y. Lam, W. Ni, Uniqueness and complete dynamics in heterogeneous competition diffusion systems,SIAM J. Appl. Math.72(2012), No. 6, 1695–1712.https://doi.org/10.
1137/120869481;MR3022283
[16] Y. Lou, D. Xiao, P. Zhou, Qualitative analasis for a Lotka–Volterra competition system in advective homogeneous environment, Discrete Contin. Dyn. Syst. A. 36(2016), No. 2, 953–969.https://doi.org/10.3934/dcds.2016.36.953;MR3392913
[17] Y. Lou, P. Zhou, Evolution of dispersal in advective homogeneous environment: the effect of boundary conditions, J. Differential Equations 259(2015), No. 1, 141–171. https:
//doi.org/10.1016/j.jde.2015.02.004;MR3335923
[18] M. H. Protter, H. F. Weinberger, Maximum principles in differential equations, Springer-Verlag, New York, 1984.https://doi.org/10.1007/978-1-4612-5282-5;MR0762825 [19] C. V. Pao, Nonlinear parabolic and elliptic equations, Springer Science & Business Media,
New York, 1992.MR1212084
[20] J. Ren, X. Li, Bifurcations in a seasonally forced predator–prey model with generalized Holling type IV functional response, Int. J. Bifurcation Chaos 26(2016), No. 12, 1650203.
https://doi.org/10.1142/S0218127416502035;MR3574812
[21] J. Ren, X. Li, How seasonal forcing influences the complexity of a predator–prey system, Discrete Contin. Dyn. Syst. B.23(2018), No. 2, 785–807.https://doi.org/10.3934/dcdsb.
2018043
[22] J. Ren, L. Yu, Codimension-two bifurcation, chaos and control in a discrete-time informa-tion diffusion model,J. Nonlinear Sci. 26(2016), No. 6, 1895–1931. https://doi.org/10.
1007/s00332-016-9323-8;MR3562398
[23] H. L. Smith, Monotone dynamical systems. An introduction to the theory of competitive and cooperative systems,Mathematical Surveys and Monographs, Vol. 41 American Mathemat-ical Society, Providence, 1995.MR1319817
[24] M. Wang, On some free boundary problems of the prey-predator model, J. Differential Equations256(2014), No. 10, 3365–3394. https://doi.org/10.1016/j.jde.2014.02.013;
MR3177899
[25] F. Wang, H. Wang, K. Xu, Diffusive logistic model towards predicting information diffu-sion in online social networks, in: 32nd International Conference on Distributed Computing Systems Workshops (ICDCS Workshops)2012, 133–139.https://doi.org/10.1109/ICDCSW.
2012.16
[26] M. Wang, J. Zhao, Free boundary problems for a Lotka–Volterra competition sys-tem, J. Dyn. Differential Equations 26(2014), No. 3, 655–672. https://doi.org/10.1007/
s10884-014-9363-4;MR3274436
[27] M. Wang, Y. Zhang, Two kinds of free boundary problems for the diffusive prey–
predator model,Nonlinear Anal. Real World Appl. 24(2015), 73–82.https://doi.org/10.
1016/j.nonrwa.2015.01.004;MR3332883
[28] X. Zhao,Dynamical systems in population biology, Springer, New York, 2003.MR1980821 [29] P. Zhou, On a Lotka–Volterra competition system: diffusion vs advection, Calc. Var.
55(2016), No. 6.https://doi.org/10.1007/s00526-016-1082-8;MR3566937
[30] L. Zhou, S. Zhang, Z. Liu, A free boundary problem of a predator–prey model with advection in heterogeneous environment, Appl. Math. Comput. 289(2016), 22–36. https:
//doi.org/10.1016/j.amc.2016.05.008;MR3515835