The roles of nonlinear diffusion, haptotaxis and ECM remodelling in determining the global solvability of a cancer invasion model

Chunhua Jin

doi:10.1017/prm.2024.70

The roles of nonlinear diffusion, haptotaxis and ECM remodelling in determining the global solvability of a cancer invasion model

Part of: Parabolic equations and systems Equations and systems of special type Physiological, cellular and medical topics

Published online by Cambridge University Press: 13 September 2024

Chunhua Jin

Show author details

Chunhua Jin*: Affiliation:
School of Mathematical Sciences, South China Normal University, Guangzhou, 510631, China (jinchhua@126.com)
*: *Corresponding author

Article contents

Abstract
Introduction
Preliminaries
Uniform energy estimations and global solvability of chemotaxis-haptotaxis system
Global solution to the haptotaxis-only system
Statements and declarations
References

Rights & Permissions

Abstract

In this paper, we consider the following PDE-ODE system modelling cancer invasion with slow diffusion and ECM remodelling,

\[ \begin{cases} u_t=\Delta u^m-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega), \\ v_t=\Delta v+u-v, \\ \omega_t={-}v\omega+\eta \omega(1-u-\omega). \end{cases} \]

For the special case $\eta =0$

, fruitful results have been achieved since Tao and Winkler's work in 2011. However, there is no any progress for the general case $\eta >0$

in the past ten years. In this paper, we analysed some commonly used research methods when $\eta =0$

, and found that these methods are completely unsuitable for situations where $\eta >0$

. By introducing some new forms of functionals, we reconstruct the relationship between the haptotactic term and the nonlinear diffusion term, and ultimately prove the global existence of weak solutions. This result improves and perfects a series of works previously presented in the literature.

Keywords

haptotaxis nonlinear diffusion ECM remodelling global solution

MSC classification

Primary: 35M13: Initial-boundary value problems for equations of mixed type

Secondary: 35K65: Degenerate parabolic equations 92C17: Cell movement (chemotaxis, etc.)

Type: Research Article
Information: Proceedings of the Royal Society of Edinburgh Section A: Mathematics , First View , pp. 1 - 32

DOI: https://doi.org/10.1017/prm.2024.70 [Opens in a new window]
Copyright: Copyright © The Author(s), 2024. Published by Cambridge University Press on behalf of The Royal Society of Edinburgh

1. Introduction

As a major disease threatening human life safety, cancer has always been the focus of social attention. The combination of mathematical models and various medical data will help people turn cancer research into a quantitative and predictable science. The mathematical modelling of tumour growth and its theoretical research have always been the concern of biologists and mathematicians. Generally speaking, the growth of a tumour usually goes through two stages: vascular phase and vascularization. That is, if there are no blood vessels to provide enough nutrients, the tumour will stop growing when it grows to a certain size. In order to continue to grow, some tumours will secrete a chemical substance, which is called urokinase Plasminogen Activator (uPA), to recruit vascular factors, and to build a vascular network around themselves to meet the growing nutritional needs. At the same time, blood vessels also provide a way for tumours to travel to other parts of the body, which is called metastasis. In recent years, some reaction-diffusion-taxis models are proposed to characterize the process of tumour growth and invasion [Reference Byrne, Chaplain, J.Pettet and Mcelwain1, Reference Chaplain and Anderson3, Reference Chaplain and Lolas4, Reference Stinner and Surulescu17]. In particular, the following model proposed by Chaplain and Lolas [Reference Chaplain and Lolas4] has attracted extensive attention of mathematicians in recent years.

(1.1)\begin{equation} \left\{ \begin{aligned} & u_t=\nabla\cdot(D(u)\nabla u)-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega), \\ & v_t=\Delta v+u-v, \\ & \omega_t={-}v\omega+\eta \omega(1-u-\omega), \end{aligned}\right. \end{equation}

which describes the invasion and diffusion process of solid tumours during the vascular growth stage. In this system, $u$, $v$, $\omega$ represent the cancer cell density, urokinase Plasminogen Activator (uPA) concentration, and the extracellular matrix (ECM) density respectively. $D(u)$ denotes the mobility of cells, the positive constants $\chi$, $\xi$ denote the uPA-mediated chemotaxis, ECM-mediated haptotaxis coefficients respectively, $\mu u(1-u-\omega )$ denotes the proliferation or death of cancer cells. In the second equation, $+u$ represents the production of uPA by cancer cells, $-v$ represents decay of uPA. As for the ECM density $\omega$, it is known that it does not diffuse, therefore, one can omit the random motion. $-v\omega$ denotes the degradation of ECM, $\eta \omega (1-u-\omega )$ represents the remodelling of ECM components.

In the past ten years, this model has been widely studied. When $D(u)\equiv 1$, $\eta =0$, we refer to [Reference Cao5, Reference Marciniak-Czochra and Ptashnyk13, Reference Tao and Wang18, Reference Walker and Webb23] for the study of global existence of bounded solutions, and refer to [Reference Hillen, Painter and Winkler6, Reference Tao and Winkler20, Reference Wang and Ke22] for the study of large time behaviour. While, if the remodelling of ECM is considered, the calculation of the haptotaxis term will bring some essential difficulties. Therefore, the following transformation is introduced to avoid the estimation of the haptotaxis term

\[ \rho=ue^{-\xi\omega}, \]

and the first equation of (1.1) is transformed into

(1.2)\begin{align} e^{\xi\omega}\rho_t-{\rm div}(e^{\xi\omega}\nabla\rho)={-}\chi{\rm div}(\rho e^{\xi\omega}\nabla v) +\xi\rho v\omega e^{\xi\omega}+\rho e^{\xi\omega}(1-\omega-\rho e^{\xi\omega})(\mu -\xi\eta\omega). \end{align}

Noticing that $\omega$ is bounded, then $\|u\|_{L^p}\sim \|\rho \|_{L^p}$. Using this method, the global existence and boundedness of classical solution for any initial datum in two dimensional space is proved [Reference Jin10, Reference Pang and Wang15]. While in three dimensional space, only a small global classical solution is established [Reference Jin11, Reference Pang and Wang16]. In fact, even for the haptotaxis-only system, the global bounded solution for any $\mu,\, \eta >0$ in dimension 3 is still open.

Considering the mechanism of avoiding crowding between individual cells, the mobility of cells should be related to density, so a model with nonlinear diffusion is also very practical, and the most representative one is the porous medium diffusion model. For example, in [Reference Byrne, Chaplain, J.Pettet and Mcelwain1, Reference Byrne, Chaplain, J.Pettet and Mcelwain2], the authors use porous medium diffusion to describe the process of trophoblast cells invading the uterine tissue, that is,

(1.3)\begin{equation} \left\{ \begin{aligned} & u_t=D_u\nabla\cdot(u^{m-1}\nabla u)-\xi\nabla\cdot(u\nabla v)+k_1 u(1-u-\omega), \\ & n_t=\Delta n+k_2nu(1-u)-k_3nv, \\ & v_t=D_v\Delta v+k_4n\omega-k_3nv, \\ & \omega_t={-}\delta n\omega+\eta\omega(1-u-\omega). \end{aligned}\right. \end{equation}

In 2011, Tao and Winkler [Reference Tao and Winkler19] first studied the global solvability of (1.1) by ignoring the remodelling of ECM, that is, for

(1.4)\begin{equation} \left\{ \begin{aligned} & u_t=\Delta u^m-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega), \\ & v_t=\Delta v+u-v, \\ & \omega_t={-}v\omega, \end{aligned}\right. \end{equation}

they proved the global existence of weak solutions when

\[ m> m^*=\begin{cases} \dfrac{2N^2+4N-4}{N(N+4)}, & \text{if} \ N\le 8, \\ \dfrac{2N^2+3N+2-\sqrt{8N(N+1)}}{N(N+1)}, & \text{if} \ N\ge 9. \end{cases} \]

After that, the global existence of weak solutions to this problem began to be widely studied, see for example [Reference Jin7, Reference Li and Lankeit12, Reference Wang21, Reference Zheng24] etc. For this case, from the third equation of (1.4), one observed that

(1.5)\begin{equation} \omega=\omega_0 e^{-\int_0^t v(x,s){\rm d}s}. \end{equation}

At this time, there are two commonly used methods to assist in dealing with the haptotactic term.

Method I: A direct calculation from (1.5) leads to

\[ |\nabla \omega(x,t)|^2\le2|\nabla \omega_0(x)|^2e^{{-}2\int_0^t v(x,s){\rm d}s} +2|w_0(x)|^2e^{{-}2\int_0^tv(x,s){\rm d}s}\Big|\int_0^t\nabla v(x,s){\rm d}s\Big|^2. \]

Therefore,

\begin{align*} \int_\Omega|\nabla \omega(x,t)|^2{\rm d}x & \le C+C\int_\Omega e^{{-}2\int_0^tv(x,s){\rm d}s}\Big|\int_0^t\nabla v(x,s){\rm d}s\Big|^2{\rm d}x\\ & = C-\frac{C}{2}\int_\Omega\nabla e^{{-}2\int_0^tv(x,s){\rm d}s}\cdot\Big(\int_0^t\nabla v(x,s){\rm d}s\Big){\rm d}x\\ & = C+\frac{C}{2}\int_\Omega e^{{-}2\int_0^tv(x,s){\rm d}s}\cdot\Big(\int_0^t\Delta v(x,s){\rm d}s\Big){\rm d}x\\ & = C+\frac{C}{2}\int_\Omega e^{{-}2\int_0^tv(x,s){\rm d}s}\cdot\Big(\int_0^t(v_t+v-u){\rm d}s\Big){\rm d}x\\ & \le C+\frac{C}{2}\int_\Omega e^{{-}2\int_0^tv(x,s){\rm d}s}\cdot\Big(v(x,t)+\int_0^tv(x,s){\rm d}s\Big){\rm d}x \\ & \le \tilde C+\frac{C}{2}\int_\Omega v(x,t){\rm d}x. \end{align*}

Similarly, one can use iteration to sequentially obtain $\|\nabla \omega (\cdot,\,t)\|_{L^4}$, $\|\nabla \omega (\cdot,\,t)\|_{L^6}$, $\|\nabla \omega (\cdot,\,t)\|_{L^{2n}}\cdots$.

Method II: A direct calculation from (1.5) leads to

\begin{align*} -\Delta\omega(x,t)& ={-}\Delta\omega_0e^{-\int_0^t v(x,s){\rm d}s}+2e^{-\int_0^t v(x,s){\rm d}s}\nabla \omega_0\cdot\int_0^t \nabla v(x,s){\rm d}s\nonumber \\ & \quad-\omega_0e^{-\int_0^t v(x,s){\rm d}s}\left|\int_0^t \nabla v(x,s){\rm d}s\right|^2+\omega_0 e^{-\int_0^t v(x,s){\rm d}s}\int_0^t \Delta v(x,s){\rm d}s \\ & ={-}\Delta\omega_0e^{-\int_0^t v(x,s){\rm d}s}-\omega_0e^{-\int_0^t v(x,s){\rm d}s}\left(\frac{\nabla\omega_0}{\omega_0}-\int_0^t \nabla v(x,s){\rm d}s\right)^2\nonumber \\ & \quad +\frac{|\nabla\omega_0|^2}{\omega_0}e^{-\int_0^t v(x,s){\rm d}s}+\omega_0 e^{-\int_0^t v(x,s){\rm d}s}\int_0^t (v_t+v-u){\rm d}s \\ & \le \|\Delta \omega_0\|_{L^\infty}+4\|\nabla\sqrt{\omega_{0}}\|_{L^\infty}^2+\|\omega_{0}\|_{L^\infty}(1+v(x,t)) \\ & \le\|\omega_0\|_{L^\infty}v(x,t)+K \end{align*}

where $K$ only depends on $\omega _0$ [Reference Tao and Winkler19]. Therefore, the above two methods can make the haptotactic term easy to handle, however, the two methods are obviously not suitable for situations where $\eta >0$. In fact, if the remodelling of ECM is considered, that is, the third equation of (1.4) is replaced with

\[ \omega_t={-}v\omega+\eta \omega(1-u-\omega), \]

the problem becomes much more complex since the regularity of $\omega$ completely depends on the regularity of $u$. Actually, due to that the ODE has no regularization effect, the regularity of $\omega$ is greatly reduced by $-\eta u\omega$ in remodelling term. Therefore, although a lot of results for the case $\eta =0$ have been achieved in the past ten years, there is no any progress in the case of $\eta >0$.

In the present paper, we consider the initial and boundary value problem for the system (1.1) with nonlinear diffusion, that is

(1.6)\begin{equation} \begin{cases} u_t=\Delta u^m-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega), & \text{in}\ Q, \\ v_t=\Delta v+u-v, & \text{in}\ Q, \\ \omega_t={-}v\omega+\eta \omega(1-u-\omega), & \text{in}\ Q, \\ \left.\displaystyle\frac{\partial u^m}{\partial{\bf n}}-\chi u\frac{\partial v}{\partial{\bf n}}-\xi u\frac{\partial \omega}{\partial{\bf n}} \right|_{\partial\Omega}=0, \left.\frac{\partial v}{\partial{\bf n}}\right|_{\partial\Omega}=0, & \\ u(x,0)=u_0(x), \quad v(x,0)=v_0(x), \quad \omega(x,0)=\omega_0(x), & x\in\Omega, \end{cases} \end{equation}

where $Q=\Omega \times \mathbb {R}^+$, $m>1$, $\Omega \subset \mathbb {R}^N$ ($N\ge 2$) is a smooth bounded domain. $\chi$, $\xi$, $\mu$, $\eta$ are positive constants.

According to previous research experience for some chemotaxis models, the larger $m$ is, the easier it is to prove the uniform boundedness of the solution, but this experience seems to be invalid for the current model. On the one hand, due to the existence of remodelling term $\eta \omega (1-u-\omega )$, the methods (Method I and Method II) used in the study of equations (1.4) is completely invalid. On the other hand, due to the lack of a good coupling structure similar to the linear diffusion case between the diffusion term and the haptotactic term, the method used for the linear diffusion case (see (1.2)) is also not suitable for the current nonlinear diffusion model. Therefore, we have to try new methods. Fortunately, we found a complex relationship between the haptotactic term and the nonlinear diffusion term. In addition to constructing classical functional such as

\[ F(u,\omega)=\int_\Omega \left(u\ln u+\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x, \]

we also construct the following new form of functionals

\begin{align*} F_1(u, \omega)& =\int_{\Omega} \left(\frac{1}{m-1}(u+\varepsilon)^m-\xi u\omega \right){\rm d}x,\\ G_k(u, \omega)& =\left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1}\\ & \quad \times \int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x, \end{align*}

which allows us to reconstruct the relationship between nonlinear diffusion term and haptotactic term. In particular, by taking the derivative of the above functional, we can use the diffusion term to neutralize the ‘bad’ effect brought by the haptotactic term. That is (see lemma 3.6)

\begin{align*} & \frac{{\rm d}}{{\rm d}t}G_k(u, \omega)+(2k+1)\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k} \\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x\cdots. \end{align*}

However, this also reduces the original good effect of the diffusion term. Fortunately, we have the logistic term, with the help of this term, we can finally get the $L^p$-norm ($\forall$ $p>1$)uniform estimation of $u$ and the $W^{1,\infty }$-norm estimation of $v$. However, it is hard to get the $L^\infty$-norm estimation of $u$ since the Moser's iteration technique is no longer applicable due to the diffusion term is not working.

Although the haptotactic term caused the main difficulties in the proof, we prove that the chemotaxis still plays a leading role in determining whether the solution can exist globally. More precisely, we prove that the weak solution will exists globally for any $m>\frac {2N}{N+2}$ if $\chi >0$, while for the haptotaxis-only model, that is the case $\chi =0$. We prove that for any $m>1$, the solution always exists globally. This work obviously improves the results in references[Reference Li and Lankeit12, Reference Tao and Winkler19, Reference Wang21, Reference Zheng24], in which, only the special case $\eta =0$ is studied.

In what follows, we give the assumptions of this paper.

\[ (H)\qquad\qquad \begin{cases} u_0, \omega_0 \in L^\infty(\Omega), \nabla\sqrt{\omega_{0}}\in L^2(\Omega),\\ \quad v_0 \in W^{1,\infty}(\Omega)\cap W^{2,p}(\Omega), & \text{for any }p>1, \\ u_0, v_0, \omega_0\ge 0 & . \end{cases} \]

We state the main results as follows.

Theorem 1.1 Assume that (H) holds, $\mu >0$ and $m>\frac {2N}{N+2}$. Then the problem (1.6) admits a global nonnegative weak solution $(u,\, v,\, \omega )\in \mathcal {X}_1\times \mathcal {X}_2\times \mathcal {X}_3$. In particular, $u,\, v,\, \omega$ are bounded uniformly in the following sense

(1.7)\begin{align} \sup_{0< t<\infty}\left(\|v({\cdot}, t)\|_{W^{1,\infty}}+\|\omega({\cdot}, t)\|_{L^{\infty}}+\|u({\cdot}, t)\|_{L^{r}}\right)\le C_r,\quad \text{for any $r>1$ }, \end{align}

where $C_r$ only depends on $r$, $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$. Here

\begin{align*} \mathcal{X}_1& =\Big\{u\in L^\infty(\mathbb{R}^+; L^r(\Omega)) \ \text{ for any}\ r>1; \nabla u^m\in L_{loc}^p(\mathbb{R}^+; L^p(\Omega)),\\ & \quad \qquad u_t\in L_{loc}^2(\mathbb{R}^+; W^{{-}1,p}(\Omega)), \\ & \quad \text{ for any}\ p\in (1,2), \nabla u^{m-\frac12}\in L_{loc}^2(\mathbb{R}^+; L^2(\Omega)) \Big\}; \\ \mathcal{X}_2& =\{v\in L^\infty(\mathbb{R}^+; W^{1,\infty}(\Omega)); v_t, D^2v\in L_{loc}^p(\mathbb{R}^+; L^p(\Omega)) \text{ for any}\ p>1 \}; \\ \mathcal{X}_3& =\{\omega\in L^\infty(Q);\nabla\sqrt\omega \in L_{loc}^\infty(\mathbb{R}^+; L^2(\Omega)), \sqrt u\nabla\sqrt\omega\in L_{loc}^2(\mathbb{R}^+; L^2(\Omega)), \\ & \quad\omega_t\in L_{loc}^r(\mathbb{R}^+; L^r(\Omega)) \ \text{ for any}\ r>1 \}. \end{align*}

Although the haptotactic term caused the main difficulties in the proof, the requirement of this index $m>\frac {N+2}{2N}$ was still caused by the chemotactic term. In fact, if $\chi =0$, we have the global existence result for any slow diffusion case $m>1$.

Theorem 1.2 Assume $\chi =0$, $\mu >0$, $m>1$ and (H) holds. Then the problem (1.6) admits a global nonnegative weak solution $(u,\, v,\, \omega )\in \mathcal {X}_1\times \mathcal {X}_2\times \mathcal {X}_3$. In particular, $u,\, v,\, \omega$ are bounded uniformly in the following sense

(1.8)\begin{equation} \sup_{0< t<\infty}\left(\|v({\cdot}, t)\|_{W^{1,\infty}}+\|\omega({\cdot}, t)\|_{L^{\infty}}+\|u({\cdot}, t)\|_{L^{r}}\right)\le C_r,\quad \text{for any }r>1, \end{equation}

where $C_r$ only depends on $r$, $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$.

Remark 1.3 As for the stability of the equilibrium points, the local stability of equilibrium point $(1,\,1,\,0)$ is easy to obtain by linear stability analysis, see for example [Reference Jin11]. However, the research on global asymptotic stability of the equilibrium point is a very challenging problem. In fact, this problem is still unsolved even for linear diffusion case $m=1$.

2. Preliminaries

We first give the definition of weak solutions.

Definition 2.1 $(u,\, v,\, \omega )$ is called a nonnegative weak solution of (1.6), if $(u,\, v,\, \omega )\in \mathcal {X}_1\times \mathcal {X}_2\times \mathcal {X}_2$, such that for any $T>0$,

\begin{align*} & -\iint_{Q_T}u\varphi_{1t} {\rm d}x{\rm d}t-\int_\Omega u(x,0)\varphi_1(x,0){\rm d}x+\iint_{Q_T}(\nabla u^m-\chi u\nabla v-\xi u\nabla \omega)\nabla\varphi_1 {\rm d}x{\rm d}t \\ & \quad=\mu\iint_{Q_T}u(1-u-\omega)\varphi_1 {\rm d}x{\rm d}t, \\ & -\iint_{Q_T}v\varphi_{2t}{\rm d}x{\rm d}t-\int_\Omega v(x,0)\varphi_2(x,0){\rm d}x+\iint_{Q_T}\nabla v\nabla\varphi_2 {\rm d}x{\rm d}t\\ & \quad =\iint_{Q_T}(u-v)\varphi_2 {\rm d}x{\rm d}t, \\ & -\iint_{Q_T}\omega\varphi_{3t} {\rm d}x{\rm d}t-\int_\Omega \omega(x,0)\varphi_3(x,0){\rm d}x+\iint_{Q_T}v\omega\varphi_3 {\rm d}x{\rm d}t\\ & \quad =\eta\iint_{Q_T}\omega(1-u-\omega)\varphi_3 {\rm d}x{\rm d}t, \end{align*}

for any $\varphi _1,\, \varphi _2,\, \varphi _3\in C^\infty (\overline Q_T)$ with $\varphi _i(x,\,T)=0$.

As a preparation, we introduce the following inequalities [Reference Jin8, Reference Jin9, Reference Mizoguchi and Souplet14], which is useful in the calculation of energy estimation.

Lemma 2.2 Assume that $\Omega$ is bounded with smooth boundary, and let $\omega \in C^2(\overline \Omega )$ satisfies $\frac {\partial \omega }{\partial \nu }\Big |_{\partial \Omega }=0$. Then we have

(2.1)\begin{equation} \frac{\partial|\nabla\omega|^2}{\partial\nu}\le 2\kappa|\nabla\omega|^2 \quad \text{on} \ \partial\Omega, \end{equation}

where $\kappa >0$ is an upper bound for the curvatures of $\Omega$. In particular, we have the following inequalities

(2.2)\begin{align} & \int_{\Omega}\frac{|\nabla \omega|^4}{\omega^3}{\rm d}x\le (2+\sqrt N)^2\int_{\Omega}\omega |D^2\ln \omega|^2 {\rm d}x, \end{align}

(2.3)\begin{align} & \int_\Omega|\nabla \omega|^{r+2} {\rm d}x\le (\sqrt N+r)^2\|\omega\|_{L^\infty}^2\int_\Omega |\nabla\omega|^{r-2}|D^2\omega|^2 {\rm d}x, \quad \text{for any }r\ge 2, \end{align}

and

(2.4)\begin{equation} \int_{\partial\Omega}\frac{1}{\omega} \frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S\le \delta\int_{\Omega} \omega |D^2\ln \omega|^2 {\rm d}x +C_\delta, \end{equation}

where $\delta >0$ is an arbitrary small constant, and $C_\delta$ is a constant depending on $\delta$.

Using Neumann heat semigroup theory, we prove the following lemma.

Lemma 2.3 Let $\Omega \subset \mathbb {R}^N$ be a bounded domain with smooth boundary. Assume that $u\in L^q(\Omega \times (0,\, T))$, $v_0\in W^{2,q}(\Omega )$ with $q\ge 2$. Then the following problem admits

\[ \begin{cases} v_t-\Delta v+v=u, & (x,t)\in \Omega\times(0, T) \\ \left.\displaystyle\frac{\partial v}{\partial{\bf n}}\right|_{\partial\Omega}=0 & \\ v(x,0)=v_0(x), & x\in\Omega \end{cases} \]

admits a unique strong solution $v\in W_q^{2,1}(\Omega \times (0,\, T))$, such that for any $r<\frac{Nq}{(N+2-q)_+}$,

(2.5)\begin{align} & \sup_{t\in (0, T)}\|\nabla v({\cdot}, t)\|_{L^r}\le C_1+C_2\left(\sup_{t\in(\tau, T)}\int_{t-\tau}^{t}\|u(s)\|_{L^q}^q{\rm d}s\right)^{\frac1q}, \end{align}

and

(2.6)\begin{align} & \sup_{t\in(\tau, T)}\int_{t-\tau}^{t}\|\nabla v\|_{L^{q+\frac{rq}N}}^{q+\frac{rq}N}{\rm d}s\le C_3\left(\sup_{t\in(\tau, T)}\int_{t-\tau}^{t}\|u(s)\|_{L^q}^q{\rm d}s\right)^{\frac{r}N+1}+C_4, \end{align}

where $\tau =\min \left \{1,\, \frac {T}2\right \}$, $C_i (i=1,\,2,\,3,\,4)$ are constants depending only on $v_0$, $\Omega$.

Proof. By Duhamel's principle, $v$ can be expressed as follows

\[ v=e^{{-}t}e^{t\Delta}v_0+\int_{0}^t e^{-(t-s)}e^{(t-s)\Delta} u(s){\rm d}s, \]

where $\{e^{t\Delta }\}_{t\ge 0}$ is the Neumann heat semigroup in $\Omega$. For any $r\in (1,\, +\infty )$, we have

\begin{align*} \|\nabla v({\cdot}, t)\|_{L^r}& \le e^{{-}t}\|\nabla v_0\|_{L^r}+\int_0^t e^{-(t-s)}\left(1+(t-s)^{-\frac N2(\frac 1q-\frac1r)-\frac12}\right)\|u(s)\|_{L^q}{\rm d}s \\ & \le e^{{-}t}\|\nabla v_0\|_{L^r}+\left(\int_0^t e^{-(t-s)}\left(1+(t-s)^{-\left(\frac{N}2(\frac 1q-\frac1r)+\frac12\right)\frac{q}{q-1}}\right){\rm d}s\right)^{\frac{q-1}q}\\ & \quad \times \left(\int_0^t e^{-(t-s)}\|u(s)\|_{L^q}^q{\rm d}s\right)^{\frac1q} \\ & \le e^{{-}t}\|\nabla v_0\|_{L^r}+\left(\int_0^t e^{{-}s}\left(1+s^{-\left(\frac{N}2(\frac 1q-\frac1r)+\frac12\right)\frac{q}{q-1}}\right){\rm d}s\right)^{\frac{q-1}q}\\ & \quad \times \left(\int_0^t e^{-(t-s)}\|u(s)\|_{L^q}^q{\rm d}s\right)^{\frac1q}. \end{align*}

A direct calculation gives

\[ \left(\frac{N}2(\frac 1q-\frac1r)+\frac12\right)\frac{q}{q-1}<1\Leftrightarrow \frac1 r>\frac{N+2-q}{Nq}. \]

Then when $r<\frac {Nq}{(N+2-q)_+}$,

\[ \|\nabla v({\cdot}, t)\|_{L^r}\le e^{{-}t}\|\nabla v_0\|_{L^r}+C\left(\int_0^t e^{-(t-s)}\|u(s)\|_{L^q}^q{\rm d}s\right)^{\frac1q}. \]

We also note that for any $t>0$, there exists an integer $N$ such that $t=N\tau +\sigma$ with $0<\sigma <\tau$. Therefore, we have

\begin{align*} \int_0^t e^{-(t-s)}\|u(s)\|_{L^q}^q{\rm d}s & =e^{{-}t}\left(\sum_{k=1}^N\int_{(k-1)\tau}^{k\tau}e^{s}\|u(s)\|_{L^q}^q{\rm d}s+\int_{N\tau}^{N\tau+\sigma}e^{s}\|u(s)\|_{L^q}^q{\rm d}s\right) \\ & \le e^{{-}t}\left(\sum_{k=1}^N e^{k\tau}\int_{(k-1)\tau}^{k\tau}\|u(s)\|_{L^q}^q{\rm d}s+e^{t}\int_{N\tau}^{t}\|u(s)\|_{L^q}^q{\rm d}s\right) \\ & \le \frac{e^{2\tau}}{e^{\tau}-1} \sup_{\tau< a< t}\int_{a-\tau}^{a}\|u(s)\|_{L^q}^q{\rm d}s, \end{align*}

which implies that (2.5). By $L^p$ theory of linear parabolic equations, we obtain

(2.7)\begin{equation} \sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\|D^2 v({\cdot}, t)\|_{L^q}^q \le C_3+C_4\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\|u(s)\|_{L^q}^q{\rm d}s, \ \text{for any} \ q>1. \end{equation}

From Gagliardo-Nirenberg interpolation inequality, we infer that

(2.8)\begin{equation} \|\nabla v\|_{L^{q+\frac{rq}N}}^{q+\frac{rq}N}\le C_5\|\nabla v\|_{L^r}^{\frac{rq}N}\|D^2 v\|_{L^q}^q+C_6\|\nabla v\|_{L^r}^{q+\frac{rq}N}. \end{equation}

Using (2.5), (2.7) and (2.8), we derive (2.6).

3. Uniform energy estimations and global solvability of chemotaxis-haptotaxis system

In order to prove the existence of the weak solution of problem (1.6), some prior estimates are necessary. Since problem (1.6) is degenerate at $u = 0$, in order to facilitate us to obtain various prior estimates later, we first use the standard method to build a framework, that is, consider its regularization problem

(3.1)\begin{align} \left\{ \begin{aligned} & u_t=m\nabla\cdot ((u+\varepsilon)^{m-1}\nabla u)-\chi\nabla\cdot(u\nabla v)-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega),\quad \text{in}\ Q, \\ & v_t=\Delta v+u-v, \quad \text{in}\ Q, \\ & \omega_t=\varepsilon\Delta\omega -v\omega+\eta \omega(1-u-\omega), \ \text{in}\ Q, \\ & \left.\frac{\partial u}{\partial{\bf n}}\right|_{\partial\Omega}=0, \left.\frac{\partial v}{\partial{\bf n}}\right|_{\partial\Omega}=0,\left.\frac{\partial \omega}{\partial{\bf n}}\right|_{\partial\Omega}=0, \\ & u(x,0)=u_{\varepsilon 0}(x), \quad v(x,0)=v_{\varepsilon 0}(x), \quad \omega(x,0)=\omega_{\varepsilon 0}(x), \quad x\in\Omega, \end{aligned}\right. \end{align}

for any $\varepsilon \in (0,\,1)$, and $(u_{\varepsilon 0},\, v_{\varepsilon 0},\, \omega _{\varepsilon 0})$ is smooth approximation of $(u_0,\, v_0,\, \omega _0)$ with $(u_{\varepsilon 0},\, v_{\varepsilon 0},\, \omega _{\varepsilon 0})$ sufficiently smooth. By standard fixed point method, it is easy to obtain the local existence of classical solutions to the above problem, see for example [Reference Tao and Winkler19]. That is

Lemma 3.1 Assume $B$, $m>1$, $\varepsilon \in (0,\,1)$. Then there exists $T_{\max }\in (0,\, +\infty ]$ such that the problem (3.1) admits a classical solution $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )\in C^{2, 1}(Q_{T_{\max }})\cap C^0(\overline \Omega \times [0,\,T_{\max }))$ with

\[ u_\varepsilon> 0, \quad v_\varepsilon> 0, \quad \omega_\varepsilon> 0 \quad \text{for all}\ (x, t)\in \Omega\times (0, T_{\max}). \]

Moreover, we have the following dichotomy: Either $T_{\max } = \infty$, or

\[ \limsup_{t\nearrow T_{\max}}\left(\|u_\varepsilon({\cdot}, t)\|_{L^\infty}+\|v_\varepsilon\|_{W^{1,\infty}}+\|\omega_\varepsilon\|_{W^{1,\infty}}\right) = \infty. \]

Throughout this paper, we denote

\[ \tau:=\min\left\{1, \frac{T_{\max}}2\right\}\le 1. \]

In what follows, we focus on the uniform energy estimates of $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )$. For simplicity, we omit the subscript $\varepsilon$ of the approximate solutions $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )$ in the subsequent energy estimate calculations.

Firstly, it is easy to obtain the following estimates.

Lemma 3.2 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Then

(3.2)\begin{align} & \sup_{t\in(0, T_{\max})}\|u({\cdot}, t)\|_{L^1}+\mu\sup_{t\in(\tau, T_{\max})}\int_{t-\tau}^t\int_\Omega (u^2+u\omega) {\rm d}x{\rm d}s\le C_1, \end{align}

(3.3)\begin{align} & \sup_{t\in(0, T_{\max})}\left(\|\omega\|_{L^\infty}+\|v({\cdot}, t)\|_{H^1}^2\right)+\sup_{t\in(\tau, T_{\max})}\int_{t-\tau}^t\|v({\cdot}, s)\|_{H^2}^2{\rm d}s\le C_2, \end{align}

where the constants $C_1$, $C_2$ only depend on $\mu$, $u_0$, $v_0$, $\omega _0$ and $\Omega$.

Proof. By comparison lemma, it is easy to see that

\[ \omega\le \max\{1, \|\omega_{\varepsilon 0}\|_{L^\infty}\}\le \max\{1, 2\|\omega_0\|_{L^\infty}\}. \]

By direct integration over $\Omega$ for the first equation of (3.1), (3.2) is easy to be obtained. By $L^2$ theory of linear parabolic equations, (3.3) is arrived.

Lemma 3.3 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Then for any $T< T_{\max }$,

(3.4)\begin{align} & \sup_{0< t< T}\int_\Omega \left(u\ln u+\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x\nonumber\\ & \quad +\int_0^T\int_\Omega\left( \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 +\frac{ |\nabla \omega|^2}{\omega} (u+v)+u^2|\ln u|\right){\rm d}x{\rm d}t \nonumber\\ & \quad+\varepsilon\int_0^T\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x{\rm d}t\le C_T, \end{align}

where $C_T$ is independent of $\varepsilon$, it depends only on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$, and $T$.

In particular, we also have

(3.5)\begin{align} & \sup_{0< t< T_{\max}}\int_\Omega \left(u\ln u+\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x\nonumber\\ & \quad + \sup_{t\in(\tau, T_{\max})}\int_{t-\tau}^t\int_\Omega \left( \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 +\frac{ |\nabla \omega|^2}{\omega} (u+v)+u^2|\ln u|\right){\rm d}x{\rm d}t \nonumber\\ & \quad+\varepsilon\sup_{t\in(\tau, T_{\max})}\int_{t-\tau}^t\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x{\rm d}t\le \frac{C}\varepsilon, \end{align}

where $C$ is independent of $\varepsilon$ and $T_{\max }$, it depends only on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$.

It should be noted that the estimation (3.4) is time-dependent, but (3.5) is not time-dependent. Although (3.5) depends on $\varepsilon$, it is useful for us to obtain uniform estimations independent of time and $\varepsilon$.

Proof. By direct calculation,

\begin{align*} & \frac12\frac{{\rm d}}{{\rm d}t}\int_\Omega\frac{ |\nabla \omega|^2}{\omega}{\rm d}x=\int_\Omega \frac{ \nabla \omega}{\omega} \nabla \omega_t {\rm d}x-\frac12\int_\Omega\frac{ |\nabla \omega|^2}{\omega^2} \omega_t{\rm d}x \\ & \quad={-}\int_\Omega \omega_t \left(\frac{ \Delta \omega}{\omega} -\frac{ |\nabla \omega|^2}{\omega^2} \right){\rm d}x-\frac12\int_\Omega\frac{ |\nabla \omega|^2}{\omega^2} \omega_t{\rm d}x \\ & \quad= \frac12\int_\Omega\frac{ |\nabla \omega|^2}{\omega^2} \omega_t{\rm d}x-\int_\Omega \frac{ \Delta \omega}{\omega} \omega_t {\rm d}x \\ & \quad= \frac12\int_\Omega\frac{ |\nabla \omega|^2}{\omega^2} (\varepsilon\Delta\omega -v\omega+\eta \omega(1-u-\omega)){\rm d}x\\ & \qquad -\int_\Omega \frac{ \Delta \omega}{\omega}(\varepsilon\Delta\omega -v\omega+\eta \omega(1-u-\omega)) {\rm d}x \\ & \quad= \frac\varepsilon 2\int_\Omega\left(\frac{ |\nabla \omega|^2}{\omega^2}\Delta\omega-2\frac{ |\Delta \omega|^2}{\omega}\right){\rm d}x-\int_\Omega\nabla\omega\nabla v {\rm d}x -\eta\int_\Omega\nabla\omega\nabla u {\rm d}x\\ & \quad +\frac12\int_\Omega\frac{ |\nabla \omega|^2}{\omega} (\eta-3\eta\omega-\eta u-v){\rm d}x. \end{align*}

Similar to the proof of lemma 2.6 in [Reference Jin8], we have

\[ \frac\varepsilon 2\int_\Omega\left(\frac{ |\nabla \omega|^2}{\omega^2}\Delta\omega-2\frac{ |\Delta \omega|^2}{\omega}\right){\rm d}x ={-}\varepsilon\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x +\frac{\varepsilon}2\int_{\partial\Omega}\frac{1}{\omega} \frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S. \]

Putting the above two equalities together, we get the following conclusion

(3.6)\begin{align} & \frac12\frac{{\rm d}}{{\rm d}t}\int_\Omega\frac{ |\nabla \omega|^2}{\omega}{\rm d}x+\varepsilon\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x +\frac12\int_\Omega\frac{ |\nabla \omega|^2}{\omega} (3\eta\omega+\eta u+v){\rm d}x \nonumber\\ & \quad={-}\int_\Omega\nabla\omega\nabla v {\rm d}x -\eta\int_\Omega\nabla\omega\nabla u {\rm d}x+\frac{\eta}2\int_\Omega\frac{ |\nabla \omega|^2}{\omega} {\rm d}x +\frac{\varepsilon}2\int_{\partial\Omega}\frac{1}{\omega}\frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S. \end{align}

Multiplying both sides of the first equation of (3.1) by $1+\ln u$, and integrating the resultant equation over $\Omega$ yields

(3.7)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_\Omega u\ln u {\rm d}x+m\int_\Omega \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 {\rm d}x \nonumber\\ & \quad=\chi\int_\Omega \nabla v\nabla u {\rm d}x+\xi\int_\Omega \nabla \omega\nabla u {\rm d}x+\mu\int_\Omega u(1+\ln u)(1-u-\omega){\rm d}x. \end{align}

The combination of (3.6) and (3.7) leads to

(3.8)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_\Omega \left(\eta u\ln u+\frac{\xi}2\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x+\eta m\int_\Omega \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 {\rm d}x+\xi\varepsilon\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x \nonumber\\ & +\frac{\xi}2\int_\Omega\frac{ |\nabla \omega|^2}{\omega} (3\eta\omega+\eta u+v){\rm d}x+\mu\eta\int_\Omega u(1+\ln u)(u+\omega){\rm d}x \nonumber\\ & \quad=\frac{\varepsilon\xi}2\int_{\partial\Omega}\frac{1}{\omega}\frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S+\frac{\xi\eta}2\int_\Omega\frac{|\nabla \omega|^2}{\omega}{\rm d}x +\chi\eta\int_\Omega \nabla v\nabla u {\rm d}x\nonumber\\ & \qquad-\xi\int_\Omega \nabla \omega\nabla v {\rm d}x+\mu\eta\int_\Omega u(1+\ln u){\rm d}x \nonumber\\ & \quad=\frac{\varepsilon\xi}2\int_{\partial\Omega}\frac{1}{\omega}\frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S+\frac{\xi\eta}2\int_\Omega\frac{|\nabla \omega|^2}{\omega}{\rm d}x -\chi\eta\int_\Omega u \Delta v {\rm d}\nonumber\\ & \qquad +\xi\int_\Omega \omega\Delta v {\rm d}x+\mu\eta\int_\Omega u(1+\ln u) \nonumber\\ & \quad\le \frac{\varepsilon\xi}2\int_{\partial\Omega}\frac{1}{\omega}\frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S+\frac{\xi\eta}2\int_\Omega\frac{|\nabla \omega|^2}{\omega}{\rm d}x +\int_\Omega (\chi^2\eta^2 u^2+\xi^2\omega^2\nonumber\\ & \qquad +|\Delta v|^2){\rm d}x+\frac{\mu\eta}4\int_\Omega u^2\ln u {\rm d}x+C. \end{align}

Recalling (2.4), we have

(3.9)\begin{equation} \frac{\varepsilon\xi}2\int_{\partial\Omega}\frac{1}{\omega}\frac{\partial}{\partial{\bf n}}|\nabla \omega|^2 {\rm d}S\le \frac{\varepsilon\xi}{2} \int_{\Omega} \omega |D^2\ln \omega|^2 {\rm d}x +C\varepsilon.\end{equation}

We substitute (3.9) into (3.8) to obtain

(3.10)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_\Omega \left(\eta u\ln u+\frac{\xi}2\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x+\eta m\int_\Omega \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 {\rm d}x+\frac{\xi\varepsilon}2\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x \nonumber\\ & +\frac{\xi}2\int_\Omega\frac{ |\nabla \omega|^2}{\omega} (3\eta\omega+\eta u+v){\rm d}x+\frac{\mu\eta}2\int_\Omega u^2(1+|\ln u|){\rm d}x \nonumber\\ & \quad\le\frac{\xi\eta}2\int_\Omega\frac{|\nabla \omega|^2}{\omega}{\rm d}x +\int_\Omega (\chi^2\eta^2 u^2+\xi^2\omega^2+|\Delta v|^2) {\rm d}x+\hat C. \end{align}

Using (3.2) and (3.3), by direct calculation, we complete the proof of (3.4).

On the other hand, we note that

(3.11)\begin{align} \frac{\xi(\eta+1)}2\int_\Omega\frac{|\nabla \omega|^2}{\omega}{\rm d}x={-}\frac{\xi(\eta+1)}2\int_\Omega\omega\Delta\ln\omega {\rm d}x \le \frac{\varepsilon\xi}4\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x+\frac{C_1}{\varepsilon},\end{align}

which together with (3.10) yields

(3.12)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_\Omega \left(\eta u\ln u+\frac{\xi}2\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x+\eta m\int_\Omega \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 {\rm d}x+\frac{\xi\varepsilon}2\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x \nonumber\\ & \qquad+\frac{\xi}2\int_\Omega\frac{ |\nabla \omega|^2}{\omega} (3\eta\omega+\eta u+v+1){\rm d}x+\frac{\mu\eta}2\int_\Omega u^2(1+|\ln u|){\rm d}x \nonumber\\ & \quad\le\frac{\xi(\eta+1)}2\int_\Omega\frac{|\nabla \omega|^2}{\omega}{\rm d}x +\int_\Omega (\chi^2\eta^2 u^2+\xi^2\omega^2+|\Delta v|^2) {\rm d}x+\hat C \nonumber\\ & \quad\le\frac{\xi\varepsilon}4\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x+\int_\Omega (\chi^2\eta^2 u^2+\xi^2\omega^2+|\Delta v|^2) {\rm d}x+\frac{\tilde C}{\varepsilon}, \end{align}

where $\tilde C$ is independent of $\varepsilon$. Denote

\[ f(t)=\int_\Omega \left(\eta u\ln u+\frac{\xi}2\frac{ |\nabla \omega|^2}{\omega} \right){\rm d}x, \]

clearly,

\begin{align*} & \eta m\int_\Omega \frac{(u+\varepsilon)^{m-1}}{u}|\nabla u|^2 {\rm d}x+\frac{\xi\varepsilon}4\int_\Omega \omega|D^2\ln \omega|^2 {\rm d}x\\& \quad+\frac{\xi}2\int_\Omega\frac{ |\nabla \omega|^2}{\omega} (3\eta\omega+\eta u+v+1){\rm d}x \\ & \quad +\frac{\mu\eta}2\int_\Omega u^2(1+|\ln u|){\rm d}x\ge\frac{\min\{\mu, \xi\}}4f(t)-C. \end{align*}

Then (3.12) implies that

\[ f'(t)+\frac{\min\{\mu, \xi\}}4f(t)\le \int_\Omega (\chi^2\eta^2 u^2+\xi^2\omega^2+|\Delta v|^2) {\rm d}x+\frac{\tilde C}{\varepsilon}+C. \]

Then the uniform boundedness of $f$ is derived from (3.2) and (3.3). Furthermore, (3.5) is obtained by integrating (3.12) from $t-\tau$ to $t$.

Next, our purpose is to continuously improve the regularity of the solution. In this process, the key is to deal with the haptotaxis term. We try to offset the influence of the haptotaxis term with the diffusion term. First, we can get the following lemma.

Lemma 3.4 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$, and assume $m>\frac {2N}{N+2}$. Then

(3.13)\begin{align} & \sup_{0< t< T_{\max}}\int_{\Omega} u^m {\rm d}x +\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}\left((u+\varepsilon)\vphantom{\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2}\right.\nonumber\\ & \quad \left.\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2 +u^{m+1}\right){\rm d}x{\rm d}s\le C, \end{align}

where $C$ is independent of $\varepsilon$ and $T_{\max }$, and it depends only on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$.

In addition, we also have

(3.14)\begin{equation} \int_0^T\int_{\Omega} (u+\varepsilon)^{2m-3} |\nabla u|^2{\rm d}x{\rm d}t\le C_T, \end{equation}

and for any $1< p<2$,

(3.15)\begin{equation} \int_0^T\int_{\Omega}|\nabla(u+\varepsilon)^{m}|^p{\rm d}x{\rm d}t\le C_{Tp}, \end{equation}

where $C_T$, $C_{Tp}$ are independent of $\varepsilon$, and $C_T$ depends on $T$, $C_{Tp}$ depends on $T$ and $p$.

Proof. Noticing that

\begin{align*} \omega|D^2\ln\omega|^2& =\frac{|D^2\omega|^2}{\omega}+\frac{|\nabla\omega|^4}{\omega^3}-2\frac{\nabla\omega D^2\omega\nabla\omega}{\omega^2} \\ & \ge \frac{|D^2\omega|^2}{\omega}+\frac{|\nabla\omega|^4}{\omega^3}-\frac12\frac{|D^2\omega|^2}{\omega}-2\frac{|\nabla\omega|^4}{\omega^3}, \end{align*}

that is

\[ \frac{|D^2\omega|^2}{\omega}\le 2\omega|D^2\ln\omega|^2+2\frac{|\nabla\omega|^4}{\omega^3}. \]

Recalling (2.2), we conclude

(3.16)\begin{align} \int_\Omega \frac{|D^2\omega|^2}{\omega}{\rm d}x& \le 2\int_\Omega\omega|D^2\ln\omega|^2{\rm d}x+2\int_\Omega\frac{|\nabla\omega|^4}{\omega^3}{\rm d}x \nonumber\\ & \le \left(2+2(2+\sqrt N)^2\right)\int_\Omega\omega|D^2\ln\omega|^2{\rm d}x.\end{align}

Noticing that

(3.17)\begin{align} & m\nabla\cdot ((u+\varepsilon)^{m-1}\nabla u)-\xi\nabla\cdot(u\nabla\omega)\nonumber\\ & \quad=\nabla\cdot\left(\frac{m}{m-1}(u+\varepsilon)\nabla (u+\varepsilon)^{m-1} -\xi(u+\varepsilon)\nabla\omega\right)+\xi\varepsilon\Delta\omega \nonumber\\ & \quad=\nabla\cdot\left((u+\varepsilon)\nabla \left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right)+\xi\varepsilon\Delta\omega, \end{align}

multiplying both sides of the first equation of (3.1) by $\frac {m}{m-1}(u+\varepsilon )^{m-1}-\xi \omega$, and using (3.17), (3.5) and lemma 3.2, we arrive at

\begin{align*} & \frac{{\rm d}}{{\rm d}t}\int_{\Omega}\left(\frac1{m-1}(u+\varepsilon)^m -\xi u\omega\right) {\rm d}x +\int_{\Omega}(u+\varepsilon)\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \quad={-}\varepsilon\xi\int_{\Omega}\nabla\omega\nabla \left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right){\rm d}x\\ & \qquad +\chi\int_{\Omega}u\nabla v\nabla\left(\frac{m}{m-1}u^{m-1}-\xi\omega\right){\rm d}x \\ & \qquad+\int_{\Omega}\mu u(1-u-\omega)\left(\frac{m}{m-1}u^{m-1}-\xi\omega\right) {\rm d}x \\ & \qquad+\xi\int_{\Omega}u(-\varepsilon\Delta\omega+v\omega-\eta\omega(1-u-\omega)){\rm d}x \\ & \quad\le\frac{\varepsilon\xi^2\|\omega\|_{L^\infty}}2\int_{\Omega}\frac{|\nabla\omega|^2}{\omega} {\rm d}x +\frac{\varepsilon}2\int_{\Omega}\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x\\ & \qquad +\frac{\chi^2}2\int_{\Omega}u|\nabla v|^2 {\rm d}x +\frac12\int_{\Omega}u\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x\\ & \qquad -\frac{3m\mu}{4(m-1)}\int_{\Omega}u^{m+1}{\rm d}x+\varepsilon^2\int_{\Omega}|\Delta\omega|^2{\rm d}x +C \\ & \quad\le \frac12\int_{\Omega}(u+\varepsilon)\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \qquad -\frac{m\mu}{2(m-1)}\int_{\Omega}u^{m+1}{\rm d}x+C\int_{\Omega}|\nabla v|^{\frac{2(m+1)}m} {\rm d}x +N\varepsilon^2\int_{\Omega}|D^2\omega|^2{\rm d}x+\tilde C, \end{align*}

namely,

(3.18)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_{\Omega}\left(\frac1{m-1}(u+\varepsilon)^m -\xi u\omega\right) {\rm d}x \nonumber\\ & \quad+\frac12\int_{\Omega}(u+\varepsilon)\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \nonumber\\ & \quad+\frac{m\mu}{2(m-1)}\int_{\Omega}u^{m+1}{\rm d}x\le C\int_{\Omega}|\nabla v|^{\frac{2(m+1)}m} {\rm d}x+N\varepsilon^2\int_{\Omega}|D^2\omega|^2{\rm d}x+\tilde C. \end{align}

Taking advantage of Gagliardo-Nirenberg interpolation inequality and (3.3), we have

(3.19)\begin{align} \|\nabla v\|_{L^{\frac{(N+2)(m+1)}{N}}}^{\frac{(N+2)(m+1)}{N}}& \le C_1\|\nabla v\|_{L^2}^{\frac{2(m+1)}{N}}\|D^2 v\|_{L^{m+1}}^{m+1}+ C_2\|\nabla v\|_{L^2}^{\frac{(N+2)(m+1)}{N}}\nonumber\\ & \le C_3\|D^2 v\|_{L^{m+1}}^{m+1}+C_4. \end{align}

Noticing that $\frac {(N+2)(m+1)}{N}>\frac {2(m+1)}m$ when $m>\frac {2N}{N+2}$, we infer from (3.18), (3.19),

(3.20)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_{\Omega}\left(\frac1{m-1}(u+\varepsilon)^m -\xi u\omega\right) {\rm d}x +\frac12\int_{\Omega}(u+\varepsilon)\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2\!{\rm d}x \nonumber\\ & \quad+\frac{m\mu}{2(m-1)}\int_{\Omega}u^{m+1}{\rm d}x\le \sigma\int_{\Omega}|D^2 v|^{m+1} {\rm d}x+N\varepsilon^2\int_{\Omega}|D^2\omega|^2{\rm d}x+C_\sigma \end{align}

for any small constant $\sigma >0$, with $C_\sigma$ depending on $\sigma$. Noticing that

\[ \int_{\Omega}u^{m+1}{\rm d}x\ge \int_{\Omega}\left(\frac1{m-1}(u+\varepsilon)^m -\xi u\omega\right) {\rm d}x-C, \]

and using $L^p$ theory of linear parabolic equations, we conclude that

\begin{align*} & \sup_{0< t< T_{\max}}\int_{\Omega} (u+\varepsilon)^m {\rm d}x +\frac12\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}(u+\varepsilon)\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2 {\rm d}x{\rm d}s \\ & \quad+\frac{m\mu}{2(m-1)}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}u^{m+1}{\rm d}x{\rm d}s \\ & \qquad\le C_5\sigma\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}|D^2 v|^{m+1} {\rm d}x{\rm d}s\\ & \qquad \quad +C_6\varepsilon^2\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}|D^2\omega|^2{\rm d}x{\rm d}s+\tilde C_\sigma \\ & \qquad\le C_7\sigma\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}|u|^{m+1} {\rm d}x{\rm d}s\\ & \qquad \quad+ C_6\varepsilon^2\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}|D^2\omega|^2{\rm d}x{\rm d}s +\hat C_\sigma. \end{align*}

Here, all these constants $C_i$, $\hat C$, $\tilde C$ are independent of $\varepsilon$. By the arbitrariness of $\sigma$, using (3.5), (3.16), and (3.13) is proved by taking $\sigma$ appropriately small in the above inequality.

In addition, using (3.3), (3.4) and (3.13), and noticing that

\begin{align*} (u+\varepsilon)\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}\right)\right|^2& \le 2(u+\varepsilon)\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2 \\ & \quad +2\xi^2(u+\varepsilon)\left|\nabla\omega\right|^2, \end{align*}

we arrive at (3.14).

For any $1< p<2$, we also notice that

\begin{align*} \|\nabla(u+\varepsilon)^{m}\|_{L^p}^p & =\|m(u+\varepsilon)^{\frac12}(u+\varepsilon)^{m-\frac32}\nabla u\|_{L^p}^p \\ & \le \|m(u+\varepsilon)^{m-\frac32}\nabla u\|_{L^2}^p\|u+\varepsilon\|_{L^{\frac{p}{2-p}}}^{\frac p2} \\ & \le \|m(u+\varepsilon)^{m-\frac32}\nabla u\|_{L^2}^2+C_p. \end{align*}

Then (3.15) is derived from (3.14) and the above inequality.

To improve the regularity of $u$, we need the following estimate on $\omega$. Although this estimation does not depend on $T_{\max }$, it depends on $\varepsilon$. In order to obtain some subsequent uniform estimates that does not depend on $\varepsilon$, we need to the exact order of its dependence on $\varepsilon$.

Lemma 3.5 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Then for any $r>2$,

(3.21)\begin{align} & \sup_{0< t< T_{\max}}\int_\Omega |\nabla \omega|^{r} {\rm d}x+ \sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}\int_\Omega \varepsilon|\nabla \omega|^{r+2} {\rm d}x{\rm d}s\nonumber\\ & \quad \le \frac{\tilde C}{\varepsilon^{\frac r2}}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u^r+v^r+1){\rm d}x{\rm d}s,\end{align}

where $\tilde C$ is independent of $\varepsilon$ and $T_{\max }$, it depends only on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $r$, $\Omega$.

Proof. Applying $\nabla$ to the third equation of (3.1), multiplying both sides of the resultant equation by $|\nabla \omega |^{r-2}\nabla \omega$ for any $r>2$, and using lemma 2.2, we obtain

(3.22)\begin{align} & \frac1r\frac{{\rm d}}{{\rm d}t}\int_\Omega |\nabla \omega|^{r} {\rm d}x +\varepsilon\int_{\Omega}|\nabla \omega|^{r-2}|\nabla^2 \omega|^2 {\rm d}x\nonumber\\ & \quad +(r-2)\varepsilon\int_{\Omega}|\nabla\omega|^{r-2}(\nabla|\nabla \omega|)^2 {\rm d}x +\int_{\Omega} |\nabla\omega|^{r}{\rm d}x \nonumber\\ & \qquad=\int_{\Omega}(v\omega-\eta\omega(1-u-\omega)) \nabla\cdot( |\nabla \omega|^{r-2}\nabla \omega) {\rm d}x\nonumber\\ & \quad \qquad +\frac{\varepsilon}2\int_{\partial\Omega}\frac{\partial(|\nabla \omega|^2)}{\partial{\bf n}}|\nabla \omega|^{r-2} {\rm d}S +\int_{\Omega} |\nabla \omega|^{r}{\rm d}x \nonumber\\ & \qquad\le \frac{\varepsilon}4\int_{\Omega}|\nabla \omega|^{r-2}|\nabla^2 \omega|^2 {\rm d}x+\frac{(r-2)\varepsilon}4\int_{\Omega}|\nabla\omega|^{r-2}(\nabla|\nabla \omega|)^2 {\rm d}x\nonumber\\ & \quad \qquad +\frac{C_1}{\varepsilon}\int_{\Omega}|\nabla \omega|^{r-2}(u^2+v^2){\rm d}x +\int_{\Omega} |\nabla \omega|^{r}{\rm d}x+\kappa\varepsilon\int_{\partial\Omega}|\nabla \omega|^r {\rm d}S. \end{align}

Recalling (2.3), and noticing that $\omega$ is bounded, then there exists a constant $\rho$ such that

(3.23)\begin{equation} \frac{\varepsilon}4\int_{\Omega}|\nabla \omega|^{r-2}|\nabla^2 \omega|^2 {\rm d}x\ge \rho\varepsilon\int_{\Omega}|\nabla\omega|^{r+2}{\rm d}x.\end{equation}

By the boundary trace embedding inequalities, we conclude that for any small $\delta >0$,

(3.24)\begin{equation} \kappa\varepsilon\int_{\partial\Omega}|\nabla \omega|^r {\rm d}S\le \frac{(r-2)\varepsilon}4 \int_{\Omega}|\nabla\omega|^{r-2}(\nabla|\nabla \omega|)^2 {\rm d}x+C_2\int_{\Omega} |\nabla\omega|^{r}{\rm d}x.\end{equation}

The combination of (3.22)–(3.24) leads to

(3.25)\begin{align} & \frac1r\frac{{\rm d}}{{\rm d}t}\int_\Omega |\nabla \omega|^{r} {\rm d}x +\frac{\varepsilon}2\int_{\Omega}|\nabla \omega|^{r-2}|\nabla^2 \omega|^2 {\rm d}x+\frac{(r-2)\varepsilon}2\int_{\Omega}|\nabla\omega|^{r-2}(\nabla|\nabla \omega|)^2 {\rm d}x\nonumber\\ & \quad +\int_{\Omega} |\nabla\omega|^{r}{\rm d}x +\rho\varepsilon\int_{\Omega}|\nabla\omega|^{r+2}{\rm d}x \nonumber\\ & \qquad\le \frac{C_1}{\varepsilon}\int_{\Omega}|\nabla \omega|^{r-2}(u^2+v^2){\rm d}x+(C_2+1)\int_{\Omega} |\nabla \omega|^{r}{\rm d}x \nonumber\\ & \qquad \le C_3\int_{\Omega} |\nabla \omega|^{r}{\rm d}x+\frac{C_4}{\varepsilon^{\frac r2}}\int_{\Omega}(u^r+v^r){\rm d}x \nonumber\\ & \qquad\le \frac{\rho}2\varepsilon\int_{\Omega}|\nabla\omega|^{r+2}{\rm d}x+\frac{C_4}{\varepsilon^{\frac r2}}\int_{\Omega}(u^r+v^r){\rm d}x+\frac{C_5}{\varepsilon^{\frac r2}}. \end{align}

It implies that

(3.26)\begin{equation} \sup_{0< t< T_{\max}}\int_\Omega |\nabla \omega|^{r} {\rm d}x\le \frac{C_6}{\varepsilon^{\frac r2}}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u^r+v^r+1){\rm d}x{\rm d}s.\end{equation}

Using (3.26), and integrating (3.25) directly, we arrive at

(3.27)\begin{equation} \sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} \varepsilon|\nabla \omega|^{r+2} {\rm d}x{\rm d}s\le \frac{C_7}{\varepsilon^{\frac r2}}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u^r+v^r+1){\rm d}x{\rm d}s.\end{equation}

Here these constants $C_i$ are independent of $\varepsilon$. Then (3.21) is derived from (3.26) and (3.27).

Based on the above lemmas, we can improve the regularity of the solution, and obtain the following result.

Lemma 3.6 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Assume that $m>\frac {2N}{N+2}$. Then for any positive integer $k$, we have

(3.28)\begin{align} & \sup_{0< t< T_{\max}}\int_{\Omega}(u+\varepsilon)^{(m-1)(2k+1)+1} {\rm d}x+\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s \nonumber\\ & \quad+\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\nonumber\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x{\rm d}s\le C, \end{align}

where $C$ is independent of $\varepsilon$ and $T_{\max }$, it depends only on $k$, $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$.

Proof. For any given positive integer $k$, multiplying both sides of the first equation of (3.1) by $(\frac {m}{m-1}(u+\varepsilon )^{m-1}-\xi \omega )^{2k+1}$, recalling (3.17), and noticing that $\omega$ is bounded, then we arrive at

\begin{align*} & \int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k+1}u_t {\rm d}x \\ & \quad+(2k+1)\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \quad={-}\varepsilon\xi(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \qquad \times \nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\nabla \omega {\rm d}x \\ & \qquad+\chi(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \qquad \times u\nabla v\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right){\rm d}x \\ & \qquad+\int_{\Omega}\mu u(1-u-\omega)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k+1} {\rm d}x \\ & \quad\le \frac{2k+1}2\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \qquad \times\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \qquad+\varepsilon\xi^2(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}|\nabla \omega|^2{\rm d}x \\ & \qquad+\chi^2(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}u|\nabla v|^2{\rm d}x \\ & \qquad -\frac{3\mu}4\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x+C_1. \end{align*}

That is

(3.29)\begin{align} & \int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k+1}u_t {\rm d}x+\frac{3\mu}4\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x \nonumber\\ & \quad+\frac{2k+1}2\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\nonumber\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \nonumber\\ & \qquad\le \varepsilon\xi^2(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}|\nabla \omega|^2{\rm d}x \nonumber\\ & \quad \qquad+\chi^2(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}u|\nabla v|^2{\rm d}x +C_1 \nonumber\\ & \qquad=I+II+C_1. \end{align}

We first calculate the first term on the left of the above inequality.

(3.30)\begin{align} & \int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k+1}u_t {\rm d}x\nonumber\\ & \quad =\int_{\Omega}\sum_{i=0}^{2k+1}(-\xi)^i C_{2k+1}^i\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}\right)^{2k+1-i}\omega^i u_t{\rm d}x \nonumber\\ & \quad=\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1}\left((u+\varepsilon)^{(m-1)(2k+1-i)+1}\right)_t\omega^i {\rm d}x \nonumber\\ & \quad=\left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1}\left(\frac{{\rm d}}{{\rm d}t}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x\right. \nonumber\\ & \qquad\left.-i\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^{i-1}\omega_t {\rm d}x\right) \nonumber\\ & \quad=\left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1} \frac{{\rm d}}{{\rm d}t}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x \nonumber\\ & \qquad-\left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=1}^{2k+1}\frac{i(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^{i-1}\omega_t {\rm d}x \nonumber\\ & \quad \ge\left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1} \frac{{\rm d}}{{\rm d}t}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x \nonumber\\ & \quad-C\int_{\Omega}\left(u^{2k(m-1)+2}+|\omega_t|^{2k(m-1)+2}+1\right) {\rm d}x \nonumber\\ \ & \quad \ge\left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1} \frac{{\rm d}}{{\rm d}t}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x \nonumber\\ & \qquad-\frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x -C_2\int_{\Omega}\left(|\omega_t|^{2k(m-1)+2}+1\right) {\rm d}x. \end{align}

In what follows, we calculate $I$, $II$ respectively. Using Young's inequality and recalling (3.5), it is not difficult to see that

(3.31)\begin{align} I& \le C_3\int_{\Omega}u^{2k(m-1)+2}+\varepsilon^{k(m-1)+1}\int_{\Omega}|\nabla \omega|^{2k(m-1)+2}{\rm d}x+C_4\varepsilon\int_{\Omega}|\nabla \omega|^2{\rm d}x \nonumber\\ & \le \frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x+\varepsilon^{k(m-1)+1}\int_{\Omega}|\nabla \omega|^{2k(m-1)+2}{\rm d}x +C_5. \end{align}

Noticing that $\omega$ is bounded, and using (3.3) gives

(3.32)\begin{equation} II\le \frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x+C_6\int_{\Omega}|\nabla v|^{\frac{4k(m-1)+2(m+1)}{m}} {\rm d}x+C_7.\end{equation}

Next, we estimate the second term $\displaystyle C_6\int _{\Omega }|\nabla v|^{\frac {4k(m-1)+2(m+1)}{m}} {\rm d}x$ in (3.32). From Gagliardo-Nirenberg inequality, we infer that

(3.33)\begin{equation} \|\nabla v\|_{L^p}^p\le C\|D^2 v\|_{L^{q_1}}^{p\alpha}\|\nabla v\|_{L^{q_2}}^{p(1-\alpha)}+C\|\nabla v\|_{L^2}^2, \ \text{with}\ \frac1p=\left(\frac1{q_1}-\frac1N\right)\alpha+\frac{1-\alpha}{q_2}.\end{equation}

By direct calculation, it is easy to obtain

\[ p\alpha< q_1 \Leftrightarrow \frac{p}{q_1}<1+\frac{q_2}N. \]

When $N=2$, take

\[ p=\frac{4k(m-1)+2(m+1)}{m}, \quad q_1=2k(m-1)+m+1, \quad q_2=2 \]

in (3.33). Then $p\alpha < q_1$. Using (3.3), we infer from (3.33) that for any small constant $\rho >0$,

(3.34)\begin{align} C_6\int_{\Omega}|\nabla v|^{\frac{4k(m-1)+2(m+1)}{m}} {\rm d}x & \le C_8\|D^2 v\|_{L^{{2k(m-1)+m+1}}}^{p\alpha}+C_9\nonumber\\ & \le \rho\int_{\Omega}|D^2 v|^{2k(m-1)+m+1} {\rm d}x+C_\rho, \quad \text{when $N=2$}, \end{align}

since $p\alpha < q_1$.

When $N\ge 3$. By (2.5) and (3.13), we see that

(3.35)\begin{equation} \sup_{0< t< T_{\max}}\|\nabla v\|_{L^r}^r\le C, \quad \forall r<\frac{N(m+1)}{(N+1-m)_+}.\end{equation}

We take

\begin{align*} & p=\frac{4k(m-1)+2(m+1)}{m}, \quad q_1=2k(m-1)+m+1, \\ & q_2=\frac{N(m+\frac{5}{2N})}{(N+1-m)_+} \quad (N\ge 3) \end{align*}

in (3.33). Then $\|\nabla v\|_{L^{q_2}}$ is uniformly bounded by (3.35), and

\[ p\alpha< q_1 \Leftrightarrow\frac{p}{q_1}<1+\frac{q_2}N\Leftrightarrow \frac{2}m<1+\frac{m+\frac{5}{2N}}{(N+1-m)_+} \quad (N\ge 3). \]

By direct verification, the above inequality $\frac {2}m<1+\frac {m+\frac {5}{2N}}{(N+1-m)_+}$ holds when $m>\frac {2N}{N+2}$. Then

(3.36)\begin{align} C_6\int_{\Omega}|\nabla v|^{\frac{4k(m-1)+2(m+1)}{m}} {\rm d}x\le \rho\int_{\Omega}|D^2 v|^{2k(m-1)+m+1} {\rm d}x+\tilde C_\rho, \quad \text{when }N\ge 3\end{align}

for any small constant $\rho >0$.

Substituting (3.34) and (3.36) into (3.32) yields

(3.37)\begin{align} II\le \frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x+\rho\int_{\Omega}|D^2 v|^{2k(m-1)+m+1} {\rm d}x+\hat C_\rho \end{align}

for any small constant $\rho >0$ with $\hat C_\rho$ depending on $\rho$.

By substituting the inequalities (3.30), (3.31) and (3.37) into (3.29) yields

(3.38)\begin{align} & \left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1} \frac{{\rm d}}{{\rm d}t}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x \nonumber\\ & \quad+\frac{2k+1}2\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\nonumber\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \nonumber\\ & \quad+\frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x \nonumber\\ & \qquad\le C_2\int_{\Omega}|\omega_t|^{2k(m-1)+2} {\rm d}x+\varepsilon^{k(m-1)+1}\int_{\Omega}|\nabla \omega|^{2k(m-1)+2}{\rm d}x\nonumber\\ & \qquad\quad +\rho\int_{\Omega}|D^2 v|^{2k(m-1)+m+1} {\rm d}x+\widehat C_{\rho}. \end{align}

Noticing that

\begin{align*} & \sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x\\ & \quad \le \frac{\mu}8 \int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x +C \end{align*}

and for any small constant $\rho >0$, there exists a constant $C_\rho$ such that

\begin{align*} \int_{\Omega}|v|^{2k(m-1)+2}{\rm d}x& \le \rho\int_{\Omega}|D^2v|^{2k(m-1)+2}{\rm d}x+C_\rho\int_{\Omega}|v|^2{\rm d}x\\ & \le \tilde \rho \int_{\Omega}|D^2v|^{2k(m-1)+m+1}{\rm d}x+\tilde C_\rho, \end{align*}

combining (3.21), and using $L^p$ theory of linear parabolic equations, from (3.38) we derive

\begin{align*} & \left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1}\frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1}\sup_{0< t< T_{\max}}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x \\ & \quad+\frac{2k+1}2\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x{\rm d}s \\ & \quad+\frac{3\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s \\ & \qquad\le C_{10}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}|\omega_t|^{2k(m-1)+2} {\rm d}x{\rm d}s+C_{11}\varepsilon^{k(m-1)+1}\\ & \quad \qquad \times \sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}|\nabla \omega|^{2k(m-1)+2}{\rm d}x{\rm d}s \\ & \quad\qquad+C_{12}\rho\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}|D^2 v|^{2k(m-1)+m+1} {\rm d}x{\rm d}s+M_{\rho} \\ & \qquad\le C_{13}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}\left(|u|^{2k(m-1)+2} +|v|^{2k(m-1)+2}+1\right){\rm d}x{\rm d}s \\ & \quad\qquad+C_{14}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}\left(|u|^{2k(m-1)} +|v|^{2k(m-1)}+1\right){\rm d}x{\rm d}s\\ & \qquad \quad +C_{12}\rho\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}|D^2 v|^{2k(m-1)+m+1} {\rm d}x{\rm d}s+M_{\rho} \\ & \qquad\le\frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s\\ & \qquad \quad +C_{15}\rho\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}|D^2 v|^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s+\widetilde M_{\rho} \\ & \qquad\le \frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s\\ & \qquad \quad +C_{16}\rho\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}|u|^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s+\widetilde M_{\rho}. \end{align*}

Take $\rho$ appropriately small in the above inequality, such that $C_{16}\rho <\frac {\mu }8(\frac {m}{m-1})^{2k+1}$. Then we finally arrive that

which implies (3.28).

By (3.28), there exists $k$ sufficiently large such that $(m-1)(2k+1)+1>N$, then by Neumann heat semigroup theory, it is easy to obtain the $L^\infty$ estimation of $\nabla v$. By $L^p$ theory of linear parabolic equations, we also have $W_p^{2,1}$ estimation of $v$ and $L^p$ estimation of $\omega _t$.

Lemma 3.7 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Assume that $m>\frac {2N}{N+2}$. Then

\[ \sup_{0< t< T_{\max}}\|v({\cdot}, t)\|_{W^{1,\infty}}< C, \]

and for any $p>1$,

\[ \sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\left(\|v_t\|_{L^p}^p+\|\omega_t\|_{L^p}^p+\|D^2v\|_{L^p}^p\right){\rm d}x{\rm d}s\le \tilde C, \]

where $C$, $\tilde C$ are independent of $\varepsilon$ and $T_{\max }$, it depends only on $k$, $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$.

In addition, by (3.28), and combining with Neumann heat semigroup theory, it is also easy to obtain the $L^\infty$ estimation of $\nabla \omega$ (depending on $\varepsilon$). Furthermore, by (3.28) and Lemma 3.7, using standard Morser iterative technique, the $L^\infty$ estimation of $u$ is also easy to obtain.

Lemma 3.8 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Assume that $m>\frac {2N}{N+2}$. Then

\[ \sup_{0< t< T_{\max}}(\|u({\cdot}, t)\|_{L^\infty}+\|\omega({\cdot}, t)\|_{W^{1,\infty}})< C_\varepsilon, \]

where $C_\varepsilon$ depends on $\varepsilon$, and is independent of $T_{\max }$.

Using lemmas 3.1, 3.7, and 3.8, the solution of the regularized problem (3.1) exists globally for any $\varepsilon >0$, and these estimations in lemma 3.2-Lemma 3.7 are independent of $\varepsilon$. To obtain the global solution to the original problem (1.6), we also need to make some a priori estimates for $u_t$.

Lemma 3.9 Let $(u,\, v,\, \omega )$ be the classical solution of (3.1) in $[0,\, T_{\max })$. Assume that $m>\frac {2N}{N+2}$. Then for any $1< r<2$, and any $T>0$,

(3.39)\begin{equation} \|u_t\|_{L^2\left((0, T);W^{{-}1,r}(\Omega)\right)}\le C_{T, r}, \end{equation}

where $C_{T, r}$ depends on $T$, $r$, and it is independent of $\varepsilon$.

Proof. For any $\varphi \in C^{\infty }(Q_T)$, we see that

\begin{align*} \int_0^T\int_{\Omega}u_t\varphi {\rm d}x{\rm d}t& ={-}\int_0^T\int_{\Omega} (u+\varepsilon) \nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right) \nabla\varphi {\rm d}x \\ & \quad+\chi\int_0^T\int_{\Omega}u\nabla v\nabla\varphi {\rm d}x{\rm d}t+\mu\int_0^T\int_{\Omega}u(1-u-\omega)\varphi {\rm d}x{\rm d}t. \end{align*}

Then for any $q>2$, using (3.13), (3.28) and lemma 3.7, we derive that

\begin{align*} & \left|\int_0^T\int_{\Omega}u_t\varphi {\rm d}x{\rm d}t\right|\le \int_0^T \left\|\sqrt{u+\varepsilon} \nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right\|_{L^2}\\ & \quad \times \|u+\varepsilon\|_{L^{\frac q{q-2}}}^{\frac12}\|\nabla\varphi\|_{L^q}{\rm d}t \\ & \quad+\chi \int_0^T \|u\|_{L^2}\|\nabla v\|_{L^\infty}\|\nabla\varphi\|_{L^2}{\rm d}t +\mu \int_0^T\|u(1-u-\omega)\|_{L^2(Q_T)}\|\varphi\|_{L^2}{\rm d}t \\ & \le C\left(\int_0^T\left\|\sqrt{u+\varepsilon} \nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right\|_{L^2}^2{\rm d}t\right)^{\frac12}\left(\int_0^T\|\nabla\varphi\|_{L^q}^2{\rm d}t\right)^{\frac12} \\ & \quad+\chi \left(\int_0^T \|u\|_{L^2}^2\|\nabla v\|_{L^\infty}^2{\rm d}t\right)^{\frac12}\left(\int_0^T\|\nabla\varphi\|_{L^2}^2{\rm d}t\right)^{\frac12}\\ & \quad +\mu \left(\int_0^T\|u(1-u-\omega)\|_{L^2}^2{\rm d}t\right)^{\frac12}\left(\int_0^T\|\varphi\|_{L^2}^2{\rm d}t\right)^{\frac12} \\ & \le C\left(\int_0^T\|\varphi\|_{W^{1,q}}^2{\rm d}t\right)^{\frac12}, \end{align*}

namely for any $q>2$,

\[ \|u_t\|_{L^2\left((0, T);W^{{-}1,\frac{q}{q-1}}(\Omega)\right)}\le C_q, \]

which implies (3.39).

Combining lemmas 3.1 -3.8, and noticing (3.16), we have the following proposition.

Proposition 3.10 Assume that $m>\frac {2N}{N+2}$. For any $\varepsilon \in (0,\,1)$, the problem (3.1) in $[0,\, T_{\max })$ admits a unique positive global classical solution $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )\in C^{2, 1}(Q)\cap C^0(\overline \Omega \times [0,\, \infty ))$ such that

(3.40)\begin{align} \sup_{0< t<\infty}\left(\|v_\varepsilon({\cdot}, t)\|_{W^{1,\infty}}+\|\omega_\varepsilon({\cdot}, t)\|_{L^{\infty}}+\|u_\varepsilon({\cdot}, t)\|_{L^{r}}\right)\le C_r,\quad \text{for any }r>1, \end{align}

and for any $T>0$,

(3.41)\begin{align} & \sup_{0< t< T}\int_\Omega \left(\frac{ |\nabla \omega_\varepsilon|^2}{\omega_\varepsilon} \right){\rm d}x+\int_0^T\int_\Omega \left((u_\varepsilon+\varepsilon)^{2m-3}|\nabla u_\varepsilon|^2\vphantom{\frac{ u_\varepsilon|\nabla \omega_\varepsilon|^2}{\omega_\varepsilon}}\right.\nonumber\\ & \quad +\left.\frac{ u_\varepsilon|\nabla \omega_\varepsilon|^2}{\omega_\varepsilon}+\varepsilon |D^2\omega_\varepsilon|^2 \right){\rm d}x{\rm d}t \le C_T, \end{align}

(3.42)\begin{align} & \int_0^T\int_\Omega\left(|D^2 v_\varepsilon|^r+|\partial_t v_\varepsilon|^r+|\partial_t \omega_\varepsilon|^r\right){\rm d}x{\rm d}t\le C_{Tr}, \quad \text{for any }r>1, \end{align}

(3.43)\begin{align} & \|\nabla(u_\varepsilon+\varepsilon)^m\|_{L^p(Q_T)}+ \|\partial_t u_\varepsilon\|_{L^2\left((0, T);W^{{-}1,p}(\Omega)\right)}\le C_{Tp}, \quad \text{for any }1< p<2, \end{align}

where $C_r$, $C_T$, $C_{Tp}$, and $C_{Tr}$ are independent of $\varepsilon$, they depend on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$, $C_r$ depends on $r$, $C_T$ depends on $T$, $C_{Tp}$ depends on $T$, $p$, and $C_{Tr}$ depends on $T$, $r$.

Proof of theorem 1.1 Since $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )$ are classical solutions of (3.1), they obviously satisfy the equations (3.1) in the sense of distribution. By proposition 3.10, using Aubin-Lions lemma, and Sobolev compact embedding theorem, there exists a subsequence of $\{(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )\}$, which we still remember as itself for convenience, such that as $\varepsilon \to 0$,

\begin{align*} & u_\varepsilon, \omega_\varepsilon \rightarrow u, \omega\ \text{in }L^r(Q_T)\text{ for any }r>1, \\ & \nabla (u_\varepsilon+\varepsilon)^m \rightharpoonup \nabla u^m, \text{ in }L^p(Q_T) \text{ for any }1< p<2, \\ & \varepsilon\Delta\omega\rightarrow 0, \ \text{in }L^2(Q_T), \\ & v_\varepsilon \rightharpoonup v, \ \text{in }W_r^{2,1}(Q_T)\text{ for any }r>1, \\ & \nabla v_\varepsilon \rightarrow \nabla v, \ \text{in }L^r(Q_T)\text{ for any }r>1, \\ & v_\varepsilon \rightarrow v, \ \text{uniformly}. \end{align*}

Letting $\varepsilon \to 0$ and utilizing the above convergence results, we can ultimately verify that $(u,\, v,\, \omega )$ satisfies definition 2.1, that is, $(u,\, v,\, \omega )\in \mathcal {X}_1\times \mathcal {X}_2\times \mathcal {X}_3$ is a weak solution. Theorem 1.1 is proved.

4. Global solution to the haptotaxis-only system

In this section, we consider the haptotaxis-only system. Similar to § 2, to prove the existence of weak solutions, we consider the following regularized problem

(4.1)\begin{equation} \begin{cases} u_t=m\nabla\cdot (u+\varepsilon)^{m-1}\nabla u-\xi\nabla\cdot(u\nabla\omega)+\mu u(1-u-\omega),\quad \text{in}\ Q, \\ v_t=\Delta v+u-v, \quad \text{in}\ Q, \\ \omega_t=\varepsilon\Delta\omega -v\omega+\eta \omega(1-u-\omega), \ \text{in}\ Q, \\ \displaystyle\left.\frac{\partial u}{\partial{\bf n}}\right|_{\partial\Omega}=0, \left.\frac{\partial v}{\partial{\bf n}}\right|_{\partial\Omega}=0, \left.\frac{\partial \omega}{\partial{\bf n}}\right|_{\partial\Omega}=0, \\ u(x,0)=u_{\varepsilon 0}(x), \quad v(x,0)=v_{\varepsilon 0}(x), \quad \omega(x,0)=\omega_{\varepsilon 0}(x), \quad x\in\Omega. \end{cases} \end{equation}

Similar to § 3, in the subsequent energy estimate calculations, we omit the subscript $\varepsilon$ of the approximate solutions $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )$ in the subsequent energy estimate calculations.

It is clear that lemmas 3.1, 3.2, and 3.3 also hold for $\chi =0$. Similar to lemma 3.4, we also have

Lemma 4.1 Let $(u,\, v,\, \omega )$ be the classical solution of (4.1) in $[0,\, T_{\max })$. Then

(4.2)\begin{align} & \sup_{0< t< T_{\max}}\int_{\Omega} u^m {\rm d}x +\sup_{\tau< t< T_{\max}}\int_{t-\tau}^{t}\int_{\Omega}\left((u+\varepsilon)\vphantom{\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2}\right.\nonumber\\ & \quad \times \left.\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2 +u^{m+1}\right){\rm d}x{\rm d}s\le C, \end{align}

where $C$ is independent of $\varepsilon$ and $T_{\max }$, it depends only on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$. In addition, we also have

(4.3)\begin{equation} \int_0^T\int_{\Omega} (u+\varepsilon)^{2m-3} |\nabla u|^2{\rm d}x{\rm d}t\le C_T, \end{equation}

and for any $1< p<2$,

(4.4)\begin{equation} \int_0^T\int_{\Omega}|\nabla(u+\varepsilon)^{m}|^p{\rm d}x{\rm d}t\le C_{Tp}, \end{equation}

where $C_T$, $C_{Tp}$ are independent of $\varepsilon$, and $C_T$ depends on $T$, $C_{Tp}$ depends on $T$ and $p$.

Proof. Similar to the proof of (3.18), we see that

\begin{align*} & \frac{{\rm d}}{{\rm d}t}\int_{\Omega}\left(\frac1{m-1}(u+\varepsilon)^m -\xi u\omega\right) {\rm d}x +\int_{\Omega}(u+\varepsilon)\left|\nabla\left(\frac{m}{m-1} (u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \quad={-}\varepsilon\xi\int_{\Omega}\nabla\omega\nabla \left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right){\rm d}x\\ & \qquad +\int_{\Omega}\mu u(1-u-\omega)\left(\frac{m}{m-1}u^{m-1}-\xi\omega\right) {\rm d}x \\ & \qquad+\xi\int_{\Omega}u(-\varepsilon\Delta\omega+v\omega-\eta\omega(1-u-\omega)){\rm d}x \\ & \quad\le\frac{\varepsilon\xi^2\|\omega\|_{L^\infty}}2\int_{\Omega}\frac{|\nabla\omega|^2}{\omega} {\rm d}x+\frac{\varepsilon}2\int_{\Omega}\left|\nabla \left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \qquad -\frac{3m\mu}{4(m-1)}\int_{\Omega}u^{m+1}{\rm d}x +\varepsilon^2\int_{\Omega}|\Delta\omega|^2{\rm d}x +C \\ & \quad\le \frac{\varepsilon}2\int_{\Omega}\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x -\frac{3m\mu}{4(m-1)}\int_{\Omega}u^{m+1}{\rm d}x \\ & \qquad +N\varepsilon^2\int_{\Omega}|D^2\omega|^2{\rm d}x+\tilde C, \end{align*}

that is

(4.5)\begin{align} & \frac{{\rm d}}{{\rm d}t}\int_{\Omega}\left(\frac1{m-1}(u+\varepsilon)^m -\xi u\omega\right) {\rm d}x +\frac12\int_{\Omega}(u+\varepsilon)\left|\nabla\left(\!\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\!\right)\right|^2{\rm d}x \nonumber\\ & \quad+\frac{3m\mu}{4(m-1)}\int_{\Omega}u^{m+1}{\rm d}x\le N\varepsilon^2\int_{\Omega}|D^2\omega|^2{\rm d}x+\tilde C. \end{align}

Then completely similar to the proof of lemma 3.4, we complete the proof.

Completely similar to lemma 3.5, we also have that

Lemma 4.2 Let $(u,\, v,\, \omega )$ be the classical solution of (4.1) in $[0,\, T_{\max })$. Then for any $r>2$,

(4.6)\begin{align} & \sup_{0< t< T_{\max}}\int_\Omega |\nabla \omega|^{r} {\rm d}x+ \sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} \varepsilon|\nabla \omega|^{r+2} {\rm d}x{\rm d}s\nonumber\\ & \quad \le \frac{\tilde C}{\varepsilon^{\frac r2}}\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u^r+v^r+1){\rm d}x{\rm d}s, \end{align}

where $\tilde C$ is independent of $\varepsilon$ and $T_{\max }$, it depends only on $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $r$, $\Omega$.

Then similar to lemma 3.6, we also have that

Lemma 4.3 Let $(u,\, v,\, \omega )$ be the classical solution of (4.1) in $[0,\, T_{\max })$. Then for any positive integer $k$, we have that

(4.7)\begin{align} & \sup_{0< t< T_{\max}}\int_{\Omega}(u+\varepsilon)^{(m-1)(2k+1)+1} {\rm d}x+\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x{\rm d}s \nonumber\\ & \quad+\sup_{\tau< t< T_{\max}}\int_{t-\tau}^t\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\nonumber\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x{\rm d}s\le C, \end{align}

where $C$ is independent of $\varepsilon$ and $T_{\max }$, it depends only on $k$, $\chi$, $\xi$, $m$, $\mu$, $\eta$, $u_0$, $v_0$, $\omega _0$, $\Omega$.

Proof. For any given positive integer $k$, multiplying both sides of the first equation of (4.1) by $(\frac {m}{m-1}(u+\varepsilon )^{m-1}-\xi \omega )^{2k+1}$, and noticing that $\omega$ is bounded, then we arrive at

\begin{align*} & \int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k+1}u_t {\rm d}x \\ & \quad+(2k+1)\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \qquad={-}\varepsilon\xi(2k+1)\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\\ & \quad \qquad \times \nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\nabla \omega {\rm d}x \\ & \quad \qquad+\int_{\Omega}\mu u(1-u-\omega)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k+1} {\rm d}x \\ & \qquad\le \frac{2k+1}2\int_{\Omega}\varepsilon\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x \\ & \quad \qquad+\frac{\varepsilon\xi^2(2k+1)}2\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}|\nabla \omega|^2{\rm d}x \\ & \qquad \quad-\frac{3\mu}4\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x+C_1. \end{align*}

Combining with (3.30), and using (3.31) yields

(4.8)\begin{align} & \left(\frac{m}{m-1}\right)^{2k+1}\sum_{i=0}^{2k+1} \frac{(-\frac{m-1}{m}\xi)^iC_{2k+1}^i}{(m-1)(2k+1-i)+1} \frac{{\rm d}}{{\rm d}t}\int_{\Omega} (u+\varepsilon)^{(m-1)(2k+1-i)+1}\omega^i {\rm d}x \nonumber\\ & \quad+\frac{2k+1}2\int_{\Omega}(u+\varepsilon)\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}\nonumber\\ & \quad \times \left|\nabla\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)\right|^2{\rm d}x +\frac{5\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x \nonumber\\ & \qquad\le \frac{\varepsilon\xi^2(2k+1)}2\int_{\Omega}\left(\frac{m}{m-1}(u+\varepsilon)^{m-1}-\xi\omega\right)^{2k}|\nabla \omega|^2{\rm d}x\nonumber\\ & \qquad \quad +C_2\int_{\Omega} |\omega_t|^{2k(m-1)+2} {\rm d}x +C_3 \nonumber\\ & \qquad\le \frac{\mu}8\left(\frac{m}{m-1}\right)^{2k+1}\int_{\Omega} u^{(2k+1)(m-1)+2} {\rm d}x+\varepsilon^{k(m-1)+1}\int_{\Omega}|\nabla \omega|^{2k(m-1)+2}{\rm d}x\nonumber\\ & \qquad \quad +C_2\int_{\Omega} |\omega_t|^{2k(m-1)+2} {\rm d}x+C_4. \end{align}

By (4.6), (2.5), and using $L^p$ theory of linear parabolic equations, we arrive at

Then (4.7) is proved.

Completely similar to lemmas 3.7– 3.9. We have the following results.

Proposition 4.4 For any $\varepsilon \in (0,\,1)$, the problem (4.1) admits a unique positive global classical solution $(u_\varepsilon,\, v_\varepsilon,\, \omega _\varepsilon )\in C^{2, 1}(Q)\cap C^0(\overline \Omega \times [0,\, \infty ))$ such that

(4.9)\begin{equation} \sup_{0< t<\infty}\left(\|v_\varepsilon({\cdot}, t)\|_{W^{1,\infty}}+\|\omega_\varepsilon({\cdot}, t)\|_{L^{\infty}}+\|u_\varepsilon({\cdot}, t)\|_{L^{r}}\right)\le C_r,\quad \text{for any }r>1, \end{equation}

and for any $T>0$,

(4.10)\begin{align} & \sup_{0< t< T}\int_\Omega \left(\frac{ |\nabla \omega_\varepsilon|^2}{\omega_\varepsilon} \right){\rm d}x+\int_0^T\int_\Omega \left((u_\varepsilon+\varepsilon)^{2m-3}|\nabla u_\varepsilon|^2\vphantom{\frac{ u_\varepsilon|\nabla \omega_\varepsilon|^2}{\omega_\varepsilon}}\right.\nonumber\\ & \quad +\left.\frac{ u_\varepsilon|\nabla \omega_\varepsilon|^2}{\omega_\varepsilon}+\varepsilon |D^2\omega_\varepsilon|^2 \right){\rm d}x{\rm d}t \le C_T, \end{align}

(4.11)\begin{align} & \int_0^T\int_\Omega\left(|D^2 v_\varepsilon|^r+|\partial_t v_\varepsilon|^r+|\partial_t \omega_\varepsilon|^r\right){\rm d}x{\rm d}t\le C_{Tr}, \quad \text{for any }r>1, \end{align}

(4.12)\begin{align} & \|\nabla(u_\varepsilon+\varepsilon)^m\|_{L^p(Q_T)}+ \|\partial_t u_\varepsilon\|_{L^2\left((0, T);W^{{-}1,p}(\Omega)\right)}\le C_{Tp}, \quad \text{for any }1< p<2, \end{align}

Letting $\varepsilon \to 0$, we complete the proof of Theorem 1.2.

Acknowledgements

The author is supported by NSFC(12271186,12171166).

Statements and declarations

The author declares that there is no competing interest, and confirms that the data supporting the findings of this study are available within the paper.

References

Byrne, H. M., Chaplain, M. A. J., J.Pettet, G. and Mcelwain, D. L. S.. Mathematical model of trophoblast invasion. J. Theor. Med. 1 (1999), 275–286.Google Scholar

Byrne, H. M., Chaplain, M. A. J., J.Pettet, G. and Mcelwain, D. L. S.. An analysis of a mathematical model of trophoblast invasion. Appl. Math. Lett. 14 (2001), 1005–1010.CrossRef Google Scholar

Chaplain, M. A. J. and Anderson, A. R. A., Mathematical modelling of tissue invasion, In: Preziosi, L., Ed., Cancer Modelling and Simulation (Chapman Hall/CRC, 2003), pp. 269–297.CrossRef Google Scholar

Chaplain, M. A. J. and Lolas, G.. Mathematical modelling of cancer invasion of tissue: dynamic heterogeneity. Netw. Heterogeneous Media 1 (2006), 399–439.CrossRef Google Scholar

Cao, X.. Boundedness in a three-dimensional chemotaxis-haptotaxis model. Z. Angew. Math. Phys. 67 (2016), 11.CrossRef Google Scholar

Hillen, T., Painter, K. J. and Winkler, M.. Convergence of a cancer invasion model to a logistic chemotaxis model. Math. Models Methods Appl. Sci. 25 (2013), 165–198.CrossRef Google Scholar

Jin, C.. Boundedness and global solvability to a chemotaxis-haptotaxis model with slow and fast diffusion. Discrete Contin. Dyn. Syst. Ser. B 23 (2018), 1675–1688.Google Scholar

Jin, C.. Global classical solutions and convergence to a mathematical model for cancer cells invasion and metastatic spread. J. Differ. Equ. 269 (2020), 3987–4021.CrossRef Google Scholar

Jin, C.. Global bounded weak solutions and asymptotic behavior to a chemotaxis-Stokes model with non-Newtonian filtration slow diffusion. J. Differ. Equ. 28 (2021), 148–184.CrossRef Google Scholar

Jin, C.. Global classical solution and boundedness to a chemotaxis-haptotaxis model with re-establishment mechanisms. Bull. Lond. Math. Soc. 50 (2018), 598–618.CrossRef Google Scholar

Jin, C.. Global strong solution and periodic dynamic behavior to Chaplain-Lolas model. J. Dyna. Diff. Eqns. (2022). doi:10.1007/s10884-022-10210-wGoogle Scholar

Li, Y. and Lankeit, J.. Boundedness in a chemotaxis–Chaptotaxis model with nonlinear diffusion. Nonlinearity 29 (2016), 1564–1595.CrossRef Google Scholar

Marciniak-Czochra, A. and Ptashnyk, M.. Boundedness of solutions of a haptotaxis model. Math. Models Methods Appl. Sci. 20 (), –.Google Scholar

Mizoguchi, N. and Souplet, P.. Nondegeneracy of blow-up points for the parabolic Keller-Segel system. Ann. I. H. Poincaré - AN 31 (2014), 851–875.Google Scholar

Pang, P. H. and Wang, Y.. Global existence of a two-dimensional chemotaxis-haptotaxis model with remodeling of non-diffusible attractant. J. Differ. Equ. 263 (2017), 1269–1292.CrossRef Google Scholar

Pang, P. H. and Wang, Y.. Global boundedness of solutions to a chemotaxis-haptotaxis model with tissue remodeling. Math. Models Methods Appl. Sci. 28 (2018), 2211–2235.CrossRef Google Scholar

Stinner, C. and Surulescu, C.. A multiscale model for pH-tactic invasion with time-varying carrying capacities. IMA J. Appl. Math. 80 (2015), 1300–1321.CrossRef Google Scholar

Tao, Y. and Wang, M.. Global solution for a chemotactic-haptotactic model of cancer invasion. Nonlinearity 21 (2008), 2221–2238.CrossRef Google Scholar

Tao, Y. and Winkler, M.. A chemotaxis-haptotaxis model: the roles of nonlinear diffusion and logistic source. SIAM J. Math. Anal. 43 (2011), 685–704.CrossRef Google Scholar

Tao, Y. and Winkler, M.. Large time behavior in a multidimensional chemotaxis-haptotaxis model with slow signal diffusion. SIAM J. Math. Anal. 47 (2015), 4229–4250.CrossRef Google Scholar

Wang, Y.. Boundedness in the higher-dimensional chemotaxis-haptotaxis model with nonlinear diffusion. J. Differ. Equ. 260 (2016), 1975–1989.CrossRef Google Scholar

Wang, Y. and Ke, Y.. Large time behavior of solution to a fully parabolic chemotaxis-haptotaxis model in higher dimensions. J. Differ. Equ. 260 (2016), 6960–6988.CrossRef Google Scholar

Walker, C. and Webb, G. F.. Global existence of classical solutions for a haptotaxis model. SIAM J. Math. Anal. 38 (2007), 1694–1713.CrossRef Google Scholar

Zheng, J.. Boundedness of solutions to a quasilinear higher-dimensional chemotaxis-haptotaxis model with nonlinear diffusion. Discrete Contin. Dyn. Syst. 37 (2017), 627–643.CrossRef Google Scholar

Article contents

The roles of nonlinear diffusion, haptotaxis and ECM remodelling in determining the global solvability of a cancer invasion model

Abstract

Keywords

MSC classification

1. Introduction

2. Preliminaries

3. Uniform energy estimations and global solvability of chemotaxis-haptotaxis system

4. Global solution to the haptotaxis-only system

Acknowledgements

Statements and declarations

References

Save article to Kindle

Save article to Dropbox

Save article to Google Drive

Reply to: Submit a response

Your details

You have entered the maximum number of contributors

Conflicting interests