Cosmin Burtea

Boris Haspot
email: haspot@ceremade.dauphine.fr

Vanishing capillarity limit of the Navier-Stokes-Korteweg system in one dimension with degenerate viscosity coefficient and discontinuous initial density

Vanishing capillarity limit of the Navier-Stokes-Korteweg system in one dimension with degenerate viscosity coefficient and discontinuous initial density

Cosmin Burtea, Boris Haspot

To cite this version:

Cosmin Burtea, Boris Haspot. Vanishing capillarity limit of the Navier-Stokes-Korteweg system in one dimension with degenerate viscosity coefficient

Introduction 1.Presentation of the models

One classical model used to study mixtures of two or more compressible fluids with different densities is the Navier-Stokes system (referred as NS in the sequel):

   ∂ t ρ + ∂ x (ρu) = 0, ∂ t (ρu) + ∂ x ρu 2 -∂ x (µ (ρ) ∂ x u) + ∂ x (aρ γ ) = 0, (ρ, u) |t=0 = (ρ 0 , u 0 ) , (1.1)
The mixture is supposed to be itself a compressible barotropic fluid and, in order to simplify matters, we suppose that it has constant temperature. Practical examples include gas bubbles in water flowing in pipes, fluids containing a melted substance, polluted air/water ..., see [START_REF] Ishii | Thermo-fluid dynamics of two-phase flow[END_REF] for more relevant examples. The unknowns are the velocity of the fluid u = u(t, x) ∈ R and the density of the fluid ρ = ρ(t, x) ∈ R + . We consider here that the pressure P (ρ) = aρ γ is defined throughout a γ state-law with γ > 1 and a > 0. In this paper, we consider that the viscosity coefficient µ = µ(ρ) ≥ 0 is density-dependent, with µ (0) = 0.

The fact that the dissipative operator depends on the density itself is motivated by physical considerations. Keeping in mind the fact that (1.1) governs flows of mixtures, it is reasonable to suppose that each of the components has its own viscosity. We also mention that system (1.1) can be obtained via the Boltzmann equation through the Chapman-Enskog expansion to the second order (see [START_REF] Chapman | The mathematical theory of nonuniform gases: an account of the kinetic theory of viscosity[END_REF]) which ensures that the viscosity coefficient is then a function of the temperature. If we consider the case of isentropic fluids, this dependence is expressed by a dependence on the density function (we refer in particular to [START_REF] Hoff | The failure of continuous dependence on initial data for the Navier-Stokes equations of compressible flow[END_REF]). Let us also mention that the case µ(ρ) = ρ is related to the so called viscous shallow water system. This system with friction has been derived by Gerbeau and Perthame in [START_REF] Gerbeau | Derivation of Saint-Venant system for laminar shallow water (numerical validation[END_REF] from the Navier-Stokes system with a free moving boundary in the shallow water regime at the first order. This derivation relies on the hydrostatic approximation where the authors follow the role of viscosity and friction on the bottom. It transpires that the mathematical analysis of this situation is more involved than the constant viscosity case since the operator -∂ x (µ (ρ) ∂ x u) is not strongly elliptic.

The system (1.1) is a Sharp Interface Model : assuming that at initial time the constituents of the mixture are separated by a jump in the initial density ρ 0 , then this discontinuity persists for all later times, and it is advected by the fluid's velocity see the early work of D. Hoff and J. Smoller [START_REF] Hoff | Solutions in the large for certain nonlinear parabolic systems[END_REF]. This allows one to distinguish at any time the zones occupied by the "pure" constituents.

This picture changes if one takes in consideration capillary forces. In this case the two fluids are not anymore separated by a sharp interface but rather by a thin layer where the density, although passing continuously from one fluid to another, it experiences large variations. This type of models are called Diffusive Interface Models. One classical model used to describe this situation is the Navier-Stokes-Korteweg (NSK in the sequel) system which, in its 1d version reads   

∂ t ρ + ∂ x (ρu) = 0, ∂ t (ρu) + ∂ x ρu 2 -∂ x (µ (ρ) ∂ x u) + ∂ x (aρ γ ) = c∂ x K, (ρ, u) |t=0 = (ρ 0 , u 0 ) .
(1.2)

In the above system K stands for the Korteweg capillarity tensor and its general formula is

K = ρκ (ρ) ∂ xx ρ + 1 2 (ρκ (ρ) -κ (ρ)) (∂ x ρ) 2 , (1.3)
for some positive function κ = κ (ρ). Moreover, c > 0 is a positive constant measuring the "strength of the capillarity" or in an other way "the thickness of the diffusive interface". Indeed the size of c is directly linked with the thickness of the transition zone separating the pure phases. For the underlying physical principles behind the governing equations, see the pioneering paper by J.-E. Dunn and J. Serrin in [START_REF] Serrin | On the thermodynamics of interstitial working[END_REF] but also [VdW79, Kor01, SD85, NT92, AMW98, CH58].

Let us observe that formally, when c → 0, we recover the NS system (1.1). Justifying mathematically such a limit in a physically interesting setting, namely, where interfaces could be tracked down, raises the following questions: Q1 Can we construct solutions for the system (1.1) with initial data ρ 0 allowing discontinuities? How "wildly" can these discontinuities behave?

Q2 Can we recover solutions of (1.1) with discontinuous initial data ρ 0 as the limit of solutions of (1.2)?

In this paper, we give positive answers to the above questions. Of course, we cannot solve the problem for general µ (ρ) and κ (ρ), additional hypothesis which ensure some extra algebraic structure to the underlying PDEs systems are needed. One typical example we can treat is µ (ρ) = ρ α with α ∈ (0, 1 2 ) and κ (ρ) = µ 2 (ρ) ρ 3 .

(1.4)

In order to obtain such results, we propose two new technical features in the context of the NSK system. Before, explaining these features, let us mention that, it is by now well-established, that the crucial point for obtaining global existence results is to be able to obtain apriori estimates which ensure that the density ρ stays bounded and bounded away from vacuum (i.e. ρ > 0) all along the time. Indeed, explosion criteria shows that this is the possible cause for the breakdown of strong solutions, see for instance Theorem 1.1. from [START_REF] Constantin | Compressible fluids and active potentials[END_REF] and Theorem 3.3. from [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF][START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF]. Even, for less regular solutions in the spirit of Hoff-Serre [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF][START_REF] Serre | Solutions faibles globales des équations de Navier-Stokes pour un fluide compressible[END_REF][START_REF] Serre | Sur l'équation monodimensionnelle d'un fluide visqueux, compressible et conducteur de chaleur[END_REF], the former qualitative information for the densities ensure that one can carry out estimation in the Hoff-class of regularity (referred to as intermediate-regularity solutions).

Regarding the second question Q2: in the case of the NSK system, first of all, we are able to use a maximum principle in order to obtain an L ∞ -bound on the density that does not degenerate in the vanishing capillarity coefficient limit c → 0. Secondly, we show that in order to obtain a lower bound for the density, it is sufficient to assume only a one sided inequality on the so-called effective velocity. At the level of the density, loosely speaking, this amounts to ask a one sided inequality:

∂ x ϕ(ρ 0 ) ≤ M 0 , (1.5)
for some M 0 ∈ R and where ϕ (ρ) is a primitive of µ(ρ) ρ 2 . Obviously, the expression in (1.5) has a meaning in the sense of measures. As these estimates are uniform w.r.t. capillarity parameter c, we are able to recover in the vanishing capillarity limit c → 0, solutions for the NS system where the initial data needs only to verify (1.5). We can observe that functions ϕ(ρ 0 ) ∈ BV loc with discontinuities with negative jump enter in the scope of (1.5). As far as the authors are aware, this is the first result in this direction, the only somehow related results known are obtained for the zero viscosity-capillarity limit by F. Charve and the second author in [START_REF] Charve | Existence of a global strong solution and vanishing capillarityviscosity limit in one dimension for the Korteweg system[END_REF] and by P. Germain and P.G. LeFloch in [START_REF] Germain | Finite energy method for compressible fluids: The Navier-Stokes-Korteweg model[END_REF]. The relative weak assumptions that are required on the initial density are in sharp contrast with what is usually assumed in the literature where the so-called Bresch-Desjardin entropy [START_REF] Bresch | On the existence of global weak solutions to the Navier-Stokes equations for viscous compressible and heat conducting fluids[END_REF] is used in order to ensure that the density stays away from vacuum (see [START_REF] Mellet | Existence and uniqueness of global strong solutions for one-dimensional compressible Navier-Stokes equations[END_REF]). However the assumption that the BD-entropy is finite at initial time allows only to consider continuous initial densities.

The method of proof reveals that in the case of the NS system (1.1), the situation is much more better: we can ensure that the density stays away from vacuum using only the fact that this is the case at the initial time t = 0. More precisely, under hypothesis (1.4) for µ (ρ) we are able to construct global strong solutions assuming only that ρ 0 and 1 ρ0 ∈ L ∞ . This gives an affirmative answer to Q1 at least for mildly degenerate density-dependent viscosities. Let us also point out that this result is also relevant with respect to recent results obtained by Bresch and Hillairet in [START_REF] Bresch | A compressible multifluid system with new physical relaxation terms[END_REF] in the context of multiphase modeling. Indeed the authors show that solutions for a Baer-Nunziato type-system can be obtained as week limits of highly-oscilating solutions of the NS system (1.1). The only norm that one can hope to keep uniformly bounded in the presence of high-oscilations is the L ∞ -norm. Our existence result stated in Theorem 1.1 thus enlarges the class of viscosity coefficients for which Theorem 1 from [START_REF] Bresch | A compressible multifluid system with new physical relaxation terms[END_REF] holds true. In particular, the aforementioned result holds true for viscosities that degenerate close to vacuum.

A short review of known results

Before, giving the formal statement of our main results we propose to the reader a short review of the results concerning the existence of global strong solutions for the NS and NSK systems. We mention here that we use the term weak solution to designate a distributional solution for which uniqueness is not know to hold true. The term strong solution will be used whenever we may ensure the uniqueness property even thought the derivatives appearing in the equations may not have a meaning a.e..

The Navier-Stokes system

We start mentioning some known results for the system (1.1).

The constant viscosity case

The study of the well-posedness of the Cauchy problem for the compressible Navier-Stokes equations with constant viscosity coefficients can be tracked back to the pioneering works of the Russian school of PDEs, see Ya. I. Kanel in [START_REF] Ya | On a model system of equations of one-dimensional gas motion[END_REF] and Kazhikhov and Shelukhin [START_REF] Vasil | Unique global solution with respect to time of initial-boundary value problems for one-dimensional equations of a viscous gas[END_REF] where they considered the case of regular initial data (in particular, the densities are always continuous). In the 80 s D. Hoff and J. Smoller started a program investigating the well-posedness issues for parabolic equations with rough initial data. Up to the knowledge of the authors, they were the first to observe that the discontinuities of the density are advected by the flow and they persist all along the time owing to the regularization of the so called effective viscous flux [START_REF] Hoff | Solutions in the large for certain nonlinear parabolic systems[END_REF]. The first results tackling the question of rough initial densities for the NS system were obtained by D. Hoff [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF] and D. Serre [START_REF] Serre | Solutions faibles globales des équations de Navier-Stokes pour un fluide compressible[END_REF][START_REF] Serre | Sur l'équation monodimensionnelle d'un fluide visqueux, compressible et conducteur de chaleur[END_REF]. Qualitative properties for weak-solutions like the non-formation of vacuum states were adressed in [START_REF] Hoff | Non-formation of vacuum states for compressible Navier-Stokes equations[END_REF]. We mention also the result of D. Hoff and D. Serre [START_REF] Hoff | The failure of continuous dependence on initial data for the Navier-Stokes equations of compressible flow[END_REF] where they proved that the assumption of constant viscosity leads to a failure of continuous dependence on the initial data. Finally D. Hoff extended the result of [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF] by showing in [START_REF] Hoff | Global solutions of the equations of one-dimensional, compressible flow with large data and forces, and with differing end states[END_REF] the existence of global weak solution with initial density admitting shocks (roughly speaking the density ρ 0 satisfies ρ 0 , 1 ρ0 ∈ L ∞ ) and with regularizing effects on the velocity, The density dependent case For initial densities that feature vacuum zones, existence of global weak solution has been obtained by Q.Jiu and Z. Xin in [START_REF] Jiu | The Cauchy problem for 1d compressible flows with density-dependent viscosity coefficients[END_REF] for viscosity coefficients verifying µ (ρ) = ρ α with α > 1 2 . In [START_REF] Haspot | Vortex solutions for the compressible Navier-Stokes equations with general viscosity coefficients in 1d: regularizing effects or not on the density[END_REF], the second author constructed global weak solutions for general viscosity coefficients with initial density admitting shock type discontinuities and with initial velocity belonging to the set of finite measures. In opposite to [START_REF] Hoff | Global solutions of the equations of one-dimensional, compressible flow with large data and forces, and with differing end states[END_REF], the initial data satisfy the BD entropy but not the classical energy which allows in particular to obtain regularizing effects for the density inasmuch as the density becomes instantaneously continuous.

We recall that the basic energy estimate, obtained formally by multiplying the velocity's equation with u and integrating by parts states that

R 1 2 ρu 2 (t, x) + ρe (ρ) (t, x) dx + t 0 R µ(ρ)(∂ x u) 2 (s, x)dsdx ≤ R 1 2 ρ 0 u 2 0 (t, x) + ρ 0 e (ρ 0 ) (x) dx, (1.6) with ρe (ρ) -e (ρ) = ρ γ .
In the context of the multidimensional NS system, D. Bresch and B. Desjardins [START_REF] Bresch | On the existence of global weak solutions to the Navier-Stokes equations for viscous compressible and heat conducting fluids[END_REF] discovered that (1.1) has some "hidden" algebraic structure. We can observe (see [START_REF] Haspot | Existence of global strong solution for the compressible Navier-Stokes equations with degenerate viscosity coefficients in 1d[END_REF]) that the so-called effective velocity:

v = u + ∂ x ϕ (ρ)
verifies the equation

∂ t (ρv) + ∂ x (ρuv) + ∂ x ρ γ = 0,
and that multiplying with v the previous relation the following functional, named BD entropy, is formally controlled:

R 1 2 ρv 2 (t, x) + ρe (ρ) (t, x) dx + γ t 0 R µ (ρ) ρ γ-3 (∂ x ρ) 2 ≤ R 1 2 ρ 0 v 2 0 (t, x) + ρ 0 e (ρ 0 ) (x) dx, (1.7)
We observe that if µ (ρ) = ρ α then the BD entropy shows that

R 1 2 ρ(v -u) 2 (t, x)dx = C α R (∂ x ρ α-1 2 ) 2 (t, x)dx Thus, supposing that ∂ x ρ α-1 2 0 ∈ L 2
, and taking 0 < α < 1 2 , the BD entropy yields that the L ∞ -norm of 1 ρ remains bounded. This was the main argument that A. Mellet and A. Vasseur [START_REF] Mellet | On the barotropic compressible Navier-Stokes equations[END_REF] used in order to construct global strong solutions for the case of degenerate viscosity coefficients with α satisfying 0 < α < 1 2 . In [START_REF] Haspot | Existence of global strong solution for the compressible Navier-Stokes equations with degenerate viscosity coefficients in 1d[END_REF], the second author has proved similar results for the case 1/2 < α ≤ 1 where he exploited the fact that the effective velocity v satisfies a damped transport equation. This allows him to obtain an L ∞ -estimates for v and using a maximum principle in order to obtain a L ∞ control on 1 ρ .

More recently Constantin et al in [START_REF] Constantin | Compressible fluids and active potentials[END_REF] have extended the previous result to the case α > 1 with γ belonging to [α, α + 1] provided that the initial data satisfy:

∂ x u 0 ≤ ρ γ-α 0 .
(1.8)

We point out that the condition (1.8) amounts to considering a negative effective flux (see for example [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF][START_REF] Lions | Compressible Models[END_REF]) at initial time. The main idea of their proof consists in proving via a maximum principle that the effective flux i.e. the function µ (ρ) ∂ x u -ρ γ , remains negative for all time t ≥ 0 if this is the case initially. This is sufficient to control the L ∞ norm of 1 ρ . In [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF] we proved that if α > 1 2 , γ ≥ max{1, α} and if the so-called effective velocity satisfies initially an Oleinik type inequality (see [START_REF] Arsen | Uniqueness and stability of the generalized solution of the Cauchy problem for a quasi-linear equation[END_REF]) then we have the existence of global strong solution provided that (ρ 0 , 1 ρ0 ) are bounded. In particular there is no restriction on the sign of the effective flux.

Let us mention that all the above results, assume that the BD entropy is bounded initially, it allow sat least when α > 1 2 to control the L ∞ norm of the density ρ all along the time by using Sobolev embedding on the quantity ∂ x ρ α-1 2 ∈ L 2 . Since this involves an L 2 information for the derivative of the density, all the above results deal with continuous initial densities.

Concerning discontinious initial densities we are only aware of two results: Fang and Zhang [START_REF] Fang | Global solutions of the Navier-Stokes equations for compressible flow with density-dependent viscosity and discontinuous initial data[END_REF] dealing with the free boundary problem and Ruxu et al [START_REF] Ruxu | Cauchy problem for the one-dimensional compressible navier-stokes equations[END_REF] dealing with the problem on the whole real line but assuming small initial energy. The estimates that the authors obtain depend on the number of points of discontinuity and, in particular, these estimates blow-up as the number of discontinuity points tends to +∞. For this reason, they are not appropiate in order to study 1d-homogeneisaition as in [START_REF] Bresch | A compressible multifluid system with new physical relaxation terms[END_REF].

As it was mentioned previously, in this paper we will provide a new method to control the L ∞ norm of the density, without using the BD entropy. By doing so, our initial density ρ 0 is just

L 2 (R) ∩ L ∞ (R)

The Navier-Stokes-Korteweg system

Let us pass now in review some results for system (1.2). F. Charve and the second author in [START_REF] Charve | Existence of a global strong solution and vanishing capillarityviscosity limit in one dimension for the Korteweg system[END_REF] proved the global existence of strong solution for the system (1.2) when µ(ρ) = ερ and κ(ρ) = ε 2 ρ . In addition, they show that the global strong solutions converge when ε goes to 0 to a global weak entropy solution of the compressible Euler system with initial data of finite energy. P. Germain and P.G. LeFloch in [START_REF] Germain | Finite energy method for compressible fluids: The Navier-Stokes-Korteweg model[END_REF] studied recently the global existence of vacuum and non-vacuum weak solutions for the Korteweg system. It is important to point out that they need to impose a tame condition on the viscosity and capillary coefficients which takes the form:

κ(ρ) µ(ρ) 2 ρ 3 and δ(ε) ε 2 , (1.9)
if we consider the vanishing viscosity capillary limit for viscosity and capillary coefficients µ ε (ρ) = εµ(ρ), κ ε (ρ) = δ(e)κ(ρ) when ε > 0 goes to 0. Roughly speaking, the previous tame condition implies in some sense that the parabolic behavior, governed by the viscosity tensor dominates the dispersive effects that are induced by the capillarity tensor. We mention that this is important because the dispersive effects tend to create strong oscillations which prevent obtaining strong convergence informations (they allow only to obtain weak-convergence results which are not well-suited to treat the compressible setting). In particular, when ε goes to 0, by assuming (1.9) we can expect recovering strong convergence in suitable functional spaces. In [START_REF] Germain | Finite energy method for compressible fluids: The Navier-Stokes-Korteweg model[END_REF] the authors study then the zero viscosity-capillarity limit associated with the Navier-Stokes-Korteweg system generalizing the results of [START_REF] Charve | Existence of a global strong solution and vanishing capillarityviscosity limit in one dimension for the Korteweg system[END_REF]. They need in particular for this analysis to assume that κ ε (ρ) = µε(ρ) 2 ρ 3

when the viscosity coefficient is degenerate with µ ε (ρ) = ερ α . We recall that we will also consider this algebraic relation in the remaining part of this paper. Recently, Chen et al. in [START_REF] Chen | Global classical solutions to the one-dimensional compressible fluid models of Korteweg type with large initial data[END_REF] and Chen in [START_REF] Chen | Asymptotic stability of strong rarefaction waves for the compressible fluid models of Korteweg type[END_REF] have proved for the first time some results of existence of global strong solutions for initial density far away from the vacuum in Lagrangian coordinates. More precisely they consider viscosity and capillary coefficients of the form µ(ρ) = ρ α1 and κ(ρ) = ρ β1 with (α 1 , β 1 ) ∈ R 2 . In comparison with the present work, there is no relation a priori between α 1 and β 1 (furthermore there is no restriction on c > 0). They manage essentially to show such result when β 1 < -2 which allows in a direct way to control the L ∞ norm of 1 ρ by using the energy estimate (indeed roughly speaking the energy estimate ensure that ∂ x ρ

β 1 2 +1 is bounded in L ∞
T (L 2 ) for any T > 0). They deal also with the case β 1 ≥ -2 but in this case α 1 < 0, in particular the viscosity coefficient explodes near vacuum. The main ideas of the proof is to obtain L 2 -estimates for the effective velocity v = u + µ(ρ) ρ 2 ∂ x ρ by using an energy method combined with Sobolev embedding in the spirit of Kanel (see [START_REF] Ya | On a model system of equations of one-dimensional gas motion[END_REF]). Furthermore the authors show also the existence of global strong solution when the initial data is a perturbation of a Riemann problem associated to a rarefaction wave for the compressible Euler problem.

In [START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF] the authors extend the results of [START_REF] Chen | Asymptotic stability of strong rarefaction waves for the compressible fluid models of Korteweg type[END_REF][START_REF] Chen | Global classical solutions to the one-dimensional compressible fluid models of Korteweg type with large initial data[END_REF] to the case of the NSK system with strongly degenerate viscosity coefficients with initial density far away from vacuum. More precisely, the viscosity coefficients take the form µ(ρ) = ρ α with α > 1 and the capillarity coefficient satisfies the algebraic relation (1.10). The main difficulty of the proof consists in estimating globally in time the L ∞ norm of 1 ρ . The method of proof relies on fine algebraic properties of the NSK system. In [START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF] we introduce two new effective velocities endowed with weight functions depending both on the viscosity and the capillarity coefficients as some power laws of the density. For these two quantities we show some Oleinik-type estimate which provide the control of the L ∞ norm of 1 ρ by applying a maximum principle. It is interesting to point out that the two effective pressure introduced in this paper depending on the capillary coefficient generalize to the NSK systems those introduced for the NS system in [BH20, CDNP20].

Main results

The algebraic structure of the NSK system

For the rest of the paper, as in [GL16, BGLV19], we will assume that:

κ (ρ) = µ 2 (ρ) ρ 3 . (1.10)
With such a choice for κ we can rewrite the capillarity tensor K as

K = µ (ρ) ∂ xx ϕ (ρ) with ϕ (ρ) = µ (ρ) ρ 2 . (1.11)
We observe also that:

∂ x K = ρ∂ x G (ρ) ∂ 2 xx G (ρ) with G (ρ) = µ (ρ) ρ 3/2 .
(1.12)

We will study the system (1.2) on the real line R with the following far field assumption:

ρ (t, x) → 1 and u (t, x) → 0 when |x| → ∞, (1.13)
for all t ≥ 0. Let r ∈ (0, 1) and observe that using the second equation of the system (1.2) we can write that:

ρ∂ t (u + r∂ x ϕ (ρ)) + ρu∂ x (u + r∂ x ϕ (ρ)) -(1 -r) ∂ x (µ(ρ)∂ x (u + r∂ x ϕ (ρ))) + a∂ x ρ γ = r 2 -r + c ∂ x K.
For c ∈ (0, 1 4 ] the equation r 2 -r + c = 0 admits the two positive roots

r 1 (c) = 1 + √ 1 -4c 2 ≤ 1 and r 0 (c) = 1 - √ 1 -4c 2 ≤ 1.
In the sequel we will assume that c ∈ (0, 1 4 ] and we observe that

r 1 (c) -→ c→0 1 and r 0 (c) -→ c→0 0. (1.14)
We introduce the following two extra variables

v i = u + r i ∂ x ϕ (ρ) for i ∈ {0, 1} ,
which will be referred as effective velocities. The effective velocities verify the following equations:

ρ∂ t v i + ρu∂ x v i -(1 -r i (c)) ∂ x (µ(ρ)∂ x v i ) + a∂ x ρ γ = 0 for i ∈ {0, 1} .
(1.15)

Let us discuss now the dissipation of energy. Using (1.12) the natural energy associated to

(1.2) is R 1 2 ρu 2 (t, x) + ρe (ρ) (t, x) + 1 2 c (∂ x G (ρ)) 2 (t, x) dx + t 0 R µ(ρ)(∂ x u) 2 (s, x)dsdx ≤ E c (ρ 0 , u 0 ) , (1.16)
where

E c (ρ 0 , u 0 ) = 1 2 R ρ 0 u 2 0 (x) + ρ 0 e (ρ 0 ) (x) + cE cap (∂ x ρ 0 ) dx, (1.17) with E cap (∂ x ρ 0 ) = 1 2 R (∂ x G (ρ 0 )) 2 (x)dx = 1 2 R ρ 0 (∂ x ϕ (ρ 0 )) 2 (x)dx and e (ρ) = (ρ γ -1 -γ (ρ -1)) (γ -1) ρ = ρ γ-1 γ -1 + 1 ρ - γ γ -1 . (1.18)
In the context of the Navier-Stokes-Korteweg system (1.2), since there are two effective-velocities we have two BD-type entropies obtained by multiplying the equation (1.15) with v i and integrating by parts:

R 1 2 ρv 2 i (t, x) + ρe (ρ) (t, x) dx + (1 -r i (c)) t 0 R µ(ρ)(∂ x v i ) 2 (s, x)dsdx + r i (c) γ t 0 t 0 µ (ρ) ρ γ-3 (∂ x ρ) 2 (s, x)dsdx (1.19) ≤ R 1 2 ρv 2 i,0 (x) + ρ 0 e (ρ 0 ) (x) dx = R 1 2 ρ 0 u 2 0 (x) + ρ 0 e (ρ 0 ) (x) dx + r i R ρ 0 u 0 ∂ x ϕ (ρ 0 ) (x)dx + r 2 i E cap (∂ x ρ 0 ) . (1.20)
If one ignores the relative size with respect to c of the coefficients r i (c) , both estimates (1.20) give the same qualitative information. Moreover, the spaces where the initial data has to be drawn from, in order to give a mathematical meaning to these functionals are such that E c (ρ 0 , u 0 ) defined in (1.17) is finite. In particular, assuming finite energy, the density in the NSK system is always a continuous function. Of course, this is normal since, as discussed in the Introduction, NSK is a Diffuse Interface Model. However, observe that as c → 0, the estimate for v 0 degenerates in the basic energy estimate for Navier-Stokes (1.6), while the estimate for v 1 degenerates to the BD-entropy estimate (1.7). This observation reveals that if one is interested in recovering a "rough" solution of the NS system by a vanishing capillarity limit, the BD entropy estimate for v 1 must explode as c → 0.

Assumptions for the viscosity coefficient

We wish now to specify the form of the viscosity coefficient with which we will work. We recall that

ϕ (ρ) = µ (ρ) ρ 2 and ψ (ρ) = µ (ρ) ρ • (1.21) Also, we introduce Ξ (ρ) = (ρ γ µ (ρ)) 1 2 ρ 2 • (1.22)
We suppose that µ is a C 1 positive function on (0, +∞) such that:

lim ρ→+∞ µ (ρ) = ∞ (µ 1 )
We assume that there exists positive constants d 1 , d 2 , d 3 , d 4 , d 5 such that the following relations hold true. First of all, for all ρ ≥ 0 we suppose that 1

d 1 µ(ρ) ≤ ψ (ρ) ≤ d 1 µ (ρ) . (µ 2 )
Next, we suppose that

µ (ρ) ≤ d 2 (1 + ρ γ ) for all ρ ≥ 0. (µ 3 )
We also require that ϕ : (0, ∞) → Imϕ is invertible and that

lim ρ→0 ϕ (ρ) = -∞. (µ 4 ) We impose that lim ρ→+∞ Ξ (ρ) = +∞. (µ 5 )
Finally, we denote

Φ (τ ) = -ϕ 1 τ > 0 , Λ (τ ) = 1 µ 1 τ > 0 (1.23)
and we ask that

Λ (τ ) ≤ d 4 (1 + Φ (τ )) 2 (µ 6 ) along with τ 1 2 ≤ d 5 (1 + Φ (τ )) 1-η (µ 7 )
for some positive η ∈ (0, 1).

Remark 1 Let us investigate how the above hypothesis translate in the case of viscosity coefficient of the form µ (ρ) = ρ α . Observe that in this case we have that

ϕ(ρ) = ρ α-1 α -1 and ψ(ρ) = ρ α α if α = 1 while ϕ(ρ) = ln ρ and ψ(ρ) = ρ if α = 1.
We see that, (µ 1 ), (µ 2 ), (µ 3 ) imply that 0 < α ≤ γ.

Hypothesis (µ 4 ) implies that α ≤ 1.

Hypothesis (µ 5 ) implies that α + γ ≥ 2.
Finally, since

Λ (τ ) = τ α and Φ (τ ) = 1 1-α τ 1-α if α = 1, the hypothesis (µ 6 ) implies that α ≤ 2 3 while (µ 7 ) implies that α < 1 2 .

Statement of the main results

Finally, after a rather long introduction, we are in the position of stating our first main result which addresses the question of global existence of strong solutions with rough initial densities.

Theorem 1.1 Consider γ > 1 and µ verifying (µ 1 ), (µ 2 ), (µ 3 ), (µ 4 ),(µ 6 ) and (µ 7 ).

Consider (ρ 0 , u 0 ) ∈ L ∞ (R) × L 2 (R) such that ρ 0 L ∞ (R) + 1 ρ 0 L ∞ (R) + E 0 (ρ 0 , u 0 ) ≤ M.
Then, there exists a unique global solution (ρ, u) for the Navier-Stokes system (1.1). Moreover, there exists

C = C T, E 0 (ρ 0 , u 0 ) , ρ 0 , 1 ρ0 L ∞
such that the pair (ρ, u) satisfies for all T > 0:

E 0 (ρ (T ) , u (T )) = R ρu 2 + ρe (ρ) (T ) + T 0 R µ(ρ)(∂ x u) 2 ≤ E 0 (ρ 0 , u 0 ) , (1.24) 1 C ≤ ρ(T, •) ≤ C, (1.25) T 0 R σρ | u| 2 + 1 2 σ (T ) R µ (ρ (T )) (∂ x u (T )) 2 ≤ C, (1.26) 1 2 T 0 R σ 2 ρ | u| 2 + (1 -r 0 ) T 0 R σ 2 µ (ρ) |∂ x u| 2 ≤ C, (1.27) T 0 σ 1 2 (τ ) ∂ x u (τ ) 2 L ∞ dτ + sup 0<t≤T σ(t) ∂ x u(t, •) L ∞ ≤ C, (1.28)
where σ (t) = min {1, t} and u = ∂ t u + u∂ x u.

Remark 2 Up to your knowledge this is the first result assuring the existence of global unique solutions for the compressible Navier-Stokes system in one dimension with degenerate viscosity coefficient and with rough initial density ρ 0 ∈ L ∞ . Of course, this includes the case of discontinuous initial data, with an arbitrary number of discontinuities and arbitrary jumps. In particular, this extends the works [BH20, CDNP20, MV08] where the authors consider initial densities which are continuous.

Remark 3 According to the Remark 1, we recall that the choice µ(ρ) = ρ α with α ∈ (0, 1/2) enters in the framework of Theorem 1.1. Since we do not need to assume (µ 5 ), any p (ρ) = ρ γ with γ > 1 fits in the framework of our result.

Remark 4 At this point it is important to remark that, using Proposition 2.4.1 proved below, the result obtained by Constantin et al. [START_REF] Constantin | Compressible fluids and active potentials[END_REF] can be adapted to obtain an existence result with discontinuous initial densities ρ 0 such that

∂ x u 0 -ρ γ-α 0 ≤ 0. (1.29)
Let us remark however that, at least if ∂ x u 0 is continuous, condition that although cannot be ensured for later times, it nevertheless implies that the jump of ρ γ-α 0 at a discontinuity must be positive. This observation shows that general initial data relevant for multifluids at a mesoscopic scale, as explained in [BDG + 18], Section 3.2. cannot be considered assuming (1.29) initially.

The second main result states that, if we impose a sign condition on the effective velocity v 1 then, we can construct solutions for the NSK system (1.2) which verify uniform estimates for the density and the effective velocity v 0 that do not degenerate as c → 0. Recall that formally, when c → 0, v 0 degenerates to u, the velocity of the NS system.

Theorem 1.2 Let c ∈ (0, 1 4 ), γ > 1 and µ, ϕ, ψ verifying (1.21), (µ 1 )-(µ 7 ). Consider the initial data (ρ 0 , u 0 ) which is uniformly bounded in c in the following space

u 0 , √ cρ 1 2 0 ∂ x ϕ (ρ 0 ) ∈ L 2 (R) , ρ 0 e (ρ 0 ) ∈ L 1 (R)
. Moreover, suppose that there exists M 0 ∈ R such that for all c ∈ 0, 1 4 we have a.e. x ∈ R :

v 1|t=0 = u 0 (x) + r 1 (c) ∂ x ϕ(ρ 0 (x)) ≤ M 0 .
Then, there exists a unique global solution for (1.2) that satisfies the following estimates, uniformly with respect to the parameter c:

R 1 2 ρu 2 + ρe (ρ) + cρ (∂ x ϕ (ρ)) 2 + t 0 R µ(ρ)(∂ x u) 2 ≤ E c (ρ 0 , u 0 ) , R 1 2 ρv 2 0 + ρe (ρ) t 0 R + (1 -r 0 (c)) t 0 R µ(ρ)(∂ x v i ) 2 (s, x)dsdx + r 0 (c) γ t 0 R µ (ρ) ρ γ-3 (∂ x ρ) 2 (s, x)dsdx ≤ 2E c (ρ 0 , u 0 ) , ρ (t, x) ≤ C (t, E c , ρ 0 L ∞ ) , u (t, x) + r 1 (c) ∂ x ϕ (ρ (t, x)) ≤ M 0 + C (t, E c , ρ 0 L ∞ ) , 1 ρ (t, •) L ∞ ≤ C t, E c , M 0 , ρ 0 , 1 ρ0 L ∞ ϕ (ρ) BV ([0,T ]×[-L,L]) + ρ BV ([0,T ]×[-L,L]) ≤ C t, L, E c , M 0 , ρ 0 , 1 ρ0 L ∞ ∀L > 0. T 0 R σρ | v0 | 2 + 1 -r 0 (c) 2 σ (T ) R µ (ρ (T )) (∂ x v 0 (T )) 2 ≤ C t, E c , M 0 , ρ 0 , 1 ρ0 L ∞ , 1 2 T 0 R σ 2 ρ | v0 | 2 + (1 -r 0 (c)) T 0 R σ 2 µ (ρ) |∂ x v0 | 2 ≤ C t, E c , M 0 , ρ 0 , 1 ρ0 L ∞ , (1.30)
where σ (t) = min {1, t} and v0 = ∂ t v 0 + u∂ x v 0 . Above, C are generic functions depending continuously on their arguments and they are increasing in t.

As we pointed out in the introduction, comparing with [START_REF] Germain | Finite energy method for compressible fluids: The Navier-Stokes-Korteweg model[END_REF] and [START_REF] Chen | Asymptotic stability of strong rarefaction waves for the compressible fluid models of Korteweg type[END_REF][START_REF] Chen | Global classical solutions to the one-dimensional compressible fluid models of Korteweg type with large initial data[END_REF], our contribution consists in the uniform estimates with respect to c for ρ L ∞ and the fact that we only use the one sided bound for v 1|t=0 = u 0 + r 1 ∂ x ϕ (ρ 0 ) in order to show that ρ is bounded from bellow. Of course this is crucial if one is interested in the capillarity-vanishing limit problem (the viscosity remaining constant).

As a corollary to Theorem 1.2 we immediately obtain

Theorem 1.3 Let c ∈ (0, 1 4 ), γ > 1 and µ, ϕ, ψ verifying (1.21), (µ 1 )-(µ 6 ). Let M 0 ∈ R, M > 0 and consider (ρ c 0 , u c 0 ) such that ρ c 0 L ∞ (R) + 1 ρ c 0 L ∞ (R) + E c (ρ c 0 , u c 0 ) ≤ M, and
u c 0 (x) + r 1 (c) ∂ x ϕ(ρ c 0 ) (x) ≤ M 0 a.e. on R. Moreover, suppose that (ρ 0 , u 0 ) is such that ρ 0 e(ρ 0 ) ∈ L 1 (R), u 0 ∈ L 2 (R), (ρ 0 , 1/ρ 0 ) ∈ L ∞ (R), ϕ(ρ 0 ) ∈ BV loc (R) such that u 0 + ∂ x ϕ(ρ 0 ) ≤ M 0 in the sense of measures (1.31)
and that:

ρ c 0 -1 → ρ 0 -1, u c 0 → u 0 in L 2 (R) , u c 0 + r 1 ∂ x ϕ(ρ c 0 ) u 0 + ∂ x ϕ(ρ 0 )
weakly in the sense of measures.

(1.32)

Then, there exists (ρ, u) such that for all T > 0,

lim c→0 (ρ c , u c ) c>0 = (ρ, u) weakly -in L ∞ ([0, T ] × R)
and (ρ, u) is the unique global solution for the NS system (1.1) with initial data (ρ 0 , u 0 ) that verifies for all T > 0:

E 0 (ρ (T ) , u (T )) = R ρu 2 + ρe (ρ) (T ) + T 0 R µ(ρ)(∂ x u) 2 ≤ E 0 (ρ 0 , u 0 ) , (1.33) C t, E 0 (ρ 0 , u 0 ) , ρ 0 , 1 ρ 0 L ∞ -1 ≤ ρ(t, x) ≤ C t, E 0 (ρ 0 , u 0 ) , ρ 0 , 1 ρ 0 L ∞ for any x ∈ R. (1.34) ϕ (ρ) BV ([0,T ]×[-L,L]) + ρ BV ([0,T ]×[-L,L]) ≤ C t, L, E c , M 0 , ρ 0 , 1 ρ 0 L ∞ ∀L > 0. (1.35) u + ∂ x ϕ (ρ) ≤ M 0 + C (t, E 0 , ρ 0 L ∞ ) in the sense of measures (1.36) T 0 R σρ | u| 2 + 1 2 σ (T ) R µ (ρ (T )) (∂ x u (T )) 2 ≤ C = C T, E 0 (ρ 0 , u 0 ) , ρ 0 , 1 ρ 0 L ∞ , (1.37) 1 2 T 0 R σ 2 ρ | u| 2 T 0 R σ 2 µ (ρ) |∂ x u| 2 ≤ C = C T, E 0 (ρ 0 , u 0 ) , ρ 0 , 1 ρ 0 L ∞ , (1.38) T 0 σ 1 2 (τ ) ∂ x u (τ ) 2 L ∞ dτ + sup 0<t≤T σ(t) ∂ x u(t, •) L ∞ ≤ C = C T, E 0 (ρ 0 , u 0 ) , ρ 0 , 1 ρ 0 L ∞ , (1.39)
.

As far as we are aware, this is the first result treating the vanishing capillarity limit. As a corollary of the previous Theorem we obtain via a vanishing capillary process the existence of global strong solutions for the 1D Navier-Stokes system with discontinuous initial densities . However, besides the fact that ρ 0 and 1 ρ0 must be in L ∞ (R) we also require that

u 0 + ∂ x ϕ(ρ 0 ) ≤ M 0 (1.40)
for some M 0 ∈ R, this is obviously an additional condition compared to Theorem 1.1. In particular (1.40) implies that ϕ(ρ 0 ) belongs to BV loc and that the only discontinuities in x on ρ 0 are such as the jump of ϕ (ρ 0 ) at x is negative:

[ϕ (ρ 0 )] (x) = lim h→0,h>0 ϕ (ρ 0 (x + h)) -lim h→0,h>0 ϕ (ρ 0 (x -h)) ≤ 0.
Remark 5 We would like also emphasize that it is possible to obtain a similar theorem if we replace the condition (1.31) by:

u 0 + ∂ x ϕ (ρ 0 ) ≥ M 0 in the sense of measures, (1.41)
with M 0 ∈ R. We refer to Section 2.6 and Remark 6 for more details on these questions.

Proof of the main results

Let us give the general plan for the proof of the main results. The main difficulty is to obtain a priori estimates assuring that the density is bounded and bounded by below. Once this is achieved, one may follow the approach of D. Hoff in order to obtain the estimates necessary to prove existence and uniqueness, see [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF], [START_REF] Hoff | Global solutions of the equations of one-dimensional, compressible flow with large data and forces, and with differing end states[END_REF] for the original approach or our more recent contributions [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF], [START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF]. Thus, we will not insist on these, by now well-understood points. First, we recall the finite-time existence results of strong solutions and some explosion criterion stating that the only way in which a classical solution might blow-up is because of the appearance of vacuum regions or because the density does not remain bounded, more precisely, the L ∞ -norm of 1 ρ or of ρ blows up. These results were stated and proved in [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF] and [START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF] (see Theorem 3.1 from these papers).

We consider initial data as in Theorem 1.1 and Theorem 1.2 and we regularize it with a family of molifiers such as to fit in the scope of the local existence, see Theorem 2.1 below. For each mollified initial data we consider the regular solution defined on a maximal time of existence. The regularity of the solutions ensures that all the computations presented below are justified. In particular, the regularity of the local-in-time solution justifies the passage from the Eulerian formulation (1.2) to the so-called mass-Lagrangian formulation, see Section 2.2. The fact that this later formulation is more adapted to obtain apriori estimates is known since the work of Kazhikhov and Shelukhin [START_REF] Vasil | Unique global solution with respect to time of initial-boundary value problems for one-dimensional equations of a viscous gas[END_REF] and was used more recently by Germain and Le Floch in [START_REF] Germain | Finite energy method for compressible fluids: The Navier-Stokes-Korteweg model[END_REF].

A first important step is to show that the density is uniformly bounded using only the uniform bound on the energy and the L ∞ norm of ρ 0 . This is the objective of Section 2.4 . More precisely, we show that the local solution (or any regular enough solution) for the NSK system verifies

ρ (t, x) ≤ C (t, E c , ρ 0 L ∞ ) for all x ∈ R
where C is continuous and increasing with respect to t. This estimate is crucial if one wants to consider the vanishing capillarity limit. The proof is rather technical and uses a lot of tricks inspired by earlier works of the Russian school of PDEs, see [START_REF] Ya | On a model system of equations of one-dimensional gas motion[END_REF], [START_REF] Vasil | Unique global solution with respect to time of initial-boundary value problems for one-dimensional equations of a viscous gas[END_REF]. We refer the reader to the monographic [START_REF] Stanislav Nikolaevich Antontsev | Boundary value problems in mechanics of nonhomogeneous fluids[END_REF] for a systematic treatment of the 1d NS system with various boundary conditions. Loosely speaking, the idea is to work with a primitive of the function defining the velocity's equation. The fact that we work with functions defined on the whole real line adds another difficulty to the proof because we need to localize the arguments. We point out that Proposition 2.3.2 below turns out to be crucial in order to carry out our proof.

In Sections 2.5 and 2.6 we show that the density part of such a regular solution remains bounded and bounded below. In particular, this will ensure that its maximal time of existence is +∞. In the same way as working with a primitive for the momentum equation leads to an upper bound for the density under not to restrictive hypothesis, working with a primitive of the equation of the effective velocity will enable us to obtain a lower bound for the density with minimal regularity assumptions. More precisely, in the NSK case we show that the density is lower bounded by a constant that depends continuously on time, initial energy, ρ 0 L ∞ , 1 ρ0 L ∞ and on M 0 ∈ R defined as the smallest constant such that

v 1|t=0 = u 0 + r 1 ∂ x ϕ (ρ 0 ) ≤ M 0 , a.e. in R.
(2.1) Relation (2.1) translates the fact the the positive part of the effective velocity should be bounded but does not offer any information on the negative part. In particular, this explains why densities having negative jumps enter the framework of (2.1). In the case of the NS system, the information we obtain is better. Indeed, we can say that the effective velocity is controlled pointwise by the initial effective velocity plus a term that we control. This allows to show that ρ is lower bounded without assuming (2.1). Finally, in Section 3, armed with the uniform estimates for the density, we show how to carry-out the Hoff-type estimates program. With respect to our recent contribution [START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF], the Hoff-estimates for the effective velocity that degenerates to u, are shown to hold true uniformly with respect to the capillarity coefficient c which is important for the proof of Theorem 1.3.

Existence of strong solution in finite time

In order to prove the existence of global strong solutions for the NS and NSK systems, we start with recalling the following result.

Theorem 2.1 Assume that the viscosity coefficients satisfy the assumptions (µ 1 )-(µ 4 ), s ≥ 3 and

(ρ 0 -1, u 0 ) ∈ H s+1 × H s (R) with 1 ρ0 ∈ L ∞ (R).
Then there exists T * > 0 such that there exists a strong solution (ρ, u) of the system (1.2) on (0, T * ) with ∀T ∈ (0, T * ):

(ρ -1) ∈ C(0, T, H s+1 (R)) ∩ L 2 (0, T, H s+2 (R)), u ∈ C(0, T, H s (R)) ∩ L 2 (0, T, H s+1 (R)),
and for all t ∈ (0, T * ):

1 ρ (t, •) L ∞ ≤ C(t),
where C(t) < +∞ if t ∈ (0, T * ). In addition, if:

sup t∈(0,T * ) [ 1 ρ (t, •) L ∞ + ρ(t, •) L ∞ ] ≤ C < +∞,
then the solution can be continued beyond (0, T * ).

The above result claims that the only way a strong solution might blow-up in finite time T is if the L ∞ -norm of 1 ρ or ρ blows-up at time T * . We refer to [START_REF] Chen | Asymptotic stability of strong rarefaction waves for the compressible fluid models of Korteweg type[END_REF][START_REF] Chen | Global classical solutions to the one-dimensional compressible fluid models of Korteweg type with large initial data[END_REF] for the proof of existence of a strong solution in finite time and we refer to [START_REF] Burtea | Existence of global strong solution for Korteweg system in one dimension for strongly degenerate viscosity coefficients[END_REF] for a proof of the blow-up criterion. Since the initial data in Theorem 1.1 and Theorem 1.2 are less regular than in the Theorem 2.1, we cannot directly use Theorem 2.1.

Consider (ρ 0 , u 0 ) ∈ L ∞ (R) × L 2 (R) such that E c (ρ 0 , u 0 ) < ∞ and, in the case of the NSK system, u 0 + r 1 (c) ∂ x ϕ (ρ 0 ) ≤ M 0 a.e. on R.
We regularize the initial data as follows:

   ρ n 0 = j n * ρ 0 , v n 1|t=0 = j n * v 1|t=0 = j n * (u 0 + r 1 (c) ∂ x ϕ (ρ 0 )) , u n 0 = v n 1|t=0 -r 1 (c) ∂ x ϕ (ρ n 0 ) .
(2.2) with j n a family of mollifiers, j n (y) = nj(ny)

with j ∈ C ∞ (R) such that 0 ≤ j ≤ 1, R j(y)dy = 1 and Supp j ⊂ [-2, 2].
We deduce that (ρ n 0 -1, v n 1|t=0 ) belong to all Sobolev spaces H s (R) with s ≥ 0. Furthermore, by the composition theorem we can prove that ϕ(ρ n 0 )-ϕ(1) belongs to H k (R) for any k ≥ 0 and consequently we obtain that u n 0 ∈ H k (R) for k ≥ 3. However, the higher order Sobolev norms explode when n → ∞. The only informations that are uniform in n ∈ N, the regularization parameter, are the following:

     0 < 1 ρ0 L ∞ ≤ ρ n 0 ≤ ρ 0 L ∞ < +∞, E c (ρ n 0 , u n 0 ) ≤ 2E c (ρ 0 , u 0 ) , v n 1|t=0 = u n 0 + r 1 (c) ∂ x ϕ (ρ n 0 ) ≤ M 0 .
(2.3)

We can now apply the Theorem 2.1 to the sequence of initial data (ρ n 0 , u n 0 ) n≥0 which provides us the existence of a strong solution (ρ n , u n ) of the system (1.2) on some finite time interval (0, T n ) with T n > 0. We denote by

v n 0 = u n + r 0 (c) ∂ x ϕ(ρ n ) and by v n 1 = u n + r 1 (c) ∂ x ϕ(ρ n
) the effective velocities introduced in Section 1.3.1. Using (2.3) and the regularity of the constructed solution we see that for any t ∈ (0, T n ) and sufficiently large n ∈ N we have that

R 1 2 ρ n |u n | 2 (t, x) + ρ n e (ρ n ) (t, x) + c(∂ x G (ρ n )) 2 dx + t 0 R µ(ρ n )(∂ x u n ) 2 (s, x)dsdx ≤ 2E c (ρ 0 , u 0 ) < +∞ (2.4) R 1 2 ρ n |v n 0 | 2 (t, x) + ρ n e (ρ n ) (t, x) dx + (1 -r 0 (c)) t 0 R µ(ρ n )(∂ x v n 0 ) 2 (s, x)dsdx ≤ 2E c (ρ 0 , u 0 ) < +∞. (2.5)
We anticipate that at the end of Section 2.6 when we will have shown that 1 ρn is bounded, we will be able to conclude that for any n ∈ N we have T n = +∞, see Proposition 2.6.2.

The end of the proof of Theorem 1.2 and 1.1 will consist in proving that (ρ n , u n ) n∈N converges up to a subsequence to a global strong solution (ρ, u) of the system (1.2). In the sequel for simplicity we omit the subscript n ∈ N and all the estimate are obtained on the time interval (0, T n ).

The Mass-Lagrangian formulation

It is well known that if inf x∈R ρ 0 (x) > 0 for all x ∈ R, then the change of variables:

ρ(t, m) = ρ(t, X(t, Y (m))) and u(t, m) = u(t, X(t, Y (m)))
with:

X(t, x) = x + t 0 u(s, X(s, x))ds and Y -1 (m) = m 0 ρ 0 (z)dz.
where Y -1 is the inverse of Y (we observe in particular that Y -1 is strictly increasing on R since inf x∈R ρ 0 (x) > 0) transforms the system (1.2) and (1.15) into:

   ∂ t ρ + ρ 2 ∂ m u = 0, ∂ t u -∂ m ( ρµ ( ρ) ∂ m u) + ∂ m ρ γ = c∂ m ( ρµ ( ρ) ∂ 2 mm ψ ( ρ)), ∂ t v i -(1 -r i ) ∂ m ( ρµ ( ρ) ∂ m v i ) + ∂ m ρ γ = 0, (2.6)
where we denote

ψ (ρ) = ρϕ (ρ) = µ (ρ) ρ , (2.7)
and with m ∈ R. We note that we have taken a = 1 in the pressure term only for simplifying the notations. The conditions (1.13) yield ρ (t, m) → 1 and u (t, m) → 0 when |m| → ∞.

Obviously, along as the solutions are regular, for instance the regularity assumed in Theorem 2.1 is sufficient, then for all t > 0, the application m → X (t, Y (m))

is a C 1 -diffeormorfism from R to R with inverse denoted by

x → m (t, x) .

(2.8)

Obviously, we have that

R F (ρ (t, x) , u (t, x)) = R F ( ρ (t, m) , u (t, m)) dm ρ (t, m) , (2.9)
for any function F . Also, we have that

∂ x G (t, X (t, Y (m))) = ρ∂ m G (t, m) (2.10)
We will use system (2.6) in order to obtain a priori estimates as it is easier to manipulate than (1.2). Also, using (2.10), we see that for i ∈ {0, 1} we have that

v i = (u + r i (c) ϕ (ρ) ∂ x ρ) = u + r i (c) ρϕ ( ρ) ∂ m ρ = u + r i (c) ∂ m ψ ( ρ) .
(2.11)

The mass equation rewrites in terms of the so-called specific volume:

∂ ∂t 1 ρ = ∂ m u.
(2.12)

Basic energy estimates and some properties of regular solutions

The natural energy functional associated to system (2.6) is for any t > 0

R u 2 2 + e ( ρ) + c 2 (∂ m ψ ( ρ)) 2 (t, x)dx + t 0 R ρµ ( ρ) (∂ m u) 2 dsdx ≤ E c (ρ 0 , u 0 ) , (2.13)
where E c (ρ 0 , u 0 ) is defined in (1.16) while e is defined by the relation (1.18). Obviously, the functional in (2.13) is nothing other than the Eulerian energy transformed via the formulae (2.9) and (2.10). Moreover we see that

R v 2 0 2 + e( ρ) (t, x)dx + (1 -r 0 (c)) t 0 R ρµ ( ρ) (∂ m v 0 ) 2 dsdx (2.14) ≤ R 1 2 v 2 0 + e( ρ 0 ) (x)dx ≤ 2E c (ρ 0 , u 0 ) . (2.15)
Remark 2.3.1 In the rest of the paper, in order to ease the reading we will rather use the notation E c instead of E c (ρ 0 , u 0 ) .

Using the uniform bounds ensured by the energy inequality (2.13), we deduce the following proposition.

Proposition 2.3.1 Consider ∈ R and θ ∈ (0, 1]. Then,

       1 2 (2 (γ -1) E c + γ) -θ γ-1 ≤ +1 dq ρ θ (t, q) ≤ E c + γ γ-1 1 θ , +1 ρ γ-1 (t, q) dq ≤ (γ -1) E c + γ.
(2.16)

Proof of Proposition 2.3.1:

The upper bound is easy to obtain as we can write that

+1 1 ρ ≤ +1 ρ γ-1 γ -1 + 1 ρ = +1 e ( ρ) + γ γ -1 ≤ E c + γ γ -1 .
Obviously, using Jensen's inequality, we have that

+1 1 ρ θ ≤ +1 1 ρ θ ,
from which we deduce the right hand side of the first inequality in (2.16). The lower bound is deduced in the following way: we suppose now that there exists ε > 0 such that +1 dq ρ θ (t, q) ≤ ε.

(2.17)

We set now:

A = m ∈ [ , + 1] : ρ (t, m) ≥ 1 (2ε) 1 θ .
Then using (2.17) we deduce that:

2ε| c A| ≤ c A dq ρ θ (t, q) ≤ +1 dq ρ θ (t, q) ≤ ε.
It implies that

meas m ∈ [ , + 1] : ρ (t, m) ≥ 1 (2ε) 1 θ ≥ 1 2 .
But then, we have that

E c ≥ +1 e ( ρ) ≥ 1 2 1 γ -1 1 2ε γ-1 θ - γ γ -1 , which implies that ε ≥ 1 2 (2 (γ -1) E c + γ) -θ γ-1 .
This concludes the proof of Proposition 2.3.1 for the first estimate in (2.16). The second one is a direct consequence of the energy estimate (2.13).

Proposition 2.3.2 There exists a constant C (E c ) that depends only on the initial energy E c such that the following holds true. For any t > 0 and ∈ R there exists a point m (t, )

∈ [ , + 1] such that sup s∈[0,t] Ξ ( ρ (s, m (t, ))) ≤ sup m∈[ , +1] Ξ (ρ 0 (m)) + C (E c ) t.
Proof of Proposition 2.3.2: We recall that the mass equation is

∂ t ρ (t, m) = -ρ 2 (t, m) ∂ m u
and as such, we get that

Ξ ( ρ (s, m)) = Ξ ( ρ 0 (m)) - s 0 ρ γ-1 2 ( ρµ ( ρ)) 1 2 ∂ m u.
From the last identity we get that sup

s∈[0,t] Ξ ( ρ (s, m)) ≤ Ξ ( ρ 0 (m)) + t 0 ρ γ-1 2 ( ρµ ( ρ)) 1 2 |∂ m u| .
Consider ∈ R and integrate the above relation and using the second inequality from Proposition 2.3.1 we obtain

+1 sup s∈[0,t] Ξ ( ρ (s, m)) dm ≤ +1 Ξ ( ρ 0 (m)) dm + t 0 +1 ρ γ-1 2 ( ρµ ( ρ)) 1 2 |∂ m u| ≤ +1 Ξ ( ρ 0 (m)) dm + t 0 +1 ρ γ-1 + t 0 +1 ρµ ( ρ) |∂ m u| 2 ≤ sup m∈[ , +1] Ξ (ρ 0 (m)) + t ((γ -1) E c + γ) + E c .
From this we infer that +1 sup s∈[0,t] Ξ ( ρ (s, m)) dm is controlled by a linear in time function the coefficients of which are controlled uniformly with respect to the initial data. It then follows that there is a point m (t, ) ∈ [ , + 1] where sup

s∈[0,t] Ξ ( ρ (s, m (t, ))) ≤ sup m∈[ , +1] Ξ (ρ 0 (m)) + C (E c ) t.
This concludes the proof of Proposition 2.3.2.

Corollary 2.3.1 Suppose that µ (ρ) verifies hypothesis (µ 5 ). There exists a constant C (E c , ρ 0 L ∞ , t) that depends only on the initial energy E c , ρ 0 L ∞ and time such that the following holds true. For any t > 0 and any ∈ R there exists a point m (t, )

∈ [ , + 1] such that sup s∈[0,t] ρ (s, m (t, )) ≤ C (E c , ρ 0 L ∞ , t) .
Obviously, 2.3.1 is a consequence of Proposition 2.3.2 and of the hypothesis (µ 5 ) assumed for the function Ξ.

In the next section we show how to use these two basic estimates in order to obtain upper and lower bounds for the density.

Upper bound for the density

The first step is to recover an upper bound for the density that does not degenerate when the capillarity coefficient c goes to 0. It transpires that for this purpose we cannot use the so called BD-entropy.

Proposition 2.4.1 Assume that µ verifies the hypothesis (µ 1 ), (µ 2 ), (µ 3 ). Consider ( ρ, u) regular enough solutions for the (2.6) system verifying the energy estimate:

R u 2 2 + e ( ρ) + c 2 (∂ m ψ ( ρ)) 2 (t, x)dx + t 0 R ρµ ( ρ) (∂ m u) 2 dsdx ≤ E c .
Also, we assume that

ρ |t=0 = ρ 0 ∈ L ∞ (R) .
1. The general case: Assume that µ also verifies (µ 5 ). Then, there exists a continuous function Remark 2.4.2 In the case of the Navier-Stokes system for coefficients of the form µ (ρ) = ρ α , the only restriction that we need is 0 < α ≤ γ.

C : [0, ∞) × [0, ∞) × [0, ∞) → [0, ∞)
Remark 2.4.3 Owing to the fact that ρ (t, x) = ρ (t, m (t, x)) , with (2.8), we also have that

ρ (t, x) ≤ C (t, E c , ρ 0 L ∞ ) . (2.19)
Remark 2.4.4 The proof of the above proposition is rather technical and uses a lot of tricks inspired by earlier works of the Russian school of PDEs, see [START_REF] Ya | On a model system of equations of one-dimensional gas motion[END_REF], [START_REF] Vasil | Unique global solution with respect to time of initial-boundary value problems for one-dimensional equations of a viscous gas[END_REF] or the monography [START_REF] Stanislav Nikolaevich Antontsev | Boundary value problems in mechanics of nonhomogeneous fluids[END_REF] for a systematic treatment of the NS system with various boundary conditions. Loosely speaking, the idea is to work with a primitive of the function defining the velocity's equation. Taking in account that the viscosity term can be written as the time derivative of ψ ( ρ), in the case of the NS system, one obtains a dumped differential equation for ψ ( ρ) while the source terms are integrals of the unknowns for which it turns out that we can provide estimates using only E c . The same conclusion holds true for the NSK system although things become more complicated. In this case, ψ ( ρ) verifies a parabolic equation such that we need to use a maximum principle. The fact that we work with functions defined on the whole real line adds another difficulty to the proof because we need to localize the arguments. We point out that Proposition 2.3.2 turns out to be crucial in order to carry out our proof.

Proof of Proposition 2.4.1: Consider m, q, ∈ R fixed arbitrarily. We are going to use the following identity which is obtained by integrating the mass-Lagrangian velocity's equation (2.6) between m and q :

d dt m q u + -ρµ ( ρ) ∂ m u + ρ γ -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) m q = 0 (2.20)
We put the above relation under the form

d dt m u + -ρµ ( ρ) ∂ m u + ρ γ -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) (t, m) = d dt q u + -ρµ ( ρ) ∂ m u + ρ γ -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) (t, q) ,
Next we multiply the above equation with 1 ρ (t, q) and we integrate between and + 1 w.r.t. q. We infer that

+1 dq ρ (t, q) d dt m u + -ρµ ( ρ) ∂ m u + ρ γ -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) (t, m) = +1 1 ρ (t, q) d dt q u dq + +1 -µ ( ρ) ∂ m u + ρ γ-1 -cµ ( ρ) ∂ 2 mm ψ ( ρ) (t, q) dq.
Now, using (2.12) and a simple integration by parts gives us:

+1 1 ρ (t, q) d dt q u dq = d dt +1 1 ρ (t, q) q u dq - +1 d dt 1 ρ (t, q) q u dq = d dt +1 1 ρ (t, q) q u dq - +1 ∂ m u (t, q) q u dq = d dt +1 1 ρ (t, q) q u dq -u (t, + 1) +1 u (t, q) dq + +1 u 2 (t, q) dq.
The above identity is rearranged to give

+1 dq ρ (t, q) d dt m u + -ρµ ( ρ) ∂ m u + ρ γ -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) (t, m) = d dt +1 1 ρ (t, q) q u dq -u (t, + 1) +1 u (t, q) dq + +1 u 2 (t, q) dq + +1 -µ ( ρ) ∂ m u + ρ γ-1 -cµ ( ρ) ∂ 2 mm ψ ( ρ) (t, q) dq.
Recalling that c = r 1 (1 -r 1 ), and using that

-ρµ ( ρ) ∂ m u -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) = -r 1 ρµ ( ρ) ∂ m u -(1 -r 1 ) ρµ ( ρ) ∂ m u -c ρµ ( ρ) ∂ 2 mm ψ ( ρ) = d dt r 1 ψ ( ρ) -(1 -r 1 ) ρµ ( ρ) ∂ m ( u + r 1 ∂ m ψ ( ρ)) = d dt r 1 ψ ( ρ) -(1 -r 1 ) ρµ ( ρ) ∂ mm m u + r 1 ψ ( ρ) ,
the left hand side of the above identity is arranged as to obtain

d dt m u + r 1 ψ( ρ) (t, m) -(1 -r 1 ) ρµ ( ρ) ∂ 2 mm m u + r 1 ψ( ρ) + ρ γ (t, m) = +1 dq ρ (t, q) -1 d dt +1 1 ρ (t, q) q u dq -u (t, + 1) +1 u (t, q) dq + +1 u 2 -µ ( ρ) ∂ m u + ρ γ-1 -cµ ( ρ) ∂ 2 mm ψ ( ρ) (t, q) dq (2.21)
From (2.12), we get

+1 dq ρ (t, q) -1 d dt +1 1 ρ (t, q) q u dq = d dt    +1 dq ρ (t, q) -1 +1 1 ρ (t, q) q u dq    + +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 +1 d dt dq ρ (t, q) = d dt    +1 dq ρ (t, q) -1 +1 1 ρ (t, q) q u dq    + +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 ( u (t, + 1) -u (t, )) (2.22)
Owing to (2.21) and (2.22) we get that

d dt m u + r 1 ψ( ρ) (t, m) -(1 -r 1 ) ρµ ( ρ) ∂ 2 mm m u + r 1 ψ( ρ) (t, m) + ρ γ (t, m) = d dt    +1 dq ρ (t, q) -1 +1 1 ρ (t, q) q u dq    + +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 ( u (t, + 1) -u (t, )) + +1 dq ρ (t, q) -1 -u (t, + 1) q u (t, q ) dq + +1 u 2 -µ( ρ)∂ m u + ρ γ-1 -cµ( ρ)∂ 2 mm ψ( ρ) (t, q) dq
Now, fix N ∈ Z and recalling that the above identity holds for all ∈ R, we integrate w.r.t. on [N, N + 1] in order to obtain that

d dt N +1 N m u d + r 1 ψ( ρ) (t, m) -(1 -r 1 ) ρµ ( ρ) ∂ 2 mm N +1 N m u d + r 1 ψ( ρ) (t, m) + ρ γ (t, m) = d dt    N +1 N +1 dq ρ (t, q) -1 +1 1 ρ (t, q) q u dq d    + N +1 N   +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 ( u (t, + 1) -u (t, ))   d + N +1 N +1 dq ρ (t, q) -1 -u (t, + 1) q u (t, q ) dq + +1 u 2 -µ( ρ)∂ m u + ρ γ-1 -cµ( ρ)∂ 2 mm ψ( ρ) (t, q) dq d (2.23)
Observe that the term appearing in the RHS depends only on time. In the following lines, we analyze the different terms appearing in the right hand side of the previous inequality and we want to prove that (2.23) implies that

d dt N +1 N m ud + r 1 ψ( ρ) (t, m) -(1 -r 1 ) ρµ ( ρ) ∂ 2 mm N +1 N m ud + r 1 ψ( ρ) (t, m) + ρ γ (t, m) ≤ dΓ dt (t) ,
where Γ is such that

Γ (t) ≤ C (E c ) (1 + t) , (2.24)
for some constant that depends only on E c . Let us write that

cµ( ρ)∂ 2 mm ψ( ρ) = r 1 µ( ρ)r 0 ∂ 2 mm ψ( ρ) = -r 1 µ( ρ)∂ m u + r 1 µ( ρ)∂ m v 0
such that we get using Proposition 2.3.1

N +1 N +1 u 2 -µ( ρ)∂ m u + ρ γ-1 -cµ( ρ)∂ 2 mm ψ( ρ) (t, q) dqd = N +1 N +1 u 2 + ρ γ-1 (t, q) dq -(1 -r 1 ) N +1 N +1 µ( ρ)∂ m ud -r 1 N +1 N +1 µ( ρ)∂ m v 0 d ≤ (γ + 1)E c + γ + (1 -r 1 ) N +1 N +1 |µ( ρ)∂ m u| d + r 1 N +1 N +1 |µ( ρ)∂ m v 0 | d . (2.25)
Using Proposition 2.3.1 and hypothesis (µ 3 ) we get that

+1 dq ρ (t, q) -1 +1 |µ( ρ)∂ m u| ≤ +1 dq ρ (t, q) -1 +1 ρ -1 µ ( ρ) 1 2 +1 ρµ( ρ)(∂ m u) 2 1 2 ≤ C (E c , γ) + +1 ρµ( ρ)(∂ m u) 2 .
Thus, we get that

N +1 N +1 dq ρ (t, q) -1 +1 |µ( ρ)∂ m u| d ≤ C (E c , γ) + N +1 N +1 ρµ( ρ)(∂ m u) 2 d
We may treat in the same manner the last term appearing in (2.25). Thus we get that

N +1 N    +1 dq ρ (t, q) -1 +1 u 2 -µ( ρ)∂ m u + ρ γ-1 -cµ( ρ)∂ 2 mm ψ( ρ) (t, q) dq    d ≤ C (E c , γ) + N +1 N +1 ρµ( ρ)(∂ m u) 2 d + N +1 N +1 ρµ( ρ)(∂ m v 0 ) 2 d . (2.26)
Next, using Proposition 2.3.1 we get that

+1 dq ρ (t, q) -1 u (t, + 1) +1 u (t, q ) dq ≤ 4 u (t, + 1) E 1 2 c (2 (γ -1) E c + γ) 1 γ-1 ≤ C (E c , γ) | u (t, + 1)| , (2.27) respectively +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 ( u (t, + 1) -u (t, )) ≤ C (E c , γ) (| u (t, + 1)| + | u (t, )|) .
(2.28)

We thus get that

N +1 N    +1 dq ρ (t, q) -1 -u (t, + 1) +1 u (t, q ) dq    d + + N +1 N    +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 ( u (t, + 1) -u (t, ))    d ≤ C (E c , γ) N +1 N (| u (t, + 1)| + | u (t, )|) d ≤ C (E c , γ) .
(2.29)

Obviously, for all t ≥ 0 we have that

+1 dq ρ (t, q) -1 +1 1 ρ (t, q) q udq dq ≤ C (E c , γ) .
(2.30)

Gathering the estimates (2.26), (2.29), (2.30) along with the maximum principle we see that

∂ ∂t N +1 N m ud + r 1 ψ( ρ) (t, m) -(1 -r 1 ) ρµ ( ρ) ∂ 2 mm N +1 N m ud + r 1 ψ( ρ) (t, m) + ρ γ (t, m) ≤ dΓ dt (t) (2.31)
We denote by

Θ (t, m) = N +1 N m ud + r 1 ψ( ρ) (t, m) -Γ (t)
and we see that (2.31) implies that

∂Θ ∂t (t, m) -(1 -r 1 ) ρµ ( ρ) ∂ 2 mm Θ (t, m) ≤ 0. Recall the conclusion of Corollary 2.3.1, let m 1 (t, N ) ∈ [N -1, N ] and m 2 (t, N ) ∈ [N + 1, N + 2] such that i ∈ {1, 2} ∀s ∈ [0, t] : ρ (s, m i (t, N )) ≤ C (t, E c , ρ 0 L ∞ ) . It then follows that Θ (s, m i (t, N )) ≤ C (t, E c , ρ 0 L ∞ ) (2.32)
Using a maximum-principle, we aim at showing that Θ is bounded on the whole interval ∀m

∈ [m 1 (t, N ) , m 2 (t, N )], namely ∀s ∈ [0, t] : Θ (s, m) ≤ max sup [N -1,N +2] Θ (0, m) , C (t, E c , ρ 0 L ∞ ) . consider A = s ∈ (0, t) : sup m∈[m1(t,N ),m2(t,N )] Θ (s, m) > max sup [N -1,N +2] Θ (0, m) , C (t, E c , ρ 0 L ∞ ) (2.33)
where C (t, E c , ρ 0 L ∞ ) is the same constant as in (2.32). The function

s → sup m∈[m1(t,N ),m2(t,N )] Θ (s, m)
is continuous on [0, t] and therefore the set A appearing in (2.33) is open and therefore it is the union of a at most countable union of open intervals I. Now, there are two ingredients that put together imply that the set A is empty:

In each of the endpoint of such an interval, say I = (t I,g , t I,d ) we have that sup

m∈[m1(t,N ),m2(t,N )] Θ (t I,g , m) =    sup m∈[m1(t,N ),m2(t,N )] Θ (0, m) if t I,g = 0, max sup [N -1,N +2] Θ (0, m) , C (t, E c , ρ 0 L ∞ ) if t I,g > 0 (2.34) while sup m∈[m1(t,N ),m2(t,N )] Θ (t I,d , m) = max sup [N -1,N +2] Θ (0, m) , C (t, E c , ρ 0 L ∞ ) .
This fact and (2.32) imply that for all s ∈ I the supremum over [m

1 (t, N ) , m 2 (t, N )] of Θ (s, m) is achieved in the interior of the segment [m 1 (t, N ) , m 2 (t, N )] i.e. ∀s ∈ I ∃m (s) ∈ (m 1 (t, N ) , m 2 (t, N )) such that Θ (s, m (s)) = sup m∈[m1(t,N ),m2(t,N )]
Θ (s, m) .

Classical considerations lead to the fact that the function s → Θ (s, m (s)) is differentiable almost everywhere on I w.r.t. s and that d ds

[Θ (s, m (s))] = ∂Θ ∂s (s, m (s)) .
The fact that m (s) is a maximum point achieved in the interior of [m 1 (t, N ) , m 2 (t, N )] implies that

∂ 2 mm Θ (s, m (s)) ≤ 0.
Using the above relations we see that s → Θ (s, m (s)) is non-increassing on I and therefore it is controlled by the value achieved in t I,g which is, in turn controlled by max sup

[N -1,N +2] Θ (0, m) , C (t, E c , ρ 0 L ∞ ) , see (2.34).
We thus obtain a contradiction of the very definition of I. The contradiction comes from the fact that we assumed the set A defined in (2.33) in non-empty. Therefore, for all s ∈ [0, t] and all m ∈ [m 1 (t, N ) , m 2 (t, N )] we have that

Θ (s, m) ≤ max sup [N -1,N +2] Θ (0, m) , C (t, E c , ρ 0 L ∞ ) . Since [N, N + 1] ⊂ [m 1 (t, N ) , m 2 (t, N )]
and N ∈ Z and t > 0 where chosen arbitrarily we conclude that

Θ (t, m) ≤ C (t, E c , ρ 0 L ∞ ) for all m ∈ R.
Owing to the definition of Θ (t, m), the estimate (2.24) on Γ, the energy estimates and the hypothesis (µ 1 ) and (µ 2 ) we conclude that there exists a constant such that

ρ (t, m) ≤ C (t, E c , ρ 0 L ∞ ) for all m ∈ R. (2.35)
In order to prove the second part of Proposition 2.4.1, let us recall that relation (2.23) when c = 0 reads

d dt N +1 N m u d + r 1 ψ( ρ) (t, m) + ρ γ (t, m) = d dt    N +1 N +1 dq ρ (t, q) -1 +1 1 ρ (t, q) q u dq d    + N +1 N   +1 1 ρ (t, q) q u dq +1 dq ρ (t, q) -2 ( u (t, + 1) -u (t, ))   d + N +1 N +1 dq ρ (t, q) -1 -u (t, + 1) q u (t, q ) dq + +1 u 2 -µ( ρ)∂ m u + ρ γ-1 (t, q) dq d (2.36)
for all N ∈ Z and for all m ∈ [N, N + 1]. Hypothesis (µ 5 ) is needed in order to have the validity of Proposition 2.3.2 in order to apply the maximum principle. However, as one can see from relation (2.36), the situation for the NS system is simpler: we just integrate in time 2.3.2 and proceed exactly as in estimates (2.26),(2.27), (2.30) in order to conclude that the validity of (2.35).

Bounds for the effective velocities

In the following, we will deduce a one sided inequality concerning v 1 , the equation of which can be put under the form

∂ t v 1 -(1 -r 1 ) ∂ m ( ρµ ( ρ) ∂ m v 1 ) + γ ρ γ r 1 µ ( ρ) ( v 1 -u) = 0.
(2.37) Lemma 2.5.1 Assume that µ satisfies the hypothesis (µ 1 ), (µ 2 ), (µ 3 ) and (µ 6 ). Then, the following estimates hold true:

1. In the general NSK case: we also suppose that (µ 5 ) holds true. Then, we have that for all m ∈ R

v 1 (t, m) ≤ max sup m∈R v 1 (0, m) , 0 + (1 + t) C (E c ) + C (E c , γ) t 0 sup q∈R Φ 1 ρ (s, q) +∞ -∞ ρ (s) µ ( ρ (s)) |∂ m u (s)| 2 1 2 ds 2.
In the case of the NS system: we do not need to impose (µ 5 ) and, moreover, the following more precise statement holds true for all m ∈ R:

v 1 (t, m) ≤ v 1 (0, m) + (1 + t) C (E c ) + C (E c , γ) t 0 sup q∈R Φ 1 ρ (s, q) +∞ -∞ ρ (s) µ ( ρ (s)) |∂ m u (s)| 2 1 2
ds Corollary 2.5.1 Assume that µ satisfies the hypothesis (µ 1 ), (µ 2 ), (µ 3 ) and (µ 6 ). Then, the following estimates hold true:

1. In the general NSK case: we also suppose that (µ 5 ) holds true. For all x ∈ R we have that

v 1 (t, x) = u (t, x) + r 1 ∂ x ϕ (ρ) (t, x) ≤ max sup x∈R (u 0 (x) + r 1 ∂ x ϕ (ρ 0 (x))) , 0 + C (E c ) (1 + t) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (s, z) +∞ -∞ µ (ρ (s, z)) |∂ x u (s, z)| 2 dz 1 2
ds.

(2.38)

2. In the case of the NS system: we do not need to impose (µ 5 ) and for all x ∈ R we have that

v 1 (t, x) = u (t, x) + r 1 ∂ x ϕ (ρ) (t, x) ≤ u 0 (x) + r 1 ∂ x ϕ (ρ 0 ) (x) + C (E c ) (1 + t) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (s, z) +∞ -∞ µ (ρ (s, z)) |∂ x u (s, z)| 2 dz 1 2
ds.

(2.39)

Remark 2.5.1 We mention that this is the key estimate leading to the fact that the density is bounded away from vacuum. In some sense, in the same way as working with a primitive for the momentum equation leads to an upper bound for the density, working with a primitive of the equation of the effective velocity will enable us to obtain a lower bound for the density.

Remark 2.5.2 The local nature of the estimate (2.39) with respect to the initial data allows us to show that, in the case of the NS system, the density is bounded away from vacuum without using any information on the derivatives of ρ 0 .

Proof of Lemma 2.5.1: Let us notice that for all t > 0 and m, q ∈ [N, N + 1] we have that

u (t, m) -u (t, q) = m q ∂ m u ≤ m q 1 ρµ ( ρ) 1 2 m q ρµ ( ρ) |∂ m u| 2 1 2 = m q 1 ρ Λ 1 ρ 1 2 m q ρµ ( ρ) |∂ m u| 2 1 2 ≤ m q 1 ρ 1 2 m q ρµ ( ρ) |∂ m u| 2 1 2 1 + sup ∈R Φ 1 ρ (t, ) , (2.40)
where we recall that Λ and Φ are defined in (1.23) and that hypothesis (µ 6 ) reads:

∀τ > 0 : Λ (τ ) = 1 µ 1 τ ≤ (1 + Φ (τ )) 2 .
Integrating the relation (2.40) w.r.t. q on [N, N + 1] and using Proposition 2.3.1 yields

u (t, m) ≤ N +1 N | u (t, q)| dq + 1 + sup q∈R Φ 1 ρ (t, q) N +1 N m q 1 ρ 1 2 m q ρµ ( ρ) |∂ m u| 2 1 2 dq ≤ C (E c ) + C (E c , γ) 1 + sup q∈R Φ 1 ρ (t, q) +∞ -∞ ρµ ( ρ) |∂ m u| 2 1 2
.

(2.41)

Since m ∈ [N, N + 1], the RHS of the above estimate is independent of N and N was arbitrarly fixed we conclude that (2.41) holds true for all m ∈ R. Using (2.41) in (2.37) along with the fact that ρ is bounded and (µ 3 ), we obtain that

∂ t v 1 -(1 -r 1 ) ∂ m ( ρµ ( ρ) ∂ m v 1 ) + γ ρ γ r 1 µ ( ρ) v 1 ≤ C (E c ) + C (E c , γ) 1 + sup q∈R Φ 1 ρ (t, q) +∞ -∞ ρµ ( ρ) |∂ m u| 2 1 2
.

(2.42)

We denote by

Ψ (t) = tC (E c ) + C (E c , γ) t 0 1 + sup q∈R Φ 1 ρ (s, q) +∞ -∞ ρµ ( ρ) |∂ m u| 2 1 2 ds,
the primitive vanishing at t = 0 of the RHS term appearing in (2.42). The inequality (2.42) can be put under the following form:

d dt ( v 1 -Ψ) -(1 -r 1 ) ∂ m ( ρµ ( ρ) ∂ m ( v 1 -Ψ)) + γ ρ γ r 1 µ ( ρ) ( v 1 -Ψ) ≤ - γ ρ γ r 1 µ ( ρ) Ψ ≤ 0.
We remark that owing to the fact that

R v 2 1 (t, m) dm < ∞,
we have that for all t ≥ 0: lim

|m|→∞ ( v 1 (t, m) -Ψ (t)) = -Ψ (t) ≤ 0.
The maximum principle then implies that

v 1 (t, m) -Ψ (t) ≤ max sup q∈R v 1 (0, q) , sup s∈[0,t]
(-Ψ (s)) ≤ max sup q∈R v 1 (0, q) , 0 and accordingly we get that

v 1 (t, m) ≤ max sup q∈R v 1 (0, q) , 0 + Ψ (t) ≤ max sup q∈R v 1 (0, q) , 0 + C (t, E c ) + C (t, E c ) t 0 sup q∈R Φ 1 ρ (s, q) +∞ -∞ ρ (s) µ ( ρ (s)) |∂ m u (s)| 2 1 2 ds.
(2.43)

In order to prove the second part of Lemma 2.5.1, we remark that in the case of the Navier-Stokes system i.e. r 1 (c) = 1, with c = 0 we can get a much more precise estimate since, in this case, (2.37) does not have a diffusion operator. Thus, we see that there is no need to invoke the maximum principle and the estimate is localized: for all m ∈ R we get using (2.37), (2.35), (µ 3 ), (2.41) that

v 1 (t, m) ≤ v 1 (0, m) + t 0 γ ρ γ (s, m) r 1 µ ( ρ (s, m)) u (s, m) exp - t s γ ρ γ µ ( ρ) ds ≤ v 1 (0, m) + C (t, E c ) t 0 u (s) L ∞ ds ≤ v 1 (0, m) + C (t, E c ) + C (E c ) t 0 1 + sup q∈R Φ 1 ρ (s, q) +∞ -∞ ρ (s) µ ( ρ (s)) |∂ m u (s)| 2 1 2 ds ≤ v 1 (0, m) + C (t, E c ) + C (E c , γ) t 0 sup q∈R Φ 1 ρ (s, q) +∞ -∞ ρ (s) µ ( ρ (s)) |∂ m u (s)| 2 1 2 ds.
This concludes the proof of Lemma 2.5.1

The result announced in Corollary is obtained just by returning back to the Eulerian variables, taking in consideration that for all x ∈ R it holds true that

v 1 (t, x) = v 1 (t, m (t, x))
where m (t, x) is the inverse of the function m → X(t, Y (m)).

In the next section we show how to use the estimates provided by Corollary 2.5.1 in order to obtain the fact that the density is bounded away from vacuum.

Lower bound for the density

In the following, we show how it is possible to find a lower bound for the density with the use of the inequality (2.43). At this point, we need an Eulerian-equivalent for Proposition 2.3.1. We claim that: Proposition 2.6.1 Consider ∈ R. Then, for all t ≥ 0 there exists a point y (t, ) ∈ [ , + 2E 0 ] such that ρ (t, y (t, )) ≥ π -1 (1/2) , where we recall that π (ρ) = ρe (ρ).

Proof of Proposition 2.6.1. Recall that according to the energy inequality for all t ≥ 0 we have +2E0 ρ (t, x) e (ρ (t, x)) dx ≤ E 0 .

Assume that for all x ∈ [ , + 2E 0 ] one has ρ (t, x) < ε, with 0 < ε ≤ 1 then, as the function π is strictly decreasing on [0, 1] we find that π (ρ (t, x)) ≥ π (ε) .

Integrating the previous inequality on [ , + 2E 0 ] yields

E 0 ≥ 2π (ε) E 0 , such that ε ≥ π -1 (1/2) ,
with π -1 (1/2) is here the unique element α of [0, 1] such that π(α) = 1 2 . The conclusion follows.

We are now in the position of obtaining a lower bound for the density.

Proposition 2.6.2 Assume that µ satisfies the hypothesis (µ 1 ), (µ 2 ), (µ 3 ), (µ 4 ), (µ 6 ) and (µ 7 ).

1. In the general NSK case we also suppose that (µ 5 ) holds true. Suppose that a.e. x ∈ R : u 0 (x) + ∂ x ϕ (ρ 0 (x)) ≤ M 0 , with M 0 ∈ R. Then, we have that for all x ∈ R:

C t, E c , M 0 , (ρ 0 , 1 ρ 0 ) L ∞ ≤ ρ (t, x) .
(2.44)

2. In the case when c = 0, i.e. for the NS system, we have that for all x ∈ R:

C t, E c , (ρ 0 , 1 ρ 0 ) L ∞ ≤ ρ (t, x) .
(2.45)

Remark 6 At this level, it is important to mention that we can obtain similar estimates if we assume that

M 1 ≤ u 0 + r 1 (c) ∂ x ϕ(ρ 0 )
for some M 1 ∈ R. It suffices to repeat the same argument except that we take y ∈ [N -2E 0 (c), N ] and we integrate between x and y.

Remark 7 In the case of the NSK system, we also have that

u (t, x) + r 1 (c) ∂ x ϕ (ρ (t, x)) ≤ C t, E c , M 0 , (ρ 0 , 1 ρ 0 ) L ∞ (2.46)
Proof of Proposition 2.6.2: As before, let N ∈ Z, consider any x ∈ [N, N + 1] and y ∈ [N + 1, N + 1 + 2E c ] (we have in particular x ≤ y) and integrate (2.38) between x and y in order to obtain that y x u (t, z) dz+r 1 (ϕ (ρ(t, y)) -ϕ (ρ (t, x))

≤ (y -x) max(M 0 , 0) + C (E c ) (1 + t) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (s, z) +∞ -∞ µ (ρ (s, z)) |∂ x u (s, z)| 2 dz 1 2 ds ≤ (2E c + 1) max(M 0 , 0) + C (E c ) (1 + t) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (s, z) +∞ -∞ µ (ρ (s, z)) |∂ x u (s, z)| 2 dz 1 2
ds.

(2.47) Using Proposition 2.6.1, there exists y

N (t) ∈ [N + 1, N + 1 + 2E c ] with the property that ρ (t, y N (t)) ≥ π -1 (1/2) .
We take y = y N (t) in (2.47) and recalling that

Φ (τ ) = -ϕ 1 τ we get that r 1 Φ 1 ρ (t, x) ≤r 1 Φ 1 ρ (t, y N (t)) + (2E c + 1) max(M 0 , 0) + C (E c , γ) (1 + t) - y N (t) x u (t, z) dz + C (E c , γ) t 0 sup z∈R Φ 1 ρ (t, z) +∞ -∞ µ (ρ) |∂ x u| 2 1 2 .
(2.48)

The only term that needs to be treated in (2.48) is the integral of u over [x, y N (t)]. This is done in the following lines using hypothesis (µ 7 ) we infer that for all ε > 0 we have:

-

y N (t) x u (t, z) dz = - y N (t) x 1 ρ (t, z) ρ (t, z)u (t, z) dz ≤ y N (t) x 1 ρ (t, z) ρ (t, z)u (t, z) dz ≤ y N (t) x 1 + Φ 1 ρ (t, z) 1-η ρ (t, z)u (t, z) dz ≤ 1 + sup z∈R Φ 1 ρ (t, z) 1-η y N (t) x ρ (t, z)u (t, z) dz ≤ 1 + sup z∈R Φ 1 ρ (t, z) 1-η (2E c + 1) 1 2 E c ≤ C (E c , η) 4ε + ε sup z∈R Φ 1 ρ (t, z)
.

Thus, we get that for all x ∈ R and all ε > 0 we have that

r 1 Φ 1 ρ (t, x) ≤ C t, E c , M 0 , (ρ 0 , 1 ρ 0 ) L ∞ + C (E c , η) 4ε + ε sup z∈R Φ 1 ρ (t, z) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (t, z) +∞ -∞ µ (ρ) |∂ x u| 2 1 2 .
(2.49)

Taking the supremum w.r.t. x ∈ R in the LHS, taking ε = r 1 /2 and using Grönwall's inequality we obtain that

Φ 1 ρ (t, x) ≤ C t, E c , M 0 , (ρ 0 , 1 ρ 0 ) L ∞ .
Thus using the hypothesis (µ 4 ) we get that ρ is bounded below and that we have for any (t,

x) ∈ R + × R C t, E c , M 0 , (ρ 0 , 1 ρ 0 ) L ∞ ≤ ρ (t, x) .
(2.50)

The second part of Proposition 2.6.2 follows owing to the fact that in the case of the Navier-Stokes system, a more precise information can be obtained (we recall that r 1 = 1 in this case). Recall that according to (2.39), for all x ∈ R we have that

v 1 (t, x) = u (t, x) + r 1 ∂ x ϕ (ρ) (t, x) ≤ u 0 (x) + r 1 ∂ x ϕ (ρ 0 ) (x) + C(E c , γ) (1 + t) (2.51) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (t, z) +∞ -∞ µ (ρ) |∂ x u| 2 1 2 .
(2.52)

Choosing as above x ∈ [N, N + 1] , y N (t) ∈ [N + 1, N + 1 + 2E c ] with ρ (t, y N (t)) ≥ π -1 (1/2)
and integrating we get that

y N (t) x u (t, z) dz + r 1 Φ 1 ρ (t, x) -r 1 Φ 1 ρ (t, y N (t)) ≤ y N (t) x (u 0 (z) + r 1 ∂ x ϕ (ρ 0 ) (z)) dz C (E c ) (1 + t) + C (E c , γ) t 0 sup x∈R Φ 1 ρ (t, x) +∞ -∞ µ (ρ) |∂ x u| 2 1 2 ≤ y N (t) x u 0 (z) dz + r 1 (ϕ (ρ 0 (y N (t)) -ϕ(ρ 0 (x))) + C (E c , γ) t 0 sup z∈R Φ 1 ρ (t, z) +∞ -∞ µ (ρ) |∂ x u| 2 1 2 . We deduce that r 1 Φ 1 ρ (t, x) ≤ r 1 Φ 1 ρ (t, y N (t)) + y N (t) x u 0 (z) dz + r 1 {ϕ(ρ 0 (y N (t))) -ϕ (ρ 0 (x))} - y N (t) x u (t, z) dz + C (E c , γ) t 0 sup x∈R Φ 1 ρ (t, x) +∞ -∞ µ (ρ) |∂ x u| 2 1 2 , from which it yields r 1 Φ 1 ρ (t, x) ≤ r 1 Φ 1 ρ (t, y N (t)) + (2E 0 + 1) 1 2 u 0 L 2 + 2 sup p∈[ 1 ρ 0 L ∞ , ρ0 L ∞ ] |ϕ (p)| - y N (t) x u (t, z) dz + C (E c , γ) t 0 sup x∈R Φ 1 ρ (t, x) +∞ -∞ µ (ρ) |∂ x u| 2 1 2
and since the terms on the RHS can be controlled we use Grönwall's lemma in order to obtain that

Φ 1 ρ (t, x) ≤ C t, E c , (ρ 0 , 1 ρ 0 ) L ∞ , from which we can deduce that C t, E c , (ρ 0 , 1 ρ 0 ) L ∞ ≤ ρ (t, x) .
This concludes the proof of Proposition 2.6.2. The estimate from (2.46) follows by combining the conclusion of Proposition 2.6.2 with the inequality (2.38).

Global existence of the approximate solutions

Now we can come back to the sequence of solution (ρ n , u n ) n∈N of the Navier-Stokes Korteweg system (1.2) and of the compressible Navier-Stokes system (1.1) constructed in the subsection 2.1. We recall that these are defined on a finite time interval (0, T n ) for each n ∈ N.

In view of Proposition 2.4.1 and the estimate (2.19) from Remark 2.4.3 that follows we obtain that:

ρ n (t, •) L ∞ (R) ≤ C (t, E c (ρ n 0 , u n 0 ) , ρ n 0 L ∞ ) < +∞.
(2.53) for any n ∈ N and t ∈ (0, T m ). Next, owing to Proposition 2.6.2, and (2.35) we know that for all t ∈ (0, T n ) and

x ∈ R we have that:

1 ρ n (t, •) L ∞ (R) ≤ C t, E c (ρ n 0 , u n 0 ) , M n 0 , ρ n 0 , 1 ρ n 0 L ∞ , ρ n (t, •) L ∞ (R) ≤ C (t, E c (ρ n 0 , u n 0 ) , ρ n 0 L ∞ ) (2.54)
where v n 1|t=0 ≤ M n 0 , with C a continuous function for the NSK system (there is no M n 0 in the case of the NS system). Of course, owing to subsection 2.1 we know that

E c (ρ n 0 , u n 0 ) , M n 0 , ρ n 0 , 1 ρ n 0 L ∞
are uniformly bounded in n and we thus deduce, that, for a larger C if necessary, we have that

ρ n (t, •) L ∞ (R) ≤ C (t, E c (ρ 0 , u 0 ) , ρ 0 L ∞ ) < +∞. (2.55) and that 1 ρ n (t, •) L ∞ (R) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ (2.56)
The blow-up criterion of the Theorem 2.1 and (2.55), (2.56) imply that for any n ∈ N we have T n = +∞ and that (2.55) and (2.56) hold true in fact for all t ≥ 0.

In the case of the NSK system, for all n ∈ N and all t ≥ 0 we also have

u n + r 1 (c) ∂ x ϕ (ρ n 0 ) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ (2.57)
2.8 Total variation estimates for solutions of the NSK system

Using the estimate (2.56) we deduce immediately for any n ∈ N, an estimate for the total variation of ϕ(ρ n ). Indeed recall that (2.57) states that u n (t, x) + r 1 (c) ∂ x ϕ (ρ n (t, x)) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 ,

1 ρ 0 L ∞ (2.58)
This implies that for all n ∈ N, for all t ≥ 0, the function x →

x 0 u n (t, z) dz + r 1 (c) ϕ (ρ n (t, x)) -C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ x is deacreasing and, taking in account (2.56), it is of bounded variation on any compact. Using the energy estimate along with (2.56) we obtain that:

T V ([-L, L] , ϕ (ρ n (t, x))) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ (1 + L) , (2.59) for all t ≥ 0 and L > 0 with T V ([-L, L] , ϕ (ρ n (t, •))) the total variation of ϕ n (t, •) on [-L, L]. Moreover, we have that

∂ t ϕ(ρ n ) = -u n ∂ x ϕ(ρ n ) + µ (ρ n ) ρ n ∂ x u n such that T 0 L -L |∂ t ϕ(ρ n )| ≤ u n L 1 t L ∞ x sup t∈[0,T ] T V ([-L, L] , ϕ(ρ n (t, •))) + µ 1 2 (ρ n ) ρ n L ∞ (R + ×R) µ 1 2 (ρ n )∂ x u n L 1 ([0,T ]×[-L,L])
Using (2.56), (2.55), (2.59) and the energy estimate, we deduce that ϕ(ρ n ) belongs uniformly to BV loc ([0, ∞) × R) namely ϕ(ρ n ) BV ([0,t]×[-L,L]) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ .

(2.60) From (2.56), (2.55) we deduce that for any T, L > 0 we have:

1 2 L 2 ∂ x (µ (ρ) ∂ x v -ρ γ ) 1 2 L 2 ≤ 1 2 C (t) t 0 σ ρ v 1 2 L 2 ∂ x u L 2 µ (ρ) ∂ x v -(ρ γ -1) 3 2 L 2 ≤ ε 2 t 0 σ ρ v 2 L 2 + C (t) 2ε t 0 ∂ x u 4 3 L 2 σ µ (ρ) ∂ x v 2 L 2 . (3.5)
Using in the identity (3.2) the estimates (3.3), (3.4) and (3.5) along with Gronwall's lemma we get that

t 0 R σρ | v0 | 2 + (1 -r) 2 σ (t) R µ (ρ (t)) (∂ x v 0 (t)) 2 ≤ C (t) . (3.6)
As it is well known from previous works [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF], [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF] we can already recover that

t 0 ∂ x v 0 L ∞ ≤ t 0 σ -1 2 (τ ) dτ 1 2 t 0 σ (t) ∂ x v 0 2 L ∞ 1 2
.

Furthermore for the NSK system we have similar estimate for v 1 but these estimates are not uniform in c > 0.

≤ ε t 0 R σ 2 µ (ρ) (∂ x v) 2 + C (t) .
(3.10) Also, we have for any ε > 0 by using again (2.56), (2.55, Young inequality) and Gagliardo-Nirenberg inequality

t 0 R σ 2 (ρµ (ρ) + µ (ρ)) µ (ρ) (1 -r) (µ (ρ) ∂ x v -(ρ γ -1)) ∂ x u∂ x v ≤ t 0 R σ 2 µ (ρ) ∂ x v -(ρ γ -1) 1 2 L 2 ρ v 1 2 L 2 ∂ x u L 2 ∂ x v L 2 ≤ C (t) t 0 σ 2 ∂ x u 2 L 2 µ (ρ) ∂ x v -(ρ γ -1) L 2 ρ 1 2 v L 2 + ε t 0 R σ 2 µ (ρ) (∂ x v) 2 ≤ C (t) t 0 σ 2 ∂ x u 2 L 2 ρ 1 2 v 2 L 2 + C (t) t 0 σ 2 ∂ x u 2 L 2 µ (ρ) ∂ x v -(ρ γ -1) 2 L 2 + ε t 0 R σ 2 µ (ρ) (∂ x v) 2 (3.11)
Putting together (3.6) , (3.7), (3.8), (3.9) , (3.10) , and (3.11) we get that

1 2 R σ 2 (t) ρ (t) v2 (t) + (1 -r) t 0 R σ 2 µ (ρ) (∂ x v) 2 ≤ C (t) .
The control over 1 ρ L ∞ , A (ρ, u) and B (ρ, u) gives us σ (t)

1 2 ∂ x v 0 (t) L ∞ ≤ C (t) , (3.12)
for any t ≥ 0. The same estimate are again true for v 1 for the NSK system but are not uniform in c > 0. This concludes the proof of Proposition 3.0.1

End of proof of the main results

The sequence of solutions constructed in Section 2.1 (ρ n , u n ) n∈N are globally defined owing to the intermediary results proved in Proposition 2.4.1, Proposition 2.6.2, the estimates from Proposition 3.0.1 and the remark that follows. Using this, we may extract a subsequence that converges to a solution of the NS system (1.1). The fact that the solution itself verifies the estimates announced in Theorem 1.1 is a consequence of the Fatou lemma. We skip the technical details and we refer the reader to [START_REF] Mellet | Existence and uniqueness of global strong solutions for one-dimensional compressible Navier-Stokes equations[END_REF][START_REF] Jiu | The Cauchy problem for 1d compressible flows with density-dependent viscosity coefficients[END_REF] for more details on this subject. It proves in particular the Theorem 1.1.

Regarding Theorem 1.2, we skip the details showing that the approximate solution (ρ n , u n ) n∈N converge up to a subsequence to a unique solution (ρ, u) of the system (1.2) on [0, T ]. Again, the fact that the solution itself verifies the estimates announced in Theorem 1.2 is a consequence of the Fatou lemma. One can consult in particular [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF].

In what Theorem 1.3 is concerned, we only mention that the only not so-classical estimates (2.60) and (2.61) are used in order to justify the fact that the condition u c (t, x) + r 1 (c) ∂ x ϕ (ρ c (t, x)) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ , is obtained also in the limit c → 0, namely

u + ∂ x ϕ (ρ) ≤ C t, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ ,
in the sense of measures. The arguments concerning the weak convergence up to subsequence of the sequence (ρ c , u c ) c>0 are the same as those in the proof of the Theorem 1.3. It is important to point out that in fact (ρ c , u c ) c>0 converges weakly and not only up to a subsequence because (ρ c , u c ) c>0 has a unique accumulation point, indeed the limit (ρ, u) when c goes to 0 is solution of the NS system and this solution is unique as we will see.

From the Hoff type estimates, we can deduce that ∂ x u belongs to L 2 ([0, T ], L ∞ (R)) for any T > 0. The uniqueness of the solution (ρ, u) is then a consequence of the fact that on [0, T ], ∂ x u belongs to L 1 ([0, T 1 ], L ∞ ). We refer to [START_REF] Burtea | New effective pressure and existence of global strong solution for compressible Navier-Stokes equations with general viscosity coefficient in one dimension[END_REF] for the details of the proof in the case of the compressible Navier-Stokes system where we use a Lagrangian formulation. The adaptation of the proof to the Navier-Stokes Korteweg system is direct.

µ

  ρ n BV ([0,T ]×[-L,L]) ≤ C t, L, E c (ρ 0 , u 0 ) , M 0 , ρ 0 , 1 ρ 0 L ∞ . (2.61) with C(T, L) independent on n. (ρ (t)) (∂ x v (t)) 2 + C (t) . (3.4)Finally, let us observe that combining (2.56), (2.55), energy estimate and Young inequality we get for anyε > ρ) + ρµ (ρ)) (µ (ρ)) 2 (µ (ρ) ∂ x v -(ρ γ -1)) 2 ∂ x u ρ) + ρµ (ρ)) (µ (ρ)) 2 |∂ x u| |µ (ρ) ∂ x v -(ρ γ -1)| µ (ρ) ∂ x v -(ρ γ -1) L ∞ ρ) + ρµ (ρ)) (µ (ρ)) 2 |∂ x u| |µ (ρ) ∂ x v -(ρ γ -1)| µ (ρ) ∂ x v -(ρ γ -1)

  which is increasing w.r.t. to the first variable and such thatρ (t, m) ≤ C (t, E c , ρ 0 L ∞ ) .(2.18)2. The NS system: If c = 0, i.e. we are considering the Navier-Stokes system, the same conclusion as in (2.18) holds true without assuming hypothesis (µ 5 ) for the viscosity function µ.

	Remark 2.4.1 The results of Proposition 2.3.1, Proposition 2.3.2 and Corollary 2.3.1 hold true for the pair ( ρ, u)
	of Proposition 2.4.1.

Acknowledgements

CB has been partially funded by the projects SingFlows (ANR-18-CE40-0027-01) and CRISIS (ANR-20-CE40-0020-01) both operated by the French National Research Agency (ANR). BH has been partially funded by the ANR project INFAMIE ANR-15-CE40-0011.

Estimates in the Hoff regularity class

We recall the Eulerian effective velocity's equation with v0 = ∂ t v 0 + u∂ x v 0 :

(3.1)

In the computations that follow, we drop the 0-lowerscript of v 0 and r 0 in order to ease readability. As it is by now well-known, once we have uniform bounds for the density we can obtain some higher order estimates by using the energy-type functionals first considered by Hoff in [START_REF] Hoff | Global existence for 1d, compressible, isentropic Navier-Stokes equations with large initial data[END_REF] (see also [START_REF] Hoff | Global solutions of the equations of one-dimensional, compressible flow with large data and forces, and with differing end states[END_REF]). These are gathered in the following proposition:

Proposition 3.0.1 The following estimates hold true:

where σ (t) = min {1, t}. As consequence we deduce that

Remark 3.0.1 In the case of the Navier-Stokes system (1.1) the functionals A (ρ, u) and B (ρ, u) are bounded by

The first Hoff energy

Multiplying the equation (3.1) with v and after some cumbersome yet straightforward computations we get that

We multiply the previous identity with σ (t) = min {1, t} and, by time integration we get that

Also, we write that

Combining the last two identities we get that

Owing to the basic energy estimate and using (2.56), (2.55), the following inequality holds true:

Next using again energy estimate and ,(2.56), (2.55), we have that

The second Hoff energy

In the following lines, we aim at obtaining higher order estimates. In order to do so, we apply the operator ∂ t + u∂ x to the equation of v 0 (again we drop the 0-lower script) in order to obtain that

Multiplying with v we see that

Next, we have that

Let us deal with the last term

Putting togeather the last two relations, we get that

Next, let us take a look at

We gather the last three relations in order to obtain 0 = 1 2

We end up with 1 2

We multiply the last relation with σ 2 (t) and we integrate in time in order to obtain

So, we see that owing to (3.6) we have that 1 2

The second term in (3.7), we see that


Antonio Brogi

Wolf Zimmermann

Kyriakos Kritikos

Preface

Service-oriented and cloud computing have made a huge impact both on the software industry and on the research community. Today, service and cloud technologies are applied to build large-scale software landscapes as well as to provide single software services to end users. Services today are independently developed and deployed as well as freely composed while they can be implemented in a variety of technologies, quite an important fact from a business perspective. Similarly, cloud computing aims at enabling flexibility by offering a centralized sharing of resources. The industry's need for agile and flexible software and IT systems has made cloud computing the dominating paradigm for provisioning computational resources in a scalable, on-demand fashion. Nevertheless, service developers, providers, and integrators still need to create methods, tools, and techniques to support cost-effective and secure development as well as the use of dependable devices, platforms, services, and service-oriented applications in the cloud.

The European Conference on Service-Oriented and Cloud Computing (ESOCC) is the premier conference on advances in the state of the art and practice of service-oriented computing and cloud computing in Europe. The main objectives of this conference are to facilitate the exchange between researchers and practitioners in the areas of service-oriented computing and cloud computing, as well as to explore the new trends in those areas and foster future collaborations in Europe and beyond. The 8th edition of ESOCC, ESOCC 2020, was held in the city of Heraklion in Crete, Greece, during September 28-30, 2020, under the auspices of FORTH-ICS.

ESOCC 2020 was a multi-event conference aimed to cover both an academic and industrial audience. The main event was associated with the main research track, which focused on the presentation of cutting-edge research in both the service-oriented and cloud computing areas. In conjunction, an industrial track was also held bringing together academia and industry by showcasing the application of service-oriented and cloud computing research, especially in the form of case studies, from industry. Overall, 20 submissions were received, out of which 6 outstanding full and 8 short papers were accepted.

Each submission was peer-reviewed by three main reviewers, either directly from the PC members or their colleagues. Due to the high quality of the manuscripts received, additional discussions were conducted, both among the PC members as well as between the two PC chairs, before the final selection was performed. The PC chairs would like to thank all the reviewers that participated in the reviewing process not only for enabling to increase the quality of the received manuscripts but also for sharing particular ideas on how the respective work, even if rejected in its current form in the ESOCC conference, could be substantially improved.

The attendees of ESOCC had the opportunity to follow two outstanding keynotes that were part of the conference program. The first keynote was conducted by Massimo Villari, Professor and Rector Delegate for ICT as well as Head of the Computer Science School in the University of Messina, Italy. This keynote concerned recent research advances and trends towards realizing the vision of Osmotic Computing. The second keynote was conducted by Joseph Spillner, Head of the Service Prototyping Lab and Associate Professor at Zurich University of Applied Sciences in Switzerland. This second keynote concerned the presentation of methods for developing production-ready, Function-as-a-Service applications concentrating on scalable event-driven data processing that are well-suited for highly dynamic environments with varying loads.

The additional events held in ESOCC 2020 included the PhD symposium, enabling PhD students to present their work in front of real experts, as well as the EU projects track, supplying researchers with the opportunity to present the main research results that they have achieved in the context of currently operating EU projects. Further, ESOCC 2020 included the organization of satellite workshops. All these events were accompanied by respective proceedings which were published separately.

Finally, this 8th edition of ESOCC included a novel track dedicated to the conduction of tutorials. This enabled the workshop participants to get acquainted with the latest results of specific European projects as well as of specific European research groups in a practical manner which included demonstrations of research prototypes.

The PC chairs and the general chair would like to gratefully thank all the persons involved in making ESOCC 2020 a success. This includes both the PC members and their colleagues that assisted in the reviews as well as the organizers of the industry track, the PhD symposium, the EU projects track, and the workshops. A special applause should also go to the members of the Local Organizing Committee for their devotion, willingness, and hospitality. Finally, a special thanks goes to all the authors of all the manuscripts submitted to ESOCC 2020, the presenters of the accepted papers who made interesting and fascinating presentations of their work, as well as the active attendees of the conference who initiated interesting discussions and gave fruitful feedback to the presenters. All these persons not only enabled a very successful organization and execution of ESOCC 2020, but also formulate an active and vibrant community which continuously contributes to the research in service-oriented and cloud computing. This also encourages ESOCC to continue contributing with new research outcomes to further facilitate and enlarge its community as well as have a greater impact and share in both the service-oriented and cloud computing research.

September 2020

Antonio Brogi Wolf Zimmermann Kyriakos Kritikos vi Preface

Organization

ESOCC 2020 was organized by FORTH-ICS, Greece.

Organizing Committee