Local Existence for Boussinesq Equations with Slip Boundary Condition in a Bounded Domain

Jian Lu

College of Science, University of Shanghai for Science and Technology, Shanghai, China.

DOI: 10.4236/jamp.2017.510165                

Abstract

In this paper, we are interested in the local existence for the Boussinesq equations with the slip boundary conditions. Energy method and Gerlakin approach are employed in this paper to get the main result.

Keywords

Boussinesq, Slip Boundary, Local Existence

1. Introduction

The Boussinesq equations are as follows:

$\{\begin{array}{l}{u}_t+\left(u\cdot \nabla \right)u+\nabla p=\gamma \Delta u+\theta f,\hfill \\ {\theta }_t+\left(u\cdot \nabla \right)\theta =\epsilon \Delta \theta ,\hfill \\ divu=0,\hfill \end{array}\left(x,t\right)\in \Omega \times \left(0,t\right)$ (1.1)

where $\Omega$ is a bounded smooth domain of $R^3$ , $u=\left(u^1,u^2,u^3\right)$ and $\theta$ represent density and temperature, $p$ is the pressure function, $f$ is the external force, $\gamma \ge 0,\epsilon \ge 0$ represent viscous coefficient and thermal conductivity coefficient.

In this paper, the initial data is given by

${\left(u,\theta \right)|}_{t=0}=\left(u_0,{\theta }_0\right)$ (1.2)

and the boundary condition is

${u\cdot n|}_{\partial \Omega }=0$ , ${\frac{\partial \theta }{\partial n}|}_{\partial \Omega }=\alpha \theta$ , ${{\left(Su\cdot n\right)}_{\tau }|}_{\partial \Omega }={\left(Bu\right)}_{\tau }$ (1.3)

Boussinesq equations are the classical model of fluid mechanics. There are a lot of important applications in marine ecology and weather forecasting. There are a lot of related conclusions about 2-D Boussinesq equations. Ye Zhuan [1] studied the Cauchy problem of the two-dimensional (2D) incompressible Boussinesq equations with fractional dissipation by making use of the nonlinear lower bounds for the fractional Laplacian established in Constantin and Vicol. Xiaojing Xu [2] studied the Cauchy problem of the two-dimensional (2D) incompressible Boussinesq equations with fractional Laplacian dissipation by the nonlinear lower bounds for the fractional Laplacian. In this paper, we want to study the 3-D Boussinesq equations. On the one hand, the regularity and well-posedness of Boussinesq equations is a popular problem that people study. Jishan fan [3] proved a regularity criterion for the 3D Boussinesq system with partial viscosity. T. Hmidi [4] studied the global well-posedness of the Euler-Boussinesq system with the term dissipation |D| α on the temperature equation. On the other hand, the existence of Boussinesq equations is always an important problem people are interested in. Ying Liu [5] applied the Fourier decomposition method to study the attenuation in $L^2$ of the weak existence of the Boussinesq equations. Wei Li [6] used the homogeneous balance method and travelling-wave transformation to acquire some exact solutions of the Boussinesq equations. Xianjin Li [7] studied the global stability in $L^2$ of the Boussinesq equations in three dimensional regional and unbounded regional.

However, for the initial questions and the boundary questions, the stations are more complicated and challenging. In fact, in fluid mechanics, two boundary conditions are considered mainly. One is the Diriclet boundary condition: $u=0,x\in \partial \Omega$ , and another is the famous boundary condition proposed by Navier:

$u\cdot n=0,{\left(Su\cdot n\right)}_{\tau }=-\alpha u_{\tau },x\in \partial \Omega$ (1.4)

where $n$ is the unit outward normal on $\partial \Omega$ , $u_{\tau }$ is the tangential part of $u$ , $Su$ denotes the deformation tensor: $Su=\frac{1}{2}\left(\nabla u+\nabla u^t\right)$ .

Now, what we need to do is the study for the local existence of the problem (1.1) - (1.3). Usually, the followed existences are considered. First, the smooth of the existence is so small [8] [9] [10] that the initial data and the existence are close to a constant in $H^2\left(\Omega \right)$ , however this existence does not have singularity. Second, the existence is the “large energy” proposed by Lions [11] . This existence has regularity but the analysis of the characteristic are more difficult and the uniqueness and the continuous dependence can’t be solved [12] [13] . Last, according to documents [14] [15] [16] , the initial data in $L^2\left(\Omega \right)$ is small and the initial density is positive and bounded. For example, when the initial data is piecewise smooth, the solutions shown in [17] satisfy the Rankine-Hugoniot conditions in a strict point wise sense. On the other hand, these solutions have so enough structure and regularity that the uniqueness and continuous dependence theory can be proved in [18] . The well-posedness theory is nearly complete in the whole space, but for the Dirichlet condition: on $\partial \Omega$ , the station $u=0$ may have problems. However, if considering the Navier condition described above, the global weak small-energy solutions can be proved to exist in a half space for initial data with small energy and bounded density. David Hoff [19] studied the equations

$\{\begin{array}{l}{\rho }_{\text{t}}+div\left(\rho u\right)=0,\hfill \\ {\left(\rho u^j\right)}_t+div\left(\rho u^{j}u\right)+P{\left(\rho \right)}_{x_j}=\mu \Delta u^j+\lambda div{u}_{x_j}+\rho f^j,\hfill \end{array}$ (1.5)

where $div(\cdot )$ denotes divergence.

Satisfying the bounded conditions:

$\{\begin{array}{l}u\cdot \upsilon =0,\hfill \\ \mu \prod \nabla u\upsilon =\prod Bu,\hfill \end{array}x\in \Omega ,t>0$ (1.6)

where $\upsilon \left(x\right)$ is the unit outward normal on $\partial \Omega$ , $\prod \left(x\right)$ is the projection onto the tangent plant to $\partial \Omega$ at x, then the existence of the solutions can be given. T. Hmidi [20] studied a fractional diffusion Boussinesq model which couples a Navier-Stokes type equation with fractional diffusion for the velocity and a transport equation for the temperature by establishing global well-posedness results with rough initial data.

In this paper, we consider the Boussinesq equations. Because the temperature and the fluid are coupled together, the study becomes more difficult. The key we solve the problem is how to deal with the temperature.

First, we define a function

$W=\left\{w\in {\left(H^1\left(\Omega \right)\right)}^3:w\left(x\right)\cdot n\left(x\right)=0,x\in \partial \Omega \right\}$ (1.7)

Definition 1.1. For a fixed time $T$ , $\left(u,\theta \right)$ is called the solutions for (1.1) - (1.3) on $\left[0,T\right]$ if the following holds:

1) $\left(u,\theta \right)\in C\left(\left[0,T\right];H^1\left(\Omega \right)\right)\cap L^2\left(\left[0,T\right];W^{2,q}\left(\Omega \right)\right)$ ,

$\left(u_t,{\theta }_t\right)\in C\left(\left[0,T\right];L^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T\right];H^1\left(\Omega \right)\right)$ .

2) For any times $t_1,t_2\in \left[0,T\right]$ ,

${\displaystyle {\int }_{t_1}^{t_2}{\displaystyle {\int }_{\Omega }\left(u_t+\left(u\cdot \nabla \right)u+\nabla p\right)\phi \text{d}x}\text{d}t}={\displaystyle {\int }_{t_1}^{t_2}{\displaystyle {\int }_{\Omega }\left(\gamma \Delta u+\theta f\right)\phi \text{d}x\text{d}t}}$ ,

${\displaystyle {\int }_{t_1}^{t_2}{\displaystyle {\int }_{\Omega }\left({\theta }_t+\left(u\cdot \nabla \right)u\right)\varphi \text{d}x}\text{d}t}={\displaystyle {\int }_{t_1}^{t_2}{\displaystyle {\int }_{\Omega }\left(\epsilon \Delta \theta \right)\varphi \text{d}x}\text{d}t}$ ,

for $\phi ,\varphi$ Lipschitz $\Omega \times \left[0,T\right]$ with $\varphi ,\phi \in W$ .

Remark: the system parameters $\Omega$ and $B$ will be assumed to satisfy the following conditions:

1) $\Omega$ is a bounded open set in $R^3$ with a $C^3$ boundary.

2) $B$ is a $C^1$ in a neighborhood of $\partial \Omega$ .

Theorem 1.2. Assume the hypotheses hold, let $q\in \left(3,6\right]$ be fixed constant and the $\left(u_0,{\theta }_0\right)$ , $p$ and $f$ satisfy the following conditions:

$\left(p,f\right)\in C\left(\left[0,T\right];L^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T\right];L^q\left(\Omega \right)\right)$ ;

$\left(p_t,f_t\right)\in L^2\left(\left[0,T\right];L^2\left(\Omega \right)\right)$ ;

${{\left(Su\cdot n\right)}_{\tau }|}_{\partial \Omega }={\left(Bu\right)}_{\tau }$ .

Then there is a small $T_{\ast }>0$ and a unique strong solution $\left(u,\theta \right)$ to the initial boundary value problems (1.2), (1.3) such that:

$\left(u,\theta \right)\in C\left(\left[0,T_{*}\right];H^1\left(\Omega \right)\right)\cap L^2\left(\left[0,T_{*}\right];W^{2,q}\left(\Omega \right)\right)$ ;

$\left(u_t,{\theta }_t\right)\in C\left(\left[0,T_{*}\right];L^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T_{*}\right];H^1\left(\Omega \right)\right)$ .

2. The Lame Operator and the Regularity of −D Operator

In this part, we introduce the regularity of Laplace and the $-\Delta$ . First we assume the condition (1.2) holds, $\Omega$ and $B$ satisfy the above hypotheses. Considering the following problem: find $u:\bar{\Omega }\to R^3$ such that:

$\{\begin{array}{l}Lu+\beta u=g,x\in \Omega \\ u\cdot n=0,\\ {\left(S\left(u\right)\cdot n\right)}_{\tau }={\left(Bu\right)}_{\tau },x\in \partial \Omega \end{array}$ (2.1)

here $Lu=-\gamma \Delta u$ , fixed $\beta \in R$ , given $g:\Omega \to R$ . Then the corresponding weak form is achieved.

Use $\omega \in W$ to multiply by the above differential Equation (2.1), we can get:

$A_{\beta }\left(u,\omega \right)={\displaystyle {\int }_{\Omega }g\cdot \omega \text{d}x}$ (2.2)

where $A_{\beta }:W\times W\to R$ is the bilinear form:

$A_{\beta }\left(u,\omega \right)=\gamma \left(-{\displaystyle {\int }_{\partial \Omega }Au\cdot \omega \text{d}s}+{\displaystyle {\int }_{\Omega }Su:S\omega \text{d}x}\right)$ (2.3)

Obviously, $A_{\beta }$ is continuous on $W\times W$ if $B$ is bounded, and we can use trace theorem to show that $A_{\beta }$ is coercive if $\beta$ is enough large depending on $\gamma$ , $\Omega$ , $B$ . In this case, there is a bounded operator $S_{\beta }:L^2\to W$ satisfying

$A_{\beta }\left(S_{\beta }g,\omega \right)={\displaystyle {\int }_{\Omega }g\cdot \omega \text{d}x}$ (2.4)

for all $\omega \in W$ . Furthermore, because the embedding $W\to L^2$ is compact, $S_{\beta }$ is a compact operator from $L^2$ to $L^2$ and the symmetry condition guarantees $S_{\beta }$ is self-adjoint.

The following lemmas are taken from document [21] .

Lemma 2.1. Assume that $\Omega \subset R^3$ is abounded open set with a $C^3$ boundary and that $B\in L^{\infty }\left(\partial \Omega \right)$ is a symmetric matrix. For $\beta$ enough large, there is a compact self-adjoint operator $S_{\beta }:L^2\left(\Omega \right)\to L^2\left(\Omega \right)$ , whose range is contained in $W$ and for which (2.2) holds for $g\in L^2\left(\Omega \right)$ and $\omega \in W$ . At the same time, there is an orthogonal basis ${\left\{{\omega }_k\right\}}_k$ for $L^2$ whose elements are in $W\cap C^{\infty }\left(\Omega \right)$ , and which are eigenfunctions of $S_{\beta }$

$S_{\beta }{\omega }_k={\gamma }_k^{-1}{\omega }_k$

where ${\gamma }_k\to \infty$ .

Lemma 2.2. Assume that $m\ge 0$ , $\Omega \subset R^3$ is abounded open set with a $C^{m+1}$ boundary and that B is a $C^{m+1}\left(\partial \Omega \right)$ mapping from a neighborhood of $\partial \Omega$ into the set of $3\times 3$ matrices, then there exists a constant $C_m=C\left(\gamma ,\Omega ,B\right)$ such that if $u$ is a solution of (2.1) in the sense of (2.2) where $\beta$ and $g\in H^m\left(\Omega \right)$ , then $u\in H^{m+2}\left(\Omega \right)$ and

${\Vert u\Vert }_{H^{m+2}\left(\Omega \right)}\le C_m\left({\Vert u\Vert }_{L^2\left(\Omega \right)}+{\Vert g\Vert }_{H^m\left(\Omega \right)}\right)$ (2.5)

For the operator $-\Delta$ , consider the following problem: find $\theta :\bar{\Omega }\to R^3$ such that

$\{\begin{array}{l}-\Delta \theta +\delta \theta =\bar{g},x\in \Omega \hfill \\ {\frac{\partial \theta }{\partial n}|}_{\partial \Omega }=\alpha \theta \hfill \end{array}$ (2.6)

where $\delta \in R$ is fixed, $\bar{g}:\Omega \to R$ is given.

Use $\varphi \in H^1\left(\Omega \right)$ to multiply by the above differential equation, and integration by parts, then we get

$B_{\delta }\left(\theta ,\varphi \right)={\displaystyle {\int }_{\Omega }\bar{g}\cdot \varphi \text{d}x}$ (2.7)

where $B_{\delta }:H^1\times H^1\to R$ is the bilinear form:

$B_{\delta }\left(\theta ,\varphi \right)={\displaystyle {\int }_{\Omega }\nabla \theta \cdot \nabla \phi +\delta \theta \phi \text{d}x}$ (2.8)

Similar to the Lame operator, there exist a bounded operator ${\Gamma }_{\delta }:L^2\left(\Omega \right)\to H^1$ meeting

$B_{\delta }\left({\Gamma }_{\delta }g,\varphi \right)={\displaystyle {\int }_{\Omega }\bar{g}\cdot \varphi \text{d}x}$ (2.9)

for all $\phi \in H^1$ and for $\delta$ large enough. Parallel to Lemmas 2.1 - 2.2, we can get the following lemmas.

Lemma 2.3. Assume that $\Omega \subset R^3$ is abounded open set with a $C^3$ boundary. For $\gamma$ large enough, there exist a compact self-adjoint operator ${\Gamma }_{\gamma }:L^2\left(\Omega \right)\to L^2\left(\Omega \right)$ , whose range is contained in $H^1$ and for which (2.9) holds for $\bar{g}\in L^2\left(\Omega \right)$ and $\varphi \in H^1$ . At the same time, there is an orthogonal basis ${\left\{{\varphi }_k\right\}}_k$ for $L^2$ whose elements are in $H^1\cap C^{\infty }\left(\Omega \right)$ and which are eigenfunctions of ${\Gamma }_{\gamma }$ , ${\Gamma }_{\gamma }{\varphi }_k=\overline{{\gamma }_k^{-1}}{\omega }_k$ , where $\overline{{\gamma }_k}\to \infty$ .

Lemma 2.4. Assume that $m\ge 0$ , $\Omega \subset R^3$ is abounded open set with a $C^{m+1}$ boundary and that B is a $C^{m+1}\left(\partial \Omega \right)$ mapping from a neighborhood of $\partial \Omega$ into the set of $3\times 3$ matrices, then there exists a constant $\overline{C_m}=\bar{C}\left(\Omega \right)$ such that if $\theta$ is a solution of (2.6) in the sense of (2.7) where $\bar{g}\in H^m\left(\Omega \right)$ , then $\theta \in H^{m+2}\left(\Omega \right)$ and

${\Vert \theta \Vert }_{H^{m+2}\left(\Omega \right)}\le \overline{C_m}\left({\Vert \theta \Vert }_{L^2\left(\Omega \right)}+{\Vert \bar{g}\Vert }_{H^m\left(\Omega \right)}\right)$ (2.10)

3. A Prior Estimates for Higher Regularity

In this part, we need the following prior estimates to prove the local existence of the solution. Assume that the following inequalities hold:

$\Vert {\left(u_0,{\theta }_0\right)}_{H^2}\Vert +1\le C_0$ (3.1)

$\underset{0\le t\le T_{\ast }}{\sup }\left({\Vert u\left(t\right)\Vert }_{H^2}+{\Vert u_t\left(t\right)\Vert }_{L^2}\right)+{\displaystyle {\int }_0^{T_{\ast }}\left({\Vert \nabla u_t\Vert }_{L^2}^2+{\Vert {\nabla }^{2}u\Vert }_{L^q}^2\right)\text{d}t}+1\le C_1$ (3.2)

where $1\le C_0\le C_1$ , $0<T_{\ast }\le T$ .

Remark: C is a constant if be not added.

3.1. Estimate for Temperature

Lemma 3.1.1. ${\Vert \theta \Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla \theta \Vert }_{L^2}^2\text{d}s}\le C$

Proof: Multiplying the second equation of (1.1) by $\theta$ and integrating over $\Omega$ , one has

$\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }{\theta }^2\text{d}x}={\displaystyle {\int }_{\Omega }\epsilon \Delta \theta \cdot \theta -\left(u\cdot \nabla \right)\theta \cdot \theta \text{d}x}$ (3.1.1)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\epsilon \Delta \theta \cdot \theta \text{d}x}\le C\left({\displaystyle {\int }_{\partial \Omega }\frac{\partial \theta }{\partial n}\theta \text{d}s}-{\displaystyle {\int }_{\Omega }{\left(\nabla \theta \right)}^2\text{d}x}\right)\\ =C\left({\displaystyle {\int }_{\partial \Omega }\alpha {\theta }^2\text{d}s}-{\displaystyle {\int }_{\Omega }{\left(\nabla \theta \right)}^2\text{d}x}\right)\\ \le C{\displaystyle {\int }_{\Omega }\left({\theta }^2+\left|\theta \right|\cdot \left|\nabla \theta \right|\right)\text{d}x}+C{\displaystyle {\int }_{\Omega }{\left(\nabla \theta \right)}^2\text{d}x}\\ \le C\left({\Vert \theta \Vert }_{L^2}+{\Vert \nabla \theta \Vert }_{L^2}+C_{\delta }{\Vert \theta \Vert }_{L^2}+\delta {\Vert \nabla \theta \Vert }_{L^2}\right)\end{array}$ (3.1.2)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\left(u\cdot \nabla \right)\theta \cdot \theta \text{d}x}={\displaystyle {\int }_{\Omega }u\cdot \nabla \theta \cdot \theta }\text{d}x\le {\Vert u\Vert }_{L^3}{\Vert \theta \Vert }_{H^1}{\Vert \nabla \theta \Vert }_{L^2}\\ \le {\Vert u\Vert }_{L^2}^{\frac{1}{2}}{\Vert u\Vert }_{H^1}^{\frac{1}{2}}\left({\Vert \theta \Vert }_{L^2}+{\Vert \nabla \theta \Vert }_{L^2}\right){\Vert \nabla \theta \Vert }_{L^2}\\ \le C\left({\Vert \theta \Vert }_{L^2}{\Vert \nabla \theta \Vert }_{L^2}+{\Vert \nabla \theta \Vert }_{L^2}^2\right)\\ \le C\left(\delta {\Vert \nabla \theta \Vert }_{L^2}^2+C_{\delta }{\Vert \theta \Vert }_{L^2}^2+{\Vert \nabla \theta \Vert }_{L^2}^2\right)\end{array}$ (3.1.3)

Substituting (3.1.2) - (3.1.3) into (3.1.1), letting $\delta$ small enough and using Gronwall’s inequality:

${\Vert \theta \Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla \theta \Vert }_{L^2}^2\text{d}s}\le C$

Lemma 3.1.2. ${\Vert \nabla \theta \Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert {\theta }_t\Vert }_{L^2}^2\text{d}s}\le C$

Proof: Multiplying the second equation of (1.1) by ${\theta }_t$ and integrating over $\Omega$ , one has

${\displaystyle {\int }_{\Omega }{\theta }_t^2\text{d}x}={\displaystyle {\int }_{\Omega }\epsilon \Delta \theta \cdot {\theta }_t-\left(u\cdot \nabla \right)\theta \cdot {\theta }_t\text{d}x}$ (3.1.4)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\epsilon \Delta \theta \cdot {\theta }_t\text{d}x}={\displaystyle {\int }_{\partial \Omega }\epsilon \frac{\partial \theta }{\partial n}{\theta }_t\text{d}s}-{\displaystyle {\int }_{\Omega }\epsilon \nabla \theta \cdot \nabla {\theta }_t\text{d}x}\\ =\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\partial \Omega }\epsilon \alpha {\theta }^2}\text{d}s-\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }\epsilon {\left(\nabla \theta \right)}^2\text{d}x}\\ \le C\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }\left({\left|\theta \right|}^2+\left|\theta \right|\cdot \left|\nabla \theta \right|\right)+{\left(\nabla \theta \right)}^2}\text{d}x\\ \le C\frac{\text{d}}{\text{d}t}\left(C_{\delta }{\Vert \nabla \theta \Vert }_{L^2}^2+\delta {\Vert \theta \Vert }_{L^2}^2+{\Vert \nabla \theta \Vert }_{L^2}^2\right)\end{array}$ (3.1.5)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\left(u\cdot \nabla \right)\theta \cdot {\theta }_t\text{d}x}={\displaystyle {\int }_{\Omega }u\cdot \nabla \theta \cdot {\theta }_t\text{d}x}\le {\Vert u\Vert }_{L^{\infty }}{\Vert {\theta }_t\Vert }_{L^2}{\Vert \nabla \theta \Vert }_{L^2}\\ \le C\left(\delta {\Vert {\theta }_t\Vert }_{L^2}^2+C_{\delta }{\Vert \nabla \theta \Vert }_{L^2}^2\right)\end{array}$ (3.1.6)

Substituting (3.1.5) - (3.1.6) into (3.1.4), letting $\delta$ small enough and using Gronwall’s inequality,

${\Vert \nabla \theta \Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert {\theta }_t\Vert }_{L^2}^2\text{d}s}\le C$

Lemma 3.1.3. ${\Vert {\theta }_t\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla {\theta }_t\Vert }_{L^2}^2\text{d}s}\le C$

Proof: Differentiating the second equation of (1.1) with respect to $t$ , multiplying the second equation of (1.1) by ${\theta }_t$ and integrating over $\Omega$ , one has

$\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }{\theta }_t^2\text{d}x}={\displaystyle {\int }_{\Omega }\epsilon \Delta {\theta }_t\cdot {\theta }_t-u_t\cdot \nabla \theta \cdot {\theta }_t-u\cdot \nabla {\theta }_t\cdot {\theta }_t\text{d}x}$ (3.1.7)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\epsilon \Delta {\theta }_t\cdot {\theta }_t\text{d}x}={\displaystyle {\int }_{\partial \Omega }\epsilon \frac{\partial {\theta }_t}{\partial n}{\theta }_t\text{d}s}-{\displaystyle {\int }_{\Omega }\epsilon {\left(\nabla {\theta }_t\right)}^2\text{d}x}\\ \le C\left({\displaystyle {\int }_{\partial \Omega }{\theta }_t^2\text{d}s}-{\displaystyle {\int }_{\Omega }{\left(\nabla {\theta }_t\right)}^2\text{d}x}\right)\\ \le C\left({\displaystyle {\int }_{\Omega }{\theta }_t^2+\nabla {\theta }_t\cdot {\theta }_t\text{d}x}+{\Vert \nabla {\theta }_t\Vert }_{L^2}^2\right)\\ \le C\left({\Vert {\theta }_t\Vert }_{L^2}^2+{\Vert \nabla {\theta }_t\Vert }_{L^2}{\Vert {\theta }_t\Vert }_{L^2}+{\Vert \nabla {\theta }_t\Vert }_{L^2}^2\right)\\ \le C\left({\Vert {\theta }_t\Vert }_{L^2}^2+C_{\delta }{\Vert \nabla {\theta }_t\Vert }_{L^2}^2+\delta {\Vert {\theta }_t\Vert }_{L^2}^2+{\Vert \nabla {\theta }_t\Vert }_{L^2}^2\right)\end{array}$ (3.1.8)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }{u}_t\cdot \nabla \theta \cdot {\theta }_t\text{d}x}\le {\Vert u_t\Vert }_{L^3}{\Vert \nabla \theta \Vert }_{L^2}{\Vert {\theta }_t\Vert }_{H^1}\le C{\Vert u_t\Vert }_{L^2}^{\frac{1}{2}}{\Vert \theta \Vert }_{L^6}^{\frac{1}{2}}{\Vert {\theta }_t\Vert }_{H^1}\\ \le C\left(C_{\delta }{\Vert {\theta }_t\Vert }_{L^2}^2+\delta {\Vert u_t\Vert }_{L^2}{\Vert \theta \Vert }_{H^1}+C_{\kappa }{\Vert \nabla {\theta }_t\Vert }_{L^2}^2+\kappa {\Vert u_t\Vert }_{L^2}{\Vert \theta \Vert }_{H^1}\right)\end{array}$ (3.1.9)

${\displaystyle {\int }_{\Omega }u\cdot \nabla {\theta }_t\cdot {\theta }_t\text{d}x}\le {\Vert u\Vert }_{L^{\infty }}{\Vert {\theta }_t\Vert }_{L^2}{\Vert \nabla {\theta }_t\Vert }_{L^2}\le C_{\delta }{\Vert {\theta }_t\Vert }_{L^2}^2+\delta {\Vert \nabla {\theta }_t\Vert }_{L^2}^2$ (3.1.10)

Substituting (3.1.8) - (3.1.10) into (3.1.7), letting $\delta$ and $\kappa$ small enough and using Gronwall’s inequality,

${\Vert {\theta }_t\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla {\theta }_t\Vert }_{L^2}^2\text{d}s}\le C$

3.2. Estimate for Velocity

Lemma 3.2.1. ${\Vert u\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla u\Vert }_{L^2}^2\text{d}s}\le C$

Proof: Multiplying the first equation of (1.1) by $u$ and integrating over $\Omega$ , one has

$\frac{1}{2}\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }{u}^2\text{d}x}={\displaystyle {\int }_{\Omega }\gamma \Delta u\cdot u+u^2\cdot \nabla u-\nabla p\cdot u+\theta f\cdot u\text{d}x}$ (3.2.1)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\gamma \Delta u\cdot u\text{d}x}={\displaystyle {\int }_{\Omega }2\gamma div\left(Su\right)\cdot u\text{d}x}=2\gamma \left({\displaystyle {\int }_{\partial \Omega }Su\cdot n\cdot u\text{d}s}-{\displaystyle {\int }_{\Omega }Su\cdot \nabla u\text{d}x}\right)\\ \le C\left({\displaystyle {\int }_{\partial \Omega }Bu\cdot u\text{d}s}+{\displaystyle {\int }_{\Omega }Su\cdot \nabla u\text{d}x}\right)\le C\left({\Vert u\Vert }_{L^2}^2+{\Vert \nabla u\Vert }_{L^2}^2\right)\end{array}$ (3.2.2)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\left(u\cdot \nabla \right)u\cdot u\text{d}x}={\displaystyle {\int }_{\Omega }{u}^2\cdot \nabla u\text{d}x}\le {\Vert u\Vert }_{L^3}{\Vert u\Vert }_{H^1}{\Vert \nabla u\Vert }_{L^2}\\ \le {\Vert u\Vert }_{L^2}^{\frac{1}{2}}{\Vert u\Vert }_{H^1}^{\frac{1}{2}}\left({\Vert u\Vert }_{L^2}+{\Vert \nabla u\Vert }_{L^2}\right){\Vert \nabla u\Vert }_{L^2}\\ \le C\left({\Vert u\Vert }_{L^2}{\Vert \nabla u\Vert }_{L^2}+{\Vert \nabla u\Vert }_{L^2}^2\right)\\ \le C\left(\delta {\Vert \nabla u\Vert }_{L^2}^2+C_{\delta }{\Vert u\Vert }_{L^2}^2+{\Vert \nabla u\Vert }_{L^2}^2\right)\end{array}$ (3.2.3)

${\displaystyle {\int }_{\Omega }\nabla p\cdot u\text{d}x}={\Vert \nabla p\Vert }_{L^2}{\Vert u\Vert }_{L^2}\le C_{\delta }{\Vert u\Vert }_{L^2}^2+\delta {\Vert \nabla p\Vert }_{L^2}^2$ (3.2.4)

${\displaystyle {\int }_{\Omega }\theta f\cdot u\text{d}x}={\Vert \theta \Vert }_{L^2}{\Vert f\Vert }_{L^{\infty }}{\Vert u\Vert }_{L^2}\le C\left(C_{\delta }{\Vert u\Vert }_{L^2}^2+\delta {\Vert \theta \Vert }_{L^2}^2\right)$ (3.2.5)

Substituting (3.2.2) - (3.2.5) into (3.2.1), letting $\delta$ small enough and using Gronwall’s inequality,

${\Vert u\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla u\Vert }_{L^2}^2\text{d}s}\le C$

Lemma 3.2.2. ${\Vert \nabla u\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert u_t\Vert }_{L^2}^2\text{d}s}\le C$

Proof: Multiplying the first equation of (1.1) by $u_t$ and integrating over $\Omega$ , one has

${\displaystyle {\int }_{\Omega }{u}_t^2\text{d}x}={\displaystyle {\int }_{\Omega }\gamma \Delta u\cdot u_t+u\cdot \nabla u\cdot u_t-\nabla p\cdot u_t+\theta f\cdot u_t\text{d}x}$ (3.2.6)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\gamma \Delta u\cdot u_t\text{d}x}={\displaystyle {\int }_{\Omega }2\gamma div\left(Su\right)\cdot u_t\text{d}x}=2\gamma \left({\displaystyle {\int }_{\partial \Omega }Su\cdot n\cdot u_t\text{d}s}-{\displaystyle {\int }_{\Omega }Su\cdot \nabla u_t\text{d}x}\right)\\ \le C\left({\displaystyle {\int }_{\partial \Omega }Bu\cdot u_t\text{d}s}+{\displaystyle {\int }_{\Omega }Su\cdot \nabla u_t\text{d}x}\right)\le C\frac{\text{d}}{\text{d}t}\left({\Vert u\Vert }_{L^2}^2+{\Vert \nabla u\Vert }_{L^2}^2\right)\end{array}$ (3.2.7)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\left(u\cdot \nabla \right)u\cdot u_t\text{d}x}={\displaystyle {\int }_{\Omega }u\cdot \nabla u\cdot u_t\text{d}x}\le {\Vert u\Vert }_{L^{\infty }}{\Vert u_t\Vert }_{L^2}{\Vert \nabla u\Vert }_{L^2}\\ \le C\left(C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert \nabla u\Vert }_{L^2}^2\right)\end{array}$ (3.2.8)

${\displaystyle {\int }_{\Omega }\nabla p\cdot u_t\text{d}x}={\Vert \nabla p\Vert }_{L^2}{\Vert u_t\Vert }_{L^2}\le C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert \nabla p\Vert }_{L^2}^2$ (3.2.9)

${\displaystyle {\int }_{\Omega }\theta f\cdot u_{t}dx}={\Vert \theta \Vert }_{L^2}{\Vert f\Vert }_{L^{\infty }}{\Vert u_t\Vert }_{L^2}\le C\left(C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert \theta \Vert }_{L^2}^2\right)$ (3.2.10)

Substituting (3.2.7) - (3.2.10) into (3.2.6), letting $\delta$ small enough and using Gronwall’s inequality,

${\Vert \nabla u\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert u_t\Vert }_{L^2}^2\text{d}s}\le C$

Lemma 3.2.3. ${\Vert u_t\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla u_t\Vert }_{L^2}^2\text{d}s}\le C$

Proof: Differentiating the first equation of (1.1) with respect to $t$ , multiplying the first equation of (1.1) by $u_t$ and integrating over $\Omega$ , one has

$\begin{array}{l}{\displaystyle {\int }_{\Omega }{u}_{tt}{u}_t\text{d}x}\\ ={\displaystyle {\int }_{\Omega }\gamma \Delta u_t\cdot u_t+u_t\cdot \nabla u\cdot u_t+u\cdot \nabla u_t\cdot u_t-\nabla p_t\cdot u_t+{\theta }_{t}f\cdot u_t+\theta f_t\cdot u_t\text{d}x}\end{array}$ (3.2.11)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }\gamma \Delta u_t\cdot u_t\text{d}x}={\displaystyle {\int }_{\Omega }2\gamma div\left(S{u}_t\right)\cdot u_t\text{d}x}=2\gamma \left({\displaystyle {\int }_{\partial \Omega }S{u}_t\cdot n\cdot u_t\text{d}s}-{\displaystyle {\int }_{\Omega }S{u}_t\cdot \nabla u_t\text{d}x}\right)\\ \le C\left({\displaystyle {\int }_{\partial \Omega }B{u}_t\cdot u_t\text{d}s}+{\displaystyle {\int }_{\Omega }S{u}_t\cdot \nabla u_t\text{d}x}\right)\le C\left({\Vert u_t\Vert }_{L^2}^2+{\Vert \nabla u_t\Vert }_{L^2}^2\right)\end{array}$ (3.2.12)

$\begin{array}{c}{\displaystyle {\int }_{\Omega }{u}_t\cdot \nabla u\cdot u_t\text{d}x}\le {\Vert u_t\Vert }_{L^3}{\Vert \nabla u\Vert }_{L^2}{\Vert u_t\Vert }_{H^1}\le C{\Vert u_t\Vert }_{L^2}^{\frac{1}{2}}{\Vert u_t\Vert }_{L^6}^{\frac{1}{2}}{\Vert u_t\Vert }_{H^1}\\ \le C\left(C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert u_t\Vert }_{L^2}{\Vert u_t\Vert }_{H^1}+C_{\kappa }{\Vert \nabla u_t\Vert }_{L^2}^2+\kappa {\Vert u_t\Vert }_{L^2}{\Vert u_t\Vert }_{H^1}\right)\end{array}$ (3.2.13)

${\displaystyle {\int }_{\Omega }u\cdot \nabla u_t\cdot u_t\text{d}x}\le {\Vert u\Vert }_{L^{\infty }}{\Vert u_t\Vert }_{L^2}{\Vert \nabla u_t\Vert }_{L^2}\le C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert \nabla u_t\Vert }_{L^2}^2$ (3.2.14)

${\displaystyle {\int }_{\Omega }\nabla p_t\cdot u_t\text{d}x}={\Vert \nabla p_t\Vert }_{L^2}{\Vert u_t\Vert }_{L^2}\le C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert \nabla p_t\Vert }_{L^2}^2$ (3.2.15)

${\displaystyle {\int }_{\Omega }{\theta }_{t}f\cdot u_t\text{d}x}={\Vert {\theta }_t\Vert }_{L^2}{\Vert f\Vert }_{L^{\infty }}{\Vert u_t\Vert }_{L^2}\le C\left(C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert {\theta }_t\Vert }_{L^2}^2\right)$ (3.2.16)

${\displaystyle {\int }_{\Omega }\theta f_t\cdot u_t\text{d}x}={\Vert \theta \Vert }_{L^2}{\Vert f_t\Vert }_{L^{\infty }}{\Vert u_t\Vert }_{L^2}\le C\left(C_{\delta }{\Vert u_t\Vert }_{L^2}^2+\delta {\Vert \theta \Vert }_{L^2}^2\right)$ (3.2.17)

Substituting (3.2.12)-(3.2.17) into (3.2.11), letting $\delta$ and $\kappa$ small enough and using Gronwall’s inequality,

${\Vert u_t\Vert }_{L^2}^2+{\displaystyle {\int }_0^t{\Vert \nabla u_t\Vert }_{L^2}^2\text{d}s}\le C$

3.3. Elliptic Estimates for Velocity and Temperature

Lemma 3.3.1. ${\Vert \nabla \theta \Vert }_{H^1}\le C$

Proof: According to the second equation of (1.1), we have $-\Delta \theta =-\frac{{\theta }_t}{\epsilon }-\frac{\left(u\cdot \nabla \right)\theta }{\epsilon }$ , then use the elliptic regularity:

${\Vert \theta \Vert }_{H^2}\le C\left({\Vert -{\theta }_t-\left(u\cdot \nabla \right)\theta \Vert }_{L^2}+{\Vert \theta \Vert }_{L^2}\right)\le C\left({\Vert {\theta }_t\Vert }_{L^2}+{\Vert \nabla \theta \cdot u\Vert }_{L^2}+{\Vert \theta \Vert }_{L^2}\right)$

$\begin{array}{c}{\Vert \nabla \theta \Vert }_{H^1}\le {\Vert \theta \Vert }_{H^2}\le C\left({\Vert {\theta }_t\Vert }_{L^2}+{\Vert \nabla \theta \cdot u\Vert }_{L^2}+{\Vert \theta \Vert }_{L^2}\right)\\ \le C\left({\Vert {\theta }_t\Vert }_{L^2}+{\Vert u\Vert }_{L^{\infty }}{\Vert \nabla \theta \Vert }_{L^2}+{\Vert \theta \Vert }_{L^2}\right)\le C\end{array}$

Lemma 3.3.2. ${\displaystyle {\int }_0^t{\Vert \nabla \theta \Vert }_{W^{1,q}}^2\text{d}s}\le C$

Proof: According to the second equation of (1.1), we have $-\Delta \theta =-\frac{{\theta }_t}{\epsilon }-\frac{\left(u\cdot \nabla \right)\theta }{\epsilon }$ , then use the elliptic regularity:

${\Vert \theta \Vert }_{W^{2,q}}\le C\left({\Vert -{\theta }_t-\left(u\cdot \nabla \right)\theta \Vert }_{L^q}+{\Vert \theta \Vert }_{L^q}\right)\le C\left({\Vert {\theta }_t\Vert }_{L^q}+{\Vert \nabla \theta \cdot u\Vert }_{L^q}+{\Vert \theta \Vert }_{L^q}\right)$

Both sides of the above inequality multiply by itself, then

${\Vert \theta \Vert }_{W^{2,q}}^2\le C\left({\Vert {\theta }_t\Vert }_{L^q}^2+{\Vert \nabla \theta \cdot u\Vert }_{L^q}^2+{\Vert \theta \Vert }_{L^q}^2\right)$

$\begin{array}{c}{\Vert \nabla \theta \Vert }_{W^{1,q}}^2\le {\Vert \theta \Vert }_{W^{2,q}}^2\le C\left({\Vert {\theta }_t\Vert }_{L^q}^2+{\Vert \nabla \theta \cdot u\Vert }_{L^q}^2+{\Vert \theta \Vert }_{L^q}^2\right)\\ \le C\left({\Vert {\theta }_t\Vert }_{H^1}^2+{\Vert u\Vert }_{L^{\infty }}^2{\Vert \nabla \theta \Vert }_{H^1}^2+{\Vert \theta \Vert }_{H^1}^2\right)\end{array}$

And integrating over $\left(0,t\right)$ , we have

${\displaystyle {\int }_0^t{\Vert \nabla \theta \Vert }_{W^{1,q}}^2\text{d}s}\le {\displaystyle {\int }_0^{t}C\left({\Vert {\theta }_t\Vert }_{H^1}^2+{\Vert u\Vert }_{L^{\infty }}^2{\Vert \nabla \theta \Vert }_{H^1}^2+{\Vert \theta \Vert }_{H^1}^2\right)\text{d}s}\le C$

Lemma 3.3.3. ${\Vert \nabla u\Vert }_{H^1}\le C$

Proof: According to the first equation of (1.1), we have $-\Delta u=-\frac{u_t}{\gamma }-\frac{\left(u\cdot \nabla \right)u}{\gamma }-\frac{\nabla p}{\gamma }+\frac{\theta f}{\gamma }$ , then use the elliptic regularity:

$\begin{array}{c}{\Vert u\Vert }_{H^2}\le C\left({\Vert -u_t-\left(u\cdot \nabla \right)u-\nabla p+\theta f\Vert }_{L^2}+{\Vert u\Vert }_{L^2}\right)\\ \le C\left({\Vert u_t\Vert }_{L^2}+{\Vert \left(u\cdot \nabla \right)u\Vert }_{L^2}+{\Vert \nabla p\Vert }_{L^2}+{\Vert \theta f\Vert }_{L^2}+{\Vert u\Vert }_{L^2}\right)\end{array}$

$\begin{array}{c}\Vert \nabla u\Vert |{|}_{H^1}\le {\Vert u\Vert }_{H^2}\le C\left({\Vert u_t\Vert }_{L^2}+{\Vert \left(u\cdot \nabla \right)u\Vert }_{L^2}+{\Vert \nabla p\Vert }_{L^2}+{\Vert \theta f\Vert }_{L^2}+{\Vert u\Vert }_{L^2}\right)\\ \le C\left({\Vert u_t\Vert }_{L^2}+{\Vert u\Vert }_{L^3}{\Vert \nabla u\Vert }_{H^1}+{\Vert \nabla p\Vert }_{L^2}+\text{||}f{\text{||}}_{L^{\infty }}{\Vert \theta \Vert }_{L^2}+{\Vert u\Vert }_{L^2}\right)\\ \le C\left({\Vert u_t\Vert }_{L^2}+{\Vert u\Vert }_{L^3}^{\frac{\text{1}}{\text{2}}}{\Vert u\Vert }_{H^1}^{\frac{\text{1}}{\text{2}}}{\Vert \nabla u\Vert }_{H^1}+{\Vert \nabla p\Vert }_{L^2}+{\Vert f\Vert }_{L^{\infty }}{\Vert \theta \Vert }_{L^2}+{\Vert u\Vert }_{L^2}\right)\\ \le C\end{array}$

Lemma 3.3.4. ${\displaystyle {\int }_0^t{\Vert \nabla u\Vert }_{W^{1,q}}^2\text{d}s}\le C$

Proof: According to the first equation of (1.1), we have $-\Delta u=-\frac{u_t}{\gamma }-\frac{\left(u\cdot \nabla \right)u}{\gamma }-\frac{\nabla p}{\gamma }+\frac{\theta f}{\gamma }$ , then use the elliptic regularity:

$\begin{array}{c}{\Vert u\Vert }_{W^{2,q}}\le C\left({\Vert -u_t-\left(u\cdot \nabla \right)u-\nabla p+\theta f\Vert }_{L^q}+{\Vert u\Vert }_{L^q}\right)\\ \le C\left({\Vert u_t\Vert }_{L^q}+{\Vert \left(u\cdot \nabla \right)u\Vert }_{L^q}+{\Vert \nabla p\Vert }_{L^q}+{\Vert \theta f\Vert }_{L^q}+{\Vert u\Vert }_{L^q}\right)\end{array}$

Both sides of the above inequality multiply by itself, then

${\Vert u\Vert }_{W^{2,q}}^2\le C\left({\Vert u_t\Vert }_{L^q}^2+{\Vert \left(u\cdot \nabla \right)u\Vert }_{L^q}^2+{\Vert \nabla p\Vert }_{L^q}^2+{\Vert \theta f\Vert }_{L^q}^2+{\Vert u\Vert }_{L^q}^2\right)$

$\begin{array}{c}{\Vert \nabla u\Vert }_{W^{1,q}}^2\le {\Vert u\Vert }_{W^{2,q}}^2\le C\left({\Vert u_t\Vert }_{L^q}^2+{\Vert \left(u\cdot \nabla \right)u\Vert }_{L^q}^2+{\Vert \nabla p\Vert }_{L^q}^2+{\Vert \theta f\Vert }_{L^q}^2+{\Vert u\Vert }_{L^q}^2\right)\\ \le C\left({\Vert u_t\Vert }_{H^1}^2+{\Vert u\Vert }_{L^{\infty }}^2{\Vert \nabla u\Vert }_{H^1}^2+{\Vert \nabla p\Vert }_{L^q}^2+{\Vert f\Vert }_{L^{\infty }}^2{\Vert \theta \Vert }_{H^1}^2+{\Vert u\Vert }_{H^1}^2\right)\end{array}$

And integrating over $\left(0,t\right)$ , we have ${\displaystyle {\int }_0^t{\Vert \nabla u\Vert }_{W^{1,q}}^2\text{d}s}\le C$ .

4. The Local Existence of the Solution of Boussinesq Equations

First, we consider the following linearized system:

$\{\begin{array}{l}{u}_t+\left(v\cdot \nabla \right)u+\nabla p=\gamma \Delta u+\theta f,\hfill \\ {\theta }_t+\left(u\cdot \nabla \right)\theta =\epsilon \Delta \theta ,\hfill \\ divu=0,\hfill \end{array}\left(x,t\right)\in \Omega \times \left(0,t\right)$ (4.1)

Lemma 4.1. Let $\Omega$ be a bounded domain in $R^3$ with smooth boundary, when $3<q\le 6$ , we have $\left(u_0,{\theta }_0\right)\in H^1\left(\Omega \right)$ . Assume that

$v\in L^{\infty }\left(\left[0,T\right];H^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T\right];W^{2,q}\left(\Omega \right)\right)$ ,

$v_t\in L^{\infty }\left(\left[0,T\right];L^2\right)\cap L^2\left(\left[0,T\right];H^1\right)$

with the boundary conditions:

${v\cdot n|}_{\partial \Omega }=0$ , ${{\left(Sv\cdot n\right)}_{\tau }|}_{\partial \Omega }={\left(Bv\right)}_{\tau }$

Then there is a unique strong solution $\left(u,\theta \right)$ meeting (1.1) - (1.3) such tat

$\left(u,\theta \right)\in C\left(\left[0,T\right];H^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T\right];W^{2,q}\left(\Omega \right)\right)$ ;

$\left(u_t,{\theta }_t\right)\in C\left(\left[0,T\right];L^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T\right];H^1\left(\Omega \right)\right)$ (4.2)

Proof: It follows from Theorem 4 in chapter 5.9 [22] , then we obtain $v\in C\left(\left[0,T\right];H^1\right)$ .

Next, Gerlakin approach is applied to prove the local existence of the solution of Equation (4.1).

Assume that ${\left\{{\varphi }_l\right\}}_1^m$ and ${\left\{{\omega }_l\right\}}_1^m$ respectively representing the eigenvectors of the operator $-\Delta$ and the operator $L$ are smooth functions.

$V^m=span\left\{{\varphi }_1,{\varphi }_2,\cdots ,{\varphi }_m\right\}$ , $W^m=span\left\{{\omega }_1,{\omega }_2,\cdots ,{\omega }_m\right\}$

for a positive constant $m$ fixed, let

${\theta }_m\left(t\right)={\displaystyle \underset{l=1}{\overset{m}{\sum }}{a}_m^l\left(t\right){\varphi }_l,}{u}_m\left(t\right)={\displaystyle \underset{l=1}{\overset{m}{\sum }}{b}_m^l\left(t\right){\omega }_l}$ (4.3)

We hope that the coefficients $a_m^l\left(t\right),b_m^l\left(t\right)$ satisfy:

$\{\begin{array}{l}{\displaystyle {\int }_{\Omega }{u}_{tm}{\omega }_m+v_m\nabla u_m{u}_{tm}{\omega }_m+\nabla p_m{u}_{tm}{\omega }_m\text{d}x}={\displaystyle {\int }_{\Omega }\gamma \Delta u_m{u}_{tm}{\omega }_m+{\theta }_m{f}_m{u}_{tm}{\omega }_m\text{d}x}\hfill \\ {\displaystyle {\int }_{\Omega }{\theta }_{tm}{\varphi }_m+u_m\nabla {\theta }_m{\varphi }_m\text{d}x}={\displaystyle {\int }_{\Omega }\epsilon \Delta {\theta }_m{\varphi }_m\text{d}x}\hfill \end{array}$ (4.4)

Thus we seek functions $\left(u_m,{\theta }_m\right)$ that satisfy the “projection” (4.4) of problem (4.1) onto the finite dimensional subspace spanned by $\left({\left\{{\omega }_l\right\}}_1^m,{\left\{{\varphi }_l\right\}}_1^m\right)$ .

It follows from Theorem 1 in chapter 7 [22] , then we obtain that for $\forall m$ , there exists unique $\left(u_m,{\theta }_m\right)$ satisfying Equation (4.4).

5. The Proof of Theme 1.2

Similar to the prior estimates in part 3, we have

$\underset{0\le t\le T_{\ast }}{\sup }\left({\Vert u_m\left(t\right)\Vert }_{H^2}+{\Vert u_t{}_m\left(t\right)\Vert }_{L^2}\right)+{\displaystyle {\int }_0^{T_{\ast }}\left({\Vert \nabla u_{tm}\Vert }_{L^2}^2+{\Vert {\nabla }^2{u}_m\Vert }_{L^q}^2\right)\text{d}t}+1\le C_1$

where $C$ has no connection with $m$ , then we have

$\left(u_m,{\theta }_m\right)\to \left(u,\theta \right)\text{in}{L}^{\infty }\left(0,T;H^1\right)$

$\left(u_m,{\theta }_m\right)\stackrel{\text{weakly}}{\to }\left(u,\theta \right)\text{in}{L}^2\left(0,T;W^{2,q}\right)$

$\left(u_m,{\theta }_m\right)\stackrel{{\text{weakly}}^{\text{*}}}{\to }\left(u,\theta \right)\text{in}{L}^{\infty }\left(0,T;H^2\right)$

$\left({\partial }_t{u}_m,{\partial }_t{\theta }_m\right)\stackrel{\text{weak}}{\to }\left(u_t,{\theta }_t\right)\text{in}{L}^2\left(0,T;H^1\right)$

It follows from Theorem 3 in chapter 7 [22] , let $m\to \infty$ then we obtain that $\left(u,\theta \right)$ is the solution of (4.1).

The proof of Lemma 4.1 is completed.

Next, the iteration method is used to prove the local existence of the solution of Boussinesq equations.

Construct approximate solutions of Boussinesq equations that meet the initial and boundary problems (1.2) - (1.3).

1) define $u^0=0$ ,

2) assume that $k\ge 1$ , define $v=u^{k-1}$ ,

$\{\begin{array}{l}{u}_t^k+\left(u^{k-1}\cdot \nabla \right)u^k+\nabla p^k=\gamma \Delta u^k+{\theta }^k{f}^k,\hfill \\ {\theta }_t^k+\left(u^k\cdot \nabla \right){\theta }^k=\epsilon \Delta {\theta }^k,\hfill \\ div{u}^k=0,\hfill \end{array}\left(x,t\right)\in \Omega \times \left(0,t\right)$ (5.1)

Initial conditions:

${\left(u^k,{\theta }^k\right)|}_{t=0}=\left(u_0^k,{\theta }_0^k\right)$ (5.2)

Boundary conditions:

${u^k\cdot n|}_{\partial \Omega }=0$ , ${\frac{\partial {\theta }^k}{\partial n}|}_{\partial \Omega }=\alpha {\theta }^k$ , ${{\left(S{u}^k\cdot n\right)}_{\tau }|}_{\partial \Omega }={\left(B{u}^k\right)}_{\tau }$ (5.3)

According to Lemma 4.1, we can know that the problems (5.1) - (5.3) exist the local solutions $\left(u^k,{\theta }^k\right)$ . Furthermore, according to the prior estimates, we get

$\underset{0\le t\le T_{\ast }}{\sup }\left({\Vert u^k\Vert }_{H^2}+{\Vert u_t^k\Vert }_{L^2}\right)+{\displaystyle {\int }_0^{T_{\ast }}\left({\Vert \nabla u_t^k\Vert }_{L^2}^2+{\Vert {\nabla }^2{u}^k\Vert }_{L^q}^2\right)\text{d}t}+1\le C_1$ (5.4)

where $C$ has no connection with $k$ .

According to Aubin-Lions lemma, one has

$\left(u^k,{\theta }^k\right)\to \left(u,\theta \right)\text{in}{L}^{\infty }\left(0,T_{*};H^1\right)$

$\left(u,\theta \right)\in L^{\infty }\left(\left[0,T_{\ast }\right];H^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T_{\ast }\right];W^{2,q}\left(\Omega \right)\right)$

$\left(u_t,{\theta }_t\right)\in L^{\infty }\left(\left[0,T_{\ast }\right];L^2\left(\Omega \right)\right)\cap L^2\left(\left[0,T_{\ast }\right];H^1\left(\Omega \right)\right)$

Last, we show the continuity of $u$ and $\theta$ over time.

$\begin{array}{l}{\Vert {\left(u_t\right)}_t\Vert }_{W^{*}}=\underset{0\le t\le T_{*}}{\sup }\left|\langle {\left(u_t\right)}_t,\psi \rangle \right|\\ \le \underset{0\le t\le T_{*}}{\sup }{\displaystyle {\int }_{\Omega }\left|\gamma \Delta u_t+u_t\cdot \nabla u+u\cdot \nabla u_t+\nabla p_t+{\theta }_{t}f+\theta f_t\right|\left|\psi \right|\text{d}x}\\ \le \underset{0\le t\le T_{*}}{\sup }C(\left({\Vert u_t\Vert }_{L^2}{\Vert \nabla u\Vert }_{L^3}+{\Vert \nabla u_t\Vert }_{L^2}{\Vert u\Vert }_{L^3}+{\Vert {\theta }_t\Vert }_{L^2}{\Vert f\Vert }_{L^3}+{\Vert \theta \Vert }_{L^3}{\Vert f_t\Vert }_{L^2}\right){\Vert \psi \Vert }_{H^1}\\ +{\Vert \nabla p_t\Vert }_{L^2}{\Vert \psi \Vert }_{L^2}+{\Vert \nabla u_t\Vert }_{L^2}{\Vert \nabla \psi \Vert }_{L^2})+\underset{0\le t\le T_{*}}{\sup }{\displaystyle {\int }_{\partial \Omega }\nabla u_t\psi \text{d}s}\end{array}$

where $W^{*}$ denotes the dual space of $W$ . It follows from Theorem 3 in chapter 5.9 [22] , then we obtain $u_t\in C\left(0,T_{*};L^2\right)$ . Then according to the elliptic regularity, we have ${\nabla }^{2}u\in L^{\infty }\left(0,T_{*};L^2\right)$ .

$\begin{array}{c}{\Vert {\left({\theta }_t\right)}_t\Vert }_{W^{*}}=\underset{{\Vert \phi \Vert }_{H^1}=1}{\sup }\left|\langle {\left({\theta }_t\right)}_t,\phi \rangle \right|\\ \le \underset{{\Vert \phi \Vert }_{H^1}=1}{\sup }{\displaystyle {\int }_{\Omega }\left|\epsilon \Delta {\theta }_t+u_t\cdot \nabla \theta +u\cdot \nabla {\theta }_t\right|\left|\phi \right|\text{d}x}\\ \le \underset{{\Vert \phi \Vert }_{H^1}=1}{\sup }C(\left({\Vert u_t\Vert }_{L^2}{\Vert \nabla \theta \Vert }_{L^3}+{\Vert \nabla {\theta }_t\Vert }_{L^2}{\Vert u\Vert }_{L^3}\right){\Vert \phi \Vert }_{H^1}\\ +{\Vert \nabla {\theta }_t\Vert }_{L^2}{\Vert \nabla \phi \Vert }_{L^2})+\underset{{\Vert \phi \Vert }_{H^1}=1}{\sup }{\displaystyle {\int }_{\partial \Omega }\nabla {\theta }_t\phi \text{d}s}\end{array}$

where ${\left(H^1\left(\Omega \right)\right)}^{*}$ denotes the dual space of $H^1\left(\Omega \right)$ . This can show ${\left({\theta }_t\right)}_t\in L^2\left(0,T_{*};{\left(H^1\left(\Omega \right)\right)}^{*}\right)$ . So we have ${\theta }_t\in L^2\left(0,T_{*};H^1\left(\Omega \right)\right)$ . Similar to the proof of the continuity of $u_t$ , it is easy to know ${\theta }_t\in C\left(0,T_{*};L^2\right)$ . Then according to the elliptic regularity, we have ${\nabla }^2\theta \in C\left(0,T_{*};L^2\right)$ .

Above all, we complete the proof of the Theorem 1.2.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Ye, Z., Xu, X.J. and Xue, L.T. (2017) On the Global Regularity of the 2D Boussinesq Equations with Fractional Dissipation. Mathematische Nachrichten, 290, 1420-1439. https://doi.org/10.1002/mana.201500413
[2]	Ye, Z. and Xu, X.J. (2016) Global Well-Posedness of the 2D Boussinesq Equations with Fractional Laplacian Dissipation. Journal of Differential Equations, 260, 6716-6744. https://doi.org/10.1016/j.jde.2016.01.014
[3]	Fan, J.S., Jia, X.J., Nakamura, G. and Zhou, Y. (2015) On Well-Posedness and Blowup Criteria for the Magnetohydrodynamics with the Hall and Ion-Slip Effects. Z. Angew. Math. Phys., 66, 1695-1706.
[4]	Hmidi, T., Keraani, S. and Rousset, F. (2011) Global Well-Posedness for Euler-Boussinesq System with Critical Dissipation, Comm. Partial Differential Equations, 36, 420-445. https://doi.org/10.1080/03605302.2010.518657
[5]	Liu, Y. The Attenuation in L2 of the Weak Existence of the Boussinesq Equations. Master’s Thesis, Capital Normal University, Beijing.
[6]	Li, W. (2008) The Exact Solutions of the Boussinesq Equations. Journal of Bohai University, 29.
[7]	Li, X.J. The Global Stability in L2 of the Boussinesq Equations in Three Dimensional Regional and Unbounded Regional. Master’s Thesis, Capital Normal University, Beijing.
[8]	Tan, Z. and Wang, Y.J. (2015) Global Well-Posedness of an Initial-Boundary Value Problem for Viscous Non-Resistive MHD Systems. Mathematics, arXiv: 1509.08349.
[9]	Ohno, M. and Hirota, T. (1998) On the Initial-Boundary Value Problem for the Linearized MHD Equations. Sbrikaisekikenkybsho K^kybroku, 1045, 26-44.
[10]	Matsumura, and Nismda, T. (1983) Initial-Boundary Value Problems of the Equations of Motion of General Fluids. Comm. Math. Phys., 89, 445-446.
[11]	LIONS (1998) Mathematical Topics in Fluid Mechanics. Vol. 2, Oxford Lecture Ser. Math. 10, Oxford University Press, New York.
[12]	HOFF (2011) Existence of Solutions to a Model for Sparse, One-Dimensional Fluids. Journal of Differential Equations, 250, 1083-1113.
[13]	HOFF and D.SERRE (1991) The Failure of Continuous Dependence on Initial Data for the Navier-Stokes Equations of Compressible Flow. SIAM Journal on Applied Mathematics, 51, 887-898. https://doi.org/10.1137/0151043
[14]	HOFF (2005) Compressible Flow in a Half-Space with Navier Boundary Conditions. Journal of Mathematical Fluid Mechanics, 7, 315-338. https://doi.org/10.1007/s00021-004-0123-9
[15]	HOFF (1997) Global Solutions of the Navier-Stokes Equations for Multidimensional, Compressible Flow with Discontinuous Initial Data. Journal of Differential Equations, 120, 1742-1760.
[16]	Fan, J., Jiang, S., Nakamura, G. and Zhou, Y. (2011) Logarithmically Improved Regularity Criteria for the Navier-Stokes and MHD Equations. Journal of Mathematical Fluid Mechanics, 13, 557.
[17]	HOFF and M.SANTOS (2008) Lagrangean Structure and Propagation of Singularities in Multidimensional Compressible Flow. Archive for Rational Mechanics and Analysis, 188, 509-543. https://doi.org/10.1007/s00205-007-0099-8
[18]	HOFF (2006) Uniqueness of Weak Solutions of the Navier-Stokes Equations of Multidimensional Compressible Flow. SIAM Journal on Mathematical Analysis, 37, 1742-1760. https://doi.org/10.1137/040618059
[19]	Hoff, D. (2012) Local Solutions of a Compressible Flow Problem with Navier Boundary Conditions in General Three Dimensional Domains Department of Mathematics. Indiana University, Bloomington, 633-650.
[20]	Hmidi, T., Keraani, S. and Rousset, F. (2010) Global Well-Posedness for a Boussinesq-Navier-Stokes System with Critical Dissipation. Journal of Differential Equations, 249, 2147-2174.
[21]	Hoff, D. (2011) Local Solutions of a Compressible Flow Problem with Navier Boundary Conditions in General Three-Dimensional Domains. SIAM Journal on Mathematical Analysis, 44, 1872-1886.
[22]	Evans (2010) Partial Differential Equations, Grad. Texts in Math. 19, American Mathematical Society, Providence.