REVIEW 3 major objections 6 minor 37 references
Asymptotic behavior toward viscous shocks for the outflow problem of barotropic Navier-Stokes equations
T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Weak viscous shocks are stable in the outflow problem for barotropic Navier-Stokes equations on a half line.
desk verdict First viscous-shock stability result for the outflow problem with a genuine Eulerian a-contraction adaptation, but the shift ODE (3.3) doesn't match the Y_g used in the proof and the referee must get that fixed. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the $a$-contraction with a dynamical shift $X(t)$: a weight function $a(\xi)=1+(\sqrt{\delta}+u_+-\tilde u(\xi))/\sqrt{\delta}$ multiplies the relative entropy, and the shift $X(t)$ solves the ODE (3.3). The weight satisfies $a'(\xi)>0$, so the monotone shock profile creates good dissipation terms, while the ODE is chosen so that the leading-order bad term $\dot X\,Y$ cancels. A change of variables $y=(u_- - \tilde u)/\delta$ maps the half line to an interval $(y_0(t),1)$ with $y_0(t)<1/6$ for large $β$, and a weighted Poincaré inequality supplies the sharp $L^2$ contraction. The boundary terms become exponentially small in $δ\u03b2$, which is what lets the contraction estimate close without a density boundary condition.
What would settle it
Set up the outflow problem with $γ=2$, $p(ρ)=ρ^2$, a right state $(ρ_+,u_+)$ on the transonic curve $u_++\sqrt{p'(ρ_+)}=0$, and $u_-$ slightly larger so that $δ=|u_--u_+|$ is small. Place the initial shock at $β=0$, so that $y_0(0)=-\tilde u(-\u03b2)/\delta$ is large, and take initial data equal to the truncated shock profile plus an $H^1$ perturbation smaller than $ε_0$. Measure $E(t)=\sup_{x\in\mathbb{R}_+}|(ρ,u)(t,x)-(\tildeρ,\tilde u)(x-σt-X(t)-β)|$; Theorem 1.1 predicts $E(t)\to0$, whereas the proof's boundary control, Lemma 4.7 requiring $y_0<1/6$, would not apply. A failure of $E(t)\to0$ in this configuration would falsify the claim as stated; success would show the boundary-distance hypothesis is stronger than necessary.
Extended reading notes
Core claim
The central claim is Theorem 1.1: for a right state $(ρ_+,u_+)$ in the subsonic or transonic region, with $u_-<0$, $u_->u_+$, and small shock strength $δ=|u_--u_+|$, the viscous 2-shock profile is asymptotically stable for the outflow problem. For initial data close to the shock profile in $H^1$ and with the shock initially placed at a sufficiently large distance $β$ from the boundary, there is a unique global strong solution, and $\sup_{x\in\mathbb{R}_+}|(ρ,u)(t,x)-(\tildeρ,\tilde u)(x-σt-X(t)-β)|\to0$ as $t\to\infty$, while $|\dot X(t)|\to0$. The proof establishes this by controlling the $L^2$ perturbation through a weighted relative entropy estimate, controlling boundary terms through the large distance $β$, and closing the argument with higher-order energy estimates.
Load-bearing premise
The proof assumes the dynamical shift stays bounded by half the shock speed, $|X(t)|\le(\sigma/2)t$, so the shock never approaches the outflow boundary; if the shift drifts faster, the boundary terms lose their exponential smallness and the $a$-contraction estimate does not close.
Editorial extensions
If this is right
- The outflow problem admits a unique global strong solution near a weak outgoing viscous 2-shock when the right state is subsonic or transonic and the shock starts far enough from the boundary.
- The solution converges in supremum norm to the shifted viscous shock, so density and velocity each become pointwise indistinguishable from a travelling shock profile whose location is set by $X(t)$.
- The shift speed $|\dot X(t)|$ decays to zero, so time-asymptotically the shifted wave is again a genuine viscous shock profile moving at the original speed $σ$.
- The Eulerian $a$-contraction argument bypasses the free-boundary formulation caused by the missing density boundary condition, opening a route for other half-line problems with no density data.
Reading between the lines
- The small-shock restriction appears tied to the geometry of the 2-Hugoniot curve: a supersonic right state cannot reach a subsonic left state with a small jump, so a large-shock version of the $a$-contraction machinery would be needed for that regime.
- The proof's exponential boundary errors $e^{-C\delta\beta}$ suggest that the initial distance $β$ controls the size of boundary effects; a quantitative convergence rate in $t$ may be extractable, although the paper only proves qualitative convergence.
- The same weighted relative entropy with a shift could be tried for composite waves involving a viscous shock and a boundary layer, where the boundary term and the shift would interact; this is not addressed in the paper.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the half-line outflow problem for the 1D barotropic Navier-Stokes system with boundary condition u(t,0)=u_-<0 and a prescribed right state (rho_+,u_+) in the subsonic or transonic region. It claims that a small-amplitude 2-viscous shock profile, initially placed far from the boundary, is nonlinearly stable under small H^1 perturbations of the initial data: the solution exists globally and converges in L^∞ to the shock profile with a time-dependent shift X(t), and |dot X(t)| tends to 0. The proof uses the a-contraction method with shift in Eulerian coordinates, a relative entropy identity, a Poincaré-type inequality in the shock-adapted variable y, higher-order estimates, and a continuation argument. Boundary terms are controlled by choosing the initial shock position beta so large that the shifted shock tail at x=0 is exponentially small over the relevant time interval.
Significance. If the theorem is correct, it is the first asymptotic stability result for viscous shocks in the outflow problem, and it fills a genuine gap left by the anti-derivative approach, which cannot handle the missing density boundary condition. The adaptation of the a-contraction framework to the Eulerian half-line problem is nontrivial, and the boundary treatment through the drift bound |X(t)|≤σt/2 is a natural and promising idea. The paper is carefully structured: the relative entropy computation is explicit, the dissipation terms Gnew, GS, Gbd, Dρ, Du1, Du2 are clearly identified, and the higher-order estimates are presented in detail. The manuscript also gives reproducible algebraic details in Lemmas 4.3, 4.5, and 4.6, which is a strength.
major comments (3)
- [Section 3.2, Eq. (3.3); Section 4.3, Lemma 4.5] The shift ODE (3.3) as printed is inconsistent with the identity dot X = -M Y_g/delta used in Lemma 4.5. In (3.3) the first integrand is a p(tilde_rho)/((sigma-tilde_u)tilde_rho) tilde_rho_x (u-tilde_u), whereas Y_g in Lemma 4.5 is defined with a p'(tilde_rho)/(sigma-tilde_u) tilde_rho_x (u-tilde_u). Since p'(tilde_rho)=gamma p(tilde_rho)/tilde_rho, the printed ODE gives dot X = -(M/delta)(Y_g - (1-1/gamma)A) with A = integral a p'(tilde_rho)/(sigma-tilde_u) tilde_rho_x (u-tilde_u) dx. The exact cancellation that produces dot X Y = -delta |dot X|^2/M + dot X Y_1 + dot X Y_2, and hence the entire negativity argument including the choice M=2(gamma+1)/rho_+ and the constant (gamma+1)rho_+ G_S/8, depends on this identity. If (3.3) is a typo, the corrected ODE and its relation to Y_g must be stated; if it is not a typo, Lemma 4.5 does not follow. This is a load-bearing internal consistency issue in the central a-contraction estimate.
- [Section 3.1 and Appendix A.1] The global existence proof is not fully self-contained. Proposition 3.1 is stated without proof, and the continuation argument in Appendix A.1 relies on the existence and Lipschitz regularity of the shift X(t) solving (3.3), which is cited from [27]. These inputs are load-bearing: the drift bound |X(t)|≤σt/2, obtained from the a priori bound on |dot X| under assumption (3.4), is used to ensure (4.9), the exponential smallness of y0(t), and the boundary estimates in Lemma 4.7. Since the outflow boundary condition does not prescribe the density at x=0, the local well-posedness is not an entirely standard textbook case. The authors should either supply the proof of Proposition 3.1 and of the ODE well-posedness in this boundary setting, or state precisely which theorem in [27] applies and verify that its hypotheses are satisfied for the shifted shock on the half-line.
- [Appendix A.1] The construction of the constants epsilon_0 and epsilon_* needs clarification. The text defines epsilon_0 = epsilon_* - C delta and then states that epsilon_0 can be chosen independently of delta, for example epsilon_0 = epsilon/(4(C0+1)). This requires showing that, for every delta<delta_0, the corresponding epsilon_* (which contains C sqrt(delta)+e^{-C delta beta}) is larger than epsilon/(4(C0+1)). Since beta is chosen after delta and can be taken arbitrarily large, the claim is plausible, but the uniform-in-delta argument should be written out explicitly.
minor comments (6)
- [Proposition 3.2, Eq. (3.5)] The square root in front of the time-integrated dissipation terms appears to be a typo; Lemma 5.2 and the applications later in the paper use the same quantity without a square root.
- [Lemma 5.2] The heading 'Under the hypothesis in (3.2)' should read 'Under the hypothesis of Proposition 3.2'.
- [Appendix A.1] The text refers to 'constants defined in Theorem 3.2', but the statement is Proposition 3.2, not a theorem.
- [Appendix A.3] In the proof of (1.11), the limit should be 'as t tends to infinity', not 'as t tends to 0'.
- [Proposition 3.1] The notation kappa_0, kappa_0, kappa_1, kappa_1 repeats the same symbols; distinct symbols (for example, underlined and overlined versions) should be used for the lower and upper bounds.
- [References [29], [30]] The author name 'Shinya, N.' should be 'Nishibata, S.' in both entries; the in-text citation to Kawashima-Nishibata-Zhu is consistent with the correct spelling.
Circularity Check
No significant circularity: the outflow stability theorem is proved from explicit a priori estimates rather than assumed or fitted, and the cited prior work is methodological.
full rationale
The central claim is not an input to the derivation. Theorem 1.1 is established through the a priori estimate Proposition 3.2, whose proof is carried out in Sections 4 and 5: the relative entropy identity (Lemma 4.2), the bad/good term decomposition (Lemmas 4.3 and 4.4), the shift estimate (Lemma 4.5), and the boundary estimate (Lemma 4.7) are all derived in the paper. The shift X(t) is defined by the explicit ODE (3.3) and is not fitted to the asymptotic state; the decay of |dot X| is proved from the a priori bounds in Appendix A.3. Citations to prior work of the authors, such as [20], [23], [24], and [27], supply the a-contraction technique and standard profile or Poincare-type facts, but they do not contain the outflow-shock stability theorem and are not used as a substitute for the target conclusion. In particular, the bound |X(t)| <= (sigma/2)t used in Lemma 4.5 is justified in the text from smallness of epsilon rather than imported as an assumption. A separate, non-circular correctness concern is that the printed shift ODE (3.3) contains p(rho_tilde)/rho_tilde while Lemma 4.5's Y_g contains p'(rho_tilde), so the statement 'Notice from (3.3) that dot X = -M Y_g/delta' is not algebraically forced as printed; that is a consistency or typo issue rather than circularity, because neither formulation assumes the theorem's conclusion.
Assumptions & free parameters
free parameters (3)
- beta (shock offset)
- delta_0 (max shock strength)
- epsilon_0 (initial perturbation size)
assumptions (5)
- domain assumption Viscous 2-shock profile exists uniquely and is strictly monotone with exponential decay (Lemma 2.1), quoted from [24].
- standard math Poincare-type inequality on bounded intervals (Lemma 2.2), quoted from [20].
- domain assumption a-contraction method and shift ODE well-posedness from [23,24,27], including the a priori bound |X(t)| <= C epsilon t.
- standard math Rankine-Hugoniot and Lax entropy conditions for the 2-shock, and the shock speed formula sigma > 0.
- standard math Sobolev embedding, interpolation, and local existence via standard iteration (Proposition 3.1).
Cite this review
Pith. "Pith review of Asymptotic behavior toward viscous shocks for the outflow problem of barotropic Navier-Stokes equations." pith.science (2026). https://pith.science/paper/H5WZW7ZD
@misc{pith2026250508171,
author = {Pith},
title = {Pith review of: Asymptotic behavior toward viscous shocks for the outflow problem of barotropic Navier-Stokes equations},
year = {2026},
howpublished = {\url{https://pith.science/paper/H5WZW7ZD}},
note = {Machine review of arXiv:2505.08171}
}
abstract
We study the large-time asymptotic stability of viscous shock profile to the outflow problem of barotropic Navier-Stokes equations on a half line. We consider the case when the far-field state as a right-end state of 2-Hugoniot shock curve belongs to the subsonic region or transonic curve. We employ the method of $a$-contraction with shifts, to prove that if the strength of viscous shock wave is small and sufficiently away from the boundary, and if a initial perturbation is small, then the solution asymptotically converges to the viscous shock up to a dynamical shift. We also prove that the speed of time-dependent shift decays to zero as times goes to infinity, the shifted viscous shock still retains its original profile time-asymptotically. Since the outflow problem in the Lagrangian mass coordinate leads to a free boundary value problem due to the absence of a boundary condition for the fluid density, we consider the problem in the Eulerian coordinate instead. Although the $a$-contraction method is technically more complicated in the Eulerian coordinate than in the Lagrangian one, this provides a more favorable framework by avoiding the difficulty arising from a free boundary.
Reference graph
Works this paper leans on
-
[20]
Huang, X., Kang, M.-J., Kim, J., and Lee, H. Asymptotic behavior toward viscous shock for imper- meable wall and inflow problem of barotropic Navier-Stokes e quations. arXiv preprint arXiv:2405.03214 (2024). 32 KANG, OH, AND W ANG
work page Pith review arXiv 2024
-
[27]
Kang, M.-J., V asseur, A. F., and W ang, Y. Time-asymptotic stability of composite waves of viscous shock and rarefaction for barotropic Navier-Stokes equati ons. Adv. Math. 419 (2023), Paper No. 108963, 66
work page 2023
-
[1]
Initial-boundary value problems for nonlinear systems of c onservation laws
Amadori, D. Initial-boundary value problems for nonlinear systems of c onservation laws. NoDEA Nonlinear Differential Equations Appl. 4 , 1 (1997), 1–42
work page 1997
-
[2]
Amadori, D., and Colombo, R. M. Continuous dependence for 2 × 2 conservation laws with boundary. J. Differential Equations 138 , 2 (1997), 229–266
work page 1997
-
[3]
Ancona, F., Marson, A., and Spinolo, L. V. Existence of vanishing physical viscosity solutions of characteristic initial-boundary value problems for sys tems of conservation laws. arXiv preprint arXiv:2401.14865 (2024)
work page Pith review arXiv 2024
-
[4]
Bardos, C., le Roux, A. Y., and N ´ed´elec, J.-C. First order quasilinear equations with boundary conditions. Comm. Partial Differential Equations 4 , 9 (1979), 1017–1034
work page 1979
-
[5]
Bianchini, S., and Spinolo, L. V. The boundary Riemann solver coming from the real vanishing viscosity approximation. Arch. Ration. Mech. Anal. 191 , 1 (2009), 1–96
work page 2009
-
[6]
Bianchini, S., and Spinolo, L. V. Characteristic boundary layers for mixed hyperbolic-para bolic systems in one space dimension and applications to the Navie r-Stokes and MHD equations. Comm. Pure Appl. Math. 73 , 10 (2020), 2180–2247
work page 2020
Show all 37 references
-
[7]
From Navier-Stokes to BV solutions of the barotropic Euler equations
Chen, G., Kang, M.-J., and V asseur, A. From Navier-Stokes to BV solutions of the barotropic Euler equations. arXiv preprint arXiv:2401.09305 (2024)
2024 arXiv
-
[8]
Spectral stability of noncharac- teristic isentropic Navier-Stokes boundary layers
Costanzino, N., Humpherys, J., Nguyen, T., and Zumbrun, K. Spectral stability of noncharac- teristic isentropic Navier-Stokes boundary layers. Arch. Ration. Mech. Anal. 192 , 3 (2009), 537–587
2009
-
[9]
Dafermos, C. M. Entropy and the stability of classical solutions of hyperbo lic systems of conserva- tion laws. In Recent mathematical methods in nonlinear wave propagation (Montecatini Terme, 1994) , vol. 1640 of Lecture Notes in Math. Springer, Berlin, 1996, pp. 48–69
1994
-
[10]
DiPerna, R. J. Uniqueness of solutions to hyperbolic conservation laws. Indiana Univ. Math. J. 28 , 1 (1979), 137–188
1979
-
[11]
c., and LeFloch, P
Dubois, F. c., and LeFloch, P. Boundary conditions for nonlinear hyperbolic systems of co nservation laws. J. Differential Equations 71 , 1 (1988), 93–122
1988
-
[12]
Stability of a Riemann shock in a physical class: From Brenne r- Navier-Stokes-Fourier to Euler
Eo, S., Eun, N., and Kang, M.-J. Stability of a Riemann shock in a physical class: From Brenne r- Navier-Stokes-Fourier to Euler. arXiv preprint arXiv:2411.03613 (2024)
2024 arXiv
-
[13]
Existence and stability of noncharac- teristic boundary layers for the compressible Navier-Stok es and viscous MHD equations
Gu`es, O., M ´etivier, G., Williams, M., and Zumbrun, K. Existence and stability of noncharac- teristic boundary layers for the compressible Navier-Stok es and viscous MHD equations. Arch. Ration. Mech. Anal. 197 , 1 (2010), 1–87
2010
-
[14]
Large-time behavior of composite waves of viscous shocks fo r the barotropic Navier-Stokes equations
Han, S., Kang, M.-J., and Kim, J. Large-time behavior of composite waves of viscous shocks fo r the barotropic Navier-Stokes equations. SIAM J. Math. Anal. 55 , 5 (2023), 5526–5574
2023
-
[15]
Convergence to superposition of boundary layer, rarefaction and shock for the 1D Navier-Stokes equations
Han, S., Kang, M.-J., Kim, J., Kim, N., and Oh, H. Convergence to superposition of boundary layer, rarefaction and shock for the 1D Navier-Stokes equations. arXiv preprint arXiv:2502.08126 (2025)
2025
-
[16]
Long-time behavior towards viscous-dispersive shock for Navier-Stokes equations of Korteweg type
Han, S., Kang, M.-J., Kim, J., and Lee, H. Long-time behavior towards viscous-dispersive shock for Navier-Stokes equations of Korteweg type. J. Differential Equations 426 (2025), 317–387
2025
-
[17]
Inviscid limit problem of radially symmetric stationary so lutions for compressible Navier-Stokes equation
Hashimoto, I., and Matsumura, A. Inviscid limit problem of radially symmetric stationary so lutions for compressible Navier-Stokes equation. Comm. Math. Phys. 405 , 9 (2024), Paper No. 220, 20
2024
-
[18]
Viscous shock wave and boundary layer solution to an inflow problem for compressible viscous gas
Huang, F., Matsumura, A., and Shi, X. Viscous shock wave and boundary layer solution to an inflow problem for compressible viscous gas. Comm. Math. Phys. 239 , 1-2 (2003), 261–285
2003
-
[19]
Uniqueness of composite wave of shock and rarefaction in the inviscid limit of Navier-Stokes equatio ns
Huang, F., W ang, W., W ang, Y., and W ang, Y. Uniqueness of composite wave of shock and rarefaction in the inviscid limit of Navier-Stokes equatio ns. SIAM Math. Anal. 56 (3) (2024), 3924– 3967
2024
-
[21]
Time-asymptotic stability of composite wave of viscous sho cks and viscous contact wave for Navier-Stokes-Fourier equations
Huang, X., and Lee, H. Time-asymptotic stability of composite wave of viscous sho cks and viscous contact wave for Navier-Stokes-Fourier equations. arXiv preprint arXiv:2504.04014 (2025)
2025
-
[22]
T., Santos, M
Kan, P. T., Santos, M. M., and Xin, Z. Initial-boundary value problem for conservation laws. Comm. Math. Phys. 186 , 3 (1997), 701–730
1997
-
[23]
Kang, M.-J., and V asseur, A. F. L2-contraction for shock waves of scalar viscous conservatio n laws. Ann. Inst. H. Poincar´ e C Anal. Non Lin´ eaire 34, 1 (2017), 139–156
2017
-
[24]
Kang, M.-J., and V asseur, A. F. Contraction property for large perturbations of shocks of t he barotropic Navier-Stokes system. J. Eur. Math. Soc. 23 , 2 (2021), 585–638
2021
-
[25]
Kang, M.-J., and V asseur, A. F. Uniqueness and stability of entropy shocks to the isentropi c Euler system in a class of inviscid limits from a large family of Nav ier-Stokes systems. Invent. Math. 224 , 1 (2021), 55–146
2021
-
[26]
Kang, M.-J., and V asseur, A. F. Well-posedness of the Riemann problem with two shocks for th e isentropic Euler system in a class of vanishing physical vis cosity limits. J. Differential Equations 338 (2022), 128–226
2022
-
[28]
F., and W ang, Y
Kang, M.-J., V asseur, A. F., and W ang, Y. Time-asymptotic stability of generic Riemann solutions for compressible Navier-Stokes-Fourier equations. arXiv preprint arXiv:2306.05604 (2023)
2023 arXiv
-
[29]
Asymptotic stability of the stationary solution to the compressible Navier-Stokes equations in the half space
Kaw ashima, S., Shinya, N., and Zhu, P. Asymptotic stability of the stationary solution to the compressible Navier-Stokes equations in the half space. Comm. Math. Phys. 240 , 3 (2003), 483–500
2003
-
[30]
Asymptotic stability of nonlinear wave for the compressibl e Navier-Stokes equations in the half space
Kaw ashima, S., Shinya, N., and Zhu, P. Asymptotic stability of nonlinear wave for the compressibl e Navier-Stokes equations in the half space. J. Differential Equations 244 , 12 (2008), 3151–3179
2008
-
[31]
Asymptotic stability of rarefaction wave for the Navier-St okes equations for a compressible fluid in the half space
Kaw ashima, S., and Zhu, P. Asymptotic stability of rarefaction wave for the Navier-St okes equations for a compressible fluid in the half space. Arch. Ration. Mech. Anal. 194 , 1 (2009), 105–132
2009
-
[32]
Inflow and outflow problems in the half space for a one-dimensi onal isentropic model system of compressible viscous gas
Matsumura, A. Inflow and outflow problems in the half space for a one-dimensi onal isentropic model system of compressible viscous gas. Methods and Applications of Analysis 8 , 4 (2001), 645–666
2001
-
[33]
Convergence to travelling fronts of solutions of the p-syst em with viscosity in the presence of a boundary
Matsumura, A., and Mei, M. Convergence to travelling fronts of solutions of the p-syst em with viscosity in the presence of a boundary. Arch. Ration. Mech. Anal. 146 , 1 (1999), 1–22
1999
-
[34]
On asymptotic stability of noncharacteristic viscous boun dary layers
Nguyen, T. On asymptotic stability of noncharacteristic viscous boun dary layers. SIAM J. Math. Anal. 42, 3 (2010), 1156–1178
2010
-
[35]
Long-time stability of large-amplitude noncharacteristi c boundary layers for hyperbolic-parabolic systems
Nguyen, T., and Zumbrun, K. Long-time stability of large-amplitude noncharacteristi c boundary layers for hyperbolic-parabolic systems. J. Math. Pures Appl. (9) 92 , 6 (2009), 547–598
2009
-
[36]
Long-time stability of multi-dimensional noncharacteris tic viscous boundary layers
Nguyen, T., and Zumbrun, K. Long-time stability of multi-dimensional noncharacteris tic viscous boundary layers. Comm. Math. Phys. 299 , 1 (2010), 1–44
2010
-
[37]
F., and Yu, Q
W ang, Y., Yang, Y. F., and Yu, Q. Y. Existence of large boundary layer solutions to inflow proble m of 1D full compressible Navier-Stokes equations. arXiv preprint arXiv: 2504.21830 (2025). (Moon-Jin Kang) Department of Mathematical Sciences, Korea Advanced Institute of Mathem...
2025 arXiv
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.