Large-time behavior of the strong solution to nonhomogeneous incompressible MHD system with general initial data

Abstract

This paper investigates the large-time behavior of strong solutions to the nonhomogeneous incompressible magnetohydrodynamic equations on a bounded domain in $\mathbb{R}^{2}$. Based on uniform estimates, we prove that the velocity, the magnetic field, and their derivatives converge to zero in $L^{2}$ norm as time goes to infinity without any additional assumption on the initial data and external force by a pure energy method.

1 Introduction

Magnetohydrodynamics (MHD) studies the theory of the macroscopic interaction of electrically conducting fluids with a magnetic field. MHD has a very broad range of applications, such as, the intensely heated and ionized fluids in an electromagnetic field in astrophysics, geophysics, high-speed aerodynamics, and plasma physics. In this paper, we are concerned with the nonhomogeneous incompressible MHD equations in a domain $\Omega\subset\mathbb{R}^{2}$ as follows (see, e.g., [1–4]):

$$\begin{aligned}& \rho_{t}+\mathrm {u}\cdot\nabla\rho=0, \end{aligned}$$

(1.1)

$$\begin{aligned}& \rho \mathrm {u}_{t}+\rho \mathrm {u}\cdot\nabla \mathrm {u}+\nabla P =-\frac{1}{2} \nabla|\mathrm {H}|^{2}+\mathrm {H}\cdot \nabla \mathrm {H}+\mu\Delta \mathrm {u}+\rho \mathrm {f}, \end{aligned}$$

(1.2)

$$\begin{aligned}& \mathrm {H}_{t}-\nu\Delta \mathrm {H}+\mathrm {u}\cdot\nabla \mathrm {H}-\mathrm {H}\cdot\nabla \mathrm {u}=0, \end{aligned}$$

(1.3)

$$\begin{aligned}& \operatorname {div}\mathrm {H}=0, \qquad \operatorname {div}\mathrm {u}=0, \end{aligned}$$

(1.4)

where $t\geq0$ is time and $x\in\Omega$ is the spatial coordinate. The unknown functions are the density $\rho(x,t)\geq0$, velocity $\mathrm {u}=(\mathrm {u}^{1}(x,t),\mathrm {u}^{2}(x,t))\in \mathbb{R}^{2}$, the pressure $P(x,t)$, and the magnetic field $\mathrm {H}=(\mathrm {H}^{1}(x,t),\mathrm {H}^{2}(x,t))\in \mathbb{R}^{2}$, respectively. f stands for the external force. The positive constants μ and ν denote the viscosity of fluid and the relative strengths of advection and diffusion of H, respectively. Throughout this paper, Ω is assumed to be a bounded domain in $\mathbb{R}^{2}$ with smooth boundary. Without loss of generality, we will take $\mu=\nu=1$ for simplicity.

In this paper, we are interested in the initial-boundary value problem of (1.1)-(1.4) subject to the following initial conditions:

$$ (\rho, \mathrm {u},\mathrm {H}) (x,0)=(\rho_{0},\mathrm {u}_{0},\mathrm {H}_{0}) (x) \quad\mbox{and} \quad \underline{\rho}\leq\rho_{0}\leq\bar{\rho} \quad\mbox{for all }x\in\Omega, $$

(1.5)

and boundary conditions:

$$ \mathrm {u}|_{\partial\Omega}=(\mathrm {H}\cdot n)|_{\partial\Omega}=\nabla \times \mathrm {H}|_{\partial\Omega}=0, $$

(1.6)

where ρ̄ and $\underline{\rho}$ are two fixed positive constants. n is the unit outward normal to ∂Ω.

System (1.1)-(1.4) has drawn many attentions of engineers and applied mathematicians due to its important physical background and mathematical feature. The Cauchy problem to (1.1)-(1.4) has been much studied in the literature, including the existence, uniqueness, and regularity of solutions of the system. Zhang [5] established local classical solutions of (1.1)-(1.4) and proved that as the viscosity (μ) and resistivity (ν) went to zero, the solution of (1.1)-(1.4) converged to the solution of ideal MHD system (i.e. $\mu=\nu=0$). Gerbeau and Le Bris [3] and Desjardins and Le Bris [6] considered the global existence of weak solutions of finite energy in the whole space or in the torus. Global existence of strong solutions with small initial data in some Besov spaces was considered by Abidi and Paicu [7]. For the initial-boundary problem of (1.1)-(1.4), Huang and Wang [8] recently proved the unique global strong solution with initial vacuum in dimension two without external force f. For more related results, we refer the reader to [1, 9–12] and the references therein.

For homogeneous incompressible MHD ($\rho\equiv \mbox{const}$ in (1.1)-(1.4)), we first emphasize that the unique global classical solutions for every initial data $(\mathrm {u}_{0}, \mathrm {H}_{0})\in H^{m}$ with $m\geq2$ have been established in [13, 14]. Furthermore, some sufficient conditions for smoothness were presented for weak solutions to the MHD equations in [15, 16]. For Cauchy problem, the authors in [17] studied the long time behavior of solutions to the MHD equations in two and three dimensions with some smallness conditions. In [18], the authors considered the asymptotic behavior of the strong solutions to MHD equations in a half space. We also notice that if the partial derivatives of the viscosity and resistivity are zero, the global regularity issue has been established in [19].

However, the large-time asymptotic behavior of strong solutions to (1.1)-(1.6) has not been well understood, especially for the case with general initial data and external force. The main purpose of this paper is to investigate the influence of the magnetic fields, viscosity, and boundary effects on the behavior of the solution of (1.1)-(1.6). We show that the velocity, the magnetic field, and their derivatives will converge to zero as time tends to infinity. The main result of this paper is stated as follows.

Theorem 1.1

Let Ω be a bounded domain with smooth boundary ∂Ω in $\mathbb{R}^{2}$. Assume that the external force and the initial data satisfy

$$ \textstyle\begin{cases} \mathrm{f}\in L^{2}(0,\infty; H^{1}(\Omega)), \quad\mathrm{f}_{t}\in L^{2}(0,\infty; L^{2}(\Omega)),\\ (\rho_{0}, \mathrm {u}_{0}, \mathrm {H}_{0})\in H^{2}, \quad\underline{\rho}\leq\rho_{0}\leq \bar{\rho} \textit{ for all } x\in\Omega, \end{cases} $$

(1.7)

then there exists a unique global solution $(\rho, \mathrm {u}, \mathrm {H})$ to (1.1)-(1.6) satisfying

$$ \textstyle\begin{cases} \rho\in C([0,\infty);H^{2}(\Omega)) \quad\textit{with } \underline{\rho }\leq\rho\leq\bar{\rho},\\ (\mathrm {u},\mathrm {H})\in C([0,\infty);H^{2}(\Omega))\cap L^{2}(0,\infty;H^{3}(\Omega)). \end{cases} $$

(1.8)

In particular, the following large-time behavior holds:

$$\begin{aligned}& \lim_{t\to\infty}\|\mathrm {u}\|_{H^{1}}=0,\qquad \lim _{t\to\infty}\|\mathrm {H}\| _{H^{1}}=0, \end{aligned}$$

(1.9)

$$\begin{aligned}& \lim_{t\to\infty} \bigl(\|\mathrm {u}_{t}\|_{L^{2}}+\| \mathrm {H}_{t}\|_{L^{2}} \bigr)=0. \end{aligned}$$

(1.10)

Remark 1.1

If $B=0$ (no magnetic field), then the system (1.1)-(1.4) becomes the classical Navier-Stokes system. The global existence and large-time behavior of the strong solution to the nonhomogeneous incompressible Navier-Stokes equations has been proved by Zhao [20].

Remark 1.2

We should point out that Huang and Wang in [8] have established a global existence of the strong solution to (1.1)-(1.4) with initial vacuum and $\mathrm {f}=0$. Theorem 1.1 is still an interesting result, because we get the large-time behavior.

Our second main result gives an exponential decay rate for the solution on taking the external force $f=0$, which is stated as follows.

Theorem 1.2

Assume the conditions in Theorem  1.1 and take the external force $f=0$ in (1.2). Then there exists a unique global solution $(\rho, \mathrm {u}, \mathrm {H})$ to (1.1)-(1.6) satisfying (1.8). In particular, one has the following exponential decay:

$$ \|\mathrm {u}\|^{2}_{L^{2}}+\|\mathrm {H}\|^{2}_{L^{2}} \leq\exp\{-C_{1}t\}, $$

(1.11)

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

Remark 1.3

Theorem 1.2 establishes exponential decay. Indeed, if we imposed some conditions on the external force f but would not take $f=0$, we can also give an algebraic rate by performing energy estimates with weights in time.

We will prove Theorem 1.1 via the pure energy method, which is based on the uniform estimate for the local solution. This approach is motivated by the previous work on the nonhomogeneous incompressible Navier-Stokes equations due to Zhao [20]. Due to the strong coupling between the velocity, magnetic field and density equations, the nonhomogeneous problem (1.1)-(1.4) under consideration is much more involved. It should be pointed out that the analysis for (1.1)-(1.4) is very complicated and more efforts should be made as regards the estimates involving these coupling terms and space-time-dependent external force.

Next, we give a brief outline of the proof. To obtain the large-time behavior, we must establish the uniform estimate. As for the uniform estimate, it is difficult to estimate the higher-order derivative of u and H due to the lack of the spatial derivatives of the solution at the boundary. So it is divided into two steps: first of all applying the standard energy estimate on the temporal derivatives of the solution, and then spatial derivatives by applying regularity theory of Stokes equation repeatedly. With all the estimates in hand, we can establish the global existence by a standard method. Finally, the proof of the large-time behavior is mainly based on the following fact: if $f(t)\in W^{1,1}(0,\infty)$, then $f(t)\rightarrow0$ as $t\rightarrow\infty$. By combining these uniform estimates in Lemmas 3.1-3.3 and the fact above we prove the desired large-time behavior of the solution in Section 4.

The rest of the paper is organized as follows. In Section 2, we state some well-known inequalities, and in Section 3 we deduce some uniform estimates and the time-dependent estimate of the higher-order derivative. Finally, we will complete the proof of Theorem 1.1 in Section 4.

2 Auxiliary lemmas

In this section, we state some inequalities which will be frequently used in this paper. First of all, we recall some Sobolev embedding inequalities (see [21, 22]).

Lemma 2.1

The following inequalities hold:

$$\begin{aligned}& (1) \quad\|f\|_{L^{p}}\leq C\|f\|_{H^{1}}, \quad\forall 1\leq p< \infty, \forall f\in H^{1}, \end{aligned}$$

(2.1)

$$\begin{aligned}& (2)\quad \|f\|_{L^{\infty}}\leq C\|f\|_{W^{1,p}},\quad \forall 2< p< \infty , \forall f\in W^{1,p}, \end{aligned}$$

(2.2)

$$\begin{aligned}& (3)\quad \|f\|^{2}_{L^{4}}\leq C\|f\|_{L^{2}}\|\nabla f \|_{L^{2}}, \quad\forall f\in H_{0}^{1}, \end{aligned}$$

(2.3)

$$\begin{aligned}& (4)\quad \|f\|_{L^{p}}\leq C\|\nabla f\|_{L^{2}}, \quad\forall 1 \leq p< \infty , \forall f\in H_{0}^{1}, \end{aligned}$$

(2.4)

$$\begin{aligned}& (5)\quad \|f\|^{2}_{L^{4}}\leq C\bigl(\|f\|^{2}_{L^{2}}+ \|f\|_{L^{2}}\|\nabla f\| _{L^{2}}\bigr), \quad\forall f\in H^{1}. \end{aligned}$$

(2.5)

Next, we state the following Sobolev embedding (see [23, 24]), which will be used to deal with the estimate of H.

Lemma 2.2

Let $v \in H^{1}(\Omega)$ be a vector function satisfying $v|_{\partial\Omega}=0$ or $v\cdot n|_{\partial \Omega}=0$, then

$$\begin{aligned} \|v\|_{L^{2}}\leq C\|\nabla v\|_{L^{2}}. \end{aligned}$$

(2.6)

3 Some a priori estimates

The local existence has been established by Chen et al. [9]. Hence, it only remains to establish some necessary a priori bounds for the strong solutions $(\rho, \mathrm {u},\mathrm {H})$ to the initial-boundary value problem (1.1)-(1.4).

Let $T>0$ be a fixed time and $(\rho, \mathrm {u}, \mathrm {H})$ be the strong solution to (1.1)-(1.4) defined on $\Omega\times(0,T]$. Throughout this paper, we will denote by C the various generic positive constants, which may depend on the initial data and ρ̄, $\underline{\rho}$, f, and Ω but are independent of t. A special dependence will be pointed out explicitly in this paper if necessary.

3.1 Uniform estimate

First, by the method of characteristics and standard energy we have the following uniform estimates:

$$ \underline{\rho}\leq\rho(x,t)\leq\bar{\rho} $$

(3.1)

and

$$\begin{aligned} \sup_{0\leq t\leq \infty} \int \biggl(| \mathrm {u}|^{2}+\frac{1}{2}| \mathrm {H}|^{2} \biggr) \,\mathrm{d} x+ \int_{0}^{\infty}\bigl(\| \nabla \mathrm {u}\|_{L^{2}}^{2}+\|\nabla \mathrm {H}\| _{L^{2}}^{2} \bigr) \,\mathrm{d} t\leq C. \end{aligned}$$

(3.2)

The next lemma is the crucial estimate in this paper. Higher-order estimates of the density, velocity and magnetic field can be obtained in a standard way provided that $\|\mathrm {u}\|_{H^{1}}$ and $\|\mathrm {H}\|_{H^{1}}$ are uniformly bounded with respect to time.

Lemma 3.1

Let $(\rho, \mathrm {u}, \mathrm {H})$ be a smooth solution of (1.1)-(1.6) on $\Omega\times(0,\infty)$. Then there exists a constant C such that

$$\begin{aligned} \sup_{0\leq t\leq \infty} \bigl(\|\mathrm {u}\|_{H^{1}}^{2}+\|\mathrm {H}\|_{H^{1}}^{2} \bigr)+ \int_{0}^{\infty}\bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}^{2}+\|\mathrm {H}_{t} \|_{L^{2}}^{2}+\|\mathrm {H}\| _{H^{2}}^{2}+\|\mathrm {u}\|_{H^{2}}^{2} \bigr)\,\mathrm{d} t\leq C. \end{aligned}$$

(3.3)

Proof

Multiplying (1.2) by $\mathrm {u}_{t}$ and integrating by parts over Ω, one obtains

$$\begin{aligned} \frac{1}{2}\frac{\mathrm {d}}{\mathrm {d}t}\|\nabla \mathrm {u}\|_{L^{2}}^{2}+\bigl\| \rho^{1/2}\mathrm {u}_{t}\bigr\| _{L^{2}}^{2} = \int (\mathrm {H}\cdot\nabla \mathrm {H}\cdot \mathrm {u}_{t}-\rho \mathrm {u}\cdot\nabla \mathrm {u}\cdot \mathrm {u}_{t}+\rho \mathrm {f}\cdot u_{t} ) \,\mathrm {d}x. \end{aligned}$$

(3.4)

Similarly, it follows from (1.3) that

$$\begin{aligned} \frac{\mathrm {d}}{\mathrm {d}t}\|\nabla \mathrm {H}\|_{L^{2}}^{2}+ \bigl(\| \mathrm {H}_{t}\|_{L^{2}}^{2}+\|\triangle \mathrm {H}\|_{L^{2}} \bigr)= \int|\mathrm {H}_{t}-\Delta \mathrm {H}|^{2}\,\mathrm {d}x= \int|\mathrm {H}\cdot\nabla \mathrm {u}-\mathrm {u}\cdot\nabla \mathrm {H}|^{2}\,\mathrm {d}x. \end{aligned}$$

(3.5)

Putting (3.4) and (3.5) together leads to

$$\begin{aligned} &\frac{\mathrm {d}}{\mathrm {d}t} \biggl(\frac{1}{2}\|\nabla \mathrm {u}\|_{L^{2}}^{2}+\|\nabla \mathrm {H}\|_{L^{2}}^{2} \biggr)+ \bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}^{2}+ \|\mathrm {H}_{t}\|_{L^{2}}^{2}+\|\triangle \mathrm {H}\|_{L^{2}}^{2} \bigr) \\ & \quad\leq \int|\mathrm {H}||\nabla \mathrm {H}||\mathrm {u}_{t}|\,\mathrm {d}x+ \int\rho|\mathrm {u}||\nabla \mathrm {u}||\mathrm {u}_{t}| \,\mathrm {d}x+2 \int|\mathrm {H}|^{2}|\nabla \mathrm {u}|^{2}\,\mathrm {d}x \\ &\qquad{} +2 \int|\mathrm {u}|^{2}|\nabla \mathrm {H}|^{2}\,\mathrm {d}x+ \int\rho| \mathrm {f}||\mathrm {u}_{t}| \,\mathrm {d}x=\sum _{i=1}^{5}I_{i}. \end{aligned}$$

(3.6)

Next, we will estimate all the terms on the right-hand side of (3.6) term by term. Using the Young inequality and (3.1), we have

$$\begin{aligned} I_{1} \leq&\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| ^{2}_{L^{2}}+C_{1}\|\mathrm {H}\| ^{2}_{L^{4}} \|\nabla \mathrm {H}\|^{2}_{L^{4}} \\ \leq&\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t}\bigr\| ^{2}_{L^{2}}+ \varepsilon\|\Delta \mathrm {H}\|_{L^{2}}+C(\varepsilon)\|\nabla \mathrm {H}\|^{4}_{L^{2}}, \end{aligned}$$

where we use the fact

$$\begin{aligned} \|\mathrm {H}\|^{2}_{L^{4}}\|\nabla \mathrm {H}\|^{2}_{L^{4}} \leq& C\bigl(\|\mathrm {H}\|^{2}_{L^{2}}+ \|\mathrm {H}\|_{L^{2}}\|\nabla \mathrm {H}\|_{L^{2}}\bigr) \bigl(\|\nabla \mathrm {H}\|_{L^{2}}\|\Delta \mathrm {H}\|_{L^{2}}+\|\nabla \mathrm {H}\|^{2}_{L^{2}} \bigr) \\ \leq&\|\mathrm {H}\| ^{2}_{L^{2}}\|\nabla \mathrm {H}\|_{L^{2}}\|\Delta \mathrm {H}\|_{L^{2}}+\|\mathrm {H}\|^{2}_{L^{2}}\| \nabla \mathrm {H}\|^{2}_{L^{2}} \\ &{}+\|\mathrm {H}\|_{L^{2}}\|\nabla \mathrm {H}\| ^{2}_{L^{2}}\|\Delta \mathrm {H}\|_{L^{2}}+\|\mathrm {H}\|^{2}_{L^{2}}\|\nabla \mathrm {H}\|^{3}_{L^{2}} \\ \leq&\frac{\varepsilon}{C_{1}}\|\Delta \mathrm {H}\|_{L^{2}}+C\| \nabla \mathrm {H}\|^{4}_{L^{2}}, \end{aligned}$$

(3.7)

due to (3.2), (2.5), and (2.6). Similarly, we have

$$\begin{aligned} I_{2} \leq&\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| ^{2}_{L^{2}}+C(\varepsilon)\| \mathrm {u}\|^{2}_{L^{4}} \|\nabla \mathrm {u}\|^{2}_{L^{4}} \\ \leq&\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t}\bigr\| ^{2}_{L^{2}}+ \varepsilon\bigl\| \nabla ^{2} u\bigr\| _{L^{2}}+C\|\nabla \mathrm {u}\|^{4}_{L^{2}}. \end{aligned}$$

Next, we turn to estimating $\|\nabla^{2} u\|_{L^{2}}$. From (1.2), we know that u satisfies the following Stokes equations:

$$ \left \{ \textstyle\begin{array}{@{}l} -\triangle \mathrm {u}+\nabla P=-\rho \mathrm {u}_{t}-\rho \mathrm {u}\cdot\nabla \mathrm {u}+\frac{1}{2}\nabla|\mathrm {H}|^{2}-\mathrm {H}\cdot\nabla \mathrm {H}+\rho \mathrm {f},\\ \operatorname {div}\mathrm {u}=0, \quad\mbox{in } \Omega,\\ \mathrm {u}=0, \quad\mbox{on } \partial\Omega. \end{array}\displaystyle \right . $$

By the well-known regularity theory on the Stokes equations (see [23]) and using (2.5), (3.2), and (3.7), we have

$$\begin{aligned} \|\mathrm {u}\|_{H^{2}} \leq& C\bigl(\|\rho \mathrm {u}_{t}\|_{L^{2}}+\| \rho \mathrm {u}\cdot\nabla \mathrm {u}\|_{L^{2}}+\|\nabla|\mathrm {H}|^{2} \|_{L^{2}}+\|\mathrm {H}\cdot\nabla \mathrm {H}\| _{L^{2}}+\|\mathrm {f}\|_{L^{2}} \bigr) \\ \leq& C\bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}+\|\mathrm {u}\|_{L^{4}}\|\nabla \mathrm {u}\| _{L^{4}}+\|\mathrm {H}\|_{L^{4}}\|\nabla \mathrm {H}\|_{L^{4}}+\|\mathrm {f}\|_{L^{2}}\bigr) \\ \leq& C\bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}+\|\nabla \mathrm {u}\|_{L^{2}}\bigl\| \nabla^{2} \mathrm {u}\bigr\| ^{1/2}_{L^{2}} \\ &{}+\| \nabla \mathrm {u}\|^{2}_{L^{2}}+\|\nabla \mathrm {H}\|_{L^{2}}\| \Delta \mathrm {H}\|^{1/2}_{L^{2}}+\|\nabla \mathrm {H}\|^{2}_{L^{2}}+\|\mathrm {f}\| _{L^{2}}\bigr) \\ \leq& C\bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}+\|\nabla \mathrm {u}\|^{2}_{L^{2}}+\|\nabla \mathrm {H}\|^{2}_{L^{2}}+\| \mathrm {f}\|_{L^{2}}\bigr)+\frac{1}{2}\bigl\| \nabla^{2} \mathrm {u}\bigr\| ^{1/2}_{L^{2}}+\|\Delta \mathrm {H}\|^{1/2}_{L^{2}}, \end{aligned}$$

which immediately leads to

$$\begin{aligned} \|\mathrm {u}\|_{H^{2}}\leq C\bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t} \bigr\| _{L^{2}}+\|\nabla \mathrm {u}\| ^{2}_{L^{2}}+\|\nabla \mathrm {H}\|^{2}_{L^{2}}+\|\mathrm {f}\|_{L^{2}}\bigr)+2\|\Delta \mathrm {H}\|_{L^{2}}. \end{aligned}$$

Substituting the above inequality into $I_{2}$, we have

$$\begin{aligned} I_{2}\leq C\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| ^{2}_{L^{2}}+\varepsilon\|\Delta \mathrm {H}\|^{2}_{L^{2}}+C \|\nabla \mathrm {u}\|^{4}_{L^{2}}+C\|\nabla \mathrm {u}\|^{2}_{L^{2}} \| \nabla \mathrm {H}\|^{2}_{L^{2}}+C\|\mathrm {f}\|^{2}_{L^{2}}. \end{aligned}$$

Similarly, for $I_{3}$ and $I_{4}$, we have

$$\begin{aligned} I_{3}\leq\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| ^{2}_{L^{2}}+\varepsilon\|\Delta \mathrm {H}\|^{2}_{L^{2}}+C \|\nabla \mathrm {H}\|^{2}_{L^{2}}\|\nabla \mathrm {u}\|^{2}_{L^{2}}+C \| \nabla \mathrm {u}\|^{4}_{L^{2}}+C\|\nabla \mathrm {H}\|^{4}_{L^{2}} \end{aligned}$$

and

$$\begin{aligned} I_{4} \leq& C\|\mathrm {u}\|^{2}_{L^{4}}\|\nabla \mathrm {H}\|^{2}_{L^{4}}\leq C\|\mathrm {u}\|_{L^{2}}\| \nabla \mathrm {u}\|_{L^{2}}\bigl(\|\nabla \mathrm {H}\|^{2}_{L^{2}}+\|\nabla \mathrm {H}\|_{L^{2}}\| \Delta \mathrm {H}\|_{L^{2}}\bigr) \\ \leq& C\|\nabla \mathrm {H}\|^{2}_{L^{2}}\|\nabla \mathrm {u}\|^{2}_{L^{2}}+\varepsilon\| \Delta \mathrm {H}\|^{2}_{L^{2}}, \end{aligned}$$

due to (3.2).

Using the Young inequality immediately leads to

$$\begin{aligned} I_{5}\leq\varepsilon\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| ^{2}_{L^{2}}+C(\varepsilon)\| \mathrm {f}\|^{2}_{L^{2}}. \end{aligned}$$

Taking ε small enough and substituting $I_{1}$-$I_{5}$ into (3.6), one obtains

$$\begin{aligned} &\frac{\mathrm {d}}{\mathrm {d}t} \bigl(\|\nabla \mathrm {u}\|_{L^{2}}^{2}+\|\nabla \mathrm {H}\|_{L^{2}}^{2} \bigr)+ \bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}^{2}+\|\mathrm {H}_{t} \|_{L^{2}}^{2}+\|\triangle \mathrm {H}\|_{L^{2}}^{2} \bigr) \\ &\quad\leq C\bigl(\|\nabla \mathrm {H}\|^{2}_{L^{2}}+\|\nabla \mathrm {u}\|^{2}_{L^{2}}\bigr)^{2}+\|\mathrm {f}\|^{2}_{L^{2}}, \end{aligned}$$

which, together with Gronwall’s inequality, immediately leads to

$$\begin{aligned} \sup_{0\leq t\leq \infty} \bigl(\|\nabla \mathrm {u}\|_{L^{2}}^{2}+ \|\nabla \mathrm {H}\|_{L^{2}}^{2} \bigr)+ \int_{0}^{\infty}\bigl(\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}^{2}+\|\mathrm {H}_{t} \|_{L^{2}}^{2}+\bigl\| \nabla^{2} \mathrm {H}\bigr\| _{L^{2}}^{2} \bigr)\,\mathrm{d} t\leq C. \end{aligned}$$

(3.8)

We get the desired estimate (3.3) by (3.2) and (3.8). □

The following lemma is concerned with the $L^{2}$-estimate of $\rho^{1/2} \mathrm {u}_{t}$ and $\mathrm {H}_{t}$.

Lemma 3.2

Let $(\rho, \mathrm {u}, \mathrm {H})$ be a strong solution of (1.1)-(1.6) on $\Omega\times(0,\infty)$. Then there exists a constant C such that

$$\begin{aligned} \sup_{0\leq t\leq \infty} \bigl(\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| _{L^{2}}^{2}+\|\mathrm {H}_{t}\|_{L^{2}}^{2} \bigr)+ \int_{0}^{\infty}\|\nabla \mathrm {u}_{t} \|_{L^{2}}^{2}+\|\nabla \mathrm {H}_{t}\| _{L^{2}}^{2} \,\mathrm{d} t \leq C. \end{aligned}$$

(3.9)

Proof

Differentiating the momentum equations (1.2) with respect to t yields

$$\begin{aligned} &\rho \mathrm {u}_{tt}+\rho \mathrm {u}\cdot\nabla \mathrm {u}_{t}+\rho \mathrm {u}_{t}\cdot\nabla \mathrm {u}+\rho_{t} (\mathrm {u}_{t}+\mathrm {u}\cdot \nabla \mathrm {u})+\nabla P_{t}\\ &\quad=\triangle \mathrm {u}_{t}+\biggl(\mathrm {H}\cdot \nabla \mathrm {H}-\frac{1}{2}\nabla|\mathrm {H}|^{2}\biggr)_{t}+(\rho \mathrm {f})_{t}. \end{aligned}$$

Multiplying the equation above with $\mathrm {u}_{t}$ and integrating by parts over Ω, one gets

$$\begin{aligned} &\frac{1}{2}\frac{\mathrm {d}}{\mathrm {d}t}\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| _{L^{2}}^{2}+\|\nabla \mathrm {u}_{t} \|_{L^{2}}^{2} \\ &\quad=- \int\rho_{t}|\mathrm {u}_{t}|^{2}\,\mathrm {d}x- \int\rho(\mathrm {u}_{t}\cdot\nabla \mathrm {u})\cdot \mathrm {u}_{t}\,\mathrm{d}x- \int\rho_{t}(\mathrm {u}\cdot\nabla \mathrm {u})\cdot \mathrm {u}_{t}\,\mathrm {d}x \\ &\qquad{} + \int(\mathrm {H}_{t}\cdot\nabla \mathrm {H}+\mathrm {H}\cdot\nabla \mathrm {H}_{t}+ \rho _{t}\mathrm {f}+\rho \mathrm {f}_{t})\cdot \mathrm {u}_{t}\,\mathrm {d}x. \end{aligned}$$

(3.10)

Differentiating (1.3) with respect to t and multiplying the resulting equation by $\mathrm {H}_{t}$, we obtain after integrating by parts

$$\begin{aligned} &\frac{1}{2}\frac{\mathrm {d}}{\mathrm {d}t}\|\mathrm {H}_{t}\|_{L^{2}}^{2}+ \|\nabla \mathrm {H}_{t}\| _{L^{2}}^{2} \\ &\quad= - \int \mathrm {u}_{t}\cdot\nabla \mathrm {H}\cdot \mathrm {H}_{t} \,\mathrm {d}x + \int \mathrm {H}_{t}\cdot\nabla \mathrm {u}\cdot \mathrm {H}_{t}+\mathrm {H}\cdot \nabla \mathrm {u}_{t}\cdot \mathrm {H}_{t}\,\mathrm {d}x. \end{aligned}$$

(3.11)

Putting (3.10) and (3.11) together leads to

$$\begin{aligned} &\frac{1}{2}\frac{\mathrm {d}}{\mathrm {d}t} \bigl(\|\mathrm {H}_{t} \|_{L^{2}}^{2}+\bigl\| \rho ^{1/2}\mathrm {u}_{t} \bigr\| _{L^{2}}^{2} \bigr)+\|\nabla \mathrm {u}_{t} \|_{L^{2}}^{2}+\|\nabla \mathrm {H}_{t}\|_{L^{2}}^{2} \\ &\quad = - \int\rho_{t}|\mathrm {u}_{t}|^{2}\,\mathrm {d}x- \int\rho(\mathrm {u}_{t}\cdot\nabla \mathrm {u})\cdot \mathrm {u}_{t}\,\mathrm {d}x- \int\rho_{t}(\mathrm {u}\cdot\nabla \mathrm {u})\cdot \mathrm {u}_{t}\,\mathrm {d}x \\ &\qquad{}+ \int\rho_{t}\mathrm {f}\cdot \mathrm {u}_{t}\,\mathrm {d}x+ \int\rho \mathrm {f}_{t}\cdot \mathrm {u}_{t}\,\mathrm {d}x \\ &\qquad{} + \int \mathrm {H}_{t}\cdot\nabla \mathrm {H}\cdot \mathrm {u}_{t}- \mathrm {u}_{t}\cdot\nabla \mathrm {H}\cdot \mathrm {H}_{t} \,\mathrm {d}x + \int \mathrm {H}_{t}\cdot\nabla \mathrm {u}\cdot \mathrm {H}_{t}\,\mathrm {d}x \\ &\qquad{}+ \int \mathrm {H}\cdot\nabla \mathrm {u}_{t}\cdot \mathrm {H}_{t}+\mathrm {H}\cdot \nabla \mathrm {H}_{t}\cdot \mathrm {u}_{t}\,\mathrm {d}x \\ &\quad\triangleq\sum_{i=1}^{8}R_{i}. \end{aligned}$$

(3.12)

We now estimate each term on the right-hand side of (3.12). First, using the Hölder and Young inequalities, we obtain

$$\begin{aligned} |R_{1}| =&\biggl\vert \int\rho \mathrm {u}\cdot\nabla|\mathrm {u}_{t}|^{2}\,\mathrm {d}x\biggr\vert \leq C\|\nabla \mathrm {u}_{t}\|_{L^{2}}\|\mathrm {u}\|_{L^{4}}\| \mathrm {u}_{t}\|_{L^{4}} \\ \leq& C\|\nabla \mathrm {u}_{t}\|_{L^{2}}\| \mathrm {u}\|^{1/2}_{L^{2}} \|\nabla \mathrm {u}\| ^{1/2}_{L^{2}}\| \mathrm {u}_{t} \|^{1/2}_{L^{2}}\|\nabla \mathrm {u}_{t}\| ^{1/2}_{L^{2}} \\ \leq& C\| \mathrm {u}_{t}\|^{1/2}_{L^{2}}\|\nabla \mathrm {u}_{t}\|^{3/2}_{L^{2}}\leq \varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C(\varepsilon)\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| ^{2}_{L^{2}}, \\ |R_{2}| \leq& C \int|\mathrm {u}_{t}|^{2}|\nabla \mathrm {u}|\,\mathrm {d}x\leq C\|\nabla \mathrm {u}\| _{L^{2}}\|\mathrm {u}_{t}\|^{2}_{L^{4}} \\ \leq& C\|\mathrm {u}_{t}\|_{L^{2}}\|\nabla \mathrm {u}_{t} \|_{L^{2}}\leq\varepsilon\| \nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\bigl\| \rho^{1/2}\mathrm {u}_{t}\bigr\| ^{2}_{L^{2}}, \end{aligned}$$

where (3.1), (3.2), and (3.3) are all used.

The estimate of $R_{3}$ is given as follows:

$$\begin{aligned} |R_{3}| =& \biggl\vert \int\rho \mathrm {u}\cdot\nabla(\mathrm {u}\cdot\nabla \mathrm {u}\cdot \mathrm {u}_{t}) \,\mathrm {d}x \biggr\vert \\ \leq& C \int|\mathrm {u}||\nabla \mathrm {u}|^{2}|\mathrm {u}_{t}|+| \mathrm {u}|^{2}|\nabla^{2} \mathrm {u}||\mathrm {u}_{t}|+| \mathrm {u}|^{2}|\nabla \mathrm {u}||\nabla \mathrm {u}_{t}|\,\mathrm {d}x\triangleq\sum _{i=1}^{3}J_{i}. \end{aligned}$$

Using (2.4) and (3.3), we can make the reduction

$$\begin{aligned} |J_{1}| \leq& C \int|\mathrm {u}||\nabla \mathrm {u}|^{2}|\mathrm {u}_{t}| \,\mathrm {d}x\leq C\| \nabla \mathrm {u}\|^{2}_{L^{4}}\|\mathrm {u}\|_{L^{4}}\| \mathrm {u}_{t}\|_{L^{4}} \\ \leq& C\bigl(\|\nabla \mathrm {u}\|^{2}_{L^{2}}+\|\nabla \mathrm {u}\|_{L^{2}}\bigl\| \nabla^{2} \mathrm {u}\bigr\| _{L^{2}}\bigr)\|\nabla \mathrm {u}\|_{L^{2}}\|\nabla \mathrm {u}_{t}\|_{L^{2}} \\ \leq& C\bigl(1+\bigl\| \nabla^{2} \mathrm {u}\bigr\| _{L^{2}}\bigr)\|\nabla \mathrm {u}\|_{L^{2}}\|\nabla \mathrm {u}_{t}\|_{L^{2}} \\ \leq&\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\|\nabla \mathrm {u}\|^{2}_{L^{2}}+ C(\varepsilon)\bigl\| \nabla^{2} \mathrm {u}\bigr\| ^{2}_{L^{2}}, \\ |J_{2}| \leq& C\bigl\| \nabla^{2} \mathrm {u}\bigr\| _{L^{2}}\|\mathrm {u}\|^{2}_{L^{8}}\|\mathrm {u}_{t}\| _{L^{4}} \\ \leq& C\bigl\| \nabla^{2} \mathrm {u}\bigr\| _{L^{2}}\|\nabla \mathrm {u}\|^{2}_{L^{2}}\|\nabla \mathrm {u}_{t}\| _{L^{2}} \\ \leq&\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\bigl\| \nabla ^{2} \mathrm {u}\bigr\| ^{2}_{L^{2}}, \\ |J_{3}| \leq& C\|\nabla \mathrm {u}\|_{L^{4}}\|\mathrm {u}\|^{2}_{L^{8}} \|\nabla \mathrm {u}_{t}\| _{L^{2}}\leq\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C(\varepsilon)\bigl\| \nabla^{2} \mathrm {u}\bigr\| ^{2}_{L^{2}}+C(\varepsilon)\|\nabla \mathrm {u}\|^{2}_{L^{2}}. \end{aligned}$$

From the estimate of $J_{1}$-$J_{3}$, we obtain

$$\begin{aligned} |R_{3}|\leq3\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C \bigl\| \nabla^{2} \mathrm {u}\bigr\| ^{2}_{L^{2}}+C\|\nabla \mathrm {u}\|^{2}_{L^{2}}. \end{aligned}$$

For $R_{4}$-$R_{7}$, we have

$$\begin{aligned} |R_{4}| =&\biggl\vert \int\rho \mathrm {u}\cdot\nabla(\mathrm {f}\cdot \mathrm {u}_{t})\,\mathrm {d}x\biggr\vert \\ \leq& C\|\mathrm {u}\|_{L^{4}}\|\mathrm {u}_{t}\|_{L^{4}}\|\nabla \mathrm {f}\|_{L^{2}}+C\|\mathrm {u}\| _{L^{4}}\|\nabla \mathrm {u}_{t} \|_{L^{2}}\|\mathrm {f}\|_{L^{4}} \\ \leq& C\|\nabla \mathrm {u}_{t}\|_{L^{2}}\|\nabla \mathrm {f}\|_{L^{2}}+C\|\mathrm {f}\|_{H^{1}}\| \nabla \mathrm {u}_{t} \|_{L^{2}} \\ \leq&\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\|\mathrm {f}\| ^{2}_{H^{1}}, \\ |R_{5}| \leq&\bigl\| \rho^{1/2}\mathrm {u}_{t} \bigr\| ^{2}_{L^{2}}+\|\mathrm {f}_{t}\|^{2}_{L^{2}}, \\ |R_{6}| \leq& C \int|\mathrm {H}_{t}||\nabla \mathrm {u}||\mathrm {u}_{t}|\,\mathrm {d}x \\ \leq& C\|\nabla \mathrm {u}\|_{L^{4}}\|\mathrm {H}_{t}\|_{L^{2}}\| \mathrm {u}_{t}\|_{L^{4}} \\ \leq& C\|\nabla \mathrm {u}\|^{1/2}_{H^{1}}\|\mathrm {H}_{t} \|_{L^{2}}\|\nabla \mathrm {u}_{t}\| _{L^{2}} \\ \leq&\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\|\nabla \mathrm {u}\|_{H^{1}}\|\mathrm {H}_{t}\|^{2}_{L^{2}} \\ \leq&\varepsilon\|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\|\mathrm {H}_{t}\| ^{4}_{L^{2}}+C(\varepsilon)\| \mathrm {u}\|^{2}_{H^{2}}, \\ |R_{7}| \leq& C \int|\mathrm {H}_{t}|^{2}|\nabla \mathrm {u}|\,\mathrm {d}x\leq C\| \mathrm {H}_{t}\|^{2}_{L^{4}}\| \nabla \mathrm {u}\|_{L^{2}} \\ \leq& C\bigl(\|\mathrm {H}_{t}\|^{2}_{L^{2}}+\| \mathrm {H}_{t}\|_{L^{2}}\|\nabla \mathrm {H}_{t}\| _{L^{2}} \bigr) \\ \leq&\varepsilon\|\nabla \mathrm {H}_{t}\|^{2}_{L^{2}}+C( \varepsilon)\| \mathrm {H}_{t}\|^{2}_{L^{2}}. \end{aligned}$$

It is easy to prove that $R_{8}=0$. Taking ε small enough and substituting $R_{1}$-$R_{8}$ into (3.12), one obtains

$$\begin{aligned} &\frac{\mathrm {d}}{\mathrm {d}t} \bigl(\|\mathrm {H}_{t}\|_{L^{2}}^{2}+\bigl\| \rho^{1/2}\mathrm {u}_{t}\bigr\| _{L^{2}}^{2} \bigr)+ \bigl(\|\nabla \mathrm {u}_{t}\|_{L^{2}}^{2}+ \|\nabla \mathrm {H}_{t}\|_{L^{2}}^{2} \bigr) \\ &\quad \leq C\|\nabla \mathrm {u}\|^{2}_{L^{2}}+ C\|\mathrm {u}\|^{2}_{H^{2}}+C\| \mathrm {H}_{t}\| ^{2}_{L^{2}}+C \| \mathrm {H}_{t}\|^{4}_{L^{2}}+C\bigl\| \rho^{1/2} \mathrm {u}_{t}\bigr\| ^{2}_{L^{2}}+C\| \nabla \mathrm {f}\|^{2}_{L^{2}}+C\|\mathrm {f}_{t}\|^{2}_{L^{2}}, \end{aligned}$$

which, together with Gronwall’s inequality and (3.3), immediately leads to the desired estimate (3.9). □

3.2 Time-dependent estimate for higher derivative

It suffices to prove the large-time behavior with the help of the uniform estimates in Lemmas 3.1 and 3.2. Lemma 3.3 below deals with the higher-order estimates of the solutions which are needed to guarantee the extension of a local classical solution to a global one.

Lemma 3.3

Let $(\rho, \mathrm {u}, \mathrm {H})$ be a strong solution of (1.1)-(1.6) on $\Omega\times(0,T)$. Then there exists a constant $C(T)$ such that

$$\begin{aligned} \sup_{0\leq t\leq T} \bigl(\|\rho\|_{H^{2}}+\|\mathrm {u}\|_{H^{2}}+\| \mathrm {H}\|_{H^{2}} \bigr)+ \int _{0}^{T}\|\mathrm {u}\|_{H^{3}}^{2}+\| \mathrm {H}\|_{H^{3}}^{2}\,\mathrm{d} t \leq C(T). \end{aligned}$$

(3.13)

Proof

The proof of this lemma is standard; details are omitted. □

4 Proof of Theorems 1.1 and 1.2

With all the a priori estimates in Section 3 at hand, the existence and uniqueness of the strong solutions to (1.1)-(1.6) can be done by standard continuous arguments. Thus, this section is mainly devoted to the proof of the large-time behavior as (1.9)-(1.10) by the uniform estimates in Section 3. The proof is split into three steps as follows.

Proof of Theorem 1.1

Step 1: Decay of u and H.

Set

$$\begin{aligned} A_{1}(t)= \int\rho|\mathrm {u}|^{2}\,\mathrm {d}x,\qquad B_{1}(t)= \int|\mathrm {H}|^{2}\,\mathrm {d}x. \end{aligned}$$

From Lemmas 2.1-2.2 and the energy inequality, we have

$$\begin{aligned}& \int_{0}^{t}A_{1}(s)\,\mathrm {d}s\leq C \int_{0}^{t}\|\mathrm {u}\|^{2}_{L^{2}}\, \mathrm {d}s\leq C \int _{0}^{t}\|\nabla \mathrm {u}\|^{2}_{L^{2}} \,\mathrm {d}s\leq C, \end{aligned}$$

(4.1)

$$\begin{aligned}& \int _{0}^{t}B_{1}(s)\,\mathrm {d}s\leq C \int_{0}^{t}\|\nabla \mathrm {H}\|^{2}_{L^{2}} \,\mathrm {d}s\leq C. \end{aligned}$$

(4.2)

Next, we need to prove that $\frac{\mathrm {d}}{\mathrm {d}t}A_{1}$ and $\frac{\mathrm {d}}{\mathrm {d}t}B_{1}$ are in $L^{1}(0,\infty)$. From (1.2)-(1.3) and using integration by parts over Ω, we have

$$\begin{aligned}& \begin{aligned}[b] \frac{\mathrm {d}}{\mathrm {d}t}A_{1}(t)&= \int\rho_{t}|\mathrm {u}|^{2}+2\rho \mathrm {u}_{t}\cdot \mathrm {u}\,\mathrm {d}x\\ &=-2 \int|\nabla \mathrm {u}|^{2}\,\mathrm {d}x+2 \int H\cdot\nabla H\cdot u\,\mathrm {d}x+2 \int\rho \mathrm {f}\cdot \mathrm {u}\,\mathrm {d}x, \end{aligned} \end{aligned}$$

(4.3)

$$\begin{aligned}& \frac{\mathrm {d}}{\mathrm {d}t}B_{1}(t)= 2 \int \mathrm {H}_{t} \cdot \mathrm {H}\,\mathrm {d}x=-2 \int |\nabla \mathrm {H}|^{2}\,\mathrm {d}x+ \int \mathrm {H}\cdot\nabla \mathrm {u}\cdot \mathrm {H}\,\mathrm {d}x. \end{aligned}$$

(4.4)

For $A_{1}(t)$, we have by (4.3) and the lemmas in Section 3

$$\begin{aligned} \int_{0}^{t}\biggl\vert \frac{\mathrm {d}}{\mathrm {d}s} \bigl[A_{1}(s)\bigr]\biggr\vert \,\mathrm {d}s \leq& C \int _{0}^{t}\|\nabla \mathrm {u}\|^{2}_{L^{2}} \,\mathrm {d}s+C \int_{0}^{t}\biggl\vert \int \mathrm {H}\cdot\nabla \mathrm {H}\cdot \mathrm {u}\,\mathrm {d}x\biggr\vert \,\mathrm {d}s \\ &{}+C \int_{0}^{t}\biggl\vert \int\rho \mathrm {f}\cdot \mathrm {u}\,\mathrm {d}x\biggr\vert \,\mathrm {d}s \\ \leq& C+C \int_{0}^{t}\biggl\vert \int \mathrm {H}\cdot\nabla \mathrm {H}\cdot \mathrm {u}\,\mathrm {d}x\biggr\vert \,\mathrm {d}s \\ \leq& C, \end{aligned}$$

(4.5)

where we use the following estimate:

$$\begin{aligned} &\int_{0}^{t}\biggl\vert \int \mathrm {H}\cdot\nabla \mathrm {H}\cdot \mathrm {u}\,\mathrm {d}x\biggr\vert \,\mathrm {d}s\\ &\quad\leq C \int_{0}^{t}\|\nabla \mathrm {H}\|_{L^{2}}\|\mathrm {H}\|_{L^{4}}\|\mathrm {u}\|_{L^{4}}\,\mathrm {d}s \\ &\quad\leq C \int_{0}^{t}\|\nabla \mathrm {H}\|_{L^{2}}\bigl(\|\mathrm {H}\|_{L^{2}}+\| \mathrm {H}\|^{1/2}_{L^{2}}\|\nabla \mathrm {H}\|^{1/2}_{L^{2}}\bigr)\|\nabla \mathrm {u}\|_{L^{2}}\,\mathrm {d}s \\ &\quad\leq C \int_{0}^{t}\|\nabla \mathrm {H}\|^{2}_{L^{2}} \|\nabla \mathrm {u}\| _{L^{2}}\,\mathrm {d}s \\ &\quad\leq C \int_{0}^{t}\|\nabla \mathrm {H}\|^{2}_{L^{2}} \,\mathrm {d}s\leq C, \end{aligned}$$

due to (2.4), (2.6), and (3.3).

Similarly, for $B_{1}(t)$ we have

$$\begin{aligned} \int_{0}^{t}\biggl\vert \frac{\mathrm {d}}{\mathrm {d}t} \bigl[B_{1}(s)\bigr]\biggr\vert \,\mathrm {d}s \leq& C \int _{0}^{t}\|\nabla \mathrm {H}\|^{2}_{L^{2}} \,\mathrm {d}s+ \int_{0}^{t} \int|\mathrm {H}|^{2}|\nabla \mathrm {u}|\,\mathrm {d}x\,\mathrm{d}s \\ \leq& C+C \int_{0}^{t} \int|\mathrm {H}|^{2}|\nabla \mathrm {u}|\,\mathrm {d}x\,\mathrm {d}s \\ \leq& C. \end{aligned}$$

(4.6)

From (4.1), (4.2), (4.5), (4.6), we prove that $A_{1}(t), B_{1}(t)\in W^{1,1}(0,\infty)$. Therefore, we conclude that

$$ \lim_{t\to\infty} \bigl(\|\mathrm {u}\|_{L^{2}}+\|\mathrm {H}\|_{L^{2}} \bigr)=0, $$

due to $\rho> \underline{\rho}$.

Step 2: Decay of $\|\nabla \mathrm {u}\|_{L^{2}}$ and $\|\nabla \mathrm {H}\|_{L^{2}}$.

Set

$$\begin{aligned} A_{2}(t)= \int|\nabla \mathrm {u}|^{2}\,\mathrm {d}x,\qquad B_{2}(t)= \int |\nabla \mathrm {H}|^{2}\,\mathrm {d}x. \end{aligned}$$

By (3.2), we see that $A_{2}$ and $B_{2}$ are in $L^{1}(0,\infty)$. Next, it remains to prove that $\frac{\mathrm {d}}{\mathrm {d}t}A_{2}$ and $\frac{\mathrm {d}}{\mathrm {d}t}B_{2}$ are in $L^{1}(0,\infty)$. Using the Young inequality, (3.2) and (3.9), we have

$$\begin{aligned}& \begin{aligned}[b] \int_{0}^{t}\biggl\vert \frac{\mathrm {d}}{\mathrm {d}t}A_{2}(s) \biggr\vert \,\mathrm {d}s&\leq C \int_{0}^{t}\|\nabla \mathrm {u}\|_{L^{2}}\|\nabla \mathrm {u}_{t}\|_{L^{2}}\,\mathrm {d}s\\ &\leq C \int_{0}^{t}\|\nabla \mathrm {u}\|^{2}_{L^{2}}+ \|\nabla \mathrm {u}_{t}\|^{2}_{L^{2}}\,\mathrm {d}s\leq C, \end{aligned} \end{aligned}$$

(4.7)

$$\begin{aligned}& \begin{aligned}[b] \int_{0}^{t}\biggl\vert \frac{\mathrm {d}}{\mathrm {d}t}B_{2}(s) \biggr\vert \,\mathrm {d}s&\leq C \int_{0}^{t}\| \nabla \mathrm {H}\|_{L^{2}}\|\nabla \mathrm {H}_{t}\|_{L^{2}}\,\mathrm {d}\\ &\leq C \int_{0}^{t}\|\nabla \mathrm {H}\|^{2}_{L^{2}}+ \|\nabla \mathrm {H}_{t}\|^{2}_{L^{2}}\,\mathrm {d}s\leq C. \end{aligned} \end{aligned}$$

(4.8)

Therefore, $A_{2}(t), B_{2}(t)\in W^{1,1}(0,\infty)$. We conclude that

$$ \lim_{t\to\infty} \bigl(\|\nabla \mathrm {u}\|_{L^{2}}+\|\nabla \mathrm {H}\|_{L^{2}} \bigr)=0. $$

Step 3: Decay of $\|\mathrm {u}_{t}\|_{L^{2}}$ and $\|\mathrm {H}_{t}\| _{L^{2}}$.

Set

$$\begin{aligned} A_{3}(t)= \int\bigl|\rho^{1/2}\mathrm {u}_{t}\bigr|^{2}\,\mathrm {d}x, \qquad B_{3}(t)= \int |\mathrm {H}_{t}|^{2}\,\mathrm {d}x. \end{aligned}$$

From Lemma 3.1, we see that $A_{3}(t)$ and $B_{3}(t)$ are in $L^{1}(0,\infty)$. By the Poincaré inequality and (3.10)-(3.11), we have

$$\begin{aligned} \int_{0}^{t}\biggl\vert \frac{\mathrm {d}}{\mathrm {d}s}A_{3}(s) \biggr\vert \,\mathrm {d}s \leq& \int_{0}^{t}\|\nabla \mathrm {u}_{t} \|_{L^{2}}^{2}\,\mathrm {d}s+ \int _{0}^{t}|R_{1}|+|R_{2}|+|R_{3}|+|R_{4}|+|R_{5}| \,\mathrm {d}s \\ &{}+ \int_{0}^{t} \int|\mathrm {H}_{t}||\nabla \mathrm {H}||\mathrm {u}_{t}|+|\mathrm {H}||\nabla \mathrm {H}_{t}||\mathrm {u}_{t}|\,\mathrm {d}x\,\mathrm {d}s \\ \leq& C+C \int_{0}^{t}\|\mathrm {H}_{t}\|_{L^{4}}\| \mathrm {u}_{t}\|_{L^{4}}\| \nabla \mathrm {H}\|_{L^{2}}+\|\nabla \mathrm {H}_{t}\|_{L^{2}}\|\mathrm {H}\|_{L^{4}}\|\mathrm {u}_{t}\| _{L^{4}}\,\mathrm {d}s \\ \leq& C \end{aligned}$$

(4.9)

and

$$\begin{aligned} \int_{0}^{t}\biggl\vert \frac{\mathrm {d}}{\mathrm {d}s}B_{3}(s) \biggr\vert \,\mathrm {d}s \leq& C+ \int_{0}^{t}|R_{6}|+|R_{7}|\, \mathrm {d}s+ \int_{0}^{t} \int \mathrm {H}\cdot\nabla u_{t}\cdot \mathrm {H}_{t}\,\mathrm {d}x\,\mathrm {d}s \\ \leq& C+ \int_{0}^{t}\|\nabla \mathrm {u}_{t} \|_{L^{2}}\|\mathrm {H}\|_{L^{4}}\|\mathrm {H}_{t}\|_{L^{4}}\,\mathrm {d}s \\ \leq& C. \end{aligned}$$

(4.10)

Therefore, $A_{3}(t), B_{3}(t)\in W^{1,1}(0,\infty)$, which, together with $\rho> \underline{\rho}$, leads to

$$ \lim_{t\to\infty} \bigl(\|\mathrm {u}_{t}\|_{L^{2}}+\| \mathrm {H}_{t}\|_{L^{2}} \bigr)=0. $$

This completes the proof of Theorem 1.1. □

Proof of Theorem 1.2

Indeed, for the case $f=0$, all the estimates in Section 3 are still valid. Thus, the existence and uniqueness follow the same method. Here, we address the decay rate. First, the continuity equation (1.1) and the fact $\operatorname{div}u=0$ lead to, for all $p\in[0,\infty]$,

$$ \|\rho\|_{L^{p}}=\|\rho_{0}\|_{L^{p}}. $$

(4.11)

Particularly, one has

$$ \underline{\rho}\leq\|\rho\|_{L^{\infty}}=\|\rho_{0} \|_{L^{\infty}}\leq\bar{\rho}. $$

(4.12)

Multiplying (1.2) and (1.3) by u and H, respectively, then adding the two resulting equations, one has after integrating by parts over Ω

$$ \frac{\mathrm {d}}{\mathrm {d}t} \int\rho|\mathrm {u}|^{2}+|\mathrm {H}|^{2}\,\mathrm {d}x +2 \bigl(\| \nabla \mathrm {u}\| _{L^{2}}^{2}+\|\nabla \mathrm {H}\|_{L^{2}}^{2} \bigr)\leq0. $$

(4.13)

Because of the boundary condition $\mathrm {u}|_{\partial\Omega}=0$ and $\mathrm {H}\cdot n|_{\partial\Omega}=0$, using the Poincaré inequality, (2.4), and (2.6), we have

$$ \bigl(\|\mathrm {u}\|_{L^{2}}^{2}+\| \mathrm {H}\|_{L^{2}}^{2} \bigr)\leq C \bigl(\| \nabla \mathrm {u}\|_{L^{2}}^{2}+\|\nabla \mathrm {H}\|_{L^{2}}^{2} \bigr). $$

(4.14)

From (4.12), (4.13), and (4.14), it is easy to prove that there exists a $\theta_{0}>0$ such that

$$ \frac{\mathrm {d}}{\mathrm {d}t} \int\rho|\mathrm {u}|^{2}+|\mathrm {H}|^{2}\,\mathrm {d}x + \theta_{0} \bigl(\bigl\| \rho ^{1/2} \mathrm {u}\bigr\| _{L^{2}}^{2}+ \| \mathrm {H}\|_{L^{2}}^{2} \bigr)\leq0, $$

which immediately implies

$$ \bigl\| \rho^{1/2} \mathrm {u}\bigr\| _{L^{2}}^{2}+\| \mathrm {H}\|_{L^{2}}^{2}\leq\exp\{-\theta_{0}t\}. $$

(4.15)

Notice that ρ is bounded from below, and the decay rate estimate (1.11) is established. Thus, Theorem 1.2 is completed. □