Jincheng Shi
Huashang College Guangdong University of Finance & Economics, Guangzhou, China.
DOI: 10.4236/jamp.2020.87099   PDF    HTML   XML   523 Downloads   1,389 Views   Citations

Abstract

In this paper, we establish the structural stability for the linear differential equations of thermo-diffusion in a semi-infinite pipe flow. Using the technology of a second-order differential inequality, we prove the continuous dependence on the density ρ and the coefficient of thermal conductivity K. These results show that small changes for these coefficients can’t cause tremendous changes for the solutions.

Keywords

Differential Equations of Thermo-Diffusion, Structural Stability, Continuous Dependence

Share and Cite:

1. Introduction

The question of continuous dependence of solutions of problems in partial differential equations on coefficients in the equations has been extensively studied in recent years for a variety of problems. This is sometimes referred to as the question of structural stability and numerous references may be found, for instance, in the book of Ames and Straughan [1] and the monograph of Straughan [2]. For more papers one can see [3] - [8]. In structural stability the emphasis is on continuous dependence (convergence result) on changes in the model itself rather than on the initial data. This means changes in coefficients in the partial differential equations and changes in the equations and may be reflected physically by changes in constitutive parameters. What’s more, the inevitable error that arises in both numerical computation and the physical measurement of data can exist. It is relevant to know the magnitude of the effect of such errors on the solution.

In the 1970s, W. Nowacki in his papers [9] [10] gave the differential equations of thermodiffusion in one dimensional space and many papers in the literatures have studied this system. For example, [10] [11] [12] investigated the initial-boundary value problem for the linear system of thermodiffusion using different arguments. [13] proved the existence, uniqueness and regularity of the solution to the initial-boundary value problems for the linear system of thermodiffusion in a solid body. L^p-L^q time decay estimates for the solution of the associated linear Cauchy problem were obtained by [14]. However, in this paper, we considered the differential equations of thermodiffusion in three dimensions and we study not only the continuous dependence on the coefficients of the equations, but also the spatial decay estimates for the solution of the system. In fact, much has been written on the subject of spatial decay bounds for various systems of differential equations, e.g., for a review of such works on Saint-Venant’s principle, one can refer to [15] - [23] and the papers cited therein. Recently, there are some new results about structural stability, one could see [24] - [28].

We shall assume that a transient flow occupies the interior of a semi-infinite cylindrical pipe $R$ with boundary $\partial R$. The pipe has arbitrary cross section denoted by $D$ and the boundary $\partial D$ and the generators of the pipe are parallel to the $x_3$ axis. We introduce the notations:

$R_z=\left\{\left(x_1,x_2,x_3\right)|\left(x_1,x_2\right)\in D,x_3>z\ge 0\right\},$

$D_z=\left\{\left(x_1,x_2,x_3\right)|\left(x_1,x_2\right)\in D,x_3=z\ge 0\right\},$

where z is a running variable along the $x_3$ axis. Clearly, $R_0=R$ and $D_0=D$. Let $u_i$, $T$, and $C$ denote the displacement, temperature, and chemical potential as independent fields, respectively. These fields depend on the space variable $\left(x_1\mathrm{,}{x}_2\mathrm{,}{x}_3\right)$ and the time variable t and satisfy the following system of equations:

$\rho {\ddot{u}}_i-\nu \Delta u_i-\left(\lambda +\nu \right)u_{j\mathrm{,}ji}+{\gamma }_1{T}_{\mathrm{,}i}+{\gamma }_2{C}_{\mathrm{,}i}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}$ (1.1)

$c\stackrel{\dot{}}{T}-K\Delta T+{\gamma }_1{\stackrel{\dot{}}{u}}_{i\mathrm{,}i}+d\stackrel{\dot{}}{C}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}$ (1.2)

$n\stackrel{\dot{}}{C}-M\Delta C+{\gamma }_2{\stackrel{\dot{}}{u}}_{i\mathrm{,}i}+d\stackrel{\dot{}}{T}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}$ (1.3)

with the initial-boundary conditions

$u_i=0,T=C=0\text{on}\partial D\times \left\{t\ge 0\right\},$ (1.4)

$u_i={\stackrel{\dot{}}{u}}_i=0,T=C=0\text{in}R\times \left\{t=0\right\}.$ (1.5)

$u_i=f_i\left(x_1,x_2,t\right),T=F\left(x_1,x_2,t\right),C=G\left(x_1,x_2,t\right)\text{on}{D}_0\times \left\{t\ge 0\right\},$ (1.6)

$u_i,u_{i,j},T,T_{,i},C,C_{,i},p=o\left(x_3^{-1}\right)\text{uniformly}\text{in}{x}_1,x_2,t\text{as}{x}_3\to \infty .$ (1.7)

In Equations (1.1)-(1.3), $\Delta$ is the Laplacian operator; $\rho$ represents the density; ${\gamma }_1$ and ${\gamma }_2$ are the coefficients of thermal and diffusion dilatation; $\lambda$ and $\nu$ are the material coefficients; K is the coefficient of thermal conductivity; M is the coefficient of diffusion. $n\mathrm{,}c\mathrm{,}d$ are the coefficients of thermodiffusion. All the above constants are positive and satisfy

$\begin{array}{l}cn-d^2>0,\hfill \end{array}$ (1.8)

which implies that (1.1)-(1.3) is a hyperbolic-parabolic system of partial differential equations. In the following several sections, we may use the below inequality. Let $D$ be a plane domain $D$ with the boundary $\partial D$. If $w=0$ on $\partial D$, then

${\displaystyle {\int }_D}{w}_{\mathrm{,}\alpha }{w}_{\mathrm{,}\alpha }\text{d}A\ge {\lambda }_1{\displaystyle {\int }_D}{w}^2\text{d}x\mathrm{,}$ (1.9)

where ${\lambda }_1$ is the smallest eigenvalue of the problem

$\Delta \varphi +\lambda \varphi =0\text{in}D\mathrm{,}$

$\varphi =0\text{on}\partial D\mathrm{.}$

This inequality has been well studied (see [29] [30]). Throughout this paper, the usual summation convention is employed with repeated Latin subscripts summed from 1 to 3 and repeat Greek subscript summed from 1 to 2. The

comma is used to indicate partial differentiation, i.e. $u_{i,j}=\frac{\partial u_i}{\partial x_j}$, ${\phi }_{\alpha ,\alpha }={\displaystyle {\sum }_{\alpha =1}^2}\frac{\partial {\phi }_{\alpha }}{\partial x_{\alpha }}$ and ${\stackrel{\dot{}}{u}}_i$ denotes $\frac{\partial u_i}{\partial t}$.

The paper is structured as follows: In Section 2, we derive the continuous dependence on $\rho$. Section 3 is devoted to seeking the continuous dependence on K.

2. Continuous Dependence on the Parameter ρ

Lemma 1. The energy $E\left(z,t\right)$ defined in [31]

$\begin{array}{c}E\left(z,t\right)=\frac{1}{2}\rho {\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_i{\stackrel{\dot{}}{u}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{u}_{i,j}{u}_{i,j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{u}_{i,i}^2\text{d}x\\ +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{,j}{T}_{,j}\text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{C}_{,j}{C}_{,j}\text{d}x\text{d}\eta \\ +{\displaystyle {\int }_{R_z}}\left[\frac{c}{2}{T}^2+dCT+\frac{n}{2}{C}^2\right]\text{d}x,\end{array}$ (2.1)

satisfies the following estimates

$E\left(z,t\right)\le E\left(0,t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}.$ (2.2)

where

$m_1\left(t\right)=\frac{\sqrt{\nu }t}{\sqrt{\rho }}+\frac{\sqrt{\lambda +\nu }t}{\sqrt{\rho }}+\frac{{\gamma }_1\sqrt{t}}{2\sqrt{\rho {\lambda }_{1}K}}+\frac{{\gamma }_2\sqrt{t}}{2\sqrt{\rho {\lambda }_{1}M}}+\frac{1}{2\sqrt{{\lambda }_1}}.$ (2.3)

Proof. These results are the main results of paper [31].

Theorem 1. The energy expression $\tilde{\phi }\left(z\mathrm{,}t\right)$ satisfies the following estimates:

$\tilde{\phi }\left(z\mathrm{,}t\right)\le \frac{2{\pi }^2}{\rho {\rho }^{*}}\tilde{E}\left(0,t\right){\text{e}}^{-\frac{1}{2m_1\left(t\right)}z}+\frac{{\pi }^2{t}^2}{{\rho }^{*}\rho }\tilde{E}\left(0,t\right)\left({\text{e}}^{-\frac{1}{2m_1\left(t\right)}z}-{\text{e}}^{-\frac{3}{2m_1\left(t\right)}z}\right)\mathrm{.}$ (2.4)

Proof. To investigate continuous dependence on $\rho$, we have to seek a bound for $\frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x$. To do this, we first differentiate (1.1), and then multiply with ${\ddot{u}}_i$, integrate over the region $R_z\times \left[\mathrm{0,}t\right]$ to obtain

$0={\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left[\rho {\stackrel{⃛}{u}}_i-\nu \Delta {\stackrel{\dot{}}{u}}_i-\left(\lambda +\nu \right){\stackrel{\dot{}}{u}}_{j,ji}+{\gamma }_1{\stackrel{\dot{}}{T}}_{,i}+{\gamma }_2{\stackrel{\dot{}}{C}}_{,i}\right]{\ddot{u}}_i\text{d}x\text{d}\eta ,$ (2.5)

which follows that

$\begin{array}{c}0=\frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x+\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{u}}_{i,3}{\ddot{u}}_i\text{d}A\text{d}\eta +\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_{i,j}{\stackrel{\dot{}}{u}}_{i,j}\text{d}x\\ +\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{u}}_{j,j}{\ddot{u}}_3\text{d}A\text{d}\eta +\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_{j,j}^2\text{d}x\\ +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}_{,i}{\ddot{u}}_i\text{d}x\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{C}}_{,i}{\ddot{u}}_i\text{d}x\text{d}\eta ,\end{array}$ (2.6)

where we have supposed that ${\ddot{u}}_i$ vanish at $t=0$. Similarly, we have

$\begin{array}{c}0=c{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\ddot{T}\stackrel{\dot{}}{T}\text{d}x\text{d}\eta -K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\Delta \ddot{T}\stackrel{\dot{}}{T}\text{d}x\text{d}\eta +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_{i,i}\stackrel{\dot{}}{T}\text{d}x\text{d}\eta +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\ddot{C}\stackrel{\dot{}}{T}\text{d}x\text{d}\eta \\ =\frac{c}{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}^2\text{d}x+K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{T}}_{,3}\stackrel{\dot{}}{T}\text{d}A\text{d}\eta +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}_{,i}{\stackrel{\dot{}}{T}}_{,i}\text{d}x\text{d}\eta \\ -{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\ddot{u}}_3\stackrel{\dot{}}{T}\text{d}A\text{d}\eta -{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}_{,i}{\ddot{u}}_i\text{d}x\text{d}\eta +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\ddot{C}\stackrel{\dot{}}{T}\text{d}x\text{d}\eta ,\end{array}$ (2.7)

and

$\begin{array}{c}0=n{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\ddot{C}\stackrel{\dot{}}{C}\text{d}x\text{d}\eta -M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\Delta \stackrel{\dot{}}{C}\stackrel{\dot{}}{C}\text{d}x\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_{i,i}\stackrel{\dot{}}{C}\text{d}x\text{d}\eta +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\ddot{T}\stackrel{\dot{}}{C}\text{d}x\text{d}\eta \\ =\frac{n}{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{C}}^2\text{d}x+M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{C}}_{,3}\stackrel{\dot{}}{C}\text{d}A\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{C}}_{,i}{\stackrel{\dot{}}{C}}_{,i}\text{d}x\text{d}\eta \\ -{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\ddot{u}}_3\stackrel{\dot{}}{C}dAd\eta -{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{C}}_{,i}{\ddot{u}}_i\text{d}x\text{d}\eta +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\ddot{T}\stackrel{\dot{}}{C}\text{d}x\text{d}\eta ,\end{array}$ (2.8)

where we also have assumed that $\stackrel{\dot{}}{T}=\stackrel{\dot{}}{C}=0$ at $t=0$. Combining (6)-(8), we set

$\begin{array}{c}\tilde{E}\left(z\mathrm{,}t\right)=\frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_{i\mathrm{,}j}{\stackrel{\dot{}}{u}}_{i\mathrm{,}j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_{j\mathrm{,}j}^2\text{d}x+K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}_{\mathrm{,}i}{\stackrel{\dot{}}{T}}_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{C}}_{\mathrm{,}i}{\stackrel{\dot{}}{C}}_{\mathrm{,}i}\text{d}x\text{d}\eta +\frac{c}{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}^2\text{d}x+\frac{1}{2}{\displaystyle {\int }_{R_z}}\left[c{\stackrel{\dot{}}{T}}^2+2d\stackrel{\dot{}}{T}\stackrel{\dot{}}{C}+n{\stackrel{\dot{}}{C}}^2\right]\text{d}x\\ =-\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{u}}_{i\mathrm{,3}}{\ddot{u}}_i\text{d}A\text{d}\eta -\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{u}}_{j\mathrm{,}j}{\ddot{u}}_3\text{d}A\text{d}\eta -K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{T}}_{\mathrm{,3}}\stackrel{\dot{}}{T}\text{d}A\text{d}\eta \\ +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\ddot{u}}_3\stackrel{\dot{}}{T}\text{d}A\text{d}\eta -M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{C}}_{\mathrm{,3}}\stackrel{\dot{}}{C}\text{d}A\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\ddot{u}}_3\stackrel{\dot{}}{C}\text{d}A\text{d}\eta \mathrm{.}\end{array}$ (2.9)

Following the method used in [31], we can get

$\begin{array}{c}\tilde{E}\left(z\mathrm{,}t\right)=\frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_{i\mathrm{,}j}{\stackrel{\dot{}}{u}}_{i\mathrm{,}j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{u}}_{j\mathrm{,}j}^2\text{d}x+K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}_{\mathrm{,}i}{\stackrel{\dot{}}{T}}_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{C}}_{\mathrm{,}i}{\stackrel{\dot{}}{C}}_{\mathrm{,}i}\text{d}x\text{d}\eta +\frac{c}{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{T}}^2\text{d}x+\frac{1}{2}{\displaystyle {\int }_{R_z}}\left[c{\stackrel{\dot{}}{T}}^2+2d\stackrel{\dot{}}{T}\stackrel{\dot{}}{C}+n{\stackrel{\dot{}}{C}}^2\right]\text{d}x\\ \le \tilde{E}\left(\mathrm{0,}t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\mathrm{.}\end{array}$

Also, we employ the argument used in [31] to get that $\tilde{E}\left(\mathrm{0,}t\right)$ may be bounded by known data. Since $cn>d^2$, we note again ${\displaystyle {\int }_{R_z}}\left[c{\stackrel{\dot{}}{T}}^2+2d\stackrel{\dot{}}{T}\stackrel{\dot{}}{C}+n{\stackrel{\dot{}}{C}}^2\right]\text{d}x>0$. So, we have

${\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x\le \frac{2}{\rho }\tilde{E}\left(\mathrm{0,}t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\mathrm{.}$ (2.10)

Now, we study the continuous dependence on the parameter $\rho$. Let $\left(u_i\mathrm{,}T\mathrm{,}C\right)$ and $\left(u_i^{\mathrm{<mo>\; *\; </mo>}}\mathrm{,}{T}^{\mathrm{<mo>\; *\; </mo>}}\mathrm{,}{C}^{\mathrm{<mo>\; *\; </mo>}}\right)$ be the solutions to (1.1)-(1.3) with same initial-boundary conditions, but for different parameters $\rho$ and ${\rho }^{\mathrm{<mo>\; *\; </mo>}}$, respectively. Define the difference variables as

$w_i=u_i-u_i^{*},\theta =T-T^{*},\Sigma =C-C^{*},\pi =\rho -{\rho }^{*}.$ (2.11)

Then, $\left(w_i\mathrm{,}\theta \mathrm{,}\Sigma \right)$ satisfy

$\pi {\ddot{u}}_i+{\rho }^{\mathrm{<mo>\; *\; </mo>}}{\ddot{w}}_i-\nu \Delta w_i-\left(\lambda +\nu \right)w_{j\mathrm{,}ji}+{\gamma }_1{\theta }_{\mathrm{,}i}+{\gamma }_2{\Sigma }_{\mathrm{,}i}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}$ (2.12)

$c\stackrel{\dot{}}{\theta }-K\Delta \theta +{\gamma }_1{\stackrel{\dot{}}{w}}_{i\mathrm{,}i}+d\stackrel{\dot{}}{\Sigma }=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}$ (2.13)

$n\stackrel{\dot{}}{\Sigma }-M\Delta \Sigma +{\gamma }_2{\stackrel{\dot{}}{w}}_{i\mathrm{,}i}+d\stackrel{\dot{}}{\theta }=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}$ (2.14)

with the initial-boundary conditions

$w_i=0,\theta =\Sigma =0\text{on}\partial D\times \left\{t\ge 0\right\},$ (2.15)

$w_i={\stackrel{\dot{}}{w}}_i=0,\theta =\Sigma =0\text{in}R\times \left\{t=0\right\}.$ (2.16)

$w_i=\theta =\Sigma =0\text{on}{D}_0\times \left\{t\ge 0\right\},$ (2.17)

Multiplying (2.12) with ${\stackrel{\dot{}}{w}}_i$ and integrating by parts, we have

$\begin{array}{c}0={\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left[\pi {\ddot{u}}_i+{\rho }^{*}{\ddot{w}}_i-\nu \Delta w_i-\left(\lambda +\nu \right)w_{j,ji}+{\gamma }_1{\theta }_{,i}+{\gamma }_2{\Sigma }_{,i}\right]{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ =\frac{1}{2}{\displaystyle {\int }_{R_z}}{\rho }^{*}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_{i}dx+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i,j}{w}_{i,j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i,i}^2\text{d}x\\ +\pi {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ +\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{i,3}{\stackrel{\dot{}}{w}}_i\text{d}A\text{d}\eta +\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{j,j}{\stackrel{\dot{}}{w}}_3\text{d}A\text{d}\eta .\end{array}$ (2.18)

Similarly, we have

$\begin{array}{c}0={\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left[c\stackrel{\dot{}}{\theta }-K\Delta \theta +{\gamma }_1{w}_{i,i}+d\stackrel{\dot{}}{\Sigma }\right]\theta \text{d}x\text{d}\eta \\ =\frac{c}{2}{\displaystyle {\int }_{R_z}}{\theta }^2\text{d}x\text{d}\eta +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta -{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\stackrel{\dot{}}{\Sigma }\theta \text{d}x\text{d}\eta +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\theta }_{,3}\theta \text{d}A\text{d}\eta -{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\theta \text{d}A\text{d}\eta ,\end{array}$ (2.19)

and

$\begin{array}{c}0={\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left[n\stackrel{\dot{}}{\Sigma }-M\Delta \Sigma +{\gamma }_2{w}_{i,i}+d\stackrel{\dot{}}{\theta }\right]\Sigma \text{d}x\\ =\frac{n}{2}{\displaystyle {\int }_{R_z}}{\Sigma }^2\text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{,i}{\Sigma }_{,i}\text{d}x\text{d}\eta -{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\stackrel{\dot{}}{\theta }\Sigma \text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\Sigma }_{,3}\Sigma \text{d}A\text{d}\eta -{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\Sigma \text{d}A\text{d}\eta .\end{array}$ (2.20)

We define a new function

$\begin{array}{c}\phi \left(z\mathrm{,}t\right)=\frac{1}{2}{\displaystyle {\int }_{R_z}}{\rho }^{*}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i,j}{w}_{i,j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i,i}^2\text{d}x+K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta \\ +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{,i}{\Sigma }_{,i}\text{d}x\text{d}\eta +\frac{1}{2}{\displaystyle {\int }_{R_z}}\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x\text{d}\eta \\ =\pi {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta -\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{i,3}{\stackrel{\dot{}}{w}}_i\text{d}A\text{d}\eta -\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{j,j}{\stackrel{\dot{}}{w}}_3\text{d}A\text{d}\eta \\ -K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\theta }_{,3}\theta \text{d}A\text{d}\eta +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\theta \text{d}A\text{d}\eta \\ -M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\Sigma }_{,3}\Sigma \text{d}A\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\Sigma \text{d}A\text{d}\eta .\end{array}$ (2.21)

Thus, we have

$\begin{array}{c}-\frac{\partial \phi \left(z\mathrm{,}t\right)}{\partial z}=\frac{1}{2}{\displaystyle {\int }_{D_z}}{\rho }^{\mathrm{<mo>\; *\; </mo>}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}A+\frac{\nu }{2}{\displaystyle {\int }_{D_z}}{w}_{i\mathrm{,}j}{w}_{i\mathrm{,}j}\text{d}A+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{D_z}}{w}_{i\mathrm{,}i}^2\text{d}A\\ +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\theta }_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}A\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\Sigma }_{\mathrm{,}i}{\Sigma }_{\mathrm{,}i}\text{d}A\text{d}\eta \\ +\frac{1}{2}{\displaystyle {\int }_{D_z}}\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}A\text{d}\eta \mathrm{.}\end{array}$ (2.22)

Similar to [31], we have

$\phi \left(z\mathrm{,}t\right)\le \pi {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta +m_1\left(t\right)\left(-\frac{\partial \phi \left(z\mathrm{,}t\right)}{\partial z}\right)\mathrm{,}$ (2.23)

where $m_1\left(t\right)$ have been defined in (2.3), but $\rho$ in $m_1\left(t\right)$ may be replaced by ${\rho }^{\mathrm{<mo>\; *\; </mo>}}$ here. By Hölder and the AG mean inequalities, from (2.23) we have

$\begin{array}{c}\phi \left(z\mathrm{,}t\right)\le \pi {\left({\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x\text{d}\eta \right)}^{\frac{1}{2}}{\left({\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \right)}^{\frac{1}{2}}+m_1\left(t\right)\left(-\frac{\partial \phi \left(z\mathrm{,}t\right)}{\partial z}\right)\\ \le \frac{{\pi }^{2}t}{{\rho }^{\mathrm{<mo>\; *\; </mo>}}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x\text{d}\eta +\frac{{\rho }^{\mathrm{<mo>\; *\; </mo>}}}{4}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+m_1\left(t\right)\left(-\frac{\partial \phi \left(z\mathrm{,}t\right)}{\partial z}\right)\mathrm{.}\end{array}$ (2.24)

Now, we let

$\begin{array}{c}\tilde{\phi }\left(z\mathrm{,}t\right)=\frac{1}{4}{\displaystyle {\int }_{R_z}}{\rho }^{\mathrm{<mo>\; *\; </mo>}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i\mathrm{,}j}{w}_{i\mathrm{,}j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i\mathrm{,}i}^2\text{d}x\\ +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{\mathrm{,}i}{\Sigma }_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +\frac{1}{2}{\displaystyle {\int }_{R_z}}\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x\text{d}\eta \mathrm{.}\end{array}$ (2.25)

From (2.24) and (2.10), we have

$\tilde{\phi }\left(z\mathrm{,}t\right)\le \frac{{\pi }^{2}t}{{\rho }^{\mathrm{<mo>\; *\; </mo>}}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x\text{d}\eta +2m_1\left(t\right)\left(-\frac{\partial \tilde{\phi }\left(z\mathrm{,}t\right)}{\partial z}\right)\mathrm{,}$ (2.26)

which follows that

$\tilde{\phi }\left(z\mathrm{,}t\right)\le \tilde{\phi }\left(\mathrm{0,}t\right){\text{e}}^{-\frac{1}{2m_1\left(t\right)}z}+\frac{{\pi }^2{t}^2}{{\rho }^{\mathrm{<mo>\; *\; </mo>}}\rho }\tilde{E}\left(0,t\right)\left(1-{\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\right){\text{e}}^{-\frac{1}{2m_1\left(t\right)}z}\mathrm{.}$ (2.27)

In order to make inequality (2.26) explicit, we need bound for $\tilde{\phi }\left(\mathrm{0,}t\right)$. So, we write (2.21) at $z=0$ and use the initial-boundary conditions to obtain

$\begin{array}{c}\phi \left(\mathrm{0,}t\right)=\frac{1}{2}{\displaystyle {\int }_{R_0}}{\rho }^{*}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_0}}{w}_{i,j}{w}_{i,j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_0}}{w}_{i,i}^2\text{d}x\\ +K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_0}}{\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_0}}{\Sigma }_{,i}{\Sigma }_{,i}\text{d}x\text{d}\eta \\ +\frac{1}{2}{\displaystyle {\int }_{R_0}}\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x\text{d}\eta \\ =\pi {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_0}}{\ddot{u}}_i{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ \le \frac{{\pi }^{2}t}{{\rho }^{*}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_0}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x\text{d}\eta +\frac{{\rho }^{*}}{4}{\displaystyle {\int }_{R_0}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x,\end{array}$ (2.28)

which results in

$\tilde{\phi }\left(\mathrm{0,}t\right)\le \frac{{\pi }^{2}t}{{\rho }^{\mathrm{<mo>\; *\; </mo>}}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_0}}{\ddot{u}}_i{\ddot{u}}_i\text{d}x\text{d}\eta \le \frac{2{\pi }^2}{\rho {\rho }^{\mathrm{<mo>\; *\; </mo>}}}\tilde{E}\left(0,t\right)\mathrm{,}$ (2.29)

where we have used (2.10). Combining (2.27) and (2.29), we have

$\tilde{\phi }\left(z\mathrm{,}t\right)\le \frac{2{\pi }^2}{\rho {\rho }^{\mathrm{<mo>\; *\; </mo>}}}\tilde{E}\left(0,t\right){\text{e}}^{-\frac{1}{2m_1\left(t\right)}z}+\frac{{\pi }^2{t}^2}{{\rho }^{\mathrm{<mo>\; *\; </mo>}}\rho }\tilde{E}\left(0,t\right)\left({\text{e}}^{-\frac{1}{2m_1\left(t\right)}z}-{\text{e}}^{-\frac{3}{2m_1\left(t\right)}z}\right)\mathrm{.}$ (2.30)

Inequality (2.30) shows that the amplitude terms in (2.25) become small as $\rho \to {\rho }^{\mathrm{<mo>\; *\; </mo>}}$ and the continuous dependence on $\rho$ is obtained.

3. Continuous Dependence on the Parameter K

Theorem 2. The energy expression $\tilde{\Phi }\left(z,t_1\right)$ satisfies the following estimates:

If $h_0-k_1-k_2=0$, we have

$\tilde{\Phi }\left(z,t_1\right)\le {\beta }^2{m}_4{\text{e}}^{-\left(h_0-k_1\right)z}+{\beta }^2{m}_3{\text{e}}^{-\left(h_0-k_1\right)z}\mathrm{.}$ (3.1)

If $h_0-k_1-k_2\ne 0$, we have

$\tilde{\Phi }\left(z,t_1\right)\le {\beta }^2\left[m_4-\frac{m_3}{h_0-k_1-k_2}\right]{\text{e}}^{-\left(h_0-k_1\right)z}+\frac{{\beta }^2{m}_3}{h_0-k_1-k_2}{\text{e}}^{-k_{2}z}\mathrm{.}$ (3.2)

Proof. In this section we compare the solutions of the following two problems

$\begin{array}{l}\rho {\ddot{u}}_i-\nu \Delta u_i-\left(\lambda +\nu \right)u_{j\mathrm{,}ji}+{\gamma }_1{T}_{\mathrm{,}i}+{\gamma }_2{C}_{\mathrm{,}i}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}\\ c\stackrel{\dot{}}{T}-K\Delta T+{\gamma }_1{\stackrel{\dot{}}{u}}_{i\mathrm{,}i}+d\stackrel{\dot{}}{C}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}\\ n\stackrel{\dot{}}{C}-M\Delta C+{\gamma }_2{\stackrel{\dot{}}{u}}_{i\mathrm{,}i}+d\stackrel{\dot{}}{T}=\mathrm{0,}\text{in}R\times \left\{t\ge 0\right\}\mathrm{,}\end{array}$ (3.3)

with the initial-boundary conditions

$\begin{array}{l}{u}_i=0,T=C=0\text{on}\partial D\times \left\{t\ge 0\right\},\\ u_i={\stackrel{\dot{}}{u}}_i=0,T=C=0\text{in}R\times \left\{t=0\right\},\\ u_i=f_i\left(x_1,x_2,t\right),T=F\left(x_1,x_2,t\right),C=G\left(x_1,x_2,t\right)\text{on}{D}_0\times \left\{t\ge 0\right\},\end{array}$ (3.4)

and

$\begin{array}{l}\rho {\ddot{u}}_i^{*}-\nu \Delta u_i^{*}-\left(\lambda +\nu \right)u_{j,ji}^{*}+{\gamma }_1{T}_{,i}^{*}+{\gamma }_2{C}_{,i}^{*}=0,\text{in}R\times \left\{t\ge 0\right\},\\ c{\stackrel{\dot{}}{T}}^{*}-K^{*}\Delta T^{*}+{\gamma }_1{\stackrel{\dot{}}{u}}_{i,i}^{*}+d{\stackrel{\dot{}}{C}}^{*}=0,\text{in}R\times \left\{t\ge 0\right\},\\ n{\stackrel{\dot{}}{C}}^{*}-M\Delta C^{*}+{\gamma }_2{\stackrel{\dot{}}{u}}_{i,i}^{*}+d{\stackrel{\dot{}}{T}}^{*}=0,\text{in}R\times \left\{t\ge 0\right\},\end{array}$ (3.5)

with the same initial-boundary conditions (3.4).

Our goal in this section is to derive the continuous dependence on the parameter K. If we define

$w_i=u_i-u_i^{*},\theta =T-T^{*},\Sigma =C-C^{*},\beta =K-K^{*},$ (3.6)

then, $\left(w_i\mathrm{,}\theta \mathrm{,}\Sigma \right)$ satisfy the system

$\begin{array}{l}\rho {\ddot{w}}_i-\nu \Delta w_i-\left(\lambda +\nu \right)w_{j,ji}+{\gamma }_1{\theta }_{,i}+{\gamma }_2{\Sigma }_{,i}=0,\text{in}R\times \left\{t\ge 0\right\},\\ c\stackrel{\dot{}}{\theta }-\beta \Delta T-K^{*}\Delta \theta +{\gamma }_1{\stackrel{\dot{}}{w}}_{i,i}+d\stackrel{\dot{}}{\Sigma }=0,\text{in}R\times \left\{t\ge 0\right\},\\ n\stackrel{\dot{}}{\Sigma }-M\Delta \Sigma +{\gamma }_2{\stackrel{\dot{}}{w}}_{i,i}+d\stackrel{\dot{}}{\theta }=0,\text{in}R\times \left\{t\ge 0\right\},\end{array}$ (3.7)

with the initial-boundary conditions

$\begin{array}{l}{w}_i=0,\theta =\Sigma =0\text{on}\partial D\times \left\{t\ge 0\right\},\\ w_i={\stackrel{\dot{}}{w}}_i=0,\theta =\Sigma =0\text{in}R\times \left\{t=0\right\},\\ w_i=0,\theta =\Sigma =0\text{on}{D}_0\times \left\{t\ge 0\right\}.\end{array}$ (3.8)

We multiply (3.7)₁ with $w_i$ and integrate by parts to have

$\begin{array}{c}0={\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left[\rho {\ddot{w}}_i-\nu \Delta w_i-\left(\lambda +\nu \right)w_{j,ji}+{\gamma }_1{\theta }_{,i}+{\gamma }_2{\Sigma }_{,i}\right]{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ =\frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i,j}{w}_{i,j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i,i}^2\text{d}x+{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta \\ +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta +\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{i,3}{\stackrel{\dot{}}{w}}_i\text{d}A\text{d}\eta +\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{j,j}{\stackrel{\dot{}}{w}}_3\text{d}A\text{d}\eta .\end{array}$ (3.9)

Similarly, we have

$\begin{array}{c}0={\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left[c\stackrel{\dot{}}{\theta }-\beta \Delta T-K^{*}\Delta \theta +{\gamma }_1{w}_{i,i}+d\stackrel{\dot{}}{\Sigma }\right]\theta \text{d}x\text{d}\eta \\ =\frac{c}{2}{\displaystyle {\int }_{R_z}}{\theta }^2\text{d}x\text{d}\eta +\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{,i}{\theta }_{,i}\text{d}x\text{d}\eta +K^{*}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta \\ -{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta +d{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\stackrel{\dot{}}{\Sigma }\theta \text{d}x\text{d}\eta +\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{T}_{,3}\theta \text{d}A\text{d}\eta \\ +K^{*}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\theta }_{,3}\theta \text{d}A\text{d}\eta -{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\theta \text{d}A\text{d}\eta .\end{array}$ (3.10)

Combining (2.20), (3.9) and (3.10), we have

$\begin{array}{c}\Phi \left(z\mathrm{,}t\right)\doteq \frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i\mathrm{,}j}{w}_{i\mathrm{,}j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i\mathrm{,}i}^2\text{d}x+K^{\mathrm{<mo>\; *\; </mo>}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{\mathrm{,}i}{\Sigma }_{\mathrm{,}i}\text{d}x\text{d}\eta +\frac{1}{2}{\displaystyle {\int }_{R_z}}\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x\text{d}\eta \\ =-\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{i\mathrm{,3}}{\stackrel{\dot{}}{w}}_i\text{d}A\text{d}\eta -\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{w}_{j\mathrm{,}j}{\stackrel{\dot{}}{w}}_3\text{d}A\text{d}\eta -K^{\mathrm{<mo>\; *\; </mo>}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\theta }_{\mathrm{,3}}\theta \text{d}A\text{d}\eta \\ +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\theta \text{d}A\text{d}\eta -M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\Sigma }_{\mathrm{,3}}\Sigma \text{d}A\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{\stackrel{\dot{}}{w}}_3\Sigma \text{d}A\text{d}\eta \\ -\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{D_z}}{T}_{\mathrm{,3}}\theta \text{d}A\text{d}\eta -\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta \mathrm{.}\end{array}$ (3.11)

Now, we define a new function

$\begin{array}{c}\Gamma \left(z,t\right)={\displaystyle {\int }_z^{\infty }}\Phi \left(\xi ,t\right)\text{d}\xi \\ =\frac{\rho }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)w_{i,j}{w}_{i,j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)w_{i,i}^2\text{d}x\\ +K^{*}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\Sigma }_{,i}{\Sigma }_{,i}\text{d}x\text{d}\eta \\ +\frac{1}{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x.\end{array}$

(3.12)

From (3.11), we have

$\begin{array}{c}\Gamma \left(z,t\right)=-\nu {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{w}_{i,3}{\stackrel{\dot{}}{w}}_i\text{d}x\text{d}\eta +\left(\lambda +\nu \right){\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{w}_{j,j}{\stackrel{\dot{}}{w}}_3\text{d}x\text{d}\eta -K^{*}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,3}\theta \text{d}x\text{d}\eta \\ +{\gamma }_1{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_3\theta \text{d}x\text{d}\eta -M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{,3}\Sigma \text{d}x\text{d}\eta +{\gamma }_2{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_3\Sigma \text{d}x\text{d}\eta \\ -\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{,3}\theta \text{d}x\text{d}\eta -\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)T_{,i}{\theta }_{,i}\text{d}x\text{d}\eta .\end{array}$ (3.13)

Following the same procedures which have been used in Section 3 (see (3.6)-(3.10)), we have

$\Gamma \left(z,t\right)\le m_1\left(t\right)\left(-\frac{\partial \Gamma }{\partial z}\left(z\mathrm{,}t\right)\right)-\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{\mathrm{,3}}\theta \text{d}x\text{d}\eta -\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)T_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta \mathrm{.}$ (3.14)

By Hölder and the AG mean inequalities, from (3.14) we have

$\begin{array}{c}\Gamma \left(z,t\right)\le m_1\left(t\right)\left(-\frac{\partial \Gamma }{\partial z}\left(z,t\right)\right)+\frac{{\beta }^2}{2K^{*}{m}_1\left(t\right){\lambda }_1}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{,i}{T}_{,i}\text{d}x\text{d}\eta \\ +\frac{m_1\left(t\right)K^{*}}{2}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta +\frac{{\beta }^2}{2K^{*}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)T_{,i}{T}_{,i}\text{d}x\text{d}\eta \\ +\frac{K^{*}}{2}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\theta }_{,i}{\theta }_{,i}\text{d}x\text{d}\eta .\end{array}$ (3.15)

In view of the definitions of the functions $\Phi \left(z,t\right)$ and $\Psi \left(z,t\right)$, we introduce the functions

$\begin{array}{c}\tilde{\Phi }\left(z\mathrm{,}t\right)=\frac{\rho }{2}{\displaystyle {\int }_{R_z}}{\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i\mathrm{,}j}{w}_{i\mathrm{,}j}\text{d}x+\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}{w}_{i\mathrm{,}i}^2\text{d}x\\ +\frac{K^{\mathrm{<mo>\; *\; </mo>}}}{2}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\theta }_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{\Sigma }_{\mathrm{,}i}{\Sigma }_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +\frac{1}{2}{\displaystyle {\int }_{R_z}}\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x\text{d}\eta \mathrm{,}\end{array}$ (3.16)

and

$\begin{array}{c}\tilde{\Gamma }\left(z\mathrm{,}t\right)={\displaystyle {\int }_z^{\infty }}\tilde{\Phi }\left(\xi ,t\right)\text{d}\xi \\ =\frac{\rho }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\stackrel{\dot{}}{w}}_i{\stackrel{\dot{}}{w}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)w_{i\mathrm{,}j}{w}_{i\mathrm{,}j}\text{d}x\\ +\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)w_{i\mathrm{,}i}^2\text{d}x+\frac{K^{\mathrm{<mo>\; *\; </mo>}}}{2}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\theta }_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right){\Sigma }_{\mathrm{,}i}{\Sigma }_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +\frac{1}{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)\left[c{\theta }^2+2d\theta \Sigma +n{\Sigma }^2\right]\text{d}x\text{d}\eta \mathrm{.}\end{array}$ (3.17)

Then, inequality (3.15) may be rewritten as

$\begin{array}{c}\tilde{\Gamma }\left(z\mathrm{,}t\right)\le 3m_1\left(t\right)\left(-\frac{\partial \tilde{\Gamma }}{\partial z}\left(z\mathrm{,}t\right)\right)+\frac{{\beta }^2}{2K^{\mathrm{<mo>\; *\; </mo>}}{m}_1\left(t\right){\lambda }_1}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{\mathrm{,}i}{T}_{\mathrm{,}i}\text{d}x\text{d}\eta \\ +\frac{{\beta }^2}{2K^{\mathrm{<mo>\; *\; </mo>}}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)T_{\mathrm{,}i}{T}_{\mathrm{,}i}\text{d}x\text{d}\eta \mathrm{.}\end{array}$ (3.18)

Combining (2.1) and (2.2), we know

${\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}{T}_{\mathrm{,}i}{T}_{\mathrm{,}i}\text{d}x\text{d}\eta \le \frac{1}{K}E\left(\mathrm{0,}t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\mathrm{.}$ (3.19)

If we set

$\hat{E}\left(z\mathrm{,}t\right)={\displaystyle {\int }_z^{\infty }}E\left(\xi \mathrm{,}t\right)\text{d}\xi \mathrm{,}$ (3.20)

from (2.2), we have

$\hat{E}\left(z\mathrm{,}t\right)\le m_1\left(t\right)E\left(\mathrm{0,}t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\mathrm{.}$ (3.21)

From the definition of $E\left(z\mathrm{,}t\right)$ and (3.20), we have

$\begin{array}{c}\hat{E}\left(z\mathrm{,}t\right)=\frac{1}{2}\rho {\displaystyle {\int }_{R_z}}\left(\xi -z\right){\stackrel{\dot{}}{u}}_i{\stackrel{\dot{}}{u}}_i\text{d}x+\frac{\nu }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)u_{i,j}{u}_{i,j}\text{d}x\\ +\frac{\lambda +\nu }{2}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)u_{i,i}^2\text{d}x+K{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)T_{,j}{T}_{,j}\text{d}A\text{d}\eta \\ +M{\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)C_{,j}{C}_{,j}\text{d}A\text{d}\eta +{\displaystyle {\int }_{R_z}}\left(\xi -z\right)\left[\frac{c}{2}{T}^2+dCT+\frac{n}{2}{C}^2\right]\text{d}x.\end{array}$ (3.22)

So, we have the following inequality

${\displaystyle {\int }_0^t}{\displaystyle {\int }_{R_z}}\left(\xi -z\right)T_{\mathrm{,}j}{T}_{\mathrm{,}j}\text{d}A\text{d}\eta \le \frac{1}{K}{m}_1\left(t\right)E\left(\mathrm{0,}t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\mathrm{.}$ (3.23)

Inserting (3.19) and (3.23) back into (3.18), we have

$\tilde{\Gamma }\left(z\mathrm{,}t\right)\le m_1\left(t\right)\left(-\frac{\partial \tilde{\Gamma }}{\partial z}\left(z\mathrm{,}t\right)\right)+{\beta }^2{m}_3\left(t\right){\text{e}}^{-\frac{1}{m_1\left(t\right)}z}\mathrm{,}$ (3.24)

where

$m_3\left(t\right)=\frac{E\left(\mathrm{0,}t\right)}{2K{K}^{\mathrm{<mo>\; *\; </mo>}}{m}_1\left(t\right){\lambda }_1}+\frac{E\left(\mathrm{0,}t\right)}{2K{K}^{\mathrm{<mo>\; *\; </mo>}}}{m}_1\left(t\right)\mathrm{.}$ (3.25)

For any fixed $t_1>0$, setting $k_1=3m_1\left(t\right)$ and $k_2=\frac{1}{m_1\left(t\right)}$, from (3.24) we obtain

$\frac{\partial \tilde{\Phi }}{\partial z}\left(z\mathrm{,}{t}_1\right)+{\displaystyle {\int }_z^{\infty }}\tilde{\Phi }\left(\xi \mathrm{,}{t}_1\right)\text{d}\xi \le k_1\tilde{\Phi }\left(z\mathrm{,}{t}_1\right)+{\beta }^2{m}_3{\text{e}}^{-k_{2}z}\mathrm{,}$ (3.26)

where we have used the fact $\frac{\partial \tilde{\Phi }}{\partial z}\left(z\mathrm{,}{t}_1\right)<0$. To get the result we want, we let

$\Pi \left(z\mathrm{,}{t}_1\right)={\text{e}}^{-k_{1}z}\tilde{\Phi }\left(z\mathrm{,}{t}_1\right)+h{\displaystyle {\int }_z^{\infty }}{\text{e}}^{-k_1\xi }\tilde{\Phi }\left(\xi \mathrm{,}{t}_1\right)\text{d}\xi \mathrm{.}$ (3.27)

Thus, inequality (3.26) may be rewritten as

$\frac{\partial \Pi }{\partial z}\left(z\mathrm{,}{t}_1\right)+h\Pi \left(z\mathrm{,}{t}_1\right)\le {\beta }^2{m}_3{\text{e}}^{-\left(k_1+k_2\right)z}\mathrm{,}$ (3.28)

provided h satisfies the quadratic equation

$h^2-k_{1}h-1=0.$ (3.29)

We make the choice of

$h=h_0=\frac{k_1+\sqrt{k_1^2+4}}{2}.$ (3.30)

For this choice of h, to integrate (3.28), we have to consider the following two cases:

1) If $h_0-k_1-k_2=0$, an integration of (3.28) leads to

$\Pi \left(z\mathrm{,}{t}_1\right)\le \Pi \left(\mathrm{0,}{t}_1\right){\text{e}}^{-h_{0}z}+{\beta }^2{m}_{3}z{\text{e}}^{-h_{0}z}\mathrm{.}$ (3.31)

In light of (3.27), we have

$\tilde{\Phi }\left(z\mathrm{,}{t}_1\right)\le \Pi \left(\mathrm{0,}{t}_1\right){\text{e}}^{-\left(h_0-k_1\right)z}+{\beta }^2{m}_3{\text{e}}^{-\left(h_0-k_1\right)z}\mathrm{.}$ (3.32)

2) If $h_0-k_1-k_2\ne 0$, an integration of (3.28) leads to

$\Pi \left(z\mathrm{,}{t}_1\right)\le \Pi \left(\mathrm{0,}{t}_1\right){\text{e}}^{-h_{0}z}+\frac{{\beta }^2{m}_3}{h_0-k_1-k_2}\left[{\text{e}}^{\left(h_0-k_1-k_2\right)z}-1\right]{\text{e}}^{-h_{0}z}\mathrm{.}$ (3.33)

It is easy to proof that the second term on the right of (3.33) is positive either $h_0-k_1-k_2>0$ or $h_0-k_1-k_2<0$. In view of (3.27), we have

$\tilde{\Phi }\left(z\mathrm{,}{t}_1\right)\le \left[\Pi \left(\mathrm{0,}{t}_1\right)-\frac{{\beta }^2{m}_3}{h_0-k_1-k_2}\right]{\text{e}}^{-\left(h_0-k_1\right)z}+\frac{{\beta }^2{m}_3}{h_0-k_1-k_2}{\text{e}}^{-k_{2}z}\mathrm{.}$ (3.34)

In order to make inequalities (3.32) and (3.34) explicit, we need bound for $\Pi \left(\mathrm{0,}{t}_1\right)$. From the definition of $\Pi \left(z\mathrm{,}{t}_1\right)$ in (3.27), we may write

$\Pi \left(\mathrm{0,}{t}_1\right)=\tilde{\Phi }\left(\mathrm{0,}{t}_1\right)+h{\displaystyle {\int }_0^{\infty }}{\text{e}}^{-k_1\xi }\tilde{\Phi }\left(\xi \mathrm{,}{t}_1\right)\text{d}\xi \le \tilde{\Phi }\left(\mathrm{0,}{t}_1\right)+h\tilde{\Gamma }\left(\mathrm{0,}{t}_1\right)\mathrm{.}$ (3.35)

From (3.35), to bound $\Pi \left(\mathrm{0,}{t}_1\right)$ we only need to bound $\tilde{\Phi }\left(\mathrm{0,}{t}_1\right)$ and $\tilde{\Psi }\left(\mathrm{0,}{t}_1\right)$. From (3.11), we have

$\Phi \left(\mathrm{0,}{t}_1\right)=-\beta {\displaystyle {\int }_0^t}{\displaystyle {\int }_R{T}_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta }\le \frac{{\beta }^2}{2K^{\mathrm{<mo>\; *\; </mo>}}}{\displaystyle {\int }_0^t}{\displaystyle {\int }_R{T}_{\mathrm{,}i}{T}_{\mathrm{,}i}\text{d}x\text{d}\eta }+\frac{K^{\mathrm{<mo>\; *\; </mo>}}}{2}{\displaystyle {\int }_0^t}{\displaystyle {\int }_R{\theta }_{\mathrm{,}i}{\theta }_{\mathrm{,}i}\text{d}x\text{d}\eta }$ (3.36)

Using (3.12) and (3.16), we have

$\tilde{\Phi }\left(\mathrm{0,}{t}_1\right)\le \frac{{\beta }^2}{2K{K}^{\mathrm{<mo>\; *\; </mo>}}}E\left(\mathrm{0,}{t}_1\right)\mathrm{.}$ (3.37)

From (3.24) and using (3.37), we can get

$\tilde{\Psi }\left(\mathrm{0,}{t}_1\right)\le \frac{m_1{\beta }^2}{2K{K}^{\mathrm{<mo>\; *\; </mo>}}}E\left(\mathrm{0,}{t}_1\right)\mathrm{.}$ (3.38)

Combining (3.35), (3.37) and (3.38), we have

$\Pi \left(\mathrm{0,}{t}_1\right)\le {\beta }^2{m}_4\mathrm{,}$ (3.39)

where $m_4=\frac{E\left(\mathrm{0,}{t}_1\right)}{2K{K}^{\mathrm{<mo>\; *\; </mo>}}}+\frac{h{m}_{1}E\left(\mathrm{0,}{t}_1\right)}{2K{K}^{\mathrm{<mo>\; *\; </mo>}}}$. Combining the above discussions, we can conclude:

If $h_0-k_1-k_2=0$, we have

$\tilde{\Phi }\left(z\mathrm{,}{t}_1\right)\le {\beta }^2{m}_4{\text{e}}^{-\left(h_0-k_1\right)z}+{\beta }^2{m}_3{\text{e}}^{-\left(h_0-k_1\right)z}\mathrm{.}$ (3.40)

If $h_0-k_1-k_2\ne 0$, we have

$\tilde{\Phi }\left(z\mathrm{,}{t}_1\right)\le {\beta }^2\left[m_4-\frac{m_3}{h_0-k_1-k_2}\right]{\text{e}}^{-\left(h_0-k_1\right)z}+\frac{{\beta }^2{m}_3}{h_0-k_1-k_2}{\text{e}}^{-k_{2}z}\mathrm{.}$ (3.41)

Inequalities (3.40) and (3.41) exhibit not only exponential decay in z, but also show that the amplitude terms in (3.40) and (3.41) become small as $K\to K^{\mathrm{<mo>\; *\; </mo>}}$.

4. Conclusion

In view of the Equations (1.2) and (1.3), we may also obtain the continuous dependence on the coefficient M by employing the methods which have been used in Section 3. Our method is also valid to study other equations. In the future, we will use the method proposed in this paper to study the structural stability for the fluid flow in porous media. We think we will get some interesting results.

Funding

The work was supported by the science foundation of Huashang College Guangdong University of Finance & Economics (Grant # 2019HSDS28).