趋向性描述了微生物在外界条件的刺激下的应激反应,当外界刺激是化学物质时称趋化性. 我们以趋化性为例,讨论微生物应激反应的机理. 本文研究了一类经典的趋化性模型 Keller-Segel 模型[1] 平面行波解的存在性. 此模型的数学描述为
对于上述趋化性模型,很多专家和学者研究了它的行波解. 当灵敏度函数 H(v)=log(v) 时, Othmer 和 Stevens[2] 建立了描述粘细菌因化学物质而聚集的模型. Keller 和 Segel[3] 在实验中观察到某些细菌沿行波带运动, 进而研究了描述这类现象的模型行波解的存在性. 在灵敏度函数 H(v)=log(v) 下, 不考虑化学物质自身的扩散 (u=0), 他们得到了行波解. Rosen[4] 则推广了反馈函数 g 的范围, Keller和 Odell[5] 也将灵敏度函数推广到 H(v)=−v−p. 然而,此时反馈函数的类型依然受到了严格的限制, 并且他们的方法不考虑化学物质的扩散. 后来,Nagai和Ikeda[6] 讨论了 g=−u 时的考虑化学物质扩散性的模型行波解的存在性. Horstmann 和 Stevens[7] 针对几种灵敏度函数和反馈函数提出了一种构造性方法, 以保证行波解的存在性. Schwetlick 和 Hartmut[8] 定性分析了 行波解存在对灵敏度函数和反馈函数的要求,证明了灵敏度函数的奇性是存在 有界行波的必要条件. 当模型考虑到细胞的繁殖和死亡因素时,即使灵敏度函数 非奇性,行波解仍然可能存在. 黎勇[9]在 H(v)=χv1−p 和 g(u,v)=−κ0vαu 时, 就 D=0,D≠0 讨论了行波解的存在性. Xue等[10] 构建了新的细菌趋化性模型,该模型在灵敏度函数非奇的情况 下考虑了细胞间的相关性. 他们证明了这个模型解的全局存在性,从数值方 法和解析方法上证明了行波解的存在性. 陈学勇等[11]讨论了一类 基于趋化性现象的强耦合非线性偏微分方程组,利用相轨分析法,得到了该模型 行波解存在的充分条件和必要条件. 此外,趋化性模型行波解的 稳定性也引起了大量学者的关注. Li 和 Wang[12] 证明了一类双曲-抛物型趋化性模型任意幅度行波 解的非线性稳定性,同时,他们也得到了一类由 Keller-Segel 模型衍变的遵守守恒律的系统在不假设波强度很小的前提下,行波解的存在性 和非线性稳定性[13]. Martin Meyries[14] 研究了奇性灵敏度函数下 带化学物质非线性扩散的 Keller-Segel 模型,得到了行波解的存在性, 在波的一个指数权重邻域内的局部适定性和一定条件下的非线性不稳定性. 最近 Wang[15] 总结了趋化性模型的行波解情况,包括存在性,波速, 渐近衰减率和稳定性等,并且提出了趋化性模型存在的一些开放性问题.
我们研究带动力项的 Othmer-Stevens 模型. 针对当灵敏度函数 H(v)=Λτlog(v) 和 g(u,v)=uγ−vαuβ 时,分别研究模型 (2.1) 在 D=0 和 D≠0 两种不同情况下的行波解的存在性
下面是我们对模型 (2.1) 研究的主要结果.
定理2.1 若常数 α,β,γ 满足条件
注2.1 定理 2.1 中若将 (2.3) 式 改为
定理2.2 当 D=0 时,若常数 α,β,γ 满足
我们将在第三章中证明定理 2.1,而定理 2.2 中行波解的存在性证明则是第四章的内容.
本定理证明分为两个部分,分别对应 D>0,D=0 的不同情形进行讨论. 我们先考虑模型 (2.1) 在 D>0 下行波解的存在性. 我们的证明基于如下事实: 对于问题
首先我们将问题 (2.1) 的求解转换成讨论问题 (3.1) 的形式,然后运用下面的已知结果.
引理3.1[16] 若 f∈C1[0,p∗],f(0)=f(p∗)=0,f′(0)>0,f′(p∗)<0 且对于任意的 p∈(0,p∗),f(p)>0,则存在一个唯一的常数 c∗ 满足
可以证明(详细见附录),若 (B(ξ),S(ξ)) 是问题 (2.2) 的解,则
为了简化 (3.7)式,我们引进新变量 y 和函数 p(y). 令 ε=Dτ,S(ξ)=1mel1ξp(y),ξ=ky, 其中,正常数 m,k,l1 满足 {mλγ−1=Cγ(ξ0),l1(λγ−1)=τcμγ,l1(α+λβ−1)=τcμβ. 于是等式 (3.7) 转化为
由 ξ=ky,k>0 知,ξ 趋向 +∞ 当且仅当 y 趋向 +∞.
可以证明(详见引理 3.3)当条件 (2.3) 成立,且对某正常数 p∗, 如果 p(y) 是下面问题
引理3.2 若 D≠0,(2.3)式 成立, 对满足不等式 (2.4) 的正常数 c,存在一个正常数 p∗ 使得问题 (3.12) 有古典解.
证 若 0<α<1, 由 γ>β 可得
注3.1 当 p>p0 时,方程 (3.11) 不存在有界行波解. 因为 limp→+∞f(p)=+∞,对任意的 p>p0,r′(p)>0, 存在一个正常数 ˆp∈(p0,+∞) 使得 f(ˆp)=0. 可以证明 {d2pdy2+kε(2l1ε+c)dpdy+k2εp(y)[(l21ε+cl1)+pτcl1μγ(y)−pτcl1μβ(y)]=0,y∈R,p(−∞)=ˆp,p(+∞)=p∗ 有严格递减的解. 然后,由函数 S(ξ)=1mel1ξp(y(ξ)) 的定义,可得 limξ→+∞S(ξ)=limξ→+∞1mel1ξp(ξk)=+∞. 因此,当 p>p0 时,方程 (3.11) 不存在有界行波解.
引理3.3 若引理 3.2 中的条件成立,则问题 (2.1) 有正行波解 (u(x,t),v(x,t)) 使得 (B(ξ),S(ξ)) 满足问题 (2.2). 同时, S(ξ) 关于 ξ 单调递增.
证 由上面的推理可知,p(y) 是问题 (3.12) 的解,则 (B(ξ),S(ξ)) 是问题 (2.2) 中方程组的解. 下面我们只需证明 (B(ξ),S(ξ)) 满足问题 (2.2) 的边界条件.
方程 d2pdy2+kε(2l1ε+c)dpdy+f(p)=0 等价于
由此可知,如果引理 3.2 中的条件成立,问题 (2.1) 存在正行波解 (u(x,t),v(x,t)) 使得 (B(\xi),S(\xi)) 是问题 (2.2) 的解.
下面我们用反证法证明 S(\xi) 的单调性. 假设 S(\xi) 关于 \xi 不是单调递增的, 则存在两个实数 \xi _1 < \xi _2 ,但是 S(\xi _1 ) > S(\xi _2) . 令 T =\{( \xi _1 ,\xi _2 ) | \xi _1 ,\xi _2 \in R,S(\xi _1 ) > S( \xi _2 )\}, 则 T 非空. 记 \xi _1^* = \min \limits_ { ( \xi _1 ,\xi _2 ) \in T } \{\xi _1\},\quad \xi _2^* = \max\limits _ { ( \xi _1^* ,\xi _2 ) \in T } \{\xi _2\} . 因为 S(\xi) \in C^1 [\xi _1^* ,\xi _2^*] ,根据闭区间上连续函数的最值性定理, S(\xi) 的最大值可在点 \xi^* 处取得. 由费马定理得 S'(\xi^*)=0 . 记 \xi^* = ky^*, p_0 =mS( \xi^* )e^{-l_1\xi^*} ,则下面三种情形等价.
综上可知,引理 3.3 得证.
现在我们考虑 D=0 的情形. 当 D=0 时,问题 (2.1) 变为
运用推导 (3.7) 式 类似的方法,我们得到
引理3.4n 假设条件 (2.3) 成立,对满足 \frac{\tau \gamma }{\gamma\Lambda - \mu} c^2 =\alpha(1-\alpha)^{\frac{1-\alpha}{\alpha}} 的正常数 c ,则存在常数 S_0 > 0 使得 (3.22) 式 存在单调递增行波解且波速为 c .
证 由 \gamma > \beta > 0,\quad \quad 0 < \alpha < 1,\quad \quad \sigma < 0 和满足条件的正常数 c ,对任意 \bar{X} \geq 0 ,函数 g(\bar X) = \bar{X}(\frac{\eta }{c}\bar X^\beta - \frac{\eta }{c}\bar X^\gamma + \sigma) \leq 0 . 易证函数 g 存在唯一正零点 \bar{X}^*
同理,任意取定 \bar{X}_0 \in (0,\bar{X}^*) ,初值问题 (3.31) 存在唯一的全局严格递增的正解. 事实上,
引理3.5 假设条件 (2.3) 成立,若存在正常数 c 满足 \frac{\tau \gamma}{\gamma\Lambda - \mu} c^2 <\alpha(1-\alpha)^{\frac{1-\alpha}{\alpha}}, 则存在常数 S_0 >0 使得问题 (3.22) 有波速为 c 的单调递增行波解.
证 由 0 < \alpha < 1,\quad \quad \gamma > \beta ,\quad \quad \sigma < 0 和满足上述条件的常数 c , 函数 g\left( {\bar X} \right) = \frac{\eta }{c}\bar X^\beta - \frac{\eta }{c}\bar X^\gamma + \sigma 有两个不同的正零点, 记较小的正零点为 \bar{X}^{*}. 任取 \bar X_0 \in ( 0,\bar {X}^* ) ,考虑初值问题
根据问题 (3.31) 在点 \bar {X} = 0 处的线性化方程,可得
由此可知 (b(x,t),s(x,t)) 是问题 (3.22) 的正行波解,从而 (B(\xi),S(\xi)) 是问题 (3.23) 的解.
S( \xi ) 的单调性证明类似于引理 3.3 中的证明,可根据下列初值问题 (3.43) 解的存在唯一性定理,利用反证法来证明.
综上,引理 3.5 得证.
引理3.6 若 D = 0 ,条件 (2.3) 成立,当常数 c > 0 ,满足条件
由引理 3.4 和引理 3.5,可得引理 3.6
综上,由引理 3.2,3.3 和 3.6,定理 2.1 得证.
值得指出的是,定理 2.1 中的条件 (2.3) 不能被 (2.3') 替代. 下面我们来证明注 2.1,其中引理 3.7 是当 D \neq 0 的情况, 而引理 3.8 是当 D = 0 的情况.
引理3.7 若 D \neq 0 , (2.3') 成立, 则对任意常数 p^* ,问题 (3.12) 无解.
证 由 \alpha = 0 , \beta = \gamma ,问题 (3.12) 中的方程变为
引理3.8 若 (2.3') 成立,问题 (3.22) 不存在正行波解.
证 由 \beta=\gamma ,方程 (3.27) 存在正解 \begin{array}{l} \bar {X}(\xi) = a_0 e^{\sigma \xi},\\xi \in R,\\ \lim\limits_{\xi\to-\infty} \bar {X}(\xi) =+\infty,\quad \quad \lim\limits_{\xi\to+\infty} \bar {X}(\xi) =0. \end{array} 于是
注 2.1 得证.
最后,我们来证明定理 2.2.
引理4.1 若 \beta > 0 > \alpha , \beta=(1-\alpha)\gamma ,则对任意常数 c > 0 ,存在正常数 S_0 使得问题 (3.22) 有单调递增的行波解,波速为 c .
证 由于 \beta > \gamma > 0,\quad c > 0,\quad \eta > 0 , \sigma < 0 ,函数 g(\bar{X})=\frac{\eta }{c}\bar X^\beta - \frac{\eta }{c}\bar X^\gamma + \sigma 有唯一的正零点,记作 \bar{X}^{*} . 取 \bar{X}_{0}\in (0,\bar{X}^{*}) , 初值问题
由此可知 (B(\xi),S(\xi)) 是问题 (3.23) 的解,进而 (b(x,t),s(x,t)) 是问题 (3.22) 的正行波解.
类似于引理 3.5 中的证明可证 S ( \xi ) 的单调性.
综上,定理 2.2 得证.
引理5.1 若 ( B(\xi),S(\xi) ) 是问题 (2.2) 的解,则 \lim\limits_{\xi\to\pm\infty }\frac{dB}{d\xi}(\xi) = 0,\quad \quad \lim\limits_{\xi\to\pm\infty }\frac{dS}{d\xi}(\xi) = 0.
证 (i) 若 \lim\limits_{\xi\to+\infty}\frac{dS}{d\xi} 存在,则它一定为 0. 否则,存在实数 a \neq 0 使得 \lim\limits_{\xi\to+\infty}\frac{dS}{d\xi} = a , \lim\limits_{\xi\to+\infty} \mid S(\xi)\mid = +\infty . 这与条件 \lim\limits_{\xi\to+\infty} S(\xi) = S_0 矛盾,因此若 \lim\limits_{\xi\to+\infty}\frac{dS}{d\xi} 存在,一定为 0.
(ii) 证明极限 \lim\limits_{\xi\to+\infty}\frac{dS}{d\xi} 的存在性. 假设极限 \lim\limits_{\xi\to+\infty}\frac{dS}{d\xi} 不存在,则存在 -\infty \leq a < b \leq+\infty 使得 \limsup_{\xi\rightarrow +\infty}\frac{dS}{d\xi}=b,\hspace{0.5cm} \liminf_{\xi\rightarrow +\infty}\frac{dS}{d\xi}=a. 由扰动引理,存在两个子列 \{\xi_{n}\},\{\eta_{n}\} 满足 \lim\limits_{n\to+\infty} \xi_{n} = + \infty,\quad \lim\limits_{n\to+\infty} \eta_{n} = + \infty, 并且 \lim\limits_{n\to+\infty}\frac{dS}{d\xi}(\xi_{n})= a,\quad \lim\limits_{n\to+\infty}\frac{d^{2}S}{d\xi^2}(\xi_{n}) = 0; \lim\limits_{n\to+\infty}\frac{dS}{d\xi}(\eta_{n})= b,\quad \lim\limits_{n\to+\infty}\frac{d^{2}S}{d\xi^2}(\eta_{n})= 0. 对子列 \{\xi_{n}\} ,因为 S(+\infty)=S_0 >0,B(+\infty)=0 , 由 问题 (2.2) 的第二个方程可得
同理在 (5.1) 式 中,用 \eta_{n} 代替 \xi_{n} , 可证 b = 0 . 这与 a < b 矛盾,所以极限 \lim\limits_{\xi\to+\infty}\frac{dS}{d\xi} 存在且为 0.
同理可证 \lim\limits_{\xi\to-\infty}\frac{dS}{d\xi}=0 . 进一步可得 \lim\limits_{\xi\to\pm\infty}\frac{d^2S}{d\xi^2}=0.
(iii) 类似于 \lim\limits_{\xi\to\pm\infty}\frac{dS}{d\xi}=0,\quad \quad \lim\limits_{\xi\to\pm\infty}\frac{d^2S}{d\xi^2}=0 的证明,从问题 (2.2) 的第一个方程,可得 \lim\limits_{\xi\to\pm\infty}\frac{dB}{d\xi}=0.
引理 5.1 得证.