Stability of strong viscous shock wave under periodic perturbation for 1-D isentropic navier-stokes system in the half space
Abstract.
In this paper, a viscous shock wave under space-periodic perturbation of 1-D isentropic Navier-Stokes system in the half space is investigated. It is shown that if the initial periodic perturbation around the viscous shock wave is small, then the solution time asymptotically tends to a viscous shock wave with a shift partially determined by the periodic oscillations. Moreover, the strength of the shock wave could be arbitrarily large. This result essentially improves the previous work ”A. Matsumura, M. Mei, Convergence to travelling fronts of solutions of the p-system with viscosity in the presence of a boundary. Arch. Ration. Mech. Anal. 146 (1999), no. 1, 1-22.” where the strength of shock wave is sufficiently small and the initial periodic oscillations vanish.
Key words and phrases:
Impermeable wall problem, Large amplitude shock, Space-periodic perturbation, Asymptotic stability1. Introduction
We consider a one-dimensional isentropic Navier-Stokes system for a general viscous gas, i.e.,
(1.1) |
where is the specific volume, the fluid velocity and is the pressure. Constant , are adiabatic constants. is the viscosity coefficient with . Without loss of generality, we assume in what follows.
The system (1.1) is a basic system of hydrodynamic equations, it has a variety of wave phenomena, such as viscous shock waves and rarefaction waves. So it is important to study the stability of the viscous shock wave for system (1.1). The stability of viscous shock wave for the Cauchy problem has been extensively studied in a large literature since the pioneer works of [5, 22], see the other interesting works [4, 7, 9, 12, 13, 16, 17, 18, 20, 23, 27, 15] as the shock wave is weak.
Physicists and engineers are more concerned with the stability of large amplitude shock (strong shock). However, the stability of large amplitude shock (strong shock) is challenging in mathematics. There have been no research results of this area until the last few years.
In 2010, Matsumura-Wang [25] proved that the large amplitude shock wave is asymptotically stable by a clever weighted energy method as . In 2016, Vasseur-Yao [28] successfully removed the condition by introducing a new variable called “effective velocity”. Recently, He-Huang [6] extended the result of [28] to general pressure and general viscosity , where could be any positive smooth function.
On the other hand, it is also interesting to investigate the stability of viscous shock waves for the initial-boundary value problem. In this paper, we considered an impermeable wall problem of (1.1) in the half space , i.e.,
(1.2) |
where . And are periodic functions with period and satisfy
(1.3) |
When the periodic functions vanish, Matsumura-Mei [21] considered the impermeable wall problem (1.1), (1.2) in 1999. And recently, an interesting result by [3] considering the multi-dimensional case of this problem.
The impermeable wall means that the velocity at the boundary must be zero because there is no flow across the boundary. They showed in [21] that when the solution of (1.1),(1.2) tends to a 2-viscous shock wave connecting the left state and the right one provided that both the strength of shock and the initial perturbation are small and the 2-viscous shock is initially far away from the boundary, where is determined by the RH condition, i.e.,
(1.6) |
Moreover, we assume that . The condition that the strength of shock is small was removed from [2]. The condition that the shock is initially far away from the boundary was removed from [24]. How to remove both these two conditions mentioned above at the same time is still open. Let us briefly recall the idea of [24].
Since at the boundary, we can exchange the impermeable wall problem (1.1) and (1.2) in the half space to the Cauchy problem in the whole space by defining as so that still satisfies the system (1.1) in the whole space, i.e.,
(1.7) |
equipped with the initial data
(1.8) |
satisfying
(1.9) |


It is obvious that the solution of the Cauchy problem (1.7)-(1.9) confined in the half line is exactly the one of the impermeable wall problem (1.1),(1.2). In view of the far field states at given by (1.9), it is expected that the solution to (1.7)-(1.9) asymptotically tends to a composite wave consisting of 1-viscous shock wave connecting at the left and an intermediate state at the right, and 2-viscous shock wave connecting at the left and at the right. Fortunately by the principle of RH condition and (1.9), see Fig (A), where means the 1-shock curve in the phase plane starting from the left state and means the 2-shock curve in the phase plane starting from the left state . The Figure B contains the graphs of the shock waves in the planes and . The wall can be regarded as a mirror and the 1-viscous shock is a mirror image of the 2-viscous shock in the plane , and the interaction between the 2-shock and the boundary for the impermeable wall problem (1.1)-(1.2) is replaced to consider the one between the 2-shock and its mirrored shock for the Cauchy problem (1.7)-(1.9).
In this paper, we want to improve the work of [2] where . Motivated by [24], the extended initial data in (1.8) satisfies
(1.10) |
We outline the strategy as follows. We apply the anti-derivative method to study the stability of the traveling wave solution , in which the anti-derivative of the perturbation (), namely, , “should” belong to some Sobolev spaces like . However, the method above can not be applicable directly in this paper since () oscillates at the far field and hence does not belong to any space for . Motivated by [30], we introduce a suitable ansatz , which has the same oscillations as the solution at the far field, so that belongs to some Sobolev spaces and the anti-derivative method is still available.
The rest of the paper will be arranged as follows. In Section 2, a suitable ansatz is constructed and the main results are stated. In Section 3, the stability problem is reformulated to a perturbation equation around the ansatz. In Section 4, the a priori estimates are established. In Section 5, the main results are proved. In Section 6, some complementary proofs are provided.
Notation. The functional defined by . When , the symbol is often omitted. As , we denote for simplicity,
In addition, denotes the -th order Sobolev space of functions defined by
2. Preliminaries and the Main Theorem
2.1. Preliminaries
As pointed out by [21, 2], when perturbation functions vanish, the solution of the impermeable wall problem (1.1)-(1.2) is expected to tend toward the outgoing viscous shock satisfying
(2.1) |
where , , is the shock speed determined by the R-H condition (1.6) and are given constants. Using and , it follows that
(2.2) |
Integrate (2.2) over . one has
(2.3) |
where . For abbreviation, we denote by . We have the following lemma.
Lemma 2.1 ([21]).
There exists a unique viscous shock up to a shift satisfying
(2.4) |
as where , ,
The initial data are assumed to satisfied
(2.5) |
and
(2.6) |
as compatibility condition, where is a constant. Set
We further assume that
(2.7) |
Borrowing from the idea of [24], we construct a composite wave. By [24], the mirrored shock , satisfies
(2.8) |
Thanks [21], one has
(2.9) |
The composite wave by two viscous shock weaves is defined as
(2.10) |
where is a constant. Motivated by [10, 8, 14, 11], we need two periodic solutions to (1.1) to establish the ansatz. Some properties of the solution are listed.
2.2. Ansatz
In order to make the anti-derivative method is available, we choose a suitable pair of ansatz such that for any . Motivated by [29], we define that the periodic solutions of (1.1) as for all which have the periodic initial data:
For the viscous shocks and define
(2.11) |
where we have used the R-H condition (1.6). It is straightforward to check that and for With functions in hand, we are ready to construct the ansatz. Let are two curves on which will be determined later. Set
Plugging the ansatz into (1.1), we have
(2.12) |
where
and
2.3. Location of The Shift and .
To apply the anti-derivative method which is always used to study the stability of viscous shock, introduced in [26], we expect that
When , the shifts and should satisfy
(2.13) |
Our next task is to show when . To make the system (2.12) as a conservative form, the curves and should satisfy
(2.14) |
With the aid of (2.2), we know , provided that the initial periodic perturbations are small. Due to (1.8) and (2.9), are odd functions and are even functions, thus , i.e, we can choose any to guarantee that . For , using (2.4) and (2.9), one gets that
(2.15) |
where
(2.16) |
By directly calculate, we have , .
(2.17) |
Moreover, choosing suitable small, we have Thus there exists a unique constant such that . Moreover, using , the constant is between and , where
(2.18) |
where we have used the following inequality
(2.19) |
By (2.5), we know exists. Thus we can obtain the curves and . More precisely, it holds that
Lemma 2.3.
Assume that (1.3), (1.6) hold. Then there exists an such that if
there exists a constant pair satisfying (2.13) where is uniquely determined and can take any constant. Moreover, there exists a unique solution to the system (2.14) with the fixed initial data ) satisfying
Moreover, the corresponding constant locations as follows,
(2.20) | ||||
and
(2.21) | ||||
where if ; , if .
2.4. the Main Result
We define
In view of (2.13), we further assume that
(2.22) |
Using the arbitrariness of , one can find a suitable , such that . From now on, we denote for simple.
Lemma 2.4.
Now, we turn to the original initial-value problem. Our main theorem is:
3. Reformulation of the Original Problem
Set
Thus satisfy
From , we know the ansazt satisfies
(3.1) |
where
(3.2) | ||||
Motivated by [21] and [1], with the help of (3.1) and (1.1), it follows that
(3.3) |
The initial condition satisfies
(3.4) |
where
(3.5) |
(3.6) |
Lemma 3.1.
We will seek the solution in the functional space for any ,
where is small.
Remark 3.1.
The function space is well defined because the Dirac function will not appear in , which can be guaranteed by .
Proposition 3.1.
Once Proposition 3.1 is obtained, the local solution can be extend to See the following lemma.
4. A Priori Estimate
For some , the problem is assumed that has a solution in this section.
(4.1) |
The Sobolev inequality gives that , and
Motivated by [28], we introduce the new effective velocity . It holds that
(4.2) |
Similarly, we define , then (3.1) becomes
(4.3) |
We define
(4.4) |
Substitute (4.3) from (4.2) and integrate the resulting system with respect to . Using (4.4), we have
(4.5) |
where
Now we give some lemmas that are useful in energy estimate.
Lemma 4.2.
The error terms
(4.8) |
satisfy
Proof.
4.1. Low Order Estimates.
Lemma 4.3.
Under the same assumptions of Proposition 3.1, we have
Proof.
We multiply and by and , respectively, sum them up, and intergrading result with respect to and over , we have
(4.13) | ||||
(4.14) |
By direct calculate, one gets that
(4.15) |
Due to (4.6), (4.15), we can get
(4.16) |
With the aid of Lemma 3.1, Hlder inequality, we have
(4.17) |
Using Hlder inequality, Sobolev inequality, combining Lemma 3.1, Lemma 4.2, one gets
(4.18) |
Inserting (4.16)-(4.18) into (4.13), using the smallness of , we obtain the proof of Lemma 4.3. ∎
Lemma 4.4.
Under the same assumptions of Proposition 3.1, we have
Proof.
We multiply and by and , respectively and sum over the result, intergrade the result with respect to and over , we have
(4.19) | ||||
With the aid of the Cauchy inequality, we have
(4.20) |
The last inequality is based on the following inequality
where we have used , , .
Lemma 4.5.
Under the same assumptions of Proposition 3.1, we have
Proof.
We multiply by and make use of , we get
(4.22) |
Intergrade with respect to and over , we have
We estimate term by term. By the Cauchy inequality, it follows that
(4.23) |
In addition, it is straightforward to imply that
(4.24) |
(4.25) |
4.2. High Order Estimates.
If we continue to get the estimates of second order derivative , new difficulties arise. In fact, in order to close the a priori estimate, should be sufficiently small. Unfortunately, it means that we have to add an additional condition “” which can guarantee that the Dirac function will not appear. Next, we need change variables to .
Lemma 4.6.
Under the same assumptions of Proposition 3.1, for , it holds that:
Proof.
This lemma is similar like [1] and the proof is omitted. ∎
Using this lemma, low order estimate (4.27) can be rewritten as
Lemma 4.7.
Under the same assumptions of Proposition 3.1, it holds that
Next, we turn to the original equation (3.3) to study the higher order estimates.
Lemma 4.8.
Under the same assumptions of Proposition 3.1, it holds that
(4.28) |
Proof.
Lemma 4.9.
Under the same assumptions of Proposition 3.1, it holds that
(4.33) |
Proof.
Differentiating with respect to , using , we have
(4.34) |
Differentiating in respect of and multiplying the derivative by , integrating the result in respect of and over , one has
(4.35) | ||||
By , one has
(4.36) |
The Cauchy inequality yields
(4.37) |
Similar to (4.30), we get
(4.38) |
can be controlled by (4.28). Using , and Cauchy inequality, we have
The Cauchy inequality yields
(4.39) |
With the help of
one gets
(4.40) |
Similar like (4.18), one gets that
(4.41) |
On the other hand, differentiating the second equation of (3.3) with respect to , multiplying the derivative by , integrating the resulting equality over , using Lemma 4.7 - Lemma 4.9, we can get the highest order estimate in the same way, which is listed as follows and the proof is omitted.
Lemma 4.10.
Under the same assumptions of Proposition 3.1, it holds that
5. Proof of Theorem 2.1
It is straightforward to imply (2.23) from Lemma 3.2. It remains to show (2.24). The following useful lemma will be used.
Lemma 5.1.
([22]) Assume that the function , then it holds that as .
Let us turn to the system (3.3). Differentiating (3.3)1 with respect to , multiplying the resulting equation by and integrating it with respect to on , we have
With the aid of Lemma 3.2, we have
which implies . By Lemma 5.1, we have
Since is bounded, the Sobolev inequality implies that
Similarly, we have
Therefore, the proof of Lemma 2.4 is completed.
5.1. Proof of Theorem 2.1
Proof.
With the aid of (see Lemma 2.1) and , it follows that . Thus if and , using (2.18), we obtain
Similar, with the help of (2.20), we have
Thus, it follows that
Set
(5.1) |
Make full use of (2.4), when , we have
Thus, we have
where is independent of and . Similarly, we can prove that and . Thus, we proved . In the same way, we have that
6. Proof of Lemma 2.3 and Lemma 3.1
6.1. Proof of Lemma 2.3
Proof.
By Lemma 2.2, we have for all . Thus and are all exist. In the following part of this subsection, we compute the two limits. Motivated by [11], we define the domain
(6.1) |
where . Using , we have
With the aid of Green formula, one gets
(6.2) |
where
We rewrite as:
where
Here
Moreover, can be rewrite as
Since , then
So we obtain
(6.3) |
where we have used (2.9),(2.11)in the last equality. With the aid of Lemma 2.2, one gets that
(6.4) |
By directly calculate, we have
Using (2.8) (2.9), one gets that
(6.5) |
The integral on in (6.2) satisfies that
(6.6) |
Here we have used Since as . By same method, we obtain
(6.7) |
Collecting (6.2)-(6.7), it follows that
Thus we obtain (2.20) where we have used R-H conditions (1.6)1. We omit the proof of (2.21), since it is similar with (2.20). ∎
6.2. Proof of Lemma 3.1
We only give the proof of , due to the fact that the proof of is similar.
Case 1. For we rewrite as follows.
where
It follows that
(6.8) |
Here we have used Lemma 2.2, (2.4), (2.9), (2.11). Moreover, can be rewritten as as follows.
(6.9) |
where
(6.10) |
and
(6.11) |
By , we get
(6.12) |
On the other hand, in the same way, it is still true to replace () with () in (6.12). We get If we choose sufficiently large, for , it follows that:
where we have used Lemma 2.1 in the second inequality and Lemma 5.2 in the last inequality. Thus, we obtain that
Similar like (6.8), one can get that
(6.13) |
Case 2. If , using (2.14), one can decompose as
(6.14) |
(6.15) |
and using similar arguments as in the case 1 to obtain that
(6.16) |
Remark 6.1.
If (isothermal gas) in our equations, we can get the same result by the same method.
Remark 6.2.
In our proof, we make the position of the shock is far away from the wall, is this necessary?
References
- [1] L. Chang, Stability of a Composite Wave of Two Seperate Strong Viscous Shock Waves for 1-D Isentropic Navier-Stokes System. arXiv:2103.15133, (2021),1-16.
- [2] L. Chang, Stability of Large Amplitude Viscous Shock Wave for 1-D Isentropic Navier-Stokes System in the Half Space. arXiv:2103.15133, (2021),1-14.
- [3] L. Chang, L. Liu, L. Xu, Nonlinear stability of planar shock wave to 3-D compressible Navier-Stokes equations in half space with Navier Boundary conditions. arXiv:2312.05565.
- [4] H. Freistuhler, D. Serre, stability of shock waves in scalar viscous conservation laws. Comm. Pure Appl. Math. 51 (1998), no. 3, 291-301.
- [5] J. Goodman, Nonlinear asymptotic stability of viscous shock profiles for conservation laws. Arch. Rational Mech. Anal. 95 (1986), no. 4, 325-344.
- [6] L. He, F. Huang, Nonlinear stability of large amplitude viscous shock wave for general viscous gas. J. Differential Equations 269 (2020), no. 2, 1226-1242.
- [7] F. Huang, A. Matsumura, Stability of a composite wave of two viscous shock waves for full compressible Navier-Stokes equation. Comm. Math. Phys. 289 (2009), no. 3, 841-861.
- [8] F. Huang, Z. Xin, L. Xu, Q. Yuan, Nonlinear asymptotic stability of compressible vortex sheets with viscosity effects. arXiv:2308.06180
- [9] F. Huang, L. Xu, Decay rate toward the traveling wave for scalar viscous conservation law. Commun. Math. Anal. Appl. 1, No. 3, 395-409 (2022).
- [10] F. Huang, L. Xu, Q. Yuan, Asymptotic stability of planar rarefaction waves under periodic perturbations for 3-d Navier-Stokes equations. Adv. Math. 404, Part B, Article ID 108452, 27 p. (2022).
- [11] F. Huang, Q. Yuan, Stability of large-amplitude viscous shock under periodic perturbation for 1-d isentropic Navier-Stokes equations, Commun. Math. Phys. 387 (2021), 1655-1679.
- [12] J. Humpherys, G. Lyng, K. Zumbrun, Multidimensional stability of large-amplitude Navier-Stokes shocks. Arch. Ration. Mech. Anal. 226 (2017), no. 3, 923-973.
- [13] S. Kawashima, A. Matsumura, Asymptotic stability of traveling wave solutions of systems for one-dimensional gas motion. Commun. Math. Phys. 101. (1985), no. 1, 97-127.
- [14] L. Liu, D. Wang, L. Xu, Asymptotic stability of the combination of a viscous contact wave with two rarefaction waves for 1-D Navier-Stokes equations under periodic perturbations. J. Differ. Equations 346, 254-276 (2023).
- [15] L. Liu, S. Wang, L. Xu, Decay rate to the planar viscous shock wave for multi-dimensional scalar conservation laws. arXiv:2312.03553.
- [16] T. Liu, Pointwise convergence to shock waves for viscous conservation laws. Comm. Pure Appl. Math. 50 (1997), no. 11, 1113-1182.
- [17] T. Liu, Y. Zeng, Time-asymptotic behavior of wave propagation around a viscous shock profile. Comm. Math. Phys. 290 (2009), no. 1, 23-82.
- [18] T. Liu, Y. Zeng, Shock waves in conservation laws with physical viscosity. Mem. Amer. Math. Soc. 234, (2015), no. 1105.
- [19] C. Mascia, K. Zumbrun, Stability of large-amplitude viscous shock profiles of hyperbolic-parabolic system. Arch. Ration. Mech. Anal. 172 (2004), no. 1, 93-131.
- [20] A. Matsumura, Waves in compressible fluids: viscous shock, rarefaction, and contact waves. Handbook of mathematical analysis in mechanics of viscous fluids. 2495-2548, Springer, Cham, 2018.
- [21] A. Matsumura, M. Mei, Convergence to travelling fronts of solutions of the -system with viscosity in the presence of a boundary.Arch. Ration. Mech. Anal. 146 (1999), no. 1, 1-22.
- [22] A. Matsumura, K. Nishihara, On the stability of travelling wave solutions of a one-dimensional model system for compressible viscous gas. Japan. J. Appl. Math. 2 (1985), no. 1, 17-25.
- [23] A. Matsumura, K. Nishihara, Asymptotic stability of traveling waves for scalar viscous conservation laws with non-convex nonlinearity. Comm. Math. Phys. 165 (1994), no. 1, 83-96.
- [24] A. Matsumura, K. Nishihara, Global Solutions for Nonlinear Differential Equations-Mathematical Analysis on Compressible Viscous Fluids (In Japanese). Nippon Hyoronsha, 2004.
- [25] A. Matsumura, Y. Wang, Asymptotic stability of viscous shock wave for a one-dimensional isentropic model of viscous gas with density dependent viscosity. Methods Appl. Anal. 17 (2010), no. 3, 279-290.
- [26] A. M. Il’in, O. A. Oleinik, Asymptotic behavior of solutions of the Cauchy problem for some quasi-linear equations for large values of the time, Mat. Sb. (N.S.) 51(93), (1960), no. 2, 191–216
- [27] A. Szepessy, Z. Xin, Nonlinear stability of viscous shock waves. Arch. Ration. Mech. Anal. 122 (1993), no. 1, 53-103.
- [28] A. Vasseur, L. Yao, Nonlinear stability of viscous shock wave to one-dimensional compressible isentropic Navier-Stokes equations with density dependent viscous coefficient. Commun. Math. Sci. 14 (2016), no. 8, 2215-2228.
- [29] Z. Xin, Q. Yuan, and Y. Yuan, Asymptotic stability of shock profiles and rarefaction waves under periodic perturbations for 1-D convex scalar viscous conservation laws, Indiana Univ. Math. J. 70 (2021), no. 6, 2295-2349.
- [30] Z. Xin, Q. Yuan, Y. Yuan, Asymptotic stability of shock waves and rarefaction waves under periodic perturbations for 1-D convex scalar conservation laws. SIAM J. Math. Anal. 51 (2019), no. 4, 2971–2994.
- [31] K. Zumbrun, Stability of large-amplitude shock waves of compressible Navier-Stokes equations. with an appendix by Helge Kristian Jenssen and Gregory Lyng. Handbook