Home Stability of combination of rarefaction waves with viscous contact wave for compressible Navier-Stokes equations with temperature-dependent transport coefficients and large data
Article Open Access

Stability of combination of rarefaction waves with viscous contact wave for compressible Navier-Stokes equations with temperature-dependent transport coefficients and large data

  • Wenchao Dong and Zhenhua Guo EMAIL logo
Published/Copyright: August 26, 2022

Abstract

In this article, we study the large-time behavior of combination of the rarefaction waves with viscous contact wave for a one-dimensional compressible Navier-Stokes system whose transport coefficients depend on the temperature. It is shown that if the adiabatic exponent γ is suitably close to 1, the unique solution global in time to ideal polytropic gas exists and asymptotically tends toward the combination of a viscous contact wave with rarefaction waves under large initial perturbation. New and subtle analysis is developed to overcome difficulties due to the smallness of γ – 1 to derive heat kernel estimates. Moreover, our results extend the studies in a previous work [F. M. Huang, J. Li, and A. Matsumura, Arch. Ration. Mech. Anal. 197 (2010), no. 1, 89–116].

MSC 2010: 35B35; 35B40; 76N06; 76N30

1 Introduction

This article is concerned with the large-time behavior of global smooth solution to the Cauchy problem of the following one-dimensional compressible Navier-Stokes equations in the Lagrangian coordinates

(1.1) v t u x = 0 , u t + p x = μ u x v x , e + 1 2 u 2 t + ( p u ) x = μ u u x v + κ θ x v x ,

supplemented with the initial data and far field conditions:

(1.2) ( v , u , θ ) ( x , 0 ) = ( v 0 , u 0 , θ 0 ) ( x ) , ( v , u , θ ) ( ± , t ) = ( v ± , u ± , θ ± ) ,

where x R , t > 0 , v ± > 0 , θ ± > 0 , and u ± R are given constants. v ( x , t ) > 0 , u ( x , t ) , e ( x , t ) > 0 , θ ( x , t ) > 0 , and p ( x , t ) are specific volume, fluid velocity, specific internal energy, absolute temperature and pressure, respectively. The viscosity coefficient μ and heat conductivity coefficient κ are prescribed through constitutive relations as smooth functions of θ which are assumed to satisfy

(1.3) μ = μ ( θ ) > 0 , κ = κ ( θ ) > 0 , θ > 0 .

In this article, we focus on the ideal polytropic gas in which constitutive relations read

e = c v θ = R γ 1 θ , p = R θ v = A v γ e s c v ,

where c v , R , A are positive constants, γ > 1 is the adiabatic constant, and s is the entropy. The study on the global solvability and the large-time behaviors of the corresponding global solution to (1.1) and (1.2) is one of the hottest topics in the field of nonlinear partial differential equations and many results have been obtained.

If μ , κ 0 , then system (1.1) becomes the classical compressible Euler system, which is one of the most important nonlinear strictly hyperbolic systems of conservation laws. It is well known that the Euler system has three basic wave patterns: a linearly degenerate wave (contact discontinuity) and two nonlinear waves (rarefaction wave or shock wave) [35]. The above three dilation invariant solutions and their linear superpositions, called Riemann solutions, govern both local and large-time behavior of solutions to the Euler system. Since the Euler system is an idealization if we neglect the dissipative effect, it is very important to study the large-time asymptotic behavior of solutions to (1.1), toward the viscous versions of these basic waves. Indeed, there have been many results on the global solvability and the precise descriptions of the large-time behaviors of the corresponding global solutions to the Cauchy problem (1.1) and (1.2) since 1985 [23,30].

For the case μ , κ are positive constants, the results are quite complete for (1.1) and (1.2). There are many meaningful studies on stability analysis of the viscous wave patterns, see [7,8,10,12,13,14, 15,17,19,20,23,27,28,30,34,39,41] and references cited therein. Particularly, in different ways, the nonlinear stability of viscous shock waves has been established in [8,13,15,23,30], etc. In the case where the Riemann solution to the Euler system consists of contact discontinuity, the stability results were obtained in [12,14,17,20,28,39]. In particular, the asymptotics toward the rarefaction waves for (1.1) and (1.2) is given by Kawashima et al. in [24] according to an energy form associated with the physical total energy and a monotone property of rarefaction waves, and more in [7,12,17,27,34].

When the transport coefficients depend on temperature, the only results available now are obtained in [6,11,16, 18,29,38]. In 2014, Liu et al. [29] primely obtained that nonvacuum constants for (1.1) and (1.3) are asymptotically stable with large perturbation, given the condition that the adiabatic exponent γ is close enough to 1. When μ = c 1 h ( v ) θ b , κ = c 2 h ( v ) θ b , b 1 , Wang and Zhao [38] and Huang and Liao [11] studied the nonlinear stability of nonvacuum constants and the combination of viscous contact wave with rarefaction waves, respectively, under large perturbation. But the result of [11] excludes case (1.4). In [16], Huang et al. discussed the stability of viscous shock profiles of the Cauchy problem (1.1) and (1.2) while μ = v k μ ( θ ) , κ = v k κ ( θ ) , 0 < k < 1 2 provided γ ( 1 , 2 ] . However, the results in [11,16,29,38] where not only δ = |θ{+}{-}|\ll1 but also the initial perturbation belongs to H 3 is needed.

Our study of the dependence of μ , κ on temperature is motivated by the following consideration. It is well known that, if the intermolecule potential is proportional to r α with r being the molecule distance and α > 0 , we can derive the compressible Navier-Stokes equations (1.1) from the Boltzmann equation with slab symmetry by employing the celebrated Chapman-Enskog expansion [5,4,9, 29], then the corresponding transport coefficients μ , κ depend on temperature, satisfy

(1.4) μ = μ ¯ θ b , κ = κ ¯ θ b , b 1 2 , + ,

where μ ¯ > 0 , κ ¯ > 0 , b = α + 4 2 α are constants. Note that for Maxwellian molecules b = 1 , while for elastic spheres b = 1 2 , see [5]. The above dependence has a strong influence on the solution behavior and leads to difficulty in analysis for global existence with large data.

Recently, Dong and Guo reported that the nonlinear stability of viscous contact wave for μ , κ satisfies (1.2) in [6], and they have obtained not only a Nishida-Smaller-type large data result, but also a decay rate ( 1 + t ) 3 8 + C δ 1 4 for the L -norm of perturbation and its derivative. From all the aforementioned studies, a natural question is to see whether similar global nonlinear stability results of the rarefaction waves corresponding to (1.1)–(1.3) hold or not. The main purpose of this article is devoted to this problem. We now focus on the case when the large-time behavior of (1.1) and (1.2) is described by the superposition of 1-rarefaction wave, viscous contact wave, and 3-rarefaction wave, then the key point to this problem with large data is to obtain the positive lower and upper bounds for v ( x , t ) , θ ( x , t ) . When the viscosity μ depends on the temperature θ , the identity corresponding to momentum equation (1.1) 2 becomes

μ ( θ ) v x v t = u t + p x + μ ( θ ) v ( θ t v x θ x u x ) .

The temperature dependence of the viscosity μ ( μ ( θ ) 0 ) has a strong influence on the solution and leads to difficulty in mathematical analysis for global solvability with large data. Motivated by [33,29], since the energy equation (1.1) 3 gives

θ t = γ 1 R p u x + μ ( θ ) u x 2 v + κ ( θ ) θ x v x ,

it is possible to use “smallness of γ 1 ” to control the nonlinear term and to yield some estimates on the positive lower and upper bounds for v , θ . Moreover, the stability of the superposition of several wave patterns is more complicated and challenging because of the fact that the stability analysis essentially depends on the underlying properties of a basic wave pattern, and these frameworks are not compatible with each other. In addition, the wave interaction between different families of wave patterns is complicated due to μ = μ ( θ ) , κ = κ ( θ ) . In order to overcome these difficulties, our idea is to use the smallness of γ 1 to obtain the energy estimates and then to obtain the stability of waves.

In the present article, we establish asymptotic stability of combination of the rarefaction waves with viscous contact wave for (1.1). According to [12], we give a new elementary inequality concerning the estimate of the term 0 t ( γ 1 ) 1 ( 1 + s ) 1 e c x 2 ( γ 1 ) ( 1 + s ) ( ϕ , ψ , ζ ) 2 d x d s , which can be extended to the general gas (Lemma 4.3). Then we derive a Nishida-Smoller-type [33] result for the Cauchy problem (1.1)–(1.3). Indeed, it is shown that the superposition of 1-rarefaction wave, viscous contact wave and 3-rarefaction wave is stable with large initial perturbation provided γ 1 is suitably small. Particularly, in [11,16,18,29,38], their method requires the initial perturbation ϕ 0 , ψ 0 , ζ 0 γ 1 H 3 ( R ) . We improve their results through this article, only asking that ϕ 0 , ψ 0 , ζ 0 γ 1 in H 2 ( R ) (see Theorem 2.1 below for the details).

This article is organized as follows. In Section 2, we recall the definition of viscous contact wave and rarefaction wave for (1.1), then state our main result. Section 3 gives four lemmas which will be used later. Finally, Section 4 is devoted to giving the a priori estimates and complete the proof of our main result (Theorem 2.1).

Notations: In the following, O ( 1 ) , c , or C from line to line denote the generic positive constants which are independent of time t unless otherwise stated. And for two functions f ( x ) , g ( x ) , f ( x ) g ( x ) , as x a means that there exists a generic positive constant C , which is independent of x such that C 1 g ( x ) f ( x ) C g ( x ) in a neighborhood of a . For function space, L p ( Ω ) ( 1 p ) represents the usual Lebesgue space on Ω with the norm L p ( Ω ) :

u L p ( Ω ) = Ω u ( x ) p d x 1 p , 1 p < , u L ( Ω ) = ess sup Ω u ( x ) .

W k , p ( Ω ) ( k 0 ) denotes the kth order Sobolev space with its norm

u W k , p ( Ω ) = j k d j u ( x ) d x j L p ( Ω ) p 1 p , 1 p < , j k d j u ( x ) d x j L ( Ω ) , p = .

Particularly, for notational simplicity, we denote H k W k , 2 and let

= L 2 ( R ) , H k = H k ( R ) , L p = L p ( R ) , d x = R d x .

2 Main result

In this section, we first revisit the definition of viscous contact wave and rarefaction wave, then state our main result.

We are interested in the global solutions at time of the Cauchy problem (1.1)–(1.3) and their large-time behaviors in the relations with the spatial asymptotic states ( v ± , u ± , θ ± ) . It is well known that these asymptotic behaviors are well characterized by those of the solutions of the corresponding Riemann problem for the hyperbolic part of (1.1) (Euler system):

(2.1) v t u x = 0 , u t + p x = 0 , e + 1 2 u 2 t + ( p u ) x = 0 ,

with the initial Riemann data

(2.2) ( v , u , θ ) ( x , 0 ) = ( v , u , θ ) , x < 0 , ( v + , u + , θ + ) , x > 0 .

The basic theory of hyperbolic systems of conservation laws [35] implies that, for any given constant state ( v , u , θ ) , there exists a suitable neighborhood Ω ( v , u , θ ) of ( v , u , θ ) such that, for any ( v + , u + , θ + ) Ω ( v , u , θ ) , the Riemann problem (2.1) and (2.2) has a unique solution. In this article, we only consider the case of the superposition of the rarefaction waves and viscous contact wave with

(2.3) ( v + , u + , θ + ) R 1 C R 3 ( v , u , θ ) Ω ( v , u , θ ) ,

where

R 1 C R 3 ( v , u , θ ) ( v , u , θ ) Ω ( v , u , θ ) s s , u u v e s s c v γ v λ ( η , s ) d η , u u e s s c v γ v v λ + ( η , s ) d η ,

s = c v ln R θ A + R ln v , s ± = c v ln R θ ± A + R ln v ± , λ ± ( v , s ) = ± A γ v γ 1 e s c v .

By the standard argument, there exists a unique pair of points ( v ± m , u ± m , θ ± m ) Ω ( v , u , θ ) that satisfy

(2.4) p m R θ m v m = R θ + m v + m , u m = u + m ,

such that ( v ± m , u ± m , θ ± m ) are connected by a contact discontinuity and the points ( v m , u m , θ m ) and ( v + m , u + m , θ + m ) belong to the 1-rarefaction wave curve R ( v , u , θ ) and the 3-rarefaction wave curve R + ( v + , u + , θ + ) , respectively, where

(2.5) R ± ( v ± , u ± , θ ± ) ( v , u , θ ) Ω ( v , u , θ ) s = s ± , u = u ± v ± v λ ± ( η , s ± ) d η , v > v ± .

Without loss of generality, we may assume u ± m = 0 in what follows. In the setting of the compressible Navier-Stokes equations (1.1), the corresponding wave ( V c , U c , Θ c ) to the contact discontinuity becomes smooth and behaves as a diffusion wave due to the dissipation effect. We call this wave a viscous contact wave. The viscous contact wave ( V c , U c , Θ c ) can be constructed as follows. Since the pressure for the profile ( V c , U c , Θ c ) is expected to be almost constant, that is, we set

R Θ c V c = p m ,

which indicates that the leading part of the energy equation (1.1) 3 is

c v Θ t c + p m U x c = κ ( Θ c ) Θ x c V c x .

The last equation and (1.1) 1 lead to a nonlinear diffusion equation

(2.6) Θ t c = ( a ( Θ c ) Θ x c ) x , Θ c ( ± , t ) = θ ± m , a ( Θ c ) = ( γ 1 ) p + κ ( Θ c ) γ R 2 Θ c > 0 ,

which has a unique self-similarity solution Θ c ( x , t ) = Θ c ( ξ ) , ξ = x ( γ 1 ) ( 1 + t ) due to [3] (Theorem 6), [2,37]. Furthermore, on the one hand, Θ c ( ξ ) is a monotone function, increasing if θ + m > θ m and decreasing if θ + m < θ m ; on the other hand, there exists positive constant δ ¯ , such that if δ c = θ + m θ m < δ ¯ , Θ c satisfies

(2.7) ξ k ( Θ c θ ± m ) c 1 δ c e c 2 ξ 2 , k 0 ,

where c 1 , c 2 are positive constants depending only on θ ± m , but independently of γ 1 . Once Θ c is determined, we define V c and U c by

(2.8) V c = R Θ c p m , U c = R a ( Θ c ) p m Θ x c .

Then the 1-rarefaction wave ( v r , u r , θ r ) x t (respectively, the 3-rarefaction wave ( v + r , u + r , θ + r ) x t ) connecting ( v , u , θ ) and ( v m , 0 , θ m ) (respectively, connecting ( v + m , 0 , θ + m ) and ( v + , u + , θ + ) ) is a weak solution of the Riemann problem of Euler system (2.1) with the following initial Riemann data:

(2.9) ( v ± r , u ± r , θ ± r ) ( x , 0 ) = ( v ± m , 0 , θ ± m ) , ± x < 0 , ( v ± , u ± , θ ± ) , ± x > 0 .

Since the rarefaction wave represents weak solutions, it is convenient to construct an approximate rarefaction wave that is smooth. Motivated by [17,31], the smooth solutions of Euler system (2.1), ( V ± r , U ± r , Θ ± r ) ( x , t ) , are defined by

(2.10) λ ± ( V ± r ( x , t ) , s ± ) = w ± ( x , t ) , U ± r ( x , t ) = u ± v ± V ± r ( x , t ) λ ± ( η , s ± ) d η , Θ ± r ( x , t ) = θ ± ( v ± ) γ 1 ( V ± r ( x , t ) ) 1 γ ,

where w ( x , t ) (respectively, w + ( x , t ) ) is the unique smooth solution of the initial problem for Burgers equation:

w t + w w x = 0 , w ( x , 0 ) = w r + w l 2 + w r w l 2 tanh x , ( x , t ) R × ( 0 , + ) ,

with w l = λ ( v , s ) , w r = λ ( v m , s ) (respectively, w l = λ + ( v + m , s + ) , w r = λ + ( v + , s + ) ).

To describe the strengths of the rarefaction wave and viscous contact wave, we let

(2.11) δ ± = v ± v ± m + u ± 0 + θ ± θ ± m , δ c = θ + m θ m , δ = max { δ ± , δ c } ,

then define

(2.12) V U Θ = V r + V c + V + r U r + U c + U + r Θ r + Θ c + Θ + r v m + v + m 0 θ m + θ + m .

More importantly, the wave ( V , U , Θ ) solves the compressible Navier-Stokes equations (1.1) time asymptotically, that is,

(2.13) V t U x = 0 , U t + P x = μ ( Θ ) U x V x + R ˜ 1 , c v Θ t + P U x = μ ( Θ ) U x 2 V + κ ( Θ ) Θ x V x + R ˜ 2 ,

where

R ˜ 1 = μ ( Θ ) U x V x + ( P P r P + r ) x + U t c , R ˜ 2 = μ ( Θ ) U x 2 V κ ( Θ ) Θ x V x + κ ( Θ c ) Θ x c V c x + ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x , P = R Θ V , P ± r = R Θ ± r V ± r .

We are now in a position to state our main result. First of all, let

(2.14) ( ϕ , ψ , ζ ) ( x , t ) = ( v V , u U , θ Θ ) ( x , t ) ,

and for any interval I [ 0 , + ) , we define a function space X ( I ) as

X ( I ) = { ( ϕ , ψ , ζ ) C ( I ; H 0 2 ) ϕ x L 2 ( I ; H 1 ) , ( ψ x , ζ x ) L 2 ( I ; H 2 ) } .

Then our main result is as follows:

Theorem 2.1

For any given ( v , u , θ ) , suppose that ( v + , u + , θ + ) satisfies (2.3), μ , κ satisfy (1.3), let ( V , U , Θ ) be defined in (2.12), assume that θ + m θ m C ( γ 1 ) α , θ ± θ ± m C ( γ 1 ) β for some constants α > 5 4 , β > 0 , and there exists ( γ 1 ) -independent positive constants v ̲ , Θ ̲ , Θ ¯ , M 0 such that

(2.15) v 0 ( x ) v ̲ , 2 Θ ̲ θ 0 ( x ) 1 2 Θ ¯ , ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 M 0 .

Then, there exists a small positive constant δ 0 such that if

(2.16) γ 1 δ 0 ,

then the Cauchy problem (1.1)–(1.3) admits a unique global solution ( v , u , θ ) satisfying ( ϕ , ψ , ζ ) X ( [ 0 , + ) ) and

(2.17) lim t + sup x R ( v V , u U , θ Θ ) ( x , t ) = 0 .

Remark 1

The difference θ + m θ m , θ ± θ ± m are naturally bounded by C ( γ 1 ) from the physical point of view (see [10,21]). If μ , κ are constants, we only need that α > 1 2 . And we can check from (2.5) and (2.16) that

u ± 0 = v ± v ± m λ ± ( η , s ± ) d η C v ± v ± m .

Moreover, without loss of generality, we assume that θ ± θ ± m since θ ± θ ± m can be discussed in the same way. We deduce from the mean value theorem that

v ± v ± m = R θ ± A 1 1 γ R θ ± m A 1 1 γ e s ± R = 1 γ 1 η γ 1 γ R A θ ± θ ± m e 1 γ 1 ln R θ ± A + ln v ± = θ ± θ ± m γ 1 R A η 1 v ± e 1 γ 1 ln R θ ± A η C e 1 γ 1 ln R θ ± A η C ( γ 1 ) β ,

where R θ ± A < η < R θ ± m A . Thus, from (2.11), we have δ ± C ( γ 1 ) β C δ 0 β 1 .

Remark 2

Theorem 2.1 means that if γ 1 is smaller, the initial perturbation can be larger, it is a generalization of the results in [17, 18,38]. Particularly, all the results obtained in [18, 38] postulate that the initial perturbation ( ϕ 0 , ψ 0 , ζ 0 ) H 3 ( R ) , but we require only that ( ϕ 0 , ψ 0 , ζ 0 ) H 2 ( R ) .

Remark 3

We can easily observe that (1.4) satisfies (1.3), for each b R . Therefore, our result in Theorem 2.1 holds for (1.4) with any b R , as in [29], can cover the case when the transport coefficients are derived from the Boltzmann equation.

Remark 4

In 2021, Sun et al. [36] studied an initial and boundary value problem of (1.1) for μ = θ α , κ = θ β , α , β 0 . And the nonlinearly exponential stability of the solution is obtained for the large data provided α 1 . In fact, we can also obtain the similar result in Theorem 2.1 for α 1 , β 0 with large perturbation. Note that our derivation of the uniform bounds on v , θ relies heavily on the assumption that γ 1 or α sufficiently small. Thus, it is an interesting and difficult problem to extend the result to the case with transport coefficients μ , κ satisfying μ = θ α , κ = θ β for α > 0 . One would conjecture that the result in Theorem 2.1 should be improved to 0 < α < c for some fixed constant c . However, this will be left for future.

Remark 5

Experimental results for gases at very high temperature show that both the viscosity μ and the heat-conductivity κ may depend on density and temperature (see [26,40]). Even for this case, similar result in Theorem 2.1 still holds if μ ( v , θ ) and κ ( v , θ ) satisfy

(2.18) μ = μ ( v , θ ) > 0 , κ = κ ( v , θ ) > 0 , v > 0 , θ > 0 ,

(2.19) limsup v 0 + , + μ ˜ v V μ ˜ ( v ) < + ,

(2.20) lim v 0 + ( v ) = , lim v + ( v ) = + ,

where μ ˜ ( v ) = min Θ ̲ θ Θ ¯ μ ( v , θ ) , ( v ) = 1 v s ln s 1 s μ ˜ ( s ) d s . Particularly, (2.18)–(2.20) involve the following cases:

μ ˜ ( v ) v a , v 0 + , v b , v + , for some a 0 , b 1 2 .

Particularly, the transport coefficients μ , κ satisfying (2.18)–(2.20) can include μ = c 1 h ( v ) θ b , κ = c 2 h ( v ) θ b in [11] as a special example. In [38], they consider μ = c 1 h ( v ) θ b , κ = c 2 h ( v ) θ b , h ( v ) c ( v l 1 + v l 2 ) , b 1 , but their result cannot cover the model satisfying (1.4) since the parameters l 1 and l 2 are assumed to be positive. Note that our result in Theorem 2.1 gives a positive answer of the problem in [38] for l 1 = l 2 = 0 and any b R if γ 1 1 .

3 Preliminaries

In this section, we give some basic properties of the viscous contact wave, rarefaction wave, the heat kernel estimate, and Gronwall’s inequality, which will be used later, whose proofs can be found in [1,3,11,12,17,37].

Let us first consider a lemma for the properties of the viscous contact wave ( V c , U c , Θ c ) for (1.1) defined in (2.6) and (2.8).

Lemma 3.1

Assume that δ c = θ + m θ m δ ¯ for a small positive constant δ ¯ . Then the viscous contact wave ( V c , U c , Θ c ) has the following properties:

x k V c + x k Θ c c 1 δ c ( γ 1 ) k 2 ( 1 + t ) k 2 e c 2 ξ 2 , k 1 , x k U c c 1 δ c ( γ 1 ) k 1 2 ( 1 + t ) k + 1 2 e c 2 ξ 2 , k 1 , t k 1 x k 2 Θ c c 1 δ c ( γ 1 ) k 2 2 ( 1 + t ) 2 k 1 + k 2 2 e c 2 ξ 2 , k 1 1 , k 2 0 , t k 1 x k 2 U c c 1 δ c ( γ 1 ) k 2 1 2 ( 1 + t ) 2 k 1 + k 2 + 1 2 e c 2 ξ 2 , k 1 1 , k 2 0 ,

where ξ = x ( γ 1 ) ( 1 + t ) and c 1 , c 2 are positive constants depending only on θ ± m , but independently of γ 1 .

We then divide R × ( 0 , t ) into three parts, i.e., R × ( 0 , t ) = Ω Ω c Ω + with

(3.1) Ω = { ( x , t ) 2 x < λ ( v m , s ) t } , Ω c = { ( x , t ) λ ( v m , s ) t 2 x λ + ( v + m , s + ) t } , Ω + = { ( x , t ) 2 x > λ + ( v + m , s + ) t } .

We have the following lemma obtained in [11,17].

Lemma 3.2

For any given ( v , u , θ ) , suppose that ( v + , u + , θ + ) satisfies (2.3). Then the viscous contact wave constructed in (2.6) and (2.8) and the smooth rarefaction waves constructed in (2.10) satisfying the following properties:

  1. ( U ± r ) x 0 .

  2. For any p [ 1 , + ] , there exists a positive constant C = C ( p , v , u , θ ) such that for small δ ,

    ( ( V ± r ) x , ( U ± r ) x , ( Θ ± r ) x ) L p C min δ , δ 1 p t 1 + 1 p , ( ( V ± r ) x x , ( U ± r ) x x , ( Θ ± r ) x x ) L p C min { δ , t 1 } , ( ( V ± r ) x x x , ( U ± r ) x x x , ( Θ ± r ) x x x ) L p C min { δ , δ 1 + 1 p t 2 + 1 p } .

  3. There exists a positive constant C = C ( v , u , θ ) such that for

    c 0 = 1 10 min { λ ( v m , s ) , λ + ( v + m , s + ) , c 2 λ 2 ( v m , s ) , c 2 λ + 2 ( v + m , s + ) } ,

    k 0 , we have

    x k ( V ± r v ± m ) + x k ( U ± r ) + x k ( Θ ± r θ ± m ) C δ e c 0 ( x + t ) , in Ω c , x k ( V c v ± m ) + x k + 1 U c + x k ( Θ c θ ± m ) C δ c ( γ 1 ) k 2 e c 0 ( γ 1 ) 1 ( x + t ) ( V r ) x + ( U r ) x + ( Θ r ) x + V r v m + Θ r θ m C δ e c 0 ( x + t ) , in Ω ± .

  4. For the rarefaction wave ( v ± r , u ± r , θ ± r ) x t determined by (2.1) and (2.9), it holds that

    lim t + sup x R ( V ± r , U ± r , Θ ± r ) ( x , t ) ( v ± r , u ± r , θ ± r ) x t = 0 .

Now, we derive an elementary inequality derived in Lemma 1 in [12].

Lemma 3.3

For 0 < T + , suppose that h ( x , t ) satisfies

h L ( 0 , T ; L 2 ) , h x L 2 ( 0 , T ; L 2 ) , h t L ( 0 , T ; H 1 ) .

Then the following estimate holds

(3.2) 0 T h 2 w 2 d x d t 4 π γ 1 h ( , 0 ) 2 + 2 π c 3 0 T h x ( , t ) 2 d t + 8 c 3 γ 1 0 T h t , h g 2 H 1 × H 1 d t ,

where for c 3 > 0 ,

(3.3) w ( x , t ) = ( γ 1 ) 1 2 ( 1 + t ) 1 2 e c 3 x 2 ( γ 1 ) ( 1 + t ) , g ( x , t ) = x w ( y , t ) d y ,

and , denotes the dual product between H 1 and H 1 .

It is easy to check that

(3.4) 4 c 3 g t = ( γ 1 ) w x , g ( , t ) L = π c 3 1 2 .

We finally give an elementary inequality derived from Lemma 6.1 in [1], which will play an essential role in the second-order energy estimate of ( ϕ , ψ , ζ ) .

Lemma 3.4

Suppose that f 1 ( t ) , f 2 ( t ) : [ 0 , T ] R are nonnegative bounded measurable function, c ( t ) : [ 0 , T ] R is a nonnegative integrable function. If f 1 ( t ) , f 2 ( t ) , c ( t ) satisfy

f 1 ( t ) f 2 ( t ) + 0 t c ( s ) f 1 ( s ) d s , t [ 0 , T ] ,

then it follows that

f 1 ( t ) f 2 ( t ) + 0 t f 2 ( s ) c ( s ) exp s t c ( τ ) d τ d s , t [ 0 , T ] .

Moreover, if f 2 ( t ) is a monotone increasing function over [ 0 , T ] , then we obtain the estimate

f 1 ( t ) f 2 ( t ) exp 0 t c ( s ) d s , t [ 0 , T ] .

4 Proof of Theorem 2.1

Now, we begin to prove Theorem 2.1. From equations (2.13), we rewrite the Cauchy problem (1.1) and (1.2) as

(4.1) ϕ t ψ x = 0 , ψ t + ( p P ) x = μ ( θ ) u x v U x V x + R 1 , c v ζ t + p u x P U x = μ ( θ ) u x 2 v U x 2 V + κ ( θ ) θ x v Θ x V x + R 2 , ( ϕ , ψ , ζ ) ( x , 0 ) = ( ϕ 0 , ψ 0 , ζ 0 ) , ( ϕ , ψ , ζ ) ( ± , t ) = ( 0 , 0 , 0 ) ,

where

R 1 = μ ( θ ) U x V x ( P P r P + r ) x U t c , R 2 = μ ( θ ) U x 2 V + κ ( θ ) Θ x V x κ ( Θ c ) Θ x c V c x ( P p m ) U x c ( P P r ) ( U r ) x ( P P + r ) ( U + r ) x .

Then we establish the existence of solution to (4.1). First of all, we introduce a function space X Θ ̲ , Θ ¯ , M ( I ) defined by

X Θ ̲ , Θ ¯ , M ( I ) = ( ϕ , ψ , ζ ) X ( I ) Θ ̲ θ Θ ¯ , v M 1 , sup I ϕ , ψ , ζ γ 1 H 2 M ,

where I [ 0 , + ) and 0 < Θ ̲ , Θ ¯ , M < + . Set

M 0 ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 , C C ( v ̲ , Θ ̲ , Θ ¯ , M 0 ) .

Since the global solution is constructed by the combination of the local existence and the a priori estimate. The proof of the existence of local solution is standard, and to prove the global existence part of Theorem 2.1, we are mainly concerned about the following a priori estimate.

Proposition 4.1

(A priori estimate) For any T > 0 , let ( ϕ , ψ , ζ ) X Θ ̲ , Θ ¯ , M ( [ 0 , T ] ) be the solution to the Cauchy problem (4.1) and μ , κ satisfy (1.3), suppose that δ c C ( γ 1 ) α , δ ± C ( γ 1 ) β for some α > 5 4 , β > 0 . Then there exists a positive constant δ 0 min { δ ¯ , 1 } depending on initial data but independent of T , M such that if γ 1 δ 0 , the following estimate holds

(4.2) sup 0 t T ϕ , ψ , ζ γ 1 ( , t ) H 2 2 + 0 T ϕ x ( , t ) H 1 2 + ( ψ x , ζ x ) ( , t ) H 2 2 d t C ( M 0 ) ,

where C ( M 0 ) is a positive constant that depends only on M 0 but independent of T , M .

Once the Proposition 4.1 is obtained, we can show the global existence of Theorem 2.1 by the definition (2.14) easily.

In what follows, the analysis is always carried out under the a priori assumptions for some positive constants Θ ̲ , Θ ¯ , M :

(4.3) Θ ̲ θ ( x , t ) Θ ¯ , v ( x , t ) M 1 , δ c C ( γ 1 ) α 1 , δ ± C ( γ 1 ) β 1 , N ( T ) sup 0 t T ϕ , ψ , ζ γ 1 H 2 M .

Thus, Proposition 4.1 is an easy consequence of lemmas below.

Lemma 4.2

Under the assumptions of Proposition 4.1, 0 < t T , ε 1 , it holds that

(4.4) Φ v V , ψ , ζ γ 1 2 + 0 t μ ( θ ) v ψ x 2 + κ ( θ ) v ζ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 2 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t ϕ x 2 d s + C ( M ) δ 1 8 ,

where Φ ( s ) = s ln s 1 , C ( M ) is a positive constant that depends only on v ̲ , Θ ̲ , Θ ¯ , M 0 , M .

Proof

Similar to [17], multiplying (4.1) 1 by R Θ ( v 1 V 1 ) , (4.1) 2 by ψ , and (4.1) 3 by θ 1 ζ , then adding the resulting equations together yield

(4.5) R Θ Φ v V + 1 2 ψ 2 + c v Θ Φ θ Θ t + μ ( θ ) Θ v θ ψ x 2 + κ ( θ ) Θ v θ 2 ζ x 2 = H x + Q + R 1 ψ + R 2 ζ θ ,

where

(4.6) H = ( P p ) ψ + μ ( θ ) u x v U x V ψ + κ ( θ ) θ x v Θ x V ζ θ , Q = Q 1 ( ( U r ) x + ( U + r ) x ) Q 2 + Q 3 , Q 1 = ( γ 1 ) P Φ v V + P ϕ 2 v V P Φ Θ θ + ( p P ) ζ θ = P Φ θ V Θ v + γ Φ v V , Q 2 = U x c P ϕ 2 v V p m Φ v V + p m γ 1 Φ Θ θ + ( p P ) ζ θ + ( γ 1 ) ( P r P ) ( U r ) x Φ v V 1 γ 1 Φ Θ θ + ( γ 1 ) ( P + r P ) ( U + r ) x Φ v V 1 γ 1 Φ Θ θ , Q 3 = μ ( θ ) 1 V 1 v U x ψ x + κ ( θ ) Θ x v θ 2 ζ ζ x + κ ( θ ) Θ Θ x v V θ 2 ϕ ζ x κ ( θ ) Θ x 2 v V θ 2 ϕ ζ + 2 μ ( θ ) U x v θ ζ ψ x μ ( θ ) U x 2 v V θ ϕ ζ ,

satisfy

(4.7) Q 1 C ( M ) 1 ( ϕ 2 + ζ 2 ) 0 , Q 2 C ( M ) U x c ϕ 2 + ζ 2 γ 1 + P r P ( U r ) x ( ϕ 2 + ζ 2 ) + P + r P ( U + r ) x ( ϕ 2 + ζ 2 ) C ( M ) ( Θ x x c + Θ x c 2 + P r P ( U r ) x + P + r P ( U + r ) x ) ( ϕ 2 + ζ 2 ) , Q 3 1 2 μ ( θ ) Θ v θ ψ x 2 + 1 2 κ ( θ ) Θ v θ 2 ζ x 2 + C ( M ) ( U x 2 + Θ x 2 ) ( ϕ 2 + ζ 2 ) .

Since

P ± r P = R Θ ± r V ± r Θ V = R Θ ± r ( V r + V c + V + r v m v + m ) V ± r ( Θ r + Θ c + Θ + r θ m θ + m ) V ± r V = R Θ ± r ( V r + V c v m v + m ) V ± r ( Θ r + Θ c θ m θ + m ) V ± r V ,

recalling (3) and (2) in Lemma 3.2, we can deduce that

(4.8) P ± r P ( U ± r ) x C ( V r v m + V c v ± m + Θ r θ m + Θ c θ ± m ) ( U ± r ) x C ( V r v m + V c v ± m + Θ r θ m + Θ c θ ± m ) Ω ± + C ( U ± r ) x Ω c Ω C δ c ( γ 1 ) 1 2 e c 0 ( γ 1 ) 1 ( x + t ) + C δ e c 0 ( x + t ) ,

(4.9) 0 t ( ( U ± r ) x 2 + ( Θ ± r ) x 2 ) ( ϕ 2 + ζ 2 ) d x d s 2 0 t ( ϕ , ζ ) ( ϕ x , ζ x ) ( ( U ± r ) x , ( Θ ± r ) x ) 2 d s C ( M ) δ 0 t ( ϕ x , ζ x ) 2 d s + 1 .

Integrating (4.5) over R × [ 0 , t ] , using (4.8) and (4.9) and taking δ small enough, we can check that

(4.10) Φ v V , ψ , ζ γ 1 2 + 0 t μ ( θ ) v ψ x 2 + κ ( θ ) v ζ x 2 d x d s + 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 2 + C ( M ) δ + C ( M ) 0 t ( Θ x x c + Θ x c 2 + U x c 2 ) ( ϕ 2 + ζ 2 ) d x d s + C ( M ) δ 0 t ϕ x 2 d s + C 0 t R 1 ψ + R 2 ζ d x d s .

Due to (3) in Lemma 3.2 and

( P r + P + r P ) x = R Θ r V r + Θ + r V + r + Θ c V c Θ V x = R ( V r ) x Θ V 2 Θ r ( V r ) 2 + R ( Θ r ) x 1 V r 1 V + R ( V c ) x Θ V 2 Θ c ( V c ) 2 + R ( Θ c ) x 1 V c 1 V + R ( V + r ) x Θ V 2 Θ + r ( V + r ) 2 + R ( Θ + r ) x 1 V + r 1 V ,

we obtain

(4.11) ( P r + P + r P ) x C δ c ( γ 1 ) 1 2 e c 0 ( γ 1 ) 1 ( x + t ) + C δ e c 0 ( x + t ) .

Thus, we can compute from Lemmas 3.1 and 3.2 that

R 1 L 1 = μ ( θ ) U x V x ( P P r P + r ) x U t c L 1 C ζ x U x , Θ x U x , V x U x , U x x ) L 1 + C δ e c 0 t + C δ ( 1 + t ) 1 C ζ x ( ( U r ) x , U x c , ( U + r ) x ) + C δ e c 0 t + C δ 1 8 ( 1 + t ) 7 8 C δ 1 2 ζ x 2 + C δ 1 8 ( 1 + t ) 7 8 ,

and then

(4.12) 0 t R 1 ψ d x d s 2 0 t ψ 1 2 ψ x 1 2 R 1 L 1 d s C ( M ) δ 1 2 0 t ζ x 2 d s + C ( M ) δ 1 8 0 t ψ x 1 2 ( 1 + s ) 7 8 d s C ( M ) δ 1 8 0 t ( ψ x , ζ x ) 2 d s + 1 .

We now go to estimate the last term in (4.10). According to (4.8), Lemma 3.2 and

(4.13) ( p m P ) U x c = R V Θ c Θ V c V V c U x c C δ c ( γ 1 ) 1 2 e c 0 ( γ 1 ) 1 ( x + t ) + C δ e c 0 ( x + t ) ,

one has

(4.14) C R 2 ζ d x C ( M ) δ e c 0 t + C ζ 1 2 ζ x 1 2 ( ( U r ) x , U x c , ( U + r ) x ) 2 + C ( κ ( θ ) κ ( Θ c ) ) Θ x V x ζ d x + C κ ( Θ c ) Θ x V Θ x c V c x ζ d x C ( M ) δ e c 0 t + ζ x 2 + ( 1 + t ) 4 3 + C ( κ ( θ ) θ x κ ( Θ c ) Θ x c ) Θ x V ζ d x + C ( κ ( θ ) κ ( Θ c ) ) Θ x x V Θ x V x V 2 ζ d x + C κ ( Θ c ) Θ x c Θ x V Θ x c V c ζ d x + C κ ( Θ c ) Θ x V Θ x c V c x ζ d x C ( M ) δ ζ x 2 + ( 1 + t ) 4 3 + i = 1 4 I i .

We then estimate I i ( i = 1 , 2 , 3 , 4 ) term by term:

I 1 C ( M ) Θ x ζ ζ x + ( Θ ± r ) x Θ x ζ + ( ζ + Θ ± r θ ± m ) Θ x c Θ x ζ d x 1 16 κ ( θ ) v ζ x 2 d x + C ( M ) Θ x c 2 ζ 2 d x + C ( M ) ζ 1 2 ζ x 1 2 ( Θ ± r ) x ζ x + C ( M ) ζ γ 1 1 2 ζ x γ 1 1 2 γ 1 ( Θ ± r ) x 2 + ( Θ ± r ) x Θ x c + Θ x c 2 Θ ± r θ ± m d x 1 8 κ ( θ ) v ζ x 2 d x + C ( M ) Θ x c 2 ζ 2 d x + C ( M ) δ ( 1 + t ) 4 3 ,

I 2 C ζ + Θ ± r θ ± m ( Θ x x c + ( Θ ± r ) x x + ( ( Θ ± r ) x + Θ x c ) ( ( V ± r ) x + V x c ) ) ζ d x C ( Θ x x c + ( Θ ± r ) x x + ( Θ ± r ) x 2 + ( V ± r ) x 2 + Θ x c 2 ) ζ 2 d x + C ( ( Θ ± r ) x x + ( Θ ± r ) x 2 + ( V ± r ) x 2 + δ e c 0 ( x + t ) + δ c ( γ 1 ) 1 e c 0 ( γ 1 ) 1 ( x + t ) ) ζ d x C ( Θ x x c + Θ x c 2 ) ζ 2 d x + C ζ 1 2 ζ x 1 2 ( ( Θ ± r ) x x L 1 + ( Θ ± r ) x 2 + ( V ± r ) x 2 + δ e c 0 t ) 1 8 κ ( θ ) v ζ x 2 d x + C ( Θ x x c + Θ x c 2 ) ζ 2 d x + C ( M ) δ 1 6 ( 1 + t ) 7 6 ,

I 3 C Θ x c ( Θ x V c Θ x c V ) ζ d x = C Θ x c ( ( Θ ± r ) x V c Θ x c ( V ± r v ± m ) ) ζ d x C ( δ e c 0 ( x + t ) + δ c ( γ 1 ) 1 2 e c 0 ( γ 1 ) 1 ( x + t ) ) ζ d x C ( M ) δ e c 0 t ,

I 4 = C κ ( Θ c ) Θ x Θ x c V + 1 V 1 V c Θ x c x ζ d x C ( ( Θ ± r ) x x + ( Θ ± r ) x V x + ( V ± r ) x Θ x c + ( V ± r v ± m ) ( Θ x x c + V x c Θ x c ) ) ζ d x C ζ 1 2 ζ x 1 2 ( ( Θ ± r ) x x L 1 + ( Θ ± r ) x 2 + ( V ± r ) x 2 ) + C ( δ e c 0 ( x + t ) + δ c ( γ 1 ) 1 e c 0 ( γ 1 ) 1 ( x + t ) ) ζ d x 1 8 κ ( θ ) v ζ x 2 d x + C ( M ) δ 1 6 ( 1 + t ) 7 6 ,

where we have used Cauchy’s inequality and Lemmas 3.1 and 3.2. Then, collecting all estimates of I i ( i = 1 , 2 , 3 , 4 ) , revisiting (4.14), we can check that

(4.15) C 0 t R 2 ζ d x d s 3 8 + C ( M ) δ 0 t κ ( θ ) v ζ x 2 d x d s + C ( M ) δ 1 6 + C ( M ) 0 t ( Θ x x c + Θ x c 2 ) ζ 2 d x d s .

Plugging (4.12) and (4.15) back into (4.10), and taking

3 8 + C ( M ) δ + C ( M ) δ 1 8 1 2 ,

we obtain

(4.16) Φ v V , ψ , ζ γ 1 2 + 0 t μ ( θ ) v ψ x 2 + κ ( θ ) v ζ x 2 d x d s + 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 2 + C ( M ) δ 1 8 + C ( M ) δ 0 t ϕ x 2 d s + C ( M ) 0 t ( Θ x x c + Θ x c 2 + U x c 2 ) ( ϕ 2 + ζ 2 ) d x d s .

Thus, this lemma follows directly from α > 1 2 and the following lemma by choosing c 3 = c 2 2 and taking ε , γ 1 suitably small.□

Lemma 4.3

Motivated by Lemma 5 in [12], for any c 3 0 , c 2 2 and w defined in (3.3), there exists a positive constant C ( M ) such that 0 < t T , 0 < ε < 1 , it holds that

(4.17) 0 t ( ϕ , ψ , ζ ) 2 w 2 d x d s C ( M ) ( γ 1 ) 1 2 + ε 0 t ( ϕ x , ψ x , ζ x ) 2 d s + 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s + C ε .

Proof

The proof of this lemma is divided into the following two parts:

(4.18) 0 t ( p P ) 2 w 2 + ψ 2 w 2 d x d s C ( M ) ( γ 1 ) 1 2 + ε 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ε ,

(4.19) 0 t ζ + P c v ϕ 2 w 2 d x d s ( ε + C ( M ) δ c ) 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C ( M ) + C ε ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s .

In fact, adding (4.19) to (4.18) and choosing ε , γ 1 suitably small thus implies (4.17) easily.

Step 1 ̲ . We first prove (4.18). Denoting by f ( x , t ) = x w 2 ( y , t ) d y , we have

(4.20) f L C ( γ 1 ) 1 2 ( 1 + t ) 1 2 , f t L C ( γ 1 ) 1 2 ( 1 + t ) 3 2 .

Rewriting (4.1) 2 leads to

(4.21) ψ t + ( p P ) x = μ ( θ ) u x v x ( P P r P + r ) x U t c ,

then multiplying (4.21) by ( p P ) f , integrating the resulting identity over R yields

1 2 ( p P ) 2 w 2 d x = ψ t ( p P ) f d x + μ ( θ ) ψ x + U x v ( ( p P ) f ) x d x + ( ( P P r P + r ) x + U t c ) ( p P ) f d x .

Then, integrating the last identity with respect to t , we obtain

(4.22) 0 t ( p P ) 2 w 2 d x d s = 2 0 t ( ψ ( p P ) f ) t d x d s 2 0 t ψ ( p P ) t f d x d s 2 0 t ψ ( p P ) f t d x d s + 2 0 t μ ( θ ) ψ x v ( ( p P ) f ) x d x d s + 2 0 t μ ( θ ) U x v ( ( p P ) f ) x d x d s + 2 0 t ( ( P P r P + r ) x + U t c ) ( p P ) f d x d s i = 5 10 I i .

We estimate I i ( i = 5 , , 10 ) as follows:

I 5 C ( M ) ( ϕ , ψ , ζ , ϕ 0 , ψ 0 , ζ 0 ) 2 ( γ 1 ) 1 2 C ( M ) ( γ 1 ) 1 2 , I 7 C ( M ) 0 t ψ ( ϕ + ζ ) ( γ 1 ) 1 2 ( 1 + s ) 3 2 d x d s C ( M ) ( γ 1 ) 1 2 0 t ( 1 + s ) 3 2 ψ ( ϕ + ζ ) d s C ( M ) ( γ 1 ) 1 2 ,

I 8 = 2 0 t μ ( θ ) ψ x v R ζ P ϕ v x f d x d s + 2 0 t μ ( θ ) ψ x v R ζ P ϕ v w 2 d x d s C ( M ) 0 t ψ x ( ζ x + ϕ x + ( ζ + ϕ ) ( V x + Θ x ) ) ( γ 1 ) 1 2 ( 1 + s ) 1 2 d x d s + C ( M ) 0 t ψ x ( ζ + ϕ ) w 2 d x d s C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 + ( ϕ 2 + ζ 2 ) ( V x 2 + Θ x 2 ) ( 1 + s ) 1 d x d s + C ( M ) 0 t ψ x L 1 + ε ε ( ϕ , ζ ) L w 2 L 1 + ε d s C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 d s + 1 + C ( M ) 0 t ψ x 1 + 3 ε 2 + 2 ε × ψ x x 1 ε 2 + 2 ε ( ϕ , ζ ) 1 2 ( ϕ x , ζ x ) 1 2 ( γ 1 ) 1 2 + ε ( 1 + s ) 1 2 + ε d s C ( M ) ( γ 1 ) 1 2 + ε 1 + 0 t ( ϕ x , ψ x , ζ x ) 2 + ψ x 1 + 3 ε 2 + 2 ε ( ϕ x , ζ x ) 1 2 ( 1 + s ) 1 2 + ε d s C ( M ) ( γ 1 ) 1 2 + ε C ε + 0 t ( ϕ x , ψ x , ζ x ) 2 d s ,

I 9 C ( M ) 0 t U x ( ζ x + ϕ x + ( ζ + ϕ ) ( V x + Θ x ) ) ( γ 1 ) 1 2 ( 1 + s ) 1 2 d x d s + C ( M ) 0 t U x ( ζ + ϕ ) w 2 d x d s C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ζ x ) U x ( 1 + s ) 1 2 d s + 0 t U x ( V x + Θ x ) ( 1 + s ) 1 2 d x d s + C ( M ) 0 t ( U x c , ( U ± r ) x ) L ( ϕ , ζ ) 1 2 ( ϕ x , ζ x ) 1 2 w 2 d s C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ζ x ) 2 d s + 1 , I 10 C ( M ) 0 t ( ( P P r P + r ) x + U t c ) ( γ 1 ) 1 2 ( 1 + s ) 1 2 d x d s C ( M ) ( γ 1 ) 1 2 ,

where we have used Cauchy’s inequality, Gagliardo-Nirenberg’s inequality, Lemmas 3.1, 3.2, and (4.11) and (4.20).

The estimate of I 6 is more subtle. We can obtain from (4.1) 1 that

I 6 = 2 0 t ψ ( p v ϕ t + ( p v P v ) V t + p θ ζ t + ( p θ P θ ) Θ t ) f d x d s = 0 t ( ψ 2 ) x p v f d x d s 2 0 t ψ ( p v P v ) V t f d x d s 2 0 t ψ ( p θ P θ ) Θ t f d x d s 2 0 t ψ p θ ζ t f d x d s i = 1 4 I 6 i .

We estimate I 6 i ( i = 1 , 2 , 3 , 4 ) as follows:

I 6 1 = 0 t ψ 2 ( p v f ) x d x d s = 0 t p v ψ 2 w 2 d x d s + 0 t ψ 2 ( p v v v x + p v θ θ x ) f d x d s 0 t p v ψ 2 w 2 d x d s + C ( M ) 0 t ψ 2 ( γ 1 ) 1 2 ( 1 + s ) 1 2 ( ϕ x + V x + ζ x + Θ x ) d x d s 0 t p v ψ 2 w 2 d x d s + C ( M ) 0 t ψ 2 δ c ( γ 1 ) 1 ( 1 + s ) 1 e c 2 ξ 2 d x d s + C ( M ) 0 t ( γ 1 ) 1 2 ( 1 + s ) 1 2 ψ 1 2 ψ x 1 2 ψ ( ϕ x , ζ x , ( V ± r ) x , ( Θ ± r ) x ) d s 0 t ( p v + C ( M ) δ c ) ψ 2 w 2 d x d s + C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 d s + 1 0 t 1 2 p v ψ 2 w 2 d x d s + C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 d s + 1 ,

I 6 2 C ( M ) 0 t ψ ( ϕ + ζ ) U x ( γ 1 ) 1 2 ( 1 + s ) 1 2 d x d s C ( M ) ( γ 1 ) 1 2 0 t ψ 1 2 ψ x 1 2 ( ϕ + ζ ) ( U x c , ( U ± r ) x ) ( 1 + s ) 1 2 d s C ( M ) ( γ 1 ) 1 2 0 t ψ x 2 d s + 1 , I 6 3 C ( M ) 0 t ψ ( ϕ + ζ ) Θ t ( γ 1 ) 1 2 ( 1 + s ) 1 2 d x d s C ( M ) 0 t ψ ( ϕ + ζ ) ( U x + Θ x x c + Θ x c 2 + ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) ( γ 1 ) 1 2 ( 1 + s ) 1 2 d x d s C ( M ) ( γ 1 ) 1 2 0 t ψ x 1 2 ( U x c , ( U ± r ) x , Θ x x c , Θ x c 2 ) ( 1 + s ) 1 2 d s + C ( M ) ( γ 1 ) 1 2 0 t δ e c 0 ( x + s ) + δ c ( γ 1 ) 1 2 e c 0 ( γ 1 ) 1 ( x + s ) d x d s C ( M ) ( γ 1 ) 1 2 0 t ψ x 2 d s + 1 ,

I 6 4 = 2 c v 0 t ψ p θ P U x p u x + μ ( θ ) u x 2 v + κ ( θ ) θ x v x κ ( Θ c ) Θ x c V c x ( P p m ) U x c ( P P r ) ( U r ) x ( P P + r ) ( U + r ) x ) f d x d s C ( M ) ( γ 1 ) 1 2 0 t ψ ( ( ϕ + ζ ) U x + U x 2 + ψ x 2 ) ( 1 + s ) 1 2 d x d s + 1 c v 0 t ψ 2 ( p θ p f ) x d x d s + 2 c v 0 t ( ψ p θ f ) x κ ( θ ) θ x v κ ( Θ c ) Θ x c V c d x d s + C ( M ) C ( M ) 0 t ψ x 2 d s + C ( M ) + C ( M ) ( γ 1 ) 1 2 0 t ( U x , U x 2 ) L ( 1 + s ) 1 2 d s + C ( M ) c v 0 t ψ 2 w 2 + ( γ 1 ) 1 2 ( 1 + s ) 1 2 ( v x + θ x ) d x d s + C ( M ) c v 0 t ψ x ( γ 1 ) 1 2 ( 1 + s ) 1 2 + ψ ( v x + θ x ) ( γ 1 ) 1 2 ( 1 + s ) 1 2 + ψ w 2 × ( ζ x + ( Θ ± r ) x + ( ϕ + V ± r v ± m + ζ + Θ ± r θ ± m ) Θ x c ) d x d s C ( M ) ( γ 1 ) 1 2 0 t ( ψ x , ζ x ) 2 d s + 1 + C ( M ) ( γ 1 ) 0 t ψ 2 w 2 d x d s ,

where we have used Cauchy’s inequality, Lemmas 3.1, 3.2, c 3 c 2 2 , C ( M ) δ c 1 , and (4.20), (2.13), (4.8), (4.13), and (4.1) 3 . Then, collecting the estimates of I i ( 5 , , 10 ) , recalling (4.22), and taking

C ( M ) ( γ 1 ) 1 4 sup p v ( v , θ ) ,

we obtain

0 t ( p P ) 2 w 2 1 4 p v ψ 2 w 2 d x d s C ( M ) ( γ 1 ) 1 2 + ε 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ε ,

which means (4.18).

Step 2 ̲ . Next, we prove (4.19). Rewriting (4.1) 3 leads to

(4.23) c v ζ t + P ϕ t = ( P p ) u x + μ ( θ ) u x 2 v + κ ( θ ) θ x v x κ ( Θ c ) Θ x c V c x ( P p m ) U x c ( P P r ) ( U r ) x ( P P + r ) ( U + r ) x .

We are supposed to use Lemma 3.3 to finish the proof of (4.19). Taking h = ζ + P c v ϕ in Lemma 3.3 and using (4.23), we know that

(4.24) c v 0 t h t , h g 2 H 1 × H 1 d s = 0 t ( c v ζ t + P ϕ t ) h g 2 d x d s + 0 t P t ϕ h g 2 d x d s = 0 t ( P p ) U x h g 2 d x d s + 0 t ( P p ) ψ x h g 2 d x d s + 0 t μ ( θ ) u x 2 v h g 2 d x d s 0 t κ ( θ ) ζ x + Θ x c v κ ( Θ c ) Θ x c V c ( h g 2 ) x d x d s + 0 t κ ( θ ) ( Θ r ) x + ( Θ + r ) x v x h g 2 d x d s 0 t ( ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x ) h g 2 d x d s + 0 t P t ϕ h g 2 d x d s i = 11 17 .

We estimate I i ( i = 11 , , 17 ) in (4.24) term by term:

I 11 C ( M ) 0 t ( ϕ 2 + ζ 2 ) ( U x c + ( U r ) x + ( U + r ) x ) d x d s C ( M ) 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s + C ( M ) δ c 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s , I 13 C ( M ) 0 t ( ϕ + ζ ) ( U x 2 + ψ x 2 ) d x d s C ( M ) 0 t ψ x 2 d s + C ( M ) 0 t ( ϕ , ζ ) L ( U x c , ( U ± r ) x ) 2 d s C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) ,

I 14 = 0 t κ ( θ ) ζ x + Θ x c v κ ( Θ c ) Θ x c V c ( h x g 2 + 2 h g g x ) d x d s C ( M ) 0 t ( ζ x + ( ϕ + V ± r v ± m + ζ + Θ ± r θ ± m ) Θ x c ) ( ϕ x + ζ x + ( γ 1 ) P x ϕ + ( ϕ + ζ ) w ) d x d s C ε ( M ) 0 t ( ϕ x , ζ x ) 2 d s + ε 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C ( M ) , I 15 = 0 t κ ( θ ) v x ( ( Θ r ) x + ( Θ + r ) x ) h g 2 d x d s + 0 t κ ( θ ) v ( ( Θ r ) x x + ( Θ + r ) x x ) h g 2 d x d s C ( M ) 0 t ( ( ϕ x , ( V ± r ) x , ζ x , ( Θ ± r ) x ) ( Θ ± r ) x + ( Θ ± r ) x x L 1 ) ( ϕ , ζ ) L d s + C ( M ) 0 t Θ x c ( Θ r ) x + ( Θ + r ) x ( ϕ + ζ ) d x d s C ( M ) 0 t ( ϕ x , ζ x ) 2 d s + C ( M ) ,

I 16 C ( M ) 0 t ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x d x d s C ( M ) 0 t δ c ( γ 1 ) 1 2 e c 0 ( γ 1 ) 1 ( x + s ) + δ e c 0 ( x + s ) d x d s C ( M ) , I 17 = 0 t ( P v V t + P θ Θ t ) ϕ h g 2 d x d s C ( M ) 0 t ( ϕ 2 + ζ 2 ) ( U x + Θ x c 2 + Θ x x c + ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x ) d x d s C ( M ) 0 t Q 1 ( ( U r ) x + ( U + r ) x ) + δ c ( ϕ 2 + ζ 2 ) w 2 d x d s + C ( M ) ,

where we have used Cauchy’s inequality, c 3 c 2 2 , the definition of Q 1 in (4.6), Lemmas 3.1, 3.2, and (4.7) 1 , (4.8), (4.13), and (2.13). The estimate of I 12 is more subtle. First, it is easy to see that

(4.25) ( P p ) = R Θ V θ v = R Θ ϕ V ζ v V = R v ζ + P c v ϕ + γ P v ϕ = R v h + γ P v ϕ .

Thus, we compute from (3.4), (4.25), and (4.23) that

I 12 = 0 t R v h + γ P v ϕ ϕ t h g 2 d x d s = 0 t R v ϕ t h 2 g 2 d x d s + 0 t γ P 2 v ( ϕ 2 ) t h g 2 d x d s = 0 t R v ϕ h 2 g 2 + γ P 2 v ϕ 2 h g 2 t d x d s 0 t R v 2 v t ϕ h 2 g 2 2 R v ϕ ( h h t g 2 + h 2 g g t ) + γ P t 2 v ϕ 2 h g 2 γ P 2 v 2 v t ϕ 2 h g 2 + γ P 2 v ϕ 2 ( h t g 2 + 2 h g g t ) d x d s C ( M ) ( ϕ , ζ ) 2 + C + C ( M ) 0 t ( ϕ 3 + ζ 3 ) ( ψ x + U x + P t ) d x d s + 1 c v 0 t 2 R v ϕ h γ P 2 v ϕ 2 g 2 c v h t d x d s + C ( M ) 0 t ( ϕ 3 + ζ 3 ) g t d x d s C ( M ) + C ( M ) 0 t ( ϕ 3 + ζ 3 ) ( ψ x + U x ) d x d s + C ( M ) 0 t ( ϕ 3 + ζ 3 ) Θ t d x d s + C ( M ) ( γ 1 ) 0 t ( ϕ 2 + ζ 2 ) ( u x 2 + ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) d x d s + 1 c v 0 t 2 R v ϕ h γ P 2 v ϕ 2 g 2 κ ( θ ) θ x v κ ( Θ c ) Θ x c V c x d x d s + C ( M ) ( γ 1 ) 0 t ( ϕ 3 + ζ 3 ) w x d x d s C ( M ) + i = 1 5 I 12 i .

We then estimate I 12 i ( i = 1 , , 5 ) term by term. It follows from Cauchy’s inequality, the definition of Q 1 in (4.6), Lemmas 3.1, 3.2, (4.8), (4.13), (2.13), and w x C ( γ 1 ) 1 ( 1 + t ) 1 e c 3 ξ 2 that

I 12 1 C ( M ) 0 t ( ϕ , ζ ) L ( ϕ , ζ ) ( ( ϕ , ζ ) L ψ x + ( ϕ , ζ ) ( ( U ± r ) x , U x c ) L ) d s C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) ,

I 12 2 C ( M ) ( γ 1 ) 0 t ( ϕ 3 + ζ 3 ) ( Θ x x c + Θ x c 2 + ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) d x d s + C ( M ) ( γ 1 ) 0 t ( ϕ 3 + ζ 3 ) U x d x d s C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) ,

I 12 3 C ( M ) 0 t ψ x 2 d s + C ( M ) 0 t ( ϕ , ζ ) ( ϕ x , ζ x ) ( ( U ± r ) x , U x c ) 2 d s + C ( M ) C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) ,

I 12 4 = 1 c v 0 t 2 R v ϕ h γ P 2 v ϕ 2 g 2 x κ ( θ ) θ x v κ ( Θ c ) Θ x c V c d x d s C ( M ) ( γ 1 ) 0 t ( ( ϕ 2 + ζ 2 ) ( w x + v x ) + ( ϕ + ζ ) ( ϕ x + ζ x + P x ϕ ) ) × ( ζ x + ( Θ ± r ) x + ( ϕ + V ± r v ± m + ζ + Θ ± r θ ± m ) Θ x c ) d x d s C ( M ) 0 t ( ϕ x , ζ x ) 2 d s + C ( M ) ( γ 1 ) 0 t ( ϕ 4 + ζ 4 ) ( w x 2 + V x 2 + Θ x 2 ) d x d s + C ( M ) ( γ 1 ) 0 t ( ϕ 2 + ζ 2 ) Θ x c 2 d x d s + C ( M ) + C ( M ) ( γ 1 ) 0 t ( ( ϕ 2 + ζ 2 ) ( w x + V x + Θ x ) + ( ϕ + ζ ) ( ϕ x + ζ x ) ) ( Θ ± r ) x d x d s C ( M ) 0 t ( ϕ x , ζ x ) 2 d s + C ( M ) 0 t ( ϕ , ζ ) 2 ( ϕ x , ζ x ) 2 ( γ 1 ) ( w x 2 , V x 2 , Θ x 2 ) L 1 d s + C ( M ) δ c 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C ( M ) + C ( M ) ( γ 1 ) 0 t ( ϕ x , ζ x ) ( w x , ( V ± r ) x , ( Θ ± r ) x ) ( Θ ± r ) x + ( ϕ x , ζ x ) 3 2 ( Θ ± r ) x d s C ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ζ x ) 2 d s + C ( M ) δ c 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C ( M ) ,

I 12 5 C ( M ) ( γ 1 ) 0 t ( ϕ , ζ ) 3 2 ( ϕ x , ζ x ) 3 2 w x L 1 d s C ( M ) 0 t ( ϕ x , ζ x ) 2 d s + C ( M ) .

Therefore, substituting the aforementioned estimates of I i ( i = 11 , , 17 ) into (4.24) yields

(4.26) 1 γ 1 0 t h t , h g 2 H 1 × H 1 d s ( ε + C ( M ) δ c ) 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C ε ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s + C ( M ) ,

where h = ζ + P c v ϕ . Then, from Lemma 3.3 and (4.26), we can check that

0 t h 2 w 2 d x d t 4 π γ 1 h ( , 0 ) 2 + 2 π c 3 1 0 t h x ( , s ) 2 d s + 8 c 3 γ 1 0 t h t , h g 2 H 1 × H 1 d s C ϕ 0 , ζ 0 γ 1 2 + C 0 t ( ϕ x , ζ x , P x ϕ ) 2 d s + C γ 1 0 t h t , h g 2 H 1 × H 1 d s ( ε + C ( M ) δ c ) 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C ε ( M ) ( γ 1 ) 1 2 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) 0 t Q 1 ( ( U r ) x + ( U + r ) x ) d x d s + C ( M ) ,

we obtain the desired estimate (4.19).

We now go to prove (4.17).

If ϕ ζ 0 , we obtain from the mean value theorem that

(4.27) ( p P ) 2 = ( p v ( v 1 , θ 1 ) ϕ + p θ ( v 2 , θ 2 ) ζ ) 2 = p v 2 ( v 1 , θ 1 ) ϕ 2 + p θ 2 ( v 2 , θ 2 ) ζ 2 + 2 p v ( v 1 , θ 1 ) p θ ( v 2 , θ 2 ) ϕ ζ C ( M ) 1 ( ϕ 2 + ζ 2 ) ,

where v i ( i = 1 , 2 ) are between v and V , θ i ( i = 1 , 2 ) are between θ and Θ . Thus, we obtain the desired estimate (4.17) from (4.18) immediately.

If ϕ ζ 0 , since

ζ + P c v ϕ 2 = ζ 2 + P 2 c v 2 ϕ 2 + 2 P c v ϕ ζ ζ 2 ,

according to (4.27), one has

(4.28) ( p P ) 2 + C ζ + P c v ϕ 2 p v 2 ( v 1 , θ 1 ) ϕ 2 + p θ 2 ( v 2 , θ 2 ) ζ 2 + 2 p v ( v 1 , θ 1 ) p θ ( v 2 , θ 2 ) ϕ ζ + C ζ 2 p v 2 ( v 1 , θ 1 ) ϕ 2 + p θ 2 ( v 2 , θ 2 ) ζ 2 1 2 p v 2 ( v 1 , θ 1 ) ϕ 2 + ( C 2 p θ 2 ( v 2 , θ 2 ) ) ζ 2 1 2 p v 2 ( v 1 , θ 1 ) ϕ 2 + p θ 2 ( v 2 , θ 2 ) ζ 2 C ( M ) 1 ( ϕ 2 + ζ 2 ) ,

where C = 2 sup p θ 2 ( v 2 , θ 2 ) . Therefore, multiplying (4.19) by C , then adding the resulting identity to (4.18) and choosing ε , γ 1 small enough, we can arrive at (4.17) by (4.28).□

Remark 6

Our result in Lemma 4.3 can be extended to the general gas, if e ( v , θ ) , p ( v , θ ) satisfy

e θ ( v , θ ) p v ( v , θ ) < 0 , v , θ > 0 .

Lemma 4.4

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.29) ζ x γ 1 2 + 0 t κ ( θ ) v ζ x x 2 d x d s C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) .

Proof

Multiply (4.1) 3 by ζ x x yields

(4.30) c v 2 ζ x 2 t ( c v ζ t ζ x ) x + κ ( θ ) v ζ x x 2 = p u x P U x κ ( θ ) v Θ x x κ ( θ ) v x θ x μ ( θ ) v u x 2 + R 3 ζ x x ,

where R 3 = κ ( Θ c ) Θ x c V c x + ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x , and then integrating the resulting equation (4.30) over R × [ 0 , t ] , using Cauchy’s inequality, Lemmas 3.1, 3.2, (4.8), (4.13) and the a priori assumptions (4.3), we have

ζ x γ 1 2 + 0 t κ ( θ ) v ζ x x 2 d x d s ζ 0 x γ 1 2 + C ( M ) 0 t ψ x 2 + ( ϕ 2 + ζ 2 ) U x 2 + Θ x x 2 + U x 4 + ψ x 4 + R 3 2 d x d s + C ( M ) 0 t ( V x 2 + ϕ x 2 + Θ x 2 + ζ x 2 ) ( Θ x 2 + ζ x 2 ) d x d s ζ 0 x γ 1 2 + C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) 0 t ( ϕ , ζ ) ( ϕ x , ζ x ) ( U x c , ( U ± r ) x ) 2 d s + C ( M ) 0 t ( U x c , V x c , Θ x c , ( U ± r ) x , ( V ± r ) x , ( Θ ± r ) x ) L 4 4 + ( Θ x x c , ( Θ ± r ) x x , R 3 ) 2 d s C ( M ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C ( M ) .

Consequently, we can obtain (4.29) from Lemma 3.2(2) immediately.□

Now we turn to derive the first-order energy estimate on ϕ .

Lemma 4.5

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.31) μ ( θ ) v ϕ x 2 + 0 t μ ( θ ) v 3 ϕ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 , ϕ 0 x , ζ 0 x γ 1 2 + C ψ 2 + C .

Proof

We rewrite (4.1) 2 as

(4.32) μ ( θ ) v ϕ x t = ψ t + ( p P ) x + μ ( θ ) v ( θ t ϕ x θ x u x ) + μ ( θ ) v 2 V x u x μ ( θ ) v U x x + ( P P r P + r ) x + U t c ,

and multiplying the last equation by μ ( θ ) v ϕ x , we find

μ 2 ( θ ) 2 v 2 ϕ x 2 μ ( θ ) v ϕ x ψ t + R μ ( θ ) θ v 3 ϕ x 2 = μ ( θ ) v ψ ψ x x + μ ( θ ) v ψ x 2 + ψ μ ( θ ) v ( θ x ψ x θ t ϕ x ) μ ( θ ) v 2 V x ψ x + μ ( θ ) v 2 U x ϕ x + R V x Θ V 2 θ v 2 + Θ x 1 v 1 V + ζ x v μ ( θ ) v ϕ x + μ ( θ ) v ( θ t ϕ x θ x u x ) μ ( θ ) v ϕ x + μ ( θ ) v 2 V x u x μ ( θ ) v U x x + ( P P r P + r ) x + U t c μ ( θ ) v ϕ x .

Then, integrating it over R × [ 0 , t ] , we obtain

(4.33) μ ( θ ) v ϕ x 2 + 0 t μ ( θ ) v 3 ϕ x 2 d x d s C ( ψ 0 , ϕ 0 x ) 2 + C ψ 2 + i = 18 22 I i .

We now estimate I i ( i = 18 , , 22 ) in (4.33) as follows:

I 18 C + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t μ ( θ ) v 3 ϕ x 2 d x d s ,

I 19 C ( M ) 0 t ψ ( θ x ψ x + Θ t ϕ x + V x ψ x + U x ϕ x ) + ψ ϕ x ζ t d x d s 1 32 0 t μ ( θ ) v 3 ϕ x 2 d x d s + 0 t μ ( θ ) v ψ x 2 + κ ( θ ) v ζ x 2 d x d s + C ( M ) 0 t ψ 2 ( U x 2 + Θ x 2 + Θ t 2 + V x 2 ) d x d s + C ( M ) 0 t ψ 1 2 ψ x 1 2 ζ x ψ x d s + C ( M ) ( γ 1 ) 0 t ψ ϕ x ( U x + ψ x + u x 2 + Θ x 2 + ζ x 2 + V x 2 + ϕ x 2 + Θ x x + ζ x x + Θ x c 2 + Θ x x c + ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) d x d s 1 16 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C + C ( M ) ( γ 1 ) 1 4 0 t ψ x 3 2 ζ x 1 2 ζ x γ 1 1 2 d s + C ( M ) ( γ 1 ) 0 t ϕ x 2 + ψ x 2 + ζ x 2 + ζ x x 2 d x d s + C ( M ) ( γ 1 ) 0 t ψ 2 ( U x 2 + U x 4 + Θ x 4 + V x 4 + Θ x x 2 + Θ x c 4 + Θ x x c 2 + ( P p m ) U x c 2 + ( P P ± r ) ( U ± r ) x 2 ) d x d s 1 16 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε + γ 1 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C + C ( M ) δ c 0 t ψ 2 w 2 d x d s + C ( M ) ( γ 1 ) 0 t ψ ψ x ( ( U ± r ) x , ( Θ ± r ) x , ( V ± r ) x , ( Θ ± r ) x x ) 2 d s 1 8 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε + γ 1 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ,

I 20 1 16 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ( M ) 0 t ( V x 2 + Θ x 2 ) ( ϕ 2 + ζ 2 ) d x d s + R 0 t μ ( θ ) v 2 ϕ x ζ x d x d s 1 8 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C + C 0 t κ ( θ ) v ζ x 2 d x d s + C ( M ) 0 t ( ϕ , ζ ) ( ϕ x , ζ x ) ( ( V ± r ) x , ( Θ ± r ) x ) 2 d s 1 8 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ,

I 21 1 8 + C ( M ) δ 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ( M ) 0 t ( Θ x 2 + ζ x 2 ) ( U x 2 + ψ x 2 ) d x d s + C ( M ) ( γ 1 ) 0 t ϕ x 2 ( U x + ψ x + u x 2 + Θ x 2 + ζ x 2 + V x 2 + ϕ x 2 + Θ x x + ζ x x + Θ x c 2 + Θ x x c + ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) d x d s 1 8 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C + C ( M ) ( γ 1 ) 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ( M ) ( γ 1 ) 0 t ϕ x L ϕ x ζ x x d s 1 8 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε + γ 1 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ,

I 22 1 8 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ( M ) ( δ c ( γ 1 ) 1 2 + δ ) 0 t ( ϕ x , ψ x ) 2 d s + C ( M ) 0 t V x 2 U x 2 + U x x 2 + ( P P r P + r ) x 2 + U t c 2 d x d s 1 8 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε 0 t μ ( θ ) v 3 ϕ x 2 d x d s + C ,

where we have used Cauchy’s inequality, Lemmas 3.1, 3.2, 4.2, 4.3, 4.4, (4.1) 3 , (2.13) 3 , Θ t δ , (4.8), (4.13), (4.11) and the a priori assumptions (4.3). Therefore, plugging the aforementioned estimates I i ( i = 18 , 22 ) back into (4.33), and taking

1 2 + C ( M ) δ + δ c ( γ 1 ) 1 2 + ε + γ 1 3 4 ,

we complete the proof of this lemma.□

Choosing ( γ 1 ) small enough, we now deduce from Lemmas 4.2 and 4.5 that

Lemma 4.6

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.34) Φ v V , ψ , ζ γ 1 , μ ( θ ) v ϕ x 2 + 0 t μ ( θ ) v ψ x 2 + κ ( θ ) v ζ x 2 + μ ( θ ) v 3 ϕ x 2 d x d s C ( ϕ 0 , ψ 0 , ζ 0 γ 1 , ϕ 0 x , ζ 0 x γ 1 ) 2 + C .

Now we use the above estimates to yield an estimate on the positive lower and upper bounds on the specific volume v ( x , t ) as follows:

Lemma 4.7

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.35) C 1 v ( x , t ) C ,

where C is a positive constant that depends only on initial data and v ̲ , Θ ̲ , Θ ¯ .

Proof

Denote v ˜ = v V , due to

v ˜ x v ˜ = ϕ x v ϕ V x v V ,

we can obtain from Lemmas 3.1, 3.2, and 4.6 that

(4.36) v ˜ x v ˜ 2 C μ ( θ ) v ϕ x 2 + C ( M ) ( ( V r ) x , V x c , ( V + r ) x ) 2 C ϕ 0 , ψ 0 , ζ 0 γ 1 , ϕ 0 x , ζ 0 x γ 1 2 + C .

To apply Kanel’s method [22] to the proof of (4.35), we let

A ( v ˜ ) = 1 v ˜ Φ ( s ) s d s , Φ ( s ) = s ln s 1 .

Moreover, one can deduce from (4.36) and Lemma 4.6 that

(4.37) A ( v ˜ ( x , t ) ) = x A ( v ˜ ( y , t ) ) y d y Φ ( v ˜ ) v ˜ x v ˜ d x Φ ( v ˜ ) v ˜ x v ˜ C ϕ 0 , ψ 0 , ζ 0 γ 1 , ϕ 0 x , ζ 0 x γ 1 2 + C .

Therefore, it is straightforward to obtain (4.35) from (4.37).□

Next, we estimate L 2 norm of ( ψ x , ζ x ) .

Lemma 4.8

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.38) ψ x , ζ x γ 1 2 + 0 t ( ψ x x , ζ x x ) 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + C .

Proof

Multiply (4.1) 2 by ψ x x yields

1 2 ψ x 2 t ( ψ t ψ x ) x + μ ( θ ) v ψ x x 2 = R ζ P ϕ v x ψ x x + ( P P r P + r ) x + U t c μ ( θ ) v U x x μ ( θ ) v x u x ψ x x ,

then integrating the resulting equation over R × [ 0 , t ] , we know that

(4.39) ψ x 2 + 0 t ψ x x 2 d s C ψ 0 x 2 + C 0 t ϕ x 2 + ζ x 2 + ( ϕ 2 + ζ 2 ) ( V x 2 + Θ x 2 ) + ( P P r P + r ) x 2 + U t c 2 + U x x 2 + ψ x 2 + ( V x 2 + Θ x 2 ) U x 2 + ( ϕ x 2 + ζ x 2 ) ψ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + C + 1 4 0 t ψ x x 2 d s ,

where we have used Cauchy’s inequality, Lemmas 3.1, 3.2, 4.3, 4.6, 4.7, (4.11), γ 1 1 ,

(4.40) 0 t ( ϕ 2 + ζ 2 ) ( V x 2 + U x 2 + Θ x 2 ) d x d s C δ c 0 t ( ϕ 2 + ζ 2 ) w 2 d x d s + C 0 t ( ϕ , ζ ) ( ϕ x , ζ x ) ( ( V ± r ) x , ( U ± r ) x , ( Θ ± r ) x ) 2 d s C 0 t ( ϕ x , ψ x , ζ x ) 2 d s + C

and

(4.41) C 0 t ( ϕ x 2 + ζ x 2 ) ψ x 2 d x d s C 0 t ψ x ψ x x ϕ x 2 d s + C 0 t ζ x ζ x x ψ x 2 d s 1 4 0 t ψ x x 2 d s + ( C + C ( M ) ( γ 1 ) ) 0 t ψ x 2 d s .

Next, integrating (4.30) over R × [ 0 , t ] , using Cauchy’s inequality, Lemmas 3.1, 4.3, 4.6, 4.7, and (4.39), (4.40), we have

(4.42) ζ x γ 1 2 + 0 t ζ x x 2 d s C ζ 0 x γ 1 2 + C 0 t ψ x 2 + ( ϕ 2 + ζ 2 ) U x 2 + ( Θ ± r ) x x 2 + ( V x 2 + Θ x 2 ) Θ x 2 + ϕ x 2 + ζ x 2 + ( ϕ x 2 + ζ x 2 ) ζ x 2 + U x 4 + ψ x 4 + Θ x c 4 + ( P p m ) U x c 2 + ( P P ± r ) ( U ± r ) x 2 + ( V ± r v ± m + Θ ± r θ ± m ) 2 Θ x x c 2 + ( ϕ + ζ ) Θ x x c ζ x x d x d s C ( M ) δ c ( γ 1 ) 1 2 0 t ζ x x 2 d s + C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + C ,

where we have used α > 1 2 ,

0 t ( ϕ + ζ ) Θ x x c ζ x x d x d s C 0 t ( ϕ , ζ ) 1 2 ( ϕ x , ζ x ) 1 2 Θ x x c ζ x x d s C 0 t ( ϕ x , ζ x ) 1 2 δ c ( γ 1 ) 3 4 ( 1 + s ) 3 4 ζ x x 1 2 ζ x x γ 1 1 2 ( γ 1 ) 1 4 d s C ( M ) δ c ( γ 1 ) 1 2 0 t ( ϕ x , ζ x , ζ x x ) 2 + ( 1 + s ) 3 2 d s

and

(4.43) 0 t ψ x 4 + ( ϕ x 2 + ζ x 2 ) ζ x 2 d x d s C 0 t ψ x ψ x x ψ x 2 d s + C 0 t ζ x ζ x x ( ϕ x , ζ x ) 2 d s C 0 t ψ x x 2 d s + ( C + C ( M ) ( γ 1 ) ) 0 t ( ϕ x , ψ x , ζ x ) 2 d s .

Consequently, we arrive at the conclusion of this lemma by (4.39) and (4.42).□

We now deal with the second-order energy estimate on ( ψ , ζ ) .

Lemma 4.9

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.44) ψ x x , ζ x x γ 1 2 + 0 t ( ψ x x x , ζ x x x ) 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 exp C 0 t ϕ x x 2 d s .

Proof

First, in a similar analysis to (4.11), (4.8), and (4.13), one has

(4.45) ( P r + P + r P ) x x C δ c ( γ 1 ) 1 e c 0 ( γ 1 ) 1 ( x + t ) + C δ e c 0 ( x + t )

and

(4.46) ( ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x ) x C δ c ( γ 1 ) 1 e c 0 ( γ 1 ) 1 ( x + t ) + C δ e c 0 ( x + t ) .

We differentiate (4.1) 2 with respect to x once and multiply the final result by ψ x x x to obtain that

1 2 ψ x x 2 t ( ψ x t ψ x x ) x + μ ( θ ) v ψ x x x 2 = R ζ P ϕ v x x ψ x x x + ( P P r P + r ) x x + U x t c μ ( θ ) v U x x x 2 μ ( θ ) v x u x x μ ( θ ) v x x u x ψ x x x .

By integrating this identity over R × [ 0 , t ] , using Cauchy’s inequality, Lemmas 3.1, 3.2, 4.64.8, and (4.45), one has

ψ x x 2 + 0 t ψ x x x 2 d s C ψ 0 x x 2 + C 0 t ϕ x x 2 + ζ x x 2 + ϕ x 2 + ζ x 2 + ( ϕ x 2 + ζ x 2 ) ϕ x 2 + ( ϕ 2 + ζ 2 ) ( V x x 2 + Θ x x 2 + V x 4 + Θ x 4 ) + ( P P r P + r ) x x 2 + U x t c 2 + U x x x 2 + U x x 2 + ψ x x 2 + ( ϕ x 2 + ζ x 2 ) ψ x x 2 + ( v x 4 + θ x 4 + v x x 2 + θ x x 2 ) u x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ϕ x x 2 d s + C ( M ) 0 t ( V x x , Θ x x , V x 2 , Θ x 2 , U x t c , U x x x , U x x ) 2 d s + C 0 t δ c ( γ 1 ) 1 e c 0 ( γ 1 ) 1 ( x + s ) + C δ e c 0 ( x + s ) d x d s + C 0 t ( ϕ x 2 + ζ x 2 ) ( ϕ x 2 + ψ x x 2 ) + ( ϕ x 4 + ζ x 4 + ϕ x x 2 + ζ x x 2 + V x x 2 + Θ x x 2 ) ψ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ϕ x x 2 d s + C 0 t ( ϕ x , ζ x ) 2 ( ϕ x ϕ x x + ψ x x ψ x x x ) d s + C 0 t ( ϕ x 2 ϕ x x 2 + ζ x 2 ζ x x 2 ) ψ x 2 d s + C 0 t ( ϕ x x , ζ x x , V x x , Θ x x ) 2 ψ x ψ x x d s 1 2 0 t ψ x x x 2 d s + C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ϕ x x 2 d s + C 0 t ( ϕ x x , ζ x x ) 2 ψ x x 2 d s ,

where we need that α > 3 4 . Then, applying Gronwall’s inequality (Lemma 3.4), we can check that

(4.47) ψ x x 2 + 0 t ψ x x x 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ϕ x x 2 d s exp C 0 t ϕ x x 2 d s .

We have, by differentiating (4.1) 3 with respect to x once and by multiplying the result by ζ x x x , that

c v 2 ζ x x 2 t ( c v ζ x t ζ x x ) x + κ ( θ ) v ζ x x x 2 = ( p u x P U x ) x + κ ( Θ c ) V c Θ x c x x κ ( θ ) v Θ x x x ζ x x x + ( ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x ) x ζ x x x μ ( θ ) v u x 2 x + 2 κ ( θ ) v x θ x x + κ ( θ ) v x x θ x ζ x x x .

Integrating it over R × [ 0 , t ] yields

(4.48) ζ x x γ 1 2 + 0 t ζ x x x 2 d s C ζ 0 x x γ 1 2 + C 0 t ψ x x 2 + U x x 2 + ( ϕ x 2 + V x 2 + ζ x 2 + Θ x 2 ) ( ψ x 2 + U x 2 ) + Θ x x x c 2 + Θ x c Θ x x c 2 + Θ x c 6 + Θ x x x 2 + ( ( P p m ) U x c + ( P P r ) ( U r ) x + ( P P + r ) ( U + r ) x ) x 2 + ( v x 2 + θ x 2 ) u x 4 + u x 2 u x x 2 + ( v x 2 + θ x 2 ) θ x x 2 + ( v x x 2 + θ x x 2 + v x 4 + θ x 4 ) θ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ( ϕ x 2 + ζ x 2 ) ψ x 2 + ( ϕ x 2 + ζ x 2 ) ψ x 4 + ψ x 2 ψ x x 2 + ( ϕ x 2 + ζ x 2 ) ζ x x 2 + ϕ x x 2 + ( ϕ x 4 + ζ x 4 + ϕ x x 2 + ζ x x 2 ) ζ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ϕ x x 2 d s + C 0 t ( ϕ x , ζ x ) 2 ψ x 2 ψ x x 2 + ψ x 2 ψ x x ψ x x x + ( ϕ x , ζ x ) 2 ζ x x ζ x x x d s + C 0 t ϕ x 2 ϕ x x 2 ζ x 2 d s + C ( M ) ( γ 1 ) 0 t ( ζ x , ϕ x x , ζ x x ) 2 d s 1 2 0 t ζ x x x 2 d s + C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ϕ x x 2 + ψ x x x 2 d s ,

where we have used α > 5 4 , Lemmas 3.1, 3.2, 4.64.8, (4.41), (4.45), and C ( M ) ( γ 1 ) 1 . Therefore, collecting the estimates (4.47) and (4.48), we obtain (4.44) immediately.□

The estimate on ϕ x x is more subtle. We first introduce the following lemma which is very important to estimate ϕ x x .

Lemma 4.10

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.49) 0 t ψ t , ζ t γ 1 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + 1 ,

(4.50) ψ t 2 + 0 t ψ x t 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 ,

(4.51) ζ t γ 1 2 + 0 t ζ x t 2 d s C ( M ) ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 ,

(4.52) ζ t γ 1 2 + 0 t ζ x t 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 exp C 0 t ϕ x x 2 d s .

Proof

Step 1 ̲ . We first prove (4.49). It is apparent from (4.1) 2 , (4.1) 3 , Lemmas 3.1, 3.2, 4.3, 4.64.8, and (4.8), (4.11), (4.13), (4.41), and (4.43) that

0 t ψ t 2 d s C 0 t ϕ x 2 + ζ x 2 + ( ϕ 2 + ζ 2 ) ( V x 2 + Θ x 2 + ϕ x 2 ) + ( v x 2 + θ x 2 ) ( ψ x 2 + U x 2 ) + ψ x x 2 + U x x 2 + ( P P r P + r ) x 2 + U t c 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + 1

and

0 t ζ t γ 1 2 d s C 0 t ψ x 2 + ( ϕ 2 + ζ 2 ) U x 2 + U x 4 + ψ x 4 + ( V x 2 + ϕ x 2 ) ( Θ x 2 + ζ x 2 ) + Θ x 4 + ζ x 4 + Θ x x 2 + ζ x x 2 + Θ x c 4 + Θ x x c 2 + ( P p m ) U x c 2 + ( P P ± r ) ( U ± r ) x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + 1 .

Step 2 ̲ . We differentiate (4.1) 2 with respect to t once and multiply the final result by ψ t to obtain that

1 2 ψ t 2 t = P ϕ R ζ v + μ ( θ ) v ψ x + μ ( θ ) v U x t ψ t x + R ζ P ϕ v μ ( θ ) v ψ x μ ( θ ) v U x t ψ x t ( ( P P r P + r ) x t + U t t c ) ψ t .

Integrating over R × [ 0 , t ] yields

ψ t ( , t ) 2 + 0 t ψ x t 2 d s C ψ t ( , 0 ) 2 + C 0 t ζ t 2 + ψ x 2 + ( ϕ 2 + ζ 2 ) ( U x 2 + Θ t 2 + ψ x 2 ) + U x t 2 + ( U x 2 + ψ x 2 + Θ t 2 + ζ t 2 ) ( ψ x 2 + U x 2 ) + ( P P r P + r ) x t 2 + U t t c 2 + ψ t 2 d x d s C ψ t ( , 0 ) 2 + C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + C + C 0 t ( ϕ 2 + ζ 2 ) ( Θ t 2 + U x 2 ) + U x t 2 + Θ t 2 U x 2 + ζ t 2 ψ x 2 + ( P P r P + r ) x t 2 d x d s C ψ t ( , 0 ) 2 + C ϕ 0 , ψ 0 , ζ 0 γ 1 H 1 2 + C + C 0 t ( ϕ , ζ ) ( ϕ x , ζ x ) ( Θ t , ( U ± r ) x ) 2 d s + C ( M ) ( γ 1 ) 0 t ( 1 + ζ x x 2 ) ψ x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C ,

where we have used Lemmas 3.1, 3.2, 4.64.8, (2.13), (4.43), (4.49), and (4.1) 3 ,

ψ t ( , 0 ) = ( P p ) x + μ ( θ ) u x v x ( P P r P + r ) x U t c ( , 0 ) C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 + C

and

U x t + Θ t U x + ( P P r P + r ) x t C ( U x t c + Θ x t c + Θ x c V t c + Θ t c V x c + V x t c + V x c V t c + ( U ± r ) x t + ( Θ t c + ( Θ ± r ) t ) ( U x c + ( U ± r ) x ) + ( Θ ± r ) x t + ( Θ ± r ) x ( V ± r ) t + ( Θ ± r ) t ( V ± r ) x + ( V ± r ) x t + ( V ± r ) x ( V ± r ) t ) C δ c ( γ 1 ) 1 2 ( 1 + t ) 3 2 e c 2 x 2 ( γ 1 ) ( 1 + t ) + ( U ± r ) x 2 + ( P ± r ) x x + ( P ± r ( U ± r ) x ) x + ( Θ ± r ) x ( U ± r ) x + ( U ± r ) x ( V ± r ) x + ( U ± r ) x x .

Therefore, we complete the proof of (4.50).

Step 3 ̲ . We then prove (4.52). Differentiate (4.1) 3 with respect to t once and multiply the final result by ζ t to obtain that

c v 2 ζ t 2 t = μ ( θ ) v u x 2 p ψ x + P ϕ R ζ v U x t ζ t + κ ( θ ) v θ x κ ( Θ c ) V c Θ x c t ζ t x κ ( θ ) v θ x κ ( Θ c ) V c Θ x c t ζ x t ( ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) t ζ t .

Integrating the last identity over R × [ 0 , t ] , then from Lemmas 3.1, 3.2, 4.64.9, and (4.41), (4.43), (4.49), and (4.50), it is clear that

(4.53) ζ t γ 1 ( , t ) 2 + 0 t ζ x t 2 d s C ζ t γ 1 ( , 0 ) 2 + C 0 t ( ( v t + θ t ) u x 2 + u x u x t + ( v t + θ t ) ψ x + ψ x t ) ζ t + ( ( ϕ t + ζ t + V t + Θ t ) U x + U x t ) ζ t + ( v t 2 + θ t 2 ) θ x 2 + Θ x t 2 + Θ x c Θ t c 2 + Θ x t c 2 + ζ t 2 + ( ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) t 2 d x d s C ζ t γ 1 ( , 0 ) 2 + C + C 0 t ( ζ t , ψ x , ψ x t , ζ x ) 2 d s + C 0 t u x 6 + ζ t 2 ψ x 2 + Θ t 2 u x 4 + ψ x 2 ψ x t 2 + U x 2 U x t 2 + ψ x 4 + ψ x ζ t 2 + ( V t 2 + Θ t 2 ) ( U x 2 + Θ x 2 ) + U x t 2 + ( ψ x 2 + ζ t 2 ) ζ x 2 + ( Θ ± r ) x t 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C 0 t ψ x 4 ψ x x 2 + ψ x ψ x x ψ x t 2 d s + C sup 0 t T ( ψ x L 2 + ψ x L + ζ x L 2 ) 0 t ( γ 1 ) ζ t γ 1 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C + C sup 0 t T ( ψ x 2 + ψ x x 2 ) 0 t ψ x t 2 d s + C ( M ) ( γ 1 ) 0 t ζ t γ 1 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 exp C 0 t ϕ x x 2 d s ,

where we have used

ζ t γ 1 ( , 0 ) = P U x p u x + μ ( θ ) u x 2 v + κ ( θ ) θ x v x κ ( Θ c ) Θ x c V c x ( P p m ) U x c ( P P r ) ( U r ) x ( P P + r ) ( U + r ) x ( , 0 ) C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 + C ,

( ( P p m ) U x c + ( P P ± r ) ( U ± r ) x ) t C ( U x t c + ( Θ t + U x ) U x c + ( U ± r ) x t + ( Θ t c + ( Θ ± r ) t + U x c + ( U ± r ) x ) ( U ± r ) x ) C δ c ( 1 + t ) 2 e c 2 x 2 ( γ 1 ) ( 1 + t ) + δ e c 0 ( x + t ) + ( P ± r ) x x + ( U ± r ) x 2

and

U x 6 + Θ t 2 U x 4 + U x 2 U x t 2 + ( V t 2 + Θ t 2 ) ( U x 2 + Θ x 2 ) + U x t 2 + ( Θ ± r ) x t 2 C ( U x 4 + U x t 2 + Θ t 2 U x 2 + U x 2 Θ x 2 + Θ t 2 Θ x 2 + ( Θ ± r ) x t 2 ) C δ c ( 1 + t ) 3 e c 2 x 2 ( γ 1 ) ( 1 + t ) + δ e c 0 ( x + t ) + ( U ± r ) x 4 + ( P ± r ) x x 2 + δ c ( 1 + t ) 1 e c 2 x 2 ( γ 1 ) ( 1 + t ) ( U ± r ) x 2 + ( U ± r ) x 2 ( Θ ± r ) x 2 + δ c ( 1 + t ) 1 e c 2 x 2 ( γ 1 ) ( 1 + t ) ( Θ ± r ) x 2 + ( P ± r ( U ± r ) x ) x 2 .

Thus, (4.53) means (4.52) is right.

Step 4 ̲ . We now estimate 0 t ψ x ψ x x ψ x t 2 d s in (4.53) in different way,

0 t ψ x ψ x x ψ x t 2 d s C ( M ) 0 t ψ x t 2 d s C ( M ) ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 .

Therefore, (4.51) follows directly from the above estimate and (4.53).□

We now turn to deal with the second-order energy estimate on ϕ .

Lemma 4.11

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.54) ϕ x x 2 + 0 t ϕ x x 2 d s C ( M 0 ) ,

where C ( M 0 ) is a positive constant that depends only on M 0 but independent of T , M .

Proof

Differentiating (4.1) 2 with respect to x once, we can arrive at

ψ x t + R ζ P ϕ v x x = μ ( θ ) v ψ x x x + x = μ ( θ ) v ϕ x x t μ ( θ ) v t ϕ x x + 2 μ ( θ ) v x ψ x x + μ ( θ ) v x x ψ x + x ,

where = μ ( θ ) U x v x ( P P r P + r ) x U t c . Then we multiply the last result by μ ( θ ) v ϕ x x and integrate over R × [ 0 , t ] to obtain

(4.55) ϕ x x 2 + 0 t ϕ x x 2 d s C ϕ 0 x x 2 + C 0 t ψ x t 2 + ζ ϕ x x 2 + ζ x x 2 + ϕ x 2 + P x x 2 + ( ζ x 2 + ϕ x 2 ) ϕ x 2 + ζ x 2 + P x 2 V x 2 + V x x 2 + V x 4 + ( v t + θ t ) ϕ x x 2 + ( v x 2 + θ x 2 ) ψ x x 2 + ( v x 4 + θ x 4 + v x x 2 + θ x x 2 ) ψ x 2 + x 2 d x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 + C γ 1 0 t ϕ x x 2 d s + C 0 t ( ϕ x 2 + ζ x 2 ) ϕ x ϕ x x d s + C δ 0 t ϕ x x 2 d s + C 0 t ψ x 1 2 ψ x x 1 2 + ζ t 1 2 ζ x t 1 2 ϕ x x 2 d s + C 0 t ( ϕ x ϕ x x + ζ x ζ x x ) ψ x x 2 d s + C ( M ) ( δ c ( γ 1 ) 1 + δ + γ 1 ) 0 t ψ x 2 d s + C 0 t ϕ x 2 ϕ x x 2 ψ x 2 + ( ϕ x x 2 + ζ x x 2 ) ψ x ψ x x d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 + 1 4 + C γ 1 + C δ 0 t ϕ x x 2 d s + C 0 t ψ x 2 + ψ x x 2 + ζ t γ 1 2 + ( γ 1 ) ζ x t 2 ϕ x x 2 d s ,

where we have used (4.50) and (2.13), Lemmas 3.1, 3.2, 4.64.8, C ( M ) ( δ + γ 1 ) 1 ,

R ζ P ϕ v x x = R ζ x x P ϕ x x 2 P x ϕ x P x x ϕ v 2 R ζ x P ϕ x P x ϕ v 2 v x + ( R ζ P ϕ ) v x x v 2 + 2 v x 2 v 3 P V ϕ x x v 2 + C ( ζ ϕ x x + ζ x x + P x ϕ x + P x x ϕ + ( ζ x + ϕ x + P x ϕ ) v x + ( ζ + ϕ ) ( V x x + v x 2 ) )

and

x = μ ( θ ) U x v x ( P P r P + r ) x U t c x C ( θ x 2 U x + θ x x U x + θ x v x U x + θ x U x x + v x x U x + v x 2 U x + v x U x x + U x x x + Θ x x c + Θ x c 2 + ( Θ ± r ) x x + ( V ± r ) x x + ( Θ ± r ) x 2 + ( V ± r ) x 2 + U x t c ) C ( ζ x 2 + ζ x x + ϕ x + ζ x + ϕ x ζ x + δ ϕ x x + ϕ x 2 + U x x + U x x x + Θ x x c + Θ x c 2 + ( Θ ± r ) x x + ( V ± r ) x x + ( Θ ± r ) x 2 + ( V ± r ) x 2 + U x t c ) .

Then, if we take δ , γ satisfy

(4.56) 1 4 + C γ 1 + C δ 1 2 ,

we can obtain from (4.55) that

ϕ x x 2 + 0 t ϕ x x 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 + C 0 t ψ x 2 + ψ x x 2 + ζ t γ 1 2 + ( γ 1 ) ζ x t 2 ϕ x x 2 d s .

Therefore, employing Gronwall’s inequality (Lemma 3.4), Lemmas 4.64.8, (4.49), (4.51), and C ( M ) ( γ 1 ) 1 , we can acquire

ϕ x x 2 + 0 t ϕ x x 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 exp C 0 t ψ x 2 + ψ x x 2 + ζ t γ 1 2 + ( γ 1 ) ζ x t 2 d s C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + 1 exp C ϕ 0 , ψ 0 , ζ 0 γ 1 H 2 2 + C C ( M 0 ) ,

which proves this lemma.□

Remark 7

To deduce an estimate on ϕ x x , we need to deal with the term 0 t θ t ϕ x x 2 d x d s in (4.55). In [29], the authors take advantage of θ t = γ 1 R κ ( θ ) v ζ x x + O ( 1 ) ( γ 1 ) Θ x x + u x + u x 2 + κ ( θ ) v x θ x . It is easy to see that to bind such a term ( γ 1 ) 0 t ζ x x ϕ x x 2 d x d s , they need to deduce an estimate on ζ x x L , and as a result, they had to close the energy-type estimates in H 3 ( R ) . In Lemma 4.11 of this article, using θ t = Θ t + ζ t and ζ t L C ζ t 1 2 ζ x t 1 2 , we obtain

0 t θ t ϕ x x 2 d x d s C 0 t ( Θ t L + ζ t 1 2 ζ x t 1 2 ) ϕ x x 2 d s 1 4 0 t ϕ x x 2 d s + C 0 t ζ t ζ x t ϕ x x 2 d s .

According to Lemma 4.10, one has 0 t ζ t ζ x t d s C . Thus, we can close the desired estimates by Gronwall’s inequality. In fact, from Lemmas 4.6–4.11, we know that this method only requires ϕ 0, ψ 0, ζ 0 in H 2 ( R ) .

In conclusion , collecting Lemmas 4.64.11, we can obtain the following lemma, then the Proposition 4.1 follows immediately.

Lemma 4.12

Under the assumptions of Proposition 4.1, 0 < t T , it holds that

(4.57) ϕ , ψ , ζ γ 1 H 2 2 + ψ t , ζ t γ 1 2 + 0 t ϕ x H 1 2 + ( ψ x , ζ x ) H 2 2 + ψ t , ζ t γ 1 , ψ x t , ζ x t 2 d s C ( M 0 ) .

Proof of Theorem 2.1

From Lemma 4.12, we derive

ζ ( x , t ) C γ 1 ζ γ 1 1 2 ζ x γ 1 1 2 C ( M 0 ) γ 1 .

Since

θ ( x , t ) = Θ ( x , t ) + ζ ( x , t ) , 2 Θ ̲ Θ ( x , t ) 1 2 Θ ¯ ,

then letting C ( M 0 ) γ 1 min Θ ̲ , 1 2 Θ ¯ , we obtain

Θ ̲ θ ( x , t ) Θ ¯ .

Furthermore, according to Lemma 4.12, we can check that

(4.58) 0 ( ϕ x , ψ x , ζ x ) 2 + d d t ( ϕ x , ψ x , ζ x ) 2 d t 2 0 ( ϕ x , ψ x , ζ x , ψ x x , ψ x t , ζ x t ) 2 d t C ( M 0 ) < ,

which means

( ϕ x , ψ x , ζ x ) ( , t ) 0 , as t .

Consequently, we have

(4.59) ( ϕ , ψ , ζ ) ( , t ) L 2 ( ϕ , ψ , ζ ) 1 2 ( ϕ x , ψ x , ζ x ) 1 2 0 , as t .

Thus, the proof of Theorem 2.1 is completed.□

Acknowledgments

The authors would like to thank the referees for the valuable comments and suggestions which greatly improved the presentation of the article.

  1. Funding information: Z. H. Guo is supported by the National Natural Science Foundation of China under Grant No. 11931013 and the Natural Science Foundation of Guangxi Province of China under Grant No. 2022GXNSFDA035078.

  2. Conflict of interest: The authors declare that there is no conflict of interests regarding the publication of this article.

References

[1] H. Amann and G. Metzen, Ordinary Differential Equations: An Introduction to Nonlinear Analysis, De Gruyter Studies in Mathematics, vol. 13, De Gruyter, Berlin, 1990. 10.1515/9783110853698Search in Google Scholar

[2] F. V. Atkinson and L. A. Peletier, Similarity profiles of flows through porous media, Arch. Ration. Mech. Anal. 42 (1971), 369–379. 10.1007/BF00250442Search in Google Scholar

[3] F. V. Atkinson and L. A. Peletier, Similarity solutions of the nonlinear diffusion equation, Arch. Ration. Mech. Anal. 54 (1974), 373–392. 10.1007/BF00249197Search in Google Scholar

[4] C. Cercignani, R. Illner, and M. Pulvirenti, The mathematical theory of dilute gases, Applied Mathematical Sciences, vol. 106, Springer-Verlag New York, New York, 1994. 10.1007/978-1-4419-8524-8Search in Google Scholar

[5] S. Chapman, T. G. Cowling, and C. Cercignani, The Mathematical Theory of Non-uniform Gases: An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion in Gases, 3rd ed., Cambridge Mathematical Library, Cambridge University Press, Cambridge, 1995. Search in Google Scholar

[6] W. C. Dong and Z. H. Guo, Global stability of viscous contact wave for one-dimensional compressible Navier-Stokes equations with temperature dependent transport coefficients and large data, preprint, 2022. 10.1515/anona-2022-0246Search in Google Scholar

[7] R. Duan, H. X. Liu, and H. J. Zhao, Nonlinear stability of rarefaction waves for the compressible Navier-Stokes equations with large initial perturbation, Trans. Amer. Math. Soc. 361 (2009), no. 1, 453–493. 10.1090/S0002-9947-08-04637-0Search in Google Scholar

[8] L. L. Fan and A. Matsumura, Asymptotic stability of a composite wave of two viscous shock waves for a one-dimensional system of non-viscous and heat-conductive ideal gas, J. Differr. Equ. 258 (2015), no. 4, 1129–1157. 10.1016/j.jde.2014.10.010Search in Google Scholar

[9] H. Grad, Asymptotic theory of the Boltzmann equation, Phys. Fluids 6 (1963), 147–181. 10.1063/1.1706716Search in Google Scholar

[10] H. Hong, Global stability of viscous contact wave for 1-D compressible Navier-Stokes equations, J. Differr. Equ. 252 (2012), no. 5, 3482–3505. 10.1016/j.jde.2011.11.015Search in Google Scholar

[11] B. K. Huang and Y. K. Liao, Global stability of combination of viscous contact wave with rarefaction wave for compressible Navier-Stokes equations with temperature-dependent viscosity, Math. Models Methods Appl. Sci. 27 (2017), no. 12, 2321–2379. 10.1142/S0218202517500464Search in Google Scholar

[12] F. M. Huang, J. Li, and A. Matsumura, Asymptotic stability of combination of viscous contact wave with rarefaction waves for one-dimensional compressible Navier-Stokes system, Arch. Ration. Mech. Anal. 197 (2010), no. 1, 89–116. 10.1007/s00205-009-0267-0Search in Google Scholar

[13] F. M. Huang and A. Matsumura, Stability of a composite wave of two viscous shock waves for the full compressible Navier-Stokes equation, Comm. Math. Phys. 289 (2009), no. 3, 841–861. 10.1007/s00220-009-0843-zSearch in Google Scholar

[14] F. M. Huang, A. Matsumura, and Z. P. Xin, Stability of contact discontinuities for the 1-D compressible Navier-Stokes equations, Arch. Ration. Mech. Anal. 179 (2006), no. 1, 55–77. 10.1007/s00205-005-0380-7Search in Google Scholar

[15] F. M. Huang, X. D. Shi, and Y. Wang, Stability of viscous shock wave for compressible Navier-Stokes equations with free boundary, Kinet. Relat. Models 3 (2010), no. 3, 409–425. 10.3934/krm.2010.3.409Search in Google Scholar

[16] B. K. Huang, S. J. Tang, and L. Zhang, Nonlinear stability of viscous shock profiles for compressible Navier-Stokes equations with temperature-dependent transport coefficients and large initial perturbation, Z. Angew. Math. Phys. 69 (2018), no. 6, Paper No. 136, 35pp. 10.1007/s00033-018-1026-6Search in Google Scholar

[17] F. M. Huang and T. Wang, Stability of superposition of viscous contact wave and rarefaction waves for compressible Navier-Stokes system, Indiana Univ. Math. J. 65 (2016), no. 6, 1833–1875. 10.1512/iumj.2016.65.5914Search in Google Scholar

[18] B. K. Huang, L. S. Wang, and Q. H. Xiao, Global nonlinear stability of rarefaction waves for compressible Navier-Stokes equations with temperature and density dependent transport coefficients, Kinet. Relat. Models 9 (2016), no. 3, 469–514. 10.3934/krm.2016004Search in Google Scholar

[19] F. M. Huang, Y. Wang, and T. Yang, Fluid dynamic limit to the Riemann solutions of Euler equations: I. Superposition of rarefaction waves and contact discontinuity, Kinet. Relat. Models 3 (2010), no. 4, 685–728. 10.3934/krm.2010.3.685Search in Google Scholar

[20] F. M. Huang, Z. P. Xin, and T. Yang, Contact discontinuity with general perturbations for gas motions, Adv. Math. 219 (2008), no. 4, 1246–1297. 10.1016/j.aim.2008.06.014Search in Google Scholar

[21] F. M. Huang and H. J. Zhao, On the global stability of contact discontinuity for compressible Navier-Stokes equations, Rend. Semin. Mat. Univ. Padova 109 (2003), 283–305. Search in Google Scholar

[22] Ya. I, Kanel, On a model system of equations of one-dimensional gas motion, Differ. Equ. 4 (1972), 374–380. Search in Google Scholar

[23] S. Kawashima and A. Matsumura, Asymptotic stability of traveling wave solutions of systems for one-dimensional gas motion, Comm. Math. Phys. 101 (1985), no. 1, 97–127. 10.1007/BF01212358Search in Google Scholar

[24] S. Kawashima, A. Matsumura, and K. Nishihara, Asymptotic behavior of solutions for the equations of a viscous heat-conductive gas, Proc. Japan Acad. Ser. A Math. Sci. 62 (1986), no. 7, 249–252. 10.3792/pjaa.62.249Search in Google Scholar

[25] S. Kawashima and M. Okada, Smooth global solutions for the one-dimensional equations in magnetohydrodynamics, Proc. Japan Acad. Ser. A Math. Sci. 58 (1982), no. 9, 384–387. 10.3792/pjaa.58.384Search in Google Scholar

[26] B. Kawohl, Global existence of large solutions to initial-boundary value problems for a viscous, heat-conducting, one-dimensional real gas, J. Differr. Equ. 58 (1985), no. 1, 76–103. 10.1016/0022-0396(85)90023-3Search in Google Scholar

[27] T. P. Liu and Z. P. Xin, Nonlinear stability of rarefaction waves for compressible Navier Stokes equations, Comm. Math. Phys. 118 (1988), no. 3, 451–465. 10.1007/BF01466726Search in Google Scholar

[28] T. P. Liu and Z. P. Xin, Pointwise decay to contact discontinuities for systems of viscous conservation laws, Asian J. Math. 1 (1997), no. 1, 34–84. 10.4310/AJM.1997.v1.n1.a3Search in Google Scholar

[29] H. X. Liu, T. Yang, H. J. Zhao, and Q. Y. Zou, One-dimensional compressible Navier-Stokes equations with temperature dependent transport coefficients and large data, SIAM J. Math. Anal. 46 (2014), no. 3, 2185–2228. 10.1137/130920617Search in Google Scholar

[30] A. Matsumura and K. Nishihara, On the stability of travelling wave solutions of a one-dimensional model system for compressible viscous gas, J. Appl. Math. 2 (1985), no. 1, 17–25. 10.1007/BF03167036Search in Google Scholar

[31] A. Matsumura and K. Nishihara, Asymptotics toward the rarefaction waves of the solutions of a one-dimensional model system for compressible viscous gas, Japan J. Appl. Math. 3 (1986), no. 1, 1–13. 10.1007/BF03167088Search in Google Scholar

[32] A. Matsumura and K. Nishihara, Global stability of the rarefaction wave of a one-dimensional model system for compressible viscous gas, Comm. Math. Phys. 144 (1992), no. 2, 325–335. 10.1007/BF02101095Search in Google Scholar

[33] T. Nishida and J. A. Smoller, Solutions in the large for some nonlinear hyperbolic conservation laws, Comm. Pure Appl. Math. 26 (1973), 183–200. 10.1002/cpa.3160260205Search in Google Scholar

[34] K. Nishihara, T. Yang, and H. J. Zhao, Nonlinear stability of strong rarefaction waves for compressible Navier-Stokes equations, SIAM J. Math. Anal. 35 (2004), no. 6, 1561–1597. 10.1137/S003614100342735XSearch in Google Scholar

[35] J. Smoller, Shock Waves and Reaction-Diffusion Equations, 2nd ed., Springer-Verlag New York, New York, 1994. 10.1007/978-1-4612-0873-0Search in Google Scholar

[36] Y. Sun, J. W. Zhang, and X. K. Zhao, Nonlinearly exponential stability for the compressible Navier-Stokes equations with temperature-dependent transport coefficients, J. Differr. Equ. 286 (2021), 676–709. 10.1016/j.jde.2021.03.044Search in Google Scholar

[37] C. J. van Duyn and L. A. Peletier, A class of similarity solutions of the nonlinear diffusion equation, Nonlinear Anal. 1 (1977), no. 3, 223–233. 10.1016/0362-546X(77)90032-3Search in Google Scholar

[38] T. Wang and H. J. Zhao, One-dimensional compressible heat-conducting gas with temperature-dependent viscosity, Math. Models Methods Appl. Sci. 26 (2016), no. 12, 2237–2275. 10.1142/S0218202516500524Search in Google Scholar

[39] Z. P. Xin, On nonlinear stability of contact discontinuities. Hyperbolic Problems: Theory, Numerics, Applications (Stony Brook, NY, 1994), World Scientific Publishing Co Pte Ltd, River Edge, NJ, 1996, pp. 249–257. Search in Google Scholar

[40] Y. Zel’dovich and Y. Rajzer, Physics of Shock Waves and High-temperature Hydrodynamic Phenomena, Academic Press, New York, 1967. Search in Google Scholar

[41] T. T. Zheng, Stability of a strong viscous contact discontinuity in a free boundary problem for compressible Navier-Stokes equations, Nonlinear Anal. Real World Appl. 25 (2015), 238–257. 10.1016/j.nonrwa.2015.03.001Search in Google Scholar

Received: 2021-08-21
Revised: 2022-02-13
Accepted: 2022-02-14
Published Online: 2022-08-26

© 2023 Wenchao Dong and Zhenhua Guo, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. On sufficient “local” conditions for existence results to generalized p(.)-Laplace equations involving critical growth
  3. On the critical Choquard-Kirchhoff problem on the Heisenberg group
  4. On the local behavior of local weak solutions to some singular anisotropic elliptic equations
  5. Viscosity method for random homogenization of fully nonlinear elliptic equations with highly oscillating obstacles
  6. Double-phase parabolic equations with variable growth and nonlinear sources
  7. Logistic damping effect in chemotaxis models with density-suppressed motility
  8. Bifurcation diagrams of one-dimensional Kirchhoff-type equations
  9. Standing wave solution for the generalized Jackiw-Pi model
  10. Weak-strong uniqueness and energy-variational solutions for a class of viscoelastoplastic fluid models
  11. Eigenvalue inequalities for the buckling problem of the drifting Laplacian of arbitrary order
  12. Homoclinic solutions for a differential inclusion system involving the p(t)-Laplacian
  13. Equivalence between a time-fractional and an integer-order gradient flow: The memory effect reflected in the energy
  14. Bautin bifurcation with additive noise
  15. Small solitons and multisolitons in the generalized Davey-Stewartson system
  16. Nonstationary Poiseuille flow of a non-Newtonian fluid with the shear rate-dependent viscosity
  17. A global compactness result with applications to a Hardy-Sobolev critical elliptic system involving coupled perturbation terms
  18. On a strongly damped semilinear wave equation with time-varying source and singular dissipation
  19. Singularity for a nonlinear degenerate hyperbolic-parabolic coupled system arising from nematic liquid crystals
  20. Stability of stationary solutions to the three-dimensional Navier-Stokes equations with surface tension
  21. Uniform decay estimates for the semi-linear wave equation with locally distributed mixed-type damping via arbitrary local viscoelastic versus frictional dissipative effects
  22. Symmetry and nonsymmetry of minimal action sign-changing solutions for the Choquard system
  23. Periodic and quasi-periodic solutions of a four-dimensional singular differential system describing the motion of vortices
  24. Multiple nontrivial solutions of superlinear fractional Laplace equations without (AR) condition
  25. Existence and blow-up of solutions in Hénon-type heat equation with exponential nonlinearity
  26. Existence and concentration behavior of positive solutions to Schrödinger-Poisson-Slater equations
  27. On the fractional Korn inequality in bounded domains: Counterexamples to the case ps < 1
  28. Gradient estimates for nonlinear elliptic equations involving the Witten Laplacian on smooth metric measure spaces and implications
  29. Dynamical analysis of a reaction–diffusion mosquito-borne model in a spatially heterogeneous environment
  30. Existence and multiplicity of solutions for a quasilinear system with locally superlinear condition
  31. Nontrivial solution for Klein-Gordon equation coupled with Born-Infeld theory with critical growth
  32. Modeling Wolbachia infection frequency in mosquito populations via a continuous periodic switching model
  33. Infinitely many localized semiclassical states for nonlinear Kirchhoff-type equation
  34. Fujita-type theorems for a quasilinear parabolic differential inequality with weighted nonlocal source term
  35. Approximations of center manifolds for delay stochastic differential equations with additive noise
  36. Periodic solutions to a class of distributed delay differential equations via variational methods
  37. Groundstate for the Schrödinger-Poisson-Slater equation involving the Coulomb-Sobolev critical exponent
  38. Existence of nontrivial solutions for the Klein-Gordon-Maxwell system with Berestycki-Lions conditions
  39. Global Sobolev regular solution for Boussinesq system
  40. Normalized solutions for the p-Laplacian equation with a trapping potential
  41. Nonlinear elliptic–parabolic problem involving p-Dirichlet-to-Neumann operator with critical exponent
  42. Blow-up for compressible Euler system with space-dependent damping in 1-D
  43. High energy solutions of general Kirchhoff type equations without the Ambrosetti-Rabinowitz type condition
  44. On the dynamics of grounded shallow ice sheets: Modeling and analysis
  45. A survey on some vanishing viscosity limit results
  46. Blow-up for logarithmic viscoelastic equations with delay and acoustic boundary conditions
  47. Generalized Liouville theorem for viscosity solutions to a singular Monge-Ampère equation
  48. Front propagation in a double degenerate equation with delay
  49. Positive solutions for a class of singular (pq)-equations
  50. Higher integrability for anisotropic parabolic systems of p-Laplace type
  51. The Poincaré map of degenerate monodromic singularities with Puiseux inverse integrating factor
  52. On a system of multi-component Ginzburg-Landau vortices
  53. Existence and concentration of solutions to Kirchhoff-type equations in ℝ2 with steep potential well vanishing at infinity and exponential critical nonlinearities
  54. Multiplicity results for fractional Schrödinger-Kirchhoff systems involving critical nonlinearities
  55. On double phase Kirchhoff problems with singular nonlinearity
  56. Estimates for eigenvalues of the Neumann and Steklov problems
  57. Nodal solutions with a prescribed number of nodes for the Kirchhoff-type problem with an asymptotically cubic term
  58. Dirichlet problems involving the Hardy-Leray operators with multiple polars
  59. Incompressible limit for compressible viscoelastic flows with large velocity
  60. Characterizations for the genuine Calderón-Zygmund operators and commutators on generalized Orlicz-Morrey spaces
  61. Non-local gradients in bounded domains motivated by continuum mechanics: Fundamental theorem of calculus and embeddings
  62. Noncoercive parabolic obstacle problems
  63. Touchdown solutions in general MEMS models
  64. Existence and nonexistence of nontrivial solutions for the Schrödinger-Poisson system with zero mass potential
  65. Short-time existence of a quasi-stationary fluid–structure interaction problem for plaque growth
  66. Fine bounds for best constants of fractional subcritical Sobolev embeddings and applications to nonlocal PDEs
  67. Symmetries of Ricci flows
  68. Asymptotic behavior of solutions of a second-order nonlinear discrete equation of Emden-Fowler type
  69. On the topological gradient method for an inverse problem resolution
  70. Supersolutions to nonautonomous Choquard equations in general domains
  71. Uniform complex time heat Kernel estimates without Gaussian bounds
  72. Global existence for time-dependent damped wave equations with nonlinear memory
  73. Normalized solutions for a critical fractional Choquard equation with a nonlocal perturbation
  74. Reversed Hardy-Littlewood-Sobolev inequalities with weights on the Heisenberg group
  75. Lamé system with weak damping and nonlinear time-varying delay
  76. Global well posedness for the semilinear edge-degenerate parabolic equations on singular manifolds
  77. Blowup in L1(Ω)-norm and global existence for time-fractional diffusion equations with polynomial semilinear terms
  78. Boundary regularity results for minimisers of convex functionals with (p, q)-growth
  79. Parametric singular double phase Dirichlet problems
  80. Special Issue on Nonlinear analysis: Perspectives and synergies
  81. Editorial to Special issue “Nonlinear analysis: Perspectives and synergies”
  82. Identification of discontinuous parameters in double phase obstacle problems
  83. Global boundedness to a 3D chemotaxis-Stokes system with porous medium cell diffusion and general sensitivity
  84. On regular solutions to compressible radiation hydrodynamic equations with far field vacuum
  85. On Cauchy problem for fractional parabolic-elliptic Keller-Segel model
  86. The evolution of immersed locally convex plane curves driven by anisotropic curvature flow
  87. Stability of combination of rarefaction waves with viscous contact wave for compressible Navier-Stokes equations with temperature-dependent transport coefficients and large data
  88. On concave perturbations of a periodic elliptic problem in R2 involving critical exponential growth
  89. Existence of solutions for a double-phase variable exponent equation without the Ambrosetti-Rabinowitz condition
Downloaded on 19.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/anona-2022-0246/html
Scroll to top button