Next Article in Journal
Dynamical Stability in a Non-Hermitian Kicked Rotor Model
Next Article in Special Issue
Nonlocal Pseudo-Parabolic Equation with Memory Term and Conical Singularity: Global Existence and Blowup
Previous Article in Journal
Editorial Conclusion for the Special Issue “Measurements in Quantum Mechanics”
Previous Article in Special Issue
Global Existence to Cauchy Problem for 1D Magnetohydrodynamics Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Behavior of the Solution to Compressible Navier–Stokes System with Temperature-Dependent Heat Conductivity in an Unbounded Domain

Center of Applied Mathematics, Yichun University, Yichun 336000, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2023, 15(1), 112; https://doi.org/10.3390/sym15010112
Submission received: 25 November 2022 / Revised: 25 December 2022 / Accepted: 28 December 2022 / Published: 31 December 2022

Abstract

:
This paper concerns the one-dimensional compressible Navier–Stokes system with temperature-dependent heat conductivity in R with large initial data. We prove that velocity and temperature are uniformly bounded from below and above in time and space when the heat conductivity coefficient takes κ = κ ¯ ( 1 + θ b ) for all b > 5 2 . In addition, we show that the global solution is asymptotically stable as time tends to infinity.

1. Introduction

This paper concerns the Cauchy problem of compressible fluids in one space dimensions. The motion of a perfect polytropic ideal heat-conducting fluids can be written in the following form [1]:
ρ t + ( ρ u ) y = 0 , ( ρ u ) t + ( ρ u 2 + P ) y = ( μ u y ) y , ( ρ ( e + 1 2 u 2 ) ) t + ( ρ ( e + 1 2 u 2 ) u + P u ) y = ( κ e y ) y + ( μ u u y ) y ,
where t > 0 and y R are the time variable and spatial variable, respectively, where the unknown ρ 0 denotes the density of the flow, u the velocity, and e the internal energy. Both pressure P and internal energy e are generally related to the density and temperature of the flow according to the equations of state: P = P ( ρ , θ ) and e = e ( ρ , θ ) . Parameters μ = μ ( ρ , θ ) denote the viscosity coefficients, and κ = κ ( ρ , θ ) is the heat conductivity.
To solve the Cauchy problem, we transform Problem (1) into Lagrangian variables. To this end, we introduce the Lagrangian symmetry variable
x = y ( t ) y ρ ( t , z ) d z ,
where y ( t ) is the particle path satisfying y ( t ) = u ( t , y ( t ) ) . The Lagrangian version of System (1) can be written as
{ (2a) v t = u x , (2b) u t + P x = ( μ u x v ) x , (2c) ( e + u 2 2 ) t + ( P u ) x = ( κ θ x v + μ u u x v ) x , (2d) P = R θ v , e = c v θ .
We consider a perfect gas for Navier–Stokes flow in this paper, that is,
P = R θ v , e = c v θ ,
where R is a positive constant, and c v is the heat capacity of the gas. System (2) is supplemented with the following initial condition:
( u , v , θ ) | t = 0 = ( u 0 , v 0 , θ 0 ) , x R ,
and the far-field condition:
lim | x | ( v ( x , t ) , u ( x , t ) , θ ( x , t ) ) = ( 1 , 0 , 1 ) , t > 0 .
Let us review some results on System (2) in different situations. When μ and κ were constants, the existential results in bounded domains for large initial data were obtained by Kazhikhov et al. [2,3,4]. Regarding initial boundary value problems in bounded domains, see [2,5,6,7,8,9] for a thorough discussion of System (2) with initial Condition (4) and far-field Condition (5). Furthermore, the existence and uniqueness of global solutions, and the regularity are known [2,5,6,7,8,9,10]. Moreover, the asymptotic behavior of the global solution was studied as time tended to infinity; see [11,12,13,14], among others. For the Cauchy problem, the global existence of a solution was obtained by Kazhikhov [15]; then, Li [16] gave the asymptotic behavior of solutions to System (2) with initial Condition (4) and far-field Condition (5).
We could obtain compressible Navier–Stokes Equations (2) from the celebrated Boltzmann equations for monatomic gas with a slab symmetry by using the Chapman–Enskog expansion. Then, viscosity coefficient μ and heat conductivity coefficient κ are functions of density and temperature; see Chapman and Cowling [17] or Vincenti and Kruger ([18], Chapter X) for a thorough discussion of these issues. When the coefficients depended on special volume and temperature, for the one-dimensional full compressible Navier–Stokes equations of ideal polytropic gas whose viscosity coefficient and heat conductivity coefficient satisfying μ = μ ¯ h ( v ) θ b , κ = κ ¯ h ( v ) θ b , Liu, Yang, et al. in [19] obtained the global nonvacuum classical solutions with a smallness mechanism (i.e., γ 1 small). Wang and Zhao in [20] obtained the global nonvacuum classical solutions with smallness assumptions for b. Later, in 2016, Wang and Zhao [21] gave the large-time behavior of the solutions under the assumptions that C h ( v ) v l 1 + v l 2 , h ( v ) 2 v C h ( v ) 3 and that b was small enough. Duan, Guo, et al. [22] proved the existence and uniqueness of a strong global solution for ideal polytropic gas flow, with μ = 1 + ρ α and κ = θ β . Kazhikhov [15] gave frameworks when μ and κ are constants. However, if the viscosity coefficient depends on temperature, Kazhikhov’s method is invalid. Li, Shu, et al. [23] proved the global existence of strong solutions to a compressible Navier–Stokes system with degenerate heat conductivity in unbounded domains. However, the asymptotic behavior of a solution with large initial data is still open.
When viscosity was a positive constant, and only heat conductivity depended on temperature, i.e.,
μ = μ ¯ , κ = κ ¯ θ b ,
Jenssen and Karper [24] proved the global existence of a weak solution to initial-boundary value problem (IBVP) (2) under the assumption that b [ 0 , 3 2 ) ; Pan and Zhang [25] extended it to b [ 0 , ) . Li and Guo [1] established the global existence of strong and classical solutions to free boundary Problem (2) for b [ 0 , ) , and the expanding rates of the interface were also studied. Recently, Li, Shu, et al. [23] proved the global existence of a solution to Cauchy Problem (2) for b [ 0 , ) . Chen and Zhang [26] proved global existence to free boundary problems. Cai, Chen, et al. [27] obtained the asymptotic behavior of the initial boundary value problem of System (2). However, the asymptotic behavior to the Cauchy problem is still open and our focus. The research on numerical and applications in engineering to system of (2) and it’s related models, see [28,29,30,31].
The mission of this paper is to establish the uniform bounds from below and above of velocity and temperature to the Cauchy problem, and the large-time behavior of strong solutions with κ = κ ¯ ( 1 + θ b ) .
Notations:
(1)
For p 1 , L p = L p ( R ) denotes the L p space with the norm · L p . For k 1 and p 1 , W k , p = W k , p ( R ) denotes the Sobolev space, whose norm is denoted as · W k , p , H k = W k , 2 ( R ) . For k 1 and p 1 , D k , p ( R ) denotes the homogeneous Sobolev space, the norm of f D k , p ( R ) is f k L p ( R ) . Q T = [ 0 , T ] × R .
(2)
For the sake of simplicity, we denote various positive constants independent of time T and depending on time T with C and C ( T ) , which may be different at different occurrences.
Definition 1.
(Global strong solution) For any ( x , t ) ( [ 0 , ) × R ) , ( v , u , θ ) is called a global strong solution if
v 1 C ( [ 0 , ) , H 1 ( R ) ) , θ 1 C ( [ 0 , ) , H 2 ( R ) ) L 2 ( [ 0 , ) , H 1 ( R ) ) , u L ( [ 0 , ) , H 2 ( R ) ) L 2 ( [ 0 , ) , W 2 , 2 ( R ) ) , v t , u t , θ t L 2 ( [ 0 , ) , D 1 , 2 ( R ) ) ,
and ( v , u , θ ) satisfies both System (2) almost everywhere in R × ( 0 , ) and Initial Value (4) almost everywhere in R .
The existence and uniqueness of local solution can be proven with a fixed-point theorem; see Tani [32], who proved the existence of local solution if the initial (4) and far-field Condition (5) are satisfied, and μ , κ are locally Lipschitz-continuous functions on ( v , θ ) . As a special case of the result in [32], the following theorem gives the local existence for our problem.
Theorem 1.
Assume that μ and κ satisfy (6) for some positive constants μ ¯ and κ ¯ . If the initial data ( v 0 , u 0 , θ 0 ) ( x ) are compatible with far-field Condition (5), satisfying
( v 0 1 , u 0 , θ 0 1 ) ( x ) H 1 × H 2 × H 2 ,
and there are constants v ̲ , v ¯ , θ ̲ , θ ¯ such that
0 < v ̲ v 0 ( x ) v ¯ , 0 < θ ̲ θ 0 ( x ) θ ¯ ,
then there exists a unique local strong solution ( v , u , θ ) ( x , t ) to (2) on R × [ 0 , T 1 ] for some C > 0 depending on the initial data, and T 1 satisfies
C 1 θ ( x , t ) C ( T 1 ) , C 1 v ( x , t ) C ( T 1 ) , ( v 1 , u , θ 1 ) ( · , t ) H 1 ( R ) 2 + 0 t ( v 1 , u , θ 1 ) ( · , s ) H 1 ( R ) 2 d s C ( T 1 ) , ( u , θ 1 ) ( · , t ) H 2 ( R ) 2 + 0 t ( v x t , u x t , u x x θ x t , θ x x ) ( · , s ) L 2 ( R ) 2 d s C ( T 1 ) .
if the initial data further satisfy
v 0 ( x ) C 1 + α , u 0 ( x ) C 2 + α , θ 0 C 2 + α ,
then v C 1 + α , α 2 ( R × [ 0 , T 1 ] ) , u C 2 + α , 1 + α 2 ( R × [ 0 , T 1 ] ) , and θ C 2 + α , 1 + α 2 ( R × [ 0 , T 1 ] ) .
Thanks to this local existence result, the existence of a global solution is established by extending the local solution with the help of the global a priori estimates stated in (12) (see Theorem 2). It is clear that (12) is sufficient to extend the local strong solution to global one by a standard continuity argument.
The following are the main results of this paper. Some uniform estimate results and the large-time behavior of the solutions are obtained when the heat conductivity coefficient is in nondegenerate form with the temperature.
Theorem 2.
Assume that the initial data ( v 0 , u 0 , θ 0 ) satisfy (8), (9), κ = κ ¯ ( 1 + θ b ) (nondegenerate case) for b ( 5 2 , ) . Let ( v , u , θ ) be a solution to (2)–(4) together with far-field Condition (5). For any T > 0 , there exists a unique global strong solution ( v , u , θ ) satisfying
( v 1 , u , θ 1 ) ( · , t ) H 1 ( R ) 2 + 0 t ( v 1 , u , θ 1 ) ( · , s ) H 1 ( R ) 2 d s C , ( u , θ 1 ) ( · , t ) H 2 ( R ) 2 + 0 t ( v x t , v x x , u x t , u x x , θ x t , θ x x ) ( · , s ) L 2 ( R ) 2 d s C .
Moreover, there exists a positive constant C depending only on μ , κ , R , c v , v ̲ , θ ̲ , and the initial value; the following uniform estimate holds
C 1 θ ( x , t ) C , C 1 v ( x , t ) C ,
and large-time behavior is obtained
lim t ( ( v 1 , u , θ 1 ) ( t ) L p ( R ) + ( v x , u x , θ x ) ( t ) L 2 ( R ) ) = 0 ,
for any p ( 2 , ] .
Remark 1.
Jiang [12] and Li [16] proved the results in Theorem 2 when κ was a constant. Jiang obtained the positive upper and lower bounds of v ( x , t ) , and Li proved that θ ( x , t ) was bounded from below and above, and the solution was asymptotically stable as time tended to infinity for large initial data.
Remark 2.
The global existence of a solution and large time were obtained in Theorem 2 for κ = κ ¯ ( 1 + θ b ) . The global existence for κ = κ ¯ θ b could also be obtained, but Large Time (14) failed for this case in our method.
We now outline the main ideas and difficulties in our problem compared to previous results. The existence of strong solutions can be easily obtained due to pioneering works, e.g., Tani [32], Kazhikhov [4], and Jesssen and Karper [24]. For the large-time behavior of such a solution, obtaining the uniform positive lower and upper bounds of v ( x , t ) and θ ( x , t ) is a great challenge due to the strong nonlinearity of κ = κ ¯ ( 1 + θ b ) . Jiang obtained uniform positive lower and upper bounds on v ( x , t ) in [12] with a decent localized version of the expression for v ( x , t ) when κ was a constant. Li and Liang deduced the uniform positive lower and upper bounds on temperature θ ( x , t ) in [16] with a smart test function method. However, there methods could not be applied to our case, since it is difficult to obtain the uniform bounds of the high-order estimate ( θ x L 2 ( R ) ), and bounds of θ from below and above. To overcome such a difficulty, motivated by [1,25], we obtained the high-order estimate Y ( t ) = sup 0 t T θ b θ x L 2 with an iterative method. The crucial techniques of proofs in [25] could not be adapted directly here since their arguments depend on bounded domain and boundary conditions that were different from ours, and we could not obtain L P ( p 1 ) norm of θ 1 under the far-field condition in this paper with an unbounded domain. In combination with the above methods in the literature, we discuss it withg a space separation technique and iterative method that could obtain the global existence of a solution. Then, combining the lower bound of the temperature when t [ 0 , T 0 ] induced by the comparison principle and the lower bound when t ( T 0 , ) obtained from a well-designed continuation argument for some suitable fixed T 0 [ 0 , ) , the positive pointwise boundedness of θ ( x , t ) from below and above independent of time and the large-time decay behavior of solutions could be obtained.
The rest of the paper is organized as follows. In Section 2, we give some a priori estimates, and prove the uniform positive lower and upper bounds of v ( x , t ) independent of time. In Section 3, on the basis of the local existence of the solutions and the a priori estimates in Section 2, we prove the global existence of solution with a standard continuity argument. In Section 4, we give the asymptotic behavior of the global solution for κ = κ ¯ ( 1 + θ b ) .

2. A Priori Estimates

In this section, we perform a sequence of estimates. We proved that volume v ( x , t ) was pointwise bounded from below and above independent of time. This is a key step in the proof of both the global existence and asymptotic behavior of the solution. Assume that ( v , u , θ ) ( x , t ) is the unique strong solution of (2), defined on R × [ 0 , T 0 ] for some T 0 > 0 .
Lemma 1.
There are positive constants e 0 and C independent of T, such that
sup 0 t < R 1 2 u 2 + R ( v ln v 1 ) + c v ( θ ln θ 1 ) d x + μ ¯ 0 R u x 2 v θ d x d t + κ ¯ 0 R ( 1 + θ b ) θ x 2 θ 2 v d x d t e 0 ,
Let Ω M ( t ) = { x R | θ ( x , t ) M > 1 } ; we derive from (15) that
  • Ω M | θ 1 | d x C Ω M ( θ ln θ 1 ) d x C ,
  • R / Ω M | θ 1 | 2 d x C R / Ω M ( θ ln θ 1 ) d x C .
Proof. 
By using Equation (2c) and a far-field condition, we obtain after a straightforward calculation that
c v θ t + R u x θ v = κ ¯ ( 1 + θ b ) θ x v x + μ ¯ ( u x ) 2 v .
Multiplying (2a) by R ( 1 v 1 ) , (2b) by u, (16) by ( 1 θ 1 ) , and adding them together, we obtain
1 2 u 2 + R ( v ln v 1 ) + c v ( θ ln θ 1 ) t + μ ¯ u x 2 v + κ ¯ ( 1 + θ b ) θ x 2 θ 2 v = μ ¯ ( u u x v ) x R ( u θ v ) x + R u x + κ ¯ ( 1 θ 1 ) ( 1 + θ b ) θ x 2 v x .
Using Taylor’s theorem, (8) and Sobolev’s imbedding theorem ( H 1 L ), we have
R 1 2 u 0 2 + R ( v 0 ln v 0 1 ) + c v ( θ 0 ln θ 0 1 ) d x C ( 1 + ( v 0 1 , u 0 , θ 0 1 ) H 1 2 ) .
Integrating (17) over R and using far-field Condition (5) obtains (15). The proof of Lemma 1 is finished. □
For some positive integer k, let ϕ W 1 , ( R ) be defined by
ϕ = 1 , x k + 1 ; k + 2 x , k + 1 x k + 2 ; 0 , x k + 2 .
For simplicity, we denote Ω k : = ( k , k + 1 ] for Cauchy problem. Then, the bounds of v ( x , t ) dependent on T can be obtained. We prove the pointwise bounds on a specific volume into two parts, when t [ 0 , T ] and t ( T , ) for some suitable T.
Lemma 2.
There exists a positive constant C, such that
C 1 ( T ) v ( x , t ) ,
for all ( x , t ) R × [ 0 , T ] .
Proof. 
For any x Ω k , we have the following local representation via Lemma 1:
k k + 1 [ ( v ln v 1 ) + ( θ ln θ 1 ) ] d x e 0 ,
which, together with Jensen’s inequality, yields
α 1 k k + 1 v ( x , t ) d x α 2 , α 1 k k + 1 θ ( x , t ) d x α 2 ,
where 0 < α 1 < 1 < α 2 are two roots of
x ln x 1 = e 0 .
Moreover, it follows from (19) that, for any t > 0 , there exists some b k ( t ) [ k , k + 1 ] , such that
( v ln v 1 ) + ( θ ln θ 1 ) ( b k ( t ) , t ) e 0 ,
which implies
α 1 v ( b k ( t ) , t ) α 2 , α 1 θ ( b k ( t ) , t ) α 2 .
Letting σ μ ¯ u x v R θ v = μ ¯ ( ln v ) t R θ v , we write (2b) as
u t = σ x .
Multiplying (2b) by ϕ gives
[ ϕ u ] t = μ ¯ u x v R θ v ϕ x ϕ x μ ¯ u x v R θ v .
Integrating over ( x , ) ( x Ω k ) with respect to x and recalling (2a) and the definition of ϕ , we have
x [ ϕ u ] t d y = μ ¯ u x v R θ v + x μ ¯ u x v R θ v ϕ x d y = μ ¯ [ ln v ] t R θ v k + 1 k + 2 μ ¯ u x v R θ v d y , x Ω k .
Furthermore, integrating over [ 0 , t ] , one has
x ( u 0 u ) ϕ d y = μ ¯ ( ln v ln v 0 ) 0 t R θ v d s 0 t k + 1 k + 2 μ ¯ u x v R θ v d y d s , x Ω k .
Denote
B ( x , t ) = v 0 exp 1 μ ¯ x ( u 0 ( y ) u ( y , t ) ) ϕ ( y ) d y , Y ( t ) = exp 1 μ ¯ 0 t k + 1 k + 2 μ ¯ u x v R θ v d y d s .
Taking the exponential on both sides of (24), the following relation appears:
1 B ( x , t ) Y ( t ) = 1 v ( x , t ) exp 0 t R θ μ ¯ v d s , x Ω k , t 0 .
Multiplying (26) by R θ ( x , t ) μ ¯ and integrating over ( 0 , t ) , we infer
exp R μ ¯ 0 t θ ( x , s ) v ( x , s ) d s = 1 + R μ ¯ 0 t θ ( x , s ) B ( x , s ) Y ( s ) d s .
Substituting the above identity into (26), we obtain
v ( x , t ) = B ( x , t ) Y ( t ) + R μ ¯ 0 t B ( x , t ) Y ( t ) B ( x , s ) Y ( s ) θ ( x , s ) d s , x Ω k , t 0 .
Since
| x ( u ( y , t ) u 0 ( y ) ) ϕ ( y ) d y | Ω k u 2 d y 1 2 + Ω k u 0 2 d y 1 2 C ( e 0 ) ,
we deduce from (25)
C 1 ( e 0 ) B ( x , t ) C ( e 0 ) .
Moreover, integrating (27) with respect to x over [ k , k + 1 ] gives
1 Y ( t ) k k + 1 v ( x , t ) d x = k k + 1 B ( x , t ) 1 + 0 t θ ( x , τ ) B ( x , τ ) Y ( τ ) d τ d x .
Hence, we have
1 Y ( t ) C ( α 1 ) + C ( α 1 , α 2 , e 0 ) 0 t 1 Y ( τ ) d τ ,
where (15), (20), and (28) were used, as well as the following simple fact:
k k + 1 θ ( x , τ ) B ( x , τ ) B ( x , t ) d x C ( e 0 ) k k + 1 θ ( x , τ ) d x C ( α 2 , e 0 ) .
Applying the Grönwall’s inequality to (29) gives
1 Y ( t ) C ( T , α 1 , α 2 , e 0 ) ,
which implies that, for any positive integer k and ( x , t ) [ k , k + 1 ] × [ 0 , T ] , from (27), we have
C 1 ( T , α 1 , α 2 , e 0 ) v ( x , t ) .
From (29), there exists a suitable constant N > 0 , such that v ( x , t ) C 1 ( T ) when k < N . Combining the fact v ( x , t ) 1 as | x | , one has v ( x , t ) 1 ϵ when k N and ϵ small enough. So, the bounds of v ( x , t ) from below are obtained in R × [ 0 , T ] . □
Lemma 3.
For any t [ 0 , T ] , and positive constants C ( T ) , it holds that
θ C ( T ) .
Proof. 
Let θ ¯ = 1 θ , and rewrite Equation (16) as follows:
θ ¯ t θ ¯ 2 = R u x c v v θ ¯ 1 c v κ ( θ ¯ ) θ ¯ x v θ ¯ 2 x + μ ¯ u x 2 c v v .
So,
θ ¯ t = R θ ¯ u x c v v + θ ¯ 2 c v κ ( θ ¯ ) θ ¯ x v θ ¯ 2 x μ ¯ θ ¯ 2 u x 2 c v v = κ ( θ ¯ ) θ ¯ x c v v x 2 κ ( θ ¯ ) θ ¯ x 2 c v v θ ¯ μ ¯ θ ¯ 2 u x 2 c v v + R θ ¯ u x c v v = κ ( θ ¯ ) θ ¯ x c v v x 2 κ ( θ ¯ ) θ ¯ x 2 c v v θ ¯ μ ¯ θ ¯ 2 c v v u x 2 R u x μ ¯ θ ¯ = κ ( θ ¯ ) θ ¯ x c v v x 2 κ ( θ ¯ ) θ ¯ x 2 c v v θ ¯ μ ¯ θ ¯ 2 c v v u x R 2 μ ¯ θ ¯ 2 + R 2 4 μ ¯ c v v ,
which implies
θ ¯ t κ ( θ ¯ ) θ ¯ x c v v x + R 2 4 μ ¯ c v v κ ( θ ¯ ) θ ¯ x c v v x + C ( T ) .
Define the operator L : = t + κ ( · ) v ( · ) x x and
L θ ˜ < 0 on [ 0 , ) × Ω , θ ˜ t = 0 0 on Ω , θ ˜ x 0 on [ 0 , ) ,
where θ ˜ ( x , t ) = C ( T ) t + max Ω ¯ θ ¯ 0 θ ¯ ( x , t ) ; then, with the comparison theorem, we obtain
min ( x , t ) Q ¯ T θ ˜ ( x , t ) 0 ,
which implies
θ ( x , t ) 1 C ( T ) t + max Ω ¯ θ ¯ 0 C ( T ) .
This completes the proof of Lemma 3. □
Now, in order to obtain the uniform upper and lower bounds of v ( x , t ) , we first show the exponential decay of Y ( t ) , and use Representation (27) to obtain the following uniform bounds on v ( x , t ) .
Lemma 4.
There exists a positive constant C independent of t, such that
C 1 v ( x , t ) C , x R , t 0 .
Proof. 
By using (28), one can first prove, by repeating the argument used in [12], the following estimates:
s t inf [ k + 1 , k + 2 ] θ ( s , · ) d s C t s C
for all 0 s t T . Then, one can obtain, with Jenssen’s inequality,
s t k + 1 k + 2 μ ¯ u x v R θ v d y d s C s t k + 1 k + 2 u x 2 θ v d y d s R 2 s t k + 1 k + 2 θ v d y d s C R 2 s t inf [ k + 1 , k + 2 ] θ ( s , · ) k + 1 k + 2 1 v d y d s C R 2 s t inf [ k + 1 , k + 2 ] θ ( s , · ) k + 1 k + 2 v d y 1 d s C R 2 α 2 s t inf [ k + 1 , k + 2 ] θ ( s , · ) d s C t s C , 0 s t T .
Recalling the definition of Y ( t ) yields that
Y ( t ) C e C t , Y ( t ) Y ( s ) C e C ( t s ) , t 0 , t s 0 .
On the other hand, for a point b k ( t ) [ k , k + 1 ] via Lemma 1 implies that
| θ 1 2 ( x , t ) θ 1 2 ( b k ( t ) , t ) | 1 2 b k ( t ) x θ 1 2 | θ x | d x 1 2 k k + 1 θ x 2 v θ 2 d x 1 2 k k + 1 v θ d x 1 2 α 2 2 Ω k θ x 2 v θ 2 d x 1 2 max [ k , k + 1 ] v 1 / 2 ( · , t ) ,
k = 0 , ± 1 , · · · . By using Young’s inequality, we have
C ( α 1 , b ) C ( α 2 ) Ω k θ x 2 v θ 2 d x max Ω ¯ k v ( · , t ) θ ( x , t ) C ( α 2 ) ( α 2 , b ) + C Ω k θ x 2 v θ 2 d x max Ω ¯ k v ( · , t ) .
Hence, inserting (28), (41), and (42) into (28), one has
v ( x , t ) C e C t + C 0 t e C ( t s ) C ( α 2 , b ) + C Ω k θ x 2 v θ 2 d x max Ω ¯ k v ( · , t ) d s ,
Applying Gronwall’s inequality and (15), we conclude
v ( x , t ) C , x Ω k , t 0 .
Now, we prove the lower bound of v ( x , t ) independent with T. The proof is divided into two parts, when t [ 0 , T 0 ] and t ( T 0 , ) , for some suitable fixed T 0 [ 0 , ) . When t [ 0 , T 0 ] , we know that v ( x , t ) C ( T ) via Lemma 2. Regarding t ( T 0 , ) , integrating (27) over [ k + 1 , k + 2 ] , using the estimate in Lemma 3, we have
α 1 C e C t + C 0 t Y ( t ) Y ( s ) d s .
It follows from (20), (28), and (41)–(44) that
v ( x , t ) C 4 0 t Y ( t ) Y ( s ) ( C 5 C 6 R θ x 2 v θ 2 d x ) d s C 7 C 8 e C t C 9 0 t / 2 + t / 2 t e C ( t s ) R θ x 2 v θ 2 d x d s C 10 C 11 e C t C 9 e ( C t ) / 2 0 t / 2 R θ x 2 v θ 2 d x d s C 9 t / 2 t R θ x 2 v θ 2 d x d s C 10 i = 1 3 J i C 12 .
Obviously, via (15), we have J 1 , J 2 0 as t . We also have
J 3 = C 9 0 t R θ x 2 v θ 2 d x d s C 9 0 t / 2 R θ x 2 v θ 2 d x d s 0 t .
So, we can take a large enough T 0 to ensure that C 12 > 0 .
Therefore, by combining (43), (29), and (45), one has
C 1 ( k ) v ( x , t ) C ( k ) ,
Due to the far-field condition, we obtain (39). This completes the proof of Lemma 4 (i.e., part of (13) in Theorem 2). □

3. Proof of Global Existence

In this section, we apply the results obtained in Section 2 to prove Theorem 2. Motivated by [1], we give the estimate on sup 0 t T R ( 1 + θ 2 b ) θ x 2 d x , which is the key step in the proof of Theorem 2.
In our case, nonlinearity κ on θ requires further attention on the control of θ . For this purpose, one of the main ingredients in this paper is the following lemma-refined estimates on temperature. In order to obtain the higher-order estimates to the solution, we follow the framework introduced in [1], and define the following two functionals
Z ( t ) = sup 0 t T R u x x 2 ( x , t ) d x , Y ( t ) = sup 0 t T R ( 1 + θ 2 b ) θ x 2 ( x , t ) d x .
These two functionals are useful in depicting the tangled relations of the higher order and upper bound of θ . First, we have that the following lemma holds.
Lemma 5.
For some positive constant, we have
sup R × [ 0 , T ] θ C + C Y 1 2 b + 3 δ ,
sup R × [ 0 , T ] | u x | C + C Z 3 8 .
for any ( x , t ) Q T .
Proof. 
For a small constant δ , via the Gagliardo–Nirenberg inequality, we infer that
sup R | θ 1 | 2 b + 3 δ = sup R x | θ 1 | 2 b + 3 δ x d y C R | θ 1 | 2 b + 2 δ | θ x | d x C R ( 1 + θ 2 b ) θ x 2 d x + C Ω M R / Ω M | θ 1 | 2 b + 4 2 δ d x C R ( 1 + θ 2 b ) θ x 2 d x + C sup R | θ 1 | 2 b + 3 2 δ Ω M | θ 1 | d x + C sup R | θ 1 | 2 b + 2 2 δ R / Ω M | θ 1 | 2 d x C + ϵ sup R | θ 1 | 2 b + 3 δ + C Y ,
which implies
sup R θ sup R | θ 1 | + 1 C + C Y 1 2 b + 3 δ ,
where we used the fact that | θ 1 | 2 b C ( 1 + θ 2 b ) , Lemma 1 and Young inequality.
Regarding (48), we have
sup R | u x | 2 R u x 2 d x + 2 R | u x u x x | d x R u 2 d x 1 2 R u x x 2 d x 1 2 + 2 R u x 2 d x 1 2 R u x x 2 d x 1 2 C Z 1 2 + 2 R u 2 d x 1 4 R u x x 2 d x 1 4 R u x x 2 d x 1 2 C + C Z 3 4 .
This completes the proof of Lemma 5. □
In addition, we have the following key estimate.
Lemma 6.
For a positive constant C, 0 t T and b 2 , it holds that
0 T sup R | θ 1 | b + 1 d t C
Proof. 
Through the Gagliardo–Nirenberg inequality and Lemma 1, we infer that
0 T sup R | θ 1 | b + 1 d t = 0 T sup R ( | θ 1 | b + 1 2 ) 2 d t 0 T R | θ 1 | b + 1 d x 1 2 R | θ 1 | b 1 θ x 2 d x 1 2 d t 0 T Ω M | θ 1 | b + 1 d x 1 2 sup R | θ 1 | 1 2 R | θ 1 | b 2 θ x 2 d x 1 2 d t + 0 T R / Ω M | θ 1 | b + 1 d x 1 2 sup R | θ 1 | R | θ 1 | b 3 θ x 2 d x 1 2 d t ε 0 T sup R | θ 1 | b + 1 d t + C Q T | θ 1 | b 2 θ x 2 d x d t + C Q T | θ 1 | b 3 θ x 2 d x d t ε 0 T sup R | θ 1 | b + 1 d t + C Q T ( 1 + θ b ) θ 2 θ x 2 d x d t ε 0 T sup R | θ 1 | b + 1 d t + C .
Here, we use the fact that, for b 2 and κ = κ ¯ ( 1 + θ b ) , we have
Q T θ b 2 θ x 2 d x d t Q T 1 + θ b θ 2 d x d t ,
and for b 1
Q T θ b 3 θ x 2 d x d t C Q T θ x 2 θ 2 d x d t + C Q T θ b θ x 2 θ 2 d x d t .
It yields that
0 T sup R | θ 1 | b + 1 d t C .
Then, the proof of Lemma 6 is completed. □
The following lemma gives estimates on the L 2 norm of v x .
Lemma 7.
For any t [ 0 , T ] , there exists a constant C independent of time; it holds that
sup 0 t T R v x 2 d x + Q T θ v x 2 d x d t C + C Y 1 2 b + 3 δ .
Proof. 
First, integrating (2b) multiplied by v x v over R , we obtain that, after using (2a),
μ ¯ 2 d d t R v x 2 v 2 d x = R R θ v x v x v d x + R u t v x v d x = R R θ x v x v 2 d x R R θ v x 2 v 3 d x + d d t R u v x v d x + R u x 2 v d x C R θ x 2 v θ d x R 2 R θ v x 2 v 3 d x + d d t R u v x v d x + R u x 2 v d x ,
which together with (15), which yields
R v x 2 v 2 d x + Q T θ v x 2 v 3 d x d t C + sup Q T θ Q T u x 2 v θ d x d t C + C Y 1 2 b + 3 δ .
We have the following relationship between the high-order estimates on u ( x , t ) and Y.
Lemma 8.
For any t [ 0 , T ] , we have following estimate:
sup 0 t T R u x 2 d x + Q T u x x 2 d x d t C + C Y 3 2 b + 3 δ .
Proof. 
We rewrite the momentum equation into the following form:
u t μ ¯ u x x v = μ ¯ u x v x v 2 R θ x v + R θ v x v 2 .
Multiplying (57) by u x x , and integrating in x over R , one has
sup 0 t T R u x 2 d x + Q T μ ¯ u x x 2 v d x d t C + | Q T μ ¯ u x v x v 2 u x x d x d t | + | Q T R θ x v u x x d x d t | + | Q T R θ v x v 2 u x x d x d t | C + μ ¯ 4 Q T u x x 2 v d x d t + C Q T ( u x 2 v x 2 + θ x 2 + θ 2 v x 2 ) d x d t .
According to Equation (2a) and the far-field conditions, one has
0 T sup R u x 2 d t = 2 0 T sup R x u x u x x d y d t C 0 T R u x 2 d x 1 2 R u x x 2 d x 1 2 d t .
Hence, via (15), Young’s inequality, and the uniform boundedness of v ( x , t ) , we have the term
Q T u x 2 v x 2 d x d t 0 T sup R u x 2 d t sup 0 t T R v x 2 d x C 0 T R u x 2 d x 1 2 R u x x 2 d x 1 2 d t sup 0 t T R v x 2 d x ε Q T u x x 2 v d x d t + C Q T u x 2 d x d t sup 0 t T R v x 2 d x 2 ε Q T u x x 2 v d x d t + C Y 3 2 b + 3 δ + C .
For the other two terms, we have the following estimate
Q T ( θ x 2 + θ 2 v x 2 ) d x d t sup Q T θ 2 Q T θ x 2 v θ d x d t + sup Q T θ Q T θ v x 2 d x d t C sup Q T θ 2 + C sup Q T θ Y 1 2 b + 3 δ C Y 2 2 b + 3 δ + C ,
Here, we use (15) and Lemma 7. Then, (58) shows that Lemma 8 holds. □
The relation of Y and Z is given in the following lemma.
Lemma 9.
For any ( x , t ) Q T and positive constant C, we have
Y + Q T ( 1 + θ b ) θ t 2 d x d t C + ε 1 Z .
Proof. 
Define K ( v , θ ) = θ v + θ 1 + b ( 1 + b ) v . Multiplying (16) by K t , integrating over Q T and by parts, we have
Q T θ t K t d x d t + κ ¯ c v Q T ( 1 + θ b ) θ x v K x t d x d t = Q T μ ¯ u x 2 c v v R θ u x c v v K t d x d t .
Then, we compute
K t = ( 1 + θ b ) θ t v θ v t v 2 θ 1 + b v t ( 1 + b ) v 2 ,
K x = ( 1 + θ b ) θ x v θ v x v 2 θ 1 + b v x ( 1 + b ) v 2 ,
K x t = ( ( 1 + θ b ) θ x v ) t + 2 θ v x u x v 3 + 2 θ 1 + b v x u x ( 1 + b ) v 3 ( 1 + θ b ) θ t v x v 2 θ u x x v 2 θ 1 + b u x x ( 1 + b ) v 2 ( ( 1 + θ b ) θ x v ) t + K ˜ .
Hence, we have
R ( 1 + θ 2 b ) θ x 2 v 2 d x + Q T ( 1 + θ b ) θ t 2 v d x d t
C + Q T θ 1 + b θ t u x ( 1 + b ) v 2 + θ θ t u x v 2 d x d t κ ¯ c v Q T ( 1 + θ b ) θ x v K ˜ d x d t + Q T μ ¯ u x 2 v R θ u x v K t d x d t
C + i = 1 3 I i .
Next, we give the estimate on I 1 , I 2 , I 3 . Using Lemmas 1 and 5, the first term can be estimated as follows:
I 1 = Q T θ 1 + b θ t u x ( 1 + b ) v 2 + θ θ t u x v 2 d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + C Q T ( θ b + 2 + θ 2 ) u x 2 d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + sup Q T ( θ b + 3 + θ 3 ) Q T u x 2 θ d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + C Y b + 3 2 b + 3 δ + C .
Then, we turn to estimate I 2 . We divide the proof into three parts through the lemmas proved in Section 2, Lemmas 5–8, and the interpolation inequality.
I 2 C Q T ( ( θ + θ 1 + 2 b ) | θ x u x v x | + ( 1 + θ 2 b ) | θ x θ t v x | + ( θ + θ 1 + 2 b ) | θ x u x x | ) d x d t i = 1 3 I 2 j .
We now give the estimate on I 2 j ( j = 1 , 2 , 3 ) term by term. For term I 21 , when b 1 , one has
I 21 C Q T ( ( θ + θ 1 + 2 b ) | θ x u x v x | d x d t sup Q T | u x | [ Q T θ 3 θ x 2 θ 2 d x d t 1 2 Q T θ v x 2 d x d t 1 2 + Q T θ 3 b θ x 2 d x d t 1 2 Q T θ 2 + b v x 2 d x d t 1 2 ] C ( 1 + Z 3 8 ) Y 2 2 b + 3 δ + sup Q T θ b + 1 Q T θ 2 + b v x 2 d x d t 1 2 C ( 1 + Z 3 8 ) [ Y 2 2 b + 3 δ + sup Q T θ b + 1 ( 0 T sup R | θ 1 | b + 1 R θ v x 2 d x d t + + Q T θ v x 2 d x d t ) 1 2 ] C ( 1 + Z 3 8 ) Y 2 2 b + 3 δ + sup Q T θ b + 1 C sup Q T θ R v x 2 d x + + Y 1 2 b + 3 δ 1 2 C ( 1 + Z 3 8 ) C + C Y b + 2 2 b + 3 δ C + C Z 3 ( 2 b + 3 δ ) 8 ( b + 1 δ ) + ϵ Y C + C Z 21 3 δ 24 8 δ + ϵ Y C + ε Z + ε Y
Here, we use (15), (47), (48), (51), (54), and Young’s inequality.
Regarding I 22 ,
I 22 Q T ( 1 + θ 2 b ) | θ x θ t v x | d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + C sup Q T ( 1 + θ b ) 0 T sup R ( θ x 2 + θ 2 b θ x 2 ) d t R v x 2 d x ε Q T ( 1 + θ b ) θ t 2 d x d t + C sup Q T ( 1 + θ b ) 0 T sup R ( θ x 2 + θ 2 b θ x 2 ) d t Y 1 2 b + 3 δ ε Q T ( 1 + θ b ) θ t 2 d x d t + C ( 1 + Y b + 1 2 b + 3 δ ) 0 T sup R ( 1 + θ 2 b ) θ x 2 d t .
For the last term in the right-hand side of (65), we have the following estimate:
0 T sup R ( 1 + θ 2 b ) θ x 2 d t Q T ( 1 + θ 2 b ) θ x 2 d x d t + C Q T ( 1 + θ b ) θ x v | ( 1 + θ b ) θ x v x | d x d t C + C sup Q T θ b + 2 + C ( 1 + sup Q T θ 2 ) Q T θ b 2 θ x 2 d x d t 1 2 × Q T ( 1 + θ b ) ( θ t 2 + u x 4 + θ 2 u x 2 ) d x d t 1 2 C + C Y b + 2 2 b + 3 δ + C ( 1 + Y 1 2 b + 3 δ ) ( Q T ( 1 + θ b ) θ t 2 d x d t + Q T ( 1 + θ b ) u x 4 d x d t + ( θ 2 + sup Q T θ b + 2 ) Q T u x 2 d x d t ) 1 2 C + C Y b + 2 2 b + 3 δ + Y b + 4 2 ( 2 b + 3 δ ) + ϵ Q T ( 1 + θ b ) θ t 2 d x d t + C ( 1 + Y 1 2 b + 3 δ ) Q T ( 1 + θ b ) u x 4 d x d t 1 2 .
Then, using (56) and the Gagliardo–Nirenberg inequality, we infer that
( 1 + sup Q T θ b 2 ) Q T u x 4 d x d t 1 2 C ( 1 + Y b 2 ( 2 b + 3 δ ) ) 0 T u 2 L u x x L 2 2 d t 1 2 C ( 1 + Y b 2 ( 2 b + 3 δ ) ) 0 T R u | u x | d x R u x x 2 d x d t 1 2 C ( 1 + Y b 2 ( 2 b + 3 δ ) ) 0 T R u 2 d x 1 2 R u x 2 d x 1 2 R u x x 2 d x d t 1 2 C ( 1 + Y b 2 ( 2 b + 3 δ ) ) sup 0 t T R u x 2 d x 1 4 0 T R u x x 2 d x d t 1 2 C ( 1 + Y b 2 ( 2 b + 3 δ ) ) Y 3 4 ( 2 b + 3 δ ) Y 3 2 ( 2 b + 3 δ ) C Y 2 b + 9 4 ( 2 b + 3 δ ) + C .
Then, we can obtain that
I 22 C + C Y b + 2 2 b + 3 δ + Y b + 5 2 ( 2 b + 3 δ ) + C Y 6 b + 17 4 ( 2 b + 3 δ ) + ε Q T ( 1 + θ b ) θ t 2 d x d t C + ε Z + ε Y + ε Q T ( 1 + θ b ) θ t 2 d x d t ,
Here, we use the fact that if b > 5 2 ; then, 6 b + 17 4 ( 2 b + 3 δ ) < 1 .
Regarding I 23 ,
I 23 C Q T ( θ + θ 1 + 2 b ) | θ x u x x | d x d t Q T u x x 2 d x d t 1 2 Q T θ 2 θ x 2 d x d t 1 2 + Q T θ 1 + b u x x 2 d x d t 1 2 Q T θ 3 b + 1 θ x 2 d x d t 1 2 C + C Y 7 2 ( 2 b + 3 δ ) + sup Q T θ b + 1 2 Q T u x x 2 d x d t 1 2 sup Q T θ 2 b + 3 2 Q T θ b 2 θ x 2 d x d t 1 2 C + C Y 7 2 ( 2 b + 3 δ ) + C Y 3 b + 7 2 ( 2 b + 3 δ ) C + ε Z + ε Y .
Here, we use the fact that 7 2 ( 2 b + 3 δ ) < 3 b + 7 2 ( 2 b + 3 δ ) < 1 , if b > 1 .
Then, we give the estimate on I 3 ,
I 3 = Q T R u x 2 c v v μ ¯ θ u x c v v ( 1 + θ b ) θ t v θ v t v 2 θ 1 + b v t ( 1 + b ) v 2 d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + C Q T ( ( 1 + θ b ) | u x | 4 + θ u x 3 + θ 1 + b u x 3 + ( 1 + θ 1 + b ) θ u x 2 + θ 2 u x 2 + θ 2 + b u x 2 ) d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + C ( 1 + sup Q T θ b ) Q T u x 4 d x d t + C sup Q T θ b + 2 sup Q T | u x | Q T u x 2 θ d x d t + C ( 1 + sup Q T θ b + 3 ) Q T u x 2 θ d x d t ε Q T ( 1 + θ b ) θ t 2 d x d t + C Y 4 b + 18 4 ( 2 b + 3 δ ) + C Y b + 2 2 b + 3 δ Z 3 8 + C Y b + 3 2 b + 3 δ + C C + ε Z + ε Y + ε Q T ( 1 + θ b ) θ t 2 d x d t .
Here, we use the fact that 4 b + 18 4 ( 2 b + 3 δ ) < 1 , if b > 3 2 and b + 2 2 b + 3 δ < 5 8 , if b > 1 2 .
Adding the estimations of I 1 , I 2 , I 3 and taking a suitable δ , it holds that
Y + Q T ( 1 + θ b ) θ t 2 d x d t C + ε 1 Z .
This completes the proof of Lemma 9. □
Next, we give the uniform boundedness of z.
Lemma 10.
For any b ( 5 2 , ) and ( x , t ) Q T , it holds that
sup [ 0 , T ] R u t 2 d x + Q T u x t 2 d x d t C + ε 2 Z ,
Z C , Y C .
Proof. 
Differentiating (2b) with respect to t, multiplying by u t , and integrating over Q T yields
R u t 2 d x + Q T u x t 2 v d x d t C + Q T μ ¯ u x 2 v 2 + R θ t v θ u x v 2 u x t d x d t + C ε Q T u x t 2 v d x d t + C Q T ( u x 4 + θ t 2 + u x 2 θ 2 ) d x d t + C ε Q T u x t 2 v d x d t + Q T θ t 2 d x d t + C ( sup Q T u x 2 θ + sup Q T θ 3 ) Q T u x 2 θ d x d t + C ε Q T u x t 2 v d x d t + C Z 3 4 Y 1 2 b + 3 δ + C Y 3 2 b + 3 δ + C ε Q T u x t 2 v d x d t + C + ε 2 Z ,
Here, we use (59) and the fact that 1 2 b + 3 δ < 1 4 , if b > 1 2 . This completes the proof of (68).
Rewrite Equation (2b) as
μ ¯ u x x v = u t + R θ v x + μ ¯ u x v x v 2 ,
which implies that
Z C sup 0 t T R u t 2 d x + R u x 2 v x 2 d x + R θ x 2 d x + R θ 2 v x 2 d x C 1 + ε 2 Z + sup Q T ( u x 2 + θ 2 ) R v x 2 d x + R θ x 2 d x C 1 + ε 2 Z + Z 3 4 Y 1 2 b + 3 δ + Y C + ε 3 Z .
Substituting (59) into (71), we obtain
Z C , Y C .
This completes the proof. □
By Lemmas 5–10, we have the following high-order estimates.
Lemma 11.
Via the estimations above, for any b ( 5 2 , ) and ( x , t ) Q T , it holds
sup Q T θ C , sup Q T | u x | C , Q T ( 1 + θ b ) θ t 2 d x d t C ,
sup 0 t T R ( u t 2 + v t 2 ) d x + Q T u x x 2 d x d t + Q T ( u x t 2 + v x t 2 ) d x d t C .
It remains to obtain the L 2 ( R ) -norm bound of θ t , L 2 ( Q T ) -norm bound of θ x t and θ x x .
Lemma 12.
For any b ( 5 2 , ) and ( x , t ) Q T , it holds that
sup 0 t T R θ t 2 d x + Q T θ x t 2 d x d t + Q T θ x x 2 d x d t C .
Proof. 
Differentiating the temperature equation with respect to t, multiplying the result equation by θ t , and integrating over Q T using (73) and (74) gives
d d t R θ t 2 d x + R ( 1 + θ b ) θ x t 2 d x C R | θ x t | ( | θ t | | u x | + | u x | 2 + | u x t | + | θ t | | θ x | + | θ x | | u x | ) d x ε R θ x t 2 d x + C R ( θ x 2 + u x 2 ) θ t 2 d x + C R θ x 2 u x 2 d x + C R ( u x 2 + θ t 2 + u x t 2 ) d x ε R θ x t 2 d x + C ( sup R u x 2 + sup R θ x 2 ) R θ t 2 d x + C R ( θ t 2 + u x t 2 ) d x + C .
Combining this with (66) and (72) , we have
0 T sup R θ x 2 d t C .
Then, applying the Gronwall inequality on (76) yields
sup 0 t T R θ t 2 d x + Q T ( 1 + θ b ) θ x t 2 d x d t C .
Next, we rewrite the temperature equation as follows:
κ ¯ ( 1 + θ b ) θ x x v = θ t + R θ u x v μ ¯ u x 2 v b κ ¯ θ b 1 θ x 2 v + κ ¯ θ b θ x v x v 2 ,
which gives
R ( 1 + θ 2 b ) θ x x 2 d x C R ( θ t 2 + u x 2 + u x 4 + θ x 4 + θ x 2 v x 2 ) d x C + C sup R θ x 2 C + C R | θ x θ x x | d x C + C R θ x x 2 d x 1 2 ,
and implies
sup 0 t T R θ x x 2 d x C .
We complete the proof of Lemma 12. □
It is clear that we carried out all estimates in (12) of Theorem 2. Then, Theorem 2 follows via the standard procedures. We omitted the details.

4. Proof of Asymptotic Behavior

With Lemmas 7–11 in hand, we study the asymptotic behavior as t of the solutions to System (2) with κ = κ ¯ ( 1 + θ b ) for b ( 5 2 , ) in this section. When b = 0 , see Li [16], who deduced the uniform positive lower and upper bounds on temperature θ ( x , t ) with a smart test function.
The following large-time behavior of global solutions together with Lemmas 7–11 finish the proof of Theorem 2.
Lemma 13.
For all t [ 0 , ) and b ( 5 2 , ) , it holds that
C 1 θ ( x , t ) C ,
lim t ( v 1 , u , θ 1 ) ( t ) L p ( R ) + ( v x , u x , θ x ) ( t ) L 2 ( R ) = 0 ,
for any p ( 2 , ] .
Proof. 
Integrating (2b) multiplied by u x x over R leads to
0 | d d t R u x 2 d x | d t + 0 R u x x 2 v d x d t C 0 ( sup R u x 2 ) R v x 2 d x d t + sup Q T θ 0 θ v x 2 d x d t + C 0 R θ x 2 d x d t C 0 sup R u x 2 d t + C sup R × [ 0 , ) θ 0 R θ x 2 θ d x d t + C C .
Next, multiplying the temperature equation by θ x x and integrating over R leads to
c v 2 d d t R θ x 2 d x + κ ¯ R ( 1 + θ b ) θ x x 2 v d x = b κ ¯ R θ b 1 θ x 2 θ x x v d x + κ ¯ R θ b θ x v x θ x x v 2 d x μ ¯ R u x 2 θ x x v d x + R R θ u x θ x x v d x .
Using the Cauchy inequality and Sobolev interpolation inequality gives
0 | R θ b 1 θ x 2 θ x x v d x + R θ b θ x v x θ x x v 2 d x R u x 2 θ x x v d x + R θ u x θ x x v d x | d t C 0 θ x x L 2 θ x L 2 θ x L d t + C 0 θ x x L 2 θ x L v x L 2 d t + C 0 θ x x L 2 u x L u x L 2 d t + C 0 θ x x L 2 u x L 2 d t C 0 θ x x L 2 ( θ x x L 2 1 2 θ x L 2 1 2 + u x x L 2 1 2 u x L 2 1 2 ) d t + C 0 θ x x L 2 u x L 2 d t ε 0 R θ x x 2 d x d t + C .
It follows from (81) and (83) that
0 | d d t u x ( · , t ) L 2 ( R ) 2 | + u x ( · , t ) L 2 ( R ) 2 d t + 0 | d d t θ x ( · , t ) L 2 ( R ) 2 | + θ x ( · , t ) L 2 ( R ) 2 d t C ,
which gives
lim t ( u x ( · , t ) L 2 ( R ) + θ x ( · , t ) L 2 ( R ) ) = 0 .
Thanks to the uniform lower and upper bounds of v ( x , t ) , and upper bounds of θ ( x , t ) , we have
sup 0 t R ( v 1 ) 2 d x + sup 0 t R ( θ 1 ) 2 d x C sup 0 t R ( v ln v 1 ) d x + C sup 0 t R ( θ ln θ 1 ) d x C .
On the other hand, we have the following estimates for all t 0 :
( θ 1 ) ( · , t ) C ( R ) 2 C ( θ 1 ) ( · , t ) L 2 ( R ) θ x ( · , t ) L 2 ( R ) C θ x ( · , t ) L 2 ( R ) ,
and
( v 1 ) ( · , t ) C ( R ) 2 C ( v 1 ) ( · , t ) L 2 ( R ) v x ( · , t ) L 2 ( R ) C v x ( · , t ) L 2 ( R ) .
This, combined with (85), shows
lim t θ ( · , t ) 1 C ( R ) = 0 .
Hence, there exists some T 0 > 0 such that for all ( x , t ) R × [ T 0 , )
1 2 θ ( x , t ) 3 2 .
Lastly, it follows from the proof in Lemma 3 that there exists a constant, such that, for all ( x , t ) R × [ 0 , T 0 ]
θ ( x , t ) C ( T 0 ) .
Hence, we have
C 1 θ ( x , t ) C for all t [ 0 , ) .
Then, combining with Lemma 7, one has
0 | d d t v x ( · , t ) L 2 ( R ) 2 | + v x ( · , t ) L 2 ( R ) 2 d t C ,
it yields that
lim t v x ( · , t ) L 2 ( R ) = 0 .
The pointwise bounds of v ( x , t ) and θ ( x , t ) from below and above independent of time were proven in Lemmas 4 and 13. The asymptotic behavior as t of the solutions was proven in Lemma 13. This completes the proof of Theorem 2. □

5. Conclusions

In this paper, we considered the Cauchy problems to a one-dimensional compressible Navier–Stokes system with temperature-dependent heat conductivity, and general large initial data and far-field conditions. We proved that velocity and temperature are uniformly bounded from below and above in time and space. Further, we proved that the global solution was asymptotically stable as time tended to infinity for b > 5 2 . Our approaches relied upon the maximal principle, and the iteration and energy estimate method. The conclusions in this manuscript are primitive. However, there are limitations to this conclusion because the corresponding results were not obtained when 0 < b 5 2 . We can further study the global well-posedness to System (2) with a viscosity coefficient and heat conduction coefficient that are both dependent on density and temperature. However, there is a great challenge: since we could not obtain an expression for velocity as in (27), the uniform estimates to velocity and temperature were difficult to obtain with the Gronwall inequality. Therefore, there were only some small initial value conclusions in this situation.

Author Contributions

Conceptualization, W.S.; methodology, W.S. and J.Z.; validation, J.Z.; formal analysis, W.S. and J.Z.; resources, W.S.; writing—original draft preparation, W.S.; writing—review and editing, W.S. and J.Z.; visualization, W.S. and J.Z.; funding acquisition, W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (nos. 12161087 and 11801495), the Natural Science Foundation of Jiangxi Province (no. 20212BAB211017), and the Science and Technology Project of Education Department of Jiangxi Province (nos. GJJ211601 and GJJ180833).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Z.L.; Guo, Z.H. On free boundary problem for compressible Navier-Stokes equations with temperature-dependent heat conductivity. Discret. Contin. Dyn. Syst. 2017, 22, 3903–3919. [Google Scholar]
  2. Antontsev, S.N.; Kazhikhov, A.V.; Monakhov, V.N. Boundary Value Problems in Mechanics of Nonhomogeneous Fluids; Studies in Mathematics and its Applications; North-Holland Publishing Co.: Amsterdam, The Netherlands, 1990; Volume 22. [Google Scholar]
  3. Kazhikhov, A.V.; Shelukhin, V.V. Unique global solution with respect to time of initial-boundary value problems for one-dimensional equations of a viscous gas. J. Appl. Math. Mech. 1977, 41, 273–282. [Google Scholar] [CrossRef]
  4. Kazhikhov, A.V. On the Cauchy problem for the equations of a viscous gas. Sib. Math. J. 1982, 23, 60–64. [Google Scholar] [CrossRef]
  5. Amosov, A.A.; Zlotnik, A.A. Solvability “in the large” of a system of equations for the one-dimensional motion of an inhomogeneous viscous heat-conducting gas. Math. Notes 1992, 52, 753–763. [Google Scholar] [CrossRef]
  6. Nagasawa, T. On the one-Dimensional motion of the polytropic ideal gas non-fixed on the boundary. J. Differ. Equ. 1986, 65, 49–67. [Google Scholar] [CrossRef] [Green Version]
  7. Nagasawa, T. On the One-Dimensional Free Boundary Problem for the Heat-Conductive Compressible Viscous Gas; Recent Topics in Nonlinear PDE IV, Lecture Notes in Numerical and Applied Analysis; Kinokuniya/North-Holland: Amsterdam, The Netherlands, 1989; Volume 10, pp. 83–99. [Google Scholar]
  8. Nagasawa, T. Equations of Motion of Compressible Viscous Fluids; Kinokuniya/North-Holland: Amsterdam, The Netherlands, 1986; Volume 18, pp. 97–128. [Google Scholar]
  9. Qin, Y.M. Nonlinear Parabolic-Hyperbolic Coupled Systems and Their Attractors; Operator Theory, Advances and Applications; Birkhäuser: Basel, Switzerland, 2008; Volume 184. [Google Scholar]
  10. Zlotnik, A.A.; Amosov, A.A. Global generalized solutions of the equations of the one-dimensional motion on a viscous heat-conducting gas. Soviet Math. Dokl. 1989, 38, 11–15. [Google Scholar]
  11. David, H. Global well-posedness of the Cauchy problem for the Navier-Stokes equations of nonisentropic flow with discontinuous initial data. J. Differ. Equ. 1992, 95, 33–74. [Google Scholar]
  12. Jiang, S. Large-time behavior of solutions to the equations of a one-dimensional viscous polytropic ideal gas in unbounded domains. Commun. Math. Phys. 1999, 200, 181–193. [Google Scholar] [CrossRef]
  13. Ya, I.K. Cauchy problem for the equations of gasdynamics with viscosity. Sib. Math. J. 1979, 20, 208–218. [Google Scholar]
  14. Kawashima, S. Large-time behaviour of solutions to hyperbolic–parabolic systems of conservation laws and applications. Proc. R. Soc. Edinb. Sect. A Math. 1987, 106, 169–194. [Google Scholar] [CrossRef]
  15. Kazhikhov, A.V. Cauchy problem for viscous gas equations. Sib. Math. J. 1982, 23, 44–49. [Google Scholar] [CrossRef]
  16. Li, J.; Liang, Z.L. Some uniform estimates and large-time behavior of solutions to one-dimensional compressible Navier-Stokes system in unbounded domains with large data. Arch. Ration. Mech. Anal. 2016, 200, 1195–1208. [Google Scholar] [CrossRef] [Green Version]
  17. Chapman, S.; Cowling, T.G. The Mathematical Theory of Non-Uniform Gases: An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion in Gases; Cambridge University Press: Cambridge, UK, 1970. [Google Scholar]
  18. Vincenti, W.G.; Kruger, C.H. Introduction to Physical Gas Dynamics; John Wiley and Sons: New York, NY, USA, 1965; p. 242. [Google Scholar]
  19. Liu, H.X.; Yang, T.; Zhao, H.J.; Zou, Q.Y. One-dimensional compressible Navier-Stokes equations with temperature dependent transport coefficients and large data. SIAM J. Math. Anal. 2014, 46, 2185–2228. [Google Scholar] [CrossRef]
  20. Wang, T.; Zhao, H.J. Global large solutions to a viscous heat-conducting one-dimensional gas with temperature-dependent viscosity. arXiv 2015, arXiv:1505.05252. [Google Scholar]
  21. Wang, T.; Zhao, H.J. One-dimeisional compressible heat-conducting gas with temperature-dependent viscosity. Math. Model. Methods Appl. Sci. 2016, 26, 2237–2275. [Google Scholar] [CrossRef]
  22. Duan, R.; Guo, A.; Zhu, C. Global strong solution to compressible Navier–Stokes equations with density dependent viscosity and temperature dependent heat conductivity. J. Differ. Equ. 2017, 262, 4314–4335. [Google Scholar] [CrossRef]
  23. Li, K.X.; Shu, X.L.; Xu, X.J. Global Existence of Strong Solutions to Compressible Navier-Stokes System with Degenerate Heat Conductivity in Unbounded Domains. Math. Methods Appl. Sci. 2020, 43, 1543–1554. [Google Scholar] [CrossRef] [Green Version]
  24. Jenssen, H.K.; Karper, T.K. One-dimensional compressible flow with temperature dependent transport coefficients. SIAM J. Math. Anal. 2010, 42, 904–930. [Google Scholar] [CrossRef]
  25. Pan, R.H.; Zhang, W.Z. Compressible Navier-Stokes equations with temperature dependent heat conductivity. Commun. Math. Sci. 2015, 13, 401–425. [Google Scholar] [CrossRef]
  26. Chen, T.; Zhang, Y. Free boundary problem for one-dimensional compressible Navier–Stokes equations with temperature-dependent viscosity and heat conductivity. Math. Methods Appl. Sci. 2021, 44, 13273–13286. [Google Scholar] [CrossRef]
  27. Cai, G.; Chen, Y.; Peng, Y.F.; Peng, Y. On the asymptotic behavior of the one-dimensional motion of the polytropic ideal gas with degenerate heat conductivity. J. Differ. Equ. 2022, 317, 225–263. [Google Scholar] [CrossRef]
  28. Ali, N.; Asghar, Z.; Sajid, M.; Anwar Bég, O. Biological interactions between Carreau fluid and microswimmers in a complex wavy canal with MHD effects. J. Braz. Soc. Mech. Sci. 2019, 41, 1–13. [Google Scholar]
  29. Asghar, Z.; Ali, N.; Javid, K.; Waqas, M.; Dogonchi, A.S.; Khan, W.A. Bio-inspired propulsion of micro-swimmers within a passive cervix filled with couple stress mucus. Comp. Meth. Prog. Biom. 2020, 189, 105313. [Google Scholar] [CrossRef]
  30. Asghar, Z.; Ali, N.; Waqas, M.; Nazeer, M.; Khan, W.A. Locomotion of an efficient biomechanical sperm through viscoelastic medium. Biomech. Model. Mechanobiol. 2020, 19, 2271–2284. [Google Scholar] [CrossRef]
  31. Domnich, A.A.; Baranovskii, E.S.; Artemov, M.A. A nonlinear model of the non-isothermal slip flow between two parallel plates. J. Phys. Conf. Ser. 2020, 1479, 012005. [Google Scholar] [CrossRef]
  32. Tani, A. On the first initial-boundary value problem of compressible viscous fluid motion. Publ. Res. Inst. Math. Sci. 1977, 13, 193–253. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Su, W.; Zhong, J. Asymptotic Behavior of the Solution to Compressible Navier–Stokes System with Temperature-Dependent Heat Conductivity in an Unbounded Domain. Symmetry 2023, 15, 112. https://doi.org/10.3390/sym15010112

AMA Style

Su W, Zhong J. Asymptotic Behavior of the Solution to Compressible Navier–Stokes System with Temperature-Dependent Heat Conductivity in an Unbounded Domain. Symmetry. 2023; 15(1):112. https://doi.org/10.3390/sym15010112

Chicago/Turabian Style

Su, Wenhuo, and Jianxin Zhong. 2023. "Asymptotic Behavior of the Solution to Compressible Navier–Stokes System with Temperature-Dependent Heat Conductivity in an Unbounded Domain" Symmetry 15, no. 1: 112. https://doi.org/10.3390/sym15010112

APA Style

Su, W., & Zhong, J. (2023). Asymptotic Behavior of the Solution to Compressible Navier–Stokes System with Temperature-Dependent Heat Conductivity in an Unbounded Domain. Symmetry, 15(1), 112. https://doi.org/10.3390/sym15010112

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop