Georgy Kitavtsev
email: georgy.kitavtsev@mis.mpg.de

Philippe Laurencot
email: laurenco@math.univ-toulouse.fr

Barbara Niethammer
email: niethammer@maths.ox.ac.uk

Philippe Laurençot

Weak solutions to lubrication equations in the presence of strong slippage

published or not. The documents may come

Weak solutions to lubrication equations in the presence of strong slippage

1 Introduction.

During the last thirty years lubrication theory was successfully applied to modeling of dewetting processes in micro and nanoscopic liquid films on a solid polymer substrates see e.g. Oron et al. [START_REF] Oron | Long-scale evolution of thin liquid films[END_REF], Münch et al. [START_REF] Münch | Lubrication models with small to large slip lengths[END_REF] and references therein. The influence of intermolecular interactions that are typically due to the competition between the long-range attractive van der Waals and short-range Born repulsive intermolecular forces play an important role in such processes, see Oron et al. [START_REF] Oron | Long-scale evolution of thin liquid films[END_REF], de Gennes [START_REF] De Gennes | Wetting: Statics and dynamics[END_REF].

Besides intermolecular forces and surface tension at the free surface of the film it has been shown by Fetzer et al. [START_REF] Fetzer | Quantifying hydrodynamic slip: A comprehensive analysis of dewetting profiles[END_REF] that the dewetting of polymer films on hydrophobic substrates also involves such boundary effect as slippage on a solid substrate. The measure of slip is a so-called slip length, which is defined as an extrapolated distance relative to the wall where the tangential velocity component of the liquid vanishes. Recently, it has been shown experimentally and theoretically that the early stages of the dewetting process and the evolving morphology depend markedly on the magnitude of the effective slip length, 1 which can be of the size of the height of the liquid film or even larger for nanoscale systems, see e.g. Fetzer et al. [START_REF] Fetzer | New slip regimes and the shape of dewetting thin liquid films[END_REF], Münch and Wagner [START_REF] Münch | Contact-line instability of dewetting thin films[END_REF]. Recently in Münch et al. [START_REF] Münch | Lubrication models with small to large slip lengths[END_REF], Kargupta et al. [START_REF] Kargupta | Instability, dynamics and morphology of thin slipping films[END_REF] closed-form one-dimensional lubrication equations over a wide range of slip lengths were derived from the underlying equations for conservation of mass and momentum, together with boundary conditions for the tangential and normal stresses, as well as the kinematic condition at the free boundary, impermeability and Navier-slip condition at the liquid-solid interface. Asymptotic arguments, based on the magnitude of the slip length show that within a lubrication scaling there are two distinguished limits, see Münch et al. [START_REF] Münch | Lubrication models with small to large slip lengths[END_REF].

These are the well-known weak-slip model

∂ t h = -∂ x (h 3 + b h 2 )∂ x (σ∂ xx h -Π(h)) (1.1)
with b denoting the slip length parameter, and the strong-slip model

Re (∂ t (hu) + ∂ x (hu 2 )) = 4∂ x (ν(h)∂ x u) + h∂ x (σ∂ xx h -Π(h)) - u β (1.2a) ∂ t h = -∂ x (hu) , (1.2b)
respectively. Here, u(x, t) and h(x, t) denote the average velocity in the lateral direction and the height profile for the free surface, respectively. The positive slip length parameters b and β are related by orders of magnitude via b ∼ η 2 β, where the (small) parameter η, 0 < η ≪ 1, refers to the vertical to horizontal scale separation of the thin film. The high order of the lubrication equations (1.1) and (1.2a)-(1.2b) is a result of the contribution from surface tension at the free boundary, reflected by the linearized curvature term σ∂ xx h with parameter σ ≥ 0. A further contribution to the pressure is denoted by Π(h) and represents that of the intermolecular forces, namely long-range attractive van der Waals and short-range Born repulsive intermolecular forces. A commonly used expression for it is given by Π(h) = 1 h 3 -α h 4 with α > 0.

(1.

3)

The terms Re (∂ t (hu)+∂ x (hu 2 )), with Re ≥ 0 denoting the Reynolds number, and 4∂ x (ν(h)∂ x u) in (1.2a)-(1.2b) represent inertial and Trouton viscosity terms, respectively. We assume below linear dependence of the viscosity coefficient on the height, namely ν(h) := νh, ν > 0.

For these choice of the constitutive laws, we first investigate the existence of global weak solutions to (1.2a)-(1.2b) on the interval (0, 1), supplemented with the boundary conditions u = 0 at x = 0, 1 , (1.4)

∂ x h = 0 at x = 0, 1 , (1.5)
and initial data (h 0 , u 0 ) with positive first component h 0 . Observe that, in particular, the boundary conditions (1.4) for u guarantee the conservation of mass

1 0 h(x, t) dx = h c := 1 0 h 0 (x) dx , t ≥ 0, (1.6)
where h c is the average of the height profile. We next investigate the behaviour of these weak solutions as some of the parameters in the model goes to zero or to infinity. More precisely, let us point out that the strong-slip model (1.2a)-(1.2b) includes as limiting cases three further lubrication models. One is obtained from the strong-slip model in the limit β → ∞ and describes the dynamics of suspended free films, see e.g. Brenner and Gueyffier [START_REF] Brenner | On the bursting of viscous films[END_REF]:

Re (∂ t (hu) + ∂ x (hu 2 )) = 4∂ x (ν(h)∂ x u) + h∂ x (σ∂ xx h -Π(h)), (1.7a) ∂ t h = -∂ x (hu) .
(1.7b)

In the second one, we neglect the inertial terms which corresponds to a vanishing Reynolds number. For the third limiting case which is derived in Münch et al. [START_REF] Münch | Lubrication models with small to large slip lengths[END_REF] the slip length parameter β I is of order of magnitude lying in between those that lead to the weak-slip model (1.

∂ t h = -∂ x h 2 ∂ x (σ∂ xx h -Π(h)) . (1.8)
It can be obtained by rescaling time in (1.1) by b and letting b → ∞ or by rescaling time and the horizontal velocity by β in (1.2a)-(1.2b) and taking the limit β → 0. We consider here the second limit and change variables in (1.2a)-(1.2b) as follows:

x := x , t := tβ , h := h , and u

:= ū β , (1.9)
where x, t, h, and ū denote the old variables. We show that, as β → 0, the rescaled solutions to (1.2a)-(1.2b) converge to a solution to

u = h∂ x (σ∂ xx h -Π(h)), (1.10a) ∂ t h = -∂ x (hu) , (1.10b)
satisfying boundary conditions (1.4)-(1.5) that is, a solution to the intermediate-slip equation (1.8).

Our proofs follow a strategy that has also been employed in other systems arising in fluid dynamics and phase transition theories. In fact, Bresch et al. [START_REF] Bresch | On some compressible fluid models: Korteweg, lubrication and shallow water systems[END_REF] and Bresch [START_REF] Bresch | Shallow-water equations and related topics[END_REF] considered the following Korteweg system

∂ t (hu) + div (hu ⊗ u) = div (µ(h)(∇u + (∇u) T ) + ∇(λ(h)div u) (1.11a) + σh∇∆h -∇P (h) ∂ t h = -div (hu) .
(1.11b)

In the case σ = 0, µ(h) = µ = const, λ(h) = λ = const,
the system(1.11a)-(1.11b) gives the compressible Navier-Stokes equations. For the latter, considered on a bounded domain with P (h) having a globally Lipschitz continuous derivative, existence of a unique strong solution was shown by Solonnikov [START_REF] Solonnikov | The solvability of the initial-boundary value problem for the equations of motion of a viscous compressible fluid[END_REF] in dimensions d ≥ 2. Recently, Bresch et al. [START_REF] Bresch | On some compressible fluid models: Korteweg, lubrication and shallow water systems[END_REF] showed existence of a weak solution to (1.11a)-(1.11b) on bounded domain with σ ≥ 0, µ(h) = νh, λ(h) = 0 (1.12)

and similar requirements on P (h) for dimensions d ≥ 2. In the proof the authors use a so called BD entropy equality introduced by them originally in Bresch and Desjardins [START_REF] Bresch | Existence of global weak solutions for a 2D viscous shallow water model and convergence to the quasigeostrophic model[END_REF] for the case σ = 0. In the latter case (1.11a)-(1.11b) coincides with viscous shallow-water equations. In a recent review by Bresch [START_REF] Bresch | Shallow-water equations and related topics[END_REF] on those equations it was shown that the BD entropy equality holds for (1.11a)-(1.11b) in dimensions d ≥ 2 if a special relation between viscosity coefficients is satisfied, namely

λ(h) = 2(µ ′ (h)h -µ(h)) (1.13)
that is satisfied for (1.12). We also notice that the strong-slip lubrication equation (1.2a)-(1.2b) has a similar form to the one dimensional case of (1.11a)- (1.11b). Indeed, in one space dimension (1.11a)-(1.11b) reduces to

∂ t (hu) + ∂ x (hu 2 ) = ∂ x (ν(h)∂ x u) + σh∂ 3 x h -∂ x P (h) (1.14a) ∂ t h = -∂ x (hu) . (1.14b)
Recently, in Mellet and Vasseur [START_REF] Mellet | Existence and uniqueness of global strong solutions for one-dimensional compressible Navier-Stokes equations[END_REF] an analogue of the BD entropy identity was shown for (1.14a)-(1.14b) with any C 2 smooth viscosity function ν(h) and σ = 0. Furthermore, they proved existence of strong solutions to (1.14a)-(1.14b) on the whole real line R with σ = 0, ν(h) = νh k and P (h) = h γ , where k ≤ 1/2 and γ > 1.

The main difference in (1.2a)-(1.2b) compared to (1.14) is the special form of the intermolecular pressure (1.3) which is singular when the height h vanishes and complicates the analysis. Nevertheless, it also provides a very useful positive lower bound for h. Equations (1.2a)-(1.2b) also include an additional slip term which is easily handled for the existence theory but plays a crucial role in the limit β → 0 (leading to equation (1.8)).

In this paper we start in section 2 with recalling the energy identity associated to (1.2a)-(1.2b) and derive a version of a BD entropy identity. These two estimates provide the main ingredients of our existence results. In section 3 we follow the strategy of Bresch and Desjardins [START_REF] Bresch | On the construction of approximate solutions for the 2D viscous shallow water model and for compressible Navier-Stokes models[END_REF] by setting up a higher-order regularized equation for which analogous energy and entropy inequalities are satisfied. The a priori estimates allow us to pass to the limit in the regularization parameter to obtain global weak solutions to (1.2a)-(1.2b). In section 4 we establish existence of global weak solutions for the cases Re = 0, β = ∞ and σ = 0 by passing to the respective limits in (1.2a)-(1.2b). While the first two cases are straightforward, we need to refine the estimates from section 2 in order to enable us to deal with the limit σ → 0. Finally, in section 5 we establish that solutions of (1.2a)-(1.2b) converge, after the rescaling (1.9), in the limit β → 0 to a solution of the intermediate-slip model (1.8).

A priori estimates.

In this section we state two relations that are satisfied by classical solutions of (1.2a)-(1.2b). As it was stated in the introduction both have their counterparts for viscous shallow-water equation (1.11a)-(1.11b) and lubrication equations (1.1) and (1.8). For the latter they were initially derived by Bernis and Friedman [START_REF] Bernis | Higher order nonlinear degenerate parabolic equations[END_REF]. The energy equality for shallow-water equations is known already for several decades, whereas the entropy equality was suggested recently by Bresch and Desjardins [START_REF] Bresch | Existence of global weak solutions for a 2D viscous shallow water model and convergence to the quasigeostrophic model[END_REF]. For the strong-slip model (1.2a)-(1.2b) the derivation of the energy equality is again standard, see Kitavtsev and Wagner [START_REF] Kitavtsev | Coarsening dynamics of slipping droplets[END_REF]. As a new result we derive here the entropy equality for (1.2a)-(1.2b) following the approach of Bresch and Desjardins [START_REF] Bresch | Existence of global weak solutions for a 2D viscous shallow water model and convergence to the quasigeostrophic model[END_REF]. At the end of this section we use the energy and entropy equalities to derive a priori estimates in the case σ > 0 on classical solutions to (1.2a)-(1.2b) having positive h and satisfying (1.6), (1.4)- (1.5). For consistency we start with the energy equality.

d dt 1 0 Re h u 2 2 + U(h) + σ |∂ x h| 2 2 dx = -4 1 0 νh|∂ x u| 2 dx - 1 0 u 2 β dx, (2.1)
where the potential function U is the indefinite integral of Π defined by

U(h) = - 1 2 h 2 + α 3 h 3 , h > 0 .
(2.2)

Proof: One way to show (2.1) is to use the fact that

E(u, h) := 1 0 Re h u 2 2 + U(h) + σ |∂ x h| 2 2 dx (2.3)
is a Lyapunov functional for the system (1.2a)-(1.2b), see Kitavtsev and Wagner [START_REF] Kitavtsev | Coarsening dynamics of slipping droplets[END_REF]. Here we give another standard derivation of (2.1). Multiplying (1.2a) by u and integrating in x we obtain Re

1 0 ∂ t (hu)u + ∂ x (hu 2 )u dx -4 1 0 ∂ x (ν h∂ x u)u dx - 1 0 hu∂ x (σ∂ xx h -Π(h)) dx + 1 0 u 2 β dx = 0.
Using now several integrations by parts, (1.2b) and boundary conditions (1.4)-(1.5) we obtain 0 = Re

1 0 ∂ t h u 2 2 + ∂ t h u 2 2 -hu 2 ∂ x u dx + 4 1 0 ν h(∂ x u) 2 dx - 1 0 ∂ t h(σ∂ xx h -Π(h)) dx + 1 0 u 2 β dx = Re 1 0 ∂ t h u 2 2 -∂ x (hu) u 2 2 -hu∂ x u 2 2 dx + 4 1 0 ν h(∂ x u) 2 dx + 1 0 ∂ t σ |∂ x h| 2 2 + U(h) dx + 1 0 u 2 β dx = Re 1 0 ∂ t h u 2 2 dx + 4 1 0 ν h(∂ x u) 2 dx + 1 0 ∂ t σ |∂ x h| 2 2 + U(h) dx + 1 0 u 2 β dx.
From the last expression (2.1) follows.

Lemma 2.2. For classical solutions of the system (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5) the following equality holds

1 2 d dt 1 0 |∂ x h| 2 h dx = 1 0 ∂ x ∂ x h h ∂ x u h dx. (2.4)
Proof: The statement is verified again using several integrations by parts, equation (1.2a) and boundary conditions (1.4)-(1.5) as follows.

1 2

d dt 1 0 |∂ x h| 2 h dx = 1 0 ∂ x h∂ xt h h dx - 1 2 1 0 ∂ t h|∂ x h| 2 h 2 dx = - 1 0 ∂ x h h ∂ xx (hu) dx + 1 2 1 0 ∂ x (hu) ∂ x h h 2 dx = 1 0 ∂ x ∂ x h h ∂ x (hu) dx - 1 0 u∂ x h∂ x ∂ x h h dx = 1 0 ∂ x ∂ x h h ∂ x u h dx.
Let us define a so-called entropy function by ϕ(h) := 4ν log h for h > 0. Then we have the following lemma.

d dt 1 0 1 2 h(Re u + ∂ x ϕ(h)) 2 - 1 β ϕ(h) + Re σ |∂ x h 2 2 + U(h) dx = -Re 1 0 u 2 β dx -4σν 1 0 |∂ xx h| 2 dx -4ν 1 0 Π ′ (h)|∂ x h| 2 dx. (2.5)
Proof: Let us multiply equation (1.2a) by ∂ x ϕ(h), integrate with respect to x and use an integration by parts, (1.2b) and (1.5):

4νRe 1 0 (∂ t u + u∂ x u)∂ x h dx = -16ν 2 1 0 ∂ x ∂ x h h ∂ x u h dx -4ν 1 0 u∂ x h βh dx -4σν 1 0 |∂ xx h| 2 dx -4ν 1 0 Π ′ (h)|∂ x h| 2 dx. (2.6)
On one hand, a further integration by parts in the left-hand side of (2.6), equation (1.2b), and the energy equality (2.1) give 4νRe 1 0

(∂ t u + u∂ x u)∂ x h dx = 4νRe d dt 1 0 u∂ x h dx - 1 0 u∂ xt h dx + 1 0 u∂ x u∂ x h dx = 4νRe d dt 1 0 u∂ x h dx - 1 0 ∂ x u∂ x (uh) dx + 1 0 u∂ x u∂ x h dx = 4νRe d dt 1 0 u∂ x h dx - 1 0 h(∂ x u) 2 dx = Re d dt 1 0 4νu∂ x h + Re h u 2 2 + U(h) + σ |∂ x h| 2 2 dx + 1 0 u 2 β dx .
On the other hand, using Lemma 2.2 and the definition of ϕ, one can write the first term on the right-hand side of (2.6) as

-16ν 2 1 0 ∂ x ∂ x h h ∂ x uh dx = - 1 2 d dt 1 0 h |∂ x ϕ(h)| 2 dx,
while it follows from (1.2b) and (1.4) that the second term on the right-hand side of (2.6) can be transformed into -4ν

1 0 u∂ x h βh dx = -4ν 1 0 ∂ x (uh) βh dx + 4ν 1 0 ∂ x u β dx = 1 β d dt 1 0 ϕ(h) dx.
Substituting finally the last three identities into (2.6) one obtains

d dt 1 0 Re 2 2 hu 2 + 1 2 h |∂ x ϕ(h)| 2 + 4νRe u∂ x h - 1 β ϕ(h) + Re σ |∂ x h| 2 2 + U(h) dx = -Re 1 0 u 2 β dx -4σν 1 0 |∂ xx h| 2 dx -4ν 1 0 Π ′ (h)|∂ x h| 2 dx.
Using the definition of ϕ the last expression can be easily transformed into (2.5).

Lemma 2.4. For smooth functions (h, u) with a positive first component h, we have

1 4 1 0 h|∂ x ϕ(h)| 2 dx ≤ 1 2 1 0 h (Re u + ∂ x ϕ(h)) 2 dx + Re E(u, h) + Re 6α 2 , (2.7)
where E(u, h) is defined in (2.3).

Proof: Using the elementary inequality

(y + z) 2 ≥ y 2 2 -z 2 , the fact that U(h) ≥ -1/(6α 2 ) (2.8)
for all h > 0, and the definition (2.3) of E one obtains

1 0 h (Re u + ∂ x ϕ(h)) 2 dx ≥ 1 2 1 0 h|∂ x ϕ(h)| 2 dx -Re 2 1 0 hu 2 dx ≥ 1 2 1 0 h|∂ x ϕ(h)| 2 dx -2Re E(u, h) - 1 0 U(h) + σ |∂ x h| 2 2 dx ≥ 1 2 1 0 h|∂ x ϕ(h)| 2 dx -2Re E(u, h) - Re 3α 2 ,
from which the statement of the lemma follows.

We now deduce several bounds from Lemmas 2.1 and 2.3 for classical solutions (h, u) to the system (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5) and a positive first component.

To this end, we assume that the initial data h 0 and u 0 for (1.2a)-(1.2b) satisfy h 0 ∈ H 1 (0, 1), h 0 (x) > 0 for all x ∈ [0, 1] and u 0 ∈ L 2 (0, 1).

(2.9)

As already mentioned, a first consequence of (1.2b) and (1.4) is the conservation of mass (1.6). We next combine (1.6), the energy and entropy identities to obtain the following bounds when σ > 0.

Proposition 2.5 (a priori estimates). For fixed positive σ, Re, β, T and initial data satisfying (2.9) there exists a positive constant C 0 > 1 depending only on T , α, ν, Re, σ, E(u 0 , h 0 ), and S β (u 0 , h 0 ) defined in (2.13) below such that the following terms are bounded by

C 0 in respective norms √ h, ∂ x √ h, h -3/2 , ∂ x h, √ Re √ hu ∈ L ∞ (0, T ; L 2 (0, 1)), (2.10) ∂ x (h -3/2 ), ∂ xx h, √ h∂ x u, u √ β ∈ L 2 ((0, 1) × (0, T )), (2.11)
and

C -1 0 ≤ h(x, t) ≤ C 0 (2.12)
for all x ∈ (0, 1) and t ∈ (0, T ). The constant C 0 tends to ∞ as σ → 0.

Proof: Integrating the energy equality (2.1) with respect to time over (0, t), t ∈ (0, T ), and using the inequality

U(h) ≥ α 6h 3 - 2 3α 2 ensure that √ Re h(t)u(t) 2 2 + α 3 h -3/2 (t) 2 2 + √ σ∂ x h(t) 2 2 ≤ C 1 := 2E(u 0 , h 0 ) + 4 3α 2 , t 0 4ν h(s)∂ x u(s) 2 2 + u(s) √ β 2 2 ds ≤ C 1 .
Therefore, the norms of √ σ∂ x h, h -3/2 , and

√ Re √ hu in L ∞ (0, T ; L 2 (0, 1)
) are bounded by a constant that depends only on C 1 while the norms of √ h∂ x u and u/ √ β in L 2 ((0, 1)×(0, T )) are bounded by a constant that depends only on C 1 and ν.Using the bound on ∂ x h one obtains

|h(x, t) -h(y, t)| = y x ∂ x h(z, t) dz ≤ |x -y| 1/2 ∂ x h(t) 2 ≤ C 1 √ σ |x -y| 1/2
for all (x, y) ∈ (0, 1) × (0, 1) and t ∈ (0, T ). Integrating the above inequality with respect to y ∈ (0, 1) and using (1.6) readily give the upper bound in(2.12). To establish the lower bound for h in (2.12), we combine the L ∞ (0, T ; L 2 (0, 1))-estimates on h -3/2 and √ σ∂ x h just established to obtain a bound on the norm of 1/ √ h in L ∞ (0, T ; W 1,1 (0, 1)) since

1 0 ∂ x h -1/2 dx = 1 2 1 0 |∂ x h| h 3/2 dx ≤ 1 2 √ σ √ σ∂ x h 2 h -3/2 2 .
Owing to the continuous embedding of W 1,1 (0, 1) in L ∞ (0, 1), the positive lower bound in (2.12) readily follows.

Next, let us introduce the functional S β (u, h) defined by

S β (u, h) := 1 0 1 2 h(Re u + ∂ x ϕ(h)) 2 + 1 β (4νh -ϕ(h)) + Re σ |∂ x h| 2 2 + Re U(h) dx. (2.13)
It follows from the mass conservation (1.6) and the entropy equality (2.5) that

S β (u(t), h(t)) + t 0 1 0 Re β u 2 + 4σν|∂ xx h| 2 + 4νΠ ′ (h)|∂ x h| 2 dx ds = S β (u 0 , h 0 ) (2.14)
for any t ∈ (0, T ). Since zlog(z) ≥ 1 for all z > 0, it follows from the energy equality (2.1), (2.8), and Lemma 2.4 that

S β (u, h) ≥ 1 2 1 0 h(Re u + ∂ x ϕ(h)) 2 dx + 1 β 1 0 (4νh -ϕ(h)) dx - Re 6α 2 ≥ 1 2 1 0 h(Re u + ∂ x ϕ(h)) 2 dx - Re 6α 2 (2.15) ≥ 1 4 1 0 h|∂ x ϕ(h)| 2 dx -Re E(u, h) - Re 3α 2 ≥ 1 4 1 0 h|∂ x ϕ(h)| 2 dx -Re E(u 0 , h 0 ) - Re 3α 2 .
Combining this with the previous estimate on √ σ∂ x h, (2.14), and the fact that

Π ′ (h) ≥ 2α h 5 - 6 5 5 1 2α 4 implies that 1 4 1 0 h(t)|∂ x ϕ(h(t))| 2 dx + t 0 1 0 4σν|∂ xx h(s)| 2 + 8να|∂ x h(s)| 2 h 5 (s) dx ds ≤ C 2 + 6 5 5 2ν α 4 t 0 1 0 |∂ x h(s)| 2 dx ds ≤ C 2 + 6 5 5 2ν α 4 C 2 1 T σ with C 2 := S β (u 0 , h 0 ) + Re E(u 0 , h 0 ) + (Re/3α 2 ). Since h|∂ x ϕ(h)| 2 = 4|∂ x √ h| 2
, this completes the proof of (2.10)-(2.11).

Remark 2.6. As in Theorem 1 of Bertozzi et al. [START_REF] Bertozzi | Dewetting films: bifurcations and concentrations[END_REF] that shows existence of uniform lower bound for solutions of (1.8) many estimates in the proof of Proposition 2.5 heavily rely on the positivity of the parameter α in Π and become unbounded as α → 0. The estimates (2.10), (2.11), and (2.12) are no longer valid in that case. The same remark holds for the case σ = 0 since the estimates depend on the bound for √ σ∂ x h. On the contrary the proof holds without changes for the cases Re = 0 or β = ∞ except for the fact that then one looses estimates on Re √ hu and u/ √ β in (2.10)-(2.11), respectively.

3 Existence of weak solutions.

We first define a weak formulation of the problem (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5). Consider two functions h 0 and u 0 satisfying (2.9).

h∂ t ψ dxdt + 1 0 h 0 ψ(•, 0) dx = - ∞ 0 1 0 hu∂ x ψ dxdt, (3.1)
Re ∞ 0 1 0 hu∂ t φ dxdt + Re 1 0 h 0 u 0 φ(•, 0) dx + Re ∞ 0 1 0 hu 2 ∂ x φ dx dt -4ν ∞ 0 1 0 h∂ x u∂ x φ dx dt -σ ∞ 0 1 0 ∂ x h∂ xx hφ dx dt -σ ∞ 0 1 0 h∂ xx h∂ x φ dx dt + ∞ 0 1 0 Π 1 (h)∂ x φ dxdt - 1 β ∞ 0 1 0 uφ dxdt = 0 (3.2)
for all ψ ∈ C ∞ 0 ([0, 1] × [0, ∞)) and φ ∈ C ∞ 0 ((0, 1) × [0, ∞)),
where

Π 1 (h) := - ∞ h τ Π ′ (τ ) dτ. (3.3)
To show the existence of weak solutions to (1.2a)-(1.2b)-(1.4)-(1.5) we construct approximating systems similar to those suggested by Bresch and Desjardins [START_REF] Bresch | On the construction of approximate solutions for the 2D viscous shallow water model and for compressible Navier-Stokes models[END_REF] for the Korteweg and viscous shallow-water equations (1.11a)-(1.11b). Although the pressure term Π(h) does not need a regularization as in the case of Bresch and Desjardins [START_REF] Bresch | On the construction of approximate solutions for the 2D viscous shallow water model and for compressible Navier-Stokes models[END_REF], one still needs to regularize the function h sufficiently in order to control additional higher order terms arising in the entropy equality. The approximating systems we take are given by

Re (∂ t (h ε u ε ) + ∂ x (h ε u 2 ε )) = 4∂ x (ν h ε ∂ x u ε ) + h ε ∂ x (σ∂ xx h ε -Π(h ε )) - u ε β (3.4a) + εh ε ∂ 7 x h ε -ε 2 ∂ 4 x u ε , ∂ t h ε = -∂ x (h ε u ε ) , (3.4b)
where ε > 0 is a small parameter. Consider (3.4a)-(3.4b) with boundary conditions

u ε = ∂ xx u ε = ∂ x h ε = ∂ 3 x h ε = ∂ 5 x h ε = 0, (x, t) ∈ {0, 1} × (0, T ), (3.5)
and initial conditions

u ε (x, 0) = u ε,0 (x) and h ε (x, 0) = h ε,0 (x) > 0, x ∈ (0, 1), (3.6)
where u ε,0 and h ε,0 are smooth functions such that

u ε,0 → u 0 in L 2 (0, 1), h ε,0 → h 0 in H 1 (0,
E ε (u ε , h ε ) := 1 0 Re h ε u 2 ε 2 + U(h ε ) + σ |∂ x h ε | 2 2 + ε 2 |∂ 3 x h ε | 2 dx,
and the entropy

S ε (u ε , h ε ) := 1 0 1 2 h ε (Re u ε + ∂ x ϕ(h ε )) 2 - 1 β ϕ(h ε ) dx + Re 1 0 σ |∂ x h ε | 2 2 + U(h ε ) + ε 2 |∂ 3 x h ε | 2 dx,
the corresponding energy and entropy equalities read

d dt E ε (u ε , h ε ) = -4 1 0 νh ε |∂ x u ε | 2 dx - 1 0 u 2 ε β dx -ε 2 1 0 |∂ xx u ε | 2 dx (3.8)
and

d dt S ε (u ε , h ε ) = -Re 1 0 u 2 ε β dx -4σν 1 0 |∂ xx h ε | 2 dx -4ν 1 0 Π ′ (h ε )|∂ x h ε | 2 dx - 1 0 Re ε 2 |∂ xx u ε | 2 + 4νε|∂ 4 x h ε | 2 + 4νε 2 ∂ xx u ε ∂ 3 x log h ε dx. (3.9)
Given ε > 0 equation (3.4a) is parabolic with respect to u. Arguing as in Bresch and Desjardins [START_REF] Bresch | On the construction of approximate solutions for the 2D viscous shallow water model and for compressible Navier-Stokes models[END_REF], the initial-boundary value problem system

, T ).

There exists a positive constant C 3 ≥ 1 depending only on T , α, ν, Re, β, σ, h 0 , and u 0 such that, for all sufficiently small ε > 0, the following terms are bounded by C 3 in respective norms

h ε , ∂ x h ε , h -3/2 ε , ∂ x h ε , h ε u ε , √ ε∂ 3 x h ε ∈ L ∞ (0, T ; L 2 (0, 1)), (3.10)
∂ x h -3/2 ε , ∂ xx h ε , h ε ∂ x u ε , u ε , √ ε∂ 4 x h ε , ε∂ xx u ε ∈ L 2 ((0, 1) × (0, T )), (3.11
)

and 1 C 3 ≤ h ε (x, t) ≤ C 3 , (x, t) ∈ (0, 1) × (0, T ) . ( 3

.12)

Proof: Throughout the proof C denotes a positive constant depending on the same parameters as C 3 that may vary from line to line. Proceeding as in the proof of Proposition 2.5, the mass conservation and the energy equality (3.8) imply the estimate (3.12) and all bounds in (3.10)-(3.11) except those on

∂ x √ h ε , ∂ xx h ε , ∂ x h -3/2 ε and √ ε∂ 4 x h ε .
Using the previously obtained bounds (in particular the lower bound (3.12) on h ε ) and the interpolation inequalities

∂ xx h ε L ∞ (0,1) ≤ C ∂ 3 x h ε L 2 (0,1) and ∂ x h ε L 6 (0,1) ≤ C ∂ xx h ε 1/3 L 2 (0,1) ∂ x h ε 2/3 L 2 (0,1) ,
one can estimate the last term on the right-hand side of the entropy equality (3.9) as follows:

4νε 2 T 0 1 0 ∂ xx u ε • ∂ 3 x log h ε dx dt ≤ 4νε T 0 1 0 ε 2 |∂ xx u ε | 2 dx dt 1/2 sup [0,T ] ||∂ 3 x log h ε || L 2 (0,1) ≤ C ε sup [0,T ] ∂ 3 x h ε h ε -3 ∂ x h ε ∂ xx h ε h 2 ε + 2 |∂ x h ε | 3 h 3 ε L 2 (0,1) ≤ C ε sup [0,T ] ∂ 3 x h ε L 2 (0,1) + ∂ x h ε L 2 (0,1) ∂ xx h ε || L ∞ (0,1) + ∂ x h ε 3 L 6 (0,1) ≤ C ε sup [0,T ] ∂ 3 x h ε L 2 (0,1) + C ∂ 3 x h ε L 2 (0,1) + C ∂ xx h ε L 2 (0,1) + C √ ε sup [0,T ] √ ε ∂ 3 x h ε L 2 (0,1) ≤ C √ ε. (3.13)
Therefore, taking ε sufficiently small, the bounds in (3.10)-(3.11) still to be proved follow from the entropy inequality (3.9) exactly as in the proof of Proposition 2.5.

We now show that solutions to the system (3. We investigate the compactness properties of the sequence (h εn , u εn ) and to this end use the estimates derived in Proposition 3.2. Let T > 0. First, owing to (3.10) and (3.12), ( h εn ) and ( h εn u εn ) are bounded in L ∞ ((0, 1) × (0, T )) and L ∞ (0, T ; L 2 (0, 1)), respectively. Hence, by (3.4b), (∂ t h εn ) is bounded in L ∞ (0, T ; H -1 (0, 1)). Since (h εn ) is bounded in L ∞ (0, T ; H 1 (0, 1)) and L 2 (0, T ; H 2 (0, 1)) by (3.10) and (3.11), it follows from the compactness of the embedding of H 1 (0, 1) in C([0, 1]) and Corollary 4 in Simon [START_REF] Simon | Compact sets in the space L p (0, T ; B)[END_REF] that there is h ∈ C([0, 1] × [0, T ]) such that, after possibly extracting a subsequence,

h εn → h in L 2 (0, T ; W 1,p (0, 1)) ∩ C([0, 1] × [0, T ]) for p ∈ [1, ∞), (3.14) ∂ t h εn ⋆ ⇀ ∂ t h in L ∞ (0, T ; H -1 (0, 1)).
Combining (3.12) and (3.14) we additionally obtain that h is positive and

h -1 εn → h -1 in C([0, 1] × [0, T ]). (3.15)
We next turn to compactness properties of (u εn ). For that purpose, we write (3.4a) as

Re ∂ t (h εn u εn ) = -Re ∂ x (h εn u 2 εn ) + 4∂ x (ν h εn ∂ x u εn ) + h εn ∂ x (σ∂ xx h εn -Π(h εn )) - u εn β + ε n h εn ∂ 7 x h εn -ε 2 n ∂ 4 x u εn (3.16)
and claim that the right-hand side of (3.16) is bounded in L 2 (0, T ; H -3 (0, 1)). Indeed, by (3.10)-(3.12), (h εn u 2 εn ), ( h εn ), and ( h εn ∂ x u εn ) are in L ∞ (0, T ; L 1 (0, 1)), L ∞ ((0, 1) × (0, T )), and L 2 ((0, 1) × (0, T )), which imply that (h εn u 2 εn ) and (∂ x (h εn ∂ x u εn )) are bounded in L 2 (0, T ; H -1 (0, 1)). Next, for any ψ ∈ C ∞ 0 ((0, 1) × (0, T )) one obtains, using integration by parts and (3.10)-(3.12), that

T 0 1 0 ψh εn ∂ 3 x h εn dx dt = T 0 1 0 ∂ xx h εn [h εn ∂ x ψ + ψ∂ x h εn ] dx dt ≤ T 0 ∂ xx h εn L 2 (0,1) h εn L ∞ (0,1) ∂ x ψ L 2 (0,1) + ψ L ∞ (0,1) ∂ x h εn L 2 (0,1) dt ≤ C T 0 ψ 2 H 1 (0,1) dt 1/2 and T 0 1 0 ψ h εn ∂ x Π(h εn ) dx dt = T 0 1 0 ∂ x ψ Π 1 (h εn ) dx dt ≤ Π 1 (h εn ) L ∞ ((0,1)×(0,T )) T 0 ψ 2 H 1 (0,1) dt 1/2
, where Π 1 (h) is defined in (3.3). Similarly,

ε n T 0 1 0 ψh εn ∂ 7 x h εn dx dt = ε n T 0 1 0 ∂ 4 x h εn ψ ∂ 3 x h εn + 3∂ x ψ ∂ xx h εn + 3∂ xx ψ ∂ x h εn + h εn ∂ 3 x ψ dx dt ≤ ε n T 0 ∂ 4 x h εn L 2 (0,1) ψ L ∞ (0,1) ∂ 3 x h εn L 2 (0,1) + h εn L ∞ (0,1) ∂ 3 x ψ L 2 (0,1) + 3 ∂ x ψ L ∞ (0,1) ∂ xx h εn L 2 (0,1) + 3 ∂ xx ψ L ∞ (0,1) ∂ x h εn L 2 (0,1) dt ≤ C √ ε n T 0 ∂ 4 x h εn L 2 (0,1) ψ H 3 (0,1) dt ≤ C ψ L 2 (0,T ;H 3 (0,1)) .
Finally, (u εn ) and (ε n ∂ xx u εn ) are bounded in L 2 ((0, 1) × (0, T )) and L 2 (0, T ; H -2 (0, 1)), respectively, by (3.11). Collecting the above information completes the proof of the boundedness of the right-hand side of (3.16), whence

(∂ t (h εn u εn )) is bounded in L 2 (0, T ; H -3 (0, 1)). (3.17)
Now, owing to (3.11) and (3.12), we realize that (u εn ) is bounded in L 2 (0, T ; H 1 (0, 1)). This fact along with (3.10)-(3.12) allows us to conclude that (h εn u εn ) is bounded in L 2 (0, T ; H 1 (0, 1)). Combining this with (3.17) and Corollary 4 in Simon [START_REF] Simon | Compact sets in the space L p (0, T ; B)[END_REF] ensures that (h εn u εn ) is compact in L 2 ((0, 1) × (0, T )). Then, thanks to (3.15), there exists u ∈ L 2 (0, T ; H 1 (0, 1)) such that

∂ x u εn ⇀ ∂ x u in L 2 ((0, 1) × (0, T )) and u εn → u in L 2 ((0, 1) × (0, T )). ( 3
1 0 Re h u 2 2 + U(h) + σ |∂ x h| 2 2 dx + T 0 1 0 4νh|∂ x u| 2 + u 2 β dxdt ≤ 1 0 Re h 0 u 2 0 2 + U(h 0 ) + σ |∂ x h 0 | 2 2 dx (3.19)
and entropy inequality sup 2) with a fixed Re > 0 and β = β n by (h βn , u βn ). Though the estimate on (u βn /β n ) is useless in that case, one still recovers the estimate of (u βn ) in L 2 (0, T ; H 1 0 (0, 1)) as a consequence of (2.11), (2.12), and the Poincaré inequality. Arguing as in the proof of Theorem 3.3, we conclude that, after possibly extracting a subsequence, (h βn , u βn ) converges towards a weak solution to the strong-slip model (1.7a)-(1.7b) describing the dynamics of free suspended films.

t∈(0,T ] 1 0 1 2 h (Re u + ∂ x ϕ(h)) 2 - 1 β ϕ(h) + Re σ |∂ x h| 2 2 + U(h) dx + T 0 1 0 Re u 2 β + 4σν|∂ xx h| 2 + 4νΠ ′ (h)|∂ x h| 2 dxdt ≤ 1 0 1 2 h 0 (Re u 0 + ∂ x ϕ(h 0 )) 2 - 1 β ϕ(h 0 ) + Re σ |∂ x h 0 | 2 2 + U(h 0 ) dx. (3.20)
If we now consider a sequence of positive real numbers (Re n ), Re n → 0, and denote the corresponding weak solutions to (3.1)-(3.2) with Re = Re n and a fixed β ∈ (0, ∞] by (h Re n , u Re n ), we may again proceed as in the proof of Theorem 3.3 to show that, after possibly extracting a subsequence, (h Ren , u Ren ) converges towards a weak solution to the strong-slip model (1.2a)-(1.2b) without inertial terms. There is however an important difference as the bound on (∂ t ( h Ren u Ren )) in L 2 (0, T ; H -3 (0, 1)) is no longer available and we only obtain the weaker conclusion

u Ren ⇀ u in L 2 (0, T ; H 1 (0, 1)).
Still, owing to the strong convergence of (h Ren ), this allows to pass to the limit as n → ∞ in all the terms which depends linearly on u Re n , that is, all the terms involving u Re n except Re n h Ren u 2

Ren . But the latter converges to zero as Re n → 0 since (h Ren ) and (u Ren ) are bounded in L ∞ ((0, 1) × (0, T )) and L 2 (0, T ; H 1 0 (

4.2 Case σ = 0.

As already pointed out in Remark 2.6, for weak solutions (h, u) to (3.1)-(3.2), the estimate on ∂ x h in (2.11) and both bounds in (2.12) depend on σ and explode as σ → 0. In order to investigate the limiting behaviour of weak solutions to (3.1)-(3.2) as σ → 0, we refine the a priori bounds of Proposition 2.5 in the next proposition so as to avoid their dependence on σ, although paying for it with an exponential growth with respect to time.

Proposition 4.1. For fixed positive σ, Re, β, T and initial data satisfying (2.9), let (h, u) be a weak solution to (3.1)-(3.2). There exists a positive constant C 4 > 1 depending only on T , α, ν, Re, E(u 0 , h 0 ), and S β (u 0 , h 0 ) defined in (2.13) (but not on σ ∈ (0, 1)) such that the following terms are bounded by

C 4 in respective norms √ h, ∂ x √ h, h -3/2 , √ hu ∈ L ∞ (0, T ; L 2 (0, 1)), (4.1)
∂ x (h -3/2 ), √ σ∂ xx h, √ h∂ x u, u ∈ L 2 ((0, 1) × (0, T )), (4.2)
and

C -1 4 ≤ h(x, t) ≤ C 4 , (x, t) ∈ (0, 1) × (0, T ). ( 4

.3)

Proof: First, the L ∞ (0, T ; L 2 (0, 1))-estimates on √ h, √ hu and h -3/2 and the L 2 ((0, 1) × (0, T ))-estimates on √ h∂ x u and u readily follow from the mass conservation and the energy inequality (3.19) as before. Next, we actually estimate in a different way the term 4νΠ ′ (h)|∂ x h| 2 dxdt in the entropy inequality (3.20). More precisely, using the estimate

Π ′ (h) = 4α -3h h 5 ≥ χ (0,α) (h) α h 5 + χ (α,∞) (h) - 3 h 4 ≥ α h 5 -χ (α,∞) (h) α h 5 + 3 h 4 ≥ α h 5 -χ (α,∞) (h) 4 h 4 ,
where χ A denotes the characteristic function of a set A ⊂ R and recalling that ϕ(h) = 4ν log h, one obtains

Π ′ (h)|∂ x h| 2 ≥ α h 5 |∂ x h| 2 -χ (α,∞) (h) h|∂ x ϕ(h)| 2 4ν 2 α 3 ≥ α h 5 |∂ x h| 2 - h|∂ x ϕ(h)| 2 4ν 2 α 3 . (4.4)
Next, let us consider as in the proof of Proposition 2.5 the function S β (u, h) defined in (2.13). Thanks to the mass conservation (1.6), the entropy inequality (3.20) also reads

S β (u(t), h(t)) + t 0 1 0 Re u 2 β + 4σν|∂ xx h| 2 dx ds ≤ S β (u 0 , h 0 ) -4ν t 0 1 0 Π ′ (h)|∂ x h| 2 dx ds. (4.5)
Recalling that

S β (u, h) ≥ 1 2 1 0 h(Re u + ∂ x ϕ(h)) 2 dx - Re 6α 2 (4.6)
by (2.15), it follows from (4.4), (4.6), and Lemma 2.4 that

-4ν 1 0 Π ′ (h)|∂ x h| 2 dx ≤ -4αν 1 0 |∂ x h| 2 h 5 dx + 1 να 3 1 0 h|∂ x ϕ(h)| 2 dx ≤ -4αν 1 0 |∂ x h| 2 h 5 dx + 4 να 3 S β (u, h) + Re E(u, h) + Re 3α 2 ≤ -4αν 1 0 |∂ x h| 2 h 5 dx + 4 να 3 S β (u, h) + Re E(u 0 , h 0 ) + Re 3α 2 ,
where we used the energy inequality (3.19)

x √ h in L ∞ (0, T ; L 2 (0, 1)) and ∂ x (h -3/2 ),
√ σ∂ xx h in L 2 ((0, 1) × (0, T )) and completes the proof of (4.1)-(4.2). To check (4.3), we first notice that the bound on √ h in L ∞ (0, T ; H 1 (0, 1)) from (4.1) and the embedding of H 1 (0, 1) in L ∞ (0, 1) guarantee the upper bound in (4.3). Next, integrating the following identity

1 h(x) = 1 h(y) - x y ∂ z h(z) h 2 (z) dz , x ∈ (0, 1) ,
with respect to y and using the Cauchy-Schwarz inequality give

1 h(x) ≤ 1 0 dy h(y) + 1 0 |∂ z h(z)| 2 h(z) dz 1/2 1 0 dz h 3 (z) 1/2 ≤ h -3/2 2/3 L 2 (0,1) + 4 ∂ x √ h L 2 (0,1) h -3/2 L 2 (0,1) .
The lower bound in (4. Proof: Let us take a sequence {σ n } → 0 and, for each n ≥ 1, denote the corresponding weak solution to (3.1)-(3.2) with σ = σ n by (h σn , u σn ). Owing to Proposition 4.1 we may proceed as in the proof of Theorem 3.3 to show the existence of functions h and u satisfying

h σn → h and h -1 σn → h -1 in C([0, 1] × [0, T ]), ∂ t h σn * ⇀ ∂ t h in L ∞ (0, T ; H -1 (0, 1)) and ∂ x u σn ⇀ ∂ x u in L 2 ((0, 1) × (0, T )) and u σn → u in L 2 ((0, 1) × (0, T )). In addition, σ n ∂ xx h σn → 0 in L 2 ((0, 1) × (0, T ))
by (4.2), from which we deduce that

σ n ∂ x h σn ∂ xx h σn ⇀ 0 and σ n h σn ∂ xx h σn ⇀ 0 in L 1 ((0, 1) × (0, T ))
with the help of (4.1). This allows us to pass to the limit as n → ∞ in all integrals in (3.1)-(3.2) and conclude that (h, u) satisfies (3.1)-(3.2) with σ = 0.

5 Convergence to solutions of intermediate-slip equation.

In the introduction we mentioned that performing the change of variables (1.9) and passing formally to the limit as β → 0 in ( Then there exist a positive function h and a subsequence of (h n , u n ) (not relabeled) such that, for any T > 0,

h n → h in L 2 (0, T ; H 1 (0, 1)) ∩ C([0, 1] × [0, T ]), u n ⇀ u := h∂ x (σ∂ xx h -Π(h)) in L 2 ((0, 1) × (0, T )),
and h is a global weak solution to (1.8) considered with initial condition h 0 satisfying (2.9) and boundary conditions (1.4)-(1.5).

Proof: Owing to (3.1)-(3.2), for each n ≥ 1, the functions (h n , u n ) satisfy the following scaled weak formulation

∞ 0 1 0 h n ∂ t ψ dx dt + 1 0 h 0 ψ(., 0) dx = - ∞ 0 1 0 h n u n ∂ x ψ dxdt, (5.2a)
β 2 n Re ∞ 0 1 0 h n u n ∂ t φ dx dt + Re 1 0 h 0 u 0 φ(., 0) dx + Re ∞ 0 1 0 h n u 2 n ∂ x φ dx dt -4β n ν ∞ 0 1 0 h n ∂ x u n ∂ x φ dx dt -σ ∞ 0 1 0 ∂ x h n ∂ xx h n φ dx dt - ∞ 0 1 0 (σ h n ∂ xx h n -Π 1 (h n )) ∂ x φ dx dt - ∞ 0 1 0 u n φ dx dt = 0 (5.2b) for all ψ ∈ C ∞ 0 ([0, 1] × [0, ∞)) and φ ∈ C ∞ 0 ((0, 1) × [0, ∞)).
In addition, for any T > 0, (h n , u n ) satisfies the scaled energy inequality

sup t∈[0,T ] 1 0 β 2 n Re h n u 2 n 2 + U(h n ) + σ |∂ x h n | 2 2 dx + T 0 1 0 4νβ n h n |∂ x u n | 2 + u 2 n dx dt ≤ 1 0 β 2 n Re h 0 u 2 0 2 + U(h 0 ) + σ |∂ x h 0 | 2 2 dx
and the scaled entropy inequality

sup t∈[0,T ] 1 0 β n 2 h n (β n Re u n + ∂ x ϕ(h n )) 2 + β n Re σ |∂ x h n | 2 2 + U(h n ) -ϕ(h) dx + T 0 1 0 β n Re u 2 n + 4σν|∂ xx h n | 2 + 4νΠ ′ (h n )|∂ x h n | 2 dx dt ≤ β n 1 0 1 2 h 0 (β n Re u 0 + ∂ x ϕ(h 0 )) 2 + Re σ |∂ x h 0 | 2 2 + U(h 0 ) dx - 1 0 ϕ(h 0 ) dx.
Using the fact that σ > 0 and proceeding as in the proof of Proposition 2.5, it follows from the mass conservation and the scaled energy inequality that ( h n ), (h -3/2 n ), (∂ x h n ) are bounded in L ∞ (0, T ; L 2 (0, 1)), (

(u n ), ( β n h n ∂ x u n ) are bounded in L 2 ((0, 1) × (0, T )), (5.4) and there exists C 5 > 1 such that C -1 5 ≤ h n (x, t) ≤ C 5 , (x, t) ∈ (0, 1) × (0, T ). We next use the scaled entropy inequality as at the end of the proof of Proposition 2.5 to establish that (∂ x (h -3/2 n )), (∂ xx h n ), are bounded in L 2 ((0, 1) × (0, T )).

(5.5)

Proceeding as in the proof of Theorems 3.3 and 4.2 one obtains that there exist functions h and u such that h n → h in L 2 (0, T ; H 1 (0, 1)) ∩ C([0, 1] × [0, T ]), h -1 n → h -1 in C([0, 1] × [0, T ]), u n ⇀ u in L 2 ((0, 1) × (0, T )), β n h n ∂ x u n → 0 in L 2 ((0, 1) × (0, T )).

It is then straightforward to pass to the limit as β n → 0 in (5.2a)-(5.2b) to obtain ) is then uniformly parabolic and it is likely that classical parabolic regularity results ensure that h is actually a classical solution to (1.8). As it is shown in Bertozzi et al. [START_REF] Bertozzi | Dewetting films: bifurcations and concentrations[END_REF] that there is a unique positive classical solution to (1.8), this implies the convergence of the whole sequence (h n ) towards it.

6 Discussion.

We would like to note some open questions that remain to be solved for the strong-slip system (1.2a)-(1.2b) as well as the corresponding shallow-water equations (1.11a)-(1.11b). First, it would be of interest to figure out whether the weak solutions constructed in the previous sections are more regular. Partial answers have already been given for particular cases of (1.11a)-(1.11b) with σ = 0, see Solonnikov [START_REF] Solonnikov | The solvability of the initial-boundary value problem for the equations of motion of a viscous compressible fluid[END_REF] and Mellet and Vasseur [START_REF] Mellet | Existence and uniqueness of global strong solutions for one-dimensional compressible Navier-Stokes equations[END_REF] where the existence of strong solutions is established. Furthermore, in Vaynblat et al. [START_REF] Vaynblat | Rupture of thin viscous films by van der Waals forces: Evolution and self-similarity[END_REF] the two-dimensional analogue of (1.7a)-(1.7b) describing dewetting of three-dimensional free suspended films is derived and reads ∂ t (hu) + div (hu ⊗ u) = div (h(∇u + (∇u) T )) + 2∇(hdiv u) + σh∇∆h -∇P (h) (6.1a) ∂ t h =div (hu) , (6.1b) which corresponds to the choice µ(h) = h and λ(h) = 2h in (1.11a)-(1.11b). As mentioned in the introduction the BD entropy equality has been so far only established for (1.11a)-(1.11b) if the relation (1.13) for the viscosities is satisfied. Unfortunately this relation is not fulfilled for (6.1a)-(6.1b). Analogously a two-dimensional version of (1.2a)-(1.2b) can be constructed but will present the same difficulty, so that the analysis presented herein is not likely to extend in a straightforward way to the two-dimensional model.

Lemma 2 . 1 (

 21 Energy equality). For classical solutions of the system (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5) the following equality holds

Lemma 2 . 3 (

 23 Entropy equality). For classical solutions of the system (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5) the following equality holds.

Definition 3 . 1 .

 31 A pair (h, u) is a global weak solution to (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5) and initial conditions (h 0 , u 0 ) if h and u enjoy the regularity properties stated in (2.10)-(2.11)-(2.12) and the following holds

Proposition 3 . 2 .

 32 (3.4a)-(3.4b) with (3.5)-(3.6) has a unique classical solution at least locally in time. The next proposition establishes a uniform lower bound for h ε and thus guarantees the global in time solvability for (3.4a)-(3.4b). For fixed positive ε, σ, Re, β and T > 0 let (h ε , u ε ) be the solution of the initial-boundary value problem system (3.4a)-(3.4b) with (3.5)-(3.6) in (0, 1) × (0

Theorem 3 . 3 .

 33 4a)-(3.4b) with boundary and initial conditions (3.5)-(3.6) converge to a solution to (3.1)-(3.2) as ε → 0. For any positive σ, Re, β and initial data (h 0 , u 0 ) satisfying (2.9), there exists a global weak solution (h, u) to the system (1.2a)-(1.2b) with boundary conditions (1.4)-(1.5) and initial conditions (h 0 , u 0 ) in the sense of Definition 3.1. Proof: Take a sequence {ε n } n≥1 → 0 and, for each n ≥ 1, denote the corresponding solution to the approximate system (3.4a)-(3.4b)-(3.5)-(3.6) with ε = ε n by (h εn , u εn ).

. 18 )Remark 3 . 4 (

 1834 The convergences(3.14),(3.15), and(3.18) then allow to pass to the limit as n → ∞ in the weak formulation of the approximating systems (3.4a)-(3.4b)-(3.5)-(3.6) in order to get that (h, u) satisfies (3.1)-(3.2). Energy and entropy inequalities). It follows from the proof of Theorem 3.3 together with the energy and entropy identities (3.8), (3.9), the property (3.7) and the estimate (3.13) that the weak solutions to (3.1)-(3.2) obtained above satisfy the corresponding energy inequality sup t∈(0,T ]

  Note that the right-hand side of (3.19)-(3.20) does not depend on T . Consequently, the statements of Lemma 2.4 and Proposition 2.5 hold also for the weak solutions to (3.1)-(3.2) constructed in Theorem 3.3. 4 Limiting cases.

4. 1

 1 Cases Re = 0 and β = ∞. By Remark 3.4, the statement of Proposition 2.5 is true for the weak solutions to (3.1)-(3.2) provided by Theorem 3.3. We may then investigate the behaviour of these solutions as either Re → 0 or β → ∞. Let us first consider a sequence of positive real numbers (β n ), β n → ∞, and denote the corresponding solutions to (3.1)-(3.

3 )Theorem 4 . 2 .

 342 then follows from the above inequality by (4.1). Using the a priori bounds from Proposition 4.1 one can show existence of weak solutions to (3.1)-(3.2) in the case σ = 0. Let σ = 0. For any positive Re, β and initial data (h 0 , u 0 ) satisfying (2.9) there exists a global weak solution (h, u) to (1.2a)-(1.2b) satisfying (3.1)-(3.2).

∂ 1 0Π 1 1 0∂ 1 0h 2 Remark 5 . 2 .

 1111252 xx h∂ x (hφ) dxdt + ∞ 0 (h)∂ x φ dxdt = ∞ 0 uφ dx dt. (5.6b)Using an approximation argument, we may actually takeφ = h∂ x ψ in (5.6b) with ψ ∈ C ∞ 0 ((0, 1) × [0, ∞)) and obtain xx h∂ x (h 2 ∂ x ψ) dx dt + ∞ 0 ∂ x (Π(h))ψ dx dt,20which coincides with the weak formulation for (1.8) introduced in Bernis and Friedman[START_REF] Bernis | Higher order nonlinear degenerate parabolic equations[END_REF] and Bertozzi and Pugh[START_REF] Bertozzi | Long-wave instabilities and saturation in thin film equations[END_REF]. Owing to the positivity of the weak solution h to (1.8) constructed in Theorem 5.1, equation(1.8

  [START_REF] Oron | Long-scale evolution of thin liquid films[END_REF], and εh ε,0 → 0 in H 3 (0, 1) as ε → 0. (3.7)

	Proceeding as in Lemmas 2.1 and 2.3 one derives formulas analogous to (2.1) and (2.5)
	for the system (3.4a)-(3.4b) with boundary and initial conditions (3.5)-(3.6). Namely,
	introducing the energy

  For fixed positive Re, σ, let (β n ) be a sequence of positive real numbers, β n → 0, and, for n ≥ 1, denote a global weak solution to (3.1)-(3.2) with initial data (h 0 , u 0 ) satisfying (2.9) and β = β n by ( hn , ūn ). Using (1.9), we define

			1.2a)-(1.2b) gives(1.10a)-(1.10b) and consequently that
	solutions of the strong-slip equation are expected to converge to that of (1.8). The next
	theorem makes this formal procedure rigorous. More precisely, we show that scaled weak
	solutions to (1.2a)-(1.2b) converge to that of (1.8) when σ > 0. In this section, we denote
	the unscaled variables in (1.2a)-(1.2b) with bars.
	Theorem 5.1. h n (x, t) := hn x,	t β	and u n (x, t) :=	1 β	ūn x,	t β	, (x, t) ∈ (0, 1) × (0, ∞). (5.1)

Acknowledgements

GK acknowledges the support from the Weierstrass Institute and the postdoctoral scholarship at the Max-Planck-Institute for Mathematics in the Natural Sciences, Leipzig. PhL and BN were partially supported by the CNRS/RSE project JP090230 and the EPSRC Science and Innovation award to the Oxford Centre for Nonlinear PDE (EP/E035027/1).