留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

具有非线性扩散的趋化模型弱解的有界性

上一篇

下一篇

张婕燕, 辛巧, 穆春来. 具有非线性扩散的趋化模型弱解的有界性[J]. 西南大学学报(自然科学版), 2021, 43(5): 88-95. doi: 10.13718/j.cnki.xdzk.2021.05.012
引用本文: 张婕燕, 辛巧, 穆春来. 具有非线性扩散的趋化模型弱解的有界性[J]. 西南大学学报(自然科学版), 2021, 43(5): 88-95. doi: 10.13718/j.cnki.xdzk.2021.05.012
ZHANG Jie-yan, XIN Qiao, MU Chun-lai. Boundedness of the Weak Solutions to the Chemotaxis Model with Nonlinear Diffusion[J]. Journal of Southwest University Natural Science Edition, 2021, 43(5): 88-95. doi: 10.13718/j.cnki.xdzk.2021.05.012
Citation: ZHANG Jie-yan, XIN Qiao, MU Chun-lai. Boundedness of the Weak Solutions to the Chemotaxis Model with Nonlinear Diffusion[J]. Journal of Southwest University Natural Science Edition, 2021, 43(5): 88-95. doi: 10.13718/j.cnki.xdzk.2021.05.012

具有非线性扩散的趋化模型弱解的有界性

  • 基金项目: 新疆维吾尔自治区自然科学基金面上项目(2018D01C004)
详细信息
    作者简介:

    张婕燕,硕士研究生,主要从事偏微分方程研究 .

  • 中图分类号: O175.26

Boundedness of the Weak Solutions to the Chemotaxis Model with Nonlinear Diffusion

  • 摘要: 研究间接信号产出的非线性扩散-趋化聚集模型,并给出具有非负且适当光滑的初值条件和相应的扩散系数. 通过能量估计方法,证明生物趋化模型弱解的全局有界性及唯一性.
  • 加载中
  • [1] KELLER E F, SEGEL L A. Initiation of Slime Mold Aggregation Viewed asan Instability [J]. Journal of Theoretical Biology, 1970, 26(3): 399-415. doi: 10.1016/0022-5193(70)90092-5
    [2] doi: http://www.ams.org/mathscinet-getitem?mr=2052525 PAINTER K J, HILLEN T. Volume-Filling and Quorum-Sensing in Models for Chemosensitive Movement [J]. Canadian Applied Mathematics Quarterly, 2002, 10(4): 501-543.
    [3] 张颖, 赵志新, 张优佳, 等. 一类拟线性抛物-椭圆趋化增长系统解的全局有界性[J]. 西南师范大学学报(自然科学版), 2019, 44(1): 34-39. doi: https://www.cnki.com.cn/Article/CJFDTOTAL-XNZK201901008.htm
    [4] CIEŚLAK T. Global Existence of Solutions to a Chemotaxis System with Volume Filling Effect [J]. Colloquium Mathematicum, 2008, 111(1): 117-134. doi: 10.4064/cm111-1-11
    [5] LIUD M, TAO Y S. Global Boundedness in a Fully Parabolic Attraction-Repulsion Chemotaxis Model [J]. Mathematical Methods in the Applied Sciences, 2015, 38(12): 2537-2546. doi: 10.1002/mma.3240
    [6] doi: http://www.ams.org/mathscinet-getitem?mr=1893940 OSAKI K, YAGI A. Finite Dimensional Attractors for One-Dimensional Keller-Segel Equations [J]. Funkcialaj Ekvacioj, 2001, 44(3): 441-469.
    [7] NAGAI T. Blow-up of Radially Symmtetric Solutions to a Chemotaxis Systems [J]. Advances in Applied Mathematics, 1995, 5(2): 581-601.
    [8] NAGAI T. Blowup of Nonradial Solutions to Parabolic-Elliptic Systems Modeling Chemotaxis in Two-Dimensional Domains [J]. Journal of Inequalities and Applications, 2001, 2001(1): 970292. doi: 10.1155/S1025583401000042
    [9] doi: http://www.ams.org/mathscinet-getitem?mr=1610709 NAGAI T, SENBA T, YOSHIDA K. Application of the Trudinger-Moser Inequality to a Parabolic System of Chemotaxis [J]. Funkcialaj Ekvacioj, 1997, 40(3): 411-433.
    [10] CORRIAS L, PERTHAME B, ZAAG H. Global Solutions ofsome Chemotaxis and Angiogenesis Systems in High Space Dimensions [J]. Milan Journal of Mathematics, 2004, 72(1): 1-28. doi: 10.1007/s00032-003-0026-x
    [11] KOWALCZYK R, SZYMAŃSKA Z. On the Global Existence of Solutions to an Aggregation Model [J]. Journal of Mathematical Analysis and Applications, 2008, 343(1): 379-398. doi: 10.1016/j.jmaa.2008.01.005
    [12] STROHM S, TYSON R C, POWELL J A. Pattern Formation in a Model for Mountain Pine Beetle Dispersal: Linking Model Predictions to Data [J]. Bulletin of Mathematical Biology, 2013, 75(10): 1778-1797. doi: 10.1007/s11538-013-9868-8
    [13] 王笑丹, 陶有山. 具间接信号产出的二维趋向性模型古典解的整体存在性和有界性[J]. 东华大学学报(自然科学版), 2016, 42(6): 922-930. doi: 10.3969/j.issn.1671-0444.2016.06.023
    [14] TELLO J I, WINKLER M. A Chemotaxis System with Logistic Source [J]. Communications in Partial Differential Equations, 2007, 32(6): 849-877. doi: 10.1080/03605300701319003
    [15] doi: http://www.sciencedirect.com/science/article/pii/S0898122118300580 QIUS Y, MU C L, WANG L C. Boundedness in the Higher-Dimensional Quasilinear Chemotaxis-Growth System with Indirect Attractant Production [J]. Computers & Mathematics With Applications, 2018, 75(9): 3213-3223.
    [16] HORSTMANN D, WINKLER M. BoundednessVs. Blow-up in a Chemotaxis System [J]. Journal of Differential Equations, 2005, 215(1): 52-107. doi: 10.1016/j.jde.2004.10.022
    [17] LAURENCOT P, WRZOSEK D. A Chemotaxis Model with Threshold Density and Degenerate Diffusion [M] // CHIPOT M, ESCHER J. Nonlinear Elliptic and Parabolic Problems. Berlin: Birkhäuser Basel, 2006.
    [18] KOWALCZYK R. Preventing Blow-up in a Chemotaxis Model [J]. Journal of Mathematical Analysis and Applications, 2005, 305(2): 566-588. doi: 10.1016/j.jmaa.2004.12.009
  • 加载中
计量
  • 文章访问数:  788
  • HTML全文浏览数:  788
  • PDF下载数:  90
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-05-28
  • 刊出日期:  2021-05-20

具有非线性扩散的趋化模型弱解的有界性

    作者简介: 张婕燕,硕士研究生,主要从事偏微分方程研究
  • 1. 伊犁师范大学 数学与统计学院,新疆 伊宁 835000
  • 2. 重庆大学 数学与统计学院,重庆 401331
基金项目:  新疆维吾尔自治区自然科学基金面上项目(2018D01C004)

摘要: 研究间接信号产出的非线性扩散-趋化聚集模型,并给出具有非负且适当光滑的初值条件和相应的扩散系数. 通过能量估计方法,证明生物趋化模型弱解的全局有界性及唯一性.

English Abstract

  • 趋化行为是细胞沿化学物质浓度高的方向作的一种定向运动,其研究在胚胎发育、血管生成及种群动力学行为等方面具有重要意义. 文献[1]建立了经典的生物趋化模型,主要用于刻画细胞对趋化—交叉扩散产生的奇异性反应的集体行为. 在过去的几十年中,经典的Keller-Segel模型解的存在性、唯一性、爆破行为及渐近行为等问题的研究得到了越来越多学者的关注,取得了大量的研究成果[2-6]. 最简单的Keller-Segel模型(记作KS模型)如下:

    KS模型解的性质(全局有界性或爆破行为)强依赖于空间维数:在空间维数d=1时,文献[6]证明了KS模型解的全局有界性;在空间维数d=2时,文献[7-9]证明了存在临界指数现象,若细胞的初始质量$ \int_{\varOmega} u(x, 0) \mathrm{d} x<4 {\rm{ \mathsf{ π}}}$,则KS模型的解全局有界,若$\int_{\varOmega} u(x, 0) \mathrm{d} x>4 {\rm{ \mathsf{ π}}} $,则KS模型存在解在有限时间爆破;在空间维数d≥3时,文献[10]证明了当初值$\left\|u_{0}\right\|_{L^{\frac{2}{d}}} $足够小时,KS模型的解全局有界,$\left\|u_{0}\right\|_{L^{\frac{2}{d}}} $足够大时,解在有限时间内爆破. 事实上,关于KS模型的研究到目前为止并不完善,还有很多需要解决的问题. 考虑到KS模型的临界指数现象,文献[11]考虑具有非线性扩散的KS模型(d≥2),在扩散函数足够大时,证明了生物趋化系统弱解的全局有界性,其研究的非线性扩散的生物趋化模型如下:

    KS模型中的化学物质由细胞直接产生,不同于该化学物质产生模式,文献[12]提出间接信号产出的生物趋化模型,用于刻画山松甲壳虫的扩散和聚集行为. 关于间接信号产生的生物趋化模型,主要集中研究带有Logistic源的情形,具体结果可以参考文献[13-15]. 特别需要提到的是,文献[15]考虑具有拟线性扩散的生物模型,其研究的拟线性扩散的生物趋化模型如下[15]

    文献[15]证明了当d≥3,μ>0,δ>0和τ>0且$p>1-\frac{4}{d} $时,生物趋化模型整体解的一致有界性. 在上述研究的启发下,本文主要研究如下具有非线性扩散的间接信号产生的生物趋化模型(不带有Logistic源):

    其中:$\varOmega \subset \mathbb{R}^{d}$(空间维数d≥2)表示光滑有界区域,υ表示边界∂Ω单位外法线方向;初值条件(u0v0w0)满足u0≠0,非负且适当光滑;扩散系数D(u)≥D0up$ p>1-\frac{4}{d}$D0>0;u(xt)表示飞行的甲壳虫的密度,v(xt)表示甲壳虫信息素的浓度,w(xt)表示筑巢的甲壳虫密度,不同于经典的KS模型,甲虫信息素是由筑巢的甲壳虫产生,致使飞行中的甲壳虫产生趋化行为,该甲虫信息素对筑巢的甲壳虫无任何作用. 本文通过类似于文献[11]中的方法来考虑不带有Logistic源的情形,主要得到以下结论:

    定理1  设$\varOmega \subset \mathbb{R}^{d} $(空间维数d≥2) 为具有C2边界的有界开区域,初值条件$ u_{0} \in W^{1, \bar{p}}(\varOmega), v_{0} \in W^{1, \infty}(\varOmega) $$ w_{0} \in L^{\infty}(\varOmega)$,其中pd. 扩散系数D(u)为具有C2光滑性的非减函数,且对任意的u≥0,有D(u)≥D0up,其中$D_{0}>0, p>2-\frac{4}{d} $. 那么,生物趋化模型(2)存在非负弱解(uvw)且一致有界.

    注1  关于生物趋化模型(2)的弱解的定义可参考文献[11]中定义1类似定义,本文略.

    注2  当山松甲壳虫的随机布朗运动足够强烈时,即山松甲壳虫的扩散强度超过甲虫信息素对其的趋化吸引作用,那么这个松林中山松甲壳虫不会存在大量聚集现象,不会对加拿大短叶松等松林造成了巨大的破坏. 如何增大山松甲壳虫的扩散运动,可能对于未来松林的保护具有重要的意义.

    定理2  设D(u)是具有C2光滑的非减函数,趋化模型(2)存在弱解且vL(0,TW2,∞(Ω)),则趋化模型(2)的解必唯一.

  • 引理1[16]  设v0W1,∞,对所有t∈(0,T),都有‖w(t)‖Lq(Ω)C,那么对任意t∈(0,T)和q∈[1,d)可得

    其中

    q=d,则对于任意的s<+∞,都有(3)式成立;若qd,则有‖v(t)‖W1,∞C.

    此外,考虑到扩散系数D(u)可能具有奇异性,先考虑其正则性的趋化模型:

    其中Dσ(u)=D(u+σ)且σ>0,利用文献[17]中Amann理论容易证明其局部适定性. 接下来,先证明趋化系统(5)的解的全局适定性,主要有如下引理.

    引理2  设d≥2,$p>2-\frac{4}{d} $,若对∀t∈(0,T),有

    其中$ q_{0} \in\left[\max (1, p), \min \left(p+\frac{2}{d}, d\right)\right]$C是常数,则有

    (i) $\|w(t)\|_{L^{q_{0}}} \leqslant C $$\|\nabla v(t)\|_{L^{s_{0}}} \leqslant C $,其中$s_{0}<\frac{d q_{0}}{d-q_{0}} $.

    (ii) 若d>2且pd-2,则对于任意的$ q<\frac{d}{d-2} p$,都有$\left\|u_{\varepsilon}(t)\right\|_{L^{q}(\varOmega)} \leqslant C(q) $.

    (iii) 若d=2或pd-2,则对于任意的q<∞,都有$\left\|u_{\varepsilon}(t)\right\|_{L^{q}(\varOmega)} \leqslant C(q) $.

      (i)对趋化模型(5)的第三个方程两边同乘wq0-1,然后在Ω上积分可得

    进一步,对(6)式的右端的最后一项用Young不等式可得

    因为uLq0范数有界,则联合(6)和(7)式可知

    其中C是常数,即可知

    最后,利用引理1可知$ \|\nabla v(t)\|_{L^{s_{0}}} \leqslant C$成立.

    (ii) 对于任意的qq0,设uε=u+ε,其中ε∈(0,σ],对趋化模型(5)的第一个方程两边同乘uεq-1,然后在Ω上积分可得

    因为

    那么由Young不等式可知

    另一方面,对于∀ε>0,有

    进而可知

    现在,令$g=u_{\varepsilon}^{\frac{p+q}{2}} $,则有

    首先假设qp+q0,则对(9)式右边的积分用Hölder不等式可知

    进一步,取$ q_{0}<q<p+\frac{2}{d} q_{0}$$ 1 \leqslant s<\min \left(\frac{d q_{0}}{d-q_{0}}, \frac{2 q_{0}}{p-q+q_{0}}\right)$,因为d≥3,可设$s<\frac{2 d}{d-2} $. 又因为$\|\nabla v\|_{L^{s}(\varOmega)} \leqslant C $,则由Gagliardo-Nirenberg不等式可知

    其中b∈(0,1). 由于-ΔDL2(Ω)的同构映射(证明见文献[16]),其中

    那么,由Young不等式可知

    整理可得

    其中C是常数.

    接下来,对趋化系统(5)的第二个方程两边同乘-Δv,然后在Ω上积分可得

    由Young不等式可知

    最后,整理可得

    对趋化模型(5)的第三个方程两边同乘w,然后在Ω上积分可得

    由Young不等式可得

    整理得

    联合不等式(11)-(13)可知

    接下来,关键之处在于处理$\int_{\varOmega} u_{\varepsilon }^{2} \mathrm{~d} x $. 下面分为两种情况进行讨论:

    1) 在p+q>2时,由Poincare不等式得

    2) 若p+q≤2,选取$n \in\left(\frac{2 d}{p+q}-d, \min \left(\frac{2 q_{0}}{p+q}, \frac{4}{p+q}\right)\right) $. 因为d≥3,可设$n<\frac{2 d}{d-2} $. 则由Gagliardo-Nirenberg不等式可知

    其中$ a=\frac{\frac{d}{n}-\frac{d(p+q)}{4}}{1-\frac{d}{2}+\frac{d}{n}}<\frac{p+q}{2}$等价于$n>\frac{2 d}{p+q}-d $.

    进一步,因为$ \|g\|_{L^{n}(\varOmega)}$有界和$\frac{4 a}{p+q}<2 $,则由Poincare不等式可得

    综上所述可得

    现在,联合不等式(14)和(15)可得

    其中C(q)表示仅依赖于q的常数. 此外,还有

    这样就可以得到

    其中C$\mathop C\limits^ \wedge $是常数. 那么,对所有的t∈(0,T),由Gronwall不等式可知

    其中C表示仅依赖于q的常数与时间无关.

    综合上述过程,在$q_{0} \in\left[\max (1, p), \min \left(p+\frac{2}{d}, d\right)\right] $$p>2-\frac{4}{d} $满足时,则对

    可知(17)式成立.

    q1d,则由引理1可知,对于任意的l<+∞,都有$ \|\nabla v(t)\|_{L^{l}} \leqslant C$,易知(9)式右边有界和(16)式成立,那么对所有的q=q2p+q1,亦可知(16)式成立. 然后,用q2替换(10)中的q0,对所有的q<2p+q1,都有(17)式成立. 重复上述过程k次,可知对于任意大的k,对所有的qkp+q1,使得(17)式成立. 因此,若q1d,对所有的q<+∞,可得(17)式成立.

    如果q1d,用q1替换(10)中q0,则有$ $都有(17)式成立,依次类推,可知对所有$ $,可得(17)式成立. 最后,若存在某个k使得qkd,那么根据上述证明可知,对所有的q<+∞,有(17)式成立. 如若不然,则对所有k都有qkd成立. 则序列$ $必存在上极限$ $. 因此,对每个qd都有(17)式成立. 因为qd,所以pd-2.

    (iii) 从上述证明可知,若d=2或pd-2,对所有的q<+∞,可得(17)式成立.

    接下来,可以利用文献[11]引理3证明得到引理3,利用文献[18]引理4.1证明得到引理4.

    引理3[11]  假设定理1中所有条件成立,则

    引理4[18]  设ΩRd(d≥1)是具有C1边界的有界的开区域,假设m>0,um>0,对∀uum,使得D(u)≥m. 如果0<tT$\|\nabla v(t)\|_{L^{\infty}} \leqslant C_{1} $,对0<tT,可得

    其中C1C2C3都是正数.

    最后,类似于文献[11]的证明方法,由σ的任意性,可知解的收敛性. 此外,上述推导出的所有估计值均不依赖于T,因此趋化系统(2)弱解有界性是全局的,可知定理1成立.

  • 本节利用能量方法来证明定理2即生物趋化模型(2)解的唯一性,主要证明过程如下.

    定理2证明  设方程(2)有两个满足初边值条件的解(u1v1w1)和(u2v2w2),现令(uvw)=(u1-u2v1-v2w1-w2)和ψ满足

    则可知

    其中X(u)=D(u). 对模型(21)的前两个方程分别乘ψv并在ΩT上积分得

    进而可得

    接下来,对(23)的右边用Cauchy不等式得

    此外可知

    其中

    进一步由方程(20)可知$\|\nabla \psi(0)\|_{L^{2}(\varOmega)}=0 $. 进而可得:

    最后可得

    其中C是常数. 因此,方程(24)等价于

    则由Gronwall's不等式可知,对所有t∈(0,T)且T>0,ψ(t)=u(t)=v(t)=0,进而有w(t)=0,即唯一性得证.

参考文献 (18)

目录

/

返回文章
返回