Global existence and boundedness in a quasilinear attraction-repulsion chemotaxis system of parabolic-elliptic type

Abstract

This paper deals with the global existence and boundedness of solutions to the following quasilinear attraction-repulsion chemotaxis system:

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} u_{t}=\nabla\cdot(D(u)\nabla u)-\nabla\cdot(\chi u\nabla v)+\nabla\cdot (\xi u\nabla w), & x\in\Omega, t>0,\\ 0=\Delta v+\alpha u-\beta v,& x\in\Omega, t>0,\\ 0=\Delta w+\gamma u-\delta w,& x\in\Omega, t>0, \end{array}\displaystyle \right . $$

under homogeneous Neumann boundary conditions in a bounded domain $\Omega\subset\mathbb{R}^{n}$ ($n\geq2$) with smooth boundary, where $D(u)\geq c_{D} u^{m-1}$ with $m\geq1$ and some constant $c_{D}>0$. It is proved that if $\xi\gamma-\chi\alpha>0$ or $m>2-\frac{2}{n}$, then for any sufficiently regular initial data, this system possesses a unique global bounded classical solution for the case of nondegenerate diffusion (i.e., $D(u)>0$ for all $u\geq 0$), whereas for the case of degenerate diffusion (i.e., $D(u)\geq0$ for all $u\geq0$), it is shown that there exists a global bounded weak solution under the same assumptions.

1 Introduction

Chemotaxis is widespread in nature. It describes the oriented migration of cells or bacteria toward the concentration gradient of a chemical substance. In 1970, Keller and Segel [1] derived the well-known and widely studied Keller-Segel attractive model. The most obvious feature of this system is that the solution may blow up in finite time (see [2–9] and references therein). Hillen and Painter [5] suggested the chemotaxis model with nonlinear diffusion and aggregation by considering the volume-filling effect. Therefore, there are many papers on the global existence or finite time blow-up of solutions (e.g., see [10–23]).

In many biological processes, the migration of cells or bacteria is generally influenced by a combination of attractive and repulsive chemicals [24, 25]. The scholars in [26, 27] have proposed the corresponding attraction-repulsion chemotaxis model

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} u_{t}=\Delta u-\nabla\cdot(\chi u\nabla v)+\nabla\cdot(\xi u\nabla w), & x\in\Omega, t>0,\\ \tau v_{t}=\Delta v+\alpha u-\beta v,& x\in\Omega, t>0,\\ \tau w_{t}=\Delta w+\gamma u-\delta w,& x\in\Omega, t>0, \end{array}\displaystyle \right . $$

(1.1)

under no-flux boundary conditions, where $\Omega\subset\mathbb{R}^{n}$ is a bounded domain with smooth boundary. Here $\chi\geq0$, $\xi\geq0$, $\alpha>0$, $\beta>0$, $\gamma>0$, $\delta>0$, and $\tau=0, 1$ are parameters. The unknown functions $u(x,t)$, $v(x,t)$, and $w(x,t)$ denote the cell density, the concentration of an attractive signal, and the concentration of a repulsive chemical, respectively. If we take $\xi =0$, then model (1.1) is the classical attractive Keller-Segel model. The first cross-diffusive term and the second in the first equation of (1.1) mean that the movement of the bacteria is directed toward the increasing concentration of an attractive substance and away from the increasing concentration of a repulsive chemical, respectively. The second and third equations in model (1.1) indicate that chemoattractant and chemorepellent are produced by cells and have attenuation. There are fewer results for (1.1) than the classical attractive Keller-Segel model, mainly since the latter possesses a useful Lyapunov functional whereas the former does not admit such a functional. When $n=1$ and $\tau=1$, the global existence and asymptotic dynamics of solutions of (1.1) were studied by [28, 29]. When $n=2$, $\tau=1$, and $\xi\gamma-\chi\alpha>0$, the model (1.1) possesses a unique global bounded classical solution with any sufficient regular initial data (see [30, 31]). When $n=2 \mbox{ or } 3$, $\tau=1$, and $\xi\gamma=\chi\alpha$, Lin et al. [32] proved that (1.1) admits a unique global bounded classical solution, and large time-behavior is considered. When $\tau=0$, the global solvability, critical mass phenomenon, blow-up, and asymptotic behavior were studied in [33, 34]. Recently, Jin and Wang [35] studied the boundedness, blow-up, and critical mass phenomenon of solutions to a variant of (1.1) for $n=2$. Liu et al. [36] also studied the pattern formation of model (1.1) with $\tau=1$ from both analytical and numerical aspects.

To the best of our knowledge, presently, there is no rigorous result on the attraction-repulsion chemotaxis model with nonlinear diffusion. Thus, this paper mainly aims to understand the competition among the repulsion, the attraction, and the nonlinear diffusion. Precisely, we will consider the global existence and boundedness of solutions to the following quasilinear attraction-repulsion chemotaxis system of parabolic-elliptic type:

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} u_{t}=\nabla\cdot(D(u)\nabla u)-\nabla\cdot(\chi u\nabla v)+\nabla\cdot (\xi u\nabla w), & x\in\Omega, t>0,\\ 0=\Delta v+\alpha u-\beta v, & x\in\Omega, t>0,\\ 0=\Delta w+\gamma u-\delta w, & x\in\Omega, t>0,\\ \frac{\partial u}{\partial\nu}=\frac{\partial v}{\partial\nu}=\frac {\partial w}{\partial\nu}=0, & x\in\partial\Omega, t>0,\\ u(x,0)=u_{0}(x),& x\in\Omega, \end{array}\displaystyle \right . $$

(1.2)

where $\Omega\subset\mathbb{R}^{n}$ ($n\geq2$) is a bounded domain with smooth boundary ∂Ω, and $\partial/\partial\nu$ represents the derivative with respect to the outer normal of ∂Ω. As usual, we assume that $\chi, \xi\ge0$ and that α, β, γ, and δ are positive parameters. For the diffusion coefficient D, we assume that

$$ D \in C^{2} \bigl([0,\infty) \bigr) $$

(1.3)

and there exist some constants $c_{D}>0$ and $m\ge1$ such that

$$ D(u)\geq c_{D} u^{m-1} \quad\mbox{for all } u \geq0. $$

(1.4)

In addition to (1.3) and (1.4), we will require that $D(u)$ satisfy

$$ D(u)>0 \quad\mbox{for all } u\geq0 $$

(1.5)

in some places. In particular, when $D(u)$ does not satisfy (1.5) (i.e., $D(u)\geq0$ for all $u\geq0$), equation (1.2)₁ may be degenerate at $u=0$.

We will show that we can allow for the case of attraction dominating the repulsion (i.e., $\xi\gamma-\chi\alpha<0$) and still obtain global existence results due to the nonlinear diffusion. Thus, our results confirm that the attraction-repulsion system with nonlinear diffusion can prevent blow-up of solutions in higher dimensions as mentioned before.

We now state the main results of this paper.

Theorem 1.1

Let $\Omega\subset\mathbb{R}^{n}$ ($n \geq2$) be a bounded domain with smooth boundary. Assume that $u_{0}\in W^{1,\infty}(\Omega)$ is a nonnegative function and $D(u)$ satisfies (1.3), (1.4), and (1.5). Suppose that

$$\xi\gamma-\chi\alpha>0 \quad\textit{or}\quad \xi\gamma-\chi\alpha\le0\textit{ and }m>2-\frac{2}{n}. $$

Then there exists a unique nonnegative bounded solution $(u,v,w)$ belonging to $C^{0} (\overline{\Omega} \times[0,\infty) )\cap C^{2,1} (\overline{\Omega} \times(0,\infty) )$ that solves system (1.2) classically.

Remark 1.1

Theorem 1.1 shows that the solution is still global, provided that the diffusion is strong enough even if the attraction prevails over the repulsion, which provides a supplement to the dichotomy boundedness vs. blow-up in attraction-repulsion chemotaxis equations of parabolic-elliptic type with nonlinear diffusion.

Remark 1.2

For $n=2$, Theorem 1.1 also shows that both the attraction and repulsion cannot result in blow-up when the linear diffusion is replaced by a nonlinear one.

For the case of $D(u)$ only fulfilling (1.3) and (1.4), since equation (1.2)₁ with $m >1$ may be degenerate at $u =0$, system (1.2) does not admit classical solutions in general as the porous medium equation does. However, we can prove that system (1.2) in this case possesses at least one nonnegative global bounded solution $(u, v, w)$ in the following weak sense.

Definition 1.1

Let $T>0$. Then a triple of nonnegative functions $(u,v,w)$ defined on $\Omega\times(0,T)$ is called a weak solution to (1.2) if

• (1) $u\in L^{\infty} ((0,T);L^{\infty}(\Omega) )$ and $D(u)\nabla u \in L^{2}_{\mathrm{loc}} ((0,T);L^{2} (\Omega) )$,

• (2) $v\in L^{\infty} ((0,T);W^{1,\infty}(\Omega) )$ and $w\in L^{\infty} ((0,T);W^{1,\infty}(\Omega) )$,

• (3) $(u,v,w)$ satisfies (1.2) in the distributional sense, that is, for every $\varphi\in C_{0} ^{\infty} (\Omega\times[0,T) )$,

  $$\begin{aligned}& \begin{aligned}[b] &\int_{0} ^{T} \int_{\Omega} \bigl(D(u)\nabla u\cdot\nabla\varphi-\chi u \nabla v \cdot\nabla\varphi+\xi u \nabla w\cdot\nabla\varphi-u\varphi _{t} \bigr)\,dx\,dt \\ &\quad= \int_{\Omega}u_{0}(x)\varphi(x,0)\,dx, \end{aligned}\\& \int_{0} ^{T} \int_{\Omega} (\nabla v\cdot\nabla\varphi+\beta v \varphi )\,dx\,dt= \int_{0} ^{T} \int_{\Omega} \alpha u \varphi \,dx\,dt, \\& \int_{0} ^{T} \int_{\Omega} (\nabla w\cdot\nabla\varphi+\delta w \varphi )\,dx\,dt= \int_{0} ^{T} \int_{\Omega} \gamma u \varphi \,dx\,dt. \end{aligned}$$

If $(u,v,w)$ is a weak solution to (1.2) on $\Omega\times(0,T)$ for all $T\in(0,\infty)$, then $(u,v,w)$ is called a global weak solution to (1.2).

Theorem 1.2

Let $\Omega\subset\mathbb{R}^{n}$ ($n \geq2$) be a bounded domain with smooth boundary. Assume that $u_{0}\in W^{1,\infty}(\Omega)$ is a nonnegative function and that $D(u)$ satisfies (1.3) and (1.4). Suppose that

$$\xi\gamma-\chi\alpha>0 \quad\textit{or}\quad \xi\gamma-\chi\alpha\le0\textit{ and }m>2-\frac{2}{n}. $$

Then there exists at least one nonnegative global weak solution $(u,v,w)$ to system (1.2). Moreover, $(u,v,w)$ satisfies

$$\bigl\| u(\cdot, t)\bigr\| _{L^{\infty}(\Omega)}+\bigl\| v(\cdot, t)\bigr\| _{L^{\infty}(\Omega)} +\bigl\| w(\cdot, t) \bigr\| _{L^{\infty}(\Omega)}\leq C \quad\textit{for all } t>0, $$

where $C>0$ is a constant independent of t.

The rest of this paper is organized as follows. In Section 2, we first prove the local existence and uniqueness of a solution to system (1.2) and then give mass estimates. In Section 3, we give some fundamental estimates for the solution $(u,v,w)$ to system (1.2) and then prove Theorem 1.1. In Section 4, we establish the existence of global bounded weak solutions to system (1.2).

2 Preliminaries

In this section, we first state the local well-posedness of system (1.2) and then give the mass estimates.

Lemma 2.1

Assume that $u_{0}\in W^{1,\infty}(\Omega)$ is a nonnegative function and D satisfies (1.3), (1.4), and (1.5). Then there exist $T_{\max}\in(0,\infty]$ and a unique triple $(u, v,w)$ of nonnegative functions from $C^{0} (\overline{\Omega }\times[0,T_{\max}) )\cap C^{2,1} (\overline{\Omega} \times(0, T_{\max}) )$ solving (1.2) classically in $\Omega\times (0,T_{\max})$. Moreover,

$$ \textit{if }T_{\max} < \infty, \quad\textit{then } \bigl\| u( \cdot ,t)\bigr\| _{L^{\infty}(\Omega)}\rightarrow\infty \textit{ as } t\rightarrow T_{\max}. $$

(2.1)

Proof

(i) Existence. Let $T\in(0,1)$, which is specified below. We define

$$\mathcal{S}_{T}:= \bigl\{ u\in\mathcal{X} | \bigl\| u(\cdot,t)\bigr\| _{L^{\infty}(\Omega)}\leq\| u_{0}\|_{L^{\infty}(\Omega)}+1=:R \mbox{ for all }t \in[0,T] \bigr\} , $$

which is a bounded closed convex subset of space $\mathcal{X}:=C^{0} (\bar{\Omega} \times[0,T] )$.

For any given $\tilde{u}\in S_{T}$, there exists a unique $(v,w)$ such that v and w solve the following elliptic equations

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} -\Delta v +\beta v=\alpha\tilde{u},& x\in\Omega, t\in(0,T),\\ \frac{\partial v}{\partial\nu}=0, & x\in\partial\Omega, t\in(0,T), \end{array}\displaystyle \right . $$

(2.2)

and

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} -\Delta w +\delta w=\gamma\tilde{u},& x\in\Omega, t\in(0,T),\\ \frac{\partial w}{\partial\nu}=0, & x\in\partial\Omega, t\in(0,T), \end{array}\displaystyle \right . $$

(2.3)

respectively. Then we can find a unique u solving the following parabolic equation:

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} u_{t}=\nabla\cdot (D(\tilde{u})\nabla u )+\nabla\cdot [(-\chi \nabla v+\xi\nabla w)u ],& x\in\Omega, t\in(0,T),\\ \frac{\partial u}{\partial\nu}=0, & x\in\partial\Omega, t\in (0,T),\\ u(x,0)= u_{0}(x), &x\in\Omega. \end{array}\displaystyle \right . $$

(2.4)

Thus, we can introduce a mapping $\Phi: \tilde{u}(\in\mathcal {S}_{T})\longmapsto u$ by defining $\Phi(\tilde{u})=u$.

We next show that Φ has a fixed point for T sufficiently small. The elliptic regularity [37], Theorem 8.34, implies that (2.2) admits a unique solution $v(\cdot,t)\in C^{1+\theta}(\Omega)$ for some $\theta\in(0,1)$. Similarly, (2.3) also possesses a unique solution $w(\cdot,t)\in C^{1+\theta}(\Omega)$. Moreover, the Sobolev embedding theorem and the $L^{p}$ estimates yield that

$$\|\nabla v\|_{L^{\infty} (\Omega\times(0,T) )}\leq C_{1}\|v\| _{L^{\infty} ((0,T);W^{2,p}(\Omega) )}\leq C_{2}\|\tilde{u}\| _{L^{\infty} ((0,T);L^{p}(\Omega) )} $$

and

$$\|\nabla w\|_{L^{\infty} (\Omega\times(0,T) )}\leq C_{1}\|w\| _{L^{\infty} ((0,T);W^{2,p}(\Omega) )}\leq C_{2}\| \tilde{u}\| _{L^{\infty} ((0,T);L^{p}(\Omega) )} $$

with $p>n$ and some constants $C_{1}>0$ and $C_{2}>0$. It then follows from [38], Theorem 6.1, that $u\in C^{\theta,\frac{\theta}{2}} (\Omega \times(0,T) )$ with

$$ \| u \|_{C^{\theta,\frac{\theta}{2}} (\Omega\times(0,T) )}\leq C_{3} $$

(2.5)

for some $\theta\in(0,1)$ and $C_{3}>0$, where $C_{3}$ depends on $\min_{0\leq s\leq R}D(s)$, $\|\nabla v\|_{L^{\infty}((0,T);C^{\theta}(\bar {\Omega}))}$ and $\|\nabla w\|_{L^{\infty}((0,T);C^{\theta}(\bar{\Omega }))}$. Thus, we obtain

$$\bigl\| u(\cdot,t)\bigr\| _{L^{\infty}(\Omega)} \le\| u_{0}\|_{L^{\infty}(\Omega)}+\bigl\| u( \cdot,t)-u_{0}\bigr\| _{L^{\infty}(\Omega)} \leq\| u_{0} \|_{L^{\infty}(\Omega)} +C_{3}t^{\frac{\theta}{2}}. $$

Hence if we take $T<(\frac{1}{C_{3}})^{\frac{2}{\theta}}$, then

$$ \bigl\| u(\cdot, t)\bigr\| _{L^{\infty}(\Omega)}\leq\|u_{0} \|_{L^{\infty}(\Omega )}+1=R \quad\mbox{for all } t\in[0,T], $$

(2.6)

which implies that $u\in\mathcal{S}_{T}$. Then we conclude that $\Phi (\mathcal{S}_{T})\subset\mathcal{S}_{T}$ and $\Phi(\mathcal{S}_{T})$ is compact in $\mathcal{S}_{T}$ by (2.5). Moreover, we can easily deduce that Φ is a continuous operator. Thus, the Schauder fixed point theorem gives that there exists at least one fixed point $u\in \mathcal{S}_{T}$ of Φ.

(ii) Regularity and nonnegativity. By the elliptic regularity theory we see that $v(\cdot,t)\in C^{2+\theta }(\bar{\Omega})$ and $w(\cdot,t)\in C^{2+\theta}(\bar{\Omega})$. It then follows from (2.5) that $v(x,t)\in C^{2+\theta,\frac{\theta }{2}} (\bar{\Omega}\times[\eta,T] )$ and $w(x,t)\in C^{2+\theta ,\frac{\theta}{2}} (\bar{\Omega}\times[\eta,T] )$ for all $\eta\in (0,T)$. The parabolic regularity theory [38], Theorem 6.1, entails that

$$u(x,t)\in C^{2+\theta,1+\frac{\theta}{2}} \bigl(\bar{\Omega}\times[\eta ,T] \bigr) \quad \mbox{for all } \eta\in(0,T). $$

We may prolong the solution to the interval $[0,T_{\max})$ with either $T_{\max}=\infty$ or $T_{\max}<\infty$, where, in the latter case,

$$\bigl\| u(\cdot,t)\bigr\| _{L^{\infty}(\Omega)}\rightarrow\infty \quad\mbox{as }t\rightarrow T_{\max}. $$

Finally, the parabolic and elliptic comparison principles ensure the nonnegativity of u, v, and w.

(iii) Uniqueness. The proof for the uniqueness of solutions to system (1.1) is inspired by a method in [23]. We suppose that $(u_{1},v_{1},w_{1})$ and $(u_{2},v_{2},w_{2})$ are two classical solutions to system (1.2) in $\Omega\times(0,T)$ with the same initial data. Fix $T_{1}\in(0,T)$.

It is clear that $v_{1}-v_{2}$ satisfies the equation

$$ -\Delta(v_{1}-v_{2})+\beta(v_{1}-v_{2})= \alpha(u_{1}-u_{2}). $$

(2.7)

Thus, we differentiate (2.7) on t and then take $v_{1}-v_{2}$ as a test function to have

$$\begin{aligned} &\frac{1}{2}\frac{d}{dt} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx+ \frac{\beta }{2}\frac{d}{dt} \int_{\Omega}|v_{1}-v_{2}|^{2}\,dx \\ &\quad=\alpha \int_{\Omega}(u_{1}-u_{2})_{t}(v_{1}-v_{2})\,dx \\ &\quad=-\alpha \int_{\Omega}\nabla\bigl(A(u_{1})-A(u_{2}) \bigr)\cdot\nabla(v_{1}-v_{2})\,dx+\alpha \chi \int_{\Omega}(u_{1}\nabla v_{1}-u_{2} \nabla v_{2})\cdot\nabla(v_{1}-v_{2})\,dx \\ &\qquad{}-\alpha\xi \int_{\Omega}(u_{1}\nabla w_{1}-u_{2} \nabla w_{2})\cdot\nabla(v_{1}-v_{2})\,dx \end{aligned}$$

(2.8)

for any $t\in(0,T_{1})$, where $A(s)=\int_{0} ^{s}D(s)\,ds$. For the first term on the right-hand side of (2.8), we obtain from the mean value theorem and the Young inequality that

$$\begin{aligned} &-\alpha \int_{\Omega}\nabla\bigl(A(u_{1})-A(u_{2}) \bigr)\cdot\nabla(v_{1}-v_{2})\,dx \\ &\quad= \alpha \int_{\Omega}\bigl(A(u_{1})-A(u_{2})\bigr) \cdot\Delta(v_{1}-v_{2})\,dx \\ &\quad= -\alpha^{2} \int_{\Omega}\bigl(A(u_{1})-A(u_{2})\bigr) (u_{1}-u_{2})\,dx+\alpha\beta \int _{\Omega}\bigl(A(u_{1})-A(u_{2})\bigr) (v_{1}-v_{2})\,dx \\ &\quad= \alpha\beta C_{4} \int_{\Omega}(u_{1}-u_{2}) (v_{1}-v_{2})\,dx -\alpha^{2}C_{4} \int _{\Omega}(u_{1}-u_{2})^{2}\,dx \\ &\quad\leq \frac{C_{4}\beta^{2}}{2} \int_{\Omega}|v_{1}-v_{2}|^{2}\,dx - \frac {C_{4}\alpha^{2}}{2} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx \end{aligned}$$

(2.9)

for some positive constant $C_{4}\in [D(s_{1}),D(s_{2}) ]$, where

$$\begin{aligned}& s_{1}:=\min \bigl\{ \|u_{1}\|_{L^{\infty} (\Omega\times(0,T_{1}) )},\| u_{2}\|_{L^{\infty} (\Omega\times(0,T_{1}) )} \bigr\} \quad \mbox{and}\\& s_{2}:=\max\bigl\{ \|u_{1}\|_{L^{\infty} (\Omega\times(0,T_{1}) )},\|u_{2}\|_{L^{\infty} (\Omega\times(0,T_{1}) )} \bigr\} . \end{aligned}$$

For the second integral on the right-hand side of (2.8), we can use the Hölder’s inequality to have

$$ \begin{aligned}[b] &\alpha\chi \int_{\Omega}(u_{1}\nabla v_{1}-u_{2} \nabla v_{2})\cdot\nabla(v_{1}-v_{2})\,dx \\ &\quad\leq\alpha \chi \biggl( \int_{\Omega}|u_{1}\nabla v_{1}-u_{2} \nabla v_{2}|^{2}\,dx \biggr)^{\frac{1}{2}} \biggl( \int_{\Omega}\bigl| \nabla(v_{1}- v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}}. \end{aligned} $$

(2.10)

Notice that $| u_{2}|\leq C_{5}$, $|\nabla v_{1}|\leq C_{6}$, and $|\nabla w_{1}|\leq C_{7}$ with some positive constants $C_{6}$, $C_{5}$, and $C_{7}$ in $\Omega\times(0, T_{1})$. Thus,

$$\begin{aligned} &\int_{\Omega}|u_{1}\nabla v_{1}-u_{2} \nabla v_{2}|^{2}\,dx \\ &\quad\leq2 \int_{\Omega}|u_{1}-u_{2}|^{2}| \nabla v_{1}|^{2}\,dx+2 \int_{\Omega}u_{2}^{2}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \\ &\quad\leq2 C_{6}^{2} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx +2 C_{5}^{2} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx. \end{aligned}$$

(2.11)

Inserting (2.11) into (2.10) and using Young’s inequality, we obtain

$$\begin{aligned} &\alpha\chi \int_{\Omega}(u_{1}\nabla v_{1}-u_{2} \nabla v_{2})\cdot\nabla (v_{1}-v_{2})\,dx \\ &\quad\leq \alpha\chi\sqrt{2} \biggl(C_{6}^{2} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx +C_{5}^{2} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \biggl( \int_{\Omega}\bigl|\nabla( v_{1}- v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \\ &\quad\leq \alpha\chi\sqrt{2} \biggl(C_{6} \biggl( \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx \biggr)^{\frac{1}{2}} \\ &\qquad{}+C_{5} \biggl( \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \biggr) \biggl( \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \\ &\quad\leq \alpha\chi\sqrt{2}C_{6} \biggl( \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx \biggr)^{\frac{1}{2}} \biggl( \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \\ &\qquad{}+ \alpha\chi\sqrt{2}C_{5} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \\ &\quad\leq \frac{\alpha^{2}C_{4}}{8} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx+ \biggl(\frac {4\chi^{2}C_{6}^{2}}{C_{4}}+\alpha\chi\sqrt{2} C_{5} \biggr) \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx. \end{aligned}$$

(2.12)

Similarly, we can conclude that

$$\begin{aligned} &{-}\alpha\xi \int_{\Omega}(u_{1}\nabla w_{1}-u_{2} \nabla w_{2})\cdot\nabla (v_{1}-v_{2})\,dx \\ &\quad\leq \frac{\alpha^{2}C_{4}}{8} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx+ \frac{4\xi ^{2}C_{7}^{2}}{C_{4}} \int_{\Omega}\bigl| \nabla(v_{1}-v_{2})\bigr|^{2}\,dx \\ &\qquad{}+\alpha\xi\sqrt{2} C_{5} \biggl( \int_{\Omega}\bigl| \nabla (v_{1}-v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \biggl( \int_{\Omega}\bigl| \nabla (w_{1}-w_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}}. \end{aligned}$$

(2.13)

To estimate the last integral in (2.13), we notice that $w_{1}-w_{2}$ satisfies the equation

$$-\Delta(w_{1}-w_{2})+\delta(w_{1}-w_{2})= \gamma(u_{1}-u_{2}). $$

Taking $w_{1}-w_{2}$ as a test function, we obtain

$$ \int_{\Omega}\bigl|\nabla(w_{1}-w_{2})\bigr|^{2}\,dx= \gamma \int_{\Omega}(u_{1}-u_{2}) (w_{1}-w_{2})\,dx-\delta \int_{\Omega}(w_{1}-w_{2})^{2}\,dx, $$

(2.14)

which, together with Young’s inequality, yields that

$$\begin{aligned}[b] \int_{\Omega}\bigl|\nabla(w_{1}-w_{2})\bigr|^{2}\,dx &\leq\delta \int_{\Omega}|w_{1}-w_{2}|^{2}\,dx+ \frac{\gamma^{2}}{4\delta} \int _{\Omega}|u_{1}-u_{2}|^{2}\,dx- \delta \int_{\Omega}|w_{1}-w_{2}|^{2}\,dx \\ &= \frac{\gamma^{2}}{4\delta} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx. \end{aligned} $$

Thus, the last term in (2.13) can be estimated as

$$\begin{aligned} & \alpha\xi\sqrt{2} C_{5} \biggl( \int_{\Omega}\bigl| \nabla(v_{1}-v_{2})\bigr|^{2}\,dx \biggr)^{\frac{1}{2}} \biggl( \int_{\Omega}\bigl| \nabla(w_{1}-w_{2})\bigr|^{2}\,dx \biggr)^{\frac {1}{2}} \\ &\quad\leq \frac{\delta\alpha^{2}C_{4}}{2\gamma^{2}} \int_{\Omega}\bigl|\nabla (w_{1}-w_{2})\bigr|^{2}\,dx+ \frac{\xi^{2}\gamma^{2} C_{5}^{2}}{\delta C_{4}} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx \\ &\quad\leq \frac{\alpha^{2}C_{4}}{8} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx+ \frac{\xi ^{2}\gamma^{2}C_{5}^{2}}{\delta C_{4}} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx. \end{aligned}$$

(2.15)

Summarily, combining (2.8), (2.9), (2.12), and (2.15), we obtain

$$\begin{aligned} &\frac{1}{2}\frac{d}{dt} \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx+ \frac{\beta }{2}\frac{d}{dt} \int_{\Omega}|v_{1}-v_{2}|^{2}\,dx+ \frac{\alpha^{2} C_{4}}{8} \int_{\Omega}|u_{1}-u_{2}|^{2}\,dx \\ &\quad\leq \biggl(\frac{4\chi^{2}C_{6}^{2}}{C_{4}}+\frac{4\xi ^{2}C_{6}^{2}}{C_{4}}+\frac{\xi^{2}\gamma^{2}C_{5}^{2}}{\delta C_{4}}+ \alpha\chi \sqrt{2}C_{5} \biggr) \int_{\Omega}\bigl|\nabla(v_{1}-v_{2})\bigr|^{2}\,dx\\ &\qquad{}+ \frac{\beta ^{2}C_{4}}{2} \int_{\Omega}|v_{1}-v_{2}|^{2}\,dx. \end{aligned}$$

By Gronwall’s inequality we derive that $v_{1}=v_{2}$ and $u_{1}=u_{2}$ in $\Omega\times(0,T_{1})$. By (2.14) we also have $w_{1}=w_{2}$ in $\Omega\times(0,T_{1})$. Hence, $v_{1}=v_{2}$, $u_{1}=u_{2}$, and $w_{1}=w_{2}$ in $\Omega\times(0,T)$ due to the arbitrariness of $T_{1}\in(0,T)$. This implies the uniqueness of solutions. □

The following lemma deals with the mass identities.

Lemma 2.2

Let the assumptions in Lemma 2.1 hold. Then the classical solution $(u,v,w)$ of (1.2) fulfills

$$\begin{aligned}& \bigl\| u(\cdot,t)\bigr\| _{L^{1}(\Omega)}=\|u_{0}\|_{L^{1}(\Omega)} \quad \textit{for all } t\in(0,T_{\max}), \end{aligned}$$

(2.16)

$$\begin{aligned}& \bigl\| v(\cdot,t)\bigr\| _{L^{1}(\Omega)} = \frac{\alpha}{\beta}\|u_{0}\| _{L^{1}(\Omega)} \quad\textit{for all } t\in(0,T_{\max}), \end{aligned}$$

(2.17)

$$\begin{aligned}& \bigl\| w(\cdot,t)\bigr\| _{L^{1}(\Omega)} = \frac{\gamma}{\delta}\|u_{0}\| _{L^{1}(\Omega)} \quad\textit{for all } t\in(0,T_{\max}). \end{aligned}$$

(2.18)

Moreover, we have

$$ u(x,t)>0 \quad \textit{for all } x\in\Omega, t>0, $$

(2.19)

provided that $u_{0}>0$.

Proof

We integrate each equation of (1.2) with respect to $x \in\Omega$ and then obtain

$$\frac{d}{dt} \int_{\Omega}u \,dx \equiv0,\qquad \alpha \int_{\Omega}u \,dx =\beta \int_{\Omega}v \,dx,\quad \mbox{and}\quad \gamma \int_{\Omega}u \,dx =\delta \int_{\Omega}w \,dx $$

for all $t\in(0,T_{\max})$. It is clear that (2.16)-(2.18) hold. By the maximum principle, we obtain the positivity (2.19) of u. □

3 Global bounded classical solutions in the case of nondegenerate diffusion

In this section, we mainly investigate the existence of global bounded classical solutions to system (1.2) with nondegenerate diffusion. We first consider the case that the repulsion prevails over the attraction (i.e., $\xi\gamma-\chi\alpha>0$).

Lemma 3.1

Assume that $\xi\gamma-\chi\alpha>0$. Suppose that $u_{0}\in W^{1,\infty }(\Omega)$ is a nonnegative function and D satisfies (1.3), (1.4), and (1.5). Then, for any $p>\frac{n}{2}$, there exists a constant $C>0$ independent of t such that the solution $(u,v,w)$ of (1.2) fulfills

$$ \int_{\Omega}u^{p}(x,t)\,dx\leq C \quad\textit{for all } t \in (0,T_{\max}). $$

(3.1)

Proof

We multiply the first equation in (1.2) by $u^{p-1}$ and integrate by parts over Ω to have

$$\begin{aligned} &\frac{1}{p}\frac{d}{dt} \int_{\Omega}u^{p}\,dx \\ &\quad = \int_{\Omega}u^{p-1}\nabla \cdot\bigl(D(u)\nabla u \bigr)\,dx - \int_{\Omega}u^{p-1} \nabla\cdot(\chi u \nabla v)\,dx+ \int_{\Omega}u^{p-1} \nabla\cdot(\xi u \nabla w)\,dx \\ &\quad=-(p-1) \int_{\Omega}u^{p-2}D(u)|\nabla u|^{2}\,dx + \frac{(p-1)\chi}{p} \int _{\Omega}\nabla u^{p}\cdot\nabla v \,dx\\ &\qquad{} - \frac{(p-1)\xi}{p} \int_{\Omega}\nabla u^{p}\cdot\nabla w \,dx \\ &\quad=-(p-1) \int_{\Omega}u^{p-2}D(u)|\nabla u|^{2}\,dx + \frac{p-1}{p} \int _{\Omega}u^{p}(-\chi\Delta v +\xi\Delta w)\,dx \end{aligned}$$

for all $t\in(0,T_{\max})$. Thus, from the second and third equations in (1.2) we obtain

$$\begin{aligned} \frac{d}{dt} \int_{\Omega}u^{p}\,dx ={}&{-}p(p-1) \int_{\Omega}u^{p-2}D(u)|\nabla u|^{2}\,dx \\ &{}+ (p-1) \int_{\Omega}u^{p} \bigl[\xi\delta w-(\xi\gamma-\chi \alpha)u-\chi\beta v \bigr]\,dx, \end{aligned}$$

(3.2)

which, together with $v\geq0$, yields that

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq-(\xi\gamma-\chi\alpha) (p-1) \int _{\Omega}u^{p+1}\,dx +\xi\delta(p-1) \int_{\Omega}u^{p}w \,dx. $$

(3.3)

By $\xi\gamma-\chi\alpha>0$ and Young’s inequality we deduce that

$$ \xi\delta(p-1) \int_{\Omega}u^{p}w \,dx \leq\frac{\xi\gamma-\chi\alpha }{2}(p-1) \int_{\Omega}u^{p+1}\,dx + C_{1} \int_{\Omega}w^{p+1}\,dx, $$

(3.4)

where $C_{1}:=\xi\delta\frac{p-1}{p+1} [\frac{2\xi\delta p}{(\xi\gamma -\chi\alpha)(p+1)} ]^{p}$. Substituting (3.4) into (3.3) yields

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq-\frac{\xi\gamma-\chi\alpha }{2}(p-1) \int_{\Omega}u^{p+1}\,dx + C_{1} \int_{\Omega}w^{p+1}\,dx. $$

(3.5)

Following a similar procedure as in [33], we go to estimate the term $\int_{\Omega}w^{p+1}\,dx$. Here we give a sketch for completeness. Since w solves

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} -\Delta w+\delta w=\gamma u,& x\in\Omega,\\ \frac{\partial w}{\partial\nu}=0,&x\in\partial\Omega, \end{array}\displaystyle \right . $$

(3.6)

where $\delta>0$ and $\gamma>0$, we can apply $L^{p}$ estimates [39, 40] on (3.6) to obtain

$$ \bigl\| w(\cdot,t)\bigr\| _{W^{2,p}(\Omega)}\leq C_{2} \bigl\| u(\cdot,t) \bigr\| _{L^{p}(\Omega)}\quad \quad \mbox{for all } t\in(0,T_{\max}) $$

(3.7)

with some constant $C_{2}>0$. Then by the Gagliardo-Nirenberg interpolation inequality [41] and the $L^{1}$ estimates of w (Lemma 2.2) we find that

$$\begin{aligned} \int_{\Omega}w^{p+1}\,dx &\leq C_{3} \bigl\| D^{2}w\bigr\| _{L^{p}(\Omega)}^{(p+1)\theta }\|w\|_{L^{1}(\Omega)}^{(p+1)(1-\theta)}+C_{3} \|w\|_{L^{1}(\Omega )}^{p+1} \\ & \leq C_{4} \|u\|_{L^{p}(\Omega)}^{(p+1)\theta}+C_{4} \quad\mbox{for all } t\in(0,T_{\max}) \end{aligned}$$

(3.8)

with some constants $C_{3}>0$ and $C_{4}>0$, where

$$\theta=\frac{1-\frac{1}{p+1}}{1+\frac{2}{n}-\frac{1}{p}}. $$

Since $p>\frac{n}{2}$, it is easy to check that $\theta\in(0,1)$ and $(p+1)\theta< p$. Hence, using Young’s inequality twice, we have

$$\begin{aligned} \int_{\Omega}w^{p+1}\,dx &\leq C_{4} \bigl(\|u \|_{L^{p}(\Omega)}^{p}+1 \bigr)+C_{4} \\ & \leq C_{4} \biggl(\kappa \int_{\Omega}u^{p+1}\,dx+\frac{|\Omega|}{p+1} \biggl[ \frac{p}{\kappa(p+1)} \biggr]^{p} \biggr)+2C_{4} \\ &= C_{4} \kappa \int_{\Omega}u^{p+1}\,dx +C_{5}(\kappa) \quad \mbox{for all } t\in(0,T_{\max}). \end{aligned}$$

(3.9)

Substituting (3.9) into (3.5), we obtain

$$\frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq-\frac{\xi\gamma-\chi\alpha }{2}(p-1) \int_{\Omega}u^{p+1}\,dx + C_{1}C_{4} \kappa \int_{\Omega}u^{p+1}\,dx+C_{1}C_{5}( \kappa). $$

Then by taking

$$\kappa=\frac{(\xi\gamma-\chi\alpha)(p-1)}{4C_{1} C_{4}} $$

we have

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq -\frac{\xi\gamma-\chi\alpha }{4}(p-1) \int_{\Omega}u^{p+1}\,dx + C_{6} \quad\mbox{for all } t\in(0,T_{\max}), $$

(3.10)

where $C_{6}:=C_{1} C_{5}$. By Young’s inequality again, we obtain

$$ \int_{\Omega}u^{p}\,dx \leq\frac{\xi\gamma-\chi\alpha}{4}(p-1) \int_{\Omega}u^{p+1}\,dx +\frac{|\Omega|}{p+1} \biggl[ \frac{4p}{(\xi\gamma-\chi\alpha )(p^{2}-1)} \biggr]^{p} $$

(3.11)

for $t\in(0,T_{\max})$. Thus, we conclude that

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx+ \int_{\Omega}u^{p}\,dx \leq C_{7} \quad \mbox{for all } t\in(0,T_{\max}), $$

(3.12)

where $C_{7}:=C_{6}+\frac{|\Omega|}{p+1} [\frac{4p}{(\xi\gamma-\chi \alpha)(p^{2}-1)} ]^{p}$. By Gronwall’s inequality we have

$$\int_{\Omega}u^{p}(x,t)\,dx\leq\max \biggl\{ \int_{\Omega}u_{0}^{p}\,dx, C_{7} \biggr\} \quad\mbox{for all } t\in(0,T_{\max}), $$

which implies the desired uniform estimates. □

Next, considering the case that the attraction dominates over the repulsion, we can deduce a similar uniform estimate under the assumption of $m>2-\frac{2}{n}$. We will show that the stronger diffusion plays a key role in deducing such a uniform bound.

Lemma 3.2

Assume that $\xi\gamma-\chi\alpha\le0$ and $m>2-\frac{2}{n}$. Suppose that $u_{0}\in W^{1,\infty}(\Omega)$ is a nonnegative function and D satisfies (1.3), (1.4), and (1.5). Then, for any $p>\frac{n}{2}$, there exists a constant $C>0$ independent of t such that the solution $(u,v,w)$ of system (1.2) fulfills

$$ \int_{\Omega}u^{p}(x,t)\,dx\leq C \quad\textit{for all } t\in (0,T_{\max}). $$

(3.13)

Proof

Combining (3.2) with (1.4), we derive

$$\begin{aligned} \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq{}&{-}\frac{4 c_{D}(p-1)p}{(p+m-1)^{2}} \int_{\Omega}\bigl|\nabla u^{\frac {p+m-1}{2}}\bigr|^{2}\,dx \\ &{}+ (p-1) \int_{\Omega}u^{p} \bigl[\xi\delta w + (\chi\alpha-\xi \gamma)u-\chi\beta v \bigr]\,dx \end{aligned}$$

for all $t\in(0,T_{\max})$. This, along with $v\geq0$, yields

$$\begin{aligned} \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq{}&{-}\frac{4 c_{D}(p-1)p}{(p+m-1)^{2}} \int_{\Omega}\bigl|\nabla u^{\frac{p+m-1}{2}}\bigr|^{2}\,dx \\ &{}+ (p-1)\xi \delta \int_{\Omega}u^{p} w \,dx+ (p-1) (\chi\alpha-\xi\gamma) \int_{\Omega}u^{p+1}\,dx. \end{aligned}$$

By Young’s inequality we obtain

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq-\frac{4 c_{D}(p-1)p}{(p+m-1)^{2}} \int_{\Omega}\bigl|\nabla u^{\frac {p+m-1}{2}}\bigr|^{2}\,dx +C_{1} \int_{\Omega}u^{p+1}\,dx +C_{2} \int_{\Omega}w^{p+1}\,dx, $$

(3.14)

where $C_{1}:=(p-1)(\chi\alpha-\xi\gamma+\xi\delta)$ and $C_{2}:=(p-1)\xi \delta$. Similarly to the deduction of (3.8) in Lemma 3.1, we find that there exist some constants $C_{3}>0$ and $C_{4}>0$ such that

$$\begin{aligned} \int_{\Omega}w^{p+1}\,dx &\leq C_{3} \bigl\| D^{2}w\bigr\| _{L^{p}(\Omega)}^{(p+1)\theta}\|w\|_{L^{1}(\Omega )}^{(p+1)(1-\theta)}+C_{3} \|w\|_{L^{1}(\Omega)}^{p+1} \\ &\leq C_{4} \|u\|_{L^{p}(\Omega)}^{(p+1)\theta}+C_{4} \quad\mbox{for all } t\in(0,T_{\max}), \end{aligned}$$

where $(p+1)\theta=\frac{np^{2}}{np+2p-n}\in(0,p)$ by $p>\frac{n}{2}$. Using Young’s inequality twice yields

$$\begin{aligned} \int_{\Omega}w^{p+1}\,dx &\leq C_{4} \bigl(\|u \|_{L^{p}(\Omega)}^{p}+1 \bigr)+C_{4} \\ &\leq C_{4} \biggl( \int_{\Omega}u^{p+1}\,dx +|\Omega| \biggr)+2C_{4} \\ &=C_{4} \int_{\Omega}u^{p+1}\,dx+C_{5} \quad\mbox{for all } t\in (0,T_{\max}), \end{aligned}$$

(3.15)

where $C_{5}:=C_{4}(|\Omega|+2)$. Hence, inserting (3.15) into (3.14), we obtain

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq-\frac{4 c_{D}(p-1)p}{(p+m-1)^{2}} \int_{\Omega}\bigl|\nabla u^{\frac {p+m-1}{2}}\bigr|^{2}\,dx +C_{6} \int_{\Omega}u^{p+1}\,dx +C_{7}, $$

(3.16)

where $C_{6}:=C_{1}+C_{2}C_{4}$ and $C_{7}:=C_{2}C_{5}$. Since $p>\frac{n}{2}$ and $m\geq1$, we have $p>\frac{(2-m)n}{2}>\frac{2n-mn-2}{2}$ and then obtain

$$\frac{2}{p+m-1}< \frac{2(p+1)}{p+m-1}< \frac{2n}{n-2}. $$

By the Gagliardo-Nirenberg inequality we derive that there exists $C_{8}>0$ such that

$$\begin{aligned} \int_{\Omega} u ^{p+1}\,dx &= \bigl\| u^{\frac{p+m-1}{2}} \bigr\| ^{\frac{2(p+1)}{p+m-1}}_{L^{\frac {2(p+1)}{p+m-1}}(\Omega)} \\ &\leq C_{8} \bigl(\bigl\| \nabla u^{\frac{p+m-1}{2}}\bigr\| ^{\theta_{1}}_{L^{2} (\Omega)} \bigl\| u^{\frac{p+m-1}{2}}\bigr\| ^{1-\theta_{1}}_{L^{\frac{2}{p+m-1}}(\Omega )} \bigr)^{\frac{2(p+1)}{p+m-1}}+C_{8} \bigl\| u^{\frac{p+m-1}{2}}\bigr\| _{L^{\frac {2}{p+m-1}}(\Omega)}^{\frac{2(p+1)}{p+m-1}} \\ &\leq C_{9} \bigl\| \nabla u^{\frac{p+m-1}{2}}\bigr\| ^{\frac{2\theta _{1}(p+1)}{p+m-1}}_{L^{2} (\Omega)} + C_{9} \quad\mbox{for all } t\in(0,T_{\max}), \end{aligned}$$

(3.17)

where

$$\theta_{1}=\frac{\frac{p+m-1}{2}-\frac{p+m-1}{2(p+1)}}{\frac{1}{n}-\frac {1}{2}+\frac{p+m-1}{2}}\in(0,1) $$

and

$$C_{9}:=\max \bigl\{ C_{8} \|u\|_{L^{1}(\Omega)}^{(1-\theta_{1})(p+1)}, C_{8} \|u\| _{L^{1}(\Omega)}^{p+1} \bigr\} =\max \bigl\{ C_{8}\|u_{0}\|_{L^{1}(\Omega )}^{(1-\theta_{1})(p+1)}, C_{8} \|u_{0}\|_{L^{1}(\Omega)}^{p+1} \bigr\} . $$

Since $m>2-\frac{2}{n}$, we have

$$\frac{2 \theta_{1} (p+1)}{p+m-1}=\frac{p}{\frac{1}{n}-\frac{1}{2}+\frac{p+m-1}{2}} < \frac{p}{\frac{1}{n}-\frac{1}{2}+\frac{p+2-\frac{2}{n}-1}{2}}=2. $$

Thus, we use Young’s inequality to derive

$$\begin{aligned} C_{6} \int_{\Omega}u^{p+1}\,dx & \leq C_{6}C_{9} \bigl\| \nabla u^{\frac{p+m-1}{2}}\bigr\| ^{\frac{2\theta _{1}(p+1)}{p+m-1}}_{L^{2} (\Omega)} + C_{6}C_{9} \\ & \leq\frac{2 c_{D}(p-1)p}{(p+m-1)^{2}} \int_{\Omega}\bigl|\nabla u^{\frac {p+m-1}{2}}\bigr|^{2}\,dx + C_{10}, \end{aligned}$$

(3.18)

where

$$C_{10}:=C_{6}C_{9}+(C_{6}C_{9})^{\frac{p+m-1}{p+m-1-\theta_{1}(p+1)}} \biggl( \frac {2 c_{D}(p-1)p}{(p+m-1)^{2}} \biggr)^{-\frac{\theta_{1}(p+1)}{p+m-1-\theta_{1}(p+1)}}. $$

Inserting (3.18) into (3.16) yields

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx \leq-\frac{2 c_{D}(p-1)p}{(p+m-1)^{2}} \int_{\Omega}\bigl|\nabla u^{\frac{p+m-1}{2}}\bigr|^{2}\,dx +C_{11} \quad\mbox{for all } t\in(0,T_{\max}), $$

(3.19)

where $C_{11}:=C_{7}+C_{10}$. Since $p>\frac{n}{2}$ and $m\geq1$, it is easy to check that

$$\frac{2}{p+m-1}< \frac{2p}{p+m-1}< \frac{2n}{n-2}. $$

By using the Gagliardo-Nirenberg inequality again, we can find a constant $C_{12}>0$ such that

$$\begin{aligned} \int_{\Omega} u ^{p}\,dx &=\bigl\| u^{\frac{p+m-1}{2}} \bigr\| ^{\frac{2p}{p+m-1}}_{L^{\frac {2p}{p+m-1}}(\Omega)} \\ &\leq C_{12} \bigl\| \nabla u^{\frac{p+m-1}{2}}\bigr\| ^{\frac{2\theta_{2} p }{p+m-1}}_{L^{2} (\Omega)} \bigl\| u^{\frac{p+m-1}{2}}\bigr\| ^{\frac{2p(1-\theta_{2})}{p+m-1}}_{L^{\frac {2}{p+m-1}} (\Omega)}+ C_{12} \bigl\| u^{\frac{p+m-1}{2}}\bigr\| ^{\frac{2p}{p+m-1}}_{L^{\frac{2}{p+m-1}} (\Omega )} \\ &\leq C_{13}\bigl\| \nabla u^{\frac{p+m-1}{2}}\bigr\| ^{\frac{2\theta_{2} p }{p+m-1}}_{L^{2} (\Omega)}+C_{13} \quad\mbox{for all } t\in(0,T_{\max}), \end{aligned}$$

(3.20)

where

$$\theta_{2}=\frac{p+m-1}{2}\frac{1-\frac{1}{p}}{\frac{1}{n}-\frac {1}{2}+\frac{p+m-1}{2}}\in(0,1) $$

and

$$C_{13}:=\max \bigl\{ C_{12} \|u_{0} \|_{L^{1}(\Omega)}^{(1-\theta_{2})p},C_{12} \|u_{0} \|_{L^{1}(\Omega)}^{p} \bigr\} . $$

Since $m>2-\frac{2}{n}$, we find that $\frac{2 \theta_{2} p}{p+m-1}=\frac{p-1}{\frac{1}{n}-\frac{1}{2}+\frac{p+m-1}{2}} <\frac{p-1}{\frac{1}{n}-\frac{1}{2}+\frac{p+2-\frac{2}{n}-1}{2}}=\frac {2(p-1)}{p}<2$. Thus, using Young’s inequality yields that

$$ \|u\|_{L^{p}(\Omega)}^{p}\leq\frac{2 c_{D}(p-1)p}{(p+m-1)^{2}} \int _{\Omega}\bigl|\nabla u^{\frac{p+m-1}{2}}\bigr|^{2}\,dx + C_{14}, $$

(3.21)

where

$$C_{14}:= \biggl[\frac{(p+m-1)^{2}}{2c_{D}(p-1)p} \biggr]^{\frac{\theta_{2} p}{p+m-1-\theta_{2} p}}C_{13}^{\frac{p+m-1}{p+m-1-\theta_{2} p}}+C_{13}. $$

Substituting (3.21) into (3.19) yields that

$$ \frac{d}{dt} \int_{\Omega}u^{p}\,dx+ \int_{\Omega}u^{p}\,dx \leq C_{15}\quad \mbox{for all } t\in(0,T_{\max}), $$

(3.22)

where $C_{15}:=C_{11}+C_{14}$. Thus, using Gronwall’s inequality, we have

$$\int_{\Omega}u^{p}\,dx \leq\max \biggl\{ \int_{\Omega}u_{0}^{p}\,dx,C_{15} \biggr\} \quad \mbox{for all } t\in(0,T_{\max}), $$

which implies the desired uniform $L^{p}$ estimates. □

We now turn to the existence of global bounded classical solutions.

Proof of Theorem 1.1

According to the $L^{p}$ estimates of w (see (3.7)), we obtain from Lemmas 3.1 and 3.2 that

$$\sup_{0< t< T_{\max}}\bigl\| w(\cdot,t)\bigr\| _{W^{2,p}(\Omega)}\leq C_{1} \quad\mbox{for all } t\in(0,T_{\max}) $$

with some positive constant $C_{1}$. Then, by choosing $p>n$, from the Sobolev embedding theorem we can derive that there exists a constant $C_{2}>0$ such that

$$ \sup_{0< t< T_{\max}}\bigl\| \nabla w(\cdot,t)\bigr\| _{L^{\infty}(\Omega)} \leq C_{2} \quad\mbox{for all } t\in(0,T_{\max}). $$

(3.23)

Similarly, there exists a constant $C_{3}>0$ such that

$$ \sup_{0< t< T_{\max}}\bigl\| \nabla v(\cdot,t)\bigr\| _{L^{\infty}(\Omega)} \leq C_{3} \quad \mbox{for all } t\in(0,T_{\max}). $$

(3.24)

With the aid of Lemmas 3.1 and 3.2 and using Lemma A.1 in [20] (see also [42]), we can conclude that there exists a positive constant $C_{4}>0$ such that

$$ \bigl\| u(\cdot,t)\bigr\| _{L^{\infty}(\Omega)}\leq C_{4} \quad\mbox{for all } t\in(0,T_{\max}), $$

(3.25)

which, together with the extensibility criterion (2.1), implies that $T_{\max}=+\infty$. Thus, $(u,v,w)$ is a global bounded classical solution to system (1.2). □

4 Global bounded weak solutions in the case of degenerate diffusion

In this section, we consider system (1.2) with degenerate diffusion (i.e., $D(u)\geq0$ for all $u\geq0$). We first consider the following regularized system with nondegenerate diffusion for $\varepsilon\in(0,1)$, which satisfies all the formal arguments:

$$ \left \{ \textstyle\begin{array}{@{}l@{\quad}l} u_{\varepsilon t}=\nabla\cdot(D_{\varepsilon}(u_{\varepsilon})\nabla u_{\varepsilon})-\nabla\cdot(\chi u_{\varepsilon}\nabla v_{\varepsilon})+\nabla\cdot(\xi u_{\varepsilon}\nabla w_{\varepsilon}),& x\in\Omega, t>0,\\ 0=\Delta v_{\varepsilon}+\alpha u_{\varepsilon}-\beta v_{\varepsilon},& x\in \Omega, t>0,\\ 0=\Delta w_{\varepsilon}+\gamma u_{\varepsilon}-\delta w_{\varepsilon},& x\in \Omega, t>0,\\ \frac{\partial u_{\varepsilon}}{\partial\nu}=\frac{\partial v_{\varepsilon}}{\partial\nu}=\frac{\partial w_{\varepsilon}}{\partial\nu }=0,& x\in\partial\Omega, t>0,\\ u_{\varepsilon}(x,0)=u_{0}(x), &x\in\Omega, \end{array}\displaystyle \right . $$

(4.1)

where $D_{\varepsilon}$ is defined by

$$D_{\varepsilon}(s):=D(s+\varepsilon) \quad\mbox{for all }s\geq0. $$

Thus, $D_{\varepsilon}$ satisfies (1.3), (1.4), and (1.5). The following proposition is a direct consequence of Theorem 1.1.

Proposition 4.1

Let $\varepsilon\in(0,1)$, and let $u_{0}\in W^{1,\infty}(\Omega)$ be a nonnegative function. Suppose that $\xi\gamma-\chi\alpha>0$ or $\xi\gamma-\chi\alpha\le0$ and $m>2-\frac{2}{n}$. Then system (4.1) admits a unique global bounded classical solution $(u_{\varepsilon}, v_{\varepsilon}, w_{\varepsilon})$.

Next, we go to find some estimates to $(u_{\varepsilon}, v_{\varepsilon}, w_{\varepsilon})$, which are independent of ε and used to obtain some convergence properties. By taking $\varepsilon\to0$ we will establish the existence of global bounded weak solutions. The following two lemmas based on the ideas in [18] are used to prove the existence of the limit function of $\nabla\int^{u_{\varepsilon}+\varepsilon}_{0} D(z)\,dz$.

Lemma 4.1

Let $T>0$, and let the assumptions in Proposition 4.1 hold. Let $(u_{\varepsilon},v_{\varepsilon},w_{\varepsilon})$ be a solution to system (4.1) on $(0,T)$. Then

$$ \bigl\| D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon)\nabla u_{\varepsilon}\bigr\| _{L^{2} (0,T;L^{2}(\Omega) )}^{2}\leq\frac {1}{2} \|u_{0}\|^{2}_{L^{2}(\Omega)}+C_{1} T, $$

(4.2)

where $C_{1}$ is a positive constant independent of ε.

Proof

Taking $u_{\varepsilon}$ as a test function on the first equation in (4.1) and integrating it over $\Omega\times(0,T)$, we derive

$$\begin{aligned} &\frac{1}{2} \bigl(\bigl\| u_{\varepsilon}(T)\bigr\| _{L^{2}(\Omega )}^{2}- \|u_{0}\|_{L^{2}(\Omega)}^{2} \bigr) \\ &\quad= - \int_{0}^{T} \int_{\Omega} D(u_{\varepsilon}+\varepsilon)|\nabla u_{\varepsilon}|^{2}\,dx \,dt + \int_{0}^{T} \int_{\Omega} \chi u_{\varepsilon }\nabla v_{\varepsilon}\cdot \nabla u_{\varepsilon}\,dx \,dt\\ &\qquad{} - \int _{0}^{T} \int_{\Omega} \xi u_{\varepsilon}\nabla w_{\varepsilon}\cdot \nabla u_{\varepsilon}\,dx \,dt \\ &\quad= - \int_{0}^{T} \int_{\Omega} D(u_{\varepsilon}+\varepsilon)|\nabla u_{\varepsilon}|^{2}\,dx \,dt +\frac{\chi}{2} \int_{0}^{T} \int_{\Omega} \nabla v_{\varepsilon}\cdot\nabla u_{\varepsilon}^{2}\,dx \,dt \\ &\qquad{}-\frac{\xi }{2} \int_{0}^{T} \int_{\Omega} \nabla w_{\varepsilon}\cdot\nabla u_{\varepsilon}^{2}\,dx \,dt \\ &\quad= - \int_{0}^{T} \int_{\Omega} D(u_{\varepsilon}+\varepsilon)|\nabla u_{\varepsilon}|^{2}\,dx \,dt -\frac{\chi}{2} \int_{0}^{T} \int_{\Omega} u_{\varepsilon}^{2} \Delta v_{\varepsilon}\,dx \,dt +\frac{\xi}{2} \int _{0}^{T} \int_{\Omega} u_{\varepsilon}^{2} \Delta w_{\varepsilon}\,dx \,dt. \end{aligned}$$

It then follows from the second and third equations in (4.1) that

$$\begin{aligned} &\frac{1}{2} \bigl(\bigl\| u_{\varepsilon}(T)\bigr\| _{L^{2}(\Omega )}^{2}- \|u_{0}\|_{L^{2}(\Omega)}^{2} \bigr) \\ &\quad= - \int_{0}^{T} \int_{\Omega} D(u_{\varepsilon}+\varepsilon)|\nabla u_{\varepsilon}|^{2}\,dx \,dt -\frac{\chi}{2} \int_{0}^{T} \int_{\Omega} (\beta v_{\varepsilon}-\alpha u_{\varepsilon}) u_{\varepsilon}^{2}\,dx \,dt \\ &\qquad{}+\frac{\xi}{2} \int_{0}^{T} \int_{\Omega} (\delta w_{\varepsilon}-\gamma u_{\varepsilon}) u_{\varepsilon}^{2}\,dx \,dt. \end{aligned}$$

From Proposition 4.1 we obtain that there exist some positive constants $c_{1}$, $c_{2}$, $c_{3}$, $c_{4}$, and $c_{5}$ independent of ε such that

$$ \begin{aligned} &\|u_{\varepsilon}\|_{L^{\infty}(\Omega)}< c_{1},\qquad \|v_{\varepsilon}\|_{L^{\infty}(\Omega)}< c_{2},\qquad \|w_{\varepsilon}\|_{L^{\infty}(\Omega)}< c_{3},\\ & \|\nabla v_{\varepsilon}\|_{L^{\infty}(\Omega)}< c_{4},\quad \mbox{and}\quad \|\nabla w_{\varepsilon}\|_{L^{\infty}(\Omega)}< c_{5}. \end{aligned} $$

(4.3)

Therefore, we have

$$\begin{aligned} &\frac{1}{2} \bigl(\bigl\| u_{\varepsilon}(T)\bigr\| _{L^{2}(\Omega )}^{2}- \|u_{0}\|_{L^{2}(\Omega)}^{2} \bigr) \\ &\quad\leq - \int_{0}^{T} \int_{\Omega} D(u_{\varepsilon}+\varepsilon )|\nabla u_{\varepsilon}|^{2}\,dx \,dt+\frac{1}{2}c_{1}^{2} \bigl(\chi\beta c_{2} +(\xi\gamma+\chi\alpha)c_{1}+\xi\delta c_{3} \bigr)|\Omega|T, \end{aligned}$$

which yields the desired estimate

$$\bigl\| D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon)\nabla u_{\varepsilon}\bigr\| _{L^{2} (0,T;L^{2}(\Omega) )}\leq\frac {1}{2}\|u_{0}\|^{2}_{L^{2}(\Omega)}+C_{1} T, $$

where $C_{1}:=\frac{1}{2}c_{1}^{2} (\chi\beta c_{2}+(\xi\gamma+\chi\alpha )c_{1}+\xi\delta c_{3} )|\Omega|$. □

Lemma 4.2

Let $T>0$, and let the assumptions in Proposition 4.1 hold. Let $(u_{\varepsilon},v_{\varepsilon},w_{\varepsilon})$ be a solution to system (4.1) on $(0,T)$. Then

$$\begin{aligned} &\biggl\Vert \sqrt{t}\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz\biggr\Vert _{L^{2} (0,T;L^{2}(\Omega) )}^{2} +\sup_{t\in(0,T)} t\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon }D(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} \\ &\quad\leq C_{2}+C_{2} T +C_{2} T^{2}, \end{aligned}$$

(4.4)

where $C_{2}$ is a positive constant independent of ε.

Proof

We multiply the first equation in (4.1) by $\frac{d}{dt}\int _{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz$ and then integrate it over Ω to obtain

$$\begin{aligned} & \int_{\Omega} \biggl(\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon }D^{\frac{1}{2}}(z)\,dz \biggr)^{2}\,dx \\ & \quad= \int_{\Omega}\nabla\cdot\bigl(D(u_{\varepsilon}+\varepsilon) \nabla u_{\varepsilon}\bigr)\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz \,dx - \int_{\Omega}\nabla\cdot(\chi u_{\varepsilon}\nabla v_{\varepsilon})\frac {d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz \,dx \\ &\qquad{}+ \int_{\Omega}\nabla\cdot(\xi u_{\varepsilon}\nabla w_{\varepsilon})\frac {d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz \,dx \\ &\quad= - \int_{\Omega}D(u_{\varepsilon}+\varepsilon)\nabla u_{\varepsilon}\cdot\frac{d}{dt} \bigl(D(u_{\varepsilon}+\varepsilon) \nabla u_{\varepsilon}\bigr)\,dx \\ &\qquad{}-\chi \int_{\Omega} (\nabla u_{\varepsilon}\cdot\nabla v_{\varepsilon}+ u_{\varepsilon}\triangle v_{\varepsilon}) D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon) \biggl(\frac{d}{dt} \int _{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz \biggr)\,dx \\ &\qquad{} + \xi \int_{\Omega} (\nabla u_{\varepsilon}\cdot\nabla w_{\varepsilon}+u_{\varepsilon}\triangle w_{\varepsilon}) D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon) \biggl(\frac{d}{dt} \int _{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz \biggr)\,dx, \end{aligned}$$

where we used $\frac{d}{dt}\int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz =D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon)\frac{d}{dt}\int _{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz$. By Young’s inequality we obtain

$$\begin{aligned} & \int_{\Omega} \biggl(\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon }D^{\frac{1}{2}}(z)\,dz \biggr)^{2}\,dx \\ &\quad\leq -\frac{1}{2}\frac{d}{dt}\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} +\frac{1}{4} \int_{\Omega} \biggl(\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz \biggr)^{2}\,dx \\ &\qquad{}+\frac{1}{4} \int_{\Omega} \biggl(\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz \biggr)^{2}\,dx \\ &\qquad{} +\chi^{2} \int_{\Omega} (\nabla u_{\varepsilon}\cdot\nabla v_{\varepsilon}+ u_{\varepsilon}\Delta v_{\varepsilon})^{2} D(u_{\varepsilon}+\varepsilon)\,dx \\ &\qquad{} +\xi^{2} \int_{\Omega} (\nabla u_{\varepsilon}\cdot\nabla w_{\varepsilon}+ u_{\varepsilon}\Delta w_{\varepsilon})^{2} D(u_{\varepsilon}+\varepsilon)\,dx. \end{aligned}$$

(4.5)

Since $\|u_{\varepsilon}\|_{L^{\infty}(\Omega)}< c_{1}$ and $D\in C^{2} ([0,\infty) )$, we have $\|D(u_{\varepsilon}+\varepsilon)\|_{L^{\infty}(\Omega)}< c_{\infty}$ with some constant $c_{\infty}>0$. Thus, by Young’s inequality and (4.3) we derive

$$\begin{aligned} &\int_{\Omega} (\nabla u_{\varepsilon}\cdot\nabla v_{\varepsilon}+ u_{\varepsilon}\triangle v_{\varepsilon})^{2} D(u_{\varepsilon}+\varepsilon )\,dx \\ &\quad\leq2 \int_{\Omega} \bigl(|\nabla u_{\varepsilon}|^{2} | \nabla v_{\varepsilon}|^{2} +u_{\varepsilon}^{2} | \triangle v_{\varepsilon}|^{2} \bigr)D(u_{\varepsilon}+ \varepsilon)\,dx \\ &\quad \leq2 \| \nabla v_{\varepsilon}\|_{L^{\infty}(\Omega)}^{2} \biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} +2c_{1} ^{2} c_{\infty}\int_{\Omega}|\triangle v_{\varepsilon}|^{2}\,dx \\ &\quad\leq2 \| \nabla v_{\varepsilon}\|_{L^{\infty}(\Omega)}^{2} \biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2}+4c_{1} ^{2} c_{\infty}\int_{\Omega}\bigl(\alpha^{2} u_{\varepsilon}^{2}+ \beta^{2} v_{\varepsilon}^{2} \bigr)\,dx \\ &\quad\leq2c_{4} ^{2}\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2}+4c_{1} ^{2} c_{\infty}\bigl(\alpha^{2} c_{1} ^{2} +\beta^{2} c_{2} ^{2} \bigr)|\Omega|. \end{aligned}$$

Similarly,

$$\begin{aligned} &\int_{\Omega} (\nabla u_{\varepsilon}\cdot\nabla w_{\varepsilon}+ u_{\varepsilon}\triangle w_{\varepsilon})^{2} D(u_{\varepsilon}+\varepsilon)\,dx \\ &\quad\leq2c_{5} ^{2}\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2}+4c_{1} ^{2} c_{\infty}\bigl(\gamma^{2} c_{1} ^{2} +\delta^{2} c_{3} ^{2} \bigr)|\Omega|. \end{aligned}$$

Substituting the last two inequalities into (4.5), we obtain

$$\begin{aligned} & \int_{\Omega} \biggl(\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon }D^{\frac{1}{2}}(z)\,dz \biggr)^{2}\,dx +\frac{d}{dt}\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon }D(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} \\ &\quad\leq 4 \bigl(\chi^{2} c_{4} ^{2} + \xi^{2} c_{5} ^{2} \bigr)\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2}\\ &\qquad{}+8c_{1} ^{2} c_{\infty}\chi^{2} \bigl(\alpha^{2} c_{1} ^{2} + \beta ^{2} c_{2} ^{2} \bigr)|\Omega|+8c_{1} ^{2} c_{\infty}\xi^{2} \bigl(\gamma^{2} c_{1} ^{2} +\delta^{2} c_{3} ^{2} \bigr)|\Omega|. \end{aligned}$$

Setting $C_{\max}:=4\max \{ (\chi^{2} c_{4} ^{2} + \xi^{2} c_{5} ^{2} ), [2c_{1} ^{2} c_{\infty}\chi^{2}(\alpha^{2} c_{1} ^{2} +\beta^{2} c_{2} ^{2})|\Omega |+2c_{1} ^{2} c_{\infty}\xi^{2}(\gamma^{2} c_{1} ^{2} +\delta c_{3} ^{2})|\Omega| ] \}$ yields

$$\begin{aligned} &\biggl\Vert \frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} +\frac{d}{dt}\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon }D(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} \\ &\quad\leq C_{\max} \biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} +C_{\max}. \end{aligned}$$

(4.6)

Multiplying (4.6) by t and integrating it over $(0, T)$, we obtain

$$ \begin{aligned}[b] &\biggl\Vert \sqrt{t}\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2}+ t\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2}\\ &\quad\leq \biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2} +C_{\max} T \biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2}\\ &\qquad{} + C_{\max} T. \end{aligned} $$

(4.7)

By (4.2) the integrals on the right-hand side of (4.7) can be estimated as

$$\begin{aligned} \biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2} &=\biggl\Vert D^{\frac{1}{2}}(u_{\varepsilon}+ \varepsilon)\nabla \int _{0}^{u_{\varepsilon}+\varepsilon}D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2} \\ &\leq c_{\infty}\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2} \\ &\leq c_{\infty}\biggl(\frac{1}{2}\|u_{0} \|^{2}_{L^{2}(\Omega)}+C_{1} T \biggr) \end{aligned}$$

(4.8)

Then substituting (4.8) into (4.7) and using (4.2) again, we have

$$\begin{aligned} &\biggl\Vert \sqrt{t}\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz\biggr\Vert _{L^{2} (0,T; L^{2}(\Omega) )}^{2}+ t\biggl\Vert \nabla \int_{0}^{u_{\varepsilon}+\varepsilon}D(z)\,dz\biggr\Vert _{L^{2}(\Omega)}^{2} \\ &\quad\leq c_{\infty}\biggl(\frac{1}{2}\|u_{0} \|^{2}_{L^{2}(\Omega)}+C_{1} T \biggr)+C_{\max} T \biggl(\frac{1}{2}\|u_{0}\|^{2}_{L^{2}(\Omega)}+C_{1} T \biggr)+ C_{\max} T \\ &\quad\leq C_{2}+C_{2} T +C_{2} T^{2}, \end{aligned}$$

where $C_{2}:=\max \{\frac{1}{2}c_{\infty}\|u_{0}\|^{2}_{L^{2}(\Omega)}, C_{\max}C_{1}, [\frac{1}{2}C_{\max}\|u_{0}\|^{2}_{L^{2}(\Omega)}+C_{\max }+c_{\infty}C_{1} ] \}$. By taking the supremum with respect to t on $(0, T)$ we complete the proof of (4.4). □

We now prove Theorem 1.2. Our method is also partially inspired by [18].

Proof of Theorem 1.2

For any given $T>0$, we have $\| u_{\varepsilon}\|_{L^{\infty} (0,T; L^{p}(\Omega) )}< C$ ($p\in[1,\infty]$), where C is a positive constant independent of T and ε. Then there exist a subsequence $\{u_{\varepsilon_{j}} \}_{j\in\mathbb{N}}$ and a function $u\in L^{\infty} (0,T; L^{p}(\Omega) )$ such that

$$ u_{\varepsilon_{j}}\rightharpoonup u \quad\mbox{weakly}^{*}\mbox{ in } L^{\infty} \bigl(0,T; L^{p}(\Omega) \bigr) $$

(4.9)

for any $p\in[1,\infty]$, where $\varepsilon_{j}\rightarrow0$ as $j\rightarrow\infty$. By using $D\in C^{2} ([0,\infty) )$ and $\| u_{\varepsilon}(t)\|_{L^{\infty}(\Omega)}< c_{1}$ again, from Lemma 4.1 we deduce that $\int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz$ is bounded in $L^{2} (0,T; H^{1}(\Omega) )$. Hence, there exist a subsequence (still denoted by $\{u_{\varepsilon_{j}}\} _{j\in\mathbb{N}}$) and a function $\vartheta\in L^{2} (0,T; H^{1}(\Omega) )$ such that

$$ \begin{aligned} & \int_{0}^{u_{\varepsilon_{j}}+\varepsilon_{j}} D^{\frac {1}{2}}(z)\,dz\rightharpoonup \vartheta \quad\mbox{weakly in } L^{2} \bigl(0,T; L^{2}( \Omega) \bigr), \\ &\nabla \int_{0}^{u_{\varepsilon_{j}}+\varepsilon_{j} }D^{\frac {1}{2}}(z)\,dz\rightharpoonup \nabla\vartheta \quad\mbox{weakly in } L^{2} \bigl(0,T; L^{2}( \Omega) \bigr). \end{aligned} $$

(4.10)

On the other hand, by letting $\tau>0$, from Lemma 4.2 we have

$$\begin{aligned} \tau\biggl\Vert \frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz\biggr\Vert _{L^{2} (\tau,T; L^{2}(\Omega) )}^{2} &\leq\biggl\Vert \sqrt{t}\frac{d}{dt} \int_{0}^{u_{\varepsilon}+\varepsilon }D^{\frac{1}{2}}(z)\,dz\biggr\Vert _{L^{2} (\tau,T; L^{2}(\Omega) )}^{2} \\ &\leq C_{2}+C_{2} T +C_{2} T^{2}, \end{aligned}$$

which implies that $\int_{0}^{u_{\varepsilon}+\varepsilon}D^{\frac {1}{2}}(z)\,dz$ is bounded in $H^{1} (\tau,T; L^{2}(\Omega) )$ (in particular, it is bounded in $H^{1} (\tau,T; H^{-1}(\Omega) )$. Thus, by the Aubin-Lions lemma there exists a subsequence (still denoted by $\{u_{\varepsilon_{j}}\}_{j\in\mathbb{N}}$) such that

$$\int_{0}^{u_{\varepsilon_{j}}+\varepsilon_{j}} D^{\frac {1}{2}}(z)\,dz\rightarrow \vartheta\quad \mbox{strongly in } L^{2} \bigl(\tau,T; L^{2}( \Omega) \bigr) \mbox{ and a.e. on } \Omega\times(\tau,T). $$

Set $f(r):=\int_{0} ^{r} D^{\frac{1}{2}}(z)\,dz$. We see that $f(r)$ is a strictly increasing and continuous function. Thus, the inverse function $f^{-1}(r)$ of f exists and is continuous. Moreover, we can obtain that

$$ u_{\varepsilon_{j}}\rightarrow u=f^{-1}(\vartheta) \quad \mbox{strongly in } L^{2} \bigl(\tau,T; L^{2}(\Omega) \bigr) \mbox{ and a.e. on } \Omega\times(\tau,T). $$

(4.11)

Since $\tau>0$ is arbitrary, we deduce from (4.10) and (4.11) that

$$ \vartheta= \int_{0} ^{u} D^{\frac{1}{2}}(z)\,dz\in L^{2} \bigl(0,T; H^{1}(\Omega) \bigr). $$

(4.12)

Since $\|v_{\varepsilon}(t)\|_{W^{1,\infty}(\Omega)}< c_{2}+c_{4}$, there exist a subsequence $\{v_{\varepsilon_{j}}\}_{j \in\mathbb{N}}$ (hereafter, we still denote the subscript of the subsequence by $\{ v_{\varepsilon_{j}}\}_{j\in\mathbb{N}}$ for simplicity) and functions v such that

$$ \begin{aligned} & v_{\varepsilon_{j}}\rightharpoonup v \quad\mbox{weakly}^{*}\mbox{ in } L^{\infty} \bigl(0,T; L^{\infty}( \Omega) \bigr), \\ & \nabla v_{\varepsilon_{j}}\rightharpoonup\nabla v \quad\mbox{weakly}^{*}\mbox{ in } L^{\infty} \bigl(0,T; L^{\infty}(\Omega) \bigr). \end{aligned} $$

(4.13)

Similarly, there exist subsequence $\{w_{\varepsilon_{j}}\}_{n\in\mathbb {N}}$ and functions w such that

$$ \begin{aligned} & w_{\varepsilon_{j}}\rightharpoonup w \quad\mbox{weakly}^{*}\mbox{ in } L^{\infty} \bigl(0,T; L^{\infty}( \Omega) \bigr), \\ & \nabla w_{\varepsilon_{j}}\rightharpoonup\nabla w \quad\mbox{weakly}^{*}\mbox{ in } L^{\infty} \bigl(0,T; L^{\infty}(\Omega) \bigr) \end{aligned} $$

(4.14)

due to $\|w_{\varepsilon}(t)\|_{W^{1,\infty}(\Omega)}< c_{3}+c_{5}$.

For any given $T\in(0,\infty)$, we take $\varphi\in C_{0} ^{\infty} (\Omega\times[0,T) )$. Then multiplying the first, second, and third equations in (4.1) by φ and integrating those on $\Omega \times(0,T)$ we see that

$$ \left \{ \textstyle\begin{array}{@{}l} \int_{0} ^{T}\int_{\Omega} (D(u_{\varepsilon_{j}}+\varepsilon_{j})\nabla u_{\varepsilon_{j}}\cdot\nabla\varphi-\chi u_{\varepsilon_{j}} \nabla v_{\varepsilon_{j}}\cdot\nabla\varphi+\xi u_{\varepsilon_{j}} \nabla w_{\varepsilon_{j}}\cdot\nabla\varphi-u_{\varepsilon_{j}}\varphi_{t} )\,dx\,dt\\ \quad=\int_{\Omega}u_{0}(x)\varphi(x,0)\,dx,\\ \int_{0} ^{T}\int_{\Omega} (\nabla v_{\varepsilon_{j}}\cdot\nabla\varphi +\beta v_{\varepsilon_{j}} \varphi )\,dx\,dt=\int_{0} ^{T}\int_{\Omega} \alpha u_{\varepsilon_{j}} \varphi \,dx\,dt,\\ \int_{0} ^{T}\int_{\Omega} (\nabla w_{\varepsilon_{j}}\cdot\nabla\varphi +\delta w_{\varepsilon_{j}} \varphi )\,dx\,dt=\int_{0} ^{T}\int_{\Omega} \gamma u_{\varepsilon_{j}} \varphi \,dx\,dt. \end{array}\displaystyle \right . $$

(4.15)

Noting that $D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon)|\nabla\varphi |\leq c^{\frac{1}{2}}_{\infty}\|\nabla\varphi\|_{L^{\infty}}$ and thus $D^{\frac{1}{2}}(u_{\varepsilon}+\varepsilon)|\nabla\varphi|\in L^{2} (0,T; L^{2}(\Omega) )$, we see from (4.11) that

$$D^{\frac{1}{2}}(u_{\varepsilon_{j}}+\varepsilon_{j})\nabla\varphi \rightarrow D^{\frac{1}{2}}(u)\nabla\varphi \quad\mbox{strongly in } L^{2} \bigl(0,T; L^{2}(\Omega) \bigr), $$

which, together with (4.10) and (4.12), yields that

$$\begin{aligned} &\int_{0} ^{T} \int_{\Omega} \bigl(D(u_{\varepsilon_{j}}+\varepsilon_{j}) \nabla u_{\varepsilon_{j}}\cdot\nabla\varphi \bigr)\,dx\,dt \\ &\quad= \int_{0} ^{T} \int_{\Omega} \biggl(D^{\frac{1}{2}}(u_{\varepsilon _{j}}+ \varepsilon_{j})\nabla\varphi\cdot\nabla \int_{0}^{u_{\varepsilon _{j}}+\varepsilon_{j}}D^{\frac{1}{2}}(z)\,dz \biggr)\,dx\,dt \\ &\quad\rightarrow \int_{0} ^{T} \int_{\Omega} \biggl(D^{\frac{1}{2}}(u)\nabla \varphi\cdot\nabla \int_{0}^{u}D^{\frac{1}{2}}(z)\,dz \biggr)\,dx\,dt \\ &\quad= \int_{0} ^{T} \int_{\Omega} \bigl(D(u)\nabla u\cdot\nabla\varphi \bigr)\,dx\,dt \end{aligned}$$

(4.16)

as $j\to\infty$. Similarly, since

$$u_{\varepsilon_{j}} \nabla\varphi\rightarrow u\nabla\varphi \quad\mbox{strongly in } L^{2} \bigl(0,T; L^{2}(\Omega) \bigr) $$

by (4.11), from (4.13) and (4.14) we see that

$$\begin{aligned} &\int_{0} ^{T} \int_{\Omega} (-\chi u_{\varepsilon_{j}} \nabla v_{\varepsilon _{j}} \cdot\nabla\varphi+\xi u_{\varepsilon_{j}} \nabla w_{\varepsilon _{j}}\cdot\nabla \varphi )\,dx\,dt \\ &\quad\rightarrow \int_{0} ^{T} \int_{\Omega} (-\chi u\nabla v\cdot\nabla\varphi+\xi u\nabla w\cdot \nabla\varphi )\,dx\,dt \end{aligned}$$

(4.17)

as $j\to\infty$. Summarily, by collecting (4.9), (4.13), (4.14), (4.16), and (4.17), from (4.15) we obtain that

$$ \left \{ \textstyle\begin{array}{@{}l} \int_{0} ^{T}\int_{\Omega} (D(u)\nabla u \cdot\nabla\varphi-\chi u \nabla v\cdot\nabla\varphi+\xi u \nabla w\cdot\nabla\varphi-u \varphi _{t} )\,dx\,dt\\ \quad=\int_{\Omega}u_{0}(x)\varphi(x,0)\,dx,\\ \int_{0} ^{T}\int_{\Omega} (\nabla v\cdot\nabla\varphi+\beta v \varphi )\,dx\,dt=\int_{0} ^{T}\int_{\Omega} \alpha u \varphi \,dx\,dt,\\ \int_{0} ^{T}\int_{\Omega} (\nabla w\cdot\nabla\varphi+\delta w \varphi )\,dx\,dt=\int_{0} ^{T}\int_{\Omega} \gamma u \varphi \,dx\,dt \end{array}\displaystyle \right . $$

(4.18)

upon letting $j\rightarrow\infty$. Hence, $(u,v,w)$ is a global weak solution to system (1.2). Moreover, we deduce from (4.9), (4.13), (4.14), and Theorem 1.1 that

$$\begin{aligned}& \| u\|_{L^{\infty} (0,T; L^{\infty}(\Omega) )}\leq\liminf_{j\rightarrow\infty}\| u_{\varepsilon_{j}} \|_{L^{\infty} (0,T; L^{\infty}(\Omega) )}\leq c_{1}, \\& \| v\|_{L^{\infty} (0,T; L^{\infty}(\Omega) )}\leq\liminf_{j\rightarrow\infty}\| v_{\varepsilon_{j}} \|_{L^{\infty} (0,T; L^{\infty}(\Omega) )}\leq c_{2}, \\& \|w\|_{L^{\infty} (0,T; L^{\infty}(\Omega) )}\leq\liminf_{j\rightarrow\infty}\| w_{\varepsilon_{j}} \|_{L^{\infty} (0,T; L^{\infty}(\Omega) )}\leq c_{3}, \end{aligned}$$

which implies the uniform boundedness of $(u,v,w)$. Thus, we complete the proof of Theorem 1.2. □