1. Introduction

In recent years there has been an increased focus on the mathematical modelling and analysis of tumour growth. Many new models have been proposed and numerical simulations have been carried out to provide new and important insights on cancer research, see for instance ^[8] and [13, Chap. 3]. In this work we analyse a diffuse interface model proposed in ^[20], which models a mixture of tumour cells and healthy cells in the presence of an unspecified chemical species acting as a nutrient. More precisely, for a bounded domain $\Omega \subset {\mathbb{R}^d}$ where the cells reside and T > 0, we consider the following set of equations,

$\begin{array}{*{20}{c}} {{\rm{div }}v = {\Gamma _v}}&{{\rm{in }}\Omega \times (0,T) = :Q,} \end{array}$

(1.1a)

$v = - K(\nabla p - (\mu + \chi \sigma )\nabla \varphi ){\text{ in }}Q,$

(1.1b)

${\partial _t}\varphi + {\text{div}}(\varphi v) = {\text{div}}(m(\varphi )\nabla \mu ) + {\Gamma _\varphi }{\text{ in }}Q,$

(1.1c)

$\mu = A\Psi '(\varphi ) - B\Delta \varphi - \chi \sigma {\text{ in }}Q,$

(1.1d)

${\partial _t}\sigma + {\text{div}}(\sigma v) = {\text{div}}(n(\varphi )(D\nabla \sigma - \chi \nabla \varphi )) - \mathcal{S}{\text{ in }}Q.$

(1.1e)

Here, ${v}$ denotes the volume-averaged velocity of the mixture, $p$ denotes the pressure, $\sigma$ denotes the concentration of the nutrient, $\varphi \in [-1, 1]$ denotes the difference in volume fractions, with $\{\varphi = 1\}$ representing the unmixed tumour tissue, and $\{\varphi = -1\}$ representing the surrounding healthy tissue, and $\mu$ denotes the chemical potential for $\varphi$.

The model treats the tumour and healthy cells as inertia-less fluids, leading to the appearance of a Darcy-type subsystem with a source term $\Gamma_{{v}}$. The order parameter $\varphi$ satisfies a convective Cahn--Hilliard type equation with additional source term $\Gamma_{\varphi}$, and similarly, the nutrient concentration $\sigma$ satisfies a convection-reaction-diffusion equation with a non-standard flux and a source term $\mathcal{S}$. We refer the reader to [20, x2] for the derivation from thermodynamic principles, and to [20, x2.5] for a discussion regarding the choices for the source terms $\Gamma_{\varphi}, \Gamma_{{v}}$ and $\mathcal{S}$.

The positive constants $K$ and $D$ denote the permeability of the mixture and the diffusivity of the nutrient, $m(\varphi)$ and $n(\varphi)$ are positive mobilities for $\varphi$ and $\sigma$, respectively. The parameter $\chi \geq 0$ regulates the chemotaxis effect (see ^[20] for more details), $\Psi(\cdot)$ is a potential with two equal minima at $\pm 1$, $A$ and $B$ denote two positive constants related to the thickness of the diffuse interface and the surface tension.

We supplement the above with the following boundary and initial conditions

${\partial _n}\varphi = {\partial _n}\mu = 0{\text{ on }}\partial \Omega \times (0, T) = :\Sigma ,$

(1.2a)

$v \cdot n = {\partial _n}p = 0{\text{ on }}\Sigma ,$

(1.2b)

$n(\varphi )D{\partial _n}\sigma = b({\sigma _\infty } - \sigma ){\text{ on }}\Sigma ,$

(1.2c)

$\varphi (0) = {\varphi _0}, \quad \sigma (0) = {\sigma _0}{\text{ on }}\Omega .$

(1.2d)

Here $\varphi_{0}$, $\sigma_{0}$ and $\sigma_{\infty}$ are given functions and $b > 0$ is a constant. We denote ${\partial _n} f := \nabla f \cdot {n}$ as the normal derivative of $f$ at the boundary $\partial \Omega$, where ${n}$ is the outer unit normal. Associated to (1.1) is the free energy density $N(\varphi, \sigma)$ for the nutrient, which is defined as

$N(\varphi , \sigma ): = \frac{D}{2}{\left| \sigma \right|^2} + \chi \sigma (1 - \varphi ).$

(1.3)

Note that

${N_{, \sigma }}: = \frac{{\partial N}}{{\partial \sigma }} = D\sigma + \chi (1 - \varphi ), \quad {N_{, \varphi }}: = \frac{{\partial N}}{{\partial \varphi }} = - \chi \sigma ,$

so that (1.1) may also be written as

${\text{div}} v = {\Gamma _v},$

(1.4a)

$v = - K(\nabla p - \mu \nabla \varphi + {N_{, \varphi }}\nabla \varphi ),$

(1.4b)

${\partial _t}\varphi + {\text{div}}(\varphi v) = {\text{div}}(m(\varphi )\nabla \mu ) + {\Gamma _\varphi },$

(1.4c)

$\mu = A\Psi '(\varphi ) - B\Delta \varphi + {N_{, \varphi }},$

(1.4d)

${\partial _t}\sigma + {\text{div}}(\sigma v) = {\text{div}}(n(\varphi )\nabla {N_{, \sigma }}) - \mathcal{S},$

(1.4e)

which is the general phase field model proposed in ^[20]. In this work we do not aim to analyse such a model with a general free energy density $N(\varphi, \sigma)$, but we will focus solely on the choice (1.3) and the corresponding model (1.1)-(1.2).

Our goal in this work is to prove the existence of weak solutions (see Definition 2.1 below) of (1.1)-(1.2) in two and three dimensions. Moreover, one might expect that by setting $\Gamma_{{v}} = 0$ and then sending $b \to 0$ and $K \to 0$, the weak solutions to (1.1)-(1.2) will converge (in some appropriate sense) to the weak solutions of

${\partial _t}\varphi = {\text{div}}(m(\varphi )\nabla \mu ) + {\Gamma _\varphi }{\text{ in }}Q,$

(1.5a)

$\mu = A\Psi '(\varphi ) - B\Delta \varphi - \chi \sigma {\text{ in }}Q,$

(1.5b)

${\partial _t}\sigma = {\text{div}}(n(\varphi )(D\nabla \sigma - \chi \nabla \varphi )) - \mathcal{S}{\text{ in }}Q,$

(1.5c)

$0 = {\partial _n}\varphi = {\partial _n}\mu = {\partial _n}\sigma {\text{ on }}\Sigma .$

(1.5d)

We denote (1.5) as the limit system of vanishing permeability, where the effects of the volume-averaged velocity are neglected. By substituting

$\Gamma_{\varphi} = \mathcal{S} = f(\varphi)(D \sigma + \chi(1-\varphi) - \mu)$

(1.6)

for some non-negative function $f(\varphi)$ leads to the model derived in ^[21]. The specific choices for $\Gamma_{\varphi}$ and $\mathcal{S}$ in (1.6) are motivated by linear phenomenological laws for chemical reactions. The analysis of (1.5) with the parameters

$D = 1, \quad \chi = 0, \quad n(\varphi) = m(\varphi) = 1$

has been the subject of study in ^{[5, 6, 7, 16]}, where well-posedness and long-time behaviour have been established for a large class of functions $\Psi(\varphi)$ and $f(\varphi)$. Alternatively, one may consider the following choice of source terms

$\Gamma_{\varphi} = h(\varphi)(\lambda_{p} \sigma - \lambda_{a}), \quad \mathcal{S} = \lambda_{c} h(\varphi) \sigma,$

(1.7)

where $\lambda_{p}$, $\lambda_{a}$, $\lambda_{c}$ are non-negative constants representing the tumour proliferation rate, the apoptosis rate, and the nutrient consumption rate, respectively, and $h(\varphi)$ is a non-negative interpolation function such that $h(-1) = 0$ and $h(1) = 1$. The above choices for $\Gamma_{\varphi}$ and $\mathcal{S}$ are motivated from the modelling of processes experienced by a young tumour.

The well-posedness of model (1.5) with the choice (1.7) has been studied by the authors in ^[17] and ^[18] with the boundary conditions (1.2) (neglecting (1.2b)) in the former and for non-zero Dirichlet boundary conditions in the latter. It has been noted in ^[17] that the well-posedness result with the boundary conditions (1.2) requires $\Psi$ to have at most quadratic growth, which is attributed to the presence of the source term $\Gamma_{\varphi} \mu = h(\varphi) \mu (\lambda_{p} \sigma - \lambda_{a})$ when deriving useful a priori estimates. Meanwhile in ^[18] the Dirichlet boundary conditions and the application of the Poincaré inequality allows us to overcome this restriction and allow for $\Psi$ to be a regular potential with polynomial growth of order less than 6, and by a Yosida approximation, the case where $\Psi$ is a singular potential is also covered.

We also mention the work of ^[19] that utilises a Schauder's fixed point argument to show existence of weak solutions for $\Psi$ with quartic growth and $\Gamma_{\varphi}, \mathcal{S}$ as in (1.7). This is based on first deducing that $\sigma$ is bounded by a comparison principle, leading to $\Gamma_{\varphi} \in L^{\infty}(\Omega)$. Then, the standard a priori estimates are derived for a Cahn--Hilliard equation with bounded source terms. The difference between ^[19] and ^{[17, 18]} is the absence of the chemotaxis and active transport mechanisms, i.e., $\chi = 0$, so that the comparison principle can be applied to the nutrient equation. We refer to ^[9] for the application of a similar procedure to a multi-species tumour model with logarithmic potentials.

On the other hand, by sending $b \to 0$ and $\chi \to 0$ in (1.1), we should obtain weak solutions of

${\text{div}} v = {\Gamma _v}{\text{ in }}Q,$

(1.8a)

$v = - K(\nabla p - \mu \nabla \varphi ){\text{ in }}Q,$

(1.8b)

${\partial _t}\varphi + {\text{div}} (\varphi v) = {\text{div}}(m(\varphi )\nabla \mu ) + {\Gamma _\varphi }{\text{ in }}Q,$

(1.8c)

$\mu = A\Psi '(\varphi ) - B\Delta \varphi {\text{ in }}Q,$

(1.8d)

${\partial _t}\sigma + {\text{div}}(\sigma v) = {\text{div}}(n(\varphi )D\nabla \sigma ) - \mathcal{S}{\text{ in }}Q,$

(1.8e)

$0 = {\partial _n}\varphi = {\partial _n}\mu = {\partial _n}\sigma = v \cdot n{\text{ on }}\Sigma .$

(1.8f)

We denote (1.8) as the limit system of vanishing chemotaxis. If the source terms $\Gamma_{{v}}$ and $\Gamma_{\varphi}$ are independent of $\sigma$, then (1.8) consists of an independent Cahn--Hilliard--Darcy system and an equation for $\sigma$ which is advected by the volume-averaged velocity field ${v}$. In the case where there is no nutrient and source terms, i.e., $\sigma = \Gamma_{{v}} = \Gamma_{\varphi} = 0$, global existence of weak solutions in two and three dimensions has been established in ^[14] via the convergence of a fully discrete and energy stable implicit finite element scheme. For the well-posedness and long-time behaviour of strong solutions, we refer to ^[25]. Meanwhile, in the case where $\Gamma_{{v}} = \Gamma_{\varphi}$ is prescribed, global weak existence and local strong well-posedness for (1.8) without nutrient is shown in ^[22].

We also mention the work of ^[3] on the well-posedness and long-time behaviour of a related system also used in tumour growth, known as the Cahn--Hilliard--Brinkman system, where in (1.8) without nutrient an additional viscosity term is added to the left-hand side of the velocity equation (1.8b) and the mass exchange terms $\Gamma_{{v}}$ and $\Gamma_{\varphi}$ are set to zero. The well-posedness of a nonlocal variant of the Cahn--Hilliard--Brinkman system has been investigated in ^[10]. Furthermore, when $K$ is a function depending on $\varphi$, the model (1.8) with $\sigma = \Gamma_{{v}} = \Gamma_{\varphi} = 0$ is also referred to as the Hele--Shaw--Cahn--Hilliard model (see ^{[23, 24]}). In this setting, $K(\varphi)$ represents the reciprocal of the viscosity of the fluid mixture. We refer to ^[30] concerning the strong well-posedness globally in time for two dimensions and locally in time for three dimensions when $\Omega$ is the $d$-dimensional torus. Global well-posedness in three dimensions under additional assumptions and long-time behaviour of solutions to the Hele--Shaw--Cahn--Hilliard model are investigated in ^[29].

We point out that from the derivation of (1.1) in ^[20], the source terms $\Gamma_{{v}}$ and $\Gamma_{\varphi}$ are connected in the sense that $\Gamma_{{v}}$ is related to sum of the mass exchange terms for the tumour and healthy cells, and $\Gamma_{\varphi}$ is related to the difference between the mass exchange terms. Thus, if $\Gamma_{\varphi}$ would depend on the primary variables $\varphi$, $\sigma$ or $\mu$, then one expects that $\Gamma_{{v}}$ will also depend on the primary variables. Here, we are able to prove existence of weak solutions for $\Gamma_{\varphi}$ of the form (2.1), which generalises the choices (1.6) and (1.7), but in exchange $\Gamma_{{v}}$ has to be considered as a prescribed function. This is attributed to the presence of the source term ${\Gamma _v}\left( {\varphi \mu + \frac{D}{2}{{\left| \sigma \right|}^2}} \right)$ when deriving useful a priori estimates. We see that if $\Gamma_{{v}}$ depends on the primary variables, we obtain triplet products which cannot be controlled by the usual regularity of $\varphi$, $\mu$ and $\sigma$ in the absence of a priori estimates.

In this work we attempt to generalise the weak existence results for the models studied in ^{[5, 16, 17, 18, 22, 25]} by proving that the weak solutions of (1.1) with $\Gamma_{{v}} = 0$ converge (in some appropriate sense) to the weak solutions of (1.5) as $b \to 0$ and $K \to 0$, and the weak solutions of (1.1) converge to the weak solutions of (1.8) as $b \to 0$ and $\chi \to 0$.

This paper is organised as follows. In Section 2 we state the main assumptions and the main results. In Section 3 we introduce a Galerkin procedure and derive some a priori estimates for the Galerkin ansatz in Section 4 for the case of three dimensions. We then pass to the limit in Section 5 to deduce the existence result for three dimensions, while in Section 6 we investigate the asymptotic behaviour of solutions to (1.1) as $K \to 0$ and $\chi \to 0$. In Section 7, we outline the a priori estimates for two dimensions and show that the weak solutions for two dimensions yields better temporal regularity than the weak solutions for three dimensions. In Section 8 we discuss some of the issues present in the analysis of (1.1) using different formulations of Darcy's law and the pressure, and with different boundary conditions for the velocity and the pressure.

Notation. For convenience, we will often use the notation $L^{p} := L^{p}(\Omega)$ and $W^{k, p} := W^{k, p}(\Omega)$ for any $p \in [1, \infty]$, $k > 0$ to denote the standard Lebesgue spaces and Sobolev spaces equipped with the norms ${\left\| \cdot \right\|_{{L^p}}}$ and $\left\| \cdot \right\|_{W^{k, p}}$. In the case $p = 2$ we use $H^{k} := W^{k, 2}$ and the norm $\left\| \cdot \right\|_{H^{k}}$. For the norms of Bochner spaces, we will use the notation $L^{p}(X) := L^{p}(0, T;X)$ for Banach space $X$ and $p \in [1, \infty]$. Moreover, the dual space of a Banach space $X$ will be denoted by $X^{*}$, and the duality pairing between $X$ and $X^{*}$ is denoted by ${\left\langle { \cdot , \cdot } \right\rangle _{X, {X^*}}}$. For $d = 2$ or $3$, let $\mathcal{H}^{d-1}$ denote the $(d-1)$ dimensional Hausdorff measure on $\partial \Omega$, and we denote ${\mathbb{R}^d}$-valued functions and any function spaces consisting of vector-valued/tensor-valued functions in boldface. We will use the notation D${f}$ to denote the weak derivative of the vector function ${f}$.

Useful preliminaries. For convenience, we recall the Poincaré inequality: There exists a positive constant $C_{p}$ depending only on $\Omega$ such that, for all $f \in H^{1}$,

${\left\| {f - \bar f} \right\|_{{L^2}}} \leq {C_p}{\left\| {\nabla f} \right\|_{{L^2}}},$

(1.9)

where $\bar f: = \frac{1}{{\left| \Omega \right|}}\int_\Omega f {\text{ dx}}$ denotes the mean of $f$. The Gagliardo--Nirenberg interpolation inequality in dimension $d$ is also useful (see [15, Thm. 10.1, p. 27], [11, Thm. 2.1] and [1, Thm. 5.8]): Let $\Omega$ be a bounded domain with Lipschitz boundary, and $f \in W^{m, r} \cap L^{q}$, $1 \leq q, r \leq \infty$. For any integer $j$, $0 \leq j < m$, suppose there is $\alpha \in \mathbb{R}$ such that

$\frac{1}{p} = \frac{j}{d} + \left ( \frac{1}{r} - \frac{m}{d} \right ) \alpha + \frac{1-\alpha}{q}, \quad \frac{j}{m} \leq \alpha \leq 1.$

Then, there exists a positive constant $C$ depending only on $\Omega$, $m$, $j$, $q$, $r$, and $\alpha$ such that

${\left\| {{D^j}f} \right\|_{{L^p}}} \leq C\left\| f \right\|_{{W^{m, r}}}^\alpha \left\| f \right\|_{{L^q}}^{1 - \alpha }.$

(1.10)

We will also use the following Gronwall inequality in integral form (see [17, Lem. 3.1] for a proof): Let $\alpha, \beta, u$ and $v$ be real-valued functions defined on $[0, T]$. Assume that $\alpha$ is integrable, $\beta$ is non-negative and continuous, $u$ is continuous, $v$ is non-negative and integrable. If $u$ and $v$ satisfy the integral inequality

$u(s) + \int_0^s v (t){\text{ dt}} \leq \alpha (s) + \int_0^s \beta (t)u(t){\text{ dt}}\quad {\text{ for }}s \in (0, T],$

then it holds that

$u(s) + \int_0^s v (t){\text{ dt}} \leq \alpha (s) + \int_0^s \beta (t)\alpha (t)\exp \left( {\int_0^t \beta (r){\text{ dr}}} \right){\text{ dt}}.$

(1.11)

To analyse the Darcy system, we introduce the spaces

$\begin{gathered} L_0^2: = \{ f \in {L^2}:\bar f = 0\} , \quad H_N^2: = \{ f \in {H^2}:{\partial _n}f = 0{\text{ on }}\partial \Omega \} , \hfill \\ ({H^1})_0^*: = \{ f \in {({H^1})^*}:{\left\langle {f, 1} \right\rangle _{{H^1}}} = 0\} . \hfill \\ \end{gathered}$

Then, the Neumann-Laplacian operator $-\Delta _{N} : H^{1} \cap L^{2}_{0} \to (H^{1})^{*}_{0}$ is positively defined and self-adjoint. In particular, by the Lax--Milgram theorem and the Poincaré inequality (1.9) with zero mean, the inverse operator $(-\Delta _{N})^{-1} : (H^{1})^{*}_{0} \to H^{1} \cap L^{2}_{0}$ is well-defined, and we set $u := (-\Delta _{N})^{-1}f$ for $f \in (H^{1})^{*}_{0}$ if $\bar u = 0$ and

$- \Delta u = f{\text{ in }}\Omega , \quad {\text{ }}{\partial _n}u = 0{\text{ on }}\partial \Omega .$

---

2. Main results

We make the following assumptions.

Assumption 2.1.

(A1) The constants A, B, K, D, $\chi$ and b are positive and fixed.

(A2) The mobilities m, n are continuous on $\mathbb{R}$ and satisfy

${m_0} \leq m(t) \leq {m_1}, {\text{ }}\quad {n_0} \leq n(t) \leq {n_1}\quad \forall t \in \mathbb{R},$

for positive constants m₀, m₁, n₀ and n₁.

(A3) $\Gamma_{\varphi}$ and $\mathcal{S}$ are of the form

$\begin{gathered} {\Gamma _\varphi }(\varphi , \mu , \sigma ) = {\Lambda _\varphi }(\varphi , \sigma ) - {\Theta _\varphi }(\varphi , \sigma )\mu , \hfill \\ \mathcal{S}(\varphi , \mu , \sigma ) = {\Lambda _S}(\varphi , \sigma ) - {\Theta _S}(\varphi , \sigma )\mu , \hfill \\ \end{gathered}$

(2.1)

where $\Theta_{\varphi}, \Theta_{S} : \mathbb{R}^{2} \to \mathbb{R}$ are continuous bounded functions with $\Theta_{\varphi}$ non-negative, and $\Lambda_{\varphi}, \Lambda_{S} : \mathbb{R}^{2} \to \mathbb{R}$ are continuous with linear growth

$\left| {{\Theta _i}(\varphi , \sigma )} \right| \leq {R_0}, {\text{ }}\quad {\text{ }}\left| {{\Lambda _i}(\varphi , \sigma )} \right| \leq {R_0}(1 + \left| \varphi \right| + \left| \sigma \right|)\quad {\text{ for }}i \in \{ \varphi , S\} ,$

(2.2)

so that

$\left| {{\Gamma _\varphi }} \right| + \left| \mathcal{S} \right| \leq {R_0}(1 + \left| \varphi \right| + \left| \mu \right| + \left| \sigma \right|),$

(2.3)

for some positive constant R₀.

(A4) $\Gamma_{{v}}$ is a prescribed function belonging to $L^{4}(0, T;L^{2}_{0})$.

(A5) $\Psi \in C^{2}(\mathbb{R})$ is a non-negative function satisfying

$\Psi (t) \geq {R_1}{\text{ }}{\left| t \right|^2} - {R_2}\quad \forall t \in \mathbb{R}$

(2.4)

and either one of the following,

1: if $\Theta_{\varphi}$ is non-negative and bounded, then

$\Psi (t) \leq {R_3}(1 + {\left| t \right|^2}), {\text{ }}\quad {\text{ }}\left| {\Psi '(t)} \right| \leq {R_4}(1 + \left| t \right|), \quad {\text{ }}\left| {\Psi ''(t)} \right| \leq {R_4};$

(2.5)

2: if $\Theta_{\varphi}$ is positive and bounded, that is,

$R_{0} \geq \Theta_{\varphi}(t, s) \geq R_{5} > 0 \quad \forall t, s \in \mathbb{R},$

(2.6)

then

$\left| {\Psi ''(t)} \right| \leq {R_6}(1 + {\left| t \right|^q}), \;q \in [0, 4),$

(2.7)

for some positive constants R₁, R₂, R₃, R₄, R₅, R₆. Furthermore we assume that

$A > \frac{2 \chi^{2}}{D R_{1}} .$

(2.8)

(A6) The initial and boundary data satisfy

${\sigma _\infty } \in {L^2}(0, T;{L^2}(\partial \Omega )), \quad {\sigma _0} \in {L^2}, \quad {\varphi _0} \in {H^1}.$

We point out that some of the above assumptions are based on previous works on the well-posedness of Cahn--Hilliard systems for tumour growth. For instance, (2.5) and (2.8) reflect the situation encountered in ^[17], where if $\Theta_{\varphi} = 0$, i.e., $\Gamma_{\varphi}$ is independent of $\mu$, then the derivation of the a priori estimate requires a quadratic potential. But in the case where (2.6) is satisfied, we can allow $\Psi$ to be a regular potential with polynomial growth of order less than 6, and by a Yosida approximation, we can extend our existence results to the situation where $\Psi$ is a singular potential, see for instance ^[18]. Moreover, the condition (2.8) is a technical assumption based on the fact that the second term of the nutrient free energy $\chi \sigma (1-\varphi)$ does not have a positive sign.

Meanwhile, the linearity of the source terms $\Gamma_{\varphi}$ and $\mathcal{S}$ with respect to the chemical potential $\mu$ assumed in (2.1) is a technical assumption based on the expectation that, at best, we have weak convergence for Galerkin approximation to $\mu$, which is in contrast with $\varphi$ and $\sigma$ where we might expect a.e convergence and strong convergence for the Galerkin approximations. Moreover, if we consider

$\Theta_{\varphi}(\varphi, \sigma) = \Theta_{S}(\varphi, \sigma) = f(\varphi), \quad \Lambda_{\varphi}(\varphi, \sigma) = \Lambda_{S}(\varphi, \sigma) = f(\varphi)(D \sigma + \chi(1-\varphi)),$

for a non-negative function $f(\varphi)$, then we obtain the source terms in ^{[5, 16, 21]}.

Compared to the set-up in ^[22], in (A4) we prescribe a higher temporal regularity for the prescribed source term $\Gamma_{{v}}$. This is needed when we estimate the source term ${\Gamma _v}\frac{D}{2}{\left| \sigma \right|^2}$ in the absence of a priori estimates, see Section 4.1.2 for more details. The mean zero condition is a consequence of the no-flux boundary condition ${v} \cdot {n} = 0$ on $\partial \Omega$ and the divergence equation (1.1a). In particular, we can express the Darcy subsystem (1.1a)-(1.1b) as an elliptic equation for the pressure $p$:

$- \Delta p = \frac{1}{K}{\Gamma _v} - {\text{div}}((\mu + \chi \sigma )\nabla \varphi ){\text{ in }}\Omega ,$

(2.9a)

${\partial _n}p = 0{\text{ on }}\partial \Omega .$

(2.9b)

Solutions to (2.9) are uniquely determined up to an arbitrary additive function that may only depend on time, and thus without loss of generality, we impose the condition $\bar p = \frac{1}{{\left| \Omega \right|}}\int_\Omega p {\text{ dx}} = 0$ to (2.9). We may then define $p$ as

$p = {( - {\Delta _N})^{ - 1}}\left( {\frac{1}{K}{\Gamma _v} - {\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right),$

(2.10)

if $\frac{1}{K}{\Gamma _v} - {\text{div}}((\mu + \chi \sigma )\nabla \varphi ) \in ({H^1})_0^*$.

Remark 2.1. In the case $\Gamma_{{v}} = 0$, one can also consider the assumption

$\mathcal{S} N_{, \sigma} - \Gamma_{\varphi} \mu = \mathcal{S} (D \sigma + \chi(1-\varphi)) - \Gamma_{\varphi} \mu \geq 0$

(2.11)

instead of (2.6), which holds automatically if $\Gamma_{\varphi}$ and $\mathcal{S}$ are chosen to be of the form (1.6). In fact this property is used in ^{[5, 16]}.

We make the following definition.

Definition 2.1 (Weak solutions for 3D). We call a quintuple $(\varphi, \mu, \sigma, {v}, p)$ a weak solution to (1.1)-(1.2) if

$\begin{gathered} \varphi \in {L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, \frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ \sigma \in {L^\infty }(0, T;{L^2}) \cap {L^2}(0, T;{H^1}) \cap {W^{1, \frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ \mu \in {L^2}(0, T;{H^1}), \quad p \in {L^{\frac{8}{5}}}(0, T;{H^1} \cap L_0^2), \quad v \in {L^2}(0, T;{L^2}), \hfill \\ \end{gathered}$

such that $\varphi(0) = \varphi_{0}$,

${\left\langle {{\sigma _0}, \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = {\left\langle {\sigma (0), \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}}\quad \forall \zeta \in {H^1},$

and

${\left\langle {{\partial _t}\varphi, \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \int_\Omega - m(\varphi )\nabla \mu \cdot \nabla \zeta + {\Gamma _\varphi }\zeta + \varphi v \cdot \nabla \zeta {\text{ dx}},$

(2.12a)

$\int_\Omega \mu \zeta {\text{ dx}} = \int_\Omega A \Psi '(\varphi )\zeta + B\nabla \varphi \cdot \nabla \zeta - \chi \sigma \zeta {\text{ dx}},$

(2.12b)

$\begin{gathered} {\left\langle {{\partial _t}\sigma, \phi } \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}} = \int_\Omega - n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla \phi - \mathcal{S}\phi + \sigma v \cdot \nabla \phi dx \\ + \int_{\partial \Omega } b ({\sigma _\infty } - \sigma )\phi {\text{ d}}{H^{d - 1}}, \\ \end{gathered}$

(2.12c)

$\int_\Omega \nabla p \cdot \nabla \zeta {\text{ dx}} = \int_\Omega {\frac{1}{K}} {\Gamma _v}\zeta + (\mu + \chi \sigma )\nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(2.12d)

$\int_\Omega {} v \cdot {\text{ }}\zeta {\text{ dx}} = \int_\Omega - K(\nabla p - (\mu + \chi \sigma )\nabla \varphi ) \cdot {\text{ }}\zeta {\text{ dx}},$

(2.12e)

for a.e. $t \in (0, T)$ and for all $\zeta \in H^{1}$, $\phi \in W^{1, 5}$, and ${\zeta} \in {L}^{2}$.

Neglecting the nutrient $\sigma$, we observe that our choice of function spaces for $(\varphi, \mu, p, {v})$ coincide with those in [22, Defn. 2.1(i)]. In contrast to the usual $L^{2}(0, T;(H^{1})^{*})$-regularity (see ^{[5, 16]}) we obtain a less regular time derivative $\partial _{t}\varphi$. The drop in the time regularity from $2$ to $\frac{8}{5}$ is attributed to the convection term ${\text{div}} (\varphi {v})$ belonging to $L^{\frac{8}{5}}(0, T;(H^{1})^{*})$. The same is true for the regularity for the time derivative $\partial _{t}\sigma$ in $L^{\frac{5}{4}}(0, T;(W^{1, 5})^{*})$ as the convection term ${\text{div}}(\sigma {v})$ lies in the same space. We refer the reader to the end of Section 4.3 for a calculation motivating the choice of function spaces for ${\text{div}} (\sigma {v})$ and $\partial _{t} \sigma$. Furthermore, the embedding of $L^{\infty}(0, T;H^{1}) \cap W^{1, \frac{8}{5}}(0, T;(H^{1})^{*})$ into $C^{0}([0, T];L^{2})$ from [28, x8, Cor. 4] guarantees that the initial condition for $\varphi$ is meaningful. However, for $\sigma$ we have the embedding $L^{\infty}(0, T;L^{2}) \cap W^{1, \frac{5}{4}}(0, T;(W^{1, 5})^{*}) \subset \subset C^{0}([0, T];(H^{1})^{*})$, and so $\sigma(0)$ makes sense as a function in $(H^{1})^{*}$. Thus, the initial condition $\sigma_{0}$ is attained as an equality in $(H^{1})^{*}$. We now state the existence result for (1.1)-(1.2).

Theorem 2.1 (Existence of weak solutions in 3D and energy inequality). Let $\Omega \subset \mathbb{R}^{3}$ be a bounded domain with $C^{3}$-boundary $\partial \Omega$. Suppose Assumption 2.1 is satisfied. Then, there exists a weak solution quintuple $(\varphi, \mu, \sigma, {v}, p)$ to (1.1)-(1.2) in the sense of Definition 2.1 with

$p \in L^{\frac{8}{7}}(0, T;H^{2}) , \quad {v} \in L^{\frac{8}{7}}(0, T;{H}^{1}),$

(2.13)

and in addition satisfies

$\begin{gathered} {\left\| \varphi \right\|_{{L^\infty }({H^1}) \cap {L^2}({H^3}) \cap {W^{1, \frac{8}{5}}}({{({H^1})}^*})}} + {\left\| \sigma \right\|_{{L^\infty }({L^2}) \cap {W^{1, \frac{5}{4}}}({{({W^{1, 5}})}^*})) \cap {L^2}({H^1})}} \hfill \\ + {\left\| \mu \right\|_{{L^2}({H^1})}} + {b^{\frac{1}{2}}}{\left\| \sigma \right\|_{{L^2}({L^2}(\partial \Omega ))}} + {\left\| p \right\|_{{L^{\frac{8}{5}}}({H^1}) \cap {L^{\frac{8}{7}}}({H^2})}} \hfill \\ + {K^{ - \frac{1}{2}}}\left( {{{\left\| v \right\|}_{{L^2}({L^2}) \cap {L^{\frac{8}{7}}}({H^1})}} + {{\left\| {{\text{div}}(\varphi v)} \right\|}_{{L^{\frac{8}{5}}}({{({H^1})}^*})}} + {{\left\| {{\text{div}}(\sigma v)} \right\|}_{{L^{\frac{5}{4}}}({{({W^{1, 5}})}^*})}}} \right) \hfill \\ \leq C, \hfill \\ \end{gathered}$

(2.14)

where the constant $C$ does not depend on $(\varphi, \mu, \sigma, {v}, p)$ and is uniformly bounded for $b, \chi \in (0, 1]$ and is also uniformly bounded for $K \in (0, 1]$ when $\Gamma_{{v}} = 0$.

The regularity result (2.13) is new compared to estimates for weak solutions in ^[22], which arises from a deeper study of the Darcy subsystem, and can be obtained even in the absence of the nutrient. We mention that higher regularity estimates for the pressure $p$ in $L^{2}(0, T;H^{2})$ and the velocity ${v}$ in $L^{2}(0, T;{H}^{1})$ are also established in [22, but these are for strong solutions local in time in three dimensions and global in time for two dimensions.

We now investigate the situation in two dimensions, where the Sobolev embeddings in two dimensions yields better integrability exponents.

Theorem 2.2 (Existence of weak solutions in 2D). Let $\Omega \subset \mathbb{R}^{2}$ be a bounded domain with $C^{3}$-boundary $\partial \Omega$. Suppose Assumption 2.1is satisfied. Then, there exists a quintuple $(\varphi, \mu, \sigma, {v}, p)$ to (1.1)-(1.2) with the following regularity

$\begin{gathered} \varphi \in {L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, w}}(0, T;{({H^1})^*}), \quad \mu \in {L^2}(0, T;{H^1}), \hfill \\ \sigma \in {L^2}(0, T;{H^1}) \cap {L^\infty }(0, T;{L^2}) \cap {W^{1, r}}(0, T;{({H^1})^*}), \hfill \\ p \in {L^k}(0, T;{H^1} \cap L_0^2) \cap {L^q}(0, T;{H^2}), \quad v \in {L^2}(0, T;{L^2}) \cap {L^q}(0, T;{H^1}), \hfill \\ \end{gathered}$

for

$1 \leq k < 2, \quad 1 \leq q < \frac{4}{3}, \quad 1 < r < \frac{8}{7}, \quad \frac{4}{3} \leq w < 2,$

such that (2.12a), (2.12b), (2.12d), (2.12e) and

${\left\langle {{\partial _t}\sigma, \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \int_\Omega - n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla \zeta - \mathcal{S}\zeta + \sigma v \cdot \nabla \zeta {\text{ dx}} + \int_{\partial \Omega } b ({\sigma _\infty } - \sigma )\zeta {\text{ d}}{H^{d - 1}}$

are satisfied for a.e. $t \in (0, T)$, for all $\zeta \in H^{1}$, and all ${\zeta} \in {L}^{2}$. Furthermore, the initial conditions $\varphi(0) = \varphi_{0}$ and $\sigma(0) = \sigma_{0}$ are attained as in Definition 2.1, and an analogous inequality to (2.14) also holds.

The proof of Theorem 2.2 is similar to that of Theorem 2.1, and hence the details are omitted. In Section 7 we will only present the derivation of a priori estimates. It is due to the better exponents for embeddings in two dimensions and the regularity result for the velocity that we obtain better regularities for the time derivatives $\partial _{t}\varphi$ and $\partial _{t}\sigma$, namely $\partial _{t} \sigma(t)$ belongs to the dual space $(H^{1})^{*}$ for a.e. $t \in (0, T)$. Furthermore, as mentioned in Remark 7.1 below, if we only have ${v} \in L^{2}(0, T;{L}^{2})$, then the convection term ${\text{div}} (\sigma {v})$ and the time derivative $\partial _{t}\sigma$ would only belong to the dual space $L^{\frac{4}{3}}(0, T;(W^{1, 4})^{*})$. However, even with the improved temporal regularity, as $\partial _{t}\sigma \notin L^{2}(0, T;(H^{1})^{*})$, we do not have a continuous embedding into the space $C^{0}([0, T];L^{2})$ and so $\sigma(0)$ may not be well-defined as an element of $L^{2}$.

We now state the two asymptotic limits of (1.1) for three dimensions, and note that analogous asymptotic limits also hold for two dimensions.

Theorem 2.3 (Limit of vanishing permeability). For $b, K \in (0, 1]$, we denote a weak solution to (1.1)-(1.2) with $\Gamma_{{v}} = 0$ and initial conditions $(\varphi_{0}, \sigma_{0})$ by $(\varphi^{K}, \mu^{K}, \sigma^{K}, {v}^{K}, p^{K})$. Then, as $b \to 0$ and $K \to 0$, it holds that

$\begin{gathered} {\varphi ^K} \to \varphi \quad {\text{ weakly - }}*\quad {\text{ in }}{L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, \frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\sigma ^K} \to \sigma \quad {\text{ weakly - }}*\quad {\text{ in }}{L^2}(0, T;{H^1}) \cap {L^\infty }(0, T;{L^2}) \cap {W^{1, \frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ {\mu ^K} \to \mu \quad {\text{ weakly }}\quad {\text{ in }}{L^2}(0, T;{H^1}), \hfill \\ {p^K} \to p\quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{H^1}) \cap {L^{\frac{8}{7}}}(0, T;{H^2}), \hfill \\ {v^K} \to 0\quad {\text{ strongly }}\quad {\text{ in }}{L^2}(0, T;{L^2}) \cap {L^{\frac{8}{7}}}(0, T;{H^1}), \hfill \\ \end{gathered}$

where $(\varphi, \mu, \sigma, p)$ satisfies

${\left\langle {{\partial _t}\varphi, \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \int_\Omega - m(\varphi )\nabla \mu \cdot \nabla \zeta + {\Gamma _\varphi }(\varphi , \mu , \sigma )\zeta {\text{ dx}},$

(2.16a)

$\int_\Omega \mu \zeta {\text{ dx}} = \int_\Omega A \Psi '(\varphi )\zeta + B\nabla \varphi \cdot \nabla \zeta - \chi \sigma \zeta {\text{ dx}},$

(2.16b)

${\left\langle {{\partial _t}\sigma, \phi } \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}} = \int_\Omega - n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla \phi - \mathcal{S}(\varphi , \mu , \sigma )\phi {\text{ dx}},$

(2.16c)

$\int_\Omega \nabla p \cdot \nabla \zeta {\text{ dx}} = \int_\Omega {(\mu + \chi \sigma )} \nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(2.16d)

for all $\zeta \in H^{1}$, $\phi \in W^{1, 5}$ and a.e. $t \in (0, T)$. A posteriori, it holds that

${\partial _t}\varphi , {\text{ }}{\partial _t}\sigma \in {L^2}(0, T;{({H^1})^*}),$

and thus $\varphi(0) = \varphi_{0}$ and $\sigma(0) = \sigma_{0}$.

Theorem 2.4 (Limit of vanishing chemotaxis). For $b, \chi \in (0, 1]$, we denote a weak solution to (1.1)-(1.2) with corresponding initial conditions $(\varphi_{0}, \sigma_{0})$ by $(\varphi^{\chi}, \mu^{\chi}, \sigma^{\chi}, {v}^{\chi}, p^{\chi})$. Then, as $b \to 0$ and $\chi \to 0$, it holds that

$\begin{gathered} {\varphi ^\chi } \to \varphi \quad {\text{ weakly - }}*\quad {\text{ in }}{L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, \frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\sigma ^\chi } \to \sigma \quad {\text{ weakly - }}*\quad {\text{ in }}{L^2}(0, T;{H^1}) \cap {L^\infty }(0, T;{L^2}) \cap {W^{1, \frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ {\mu ^\chi } \to \mu \quad {\text{ weakly }}\quad {\text{ in }}{L^2}(0, T;{H^1}), \hfill \\ {p^\chi } \to p\quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{H^1}) \cap {L^{\frac{8}{7}}}(0, T;{H^2}), \hfill \\ {v^\chi } \to v\quad {\text{ weakly }}\quad {\text{ in }}{L^2}(0, T;{L^2}) \cap {L^{\frac{8}{7}}}(0, T;{H^1}), \hfill \\ \end{gathered}$

and

$\begin{gathered} {\text{div}}({\varphi ^\chi }{v^\chi }) \to {\text{div}}(\varphi v)\quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\text{div}}({\sigma ^\chi }{v^\chi }) \to {\text{div}}(\sigma v)\quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ \end{gathered}$

where ($\varphi, \mu, \sigma, {v}, p)$ satisfies

${\left\langle {{\partial _t}\varphi, \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \int_\Omega - m(\varphi )\nabla \mu \cdot \nabla \zeta + {\Gamma _\varphi }(\varphi , \mu , \sigma )\zeta + \varphi v \cdot \nabla \zeta {\text{ dx}},$

(2.19a)

$\int_\Omega \mu \zeta {\text{ dx}} = \int_\Omega A \Psi '(\varphi )\zeta + B\nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(2.19b)

${\left\langle {{\partial _t}\sigma, \phi } \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}} = \int_\Omega - n(\varphi )D\nabla \sigma \cdot \nabla \phi - \mathcal{S}(\varphi , \mu , \sigma )\phi + \sigma {\text{ }}v \cdot \nabla \phi {\text{ dx}},$

(2.19c)

$\int_\Omega \nabla p \cdot \nabla \zeta {\text{ dx}} = \int_\Omega {\frac{1}{K}} {\Gamma _v}\zeta + \mu \nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(2.19d)

$\int_\Omega {} v \cdot {\text{ }}\zeta {\text{ dx}} = \int_\Omega - K(\nabla p - \mu \nabla \varphi ) \cdot \zeta {\text{ dx}},$

(2.19e)

for all $\zeta \in H^{1}$, $\phi \in W^{1, 5}$, ${\zeta} \in {L}^{2}$ and a.e. $t \in (0, T)$, with the attainment of initial conditions as in Definition 2.1.

---

3. Galerkin approximation

We will employ a Galerkin approximation similar to the one used in ^[22]. For the approximation, we use the eigenfunctions of the Neumann--Laplacian operator $\{w_{i}\}_{i \in \mathbb{N}}$. Recall that the inverse Neumann--Laplacian operator $\mathcal{L} := (-\Delta _{N})^{-1} \vert_{L^{2}_{0}} : L^{2}_{0} \to L^{2}_{0}$ is compact, positive and symmetric. Indeed, let $f, g \in L^{2}_{0}$ with $z = \mathcal{L}f$, $y = \mathcal{L}g$. Then,

${(\mathcal{L}f, f)_{{L^2}}} = \int_\Omega z f{\text{ dx}} = \int_\Omega {} {\left| {\nabla z} \right|^2}{\text{ dx}} \geq 0, \quad {(\mathcal{L}f, g)_{{L^2}}} = \int_\Omega \nabla z \cdot \nabla y{\text{ dx}} = {(f, \mathcal{L}g)_{{L^2}}}.$

Furthermore, let $\{f_{n}\}_{n \in \mathbb{N}} \subset L^{2}_{0}$ denote a sequence with corresponding solution sequence $\{z_{n} = \mathcal{L} f_{n} \}_{n \in \mathbb{N}} \subset H^{1} \cap L^{2}_{0}$. By elliptic regularity theory, we have that $z_{n} \in H^{2}_{N}$ for all $n \in \mathbb{N}$. Then, by reflexive weak compactness theorem and Rellich--Kondrachov theorem, there exists a subsequence such that $z_{n_{j}} \to z \in H^{1} \cap L^{2}_{0}$ as $j \to \infty$.

Thus, by the spectral theorem, the operator $\mathcal{L}$ admits a countable set of eigenfunctions $\{v_{n}\}_{n \in \mathbb{N}}$ that forms a complete orthonormal system in $L^{2}_{0}$. The eigenfunctions of the Neumann--Laplacian operator is then given by $w_{1} = 1$, $w_{i} = v_{i-1}$ for $i \geq 2$, and $\{w_{i}\}_{i \in \mathbb{N}}$ is a basis of $L^{2}$.

Elliptic regularity theory gives that $w_{i} \in H^{2}_{N}$ and for every $g \in H^{2}_{N}$, we obtain for $g_{k} := \sum_{i=1}^{k} (g, w_{i})_{L^{2}} w_{i}$ that

$\Delta {g_k} = \sum\limits_{i = 1}^k {{{(g, {w_i})}_{{L^2}}}} \Delta {w_i} = \sum\limits_{i = 1}^k {{{(g, {\lambda _i}{w_i})}_{{L^2}}}} {w_i} = \sum\limits_{i = 1}^k {{{(g, {\text{ }}\Delta {w_i})}_{{L^2}}}} {w_i} = \sum\limits_{i = 1}^k {{{(\Delta g, {w_i})}_{{L^2}}}} {w_i},$

where $\lambda_{i}$ is the corresponding eigenvalue to $w_{i}$. This shows that $\Delta g_{k}$ converges strongly to $\Delta g$ in $L^{2}$. Making use of elliptic regularity theory again gives that $g_{k}$ converges strongly to $g$ in $H^{2}_{N}$. Thus the eigenfunction $\{w_{i}\}_{i \in \mathbb{N}}$ of the Neumann--Laplace operator forms an orthonormal basis of $L^{2}$ and is also a basis of $H^{2}_{N}$.

Later in Section 5, we will need to use the property that $H^{2}_{N}$ is dense in $H^{1}$ and $W^{1, 5}$. We now sketch the argument for the denseness of $H^{2}_{N}$ in $W^{1, 5}$ and the argument for $H^{1}$ follows in a similar fashion.

Lemma 3.1. $H^{2}_{N}$ is dense in $W^{1, 5}$.

Proof. Take $g \in W^{1, 5}$, as $\Omega$ has a $C^{3}$-boundary, by standard results [12, Thm. 3, x5.3.3] there exists a sequence $g_{n} \in C^{\infty}(\overline{\Omega})$ such that $g_{n} \to g$ strongly in $W^{1, 5}$. Let $\varepsilon > 0$ be fixed, and define $D_{\varepsilon } := \{ x \in \Omega : \text{ dist}(x, \partial \Omega) \leq \varepsilon \}$. Let $\zeta_{\varepsilon } \in C^{\infty}_{c}(\Omega)$ be a smooth cut-off function such that $\zeta_{\varepsilon } = 1$ in $\Omega \setminus \overline{D_{\varepsilon }}$ and $\zeta_{\varepsilon } = 0$ in $\overline{D_{\frac{\varepsilon }{2}}}$.

As $g_{n} \in C^{\infty}(\overline{\Omega})$, its trace on $\partial \Omega$ is well-defined. Choosing $\varepsilon$ sufficiently small allows us to use a classical result from differential geometry about tubular neighbourhoods, i.e., for any $z \in {\text{Tu}}{{\text{b}}_\varepsilon }(\partial \Omega ): = \{ x \in {\mathbb{R}^d}:{\text{ }}\left| {{\text{dist}}(z, \partial \Omega )} \right| \leq \varepsilon \}$ there exists a unique $y \in \partial \Omega$ such that

$z = y + {\text{dist}}(z, \partial \Omega )n(y),$

where ${n}$ is the outer unit normal of $\partial \Omega$. We consider a bounded smooth function $f_{n, \varepsilon }: \mathbb{R}^{d} \to \mathbb{R}$ such that

${f_{n, \varepsilon }}(z) = {g_n}(y){\text{ for all }}z \in {\text{Tu}}{{\text{b}}_\varepsilon }(\partial \Omega ){\text{ satisfying }}z = y + {\text{ dist}}(z, \partial \Omega )n(y).$

We now define the smooth function $G_{n, \varepsilon }$ as

$G_{n, \varepsilon }(x) := \zeta_{\varepsilon }(x) g_{n}(x) + (1-\zeta_{\varepsilon }(x)) f_{n, \varepsilon }(x).$

By construction, the values of the function $f_{n, \varepsilon }$ in $D_{\varepsilon } \subset \mathrm{Tub}_{\varepsilon }(\partial \Omega)$ are constant in the normal direction, so $\nabla G_{n, \varepsilon } \cdot {n} = 0$ on $\partial \Omega$ and thus $G_{n, \varepsilon } \in H^{2}_{N}$. Furthermore, we compute that

$\begin{gathered} {\left\| {{G_{n, \varepsilon }} - {g_n}} \right\|_{{L^5}}} = {\left\| {(1 - {\zeta _\varepsilon })({f_{n, \varepsilon }} - {g_n})} \right\|_{{L^5}({D_\varepsilon })}}, \hfill \\ {\left\| {\nabla ({G_{n, \varepsilon }} - {g_n})} \right\|_{{L^5}}} = {\left\| {({g_n} - {f_{n, \varepsilon }})\nabla {\zeta _\varepsilon } + (1 - {\zeta _\varepsilon })\nabla {g_n} + (1 - {\zeta _\varepsilon })\nabla {f_{n, \varepsilon }}} \right\|_{{L^5}}}. \hfill \\ \end{gathered}$

Using that $g_{n}, f_{n, \varepsilon }$ are smooth functions on $\overline{\Omega}$ and that the Lebesgue measure of $D_{\varepsilon }$ tends to zero as $\varepsilon \to 0$ we have the strong convergence of $G_{n, \varepsilon }$ to $g_{n}$ in $L^{5}$. For the difference in the gradients, we use that $\zeta_{\varepsilon } \to 1$ a.e. in $\Omega$, Lebesgue's dominated convergence theorem and the boundedness of $\nabla g_{n}$ and $\nabla f_{n, \varepsilon }$ to deduce that

${\left\| {(1 - {\zeta _\varepsilon })\nabla {g_n}} \right\|_{{L^5}}} + {\left\| {(1 - {\zeta _\varepsilon })\nabla {f_{n, \varepsilon }}} \right\|_{{L^5}}} \to 0{\text{ as }}\varepsilon \to 0.$

For the remaining term ${\left\| {({g_n} - {f_{n, \varepsilon }})\nabla {\zeta _\varepsilon }} \right\|_{{L^5}}}$ we use that the support of $\nabla \zeta_{\varepsilon }$ lies in $D_{\varepsilon } \setminus \overline{D_{\frac{\varepsilon }{2}}}$ and for any $z \in D_{\varepsilon } \setminus \overline{D_{\frac{\varepsilon }{2}}}$,

$\begin{gathered} \left| {{f_{n, \varepsilon }}(z) - {g_n}(z)} \right| = \left| {{g_n}(y) - {g_n}(y + {\text{dist}}(z, \partial \Omega )n(y))} \right| \hfill \\ \leq \int_0^{{\text{dist}}(z, \partial \Omega )} {} \left| {\nabla {g_n}(y + \xi n(y))} \right|d\xi \leq {\left\| {\nabla {g_n}} \right\|_{{L^\infty }}}dist({\text{z}}, \partial \Omega ) \leq {\text{C}}\varepsilon . \hfill \\ \end{gathered}$

That is, $f_{n, \varepsilon }$ converges uniformly to $g_{n}$ in $D_{\varepsilon } \setminus \overline{D_{\frac{\varepsilon }{2}}}$. Furthermore, using ${\left\| {\nabla {\zeta _\varepsilon }} \right\|_{{L^\infty }}} \leq \frac{C}{\varepsilon }$ in $D_{\varepsilon } \setminus \overline{D_{\frac{\varepsilon }{2}}}$ and $\left| {{D_\varepsilon } \setminus \overline {{D_{\frac{\varepsilon }{2}}}} } \right| \leq C\varepsilon$ we obtain ${\left\| {({g_n} - {f_{\varepsilon , n}})\nabla {\zeta _\varepsilon }} \right\|_{{L^5}}} \leq C{\varepsilon ^{\frac{1}{5}}} \to 0$ as $\varepsilon \to 0$. This shows that $G_{n, \varepsilon }$ converges strongly to $g_{n}$ in $W^{1, 5}$.

We denote

$W_{k} := \mathrm{span} \{ w_{1}, \dots, w_{k} \}$

as the finite dimensional space spanned by the first k basis functions and consider

$\varphi_{k}(t, x) =\sum_{i=1}^{k} \alpha_{i}^{k}(t) w_{i}(x), \; \mu_{k}(t, x) = \sum_{i=1}^{k} \beta_{i}^{k}(t) w_{i}(x), \; \sigma_{k}(t, x) = \sum_{i=1}^{k} \gamma_{i}^{k}(t) w_{i}(x),$

(3.1a)

and the following Galerkin ansatz: For $1 \leq j \leq k$,

$\int_\Omega {{\partial _t}} {\varphi _k}{w_j}{\text{ dx}} = \int_\Omega - m({\varphi _k})\nabla {\mu _k} \cdot \nabla {w_j} + {\Gamma _{\varphi , k}}{w_j} + {\varphi _k}{v_k} \cdot \nabla {w_j}dx,$

(3.2a)

$\int_\Omega {{\mu _k}} {w_j}{\text{ dx}} = \int_\Omega A \Psi '({\varphi _k}){w_j} + B\nabla {\varphi _k} \cdot \nabla {w_j} - \chi {\sigma _k}{w_j}{\text{ dx}},$

(3.2b)

$\begin{gathered} \int_\Omega {{\partial _t}} {\sigma _k}{w_j}dx = \int_\Omega - n({\varphi _k})(D\nabla {\sigma _k} - \chi \nabla {\varphi _k}) \cdot \nabla {w_j} - {\mathcal{S}_k}{w_j} + {\sigma _k}{v_k} \cdot \nabla {w_j}dx \\ + \int_{\partial \Omega } b ({\sigma _\infty } - {\sigma _k}){w_j}{\text{ d}}{H^{d - 1}}, \\ \end{gathered}$

(3.2c)

where we define the Galerkin ansatz for the pressure $p_{k}$ and the velocity field ${v}_{k}$ by

${p_k} = {( - {\Delta _N})^{ - 1}}\left( {\frac{1}{K}{\Gamma _v} - {\text{div}} (({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k})} \right),$

(3.3)

${v_k} = - K(\nabla {p_k} - ({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k}),$

(3.4)

and we set

$\Gamma_{\varphi, k} := \Gamma_{\varphi}(\varphi_{k}, \mu_{k}, \sigma_{k}), \quad \mathcal{S}_{k} := \mathcal{S}(\varphi_{k}, \mu_{k}, \sigma_{k}).$

Note that in (3.3), the properties $\Gamma_{{v}} \in L^{2}_{0}$ and $\nabla \varphi_{k} \cdot {n} = 0$ on $\partial \Omega$ show that the term inside the bracket belongs to $L^{2}_{0}$ and hence $p_{k}$ is well-defined. Let ${M}$ and ${S}$ denote the following mass and stiffness matrices, respectively: For $1 \leq i, j \leq k$,

${M_{ij}} = \int_\Omega {{w_i}} {w_j}{\text{ dx}}, \quad {\text{ }}{S_{ij}}: = \int_\Omega \nabla {w_i} \cdot \nabla {w_j}{\text{ dx}}.$

Thanks to the orthonormality of $\{w_{i}\}_{i \in \mathbb{N}}$ in $L^{2}$, we see that ${M}$ is the identity matrix. It is convenient to define the following matrices with components

$\begin{gathered} {({C^k})_{ji}}: = \int_\Omega {{w_i}} {v_k} \cdot \nabla {w_j}dx, \quad {({M_{\partial \Omega }})_{ji}}: = \int_{\partial \Omega } {{w_i}} {w_j}dH, \hfill \\ {(S_m^k)_{ji}}: = \int_\Omega m ({\varphi _k})\nabla {w_i} \cdot \nabla {w_j}dx, \quad {(S_n^k)_{ji}}: = \int_\Omega n ({\varphi _k})\nabla {w_i} \cdot \nabla {w_j}dx, \hfill \\ \end{gathered}$

for $1 \leq i, j \leq k$. Furthermore, we introduce the notation

$R_{\varphi , j}^k: = \int_\Omega {{\Gamma _{\varphi , k}}} {w_j}{\text{ dx}}, \quad R_{S, j}^k: = \int_\Omega {{\mathcal{S}_k}} {w_j}{\text{ dx}}, \quad \psi _j^k: = \int_\Omega {\Psi '} ({\varphi _k}){w_j}{\text{ dx}}, \quad \Sigma _j^k: = \int_{\partial \Omega } {{\sigma _\infty }} {w_j}{\text{ dH}},$

for $1 \leq i, j \leq k$, and denote

$R_\varphi ^k: = {(R_{\varphi , 1}^k, \ldots , R_{\varphi , k}^k)^ \top }, \quad {\text{ }}R_S^k: = {(R_{S, 1}^k, \ldots , R_{S, k}^k)^ \top }, \quad {\text{ }}{\psi ^k}: = {(\psi _1^k, \ldots , \psi _k^k)^ \top }, \quad {\text{ }}{\Sigma ^k}: = {(\Sigma _1^k, \ldots , \Sigma _k^k)^ \top },$

as the corresponding vectors, so that we obtain an initial value problem for a system of equations for ${\alpha _k}: = {(\alpha _1^k, \ldots \alpha _k^k)^ \top }, {\beta _k}: = {(\beta _1^k, \ldots , \beta _k^k)^ \top }$, and ${\gamma _k}: = {(\gamma _1^k, \ldots , \gamma _k^k)^ \top }$ as follows,

$\frac{{\text{d}}}{{{\text{dt}}}}{\alpha _k} = - S_m^k{\beta _k} + R_\varphi ^k + {C^k}{\alpha _k},$

(3.5a)

${\beta _k} = A{\psi ^k} + BS{\alpha _k} - \chi {\gamma _k},$

(3.5b)

$\frac{{\text{d}}}{{{\text{dt}}}}{\gamma _k} = - S_n^k(D{\gamma _k} - \chi {\alpha _k}) - R_S^k + {C^k}{\gamma _k} - b{M_{\partial \Omega }}{\gamma _k} + b{\Sigma ^k},$

(3.5c)

${p_k} = {\left( { - {\Delta _N}} \right)^{ - 1}}\left( {\frac{1}{K}{\Gamma _v} - {\text{div}} (({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k})} \right)$

(3.5d)

${v_k} = - K(\nabla {p_k} - ({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k}),$

(3.5e)

and we supplement (3.5) with the initial conditions

${({\alpha _k})_j}(0) = \int_\Omega {{\varphi _0}} {w_j}{\text{ dx}}, \quad {({\gamma _k})_j}(0) = \int_\Omega {{\sigma _0}} {w_j}{\text{ dx}},$

(3.6)

for $1 \leq j \leq k$, which satisfy

${\left\| {\sum\limits_{i = 1}^k {{{({\alpha _k})}_i}} (0){w_i}} \right\|_{{H^1}}} \leq C{\left\| {{\varphi _0}} \right\|_{{H^1}}}, \quad {\text{ }}{\left\| {\sum\limits_{j = 1}^k {{{({\gamma _k})}_i}} (0){w_i}} \right\|_{{L^2}}} \leq {\left\| {{\sigma _0}} \right\|_{{L^2}}}\quad \forall k \in \mathbb{N},$

(3.7)

for some constant C not depending on k.

We can substitute (3.5b), (3.5d) and (3.5e) into (3.5a) and (3.5c), and obtain a coupled system of ordinary differential equations for ${\alpha}_{k}$ and ${\gamma}_{k}$, where ${S}_{m}^{k}$, ${C}^{k}$ and ${S}_{n}^{k}$ depend on the solutions ${\alpha}_{k}$ and ${\gamma}_{k}$ in a non-linear manner. Continuity of $m(\cdot)$, $n(\cdot)$, $\Psi'(\cdot)$ and the source terms, and the stability of $(-\Delta _{N})^{-1}$ under perturbations imply that the right-hand sides of (3.5) depend continuously on $({\alpha}_{k}, {\gamma}_{k})$. Thus, we can appeal to the theory of ODEs (via the Cauchy--Peano theorem [4, Chap. 1, Thm. 1.2]) to infer that the initial value problem (3.5)-(3.6) has at least one local solution pair $({\alpha}_{k}, {\gamma}_{k})$ defined on $[0, t_{k}]$ for each $k \in \mathbb{N}$.

We may define ${\beta}_{k}$ via the relation (3.5b) and hence the Galerkin ansatz $\varphi_{k}, \mu_{k}$ and $\sigma_{k}$ can be constructed from (3.1). Then, we can define $p_{k}$ and ${v}_{k}$ via (3.3) and (3.4), respectively. Furthermore, as the basis function $w_{j}$ belongs to $H^{2}$ for each $j \in \mathbb{N}$, by the Sobolev embedding $H^{2} \subset L^{\infty}$, we obtain that ${\text{div}} (w_{i} \nabla w_{j}) \in L^{2}$ for $i, j \in \mathbb{N}$ and hence the function ${\text{div}} ((\mu_{k} + \chi \sigma_{k}) \nabla \varphi_{k})$ belongs to $L^{2}$. Then, by elliptic regularity theory, we find that $p_{k}(t) \in H^{2}_{N} \cap L^{2}_{0}$ for all $t \in [0, t_{k}]$. This in turn implies that

${v_k}(t) \in \{ f \in {\text{ }}{H^1}:{\text{div}} f = {\Gamma _v}, \;{\text{ }}f \cdot n = 0{\text{ on }}\partial \Omega \} {\text{ for all }}t \in [0, {t_k}].$

(3.8)

Next, we show that the Galerkin ansatz can be extended to the interval [0, T] using a priori estimates.

---

4. A priori estimates

In this section, the positive constants $C$ are independent of $k$, $\Gamma_{{v}}$, $K$, $b$ and $\chi$, and may change from line to line. We will denote positive constants that are uniformly bounded for $b, \chi \in (0, 1]$ and are also uniformly bounded for $K \in (0, 1]$ when $\Gamma_{{v}} = 0$ by the symbol $\mathcal{E}$.

We first state the energy identity satisfied by the Galerkin ansatz. Let $\delta_{ij}$ denote the Kronecker delta. Multiplying (3.2a) with $\beta_{j}^{k}$, (3.2b) with $\frac{{\text{d}}}{{{\text{dt}}}}\alpha _j^k$, (3.2c) with $N_{, \sigma}^{k} := D \gamma_{j}^{k} + \chi (\delta_{1j} - \alpha_{j}^{k})$, and then summing the product from $j = 1$ to $k$ lead to

$\begin{gathered} \int_\Omega {{\partial _t}} {\varphi _k}{\mu _k}{\text{dx}} = \int_\Omega - m({\varphi _k}){\left| {\nabla {\mu _k}} \right|^2} + {\Gamma _{\varphi , k}}{\mu _k} + {\varphi _k}{v_k} \cdot \nabla {\mu _k}{\text{dx}}, \hfill \\ \int_\Omega {{\mu _k}} {\partial _t}{\varphi _k}{\text{dx}} = \frac{{\text{d}}}{{{\text{dt}}}}\int_\Omega A \Psi '({\varphi _k}) + \frac{B}{2}{\left| {\nabla {\varphi _k}} \right|^2}{\text{dx}} - \int_\Omega \chi {\sigma _k}{\partial _t}{\varphi _k}{\text{dx}}, \hfill \\ \int_\Omega {{\partial _t}} {\sigma _k}N_{, \sigma }^k{\text{dx}} = \int_\Omega - n({\varphi _k}){\left| {\nabla N_{, \sigma }^k} \right|^2} - {\mathcal{S}_k}N_{, \sigma }^k + {\sigma _k}{v_k} \cdot \nabla N_{, \sigma }^k{\text{dx}} \hfill \\ {\text{ }} + \int_{\partial \Omega } b ({\sigma _\infty } - {\sigma _k})N_{, \sigma }^k{\text{d}}{H^{d - 1}}. \hfill \\ \end{gathered}$

Here, we used that $w_{1} = 1$ and $\nabla w_{1} = {0}$. Then, summing the three equations leads to

$\begin{gathered} \frac{{\text{d}}}{{{\text{dt}}}}\int_\Omega A \Psi ({\varphi _k}) + \frac{B}{2}{\left| {\nabla {\varphi _k}} \right|^2} + N({\varphi _k}, {\sigma _k}){\text{dx}} \hfill \\ \quad + \int_\Omega m ({\varphi _k}){\left| {\nabla {\mu _k}} \right|^2} + n({\varphi _k}){\left| {\nabla N_{, \sigma }^k} \right|^2}{\text{dx}} + \int_{\partial \Omega } D b{\left| {{\sigma _k}} \right|^2}{\text{d}}{H^{d - 1}} \hfill \\ \quad = \int_\Omega {{\Gamma _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}N_{, \sigma }^k + ({\varphi _k}{v_k} \cdot \nabla {\mu _k} + {\sigma _k}{v_k} \cdot \nabla N_{, \sigma }^k)dx \hfill \\ \quad + \int_{\partial \Omega } b ({\sigma _\infty }N_{, \sigma }^k - {\sigma _k}\chi (1 - {\varphi _k})){\text{d}}{H^{d - 1}}. \hfill \\ \end{gathered}$

(4.2)

Next, multiplying (3.4) with $\frac{1}{K} {v}_{k}$, integrating over $\Omega$ and integrating by parts gives

$\begin{gathered} \int_\Omega {\frac{1}{K}} {\left| {{v_k}} \right|^2}{\text{dx}} = \int_\Omega - \nabla {p_k} \cdot {v_k} + ({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k} \cdot {v_k}{\text{dx}} \\ = \int_\Omega {{\Gamma _v}} {p_k} + ({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k} \cdot {v_k}{\text{ dx}}, \\ \end{gathered}$

where we used that ${\text{div}} {v}_{k} = \Gamma_{{v}}$ and ${v}_{k} \cdot {n} = 0$ on $\partial \Omega$. Similarly, we see that

$\begin{gathered} \int_\Omega {({\varphi _k}\nabla {\mu _k} + {\sigma _k}\nabla N_{, \sigma }^k)} \cdot {v_k}{\text{dx}} = \int_\Omega {{\varphi _k}} {v_k} \cdot \nabla {\mu _k} + {\sigma _k}{v_k} \cdot \nabla (D{\sigma _k} + \chi (1 - {\varphi _k})){\text{dx}} \\ = - \int_\Omega {{\varphi _k}} {\Gamma _v}{\mu _k} + ({\mu _k} + \chi {\sigma _k}){v_k} \cdot \nabla {\varphi _k} - \tfrac{D}{2}{v_k} \cdot \nabla {\left| {{\sigma _k}} \right|^2}{\text{dx}} \\ = - \int_\Omega {{\Gamma _v}} \left( {{\varphi _k}{\mu _k} + \tfrac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right) + ({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k} \cdot {v_k}{\text{ dx}}. \\ \end{gathered}$

In particular, we have

$\int_\Omega {\frac{1}{K}} {\left| {{v_k}} \right|^2}{\text{ dx}} = \int_\Omega {{\Gamma _v}} \left( {{p_k} - {\mu _k}{\varphi _k} - \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right) - ({\varphi _k}\nabla {\mu _k} + {\sigma _k}\nabla N_{, \sigma }^k) \cdot {v_k}{\text{ dx}}.$

(4.4)

Adding (4.3) to (4.2) leads to

$\begin{gathered} \frac{{\text{d}}}{{{\text{dt}}}}\int_\Omega A \Psi ({\varphi _k}) + \frac{B}{2}{\left| {\nabla {\varphi _k}} \right|^2} + N({\varphi _k}, {\sigma _k}){\text{dx}} \hfill \\ \quad + \int_\Omega m ({\varphi _k}){\left| {\nabla {\mu _k}} \right|^2} + n({\varphi _k}){\left| {\nabla N_{, \sigma }^k} \right|^2} + \frac{1}{K}{\left| {{v_k}} \right|^2}{\text{dx}} + \int_{\partial \Omega } D b{\left| {{\sigma _k}} \right|^2}{\text{d}}{H^{d - 1}} \hfill \\ \quad = \int_\Omega {{\Gamma _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}N_{, \sigma }^k + {\Gamma _v}\left( {{p_k} - {\mu _k}{\varphi _k} - \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right){\text{dx}} \hfill \\ \quad + \int_{\partial \Omega } b ({\sigma _\infty }(D{\sigma _k} + \chi (1 - {\varphi _k})) - {\sigma _k}\chi (1 - {\varphi _k})){\text{d}}{H^{d - 1}}. \hfill \\ \end{gathered}$

(4.3)

To derive the first a priori estimate for the Galerkin ansatz, it suffices to bring (4.4) into a form where we can apply Gronwall's inequality. We start with estimating the boundary term on the right-hand side of (4.4). By Hölder's inequality and Young's inequality,

$\begin{gathered} \left| {\int_{\partial \Omega } b ({\sigma _\infty }(D{\sigma _k} + \chi (1 - {\varphi _k})) - {\sigma _k}\chi (1 - {\varphi _k})){\text{d}}{H^{d - 1}}} \right| \hfill \\ \quad \leq b\left( {{{\left\| {{\sigma _\infty }} \right\|}_{{L^2}(\partial \Omega )}}{{\left\| {D{\sigma _k} + \chi (1 - {\varphi _k})} \right\|}_{{L^2}(\partial \Omega )}} + \chi {{\left\| {{\sigma _k}} \right\|}_{{L^2}(\partial \Omega )}}\left( {{{\left| {\partial \Omega } \right|}^{\frac{1}{2}}} + {{\left\| {{\varphi _k}} \right\|}_{{L^2}(\partial \Omega )}}} \right)} \right) \hfill \\ \quad \leq \frac{{Db}}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2 + b\left( {1 + \frac{{{\chi ^2}}}{D}} \right)\left\| {{\varphi _k}} \right\|_{{L^2}(\partial \Omega )}^2 + bC\left( {\chi + (1 + {\chi ^2})\left\| {{\sigma _\infty }} \right\|_{{L^2}(\partial \Omega )}^2} \right). \hfill \\ \end{gathered}$

By the trace theorem and the growth condition (2.4), we have

$\begin{gathered} \left\| \varphi \right\|_{{L^2}(\partial \Omega )}^2 \leq C_{{\text{tr}}}^2\left( {\left\| \varphi \right\|_{{L^2}}^2 + \left\| {\nabla \varphi } \right\|_{{L^2}}^2} \right) \\ \leq C_{{\text{tr}}}^2\left( {\frac{1}{{{R_1}}}{{\left\| {\Psi (\varphi )} \right\|}_{{L^1}}} + \left\| {\nabla \varphi } \right\|_{{L^2}}^2} \right) + C({R_2}, \left| \Omega \right|, {C_{{\text{tr}}}}), \\ \end{gathered}$

(4.5)

where the positive constant $C_{\mathrm{tr}}$ from the trace theorem only depends on $\Omega$, and so

$\begin{gathered} \int_{\partial \Omega } b ({\sigma _\infty }(D{\sigma _k} + \chi (1 - {\varphi _k})) - {\sigma _k}\chi (1 - {\varphi _k})){\text{d}}{H^{d - 1}} \hfill \\ \quad \leq \frac{{Db}}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2 + Cb\left( {1 + {\chi ^2}} \right)\left( {{{\left\| {\Psi ({\varphi _k})} \right\|}_{{L^1}}} + \left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2} \right) \hfill \\ \quad + Cb\left( {1 + {\chi ^2}} \right) + bC\left( {\chi + (1 + {\chi ^2})\left\| {{\sigma _\infty }} \right\|_{{L^2}(\partial \Omega )}^2} \right). \hfill \\ \end{gathered}$

(4.6)

---

4.1. Estimation of the source terms

For the source term

$\int_\Omega {{\Gamma _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}N_{, \sigma }^k + {\Gamma _v}\left( {{p_k} - {\mu _k}{\varphi _k} - \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right){\text{ dx}}$

that appears on the right-hand side of (4.4) we will divide its analysis into two parts. We first analyse the part involving $\Gamma_{{v}}$, which will involve a closer look at the Darcy subsystem to deduce an estimate on ${\left\| {{p_k}} \right\|_{{L^2}}}$. For the remainder $\Gamma_{\varphi, k} \mu_{k} - \mathcal{S}_{k} N_{, \sigma}^{k}$ term we will estimate it differently based on the assumptions on $\Theta_{\varphi}$.

---

4.1.1 Pressure estimates

Before we estimate the source terms involving $\Gamma_{{v}}$, we look at the Darcy subsystem, which can be expressed as an elliptic equation for the pressure (we will drop the subscript $k$ for clarity)

$- \Delta p = \frac{1}{K}{\Gamma _v} - {\text{div}} ((\mu + \chi \sigma )\nabla \varphi ){\text{ in }}\Omega , {\text{ with }}\bar p = 0,$

(4.7a)

$\partial p = 0{\text{ on }}\partial \Omega .$

(4.7b)

The following lemma is similar to [22, Lem. 3.1], and the hypothesis is fulfilled by the Galerkin ansatz.

Lemma 4.1. Let $\Omega \subset \mathbb{R}^{3}$ be a bounded domain with $C^{3}$-boundary. Given $\varphi \in H^{2}_{N}$, $\mu, \sigma \in H^{1}$, the source term $\Gamma_{{v}} \in L^{2}_{0}$, and the function $p$ satisfying the above elliptic equation (4.7). Then, the following estimate hold

${\left\| p \right\|_{{L^2}}} \leq \frac{C}{K}{\left\| {{\Gamma _v}} \right\|_{{L^2}}} + C\left( {{{\left\| {\nabla \mu } \right\|}_{{L^2}}} + \chi {{\left\| \sigma \right\|}_{{L^6}}}} \right){\left\| {\nabla \varphi } \right\|_{{L^{\frac{3}{2}}}}} + C\bar \mu {\left\| {\nabla \varphi } \right\|_{{L^2}}},$

(4.8)

for some positive constant $C$ depending only on $\Omega$.

Proof. We first recall some properties of the inverse Neumann-Laplacian operator. Suppose for $g \in L^{2}_{0}$, $f = (-\Delta _{N})^{-1} g \in H^{1} \cap L^{2}_{0}$ solves

$- \Delta f = g{\text{ in }}\Omega , \quad {\text{ }}\partial f = 0{\text{ on }}\partial \Omega .$

(4.9)

Then, testing with $f$ and integrating over $\Omega$, applying integration by parts and the Poincaré inequality (1.9) leads to

${\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}g} \right\|_{{H^1}}} = {\left\| f \right\|_{{H^1}}} \leq c{\left\| {\nabla f} \right\|_{{L^2}}} \leq C{\left\| g \right\|_{{L^2}}},$

(4.10)

for positive constants $c$ and $C$ depending only on $C_{p}$. Elliptic regularity theory then gives that $f \in H^{2}_{N}$ with

${\left\| f \right\|_{{H^2}}} \leq C\left( {{{\left\| f \right\|}_{{H^1}}} + {{\left\| g \right\|}_{{L^2}}}} \right) \leq C{\left\| g \right\|_{{L^2}}},$

(4.11)

with a positive constant $C$ depending only on $\Omega$. Returning to the pressure system, we observe from (2.10) and the above that

$\begin{gathered} {\left\| p \right\|_{{L^2}}} \leq \frac{1}{K}{\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}{\Gamma _v}} \right\|_{{L^2}}} + {\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}\left( { {\text{div}} ((\mu + \chi \sigma )\nabla \varphi )} \right)} \right\|_{{L^2}}} \\ \leq \frac{C}{K}{\left\| {{\Gamma _v}} \right\|_{{L^2}}} + {\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}\left( {div((\mu - \bar \mu + \chi \sigma )\nabla \varphi )} \right)} \right\|_{{L^2}}} \\ + {\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}(div(\bar \mu {\mkern 1mu} \nabla \varphi ))} \right\|_{{L^2}}}, \\ \end{gathered}$

(4.12)

for some positive constant $C$ depending only on $C_{p}$. Note that the third term on the right-hand side can be estimated as

$\bar \mu {\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}{\text{div}}\nabla (\varphi - \bar \varphi )} \right\|_{{L^2}}} = \bar \mu {\left\| {\varphi - \bar \varphi } \right\|_{{L^2}}} \leq {C_p}\bar \mu {\left\| {\nabla \varphi } \right\|_{{L^2}}}.$

(4.13)

We now consider estimating the second term on the right-hand side of (4.12). By assumption $\mu, \sigma \in H^{1}$ and $\varphi \in H^{2}_{N}$, we have that

${\left\| {(\mu - \bar \mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}} \leq {\left\| {\mu - \bar \mu + \chi \sigma } \right\|_{{L^6}}}{\left\| {\nabla \varphi } \right\|_{{L^3}}},$

(4.14)

and so if we consider the function $h: = {( - {\Delta _N})^{ - 1}}({\text{div}}((\mu - \bar \mu + \chi \sigma )\nabla \varphi ))$, then we obtain that

$\int_\Omega \nabla h \cdot \nabla \zeta {\text{ dx}} = \int_\Omega - (\mu - \bar \mu + \chi \sigma )\nabla \varphi \cdot \nabla \zeta {\text{ dx}}\quad \forall \zeta \in {H^1}$

(4.15)

must hold, and by (4.14) and the Poincaré inequality (1.9) with zero mean it holds that $h \in H^{1} \cap L^{2}_{0}$. We now define $f:= (-\Delta _{N})^{-1}(h) \in H^{2}_{N}$, and consider testing with $\zeta = f$ in (4.15), leading to

$\int_\Omega {} {\left| h \right|^2}{\text{ dx}} = \int_\Omega \nabla h \cdot \nabla f{\text{ dx}} = \int_\Omega - (\mu - \bar \mu + \chi \sigma )\nabla \varphi \cdot \nabla f{\text{ dx}}.$

Since $f \in H^{2}_{N}$, elliptic regularity theory and Hölder's inequality gives

$\begin{gathered} \left\| h \right\|_{{L^2}}^2 \leq {\left\| {(\mu - \bar \mu + \chi \sigma )\nabla \varphi } \right\|_{{L^{\frac{6}{5}}}}}{\left\| {\nabla f} \right\|_{{L^6}}} \leq C{\left\| {(\mu - \bar \mu + \chi \sigma )\nabla \varphi } \right\|_{{L^{\frac{6}{5}}}}}{\left\| f \right\|_{{H^2}}} \\ \leq C{\left\| {(\mu - \bar \mu + \chi \sigma )\nabla \varphi } \right\|_{{L^{\frac{6}{5}}}}}{\left\| h \right\|_{{L^2}}}, \\ \end{gathered}$

where the constant $C$ depends on $\Omega$ and the constant in (4.11). Thus we obtain

$\begin{gathered} {\left\| {{{\left( { - {\Delta _N}} \right)}^{ - 1}}\left( {{\text{div}}((\mu - \bar \mu + \chi \sigma )\nabla \varphi )} \right)} \right\|_{{L^2}}} \leq C{\left\| {(\mu - \bar \mu + \chi \sigma )\nabla \varphi } \right\|_{{L^{\frac{6}{5}}}}} \\ \leq C\left( {{{\left\| {\mu - \bar \mu } \right\|}_{{L^6}}} + \chi {{\left\| \sigma \right\|}_{{L^6}}}} \right){\left\| {\nabla \varphi } \right\|_{{L^{\frac{3}{2}}}}} \\ \end{gathered}$

(4.16)

for some constant $C$ depending only on $\Omega$. By the Sobolev embedding $H^{1} \subset L^{6}$ (with constant $C_{\mathrm{Sob}}$ that depends only on $\Omega$) and the Poincaré inequality, we find that

${\left\| {\mu - \bar \mu } \right\|_{{L^6}}} \leq {C_{{\text{Sob}}}}{\left\| {\mu - \bar \mu } \right\|_{{H^1}}} \leq c({C_{{\text{Sob}}}}, {C_p}){\left\| {\nabla \mu } \right\|_{{L^2}}}.$

(4.17)

Substituting the above elements into (4.12) yields (4.8).

Remark 4.1. We choose not to use the estimate

$c{\left\| h \right\|_{{L^2}}} \leq {\left\| {\nabla h} \right\|_{{L^2}}} \leq {\left\| {(\mu - \bar \mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}}$

(4.18)

obtained from substituting $\zeta = h$ in (4.15), where $c$ is a positive constant depending only on $C_{p}$, since by (4.14) we require control of $\nabla \varphi$ in the ${L}^{3}(\Omega)$-norm and this is not available when deriving the first a priori estimate. Thus, we make use of the auxiliary problem $f = (-\Delta _{N})^{-1}(h)$ to derive another estimate on ${\left\| h \right\|_{{L^2}}}$ that involves controlling $\nabla \varphi$ in the weaker ${L}^{\frac{3}{2}}(\Omega)$-norm.

Next, we state regularity estimates for the pressure and the velocity field. The hypothesis will be fulfilled for the Galerkin ansatz once we derived the a priori estimates in Section 4. Note that in Lemma 4.2 below, we consider a source term $\Gamma_{{v}} \in L^{2}(0, T;L^{2}_{0})$, so that our new regularity results for the pressure and the velocity is also applicable to the setting considered in ^[22].

Lemma 4.2. Let $\varphi \in L^{\infty}(0, T;H^{1}) \cap L^{2}(0, T;H^{2}_{N} \cap H^{3})$, $\sigma \in L^{2}(0, T;H^{1})$, $\mu \in L^{2}(0, T;H^{1})$, the source term $\Gamma_{{v}} \in L^{2}(0, T;L^{2}_{0})$, and the function $p$ satisfying (4.7). Then,

$\left\| p \right\|_{{L^{\frac{8}{5}}}({H^1})}^{\frac{8}{5}} \leq {C_1}\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{6}{5}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{5}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{2}{5}} + \frac{{{C_1}}}{K}{T^{\frac{1}{5}}}\left\| {{\Gamma _v}} \right\|_{{L^2}({L^2})}^{\frac{8}{5}},$

(4.19)

for some positive constant $C_{1}$ depending only on $\Omega$, and

$\begin{gathered} \left\| p \right\|_{{L^{\frac{8}{7}}}({H^2})}^{\frac{8}{7}} \leq {C_2}{T^{\frac{3}{7}}}{K^{ - \frac{8}{7}}}\left\| {{\Gamma _v}} \right\|_{{L^2}({L^2})}^{\frac{8}{7}} + {C_2}{T^{\frac{2}{7}}}\left\| p \right\|_{{L^{\frac{8}{5}}}({H^1})}^{\frac{8}{7}} \\ + {C_2}\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{2}{7}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{7}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{6}{7}}, \\ \end{gathered}$

(4.20)

for some positive constant $C_{2}$ depending only on $\Omega$. Moreover, if we have the relation

$v = - K\left( {\nabla p - (\mu + \chi \sigma )\nabla \varphi } \right),$

then

$\left\| {{\text{D}}v} \right\|_{{L^{\frac{8}{7}}}({L^2})}^{\frac{8}{7}} \leq {C_3}K\left\| p \right\|_{{L^{\frac{8}{7}}}({H^2})}^{\frac{8}{7}} + {C_3}K\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{7}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{6}{7}}\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{2}{7}},$

(4.21)

for some positive constant $C_{3}$ depending only on $\Omega$.

Proof. From (4.7) we see that $p$ satisfies $\bar p{\text{ = 0}}$ and

$\int_\Omega \nabla p \cdot \nabla \zeta {\text{ dx}} = \int_\Omega {(\mu + \chi \sigma )} \nabla \varphi \cdot \nabla \zeta + \frac{1}{K}{\Gamma _v}\zeta {\text{ dx}}\quad \forall \zeta \in {H^1}(\Omega ).$

Testing with $\zeta = p$ and applying the Hölder's inequality and the Poincaré inequality (1.9) gives

${\left\| {\nabla p} \right\|_{{L^2}}} \leq {\left\| {(\mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}} + \frac{{{C_p}}}{K}{\left\| {{\Gamma _v}} \right\|_{{L^2}}}.$

(4.22)

Applying Hölder's inequality and the Sobolev embedding $H^{1} \subset L^{6}$ yields that

${\left\| {(\mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}} \leq {\left\| {\mu + \chi \sigma } \right\|_{{L^6}}}{\left\| {\nabla \varphi } \right\|_{{L^3}}} \leq {C_{{\text{Sob}}}}{\left\| {\mu + \chi \sigma } \right\|_{{H^1}}}{\left\| {\nabla \varphi } \right\|_{{L^3}}}.$

(4.23)

By the Gagliardo--Nirenberg inequality (1.10) with parameters $j = 0$, $p = 3$, $r = 2$, $m = 2$, $d = 3$ and $q = 2$,

${\left\| {\nabla \varphi } \right\|_{{L^3}}} \leq C\left\| {\nabla \varphi } \right\|_{{H^2}}^{\frac{1}{4}}\left\| {\nabla \varphi } \right\|_{{L^2}}^{\frac{3}{4}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{1}{4}}\left\| \varphi \right\|_{{H^1}}^{\frac{3}{4}},$

(4.24)

where $C > 0$ is a constant depending only on $\Omega$. Then, the boundedness of $\mu, \sigma$ in $L^{2}(0, T;H^{1})$ and $\varphi$ in $L^{2}(0, T;H^{3}) \cap L^{\infty}(0, T;H^{1})$ leads to

$\begin{gathered} \int_0^T {} \left\| {(\mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}^{\frac{8}{5}}{\text{dt}} \leq C\int_0^T {} \left\| {\mu + \chi \sigma } \right\|_{{H^1}}^{\frac{8}{5}}\left\| \varphi \right\|_{{H^3}}^{\frac{2}{5}}\left\| \varphi \right\|_{{H^1}}^{\frac{6}{5}}{\text{dt}} \\ \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{6}{5}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{5}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{2}{5}}. \\ \end{gathered}$

By (4.22) we find that

$\begin{gathered} \int_0^T {} \left\| {\nabla p} \right\|_{{L^2}}^{\frac{8}{5}}{\text{dt}} \leq \int_0^T {} \left\| {(\mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}^{\frac{8}{5}} + \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}}^{\frac{8}{5}}{\text{dt}} \\ \leq C\left[\varphi \right]_{{L^\infty }({H^1})}^{\frac{6}{5}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{5}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{2}{5}} + \frac{C}{K}{T^{\frac{1}{5}}}\left\| {{\Gamma _v}} \right\|_{{L^2}({L^2})}^{\frac{8}{5}}, \\ \end{gathered}$

where the positive constant $C$ depends only on $\Omega$. As $\bar p{\text{ = 0}}$, by the Poincarö inequality (1.9), we see that

$\left\| p \right\|_{{L^{\frac{8}{5}}}({H^1})}^{\frac{8}{5}} \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{6}{5}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{5}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{2}{5}} + \frac{C}{K}{T^{\frac{1}{5}}}\left\| {{\Gamma _v}} \right\|_{{L^2}({L^2})}^{\frac{8}{5}},$

for some positive constant $C$ depending only on $\Omega$. Next, we see that

$\begin{gathered} {\left\| {{\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right\|_{{L^2}}} \leq {\left\| {(\mu + \chi \sigma )\Delta \varphi } \right\|_{{L^2}}} + {\left\| {\nabla (\mu + \chi \sigma ) \cdot \nabla \varphi } \right\|_{{L^2}}} \\ \leq {\left\| {\mu + \chi \sigma } \right\|_{{L^6}}}{\left\| {\Delta \varphi } \right\|_{{L^3}}} + {\left\| {\nabla (\mu + \chi \sigma )} \right\|_{{L^2}}}{\left\| {\nabla \varphi } \right\|_{{L^\infty }}}. \\ \end{gathered}$

By the Gagliardo-Nirenberg inequality (1.10), we find that

$\begin{gathered} {\left\| {{D^2}\varphi } \right\|_{{L^3}}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{3}{4}}\left\| \varphi \right\|_{{L^6}}^{\frac{1}{4}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{3}{4}}\left\| \varphi \right\|_{{H^1}}^{\frac{1}{4}}, \hfill \\ {\left\| {\nabla \varphi } \right\|_{{L^\infty }}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{3}{4}}\left\| \varphi \right\|_{{L^6}}^{\frac{1}{4}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{3}{4}}\left\| \varphi \right\|_{{H^1}}^{\frac{1}{4}}, \hfill \\ \end{gathered}$

(4.25)

and so, we have

${\left\| {{\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right\|_{{L^2}}} \leq C{\left\| {\mu + \chi \sigma } \right\|_{{H^1}}}\left\| \varphi \right\|_{{H^3}}^{\frac{3}{4}}\left\| \varphi \right\|_{{H^1}}^{\frac{1}{4}}.$

(4.26)

That is, ${\text{div}} (( \mu + \chi \sigma) \nabla \varphi) \in L^{2}$. Since by assumption $\Gamma_{{v}} \in L^{2}_{0}$, using elliptic regularity theory, we find that $p(t, \cdot) \in H^{2}$ for a.e. $t$ and there exists a constant $C$ depending only on $\Omega$, such that

${\left\| p \right\|_{{H^2}}} \leq C\left( {{{\left\| p \right\|}_{{H^1}}} + {{\left\| {{\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right\|}_{{L^2}}} + {K^{ - 1}}{{\left\| {{\Gamma _v}} \right\|}_{{L^2}}}} \right).$

(4.27)

Furthermore, from (4.26), we see that

$\begin{gathered} \int_0^T {} {\text{div}}((\mu + \chi \sigma )\nabla \varphi )_{{L^2}}^{\frac{8}{7}}{\text{dt}} \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{2}{7}}\int_0^T {} \left\| {\mu + \chi \sigma } \right\|_{{H^1}}^{\frac{8}{7}}\left\| \varphi \right\|_{{H^3}}^{\frac{6}{7}}{\text{dt}} \\ \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{2}{7}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{7}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{6}{7}}, \\ \end{gathered}$

and so for some positive constant $C$ depending only on $\Omega$, it holds that

$\begin{gathered} \int_0^T {} \left\| p \right\|_{{H^2}}^{\frac{8}{7}}{\text{dt}} \leq C{T^{\frac{3}{7}}}{K^{ - \frac{8}{7}}}\left\| {{\Gamma _v}} \right\|_{{L^2}({L^2})}^{\frac{8}{7}} + C{T^{\frac{2}{7}}}\left\| p \right\|_{{L^{\frac{8}{5}}}({H^1})}^{\frac{8}{7}} \\ + C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{2}{7}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{7}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{6}{7}}. \\ \end{gathered}$

(4.28)

For the velocity field ${v}$ we estimate as follows. Let $1 \leq i, j \leq 3$ be fixed, we obtain from (4.25),

$\begin{gathered} {\left\| {{D_i}{v_j}} \right\|_{{L^2}}} = K{\left\| {{D_i}{D_j}p - {D_i}(\mu + \chi \sigma ){D_j}\varphi - (\mu + \chi \sigma ){D_i}{D_j}\varphi } \right\|_{{L^2}}} \\ \leq K\left( {{{\left\| p \right\|}_{{H^2}}} + {{\left\| {\nabla (\mu + \chi \sigma )} \right\|}_{{L^2}}}{{\left\| {\nabla \varphi } \right\|}_{{L^\infty }}} + {{\left\| {\mu + \chi \sigma } \right\|}_{{L^6}}}{{\left\| {{D^2}\varphi } \right\|}_{{L^3}}}} \right) \\ \leq K\left( {{{\left\| p \right\|}_{{H^2}}} + C{{\left\| {\mu + \chi \sigma } \right\|}_{{H^1}}}\left\| \varphi \right\|_{{H^3}}^{\frac{3}{4}}\left\| \varphi \right\|_{{H^1}}^{\frac{1}{4}}} \right). \\ \end{gathered}$

(4.29)

Applying the same calculation as in (4.28) yields

$\begin{gathered} \int_0^T {} \left\| {{\text{D}}v} \right\|_{{L^2}}^{\frac{8}{7}}{\text{dt}} \leq CK\int_0^T {} \left\| p \right\|_{{H^2}}^{\frac{8}{7}} + \left\| {\mu + \chi \sigma } \right\|_{{H^1}}^{\frac{8}{7}}\left\| \varphi \right\|_{{H^3}}^{\frac{6}{7}}\left\| \varphi \right\|_{{H^1}}^{\frac{2}{7}}{\text{dt}} \\ \leq CK\left( {\left\| p \right\|_{{L^{\frac{8}{7}}}({H^2})}^{\frac{8}{7}} + \left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{8}{7}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{6}{7}}\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{2}{7}}} \right), \\ \end{gathered}$

for some positive constant $C$ depending only on $\Omega$.

---

4.1.2 Source term from the Darcy system

To estimate the third source term

$\int_\Omega {{\Gamma _v}} \left( {{p_k} - {\mu _k}{\varphi _k} - \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right){\text{ dx}} = \int_\Omega {{\Gamma _v}} \left( {{p_k} - {{\bar \mu }_k}{\varphi _k} + ({\mu _k} - {{\bar \mu }_k}){\varphi _k} - \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right){\text{ dx}}$

of the energy equality we use Hölder’s inequality to obtain

$\begin{gathered} \left| {\int_\Omega {{\Gamma _v}} \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2} + {\Gamma _v}{\varphi _k}({\mu _k} - {{\bar \mu }_k}){\text{dx}}} \right| \leq \frac{D}{2}{\left\| {{\Gamma _v}} \right\|_{{L^2}}}\left\| {{\sigma _k}} \right\|_{{L^4}}^2 \\ + {\left\| {{\Gamma _v}} \right\|_{{L^{\frac{3}{2}}}}}{\left\| {{\mu _k} - \bar \mu } \right\|_k}_{{L^6}}{\left\| \varphi \right\|_k}_{{L^6}}. \\ \end{gathered}$

By the Gagliardo--Nirenberg inequality (1.10) with $j = 0$, $r = 2$, $m = 1$, $p = 4$, $q = 2$ and $\alpha = \frac{3}{4}$, we have

$\left\| {{\sigma _k}} \right\|_{{L^4}}^2 \leq C\left\| {{\sigma _k}} \right\|_{{H^1}}^{\frac{3}{2}}\left\| {{\sigma _k}} \right\|_{{L^2}}^{\frac{1}{2}} = C\left( {\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + \left\| {{\sigma _k}} \right\|_{{L^2}}^{\frac{1}{2}}{{\left\| {\nabla \sigma } \right\|}_k}_{{L^2}}^{\frac{3}{2}}} \right).$

By Young's inequality with Hölder exponents (i.e., $ab \leq \frac{\varepsilon }{p} a^{p} + \frac{\varepsilon ^{-q/p}}{q} b^{q}$ for $\frac{1}{p} + \frac{1}{q} = 1$ and $\varepsilon > 0$), we find that

$\frac{D}{2}{\left\| {{\Gamma _v}} \right\|_{{L^2}}}\left\| {{\sigma _k}} \right\|_{{L^4}}^2 \leq C\left( {{{\left\| {{\Gamma _v}} \right\|}_{{L^2}}}\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4\left\| {{\sigma _k}} \right\|_{{L^2}}^2} \right) + \frac{{{n_0}{D^2}}}{4}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2,$

for some positive constant $C$ depending only in $n_{0}$, $D$ and $\Omega$. Then, by (4.17) we have

$\begin{gathered} \left| {\int_\Omega {{\Gamma _v}} \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2} + {\Gamma _v}{\varphi _k}({\mu _k} - \overline {{\mu _k}}){\text{dx}}} \right| \hfill \\ \leq \frac{{{n_0}{D^2}}}{4}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + C\left( {1 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + {\left\| {{\Gamma _v}} \right\|_{{L^{\frac{3}{2}}}}}c({C_p}, {C_{{\text{Sob}}}}){\left\| {\nabla {\mu _k}} \right\|_{{L^2}}}{\left\| {{\varphi _k}} \right\|_{{H^1}}} \hfill \\ \leq \frac{{{n_0}{D^2}}}{4}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + C\left( {1 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + \frac{{{m_0}}}{8}\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + C\left\| {{\Gamma _v}} \right\|_{{L^2}}^2{\varphi _k}_{{H^1}}^2, \hfill \\ \end{gathered}$

where the positive constant $C$ depends only on $\Omega$, $m_{0}$, $n_{0}$ and $D$. Here we point out that the assumption $\Gamma_{{v}} \in L^{4}(0, T;L^{2}_{0})$ is needed. For the remainder term $\Gamma_{{v}} (p_{k} - \overline{\mu_{k}} \varphi_{k})$, we find that

$\begin{gathered} {p_k} - \overline {{\mu _k}} {\varphi _k} \hfill \\ = \left( {{{\left( { - {\Delta _N}} \right)}^{ - 1}}\left( {\frac{1}{K}{\Gamma _v} - {\text{div}}(({\mu _k} - \overline {{\mu _k}} + \chi {\sigma _k})\nabla {\varphi _k}) - \overline {{\mu _k}} {\text{div}}\nabla ({\varphi _k} - \overline {{\varphi _k}} )} \right)} \right) - \overline {{\mu _k}} {\varphi _k} \hfill \\ = \left( {{{\left( { - {\Delta _N}} \right)}^{ - 1}}\left( {\frac{1}{K}{\Gamma _v} - {\text{div}}(({\mu _k} - \overline {{\mu _k}} + \chi {\sigma _k})\nabla {\varphi _k})} \right)} \right) - \overline {{\mu _k}} \overline {{\varphi _k}} , \hfill \\ \end{gathered}$

where we used

${( - {\Delta _N})^{ - 1}}\left( { - \overline {{\mu _k}} {\text{div}}\nabla ({\varphi _k} - \overline {{\varphi _k}} )} \right) = \overline {{\mu _k}} ({\varphi _k} - \overline {{\varphi _k}} ).$

Then, by

$\int_\Omega {{\Gamma _v}} \overline {{\varphi _k}} \overline {{\mu _k}} {\text{ dx}} = \overline {{\mu _k}} \overline {{\varphi _k}} \int_\Omega {{\Gamma _v}} {\text{ dx}} = 0,$

it holds that

$\int_\Omega {{\Gamma _v}} ({p_k} - \overline {{\mu _k}} {\varphi _k}){\text{ dx}} = \int_\Omega {{\Gamma _v}} \left( {{{\left( { - {\Delta _N}} \right)}^{ - 1}}\left( {\frac{1}{K}{\Gamma _v} - {\text{div}}(({\mu _k} - \overline {{\mu _k}} + \chi {\sigma _k})\nabla {\varphi _k})} \right)} \right).$

Applying the calculations in the proof of Lemma 4.1 (specifically (4.10), (4.16) and (4.17)), H¨older’s inequality and Young’s inequality, we find that

$\begin{gathered} \left| {\int_\Omega {{\Gamma _v}} ({p_k} - \overline {{\mu _k}} {\varphi _k}){\text{dx}}} \right| \hfill \\ \leq \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 + C{\left\| {{\Gamma _v}} \right\|_{{L^2}}}\left( {{{\left\| {\nabla {\mu _k}} \right\|}_{{L^2}}} + \chi {{\left\| {{\sigma _k}} \right\|}_{{H^1}}}} \right){\left\| {\nabla {\varphi _k}} \right\|_{{L^{\frac{3}{2}}}}} \hfill \\ \leq \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 + \frac{{{m_0}}}{8}\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + \frac{{{n_0}{D^2}}}{4}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + {\chi ^2})\left\| {{\Gamma _v}} \right\|_{{L^2}}^2\left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2, \hfill \\ \end{gathered}$

where $C$ is a positive constant depending only on $\left| \Omega \right|$, $C_{p}$, $C_{\mathrm{Sob}}$, $D$, $n_{0}$ and $m_{0}$. Here we point out that if we applied (4.18) instead of (4.8) then we obtain a term containing ${\left\| {\nabla \varphi } \right\|_{{L^3}}}$ on the right-hand side and this cannot be controlled by the left-hand side of (4.4). Using (2.4) we have

$\left\| {{\varphi _k}} \right\|_{{L^2}}^2 \leq \frac{1}{{{R_1}}}{\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + \frac{{{R_2}}}{{{R_1}}}\left| \Omega \right|.$

(4.30)

Then, we obtain the following estimate

$\begin{gathered} \left| {\int_\Omega {{\Gamma _v}} \left( {{p_k} - {\mu _k}{\varphi _k} - \frac{D}{2}{{\left| {{\sigma _k}} \right|}^2}} \right){\text{dx}}} \right| \hfill \\ \leq C\left( {\frac{1}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 + 1} \right) + \frac{{{n_0}{D^2}}}{2}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{{{m_0}}}{4}\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 \hfill \\ + C\left( {1 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C{\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + C(1 + {\chi ^2})\left\| {{\Gamma _v}} \right\|_{{L^2}}^2\left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2, \hfill \\ \end{gathered}$

(4.31)

for some positive constant $C$ depending only on $R_{1}$, $R_{2}$, $\Omega$, $m_{0}$, $n_{0}$ and $D$. Here we point out that it is crucial for the source term $\Gamma_{{v}}$ to be prescribed and is not a function of $\varphi$, $\mu$ and $\sigma$, otherwise the product term $\left\| {{\Gamma _v}} \right\|_{{L^2}}^4\left\| {{\sigma _k}} \right\|_{{L^2}}^2$ and $\left\| {{\Gamma _v}} \right\|_{{L^2}}^2\left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2$ cannot be controlled in the absence of any a priori estimates. For the remaining source term

$\int_\Omega {{\Gamma _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}N_{, \sigma }^k{\text{ dx}}$

we split the analysis into two cases and combine with (4.31) to derive an energy inequality.

---

4.1.3 Energy inequality for non-negative $\Theta_{\varphi}$

Suppose $\Theta_{\varphi}$ is non-negative and bounded, and $\Psi$ is a potential that satisfies (2.5). We will estimate the mean of $\mu_{k}$ by setting $j = 1$ in (3.2b), and using the growth condition (2.5) to obtain

$\begin{gathered} {\left| {\int_\Omega {{\mu _k}} {\text{dx}}} \right|^2} = {\left| {\int_\Omega A \Psi '({\varphi _k}) - \chi {\sigma _k}{\text{dx}}} \right|^2} \leq 2{A^2}\left\| {\Psi '({\varphi _k})} \right\|_{{L^1}}^2 + 2{\chi ^2}\left\| {{\sigma _k}} \right\|_{{L^1}}^2 \\ \leq 2{A^2}R_4^2{\left( {\left| \Omega \right| + {{\left| \Omega \right|}^{\frac{1}{2}}}{{\left\| {{\varphi _k}} \right\|}_{{L^2}}}} \right)^2} + 2{\chi ^2}\left| \Omega \right|\left\| {{\sigma _k}} \right\|_{{L^2}}^2 \\ \leq C(A, {R_4}, \left| \Omega \right|) + 4{A^2}R_4^2\left| \Omega \right|\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + 2{\chi ^2}\left| \Omega \right|\left\| {{\sigma _k}} \right\|_{{L^2}}^2. \\ \end{gathered}$

Then, by the Poincarö inequality (1.9) and the growth condition (2.4), we find that

$\begin{gathered} \left\| {{\mu _k}} \right\|_{{L^2}}^2 \leq 2C_P^2\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + 2\left| \Omega \right|{\left| {\overline {{\mu _k}} } \right|^2} \\ \leq 2C_p^2\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + 8{A^2}R_4^2\left\| {{\varphi _k}} \right\|_{{L^2}}^2 + 4{\chi ^2}\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(A, {R_4}, \left| \Omega \right|) \\ \leq 2C_p^2\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + \frac{{8{A^2}R_4^2}}{{{R_1}}}{\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + 4{\chi ^2}\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(A, {R_4}, {R_1}, {R_2}, \left| \Omega \right|). \\ \end{gathered}$

(4.32)

Note that by the specific form (2.1) for $\Gamma_{\varphi}$ we have that

${\Gamma _{\varphi , k}}{\mu _k} = {\Lambda _\varphi }({\varphi _k}, {\sigma _k}){\mu _k} - {\Theta _\varphi }({\varphi _k}, {\sigma _k}){\text{ }}{\left| {{\mu _k}} \right|^2}.$

Moving the non-negative term ${\Theta _\varphi }({\varphi _k}, {\sigma _k}){\left| {{\mu _k}} \right|^2}$ to the left-hand side of (4.4) and subsequently neglecting it, we estimate the remainder using the growth condition (2.3) and Hölder's inequality as follows (here we use the notation $\Lambda_{\varphi, k} := \Lambda_{\varphi}(\varphi_{k}, \sigma_{k})$),

$\begin{gathered} \left| {\int_\Omega {{\Lambda _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}(D{\sigma _k} + \chi (1 - {\varphi _k})){\text{dx}}} \right| \hfill \\ \leq {\left\| {{\Lambda _{\varphi , k}}} \right\|_{{L^2}}}{\left\| {{\mu _k}} \right\|_{{L^2}}} + \left( {{{\left\| {{\Lambda _{S, k}}} \right\|}_{{L^2}}} + {R_0}{{\left\| {{\mu _k}} \right\|}_{{L^2}}}} \right){\left\| {D{\sigma _k} + \chi (1 - {\varphi _k})} \right\|_{{L^2}}} \hfill \\ \leq C\left( {1 + \chi + (1 + \chi ){{\left\| {{\varphi _k}} \right\|}_{{L^2}}} + (1 + D){{\left\| {{\sigma _k}} \right\|}_{{L^2}}}} \right){\left\| {{\mu _k}} \right\|_{{L^2}}} \hfill \\ + C\left( {1 + {{\left\| {{\varphi _k}} \right\|}_{{L^2}}} + {{\left\| \sigma \right\|}_{{L^2}}}} \right)\left( {\chi {{\left| \Omega \right|}^{\frac{1}{2}}} + D{{\left\| \sigma \right\|}_{{L^2}}} + \chi {{\left\| {{\varphi _k}} \right\|}_{{L^2}}}} \right) \hfill \\ \end{gathered}$

(4.33)

where $C$ is a positive constant depending only on $R_{0}$ and $\left| \Omega \right|$. By Young's inequality, (4.32) and (4.30), we have

$\begin{gathered} \left| {\int_\Omega {{\Lambda _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}(D{\sigma _k} + \chi (1 - {\varphi _k})){\text{dx}}} \right| \hfill \\ \leq \frac{{{m_0}}}{{8C_p^2}}\left\| {{\mu _k}} \right\|_{{L^2}}^2 + C(1 + \chi + D + {\chi ^2})\left\| {{\varphi _k}} \right\|_{{L^2}}^2 + C{(1 + \chi + D)^2}\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + \chi + {\chi ^2}) \hfill \\ \leq \frac{{{m_0}}}{4}\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + C(1 + {\chi ^2})\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + {\chi ^2}){\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + C(1 + {\chi ^2}), \hfill \\ \end{gathered}$

(4.34)

for some positive constant $C$ depending only on $\left| \Omega \right|$, $R_{0}$, $R_{1}$, $R_{2}$, $R_{4}$, $A$, $D$, $C_{p}$ and $m_{0}$. Using the fact that

${\left\| {D\nabla \sigma } \right\|_{{L^2}}} \leq {\left\| {\nabla (D\sigma + \chi (1 - \varphi ))} \right\|_{{L^2}}} + {\left\| {\chi \nabla \varphi } \right\|_{{L^2}}},$

we now estimate the right-hand side of (4.4) using (4.6), (4.31) and (4.34), which leads to

$\begin{gathered} \frac{{\text{d}}}{{{\text{dt}}}}\int_\Omega A \Psi ({\varphi _k}) + \frac{B}{2}{\left| {\nabla {\varphi _k}} \right|^2} + \frac{D}{2}{\left| {{\sigma _k}} \right|^2} + \chi {\sigma _k}(1 - {\varphi _k}){\text{dx}} \hfill \\ \quad + \frac{{{m_0}}}{2}\left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + \frac{{{n_0}{D^2}}}{2}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{{Db}}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2 \hfill \\ \quad \leq C(1 + b)(1 + {\chi ^2}){\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + C\left( {\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 + b} \right)(1 + {\chi ^2})\left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2 \hfill \\ \quad + C\left( {1 + {\chi ^2} + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + b)(1 + {\chi ^2}) + \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 \hfill \\ \quad + bC(1 + {\chi ^2})\left\| {{\sigma _\infty }} \right\|_{{L^2}(\partial \Omega )}^2, \hfill \\ \end{gathered}$

(4.35)

for some positive constant $C$ not depending on $\Gamma_{{v}}$, $K$, $b$ and $\chi$. Integrating (4.35) with respect to $t$ from $0$ to $s \in (0, T]$ leads to

$\begin{gathered} A{\left\| {\Psi ({\varphi _k}(s))} \right\|_{{L^1}}} + \frac{B}{2}\left\| {\nabla {\varphi _k}(s)} \right\|_{{L^2}}^2 + \frac{D}{2}\left\| {{\sigma _k}(s)} \right\|_{{L^2}}^2 + \int_\Omega \chi {\sigma _k}(s)(1 - {\varphi _k}(s)){\text{dx}} \hfill \\ \quad + \int_0^s {\frac{{{m_0}}}{2}} \left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + \frac{{{n_0}{D^2}}}{2}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{{Db}}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2{\text{dt}} \hfill \\ \quad \leq \int_0^s C (1 + b)(1 + {\chi ^2})\left( {1 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left( {{{\left\| {\Psi ({\varphi _k})} \right\|}_{{L^1}}} + \left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2 + \left\| {{\sigma _k}} \right\|_{{L^2}}^2} \right){\text{dt}} \hfill \\ \quad + C(1 + b)(1 + {\chi ^2})T + \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}(0, T;{L^2})}^2 + Cb(1 + {\chi ^2})\left\| {{\sigma _\infty }} \right\|_{{L^2}(0, T;{L^2}(\partial \Omega ))}^2 \hfill \\ \quad + C\left( {{{\left\| {\Psi ({\varphi _0})} \right\|}_{{L^1}}} + \left\| {{\varphi _0}} \right\|_{{H^1}}^2 + \left\| {{\sigma _0}} \right\|_{{L^2}}^2} \right), \hfill \\ \end{gathered}$

(4.36)

for some positive constant $C$ independent of $\Gamma_{{v}}$, $K$, $\chi$ and $b$. Here we used $\sigma_{0} \in L^{2}$ and $\varphi_{0} \in H^{1}$, which implies by the growth condition (2.5) that $\Psi(\varphi_{0}) \in L^{1}$. Next, by Hölder's inequality and Young's inequality we have

$\begin{gathered} \left| {\int_\Omega \chi {\sigma _k}(x, s)(1 - {\varphi _k}(x, s)){\text{dx}}} \right| \leq \frac{{2D}}{8}\left\| {{\sigma _k}(s)} \right\|_{{L^2}}^2 + \frac{{2{\chi ^2}\left| \Omega \right|}}{D} + \frac{{2{\chi ^2}}}{D}\left\| {{\varphi _k}(s)} \right\|_{{L^2}}^2 \hfill \\ \quad \leq \frac{D}{4}\left\| {{\sigma _k}(s)} \right\|_{{L^2}}^2 + \frac{{2{\chi ^2}}}{{D{R_1}}}{\left\| {\Psi ({\varphi _k}(s))} \right\|_{{L^1}}} + \frac{{2{\chi ^2}\left| \Omega \right|}}{D}(1 + {R_2}). \hfill \\ \end{gathered}$

(4.37)

Substituting (4.37) into (4.36) then yields

$\begin{gathered} \min \left( {A - \frac{{2{\chi ^2}}}{{D{R_1}}}, \frac{B}{2}, \frac{D}{4}} \right)\left( {{{\left\| {\Psi ({\varphi _k}(s))} \right\|}_{{L^1}(\Omega )}} + \left\| {\nabla {\varphi _k}(s)} \right\|_{{L^2}(\Omega )}^2 + \left\| {{\sigma _k}(s)} \right\|_{{L^2}(\Omega )}^2} \right) \hfill \\ + \int_0^s {} \left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + \left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{b}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2{\text{dt}} \hfill \\ \leq \int_0^s C (1 + b)(1 + {\chi ^2})\left( {1 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left( {{{\left\| {\Psi ({\varphi _k})} \right\|}_{{L^1}}} + \left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2 + \left\| {{\sigma _k}} \right\|_{{L^2}}^2} \right){\text{dt}} \hfill \\ + C(1 + b)(1 + {\chi ^2})(1 + T) + \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}(0, T;{L^2})}^2, \hfill \\ \end{gathered}$

(4.38)

for some positive constant $C$ independent of $\Gamma_{{v}}$, $K$, $b$ and $\chi$. Setting

$\begin{gathered} \alpha : = C(1 + b)(1 + {\chi ^2})(1 + T) + \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}(0, T;{L^2})}^2, \hfill \\ \beta : = C(1 + b)(1 + {\chi ^2})\left( {1 + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right) \in {L^1}(0, T), \hfill \\ \end{gathered}$

(4.39)

and noting that

$\alpha \left( {1 + \int_0^s \beta (t)\exp \left( {\int_0^t \beta (r){\text{ dr}}} \right){\text{ dt}}} \right) \leq \alpha \left( {1 + {{\left\| \beta \right\|}_{{L^1}(0, T)}}\exp \left( {{{\left\| \beta \right\|}_{{L^1}(0, T)}}} \right)} \right) < \infty ,$

we find after applying the Gronwall inequality (1.11) to (4.38) leads to

$\begin{gathered} \mathop {\sup }\limits_{s \in (0, T]} \left( {{{\left\| {\Psi ({\varphi _k}(s))} \right\|}_{{L^1}}} + \left\| {\nabla {\varphi _k}(s)} \right\|_{{L^2}}^2 + \left\| {{\sigma _k}(s)} \right\|_{{L^2}}^2} \right) \hfill \\ + \int_0^T {} \left\| {\nabla {\mu _k}} \right\|_{{L^2}}^2 + \left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{b}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2{\text{dt}} \leq \mathcal{E}, \hfill \\ \end{gathered}$

(4.40)

where we recall that $\mathcal{E}$ denotes a constant that is uniformly bounded for $b, \chi \in (0, 1]$ and is also uniformly bounded for $K \in (0, 1]$ when $\Gamma_{{v}} = 0$.

---

4.1.4 Energy inequality for positive $\Theta_{\varphi}$

Suppose $\Theta_{\varphi}$ satisfies (2.6) and $\Psi$ is a potential satisfying the growth condition (2.7). Similar to the previous case, we see that the specific form for $\Gamma_{\varphi}$ leads to

${\Gamma _{\varphi , k}}{\mu _k} = {\Lambda _\varphi }({\varphi _k}, {\sigma _k}){\mu _k} - {\Theta _\varphi }({\varphi _k}, {\sigma _k}){\text{ }}{\left| {{\mu _k}} \right|^2}.$

We move the term ${\Theta _\varphi }({\varphi _k}, {\sigma _k}){\left| {{\mu _k}} \right|^2}$ to the left-hand side of (4.4) and estimate the remainder as in (4.33). Using Young's inequality differently and also (4.30), we have

$\begin{gathered} \left| {\int_\Omega {{\Lambda _{\varphi , k}}} {\mu _k} - {\mathcal{S}_k}(D{\sigma _k} + \chi (1 - {\varphi _k})){\text{dx}}} \right| \hfill \\ \quad \leq \frac{{{R_5}}}{2}\left\| {{\mu _k}} \right\|_{{L^2}}^2 + C(1 + \chi + D + {\chi ^2})\left\| {{\varphi _k}} \right\|_{{L^2}}^2 + C{(1 + \chi + D)^2}\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + \chi + {\chi ^2}) \hfill \\ \quad \leq \frac{{{R_5}}}{2}\left\| {{\mu _k}} \right\|_{{L^2}}^2 + C(1 + {\chi ^2})\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + {\chi ^2}){\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + C(1 + {\chi ^2}), \hfill \\ \end{gathered}$

(4.41)

for some positive constant $C$ depending only on $\left| \Omega \right|$, $R_{5}$, $R_{1}$, $R_{2}$, $A$, $D$, and $C_{p}$. Using (4.6), (4.31), (4.41) and the lower bound $\Theta_{\varphi} \geq R_{5}$, instead of (4.35) we obtain from (4.4)

$\begin{gathered} \frac{{\text{d}}}{{{\text{dt}}}}\int_\Omega A \Psi ({\varphi _k}) + \frac{B}{2}{\left| {\nabla {\varphi _k}} \right|^2} + \frac{D}{2}{\left| \sigma \right|_k}^2 + \chi {\sigma _k}(1 - {\varphi _k}){\text{dx}} \hfill \\ \quad + \frac{{{R_5}}}{2}\left\| {{\mu _k}} \right\|_{{L^2}}^2 + \frac{{{m_0}}}{2}\left\| {\nabla {\mu _k}} \right\|_{{L^2}(\Omega )}^2 + \frac{{{n_0}{D^2}}}{2}\left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{{Db}}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2 \hfill \\ \quad \leq C(1 + b)(1 + {\chi ^2}){\left\| {\Psi ({\varphi _k})} \right\|_{{L^1}}} + C\left( {\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 + b} \right)(1 + {\chi ^2})\left\| {\nabla {\varphi _k}} \right\|_{{L^2}}^2 \hfill \\ \quad + C\left( {1 + {\chi ^2} + \left\| {{\Gamma _v}} \right\|_{{L^2}}^4} \right)\left\| {{\sigma _k}} \right\|_{{L^2}}^2 + C(1 + b)(1 + {\chi ^2}) + \frac{C}{K}\left\| {{\Gamma _v}} \right\|_{{L^2}}^2 \hfill \\ \quad + Cb(1 + {\chi ^2})\left\| {{\sigma _\infty }} \right\|_{{L^2}(\partial \Omega )}^2, \hfill \\ \end{gathered}$

(4.42)

for some positive constant $C$ independent of $\Gamma_{{v}}$, $K$, $b$ and $\chi$. We point out the main difference between (4.35) and the above is the appearance of the term $\frac{{{R_5}}}{2}\left\| {{\mu _k}} \right\|_{{L^2}}^2$ on the left-hand side. The positivity of $\Theta_{\varphi}$ allows us to absorb the $\left\| {{\mu _k}} \right\|_{{L^2}}^2$ term on the right-hand side of (4.41) and thus we do not need to use (4.32), which was the main reason why $\Psi$ has to be a quadratic potential for a non-negative $\Theta_{\varphi}$. Then, applying a similar argument as in Section 4.1.3, we arrive at an analogous energy inequality to (4.40),

$\begin{gathered} \mathop {\sup }\limits_{s \in (0, T]} \left( {{{\left\| {\Psi ({\varphi _k}(s))} \right\|}_{{L^1}}} + \left\| {\nabla {\varphi _k}(s)} \right\|_{{L^2}}^2 + \left\| {{\sigma _k}(s)} \right\|_{{L^2}}^2} \right) \hfill \\ + \int_0^T {} \left\| {{\mu _k}} \right\|_{{H^1}}^2 + \left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{b}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2{\text{dt}} \leq \mathcal{E}. \hfill \\ \end{gathered}$

(4.43)

Using (4.32) and (4.30) applied to (4.40), and similarly using (4.30) applied to (4.43) we obtain

$\begin{gathered} \mathop {\sup }\limits_{s \in (0, T]} \left( {{{\left\| {\Psi ({\varphi _k}(s))} \right\|}_{{L^1}}} + \left\| {{\varphi _k}(s)} \right\|_{{H^1}}^2 + \left\| {{\sigma _k}(s)} \right\|_{{L^2}}^2} \right) \hfill \\ + \int_0^T {} \left\| {{\mu _k}} \right\|_{{H^1}}^2 + \left\| {\nabla {\sigma _k}} \right\|_{{L^2}}^2 + \frac{1}{K}\left\| {{v_k}} \right\|_{{L^2}}^2 + \frac{b}{2}\left\| {{\sigma _k}} \right\|_{{L^2}(\partial \Omega )}^2{\text{dt}} \leq \mathcal{E}. \hfill \\ \end{gathered}$

(4.44)

This a priori estimate implies that the Galerkin ansatz $\varphi_{k}$, $\mu_{k}$, $\sigma_{k}$ and ${v}_{k}$ can be extended to the interval $[0, T]$. To determine if $p_{k}$ can also be extended to the interval $[0, T]$ we require some higher order estimates for $\varphi_{k}$ in order to use (4.19).

---

4.2. Higher order estimates

Let $\Pi_{k}$ denote the orthogonal projection onto the finite-dimensional subspace $W_{k}$. From (3.2b) we may view $\varphi_{k}$ as the solution to the following elliptic equation

$- B\Delta u + u = {\mu _k} - A{\Pi _k}\left( {\Psi '(u)} \right) + \chi {\sigma _k} + u{\text{ in }}\Omega ,$

(4.45a)

${\partial _n}u = 0{\text{ on }}\partial \Omega .$

(4.45b)

For the case where $\Psi$ satisfies (2.5), as $\{\varphi_{k}\}_{k \in \mathbb{N}}$ is bounded in $L^{\infty}(0, T;H^{1})$, we have that $\{ \Psi'(\varphi_{k})\}_{k \in \mathbb{N}}$ is also bounded in $L^{\infty}(0, T;H^{1})$. Using the fact that our basis functions $\{w_{i}\}_{i \in \mathbb{N}}$ are the eigenfunctions of the inverse Neumann-Laplacian operator and is therefore orthogonal in $H^{1}$, and the Sobolev embedding $H^{1} \subset L^{r}$ for $r \in [1, 6]$, there exists a positive constant $C$ independent of $\varphi_{k}$ such that

${\left\| {{\Pi _k}(\Psi '({\varphi _k}))} \right\|_X} \leq C{\left\| {\Psi '({\varphi _k})} \right\|_X}\quad {\text{ for }}X = {H^1}{\text{ or }}{L^r}, 1 \leq r \leq 6.$

(4.46)

Then, this implies that $\{\Pi_{k} \left ( \Psi'(\varphi_{k} ) \right )\}_{k \in \mathbb{N}}$ is also bounded in $L^{\infty}(0, T;H^{1})$. As the right-hand side of (4.45a) belongs to $H^{1}$ for a.e. $t \in (0, T)$, and the boundary $\partial \Omega$ is $C^{3}$, by elliptic regularity theory, we have

${\left\| {{\varphi _k}} \right\|_{{L^2}({H^3})}} \leq C\left( {1 + {{\left\| {{\varphi _k}} \right\|}_{{L^2}({H^1})}} + {{\left\| {{\mu _k} + \chi {\sigma _k}} \right\|}_{{L^2}({H^1})}}} \right) \leq \mathcal{E},$

(4.47)

for some positive constant $C$ depending only on $\Omega$ and $R_{4}$. For the case where $\Psi$ satisfies (2.7), we employ a bootstrap argument from [18, x3.3]. The growth assumption (2.7) implies that

$\left| {\Psi '(y)} \right| \leq C\left( {1 + {{\left| y \right|}^m}} \right), \quad {\text{ }}\left| {\Psi ''(y)} \right| \leq C\left( {1 + {{\left| y \right|}^{m - 1}}} \right){\text{ for }}m \in [1, 5).$

(4.48)

For fixed $m \in [1, 5)$, we define a sequence of positive numbers $\{l_{j}\}_{j \in \mathbb{N}}$ by

$l_{1}m \leq 6, \quad l_{j+1} = \frac{6 l_{j}}{6 - (5-m)l_{j}}.$

It can be shown that $\{l_{j}\}_{j \in \mathbb{N}}$ is a strictly increasing sequence such that $l_{j} \to \infty$ as $j \to \infty$. The Gagliardo--Nirenberg inequality (1.10) then yields the following continuous embedding

$L^{2}(0, T;W^{2, l_{j}}) \cap L^{\infty}(0, T;L^{6}) \subset L^{2m}(0, T;L^{m l_{j+1}}).$

(4.49)

At the first step, the boundedness of $\{\varphi_{k}\}_{k \in \mathbb{N}}$ in $L^{\infty}(0, T;H^{1})$ yields

$\left\| {{\Pi _k}(\Psi '({\varphi _k}))} \right\|_{{L^{{l_1}}}}^2 \leq C\left( {1 + \left\| {{\varphi _k}} \right\|_{{L^6}}^{2m}} \right),$

which implies that $\{\Pi_{k}(\Psi'(\varphi_{k})) \}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T;L^{l_{1}})$. As the other terms on the right-hand side of (4.45) are bounded in $L^{2}(0, T;H^{1})$, elliptic regularity then yields that $\{\varphi_{k}\}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T;W^{2, l_{1}})$, and thus in $L^{2m}(0, T;L^{ml_{2}})$ by (4.49).

At the $j$-th step, we have $\{\varphi_{k}\}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T, W^{2, l_{j}}) \cap L^{2m}(0, T;L^{m l_{j+1}})$. Then, it holds that

$\left\| {{\Pi _k}(\Psi '({\varphi _k}))} \right\|_{{L^2}({L^{{l_{j + 1}}}})}^2 \leq C\left( {1 + \left\| {{\varphi _k}} \right\|_{{L^{2m}}({L^{m{l_{j + 1}}}})}^{2m}} \right),$

and so $\{\Pi_{k}(\Psi'(\varphi_{k})) \}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T;L^{l_{j+1}})$. Elliptic regularity then implies that $\{\varphi_{k}\}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T;W^{2, l_{j+1}})$.

We terminate the bootstrapping procedure once $l_{j} \geq 6$ for some $j \in \mathbb{N}$. This occurs after a finite number of steps as $\lim_{j \to \infty} l_{j} = \infty$. Altogether, we obtain that $\{\varphi_{k}\}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T;W^{2, 6})$. From (4.48) it holds that

${\left\| {\Psi ''({\varphi _k})\nabla {\varphi _k}} \right\|^2} \leq C\left( {1 + {{\left| {{\varphi _k}} \right|}^{2m - 2}}} \right){\left| {\nabla {\varphi _k}} \right|^2}{\text{ for }}m \in [1, 5),$

and by the following continuous embeddings obtain from the Gagliardo-Nirenberg inequality (1.10),

$L^{2}(0, T;W^{2, 6}) \cap L^{\infty}(0, T;H^{1}) \subset L^{2m}(0, T;L^{2m}) \cap L^{2m-2}(0, T;L^{\infty}) \text{ for } m \in [1, 5),$

we find that $\{\Pi_{k}(\Psi'(\varphi_{k}))\}_{k \in \mathbb{N}}$ is bounded in $L^{2}(0, T;H^{1})$. Applying elliptic regularity once more leads to the boundedness of $\{\varphi_{k}\}_{k \in \mathbb{N}}$ in $L^{2}(0, T;H^{3})$. Consequently, the hypotheses of Lemma 4.2 are satisfied and we obtain that

${\left\| {{p_k}} \right\|_{{L^{\frac{8}{5}}}({H^1})}} \leq \mathcal{E},$

which implies that the Galerkin ansatz $p_{k}$ can be extended to the interval $[0, T]$.

---

4.3. Estimates for the convection terms and the time derivatives

By the Gagliardo-Nirenberg inequality (1.10) with $j = 0$, $p = \infty$, $m = 3$, $r = 2$, $q = 2$ and $d = 3$, we have

${\left\| {{\varphi _k}} \right\|_{{L^\infty }}} \leq C\left\| {{\varphi _k}} \right\|_{{H^3}}^{\frac{1}{4}}\left\| {{\varphi _k}} \right\|_{{L^6}}^{\frac{3}{4}} \leq C\left\| {{\varphi _k}} \right\|_{{H^3}}^{\frac{1}{4}}\left\| {{\varphi _k}} \right\|_{{H^1}}^{\frac{3}{4}}.$

For any $\zeta \in L^{\frac{8}{3}}(0, T;H^{1})$ with coefficients $\{\zeta_{kj}\}_{1 \leq j \leq k} \subset \mathbb{R}^{k}$ such that $\Pi_{k} \zeta = \sum_{j=1}^{k} \zeta_{kj} w_{j}$, we can estimate

$\begin{gathered} \left| {\int_0^T {\int_\Omega {{\varphi _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{dx dt}}} \right| \leq \int_0^T {} {\left\| {{v_k}} \right\|_{{L^2}}}{\left\| {{\varphi _k}} \right\|_{{L^\infty }}}{\left\| {\nabla {\Pi _k}\zeta } \right\|_{{L^2}}}{\text{dt}} \\ \leq C\left\| {{\varphi _k}} \right\|_{{L^\infty }({H^1})}^{\frac{3}{4}}{\left\| {{v_k}} \right\|_{{L^2}({L^2})}}\left\| {{\varphi _k}} \right\|_{{L^2}({H^3})}^{\frac{1}{4}}{\left\| \zeta \right\|_{{L^{\frac{8}{3}}}({H^1})}}. \\ \end{gathered}$

(4.50)

Using (4.44) and (4.47), we find that

${\left\| {{\text{div}}({\varphi _k}{v_k})} \right\|_{{L^{\frac{8}{5}}}({{({H^1})}^*})}} \leq {K^{\frac{1}{2}}}\mathcal{E}.$

(4.51)

Next, multiplying (3.2a) by $\zeta_{kj}$, summing from $j = 1$ to $k$ and then integrating in time from $0$ to $T$ leads to

$\begin{gathered} \left| {\int_0^T {\int_\Omega {{\partial _t}} } {\varphi _k}\zeta {\text{dx dt}}} \right| \leq \int_0^T {{m_1}} {\left\| {\nabla {\mu _k}} \right\|_{{L^2}}}{\left\| {\nabla {\Pi _k}\zeta } \right\|_{{L^2}}}{\text{dt}} \\ + \int_0^T {} {\left\| {{\Gamma _{\varphi , k}}} \right\|_{{L^2}}}{\left\| {{\Pi _k}\zeta } \right\|_{{L^2}}} + {\left\| {{\text{div}}({\varphi _k}{v_k})} \right\|_{{{({H^1})}^*}}}{\left\| {{\Pi _k}\zeta } \right\|_{{H^1}}}{\text{dt}}. \\ \end{gathered}$

By (2.1), (2.2) and (4.44), we find that

${\left\| {{\Gamma _{\varphi , k}}} \right\|_{{L^2}({L^2})}} \leq C({R_0}, \left| \Omega \right|, T)\left( {1 + {{\left\| {{\varphi _k}} \right\|}_{{L^2}({L^2})}} + {{\left\| {{\mu _k}} \right\|}_{{L^2}({L^2})}} + {{\left\| {{\sigma _k}} \right\|}_{{L^2}({L^2})}}} \right) \leq \mathcal{E},$

and so, by Hölder’s inequality, we find that

$\left| {\int_0^T {\int_\Omega {{\partial _t}} } {\varphi _k}\zeta {\text{ dx dt}}} \right| \leq \left( {\mathcal{E}{T^{\frac{1}{8}}} + {{\left\| {{\text{div}}({\varphi _k}{v_k})} \right\|}_{{L^{\frac{8}{5}}}({{({H^1})}^*})}}} \right){\left\| \zeta \right\|_{{L^{\frac{8}{3}}}({H^1})}}.$

Taking the supremum over $\zeta \in L^{\frac{8}{3}}(0, T;H^{1})$ and using (4.44) and (4.51) yields that

${\left\| {{\partial _t}{\varphi _k}} \right\|_{{L^{\frac{8}{5}}}({{({H^1})}^*})}} \leq \mathcal{E}\left( {1 + {K^{\frac{1}{2}}}} \right),$

(4.52)

Similarly, by Hölder's inequality and the following Gagliardo--Nirenberg inequality (1.10) with $j = 0$, $r = 2$, $m = 1$, $p = \frac{10}{3}$, $q = 2$ and $d = 3$,

${\left\| f \right\|_{{L^{\frac{{10}}{3}}}}} \leq C\left\| f \right\|_{{H^1}}^{\frac{3}{5}}\left\| f \right\|_{{L^2}}^{\frac{2}{5}},$

which in turn implies that $\{\sigma_{k}\}_{k \in \mathbb{N}}$ is bounded uniformly in $L^{\frac{10}{3}}(Q)$. Then, we find that for any $\zeta \in L^{5}(0, T;W^{1, 5})$,

$\begin{gathered} \left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{dx dt}}} \right| \leq \int_0^T {} {\left\| {{\sigma _k}} \right\|_{{L^{\frac{{10}}{3}}}}}{\left\| {{v_k}} \right\|_{{L^2}}}{\left\| {\nabla \zeta } \right\|_{{L^5}}}{\text{dt}} \\ \leq {\left\| {{\sigma _k}} \right\|_{{L^{\frac{{10}}{3}}}(Q)}}{\left\| {{v_k}} \right\|_{{L^2}({L^2})}}{\left\| {\nabla \zeta } \right\|_{{L^5}({L^5})}}, \\ \end{gathered}$

(4.53)

and

${\left\| {{\text{div}}({\sigma _k}{v_k})} \right\|_{{L^{\frac{5}{4}}}({{({W^{1, 5}})}^*})}} \leq {K^{\frac{1}{2}}}\mathcal{E}.$

(4.54)

A similar calculation to (4.52) yields that

${\left\| {{\partial _t}{\sigma _k}} \right\|_{{L^{\frac{5}{4}}}({{({W^{1, 5}})}^*})}} \leq \mathcal{E}\left( {1 + {K^{\frac{1}{2}}}} \right).$

(4.55)

Remark 4.2. We may also use the Gagliardo-Nirenberg inequality to deduce that

${\left\| f \right\|_{{L^r}}} \leq C\left\| f \right\|_{{H^1}}^{\frac{{3(r - 2)}}{{2r}}}\left\| f \right\|_{{L^2}}^{\frac{{6 - r}}{{2r}}}{\text{ for any }}r \in (2, 6).$

Then, the computation (4.53) becomes

$\left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{ dx dt}}} \right| \leq C{\left\| {{v_k}} \right\|_{{L^2}({L^2})}}\left\| {{\sigma _k}} \right\|_{{L^\infty }({L^2})}^{\frac{{6 - r}}{{2r}}}\left\| {{\sigma _k}} \right\|_{{L^2}({H^1})}^{\frac{{3(r - 2)}}{{2r}}}{\left\| {\nabla \zeta } \right\|_{{L^{\frac{{4r}}{{6 - r}}}}({L^{\frac{{2r}}{{r - 2}}}})}},$

which implies that $\{ {\text{div}} (\sigma_{k} {v}_{k}) \}_{k \in \mathbb{N}}$ and $\{ \partial _{t}\sigma_{k}\}_{k \in \mathbb{N}}$ are bounded uniformly in

$L^{\frac{4r}{5r-6}}(0, T;(W^{1, \frac{2r}{r-2}})^{*}) \text{ for } r \in (2, 6).$

Note that the temporal exponent decreases while the spatial exponent increases as r increases, and they intersect at the point $r = \frac{10}{3}$

Here we point out that even with the improved regularity ${v}_{k} \in L^{\frac{8}{7}}(0, T;{H}^{1})$, we are unable to show ${\text{div}}(\sigma_{k} {v}_{k})$ is bounded in the dual space $(H^{1})^{*}$. Indeed, let $q$, $r > 1$ be constants yet to be determined such that $\frac{1}{q} + \frac{1}{r} = \frac{1}{2}$. Then, from Hölder's inequality we have

$\left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{ dx dt}}} \right| \leq \int_0^T {} {\left\| {{\sigma _k}} \right\|_{{L^q}}}{\left\| {{v_k}} \right\|_{{L^r}}}{\left\| {\nabla \zeta } \right\|_{{L^2}}}{\text{ dt}}.$

By the Gagliardo--Nirenberg inequality we have for $\alpha = \frac{3}{2} - \frac{3}{q} \leq 1$, $\beta = \frac{3}{2} - \frac{3}{r} \leq 1$,

$\begin{gathered} \left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{dx dt}}} \right| \leq C\int_0^T {} \left\| {{\sigma _k}} \right\|_{{L^2}}^{1 - \alpha }\left\| {{\sigma _k}} \right\|_{{H^1}}^\alpha \left\| {{v_k}} \right\|_{{H^1}}^\beta \left\| {{v_k}} \right\|_{{L^2}}^{1 - \beta }{\left\| {\nabla \zeta } \right\|_{{L^2}}}{\text{dt}} \hfill \\ \quad \leq C\left\| {{\sigma _k}} \right\|_{{L^\infty }({L^2})}^{1 - \alpha }\int_0^T {} \left\| {{\sigma _k}} \right\|_{{H^1}}^\alpha \left\| {{v_k}} \right\|_{{H^1}}^\beta \left\| {{v_k}} \right\|_{{L^2}}^{1 - \beta }{\left\| \zeta \right\|_{{H^1}}}{\text{dt}} \hfill \\ \quad \leq C\left\| {{\sigma _k}} \right\|_{{L^\infty }({L^2})}^{1 - \alpha }\left\| {{\sigma _k}} \right\|_{{L^{\alpha {x_1}}}({H^1})}^\alpha \left\| {{v_k}} \right\|_{{L^{\beta {x_2}}}({H^1})}^\beta \left\| {{v_k}} \right\|_{{L^{(1 - \beta ){x_3}}}({L^2})}^{1 - \beta }{\left\| \zeta \right\|_{{L^{{x_4}}}({H^1})}}, \hfill \\ \end{gathered}$

where

$\frac{1}{x_{1}} + \frac{1}{x_{2}} + \frac{1}{x_{3}} + \frac{1}{x_{4}} = 1, \quad \alpha x_{1} \leq 2, \quad \beta x_{2} \leq \frac{8}{7}, \quad (1-\beta) x_{3} \leq 2.$

(4.56)

Note that $\alpha = \frac{3}{2} - \frac{3}{q} = \frac{3}{r}$, and then substituting into the constraints (4.56) we find that

$\frac{1}{x_{1}} + \frac{1}{x_{2}} + \frac{1}{x_{3}} \geq \frac{\alpha}{2} + \frac{7}{8} \beta + \frac{1-\beta}{2} = \frac{3}{2r} + \frac{21}{16} - \frac{21}{8r} + \frac{3}{2r} - \frac{1}{4} = \frac{17}{16} + \frac{3}{8}r > 1.$

(4.57)

Hence, we cannot find $x_{1}$, $x_{2}$, $x_{3}$ and $x_{4}$ satisfying (4.56) and we are unable to deduce that ${\text{div}} (\sigma_{k} {v}_{k})$ lies in the dual space $(H^{1})^{*}$ even with the improved regularity ${v}_{k} \in L^{\frac{8}{7}}(0, T; {H}^{1})$.

---

5. Passing to the limit

From (4.44), (4.47), (4.19), (4.20), (4.21), (4.51), (4.52), (4.54), (4.55), we find that

$\begin{gathered} {\{ {\varphi _k}\} _{k \in \mathbb{N}}} {\text{bounded in }}{L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}), \hfill \\ {\{ {\partial _t}{\varphi _k}\} _{k \in \mathbb{N}}}, {\{ {\text{div}}({\varphi _k}{v_k})\} _{k \in \mathbb{N}}}{\text{ bounded in }}{L^{\frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\{ {\sigma _k}\} _{k \in \mathbb{N}}}{\text{ bounded in }}{L^\infty }(0, T;{L^2}) \cap {L^2}(0, T;{H^1}) \cap {L^2}(\Sigma ), \hfill \\ {\{ {\partial _t}{\sigma _k}\} _{k \in \mathbb{N}}}, {\{ {\text{div}}({\sigma _k}{v_k})\} _{k \in \mathbb{N}}}{\text{ bounded in }}{L^{\frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ {\{ {\mu _k}\} _{k \in \mathbb{N}}}{\text{ bounded in }}{L^2}(0, T;{H^1}), \hfill \\ {\{ {p_k}\} _{k \in \mathbb{N}}}{\text{ bounded in }}{L^{\frac{8}{5}}}(0, T;{H^1}) \cap {L^{\frac{8}{7}}}(0, T;{H^2}), \hfill \\ {\{ {v_k}\} _{k \in \mathbb{N}}}{\text{ bounded in }}{L^2}(0, T;{L^2}) \cap {L^{\frac{8}{7}}}(0, T;{H^1}). \hfill \\ \end{gathered}$

By standard compactness results (Banach-Alaoglu theorem and reflexive weak compactness theorem), and [28, x8, Cor. 4], and the compact embeddings in dimension 3 (see [1, Thm. 6.3] and [15, Thm. 11.2, p. 31])

$H^{j+1}(\Omega) = W^{j+1, 2}(\Omega) \subset \subset W^{j, q}(\Omega) \quad \forall j \geq 0, j \in \mathbb{Z}, \; 1 \leq q < 6,$

and the compact embedding $L^{2} \subset \subset (H^{1})^{*}$, we obtain, for a relabelled subsequence, the following weak/weak-* convergences:

$\begin{gathered} {\varphi _k} \to \varphi \quad {\text{ weakly - }}*\quad {\text{ in }}{L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, \frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\sigma _k} \to \sigma \quad {\text{ weakly - }}*\quad {\text{ in }}{L^2}(0, T;{H^1}) \cap {L^\infty }(0, T;{L^2}) \cap {L^2}(\Sigma ), \hfill \\ {\partial _t}{\sigma _k} \to {\partial _t}\sigma \quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ {\mu _k} \to \mu \quad {\text{ weakly }}\quad {\text{ in }}{L^2}(0, T;{H^1}), \hfill \\ {p_k} \to p\quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{H^1}) \cap {L^{\frac{8}{7}}}(0, T;{H^2}), \hfill \\ {v_k} \to v\quad {\text{ weakly }}\quad {\text{ in }}{L^2}(0, T;{L^2}) \cap {L^{\frac{8}{7}}}(0, T;{H^1}), \hfill \\ {\text{div}}({\varphi _k}{v_k}) \to \xi \quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\text{div}}({\sigma _k}{v_k}) \to \theta \quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ \end{gathered}$

and the following strong convergences:

$\begin{gathered} {\varphi _k} \to \varphi \quad {\text{ strongly }}\quad {\text{ in }}{C^0}([0, T];{L^r}) \cap {L^2}(0, T;{W^{2, r}}){\text{ and a}}{\text{.e}}{\text{. in }}Q, \hfill \\ {\sigma _k} \to \sigma \quad {\text{ strongly }}\quad {\text{ in }}{C^0}([0, T];{({H^1})^*}) \cap {L^2}(0, T;{L^r}){\text{ and a}}{\text{.e}}{\text{. in }}Q, \hfill \\ \end{gathered}$

for any $r \in [1, 6)$ and some functions $\xi \in L^{\frac{8}{5}}(0, T;(H^{1})^{*})$, $\theta \in L^{\frac{5}{4}}(0, T;(W^{1, 5})^{*})$.

For the rest of this section, we fix $1 \leq j \leq k$ and $\delta \in C^{\infty}_{c}(0, T)$. Then, we have $\delta(t) w_{j} \in C^{\infty}(0, T;H^{2})$. By continuity of $m(\cdot)$, we see that $m(\varphi_{k}) \to m(\varphi)$ a.e. in $Q$. Thanks to the boundedness of $m(\cdot)$, applying Lebesgue's dominated convergence theorem to ${(m({\varphi _k}) - m(\varphi ))^2}{\left| {\delta \nabla {w_j}} \right|^2}$ yields

${\left\| {m({\varphi _k})\delta \nabla {w_j} - m(\varphi )\delta \nabla {w_j}} \right\|_{{L^2}(Q)}} \to 0{\text{ as }}k \to \infty .$

Together with the weak convergence $\mu_{k} \rightharpoonup \mu$ in $L^{2}(0, T;H^{1})$, we obtain by the product of weak-strong convergence

$\int_0^T {\int_\Omega m } ({\varphi _k})\delta \nabla {w_j} \cdot \nabla {\mu _k}{\text{ dx dt}} \to \int_0^T {\int_\Omega m } (\varphi )\delta \nabla {w_j} \cdot \nabla \mu {\text{ dx dt as }}k \to \infty .$

Terms involving $n(\cdot)$ can be dealt with in a similar fashion. For the source term $\Gamma_{\varphi, k} = \Lambda_{\varphi}(\varphi_{k}, \sigma_{k}) - \Theta_{\varphi}(\varphi_{k}, \sigma_{k})\mu_{k}$, by the continuity and boundedness of $\Theta_{\varphi}$, the a.e. convergence of $\varphi_{k} \to \varphi$ and $\sigma_{k} \to \sigma$ in $Q$, we may apply Lebesgue's dominated convergence theorem to deduce that

$\int_0^T {\int_\Omega {} } {\left| {\delta {w_j}({\Theta _\varphi }({\varphi _k}, {\sigma _k}) - {\Theta _\varphi }(\varphi , \sigma ))} \right|^2}{\text{ dx dt}} \to 0{\text{ as }}k \to \infty ,$

that is, we obtain the strong convergence $\delta w_{j} \Theta_{\varphi}(\varphi_{k}, \sigma_{k}) \to \delta w_{j} \Theta_{\varphi}(\varphi, \sigma)$ in $L^{2}(Q)$. Hence, the weak convergence $\mu_{k} \rightharpoonup \mu$ in $L^{2}(0, T;H^{1})$ yields

$\int_0^T {\int_\Omega \delta } {w_j}{\Theta _\varphi }({\varphi _k}, {\sigma _k}){\mu _k}{\text{ dx dt}} \to \int_0^T {\int_\Omega \delta } {w_j}{\Theta _\varphi }(\varphi , \sigma )\mu {\text{ dx dt as }}k \to \infty .$

Meanwhile, by the triangle inequality $\left| {\left| a \right| - \left| b \right|} \right| \leq \left| {a - b} \right|$, and Hölder's inequality, we obtain

$\int_0^T {\int_\Omega {} } \left| {(\left| {{\varphi _k}} \right| - \left| \varphi \right|)(\delta {w_j})} \right|{\text{ dx dt}} \leq {\left\| {{\varphi _k} - \varphi } \right\|_{{L^2}(0, T;{L^2})}}{\left\| {\delta {w_j}} \right\|_{{L^2}(0, T;{L^2})}} \to 0$

and

$\int_0^T {\int_\Omega {} } \left| {(\left| {{\sigma _k}} \right| - \left| \sigma \right|)(\delta {w_j})} \right|{\text{ dx dt}} \leq {\left\| {{\sigma _k} - \sigma } \right\|_{{L^2}(0, T;{L^2})}}{\left\| {\delta {w_j}} \right\|_{{L^2}(0, T;{L^2})}} \to 0$

as $k \to \infty$. In particular, we have

$(1 + \left| {{\varphi _k}} \right| + \left| {{\sigma _k}} \right|)\left| {\delta {w_j}} \right| \to (1 + \left| \varphi \right| + \left| \sigma \right|)\left| {\delta {w_j}} \right|{\text{ strongly in }}{L^1}(Q){\text{ as }}k \to \infty .$

By the continuity of $\Lambda_{\varphi}$ we have

${\Lambda _\varphi }({\varphi _k}, {\sigma _k}) \to {\Lambda _\varphi }(\varphi , \sigma ){\text{ a}}{\text{.e}}{\text{. as }}k \to \infty , \quad \left| {{\Lambda _\varphi }({\varphi _k}, {\sigma _k})\delta {w_j}} \right| \leq {R_0}(1 + \left| {{\varphi _k}} \right| + \left| {{\sigma _k}} \right|)\left| {\delta {w_j}} \right|.$

which leads to

Then, the generalised Lebesgue dominated convergence theorem (see [27, Thm. 1.9, p. 89], or [2, Thm. 3.25, p. 60]) yields

$\Lambda_{\varphi}(\varphi_{k}, \sigma_{k}) \delta w_{j} \to \Lambda_{\varphi}(\varphi, \sigma) \delta w_{j} \text{ strongly in } L^{1}(Q) \text{ as } k \to \infty,$

which leads to

$\int_0^T {\int_\Omega {{\Gamma _\varphi }} } ({\varphi _k}, {\mu _k}, {\sigma _k})\delta {w_j}{\text{ dx dt}} \to \int_0^T {\int_\Omega {{\Gamma _\varphi }} } (\varphi , \mu , \sigma )\delta {w_j}{\text{ dx dt as }}k \to \infty .$

(5.1)

The same arguments can be applied for the source term $\mathcal{S}$ and for the derivative $\Psi'(\varphi)$ satisfying the linear growth condition (2.5). For potentials satisfying the growth condition (2.7), we refer to the argument in [18, x3.1.2].

To identify the limits $\xi$ and $\theta$ of the convection terms ${\text{div}} (\varphi_{k} {v}_{k})$ and ${\text{div}} (\sigma_{k} {v}_{k})$, respectively, we argue as follows. Since $\delta w_{j} \in C^{\infty}(0, T;H^{2}) \subset L^{\frac{8}{3}}(0, T;H^{1})$, by the weak convergence ${\text{div}} (\varphi_{k} {v}_{k}) \rightharpoonup \xi$ in $L^{\frac{8}{5}}(0, T;(H^{1})^{*})$, we have

$\int_0^T {\int_\Omega {{\text{div}}} } ({\varphi _k}{v_k})\delta {w_j}{\text{ dx dt}} \to \int_0^T {} {\left\langle {\xi, {w_j}} \right\rangle _{{H^1}, {{({H^1})}^*}}}\delta {\text{ dt as }}k \to \infty .$

Next, applying integrating by parts and by the boundary conditions ${v}_{k} \cdot {n} = 0$ on $\partial \Omega$ (see (3.8)), we see that

$\int_0^T {\int_\Omega {} } div({\varphi _k}{v_k})\delta {w_j}{\text{ dx dt}} = - \int_0^T {\int_\Omega \delta } {\varphi _k}{v_k} \cdot \nabla {w_j}{\text{ dx dt}}.$

(5.2)

Moreover, we claim that $\delta \varphi_{k} \nabla w_{j}$ converges strongly to $\delta \varphi \nabla w_{j}$ in $L^{2}(0, T;{L}^{2})$. Indeed, we compute

$\begin{gathered} \int_0^T {\int_\Omega {} } {\left| \delta \right|^2}{\left| {\nabla {w_j}} \right|^2}{\left| {{\varphi _k} - \varphi } \right|^2}{\text{dx dt}} \leq \int_0^T {} {\left| \delta \right|^2}\left\| {\nabla {w_j}} \right\|_{{L^6}}^2\left\| {{\varphi _k} - \varphi } \right\|_{{L^3}}^2{\text{dt}} \\ \leq \left\| {{w_j}} \right\|_{{H^2}}^2\left\| \delta \right\|_{{L^\infty }(0, T)}^2\left\| {{\varphi _k} - \varphi } \right\|_{{L^2}({L^3})}^2 \to 0 \\ \end{gathered}$

as $k \to \infty$ by the strong convergence $\varphi_{k} \to \varphi$ in $L^{2}(0, T;L^{r})$ for $r \in [1, 6)$. Together with the weak convergence ${v}_{k} \rightharpoonup {v}$ in $L^{2}(0, T;{L}^{2})$, when passing to the limit $k \to \infty$ in (5.2) we find that

$\int_0^T {} {\left\langle {\xi , {w_j}} \right\rangle _{{H^1}, {{({H^1})}^*}}}\delta {\text{ dt}} = - \int_0^T {\int_\Omega \delta } \varphi v \cdot \nabla {w_j}{\text{ dx dt}}.$

Applying integration by parts on the right-hand side shows that $\xi = {\text{div}} (\varphi {v})$ in the sense of distributions.

Now considering $\delta(t)w_{j}$ as an element in $L^{5}(0, T;W^{1, 5})$, a similar argument can be used to show $\theta = {\text{div}}(\sigma {v})$ in the sense of distributions using the strong convergence $\sigma_{k} \to \sigma$ in $L^{2}(0, T;L^{r})$ for $r \in [1, 6)$, the weak convergence ${v}_{k} \rightharpoonup {v}$ in $L^{2}(0, T;{L}^{2})$, and the weak convergence ${\text{div}}(\sigma_{k} {v}_{k}) \rightharpoonup \phi$ in $L^{\frac{5}{4}}(0, T;(W^{1, 5})^{*})$.

For the pressure and the velocity, we apply $-\Delta _{N}$ to both sides of (3.3) and test with $w_{j}$, then integrating by parts leads to

$\int_\Omega \nabla {p_k} \cdot \nabla {w_j}{\text{ dx}} = \int_\Omega {\frac{1}{K}} {\Gamma _v}{w_j} + ({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k} \cdot \nabla {w_j}{\text{ dx}}.$

Multiplying by $\delta(t)$, integrating in time and passing to the limit $k \to \infty$, keeping in mind the weak convergences $p_{k} \rightharpoonup p$ in $L^{\frac{8}{5}}(0, T;H^{1})$, $\mu_{k} \rightharpoonup \mu$ in $L^{2}(0, T;H^{1})$, $\sigma_{k} \rightharpoonup \sigma$ in $L^{2}(0, T;H^{1})$, and the strong convergence $\varphi_{k} \to \varphi$ in $L^{2}(0, T;W^{2, r})$ for $r \in [1, 6)$ leads to

$\int_0^T {\int_\Omega \delta } (t)\nabla p \cdot \nabla {w_j}{\text{ dx dt}} = \int_0^T {\int_\Omega \delta } (t)\left( {\frac{1}{K}{\Gamma _v}{w_j} + (\mu + \chi \sigma )\nabla \varphi \cdot \nabla {w_j}} \right){\text{dx dt}}.$

(5.3)

Here we used that $w_{j} \in H^{2}$, and

$\begin{gathered} \int_0^T {\int_\Omega {} } {\left| \delta \right|^2}{\left| {\nabla {\varphi _k} - \nabla \varphi } \right|^2}{\left| {\nabla {w_j}} \right|^2}{\text{dx dt}} \hfill \\ \leq \left\| \delta \right\|_{{L^\infty }(0, T)}^2\left\| {{w_j}} \right\|_{{W^{1, 6}}}^2\left\| {{\varphi _k} - \varphi } \right\|_{{L^2}({W^{1, 3}})}^2 \to 0{\text{ as }}k \to \infty , \hfill \\ \end{gathered}$

(5.4)

to deduce that $\delta \nabla \varphi_{k} \cdot \nabla w_{j} \to \delta \nabla \varphi \cdot \nabla w_{j}$ in $L^{2}(0, T;L^{2})$. Fix $1 \leq j_{1}, j_{2}, j_{3} \leq k$, and define ${\zeta} = (w_{j_{1}}, w_{j_{2}}, w_{j_{3}})^{\top}$. Then, we can consider $\delta(t) {\zeta}$ as an element in $L^{\frac{8}{3}}(0, T;{L}^{2}) \subset L^{2}(0, T;{L}^{2})$. Taking the scalar product of (3.4) with $\delta {\zeta}$, integrating over $\Omega$ and in time from $0$ to $T$ leads to

$\int_0^T {\int_\Omega \delta } (t)({v_k} + K\nabla {p_k}) \cdot \zeta {\text{ dx dt}} = \int_0^T {\int_\Omega \delta } K({\mu _k} + \chi {\sigma _k})\nabla {\varphi _k} \cdot \zeta {\text{ dx dt}}.$

(5.5)

By the weak convergences ${v}_{k} \rightharpoonup {v}$ in $L^{2}(0, T;{L}^{2})$, $\mu_{k} \rightharpoonup \mu$ in $L^{2}(0, T;H^{1})$, $\sigma_{k} \rightharpoonup \sigma$ in $L^{2}(0, T;H^{1})$, $\nabla p_{k} \rightharpoonup \nabla p$ in $L^{\frac{8}{5}}(0, T;{L}^{2})$, and the strong convergence $\delta \nabla \varphi_{k} \cdot {\zeta} \to \delta \nabla \varphi \cdot {\zeta}$ in $L^{2}(0, T; L^{2})$ (which is proved in a similar manner as (5.4)), we find that passing to the limit in (5.5) yields

$\int_0^T {\int_\Omega \delta } (t)(v + K\nabla p) \cdot \zeta {\text{ dx dt}} = \int_0^T {\int_\Omega \delta } (t)K(\mu + \chi \sigma )\nabla \varphi \cdot \zeta {\text{ dx dt}}.$

(5.6)

Then, multiplying (3.2) with $\delta \in C^{\infty}_{c}(0, T)$, integrating with respect to time from $0$ to $T$, and passing to the limit $k \to \infty$, we obtain

$\begin{gathered} \int_0^T \delta (t){\left\langle {{\partial _t}\varphi , {w_j}} \right\rangle _{{H^1}, {{({H^1})}^*}}}{\text{dt}} = \int_0^T {\int_\Omega \delta } (t)\left( { - m(\varphi )\nabla \mu \cdot \nabla {w_j} + {\Gamma _\varphi }{w_j} + \varphi v \cdot \nabla {w_j}} \right){\text{dx dt}}, \hfill \\ \int_0^T {\int_\Omega \delta } (t)\mu {w_j}{\text{dx dt}} = \int_0^T {\int_\Omega \delta } (t)\left( {A\Psi '(\varphi ){w_j} + B\nabla \varphi \cdot \nabla {w_j} - \chi \sigma {w_j}} \right){\text{dx dt}}, \hfill \\ \int_0^T \delta (t){\left\langle {{\partial _t}\sigma , {w_j}} \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}}{\text{dt}} = \int_0^T {\int_\Omega \delta } (t)\left( { - n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla {w_j} - \mathcal{S}{w_j}} \right){\text{dx dt}} \hfill \\ + \int_0^T \delta (t)\left( {\int_\Omega \sigma v \cdot \nabla {w_j}{\text{ dx}} + \int_\Gamma b ({\sigma _\infty } - \sigma ){w_j}{\text{ d}}{H^{d - 1}}} \right){\text{dt}}. \hfill \\ \end{gathered}$

Since the above, (5.3) and (5.6) hold for all $\delta \in C^{\infty}_{c}(0, T)$, we infer that $\{\varphi, \mu, \sigma, p, {v}\}$ satisfies (2.12) with $\zeta = \phi = w_{j}$ for a.e. $t \in (0, T)$ and for all $j \geq 1$. As $\{ w_{j}\}_{j \in \mathbb{N}}$ is a basis for $H^{2}_{N}$, and $H^{2}_{N}$ is dense in both $H^{1}$ and $W^{1, 5}$ (see Section 3), we see that $\{ \varphi, \mu, \sigma, p, {v}\}$ satisfy (2.12a), (2.12b), (2.12d) for all $\zeta \in H^{1}$, (2.12c) for all $\phi \in W^{1, 5}$, and (2.12e) for all ${\zeta} \in {L}^{2}$.

Attainment of initial conditions. It remains to show that $\varphi$ and $\sigma$ attain their corresponding initial conditions. Strong convergence of $\varphi_{k}$ to $\varphi$ in $C^{0}([0, T];L^{2})$, and the fact that $\varphi_{k}(0) \to \varphi_{0}$ in $L^{2}$ imply that $\varphi(0) = \varphi_{0}$. Meanwhile, as the limit function $\sigma$ belongs to the function space $C^{0}([0, T]; (H^{1})^{*})$, we see that $\sigma(0):= \sigma(\cdot, 0)$ makes sense as an element of $(H^{1})^{*}$. Let $\zeta \in H^{1}$ be arbitrary, then by the strong convergence $\sigma_{k} \to \sigma$ in $C^{0}([0, T];(H^{1})^{*})$ we see that

${\left\langle {{\sigma _k}(0), \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} \to {\left\langle {\sigma (0), \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}}{\text{ as }}k \to \infty .$

On the other hand, by (3.7), we have $\sigma_{k}(0) \to \sigma_{0}$ in $L^{2}$. This yields

${\left\langle {{\sigma _0}, \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \mathop {\lim }\limits_{k \to \infty } {\left\langle {{\sigma _k}(0), \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = {\left\langle {\sigma (0), \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}}.$

Energy inequality. For the energy inequality (2.14) we employ the weak/weak-* lower semicontinuity of the norms and dual norms to (4.44), (4.47), (4.19), (4.20), (4.21), (4.51), (4.52), (4.54), and (4.55).

---

6. Asymptotic limits

---

6.1. Limit of vanishing permeability

For $K, b \in (0, 1]$ let $(\varphi^{K}, \mu^{K}, \sigma^{K}, {v}^{K}, p^{K})$ denote a weak solution to (1.1)-(1.2) with $\Gamma_{{v}} = 0$, obtain from Theorem 2.1. By we deduce that, for a relabelled subsequence as $b \to 0$ and $K \to 0$, the following weak/weak-* convergences:

$\begin{gathered} {\varphi ^K} \to \varphi \quad {\text{ weakly - }}*\quad {\text{ in }}{L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, \frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\sigma ^K} \to \sigma \quad {\text{ weakly - }}*\quad {\text{ in }}{L^2}(0, T;{H^1}) \cap {L^\infty }(0, T;{L^2}) \cap {W^{1, \frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ {\mu ^K} \to \mu \quad {\text{ weakly }}\quad {\text{ in }}{L^2}(0, T;{H^1}), \hfill \\ {p^K} \to p\quad {\text{ weakly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{H^1}) \cap {L^{\frac{8}{7}}}(0, T;{H^2}), \hfill \\ \end{gathered}$

and the following strong convergences:

$\begin{gathered} {\varphi ^K} \to \varphi \quad {\text{ strongly }}\quad {\text{ in }}{C^0}([0, T];{L^r}) \cap {L^2}(0, T;{W^{2, r}}){\text{ and a}}{\text{.e}}{\text{. in }}Q, \hfill \\ {\sigma ^K} \to \sigma \quad {\text{ strongly }}\quad {\text{ in }}{C^0}([0, T];{({H^1})^*}) \cap {L^2}(0, T;{L^r}){\text{ and a}}{\text{.e}}{\text{. in }}Q, \hfill \\ {v^K} \to 0\quad {\text{ strongly }}\quad {\text{ in }}{L^2}(0, T;{L^2}) \cap {L^{\frac{8}{7}}}(0, T;{H^1}), \hfill \\ {\text{div}}({\varphi ^K}{v^K}) \to 0\quad {\text{ strongly }}\quad {\text{ in }}{L^{\frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\text{div}}({\sigma ^K}{v^K}) \to 0\quad {\text{ strongly }}\quad {\text{ in }}{L^{\frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ \end{gathered}$

for any $r \in [1, 6)$. The strong convergence of the velocity and the convection terms to zero follows from (2.14). Upon multiplying (2.12) by $\delta \in C^{\infty}_{c}(0, T)$ and passing to the limit $b, K \to 0$, we obtain that the limit functions $(\varphi, \mu, \sigma, p)$ satisfy

${\left\langle {{\partial _t}\varphi , \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \int_\Omega - m(\varphi )\nabla \mu \cdot \nabla \zeta + {\Gamma _\varphi }(\varphi , \mu , \sigma )\zeta {\text{ dx}},$

(6.1a)

$\int_\Omega \mu \zeta {\text{ dx}} = \int_\Omega A \Psi '(\varphi )\zeta + B\nabla \varphi \cdot \nabla \zeta - \chi \sigma \zeta {\text{ dx}},$

(6.1b)

${\left\langle {{\partial _t}\sigma , \phi } \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}} = \int_\Omega - n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla \phi - \mathcal{S}(\varphi , \mu , \sigma )\phi {\text{ dx}}$

(6.1c)

$\int_\Omega \nabla p \cdot \nabla \zeta {\text{ dx}} = \int_\Omega {(\mu + \chi \sigma )} \nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(6.1d)

for all $\zeta \in H^{1}$ and $\phi \in W^{1, 5}$ and a.e. $t \in (0, T)$.

Note that substituting any $\zeta \in L^{2}(0, T;H^{1})$ into (6.1a), integrating in time from $0$ to $T$, using Hölder's inequality and the linear growth condition for $\Gamma_{\varphi}$ leads to the deduction that $\partial _{t} \varphi \in L^{2}(0, T;(H^{1})^{*})$. To show that $\partial _{t} \sigma \in L^{2}(0, T;(H^{1})^{*})$ we argue as follows. For any $\xi \in L^{2}(0, T;H^{1})$, we can define

$F(\xi ): = \int_0^T {\int_\Omega - } n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla \xi - \mathcal{S}(\varphi , \mu , \sigma )\xi {\text{ dx dt}}.$

By Hölder's inequality and the growth condition on $\mathcal{S}$, we see that $F \in L^{2}(0, T;(H^{1})^{*})$. It is known that the set of functions that are finite linear combinations of $C^{1}_{c}(0, T) \cdot H^{2}_{N}(\Omega) := \{ \delta(t) \phi(x) : \delta \in C^{1}_{c}(0, T), \; \phi \in H^{2}_{N}(\Omega) \}$ is dense in $C^{1}_{c}(0, T;H^{1})$ (see for instance [26, p. 384], and in fact this is what we use in Section 5). Let $\zeta \in C^{1}_{c}(0, T;H^{1})$ and let $\{\zeta^{n}\}_{n \in \mathbb{N}}$ denote a sequence of functions of the above form such that $\zeta^{n} \to \zeta$ in $C^{1}_{c}(0, T;H^{1})$ as $n \to \infty$. Then, substituting $\phi = \zeta^{n}$ in (6.1c), integrating over $t$ from $0$ to $T$, and passing to the limit $n \to \infty$ yields

$\mathop {\lim }\limits_{n \to \infty } \int_0^T {} {\left\langle {{\partial _t}\sigma , {\zeta ^n}} \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}}{\text{ dt}} = \mathop {\lim }\limits_{n \to \infty } F({\zeta ^n}) = F(\zeta ).$

Moreover, by the definition of the weak time derivative, we have

$\int_0^T {} {\left\langle {{\partial _t}\sigma , {\zeta ^n}} \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}}{\text{ dt}} = - \int_0^T {\int_\Omega \sigma } {\partial _t}{\zeta ^n}{\text{ dx dt}} \to - \int_0^T {\int_\Omega \sigma } {\partial _t}\zeta {\text{ dx dt as }}n \to \infty .$

Hence, we obtain

$- \int_0^T {\int_\Omega \sigma } {\partial _t}\zeta {\text{ dx dt}} = F(\zeta ) = \int_0^T {\int_\Omega - } n(\varphi )(D\nabla \sigma - \chi \nabla \varphi ) \cdot \nabla \zeta - \mathcal{S}(\varphi , \mu , \sigma )\zeta {\text{ dx dt}}$

for all $\zeta \in C^{1}_{c}(0, T;H^{1})$. This implies that the weak time derivative $\partial _{t} \sigma$ satisfies

$\int_0^T {} {\partial _t}\sigma {\zeta _{{H^1}, {{({H^1})}^*}}}{\text{ dt}} = F(\zeta )\quad \forall \zeta \in C_c^1(0, T;{H^1}),$

and as $F$ belongs to $L^{2}(0, T;(H^{1})^{*})$, we find that $\partial _{t} \sigma$ also belongs to $L^{2}(0, T;(H^{1})^{*})$. Furthermore, due to the improved regularity $\partial _{t} \sigma \in L^{2}(0, T;(H^{1})^{*})$, we use the continuous embedding

$L^{2}(0, T;H^{1}) \cap H^{1}(0, T;(H^{1})^{*}) \subset C^{0}([0, T];L^{2})$

to deduce that $\sigma(0) = \sigma_{0}$.

---

6.2. Limit of vanishing chemotaxis

For $\chi, b \in (0, 1]$, let $(\varphi^{\chi}, \mu^{\chi}, \sigma^{\chi}, {v}^{\chi}, p^{\chi})$ denote a weak solution to (1.1)-(1.2) obtain from Theorem 2.1. By (2.14) we deduce that, for a relabelled subsequence as $b \to 0$ and $\chi \to 0$, the following weak/weak-* convergences:

$\begin{gathered} {\varphi ^\chi } \to \varphi {\mkern 1mu} {\text{ weakly - }}*{\mkern 1mu} {\text{ in }}{L^\infty }(0, T;{H^1}) \cap {L^2}(0, T;{H^3}) \cap {W^{1, \frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\sigma ^\chi } \to \sigma {\mkern 1mu} {\text{ weakly - }}*{\mkern 1mu} {\text{ in }}{L^2}(0, T;{H^1}) \cap {L^\infty }(0, T;{L^2}) \cap {W^{1, \frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ {\mu ^\chi } \to \mu {\mkern 1mu} {\text{ weakly }}{\mkern 1mu} {\text{ in }}{L^2}(0, T;{H^1}), \hfill \\ {p^\chi } \to p{\mkern 1mu} {\text{ weakly }}{\mkern 1mu} {\text{ in }}{L^{\frac{8}{5}}}(0, T;{H^1}) \cap {L^{\frac{8}{7}}}(0, T;{H^2}), \hfill \\ {v^\chi } \to v{\mkern 1mu} {\text{ weakly }}{\mkern 1mu} {\text{ in }}{L^2}(0, T;{L^2}) \cap {L^{\frac{8}{7}}}(0, T;{H^1}), \hfill \\ {\text{div}}({\varphi ^\chi }{v^\chi }) \to {\text{div}}(\varphi v){\mkern 1mu} {\text{ weakly }}{\mkern 1mu} {\text{ in }}{L^{\frac{8}{5}}}(0, T;{({H^1})^*}), \hfill \\ {\text{div}}({\sigma ^\chi }{v^\chi }) \to {\text{div}}(\sigma v){\mkern 1mu} {\text{ weakly }}{\mkern 1mu} {\text{ in }}{L^{\frac{5}{4}}}(0, T;{({W^{1, 5}})^*}), \hfill \\ \end{gathered}$

and the following strong convergences:

$\begin{gathered} {\varphi ^\chi } \to \varphi \quad {\text{ strongly }}\quad {\text{ in }}{C^0}([0, T];{L^r}) \cap {L^2}(0, T;{W^{2, r}}){\text{ and a}}{\text{.e}}{\text{. in }}Q, \hfill \\ {\sigma ^\chi } \to \sigma \quad {\text{ strongly }}\quad {\text{ in }}{C^0}([0, T];{({H^1})^*}) \cap {L^2}(0, T;{L^r}){\text{ and a}}{\text{.e}}{\text{. in }}Q, \hfill \\ \end{gathered}$

for any $r \in [1, 6)$. For any $\delta \in C^{\infty}_{c}(0, T)$ and $\zeta \in H^{1}$, we have

$\begin{gathered} \left| {\int_0^T {\int_\Omega \delta } \chi {\sigma ^\chi }\zeta {\text{dx dt}}} \right| \leq \chi {\left\| {{\sigma ^\chi }} \right\|_{{L^2}({L^2})}}{\left\| \zeta \right\|_{{L^2}}}{\left\| \delta \right\|_{{L^2}(0, T)}} \to 0, \hfill \\ \left| {\int_0^T {\int_\Omega \delta } n({\varphi ^\chi })\chi \nabla {\varphi ^\chi } \cdot \nabla \zeta {\text{dx dt}}} \right| \leq {n_1}\chi {\left\| {\nabla {\varphi ^\chi }} \right\|_{{L^2}({L^2})}}{\left\| {\nabla \zeta } \right\|_{{L^2}}}{\left\| \delta \right\|_{{L^2}(0, T)}} \to 0, \hfill \\ \left| {\int_0^T {\int_\Omega \delta } \chi {\sigma ^\chi }\nabla {\varphi ^\chi } \cdot \nabla \zeta {\text{ dx dt}}} \right| \leq \chi {\left\| {\nabla \zeta } \right\|_{{L^2}}}{\left\| {{\sigma ^\chi }} \right\|_{{L^2}({L^4})}}{\left\| {\nabla {\varphi ^\chi }} \right\|_{{L^2}({L^4})}}{\left\| \delta \right\|_{{L^\infty }(0, T)}} \to 0, \hfill \\ \end{gathered}$

as $\chi \to 0$. Thus, multiplying (2.12) with $\delta \in C^{\infty}_{c}(0, T)$, and then passing to the limit $b, \chi \to 0$, we see that $(\varphi, \mu, \sigma, {v}, p)$ satisfies

${\left\langle {{\partial _t}\varphi , \zeta } \right\rangle _{{H^1}, {{({H^1})}^*}}} = \int_\Omega - m(\varphi )\nabla \mu \cdot \nabla \zeta + {\Gamma _\varphi }(\varphi , \mu , \sigma )\zeta + \varphi v \cdot \nabla \zeta {\text{ dx}},$

(6.2a)

$\int_\Omega \mu \zeta {\text{ dx}} = \int_\Omega A \Psi '(\varphi )\zeta + B\nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(6.2b)

${\left\langle {{\partial _t}\sigma , \phi } \right\rangle _{{W^{1, 5}}, {{({W^{1, 5}})}^*}}} = \int_\Omega - n(\varphi )D\nabla \sigma \cdot \nabla \phi - \mathcal{S}(\varphi , \mu , \sigma )\phi + \sigma v \cdot \nabla \phi {\text{ dx}},$

(6.2c)

$\int_\Omega \nabla p \cdot \nabla \zeta {\text{ dx}} = \int_\Omega {\frac{1}{K}} {\Gamma _v}\zeta + \mu \nabla \varphi \cdot \nabla \zeta {\text{ dx}},$

(6.2d)

$\int_\Omega {} v \cdot \zeta {\text{ dx}} = \int_\Omega - K(\nabla p - \mu \nabla \varphi ) \cdot \zeta {\text{ dx}},$

(6.2e)

for all $\zeta \in H^{1}$, $\phi \in W^{1, 5}$, ${\zeta} \in {L}^{2}$ and a.e. $t \in (0, T)$.

---

7. Existence in two dimensions

We first derive an analogous result to Lemma 4.2 for two dimensions.

Lemma 7.1. Let $\Omega \subset \mathbb{R}^{2}$ be a bounded domain with $C^{3}$-boundary. Let $\varphi \in L^{\infty}(0, T;H^{1}) \cap L^{2}(0, T;H^{2}_{N} \cap H^{3})$, $\sigma \in L^{2}(0, T;H^{1})$, $\mu \in L^{2}(0, T;H^{1})$, the source term $\Gamma_{{v}} \in L^{2}(0, T;L^{2}_{0})$, and the function $p$ satisfying (4.7). Then,

$p \in {L^k}(0, T;{H^1}) \cap {L^q}(0, T;{H^2}), \quad {\text{ }}v \in {L^q}(0, T;{H^1}),$

for any

$1 \leq k < 2, \quad 1 \leq q < \frac{4}{3}.$

Proof. We estimate (4.23) differently than in the proof of Lemma 4.2. By Hölder’s inequality for any $1 \leq s < \infty$ we have

${\left\| {(\mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}} \leq {\left\| {\mu + \chi \sigma } \right\|_{{L^{2s}}}}{\left\| {\nabla \varphi } \right\|_{{L^{\frac{{2s}}{{s - 1}}}}}}.$

Then, by the Gagliardo--Nirenberg inequality (1.10) with $p = \frac{2s}{s-1}$, $j = 0$, $r = 2$, $m = 2$, $d = 2$, $q = 2$, and $\alpha = \frac{1}{2} - \frac{s-1}{2s} = \frac{1}{2s}$, we find that

${\left\| {\nabla \varphi } \right\|_{{L^{\frac{{2s}}{{s - 1}}}}}} \leq C\left\| {\nabla \varphi } \right\|_{{H^2}}^{\frac{1}{{2s}}}\left\| {\nabla \varphi } \right\|_{{L^2}}^{1 - \frac{1}{{2s}}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{1}{{2s}}}\left\| \varphi \right\|_{{H^1}}^{1 - \frac{1}{{2s}}}.$

Then, by Hölder's inequality and the Sobolev embedding $H^{1} \subset L^{r}$ for $1 \leq r < \infty$ in two dimensions, we have for $w, y \geq 1$,

$\begin{gathered} \int_0^T {} \left\| {(\mu + \chi \sigma )\nabla \varphi } \right\|_{{L^2}}^w{\text{dt}} \leq C\int_0^T {} \left\| {\mu + \chi \sigma } \right\|_{{L^{2s}}}^w\left\| {\nabla \varphi } \right\|_{{L^{\frac{{2s}}{{s - 1}}}}}^w{\text{dt}} \hfill \\ \quad \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{w\frac{{2s - 1}}{{2s}}}\left\| {\mu + \chi \sigma } \right\|_{{L^{wy}}({H^1})}^w\left\| \varphi \right\|_{{L^{\frac{w}{{2s}}\frac{y}{{y - 1}}}}({H^3})}^{\frac{w}{{2s}}}. \hfill \\ \end{gathered}$

As $\mu, \sigma$ belong to $L^{2}(0, T;H^{1})$ and $\varphi$ belongs to $L^{2}(0, T;H^{3}) \cap L^{\infty}(0, T;H^{1})$, we need

$wy = 2, \quad \frac{w}{{2s}}\frac{y}{{y - 1}} = 2 \Rightarrow y = \frac{{2s + 1}}{{2s}}, \quad w = \frac{{4s}}{{1 + 2s}}.$

Since $w = \frac{4s}{1+2s} < 2$ for all $s \in [1, \infty)$, and $\Gamma_{{v}} \in L^{2}(0, T;L^{2}_{0})$, the computations in the proof of Lemma 4.2 yields that

$p \in L^{k}(0, T;H^{1}) \quad \text{ for } 1 \leq k < 2$

Next, we see that

$\begin{gathered} {\left\| {{\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right\|_{{L^2}}} \leq {\left\| {(\mu + \chi \sigma )\Delta \varphi } \right\|_{{L^2}}} + {\left\| {\nabla (\mu + \chi \sigma ) \cdot \nabla \varphi } \right\|_{{L^2}}} \\ \leq {\left\| {\mu + \chi \sigma } \right\|_{{L^{2s}}}}{\left\| {\Delta \varphi } \right\|_{{L^{\frac{{2s}}{{s - 1}}}}}} + {\left\| {\nabla (\mu + \chi \sigma )} \right\|_{{L^2}}}{\left\| {\nabla \varphi } \right\|_{{L^\infty }}}. \\ \end{gathered}$

By the Gagliardo--Nirenberg inequality (1.10) with $p = \infty$, $j = 0$, $r = 2$, $m = 2$, $d = 2$, $q = 2$ and $\alpha = \frac{1}{2}$, we have

${\left\| {\nabla \varphi } \right\|_{{L^\infty }}} \leq C\left\| {\nabla \varphi } \right\|_{{H^2}}^{\frac{1}{2}}\left\| {\nabla \varphi } \right\|_{{L^2}}^{\frac{1}{2}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{1}{2}}\left\| \varphi \right\|_{{H^1}}^{\frac{1}{2}},$

(7.1)

and with $p = \frac{2s}{s-1}$, $j = 1$, $r = 2$, $m = 2$, $d = 2$, $q = 2$ and $\alpha = \frac{s+1}{2s} \in (\frac{1}{2}, 1]$ for $s \in [1, \infty)$, we have

${\left\| {\Delta \varphi } \right\|_{{L^{\frac{{2s}}{{s - 1}}}}}} \leq C\left\| {\nabla \varphi } \right\|_{{H^2}}^{\frac{{s + 1}}{{2s}}}\left\| {\nabla \varphi } \right\|_{{L^2}}^{\frac{{s - 1}}{{2s}}} \leq C\left\| \varphi \right\|_{{H^3}}^{\frac{{s + 1}}{{2s}}}\left\| \varphi \right\|_{{H^1}}^{\frac{{s - 1}}{{2s}}}.$

(7.2)

Hence, for $w, y, z \geq 1$, we find that

$\begin{gathered} \int_0^T {} \left\| {{\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right\|_{{L^2}}^w{\text{dt}} \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{w}{2}}\left\| {\mu + \chi \sigma } \right\|_{{L^{wz}}({H^1})}^w\left\| \varphi \right\|_{{L^{\frac{w}{2}\frac{z}{{z - 1}}}}({H^3})}^{\frac{w}{2}} \\ + C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{w\frac{{s - 1}}{{2s}}}\left\| {\mu + \chi \sigma } \right\|_{{L^{wy}}({H^1})}^w\left\| \varphi \right\|_{{L^{w\frac{{s + 1}}{{2s}}\frac{y}{{y - 1}}}}({H^3})}^{w\frac{{s + 1}}{{2s}}}. \\ \end{gathered}$

Since $\frac{s+1}{2s} \leq 1$ for all $s \in [1, \infty)$, we require

$wy = 2, \quad \frac{{wy(s + 1)}}{{2s(y - 1)}} = 2 \Rightarrow y = \frac{{3s + 1}}{{2s}}, \quad w = \frac{{4s}}{{1 + 3s}}.$

We choose $z = \frac{3s+1}{2s} \in (\frac{3}{2}, 2]$ so that

$wz = 2, \quad \frac{w}{2} \frac{z}{z-1} = \frac{2s}{1 + 3s} \frac{3s+1}{s+1} = \frac{2s}{s+1} \in [1, 2),$

and thus we obtain

$\begin{gathered} \int_0^T {} \left\| {{\text{div}}((\mu + \chi \sigma )\nabla \varphi )} \right\|_{{L^2}}^{\frac{{4s}}{{1 + 3s}}}{\text{dt}} \leq C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{{2s}}{{1 + 3s}}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{{4s}}{{1 + 3s}}}\left\| \varphi \right\|_{{L^{\frac{{2s}}{{s + 1}}}}({H^3})}^{\frac{{2s}}{{1 + 3s}}} \\ + C\left\| \varphi \right\|_{{L^\infty }({H^1})}^{\frac{{2s - 2}}{{1 + 3s}}}\left\| {\mu + \chi \sigma } \right\|_{{L^2}({H^1})}^{\frac{{4s}}{{1 + 3s}}}\left\| \varphi \right\|_{{L^2}({H^3})}^{\frac{{2s + 2}}{{1 + 3s}}}. \\ \end{gathered}$

From (4.27) and using the fact that $\frac{4s}{1+3s} < \frac{4s}{1+2s}$ for all $s \in [1, \infty)$, we see that

$p \in L^{q}(0, T;H^{2}) \quad \text{ for } 1 \leq q < \frac{4}{3}.$

Similarly, from (4.29), (7.1) and (7.2), we obtain for fixed $1 \leq i, j \leq 2$, and any $s \in [1, \infty)$,

$\begin{gathered} {\left\| {{D_i}{v_j}} \right\|_{{L^2}}} = K{\left\| {{D_i}{D_j}p - ({D_i}(\mu + \chi \sigma ){D_j}\varphi - (\mu + \chi \sigma ){D_i}{D_j}\varphi } \right\|_{{L^2}}} \\ \leq K\left( {{{\left\| p \right\|}_{{H^2}}} + {{\left\| {\nabla (\mu + \chi \sigma )} \right\|}_{{L^2}}}{{\left\| {\nabla \varphi } \right\|}_{{L^\infty }}} + {{\left\| {\mu + \chi \sigma } \right\|}_{{L^{2s}}}}{{\left\| {{D^2}\varphi } \right\|}_{{L^{\frac{{2s}}{{s - 1}}}}}}} \right) \\ \leq K\left( {{{\left\| p \right\|}_{{H^2}}} + C{{\left\| {\mu + \chi \sigma } \right\|}_{{H^1}}}\left( {\left\| \varphi \right\|_{{H^3}}^{\frac{1}{2}}\left\| \varphi \right\|_{{H^1}}^{\frac{1}{2}} + \left\| \varphi \right\|_{{H^3}}^{\frac{{s + 1}}{{2s}}}\left\| \varphi \right\|_{{H^1}}^{\frac{{s - 1}}{{2s}}}} \right)} \right). \\ \end{gathered}$

(7.3)

Then, a similar calculation shows that the right-hand side is bounded in $L^{\frac{4s}{1+3s}}(0, T)$, which in turn implies that

$v \in {L^q}(0, T;{H^1}){\text{ for }}1 \leq q < \frac{4}{3}.$

By the above new estimates we can show that ${\text{div}}(\varphi {v})$ and $\partial _{t} \varphi$ have improved temporal regularity, and that ${\text{div}} (\sigma {v})$ and $\partial _{t}\sigma$ belong to the dual space $(H^{1})^{*}$.

Lemma 7.2. For dimension $d = 2$, let $(\varphi_{k}, \mu_{k}, \sigma_{k}, p_{k}, {v}_{k})$ denote the Galerkin ansatz from Section 3 satisfying (4.44). Then, it holds that for $\frac{4}{3} \leq w < 2$ and $1 < r < \frac{8}{7}$,

$\begin{gathered} {\left\| {{\text{div}}({\varphi _k}{v_k})} \right\|_{{L^w}({{({H^1})}^*})}} + {\left\| {{\text{div}}({\sigma _k}{v_k})} \right\|_{{L^r}({{({H^1})}^*})}} \leq {K^{\frac{1}{2}}}\mathcal{E}, \hfill \\ {\left\| {{\partial _t}{\varphi _k}} \right\|_{{L^w}({{({H^1})}^*})}} + {\left\| {{\partial _t}{\sigma _k}} \right\|_{{L^r}({{({H^1})}^*})}} \leq \mathcal{E}\left( {1 + {K^{\frac{1}{2}}}} \right), \hfill \\ \end{gathered}$

where $\mathcal{E}$ denotes positive constants that are uniformly bounded for $b, \chi \in (0, 1]$ and are also uniformly bounded for $K \in (0, 1]$ when $\Gamma_{{v}} = 0$.

Proof. The assertions for $\partial _{t}\varphi_{k}$ and $\partial _{t}\sigma_{k}$ will follow via similar arguments in Section 4.3 once we establish the assertion for the convection terms. In dimension $d = 2$, we have the embedding $L^{2}(0, T;H^{1}) \cap L^{\infty}(0, T;L^{2}) \subset L^{4}(Q)$, and by the Gagliardo--Nirenberg inequality (1.10) with $p = 4$, $j = 0$, $r = 2$, $d = 2$, $m = 1$, $q = 2$ and $\alpha = \frac{1}{2}$,

${\left\| f \right\|_{{L^4}}} \leq C\left\| f \right\|_{{H^1}}^{\frac{1}{2}}\left\| f \right\|_{{L^2}}^{\frac{1}{2}}.$

Consider an arbitrary $\zeta \in L^{s}(0, T;H^{1})$ for some $s \geq 1$ yet to be determined. Then, we compute that

$\begin{gathered} \left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{dx dt}}} \right| \leq \int_0^T {} {\left\| {{\sigma _k}} \right\|_{{L^4}}}{\left\| {{v_k}} \right\|_{{L^4}}}{\left\| {\nabla \zeta } \right\|_{{L^2}}}{\text{dt}} \\ \leq C{\left\| {{\sigma _k}} \right\|_{{L^4}(Q)}}{\left( {\int_0^T {} \left\| {{v_k}} \right\|_{{H^1}}^{\frac{2}{3}}\left\| {{v_k}} \right\|_{{L^2}}^{\frac{2}{3}}\left\| \zeta \right\|_{{H^1}}^{\frac{4}{3}}{\text{dt}}} \right)^{\frac{3}{4}}} \\ \leq C{\left\| {{\sigma _k}} \right\|_{{L^4}(Q)}}\left\| {{v_k}} \right\|_{{L^{\frac{2}{3}{x_1}}}({H^1})}^{\frac{1}{2}}\left\| {{v_k}} \right\|_{{L^{\frac{2}{3}{x_2}}}({L^2})}^{\frac{1}{2}}{\left\| \zeta \right\|_{{L^{\frac{4}{3}{x_3}}}({H^1})}}, \\ \end{gathered}$

where $x_{1}, x_{2}, x_{3} \geq 1$ satisfy

$\frac{1}{{{x_1}}} + \frac{1}{{{x_2}}} + \frac{1}{{{x_3}}} = 1, \quad \frac{2}{3}{x_1} < \frac{4}{3}, \quad \frac{2}{3}{x_2} \leq 2 \Rightarrow {x_3} > 6.$

Then, from (4.44) and (7.3), it holds that

$\left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{ dx dt}}} \right| \leq \mathcal{E}{K^{\frac{1}{2}}}{\left\| \zeta \right\|_{{L^s}({H^1})}}\quad {\text{ for }}s = \frac{4}{3}{x_3} > 8,$

that is, $\{{\text{div}} (\sigma_{k} {v}_{k})\}_{k \in \mathbb{N}}$ is uniformly bounded in the dual space of $L^{s}(0, T;H^{1})$ for $s > 8$. Similarly, by the Gagliardo--Nirenberg inequality (1.10) with $p = \infty$, $j = 0$, $r = 2$, $d = 2$, $m = 3$, $q \in [1, \infty)$ and $\alpha = \frac{1}{q+1}$,

${\left\| {{\varphi _k}} \right\|_{{L^\infty }}} \leq C\left\| {{\varphi _k}} \right\|_{{H^3}}^{\frac{1}{{q + 1}}}\left\| {{\varphi _k}} \right\|_{{L^q}}^{\frac{q}{{q + 1}}} \leq C\left\| {{\varphi _k}} \right\|_{{H^3}}^{\frac{1}{{q + 1}}}\left\| {{\varphi _k}} \right\|_{{H^1}}^{\frac{q}{{q + 1}}}.$

Proceeding as in (4.50), we find that for an arbitrary $\zeta \in L^{s}(0, T;H^{1})$, where $s \geq 1$ is yet to be determined,

$\begin{gathered} \left| {\int_0^T {\int_\Omega {{\varphi _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{dx dt}}} \right| \leq \int_0^T {} {\left\| {{v_k}} \right\|_{{L^2}}}{\left\| {{\varphi _k}} \right\|_{{L^\infty }}}{\left\| {\nabla {\Pi _k}\zeta } \right\|_{{L^2}}}{\text{dt}} \\ \leq C\left\| {{\varphi _k}} \right\|_{{L^\infty }({H^1})}^{\frac{q}{{q + 1}}}{\left\| {{v_k}} \right\|_{{L^2}({L^2})}}\left\| {{\varphi _k}} \right\|_{{L^2}({H^3})}^{\frac{1}{{q + 1}}}{\left\| \zeta \right\|_{{L^{\frac{{2(q + 1)}}{q}}}({H^1})}} \\ \leq \mathcal{E}{K^{\frac{1}{2}}}{\left\| \zeta \right\|_{{L^{\frac{{2(q + 1)}}{q}}}({H^1})}}, \\ \end{gathered}$

and so $\{{\text{div}}(\varphi_{k} {v}_{k})\}_{k \in \mathbb{N}}$ is uniformly bounded in the dual space of $L^{s}(0, T;H^{1})$ for $s = 2 + \frac{2}{q} \in (2, 4]$.

Remark 7.1. We point out that in the absence of the regularity result ${v}_{k} \in L^{q}(0, T;{H}^{1})$ from Lemma 7.1, and if we only have ${v}_{k} \in L^{2}(0, T;{L}^{2})$, then we obtain

$\left| {\int_0^T {\int_\Omega {{\sigma _k}} } {v_k} \cdot \nabla {\Pi _k}\zeta {\text{ dx dt}}} \right| \leq {\left\| {{\sigma _k}} \right\|_{{L^4}(Q)}}{\left\| {{v_k}} \right\|_{{L^2}({L^2})}}{\left\| {\nabla \zeta } \right\|_{{L^4}({L^4})}},$

and this implies that both $\{{\text{div}}(\sigma_{k} {v}_{k})\}_{k \in \mathbb{N}}$ and $\{\partial _{t}\sigma_{k}\}_{k \in \mathbb{N}}$ are bounded uniformly only in $L^{\frac{4}{3}}(0, T;(W^{1, 4})^{*})$.

---

8. Discussion

Reformulations of Darcy’s law and the pressure. Associated to Darcy's law (1.1b) is the term ${\lambda _v}: = p - \mu \varphi - \frac{D}{2}{\left| \sigma \right|^2}$ which will contribute the source term $\Gamma_{{v}} \lambda_{{v}}$ in the energy identity (4.4). In [20, Rmk. 2.1] three other reformulations of Darcy's law (1.1b) and the pressure are considered:

(R1) Let $q: = p - A\Psi (\varphi ) - \frac{B}{2}{\left| {\nabla \varphi } \right|^2}$ so that

${\lambda _v} = q + A\Psi (\varphi ) + \frac{B}{2}{\left| {\nabla \varphi } \right|^2} - \frac{D}{2}{\left| \sigma \right|^2} + \chi \sigma (1 - \varphi ) - \mu \varphi ,$

(8.1a)

$v = K(\nabla ( - q - \tfrac{B}{2}{\left| {\nabla \varphi } \right|^2}) - B\Delta \varphi \nabla \varphi ) = - K(\nabla q + B {\text{div}} (\nabla \varphi \otimes \nabla \varphi )).$

(8.1b)

(R2) Let $\hat p: = p + \frac{D}{2}{\left| \sigma \right|^2} + \chi \sigma (1 - \varphi )$ so that

${\lambda _v} = \hat p - \mu \varphi - D{\left| \sigma \right|^2} - \chi \sigma (1 - \varphi ),$

(8.2a)

$v = - K(\nabla \hat p - \mu \nabla \varphi - (D\sigma + \chi (1 - \varphi ))\nabla \sigma ).$

(8.2b)

(R3) Let $\tilde p: = p - \frac{D}{2}{\left| \sigma \right|^2} - \mu \varphi$ so that

${\lambda _v} = \tilde p,$

(8.3a)

$v = - K(\nabla \tilde p + \varphi \nabla \mu + \sigma \nabla (D\sigma + \chi (1 - \varphi ))).$

(8.3b)

From the viewpoint of estimating the source term $\Gamma_{{v}} \lambda_{{v}}$, we see that (8.3a) has the advantage of being the simplest. Meanwhile, for (8.2a) the analysis for $\Gamma_{{v}} \lambda_{{v}}$ is similar to that performed in Section 4.1.2, but for (8.1a) the main difficulty will be to estimate the terms $(A\Psi (\varphi ) + \frac{B}{2}{\left| {\nabla \varphi } \right|^2}){\Gamma _v}$ and $( - \frac{D}{2}{\left| \sigma \right|^2} + \chi \sigma \varphi ){\Gamma _v}$, which at first glance would require the assumption that $\Gamma_{{v}} \in L^{\infty}(Q)$, and obtaining an $L^{2}$-estimate for the pressure $q$ from the Darcy law (8.1b) would be difficult due to the term ${\text{div}} (\nabla \varphi \otimes \nabla \varphi)$.

Other boundary conditions for the pressure and velocity. In [20, x2.4.4] the authors have discussed possible boundary conditions for the velocity and for the pressure. As discussed in Section 2 following Assumption 2.1, we require the source term $\Gamma_{{v}}$ to have zero mean due to the no-flux boundary condition ${v} \cdot {n} = 0$ on $\partial \Omega$. The general energy identity (with homogeneous Neumann boundary conditions for $\varphi$ and $\mu$) from [20, Equ. (2.27)] reads as

$\begin{gathered} \frac{{\text{d}}}{{{\text{dt}}}}\int_\Omega A \Psi (\varphi ) + \frac{B}{2}{\left| {\nabla \varphi } \right|^2} + \frac{D}{2}{\left| \sigma \right|^2} + \chi \sigma (1 - \varphi ){\text{dx}} \hfill \\ \quad + \int_\Omega m (\varphi ){\left| {\nabla \mu } \right|^2} + n(\varphi ){\left| {\nabla (D\sigma + \chi (1 - \varphi ))} \right|^2} + \frac{1}{K}{\left| v \right|^2}{\text{dx}} \hfill \\ \quad = \int_\Omega {{\Gamma _\varphi }} \mu - \mathcal{S}(D\sigma + \chi (1 - \varphi )) + {\Gamma _v}{\lambda _v}{\text{dx}} \hfill \\ \quad + \int_{\partial \Omega } {(D\sigma + \chi (} 1 - \varphi ))n(\varphi )(D{\partial _n}\sigma ) - (v \cdot n)\left( {\frac{D}{2}{{\left| \sigma \right|}^2} + \chi \sigma (1 - \varphi ) + p} \right){\text{ d}}{H^{d - 1}}, \hfill \\ \end{gathered}$

and we see the appearance of an extra boundary source term involving the normal component of the velocity and the pressure. Here it would be advantageous to use the rescaled pressure $\hat{p}$ and the Darcy law (8.2b), as the extra boundary source term will become

$\int_{\partial \Omega } - (v \cdot n)\hat p{\text{ d}}{H^{d - 1}},$

which motivates the consideration of a Robin-type boundary condition for $\hat{p}$

$g = a\hat p - v \cdot n = a\hat p + K{\partial _n}\hat p - K(D\sigma + \chi (1 - \varphi )){\partial _n}\sigma {\text{ on }}\partial \Omega ,$

for some given datum $g$ and positive constant $a$. On one hand, this would allow us to consider source terms $\Gamma_{{v}}$ that need not have zero mean, but on the other hand, the analysis of the Darcy system becomes more complicated. In particular, the weak formulation of the pressure system now reads as

$\begin{gathered} \int_\Omega K \nabla \hat p \cdot \nabla \zeta {\text{dx}} + \int_{\partial \Omega } a \hat p\zeta {\text{d}}{H^{d - 1}} = \int_\Omega {{\Gamma _v}} \zeta + K\left( {\mu \nabla \varphi + \left( {D\sigma + \chi (1 - \varphi )} \right)\nabla \sigma } \right) \cdot \nabla \zeta {\text{dx}} \\ + \int_{\partial \Omega } g \zeta {\text{d}}{H^{d - 1}}, \\ \end{gathered}$

and we observe that the term $D\sigma \nabla \sigma$ on the right-hand side belongs to ${L}^{1}$ as $\sigma$ has at most $H^{1}$-spatial regularity from the energy identity. Thus, it is not clear if the pressure system can be solved with the regularities stated in Lemma 4.1. A deeper study into the theory of linear elliptic equations with right-hand sides of the form ${\text{div}} {f}$ where ${f} \in {L}^{1}$ is required.

---

Conflict of Interest

All authors declare no conflicts of interest in this paper.

---