\qquad 随着现代科学技术的发展,医疗水平的提高以及人类文明的不断进步, 人类已经能够有效地预防和控制诸如天花、麻风、霍乱等曾经肆虐全球的传染性疾病, 但是一些新的、不断变异的传染病病毒却悄悄向人类袭来. 这些疾病的蔓延严重影响了人类的生活和工作. 20世纪80年代一种更为险恶的艾滋病毒正跨越国界在全球蔓延; 特别是2003年春SARS[1]病毒突袭人类,给人们的生命财产带来极大的危害. 因此,人们越来越重视对传染病的传播机理,预报和控制的研究.
为了研究如SARS 这种具有潜伏期的传染病问题,人们常用SEI模型[2, 3, 4, 5]来描述, 即根据疾病传播过程不同的状态,把所要研究的群体分成三种类别: S类,表示未被感染但有可能被感染的一类人,即易感者(Suscepible); E类,表示已被感染且具有传染力但未出现病症的一类人,即病毒携带者(Exposed)或带菌者; I 类,表示已被感染成病人具有传染力的一类人,即染病者(Infectious). 这样总人口数N=S+E+I. 如果潜伏者(E)和染病者(I)一样具有传染性, 正如2003年出现的SARS病毒,则相应的微分系统是 {dSdt=A−λ1SNE−λ2SNI−dS,[3mm]dEdt=λ1SNE+λ2SNI−γE−(d+α1)E,[3mm]dIdt=γE−kI−(d+α2)I, 其中正常数d表示自然死亡率; 非负常数α1,α2分别表示潜伏者和染病者的因病死亡率; A是易感染者的常数输出率; λ1,λ2 分别是潜伏者和染病者的有效接触率; 参数γ简称为转移率, 1/γ为平均潜伏期; 常数k表示隔离率.
对于模型(1.1),文献[4, 5]的作者证明了当基本再生数小于等于1时, 问题的无病平衡点是全局渐近稳定的,当基本再生数大于1时, 染病平衡点则是全局渐近稳定的. 在模型(1.1)的基础上,结合传染病动力学一些特征, 如种群在其所处环境中分布的不均匀性,种群内部流动的不均性等,该模型可扩展为 {u1t−d1Δu1=A−u1(λ1u2u1+u2+u3+λ2u3u1+u2+u3+d),x∈Ω, t>0,[3mm]u2t−d2Δu2=λ1u1u2u1+u2+u3+λ2u1u3u1+u2+u3−(γ+d+α1)u2,x∈Ω, t>0,[3mm]u3t−d3Δu3=γu2−(k+d+α2)u3,x∈Ω, t>0,[2mm]∂u1∂η=∂u2∂η=∂u3∂η=0,x∈∂Ω, t>0,[2mm]ui(x,0)=ui,0(i=1,2,3),x∈¯Ω, 其中Ω是Rn中的有界固定区域,u1(x,t)表示易感者在t时刻的空间分布密度, u2(x,t)表示病毒潜伏者在t时刻的空间分布密度, u3(x,t)表示染病者的分布密度; 扩散的引入允许了所有人口的迁移, 正常数d1,d2,d3分别表示易感者,潜伏者,染病者的空间扩散率. 边界∂Ω光滑,η是边界上的单位外法向量; 这里引入的齐次Neumann边界条件说明上述系统是封闭的,即在边界上没有人口移动. 初始函数ui,0在¯Ω 上是非负有界的连续函数, 根据生物意义,我们假定在¯Ω上ui,0(x)>0.
对固定区域Ω构建的模型(1.2),我们知道不管初始时刻传染源被限制在Ω 内多小的区域内,问题的解总是正的, 这表明疾病都会立即传播到整个固定的区域. 显然这样的过程不能用来描述传染病传播的实际情况, 而自由边界却能很好地解决这个问题, 它能清晰地描述所研究问题的边界随着时间的推移而变化, 而不是从疾病开始爆发就立即能传播到整个区域. 最初的自由边界问题是Josef Stefan在研究极地海洋中冰的溶解问题时所建立的数学模型[6]. 这类问题在许多领域的研究中都会碰到,如图像处理[7], 化学中热裂解碳的蒸汽渗透[8],医学中的肿瘤治愈[9]等. 由于研究问题的边界是随时间而不断变化的,因此数学上求解自由边界问题, 移动边界的刻画就显得尤为重要, 它将作为解的一部分随解一起给出. 近几年自由边界开始被用来描述动物种群的扩张[10]. 动物种群的扩张现象是Skellam[11]在1953年发现的, 他首先结合野外获取的数据,计算出麝鼠活动区域的面积,进而求平方根, 即给出扩张半径,将这些数据按照年份划分, 结果发现采集的数据点落在同一条直线上,故他认为麝鼠的扩张符合线性增长. 由于扩张的存在, 我们在其它动物种群也会发现类似的现象. 由此可见这种扩张是普遍存在的, 基于这种有趣的现象, 为了能进一步解释这种扩张,许多科学工作者致力于模拟种群扩张, 尽管这种数学模拟是比较困难的,但还是取得了很多成果(参见文献[12, 13, 14, 15, 16]).
本文尝试着引入自由边界来描述传染病的蔓延. 为了方便, 只考虑λ2=0的情形,即忽略与染病者的接触,并且假设区域是球对称的. 令 r=|x|, x∈Rn, 考虑下面问题的正解(u1(r,t),u2(r,t),u3(r,t);h(t))的性质,记 N(r,t)=u1+u2+u3, 建立如下自由边界问题 {u1t−d1Δu1=A−u1(r,t)(λ1u2(r,t)N(r,t)+d),0<r, t>0,[3mm]u2t−d2Δu2=λ1u1(r,t)u2(r,t)N(r,t)−δu2(r,t),0<r<h(t), t>0,[2mm]u3t−d3Δu3=γu2(r,t)−ωu3(r,t),0<r<h(t), t>0,u1r(0,t)=u2r(0,t)=u3r(0,t)=0,t>0,u2(r,t)=u3(r,t)=0,h(t)≤r,t≥0,h′(t)=−μu2r(h(t),t), h(0)=h0>0,t>0,ui(r,0)=ui,0(i=1,2,3),r≥0, 其中uit=∂ui∂t, Δui=(ui)rr+n−1r(ui)r (i=1,2,3), δ=γ+d+α1,ω=k+d+α2. r=h(t)是需要确定的移动边界,h0,μ均为正常数, 初始函数ui,0≥0 (i=1,2,3) 且满足 {u1,0∈C2([0,+∞))⋂L∞([0,+∞)),u2,0,u3,0∈C2[0,h0],u2,0(r)=u3,0(r)=0,r∈[h0,+∞),u2,0(r)>0,r∈[0,h0). 生态学上,该模型表示初期染病者和带菌者只在[0,h0)上, 这些人口仅能从初始区域的边界r=h(t)向外迁移, 而在自由边界r=h(t)之外只存在易感者,没有带菌者和染病者. h′(t)=−μu2r(h(t),t)是为研究种群扩散而构造的Stefan条件[17, 18]. 这儿该条件表明穿过自由边界的感染者人口数量与相应的边界移动的距离成正比, 其中正常数μ 表示病毒在新区域的传播能力.
本文我们将分别考虑在固定区域Ω和移动区域上引入扩散项的SEI模型 并研究相应问题解的渐近性态. 第2节将给出系统(1.2) 解的基本性质,讨论偏微分方程组解的局部稳定性和全局稳定性; 第3节考虑自由边界问题(1.3)的适定性,最后第4节着重研究自由边界的性质, 给出疾病蔓延或消退的条件.
由于系统(1.2)初始条件ui,0>0且其三个方程右端的反应函数在R3+中充分光滑, 根据标准的偏微分方程理论, 问题(1.2)存在唯一全局解(u1(x,t),u2(x,t),u3(x,t)),(x,t)∈¯Ω×(0,∞). 下面首先给出问题解的正性和一致有界性.
定理2.1 设(u1(x,t),u2(x,t),u3(x,t))∈[C(¯Ω×[0,∞))⋂C2,1(Ω×(0,∞))]3是问题(1.2)的解, 则 0<ui(x,t)≤Mi (i=1,2,3), 其中 M1=max{‖u1,0(x)‖L∞(Ω),Ad}, M2=max{‖u2,0(x)‖L∞(Ω),λ1+λ2δM1}, M3=max{‖u3,0(x)‖L∞(Ω),γωM2}.
证 由比较原理和强极值原理得解的正性. 下面只验证上界. 注意到u1满足 {u1t−d1Δu1≤A−du1,x∈Ω, t>0,∂u1∂η=0,x∈∂Ω, t>0,[2mm]u1(x,0)=u1,0>0,x∈Ω. 令V1(t)是问题 V′(t)=A−dV1, t>0,V1(0)=‖u1,0(x)‖L∞(Ω) 的解,显然有 V1(t)≤max{‖u1,0(x)‖L∞(Ω),Ad}:=M1, 又由标准的抛物比较原理得 u1(x,t)≤V1(t)≤M1.
同理,既然 u2t−d2Δu2≤λ1u1+λ2u1−δu2≤(λ1+λ2)M1−δu2, 由比较原理知u2(x,t)≤max{‖u2,0(x)‖L∞(Ω),λ1+λ2δM1}:=M2. 注意到u3t−d3Δu3≤γM2−ωu3,故u3(x,t)≤max{‖u3,0(x)‖L∞(Ω),γωM2}:=M3.
现在讨论系统(1.2)平衡解的局部稳定性和全局稳定性. 显然方程组(1.2)总有一个无病平衡点P0(Ad,0,0),此时系统没有带菌者和染病者. 进一步地, 若 R0:=λ1ω+λ2γωδ>1, 系统(1.2)还有一个染病平衡点P∗=(u∗1,u∗2,u∗3),其中 u∗1=N∗A[(δ−d)ω−γd](λ1ω+λ2γ)(A−dN∗)+dN∗[(δ−d)Ω−γd],u∗2=ω(A−dN∗)(δ−d)Ω−γd, u∗3=γ(A−dN∗)(δ−d)Ω−γd,N∗=u∗1+u∗2+u∗3. 传染病学上,R0被称为基本再生数,通常指一个病人在平均患病期内所传染的人数. 如果把R0拆成两项R0=λ1δ+λ2γωδ, 只考虑对染病者进行隔离,即k充分大,从而使ω充分大,所以λ2γωδ 就充分小了,因此可以忽略第二项, 这样R0取决于第一项λ1δ. 这表示当有效接触率λ1大或平均潜伏期1γ较长时(即γ很小, 从而δ小),系统(1.2)存在一个染病平衡点P∗=(u∗1,u∗2,u∗3).
考虑染病平衡点P∗的性质,利用文献[20]中有关空间分解的方法,令0=μ1<μ2<μ3<⋯是具有齐次Neumann边界条件的算子−Δ在Ω上的特征值, 且E(μi)是μi在C1(Ω)上的特征空间. 令 X={U=(u1,u2,u3)∈[C1(¯Ω)]3|∂U∂η=0,x∈∂Ω}, {ϕij,j=1,1,⋯,dimE(μi)}是E(μi)的基本正交基, Xij={c⋅ϕij|c∈R3}. Xi=dimE(μi)⨁j=1Xij, X=∞⨁i=1Xi.
下面将方程组进行线性化: 先取U∗=(u∗1,u∗2,u∗3),令U=U∗+V, 代入系统(1.2)可得Vt=LV:=DΔV+FU(u∗)U, 其中D=diag(d1,d2,d3),FU(U∗)是矩阵{aij}: a11=−N∗−u∗1N∗2(λ1u∗2+λ2u∗3)−d,a12=−u∗1[λ1N∗−(λ1u∗2+λ∗3)]N∗2, a13=−u∗1[λ2N∗−(λ1u∗2+λ2u∗3)]N∗2,a21=N∗−u∗1N∗2(λ1u∗2+λ∗2u∗3), a22=u∗1[λ1N∗−(λ1u∗2+λ2u∗3)]N∗2−δ,a23=u∗1[λ2N∗−(λ1u∗2+λ2u∗3)]N∗2, a31=0,a32=γ,a33=−ω. 对每个i≥1,Xi是算子L的不变子空间,则算子作用在Xi 上为 Vt=LV=−μiDV+FU(U∗)V. 这样φi(λ):=|λI+μiD−FU(U∗)|=λ3+Aiλ2+Biλ+Ci=0为其对应的特征方程,其中 Ai=(d1+d2+d3)μi−(a11+a22+a33), Bi=(d1d2+d2d3+d3d1)μ2i−[(a22+a33)d1+(a11+a33)d2+(a11+a22)d3]μi+(a11a22−a12a21)+(a22a33−a23a32)+a11a33, Ci=d1d2d3μ3i−(a11d2d3+a22d1d3+a33d1d2)μ2i+[(a11a22−a21a12)d3+(a22a33−a23a32)d1+a11a33d2]μi−a33(a11a22−a12a21)+a32(a23a11−a21a13). 首先,注意到a11<0,a21>0,且a33<0. 进一步地,由于 u∗1=δu∗2N∗λ1u∗2+λ2u∗3,u∗2=ωγu∗3, 由此我们可以得到 a22=u∗1[λ1N∗−(λ1u∗2+λ2u∗3)]N∗2−δ≤−λ2δγλ1ω+λ2γ<0, a11a22−a21a12=(−a21−d)(−a12−δ)−a21a12=−da22+δa21>0, a22a33−a32a23=(ω+γ)δu∗2N∗>0, (−a33)(a11a22−a21a12)+a32(a23a11−a21a13)=ωδa21+d(a22a33−a32a23)>0, 这样就使得对任意的i≥1,有Ai>0,Bi>0和Ci>0. 另一方面, Hi:=AiBi−Ci=k1μ3i+k2μ2i+k3μi+k4,其中k1>0,k2>0且k3>0, 这儿我们省略了详细的计算. 进一步,我们还有 k4=−(a11+a22+a33)(a11a22−a21a12+a22a33−a23a32+a33a11)−[(−a33)(a11a22−a21a12)+a32(a23a11−a21a13)]=(ω+d+λ2δγλ1ω+λ2γ+λ1u∗2+λ2u∗3N∗)[−da22+δa21+δ(ω+γ)u∗2N∗+ωd+ωN∗−u∗1N∗2(λ1u∗2+λ2u∗3)]−ωδa21−dδ(ω+γ)u∗2N∗>0. 因此对i≥1,Hi>0. 根据Routh-Hurwitz准则,特征方程φi(λ)=0的三个根λi,1,λi,2,λi,3都具有负实部. 利用方程根的连续性不难证明, 所有特征根的实部都小于一个负数,即存在ε>0使得
下面考虑无病平衡点P0=(A/d,0,0)的性质,我们要证明当R0>1时,P0是不稳定的. 事实上,这里 L=D△+FU(P0)=D△+{aij},其中 a11=−d,a12=−λ1,a13=−λ2,a21=0,a22=λ1−δ,a23=λ2,a31=0,a32=γ,a33=−ω. 线性化模型(1.2)得Vt=LV. 取i=1,这样μi=0,对应的特征方程为 φ1(λ):=|λI−FU(P0)|=λ3+A1λ2+B1λ+C1=0, 其中 C1=d[(δ−λ1)ω−λ2γ]=dωδ(1−R0)<0. 当λ→+∞时, φ1(λ)→+∞,而φ1(0)=C1<0,故由介值定理特征方程 φ1(λ)=0至少存在一个正根,于是算子L的谱至少包含一个实部大于0的特征值, 再依据文献[20]的定理5.1.3得,当R0>1时,P0是不稳定的.
类似地,我们可以证明当R0<1时,P0是局部渐近稳定的.
定理2.2 当R0>1时,系统(1.2)的染病平衡点P∗是局部渐近稳定的,而无病平衡点P0是不稳定的. 相反,当R0<1时,无病平衡点P0是局部渐近稳定的,此时不存在染病平衡点.
下面我们用构造Lyapunov函数的方法,证明无病平衡点P0是全局渐近稳定的.
定理2.3 当R0<1时,系统(1.2)的无病平衡点P0是全局渐近稳定的.
证 取Lyapunov函数为 V(t)=∫Ω[12k1(u1−u∗1)2+u2+k2u3]dx, 其中k1,k2待定,注意到 V(t)≥min{12k1,1,k2}∫Ω[(u1−u∗1)2+u2+u3]dx≥0, 当且仅当(u1,u2,u3)取(Ad,0,0)时,V(t)取得最小值0. 沿着系统的轨线求导得 dVdt=∫Ωk1(u1−u∗1)u1t+u2t+k2u3t]dx=∫Ωk1(u1−u∗1)d1Δu1+d2Δu2+k2d3Δu3]dx+∫Ω{k1(u1−u∗1)(A−du1)+(k2γ−δ)u2−k2ωu3}dx+∫Ωλ1u1u2+λ2u1u3N[1−k1(u1−u∗1)]dx:=E+F, 其中 E=∫Ωk1(u1−u∗1)d1Δu1+d2Δu2+k2d3Δu3]dx=−k1d1∫Ω|∇u1|2dx, F=∫Ω{k1(u1−u∗1)(A−du1)+(k2γ−δ)u2−k2ωu3}dx+∫Ωλ1u1u2+λ2u1u3N[1−k1(u1−u∗1)]dx=∫Ω{−dk1(u1−u∗1)2+(k2γ−δ)u2−k2ωu3+λ1u1u2+λ2u1u3N[1−k1(u1−u∗1)]}dx. 首先E≤0,要使dVdt≤0,只要F≤0, 即选取适当的k1和k2使得 k2γ−δ+(1+k1Ad)λ1<0,−k2ω+(1+k1Ad)λ2<0, 于是取k1和k2使得 0<k1<d(1−R0)AR0,λ2ω<k2<δ−λ1γ, 由此得到dVdt=E+F≤0,由文献[21]中的引理2.5.3得到 limt→∞∫Ω(u1−u∗1)2dx=limt→∞∫Ωu22dx=limt→∞∫Ωu23dx=0, 再利用一致估计及紧性可得无病平衡点P0是全局渐近稳定的,详见文献[20].
下面我们应用压缩映像原理证明问题(1.3)解的局部存在性及唯一性.
定理3.1 对于满足(1.4)式的初值 (u1,0,u2,0,u3,0)及任意β∈(0,1),存在T>0使得问题(1.3)存在唯一解 (u1,u2,u3;h)∈C1+β,(1+β)/2(D∞T)×[C1+β,(1+β)/2(DT)]2×C1+β/2([0,T]), 此外, ‖u1‖C1+β,(1+β)/2(D∞T)+‖u2‖C1+β,(1+β)/2(DT)+‖u3‖C1+β,(1+β)/2(DT)+‖h‖C(1+β)/2([0,T]≤C, 这里D∞T={(r,t)∈R2: r∈[0,+∞),t∈[0,T])}, DT={(r,t)∈R2: r∈[0,h(t)),t∈[0,T])}, C和T仅依赖于h0,γ,‖u1,0‖C2([0,∞)), ‖u1,0‖L∞([0,+∞)), ‖u2,0‖C2([0,h0])和‖u3,0‖C2([0,h0]).
证 类似于文献[22, 23],我们首先把自由边界拉直, 将自由边界问题转化为固定区域上的偏微分方程组. 利用压缩映像原理可以证明局部解的存在性和唯一性. 该证明与文献[23]中定理2.1 相似,这里从略.
为了将定理3.1得到的局部解推广到对所有的t>0都成立,我们需要进行下面的估计.
引理3.2 设(u1,u2,u3;h)是问题(1.3)在[0,T0]的一个解,其中T0∈(0,+∞), 则存在不依赖于T0的常数C1,C2和C3使得 0<u1(r,t)≤C1,0≤r<+∞, t∈(0,T0], 0<u2(r,t)≤C2, 0<u3(r,t)≤C3,0≤r<h(t), t∈(0,T0].
证 首先在[0,+∞)×[0,T0]内只要有解存在,则显然有u1≥0,u2≥0和u3≥0. 再对方程(1.3) 应用强极值原理得到 u1(r,t)>0, 0≤r<+∞, t∈(0,T0],u2(r,t)>0, u3(r,t)>0, 0≤r<h(t), t∈(0,T0]. 由问题(1.3)中的第一个方程可知u1(r,t)满足 {u1t−d1Δu1=A−u1(r,t)(λ1u2(r,t)N(r,t)+d),0<r, t>0,[2mm]u1(r,0)=u1,0≥0,r≥0, 利用Phragman-Lindel¨of原则,我们有这样的估计 u1(r,t)≤max{‖u1,0(r)‖L∞(0,∞),Ad}:=C1,0≤r<+∞, t∈(0,T0]. 同理可得 u2(r,t)≤max{‖u2,0(r)‖L∞(0,h0),λ1δC1}:=C2,0≤r≤h(t), t∈(0,T0], u3(r,t)≤max{‖u3,0(r)‖L∞(0,h0),γωC2}:=C3,0≤r≤h(t), t∈(0,T0]. 证毕.
下面的引理说明问题(1.3)的自由边界是严格单调增的.
引理3.3 设(u1,u2,u3;h)是问题(1.3)在[0,T0]的一个解, 则存在不依赖于T0的正常数C4使得 0<h′(t)≤C4, t∈(0,T0].
证 对问题(1.3)中关于u2 的方程应用Hopf引理可导出 (u2)r(h(t),t)<0,0<t≤T0. 因此,由Stefan条件得h′(t)>0, t∈(0,T0].
下面证明对所有的t∈(0,T0]都存在不依赖于T0的常数C4, 使h′(t)≤C4, t∈(0,T0]. 类似文献[17],定义 Γ=ΓM:={(r,t):h(t)−M−1<r<h(t), 0<t≤T0}, 同时构造一个辅助函数 w(r,t):=C2[2M(h(t)−r)−M2(h(t)−r)2], 选取 M:=max{√λ12d2, 4‖u2,0‖C1([0,h0])3C2}, 通过计算容易验证w是u2在Γ中的上解,即w(r,t)≥u2(r,t),(r,t)∈Γ, 进而可以得到 (u2)r(h(t),t)≥wr(h(t),t)=−2MC2, h′(t)=−μ(u2)r(h(t),t)≤2μMC2. 证毕.
接下来,我们给出问题(1.3)解的全局存在性.
引理3.4 对于所有的t∈(0,∞),问题(1.3)的解存在且唯一.
证 设[0,Tmax)是解存在的最大时间范围, 下面将证明Tmax=∞. 反证,假设Tmax为一有限数,则根据引理3.1和引理3.2, 当t∈[0,Tmax),r∈[0,h(t))时,存在不依赖于Tmax的常数 C1,C2,C3和C4,使得下式均成立 0<u1(r,t)≤C1,(r,t)∈[0,+∞)×[0,Tmax), 0<u2(r,t)≤C2, 0<u3(r,t)≤C3, (r,t)∈[0,h(t)]×[0,Tmax), h_0\leq h(t)\leq h_0+C_4t,\,0\leq h'(T)\leq C_4,t\in [0,T_{\max}). 现在我们固定一点\delta_0\in (0,T_{\max}),且令M>T_{\max}, 由标准抛物型偏微分方程理论,我们可以找到一个依赖于\delta_0,M,C_1,C_2,C_3和C_4 的正常数C_5,使得当t\in [\delta_0,T_{\max})时 \| u_{1}(\cdot,t)\| _{ C^{1+\beta}[0,+\infty)},\| u_{2}(\cdot,t)\| _{ C^{2}[0,h(t)]},\| u_{3}(\cdot,t)\| _{C^{2}[0,h(t)]}\leq C_5. 由定理3.1的证明可知,存在一个仅依赖于C_1,C_2,C_3,c_4和C_5的\tau, 使得问题(1.3)在初始时刻T_{\max}-\frac \tau 2 的解可以唯一地延拓到 时刻T_{\max}+\frac \tau 2,这和我们假设T_{\max}是最大存在时间矛盾.
本节考虑自由边界的解的性质,以此说明疾病的传播何时消退或者蔓延,并给出具体的判断标准. 下面分两种情型进行讨论: (I)\,R_0<1,(II)\,R_0>1.
由上节引理3.2,我们得到r=h(t)是单调递增的,因此存在h_\infty\in (0,+\infty]使得 \lim\limits_{t\to \infty} h(t)=h_\infty. 如果h_\infty<\infty且\lim\limits_{t\to \infty}\| u_{2}(\cdot,t)\| _{ C^{2}[0,h(t)]}=\lim\limits_{t\to \infty}\| u_{3}(\cdot,t)\| _{ C^{2}[0,h(t)]}=0,则称疾病消退,这表明疾病扩展的范围有限且带菌者和染病者消失; 反之,则称疾病蔓延. 下面的定理表明在情型(I)下,疾病将消退.
定理4.1 如果R_0(:=\frac{\lambda_1}\delta)<1,则\lim\limits_{t\to \infty}\| u_{2}(\cdot,t)\| _{ C^{2}[0,h(t)]}=0,且有h_\infty<\infty,进而有\lim\limits_{t\to \infty}\| u_{3}(\cdot,t)\| _{ C^{2}[0,h(t)]}=0及u_1(r,t)在[0,\infty)的 任何有界子集上一致地收敛到\frac Ad.
证 注意到u_2满足问题 \begin{eqnarray*} \left\{ \begin{array}{lll} u_{2t}-d_{2}\Delta u_2\leq(\lambda_1-\delta)U_2,\; &0<r<h(t),\; t>0,\\ u_{2r}(0,t)=u_2(r,t)=0,\; &r=h(t),\; t>0,\\ u_2(r,0)\geq 0,\; &0\leq r\leq h_0, \end{array} \right. \end{eqnarray*} 考虑定理条件R_0<1,故\lambda_1-\delta<0,于是当t\to\infty时, 有\| u_{2}(\cdot,t)\| _{ C^{2}[0,h(t)]}\to 0.
下面说明h_\infty<\infty. 事实上,直接计算可以推得 \begin{eqnarray*} \frac{\textrm{d}}{\textrm{d} t}\int_{0}^{h(t)} r^{n-1}u_2(r, t)\textrm{d}r &=&\int_{0}^{h(t)}r^{n-1}u_{2t}(r,t)\textrm{d}r+h'(t)r^{n-1}u_2(h(t),t)\\ &=&\int_{0}^{h(t)}r^{n-1}d_{2}\Delta u_2\textrm{d}r+\int_{0}^{h(t)}r^{n-1} \bigg[\lambda_1 u_2\frac{u_1}{N}-\delta u_2\bigg] \textrm{d}r\\ &=&\int_{0}^{h(t)}d_{2}(r^{n-1}u_{2r}(r,t))_{r}\textrm{d}r+\int_{0}^{h(t)}r^{n-1} \bigg[\lambda_1 u_2\frac{u_1}{N}-\delta u_2\bigg]\textrm{d}r\\ &=&-\frac{d_{2}}{\mu}h^{n-1}h'(t)+\int_{0}^{h(t)}r^{n-1} \bigg[\lambda_1 u_2\frac{u_1}{N}-\delta u_2\bigg]\textrm{d}r. \end{eqnarray*} 从T_0到t(>T_0)积分推得 \begin{eqnarray*} \int_{0}^{h(t)}r^{n-1}u_2(r,t)\textrm{d}r&=&\int^{h(T_{0})}_{0} r^{n-1}u_2(r,T_{0})\textrm{d}r +\frac {d_{2}}{n \mu}h^{n}(T_{0})-\frac {d_{2}}{n\mu}h^{n}(t)\\ &&+\int_{T_0}^t\int_0^{h(s)}r^{n-1} \bigg[\lambda_1 u_2\frac{u_1}{N}-\delta u_2\bigg]\textrm{d}r\textrm{d}s,\ t\geq T_{0}. \end{eqnarray*} 因为此时R_0(:=\frac{\lambda_1}\delta)<1,所以我们有 \lambda_1 u_2\frac{u_1}{N}-\delta u_2<0, \int_{0}^{h(t)}r^{n-1}u-2(r,t)\textrm{d}r\leq \int ^{h(T_{0})} _{0}r^{n-1}u_2(r,T_{0})\textrm{d}r +\frac{d_{2}}{n\mu}h^{n} (T_{0})-\frac {d_{2}}{n\mu}h^{n}(t),\ t\geq T_{0}, 这样我们就证明了h_\infty<\infty.
在此基础上,不难证明\lim\limits_{t\to \infty}\| u_{3}(\cdot,t)\| _{ C^{2}[0,h(t)]}=0, 而u_1(r,t)在[0,\infty)的任何有界子集上一致地收敛到\frac Ad.
现考虑情形(II). 为了后面对u_1,u_2,u_3以及自由边界r=h(t)进行估计, 下面给出比较引理,其证明类似于文献[22,引理3.5].
引理4.2 若T\in (0,+\infty),\overline{h}\in C^{1}([0,T]),\overline{u}_1\in C([0,\infty)\times (-\infty,T]) \cap C^{2,1}((0,\infty)\times (-\infty,T]), \overline{u}_2,\overline u_3\in C([0,\overline{h}(t)]\times [0,T])\cap C^{2,1}((0, \overline h(t))\times (0,T]),且 \begin{eqnarray*} \left\{ \begin{array}{lll} \overline {u}_{1t}-d_{1}\Delta\overline {u}_1\geq A-d \overline u_1,\; 0<r,\; 0<t\leq T,\\ [2mm] \overline {u}_{2t}-d_{2}\Delta\overline {u}_2\geq \lambda_1 \frac{\overline u_1 \overline u_2}{\overline u_1 +\overline u_2} -\delta \overline u_2,\; 0<r<\overline {h}(t),\; 0<t\leq T,\\ [2mm] \overline {u}_{3t}-d_{3}\Delta\overline {u}_3\geq \gamma\overline u_2-\omega\overline u_3,\; 0<r<\overline {h}(t),\ 0<t\leq T,\\ \overline u_{1r}(0,t)=\overline u_{2r}(0,t)=\overline u_{3r}(0,t)=0,0<t\leq T,\\ \overline u_2(r,t)=\overline u_3(r,t)=0,\; r\geq\overline {h}(t),\ 0\leq t\leq T,\\ \overline h'(t)\geq -\mu\overline {u}_{2r}(\overline h(t),t),\ \overline{h}(0):=\overline {h}_0\geq h_{0},\; 0<t\leq T,\\ \overline {u}_1(r,0)\geq u_{1,0}(r),0\leq r,\\ \overline{u}_2(r,0)\geq u_{2,0}(r),\; \overline{u}_3(r,0)\geq u_{3,0}(r),0\leq r\leq h_{0}, \end{array} \right. \end{eqnarray*} 则自由边界问题(1.3)的解(u_1,u_2,u_3;h)满足 u_1(r,t)\leq\overline u_1(r,t),\ h(t)\leq\overline {h}(t),\ (r,t)\in (0,\infty)\times\in(0,T]; u_2(r,t)\leq\overline u_2(r,t),\ u_3(r,t)\leq\overline u_3(r,t),\ (r,t)\in (0,h(t))\times(0,T].
接下来首先给出当h_0 和\mu足够小时,疾病消退的充分条件.
定理4.3 若R_0(:=\frac{\lambda_1}\delta)>1,当h_0\leq \min\{\sqrt{\frac{d_2}{16(\lambda_1-\delta)}},\,\sqrt{\frac{d_3}{16\gamma}}\} 且\mu\leq \frac d{8M}时,则有h_\infty<\infty,其中d=\min \{d_2,\,d_3\},M=\frac 43\max\{\| u_{2,0}\| _{L^\infty},\,\| u_{3,0}\| _{L^\infty}\}.
证 首先我们构造问题(1.3)的上解, \overline u_1(r,t)=M_{1}:=\max\bigg\{\|u_{1,0}\|_{L^\infty},\frac{A}{d}\bigg\}, \overline u_2=\overline u_3=\left\{ \begin{array}{lll} M e^{-\tau t}V(r/\theta(t)),\; 0\leq r\leq \theta(t),\\ 0,\; r> \theta(t), \end{array} \right. \theta (t)=2h_{0}(2- e^{-\tau t}),t\geq 0, V(y)=1-y^{2},0\leq y\leq 1, 其中M和\tau 是待定的正常数.
下面选择合适的M和\tau使得(\overline u_1,\overline u_2,\overline u_3; \theta(t))是 问题(1.3)的上解. 注意到\overline u_1是常数,显然满足 \overline {u}_{1t}- d_{1}\Delta\overline u_1\geq A-d \overline u_1. 又直接计算得 \begin{eqnarray*} & &\overline u_{2t}-d_{2}\Delta\overline u_2- \lambda_1 \frac{\overline u_1\overline u_2}{\overline u_1+\overline u_2}+\delta \overline u_2 \\ &=&M e^{-\tau t}\bigg[-\tau V-r\theta '\theta^{-2}V'-d_{2} \theta^{-2}V''-d_{2}\frac{n-1}{r}\theta^{-1}V'- \bigg(\lambda_1\frac{\overline u_1}{\overline u_1+\overline u_2}-\delta \bigg)V \bigg]\\ &\geq& M e^{-\tau t}\bigg[\frac {d_{2}}{8h_{0}^{2}}-\tau-(\lambda_1-\delta)\bigg], \end{eqnarray*} \overline u_{3t}-d_{3}\Delta \overline u_3-\gamma \overline u_2+\omega \overline u_3\geq M e^{-\tau t}\bigg[\frac {d_{3}}{8h_{0}^{2}}-\tau-\gamma\bigg]. 另外我们有\theta'(t)=2h_{0} \tau e^{-\tau t}和-\mu\overline u_{2r}(\theta(t),t)=2M\mu\theta^{-1}(t)e^{-\tau t}, 而当r\in [0,h_{0}] 时,我们有 \overline {u}_1(r,0)\geq u_{1,0}(r),\ \overline u_2(r, 0)=\overline u_3(r,0)=M\bigg(1-\frac{r^{2}}{4h_{0}^{2}}\bigg)\geq\frac{3}{4}M. 如果我们选择 M=\frac{4}{3}\max\{\|u_{2,0}\|_{L^\infty},\ \|u_{3,0}\|_{L^\infty} \}, \tau=\frac{d}{16h^2_0},并且假设 \mu\leq \frac{d}{8M},\ h_{0}\leq\min \bigg\{\sqrt{\frac{d_2}{16(\lambda_1-\delta)}},\ \sqrt{\frac{d_3}{16\gamma}}\bigg\}, \ d=\min\{d_2,d_3\}, 于是又有 \begin{eqnarray*} \left\{ \begin{array}{lll} \overline {u}_{1t}-d_{1}\Delta\overline {u}_1\geq A-d \overline u_1,\; 0<r,\; 0<t,\\ [2mm] \overline {u}_{2t}-d_{2}\Delta\overline {u}_2\geq \lambda_1 \frac{\overline u_1 \overline u_2}{\overline u_1 +\overline u_2} -\delta \overline u_2,\; 0<r<\theta(t),\; 0<t,\\[2mm] \overline {u}_{3t}-d_{3}\Delta\overline {u}_3\geq \gamma\overline u_2-\omega\overline u_3,\; 0<r<\theta(t),\ 0<t,\\ \overline u_{1r}(0,t)=\overline u_{2r}(0,t)=\overline u_{3r}(0,t)=0,0<t,\\ \overline u_2(r,t)=\overline u_3(r,t)=0,\; r\geq\theta(t),\ 0\leq t,\\ \theta'(t)\geq -\mu\overline {u}_{2r}(\theta(t),t),\ \theta(0)=2 {h}_0> h_{0},\; 0<t,\\ \overline {u}_1(r,0)\geq u_{1,0}(r),0\leq r,\\ \overline{u}_2(r,0)\geq u_{2,0}(r),\; \overline{u}_3(r,0)\geq u_{3,0}(r),0\leq r\leq h_{0}, \end{array} \right. \end{eqnarray*} 因此应用引理4.1得到,当t>0时,h(t)\leq \theta (t), 从而有h_\infty\leq \lim\limits_{t\to \infty} \theta(t)=4h_0<\infty.
引理4.4 若h_{\infty}<\infty,则有 \lim\limits_{t\to+\infty}\|u_2(\cdot,t)\|_ {C([0,h(t)])}=\lim\limits_{t\to+\infty}\|u_3(\cdot,t)\|_{C([0,h(t)])}=0, 并且u_1(r,t) 在[0,\infty) 的任何有界子集上一致地收敛到\frac{A}{d}.
证(反证法) 假设\delta_0:=\limsup\limits_{t\to +\infty}\|u_2(\cdot,t)\|_{C([0,h(t)])}>0. 那么存在序列(r_{k},t_{k} )\in[0,h(t))\times(0, \infty),使得k\to \infty,t_{k}\to \infty,且对所有的k\in{\mathbb N}, u_2(r_{k},t_{k})\geq \delta_0 /2. 由于0\leq r_{k}<h(t)<h_{\infty}<\infty, 于是存在序列\{r_{k}\}的子列使得该子列收敛到r_{0}\in [0,h_{\infty}). 不失一般性,我们仍记为原序列,即r_{k}\to r_{0}\in [0,h_{\infty}), k\to \infty. 定义u^{k}_i(r,t)=u_i(r,t+t_k)\,(i=1,2,3), r\in (0,h(t+t_k)),t\in(-t_{k},\infty). 由抛物方程的正则性可知, 序列\{(u^k_1,u^k_2,u^k_3)\}有一个子序列 \{(u^{k_i}_1,u^{k_i}_2,u^{k_i}_3)\}使当i\to\infty时, (u^{k_i}_1,u^{k_i}_2,u^{k_i}_3)\to (\tilde u_1,\tilde u_2,\tilde u_3),其中 (\tilde u_1,\tilde u_2,\tilde u_3) 满足方程 \begin{eqnarray*} \left\{ \begin{array}{lll} \tilde {u}_{1t}-d_{1}\Delta \tilde {u}_1= A-\tilde u_1 \bigg(\frac {\lambda_1\tilde u-2}{\tilde N}+d\bigg),\; 0<r<h_{\infty},\; t\in(-\infty,+\infty),\\ [3mm] \tilde {u}_{2t}-d_{2}\Delta \tilde {u}_2= \frac {\lambda_1\tilde u_1\tilde u-2}{\tilde N}-\delta \tilde u_2,\; 0<r<h_{\infty},\; t\in(-\infty,+\infty),\\ [2mm] \tilde {u}_{3t}-d_{3}\Delta \tilde {u}_3=\gamma\tilde u_2-\omega\tilde u_3,\; 0<r<h_{\infty},\; t\in(-\infty,+\infty), \end{array} \right. \end{eqnarray*} 其中\tilde {N}=\tilde {u}_1+\tilde u_2+\tilde u_3. 注意到 \tilde u_2(r_{0},0 )\geq \delta_0 /2,因此在[0,h_{\infty})\times(-\infty,\infty)上有 \tilde u_2>0. 由Hopf引理得,\tilde u_{2r}(h_{\infty},0 )\leq -\varepsilon_{0}<0.
另一方面,h(t)是单调递增且有界的. 进而,对任意0<\beta <1, 存在一个依赖于\beta,h_{0},\|u_{2,0}\|C^{2}[0,h_{0}] 以及h_{\infty} 的正常数\tilde C使得 \begin{eqnarray*} \|u_2\|_{C^{(1+\beta)/2,1+\beta}([0,h(t))\times [0,\infty))} +\|h\|_{C^{1+\beta/2}([0,\infty))}\leq \tilde C. \end{eqnarray*} 于是便得h'(t)\to 0,t\to \infty,再由自由边界条件知, 当k\to \infty时,有 u_{2r}(h(t_{k}),t_{k}+0)=(u^{k})_{2r}(h(t_{k}),0)\to \tilde u_{2r}(h_{\infty},0). 于是推出\tilde u_{2r}(h_{\infty},0)=0,这与我们之前得到的\tilde u_{2r}(h_{\infty},0 )\leq -\varepsilon_{0}<0相矛盾. 故假设不成立,\lim\limits_{t\to+\infty}\|u_2(\cdot,t)\|_ {C([0,h(t)])}=0. 由此知,\lim\limits_{t\to+\infty}\|u_3(\cdot,t)\|_{C([0,h(t)])}=0 并且u_1(r,t) 在[0,\infty) 的任何有界子集上一致地收敛到\frac{A}{d}.
最后讨论情型(II). 我们的结果表明当h_0适当大时,疾病总会蔓延.
定理4.5 考虑R_0(:=\frac{\lambda_1}\delta)>1. 若h_{0}>h^{*},则h_{\infty}=\infty,其中\lambda(h^{*})=\frac {\delta (R_0-1)}{2d_2}, 而\lambda(h^{*}) 是问题 \begin{eqnarray*} \left\{ \begin{array}{lll} -\Delta w=\lambda w,\; 0<r<h^*,\\ w_r(0)=0,\; w(h^*)= 0\ \end{array} \right. \end{eqnarray*} 的第一特征值.
证 (反证法) 假设h_{\infty}<+\infty. 由引理4.2,有 \lim\limits_{t\to+\infty}\|u_2(\cdot,t)\|_ {C([0,h(t)])}=\lim\limits_{t\to+\infty}\|u_3(\cdot,t)\|_{C([0,h(t)])}=0 并且u_1(r,t) 在[0,\infty) 的任何有界子集上一致地收敛到\frac{A}{d}. 因此我们可以推出, 当r\in [0,h(t))时,\lim\limits_{t\to+\infty}(\frac {\lambda_1u_1}{N}-\delta)=\lambda_1-\delta=\delta (R_0-1)一致成立, 故存在t^*>0使得当0<r<h(t),t>t^*时, \frac {\lambda_1u_1}{N}=\lambda_1-\delta>\frac {\delta}2 (R_0-1).
注意到 \begin{eqnarray*} \left\{ \begin{array}{lll} u_{2t}-d_{2}\Delta u_2\geq \frac {\delta}2 (R_0-1)u_2,\; 0<r<h_{0},\,t>t^*,\\ [2mm] u_{2r}(0,t)=0,\; u_2(h_{0},t)\geq 0,\ t>t^{*},\\ u_2(r,t^{*})>0,0\leq r<h_{0}, \end{array} \right. \end{eqnarray*} 而当h_{0}>h^{*}时,问题 \begin{eqnarray*} \left\{ \begin{array}{lll} \underline{u}_{2t}-d_{2}\Delta \underline{u}_2=\frac {\delta}2 (R_0-1)\underline{u}_2,\; 0<r<h_0,\,t>t^*,\\[2mm] \underline{u}_{2r}(0,t)=0,\; \underline{u}_2(h_0,t)=0,\ t>t^{*},\\ \underline{u}_2(r,t^{*})=u_2(r,t^{*}),0\leq r<h_{0} \end{array} \right. \end{eqnarray*} 的解是无界的,于是由比较原理得\lim\limits_{t\to +\infty}\|u_2(\cdot,t)\|_{C([0,h(t)])}=+\infty. 这就和 \lim\limits_{t\to +\infty}\|u_2(\cdot,t)\|_{C([0,h(t)])}=0 矛盾. 故h_\infty=\infty.