Boundary layers for Navier-Stokes equations with slip boundary conditions

Similar documents
On convergence criteria for incompressible Navier-Stokes equations with Navier boundary conditions and physical slip rates

Remarks on the inviscid limit for the compressible flows

Sublayer of Prandtl boundary layers

Note on the analysis of Orr-Sommerfeld equations and application to boundary layer stability

VORTICITY LAYERS OF THE 2D NAVIER-STOKES EQUATIONS WITH A SLIP TYPE BOUNDARY CONDITION

On nonlinear instability of Prandtl s boundary layers: the case of Rayleigh s stable shear flows

VANISHING VISCOSITY LIMIT FOR THE 3D NONHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATION WITH SPECIAL SLIP BOUNDARY CONDITION

On the well-posedness of the Prandtl boundary layer equation

LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS

Global well-posedness of the primitive equations of oceanic and atmospheric dynamics

REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID

OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES

arxiv: v1 [math.ap] 27 Nov 2014

Serrin Type Criterion for the Three-Dimensional Viscous Compressible Flows

Recent progress on the vanishing viscosity limit in the presence of a boundary

Séminaire Laurent Schwartz

Stability of Shear Flow

Formulation of the problem

Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems

ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN

VANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE

On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions

Separation for the stationary Prandtl equation

Homogenization of the compressible Navier Stokes equations in domains with very tiny holes

Weak-strong uniqueness for the compressible Navier Stokes equations with a hard-sphere pressure law

Viscous capillary fluids in fast rotation

Homogenization of stationary Navier-Stokes equations in domains with tiny holes

ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen

Decay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients

OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

On the local well-posedness of compressible viscous flows with bounded density

IMA Preprint Series # 2143

PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems

EXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON

Existence of minimizers for the pure displacement problem in nonlinear elasticity

ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS

THE STOKES SYSTEM R.E. SHOWALTER

Global well-posedness and decay for the viscous surface wave problem without surface tension

RESEARCH STATEMENT GUNG-MIN GIE

Existence and Continuation for Euler Equations

The incompressible Navier-Stokes equations in vacuum

Zero Viscosity Limit for Analytic Solutions of the Primitive Equations

arxiv: v1 [math.ap] 21 Dec 2016

On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations

A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem

RETRACTED. INSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC

The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method

A new regularity criterion for weak solutions to the Navier-Stokes equations

TRAVELING WAVES IN 2D REACTIVE BOUSSINESQ SYSTEMS WITH NO-SLIP BOUNDARY CONDITIONS

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction

The 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University

arxiv: v1 [math.ap] 22 Jul 2016

Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator

Compressible hydrodynamic flow of liquid crystals in 1-D

ISABELLE GALLAGHER AND MARIUS PAICU

Instability of Finite Difference Schemes for Hyperbolic Conservation Laws

α 2 )(k x ũ(t, η) + k z W (t, η)

Measure-valued - strong uniqueness for hyperbolic systems

Singular Limits of the Klein-Gordon Equation

The Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany

CONNECTIONS BETWEEN A CONJECTURE OF SCHIFFER S AND INCOMPRESSIBLE FLUID MECHANICS

Topics in Fluid Dynamics: Classical physics and recent mathematics

Stability of Small Viscosity Noncharacteristic Boundary Layers. Guy Métivier MAB Université de Bordeaux I Talence Cedex, France

Decay in Time of Incompressible Flows

Low Froude Number Limit of the Rotating Shallow Water and Euler Equations

Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System

Mathematical analysis of the stationary Navier-Stokes equations

On the ill-posedness of the Prandtl equation

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

STOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN

Navier-Stokes equations in thin domains with Navier friction boundary conditions

Relaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni

THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS

Issues for a mathematical definition of LES

Relative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system

On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1

INSTITUTE of MATHEMATICS. ACADEMY of SCIENCES of the CZECH REPUBLIC

A non-newtonian fluid with Navier boundary conditions

Low Mach number limit of the full Navier-Stokes equations

Zdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)

On the zero viscosity limit in incompressible fluids

Mixed exterior Laplace s problem

A regularity criterion for the weak solutions to the Navier-Stokes-Fourier system

Euler Equations: local existence

arxiv: v1 [math.ap] 28 Apr 2009

Presenter: Noriyoshi Fukaya

Chapter 5. The Differential Forms of the Fundamental Laws

Quantum Hydrodynamic Systems and applications to superfluidity at finite temperatures

Zdzislaw Brzeźniak. Department of Mathematics University of York. joint works with Mauro Mariani (Roma 1) and Gaurav Dhariwal (York)

Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases

Splash singularity for a free-boundary incompressible viscoelastic fluid model

Mathematical analysis of low Mach number limit of some two phase flow models

Viscous non-linear theory of Richtmyer-Meshkov Instability. Abstract

Existence of Weak Solutions to a Class of Non-Newtonian Flows

STEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel

Nonlinear Evolution Equations 1

Good reasons to study Euler equation.

Transcription:

Boundary layers for Navier-Stokes equations with slip boundary conditions Matthew Paddick Laboratoire Jacques-Louis Lions Université Pierre et Marie Curie (Paris 6) Séminaire EDP, Université Paris-Est Créteil 16 th April 2015 1 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

1 Stability and instability in the 2D incompressible model 2 2 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

1 Stability and instability in the 2D incompressible model 2 3 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

2D incompressible Navier-Stokes on the half-space (x, y) Ω = R ]0, + [: t u ε + u ε u ε ε u ε + p ε = 0 div u ε = 0 (1) u ε t=0 = u0 ε, modelling the motion of a viscous, homogeneous fluid. Usual boundary conditions (BC) for macroscopic flows: impermeability of the boundary u 2 y=0 = 0, and the fluid does not slip along the boundary - Dirichlet BC u y=0 = 0. 4 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Nondimensional viscosity ε = 1/Re, Re the Reynolds number. Inviscid limit problem: behaviour of the solutions u ε relative to v, solution of the incompressible Euler equation t v + v v + q = 0 div v = 0 (2) v t=0 = v 0 when ε goes to 0 high-reynolds approximation. The Euler equation is of order one, so only one BC is required: v 2 y=0 = 0. 5 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

As Navier-Stokes solutions can converge to Euler solutions in the whole space, a Prandtl boundary layer expansion 1 is proposed: ( ) u ε y (t, x, y) v(t, x, y) + V t, x, (3) ε 0 ε Smithsonian National Air and Space Museum 1 Prandtl, Verhand. III Intern. Math. Kongresses Heidelberg, 1904 6 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

As Navier-Stokes solutions can converge to Euler solutions in the whole space, a Prandtl boundary layer expansion 1 is proposed: ( ) u ε y (t, x, y) v(t, x, y) + V t, x, (3) ε 0 ε Smithsonian National Air and Space Museum Boundary layer stability in a normed function space (E, ) if sup u ε (t) v(t) ε 0 0. t [0,T ] 1 Prandtl, Verhand. III Intern. Math. Kongresses Heidelberg, 1904 6 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Issues with the inviscid limit problem with Dirichlet BCs: if u = (u 1, u 2 )(t, x, y) = (ũ, εṽ)(t, x, y = ε 1/2 y), the boundary layer is governed by the ill-posed 2 Prandtl equation: t ũ + ũ x ũ + ṽ y ũ y y ũ = ( tv + v v) 1 y=0 x ũ + y ṽ = 0 (ũ, ṽ) y=0 = 0 lim y + ũ = v 1 y=0 ; 2 Gérard-Varet & Dormy, J. Amer. Math. Soc., 2010 3 Grenier, Comm. Pure Appl. Math., 2000 4 Kato, Seminar on nonlinear PDEs, 1984 7 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Issues with the inviscid limit problem with Dirichlet BCs: if u = (u 1, u 2 )(t, x, y) = (ũ, εṽ)(t, x, y = ε 1/2 y), the boundary layer is governed by the ill-posed 2 Prandtl equation: t ũ + ũ x ũ + ṽ y ũ y y ũ = ( tv + v v) 1 y=0 x ũ + y ṽ = 0 (ũ, ṽ) y=0 = 0 lim y + ũ = v 1 y=0 ; even if Prandtl is well-posed, ansatz (3) can be invalid due to creation of vorticity inside the boundary layer 3, and control of the vorticity (turbulence) near the boundary is crucial 4. 2 Gérard-Varet & Dormy, J. Amer. Math. Soc., 2010 3 Grenier, Comm. Pure Appl. Math., 2000 4 Kato, Seminar on nonlinear PDEs, 1984 7 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Our boundary condition: the Navier, or slip, boundary condition 5, a mixed BC that allows the fluid to slip along the boundary at a speed proportional to the normal stress: u 2 y=0 = 0 and y u ε 1 y=0 = 2a ε u ε 1 y=0 (4) The slip phenomenon is observed at microscopic levels, for instance in capillary blood vessels. 5 Navier, Mém. Acad. Roy. Sci. Inst. France, 1823 8 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Let a ε = a/ε β for a, β R. β > 0: approximation of Dirichlet BCs as ε 0, transition between the stable β = 0 case and the unstable vorticity in the Dirichlet (a ε = ) case. With a > 0 and β = 0: stability in L 2, using energy estimates 6. 6 Iftimie & Planas, Nonlinearity, 2006 9 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Let a ε = a/ε β for a, β R. β > 0: approximation of Dirichlet BCs as ε 0, transition between the stable β = 0 case and the unstable vorticity in the Dirichlet (a ε = ) case. With a > 0 and β = 0: stability in L 2, using energy estimates 6. We extend this result: Theorem 1.1 (Diff. Integral Eqns. 2014) Let u ε 0 be a converging sequence in L2 (Ω), whose limit v 0 is in H s (Ω) for some s > 2. Then, for every time T > 0, the sequence of global weak solutions to the Navier-Stokes equation u ε with initial condition u ε 0 and Navier boundary condition as above, converges in L ([0, T ], L 2 (Ω)) towards the unique strong solution v of the Euler equation with initial condition v 0 if either a > 0 and β < 1, or a < 0 and β 1/2. 6 Iftimie & Planas, Nonlinearity, 2006 9 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Further extensions. The theorem can be extended to higher dimensions by starting with the Leray energy inequality. 7 Sueur, J. Math. Fluid Mech., 2014 10 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Further extensions. The theorem can be extended to higher dimensions by starting with the Leray energy inequality. The first point of the theorem (with a 0) has been extended to the isentropic compressible model by Sueur 7. 7 Sueur, J. Math. Fluid Mech., 2014 10 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

1 Stability and instability in the 2D incompressible model 2 11 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Set a > 0 and β = 1/2: (4) is written y u 1 y=0 = 2a ε u 1 y=0. (5) 8 Grenier, Comm. Pure Appl. Math., 2000 12 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Set a > 0 and β = 1/2: (4) is written y u 1 y=0 = 2a ε u 1 y=0. (5) Shear flow: vector field of the form u sh (y) = (u sh (y), 0); stationary solution to the incompressible Euler equation (with constant pressure). It is linearly unstable if there exist λ C with R(λ) > 0, k R and ψ H 1 0 (R+ ) such that u(t, x, y) = e λt (e ikx ψ(y)) solves the Euler equation linearised around u sh : t u + u u sh + u sh u + p = 0 div u = 0 u 2 y=0 = 0, For the Euler equation, linear instability implies nonlinear instability 8. (6) 8 Grenier, Comm. Pure Appl. Math., 2000 12 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Theorem 1.2 (Diff. Integral Eqns. 2014) Let Y = y/ ε. We consider a family of boundary-layer scale initial conditions written as u sh (Y ), with u sh a linearly unstable shear profile for the Euler equation. We generate a reference solution to the Navier-Stokes equation as follows: t u sh (t, Y ) YY 2 u sh(t, Y ) = 0 u sh (0, Y ) = u sh (Y ) Y u sh (t, 0) = 2au sh (t, 0) 13 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Theorem 1.2 (Diff. Integral Eqns. 2014) Let Y = y/ ε. We consider a family of boundary-layer scale initial conditions written as u sh (Y ), with u sh a linearly unstable shear profile for the Euler equation. We generate a reference solution to the Navier-Stokes equation as follows: t u sh (t, Y ) YY 2 u sh(t, Y ) = 0 u sh (0, Y ) = u sh (Y ) Y u sh (t, 0) = 2au sh (t, 0) Then, for any n N, there exist δ 0 > 0 and ε 0 > 0 such that, for every 0 < ε < ε 0, there exist initial data u0 ε such that u ε 0(x, y) u sh (Y ) H s (Ω) Cεn for some s > 0, which generates a solution to the Navier-Stokes equation with boundary condition (5) that satisfies, at a time T ε n ε ln(ε), u ε (T ε, x, y) u sh (T ε, Y ) L (Ω) δ 0. 13 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Main ideas of the proof. After rescaling the equation (t, x, y) = ε 1/2 (t, x, y), construct an approximate solution to the Navier-Stokes equation written as: u ap (t, x, y) = u sh ( εt, y) + N ε jn Uj i (t, x, y) + ε jn+1/4 Uj b (t, x, y/ε 1/4 ). j=1 (7) Note the ε 1/4 amplitude of the boundary layer terms, hinted by the asymptotic expansion for fixed Navier BCs 9. It will allow to bound u ap independently of ε. Each Uj i and Uj b are built using another asymptotic expansion in powers of ε 1/8, e.g. U1 i = 8n 1 m=0 εm/8 um. i 9 Iftimie & Sueur, Arch. Ration. Mech. Anal., 2011 14 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The first term u0 i is chosen as a wavepacket of unstable solutions of (6): u0(t, i x, y) = ϕ 0 (k)e λ(k)t+ikx v 0 (k, y) dk (8) R with ϕ 0 C0 (R) localising the eigenmodes around the most unstable wavenumber. 15 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The first term u0 i is chosen as a wavepacket of unstable solutions of (6): u0(t, i x, y) = ϕ 0 (k)e λ(k)t+ikx v 0 (k, y) dk (8) R with ϕ 0 C0 (R) localising the eigenmodes around the most unstable wavenumber. Denote σ 0 = max{rλ(k) k R} the most unstable eigenvalue, and let g ε (t) = ε n (1 + t) 1/2 e σ0t. Time of study: t T0 ε τ such that g ε (T0 ε ) = 1, and τ > 0 independent of ε. 15 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The first term of u ap carries the growth we are after: u i 0 (t) L 2 (Ω A (t)) C 1(1 + t) 1/4 g ε (t), for Ω A (t) Ω of measure 1 + t. 16 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The first term of u ap carries the growth we are after: u i 0 (t) L 2 (Ω A (t)) C 1(1 + t) 1/4 g ε (t), for Ω A (t) Ω of measure 1 + t. By induction, we then show that ε jn (Uj i (t) + ε 1/4 Uj b (t)) C j (1 + t) 1/4 (g ε (t)) j. L 2 (Ω) This is done by noticing that the interior terms Uj i solve (6) with a source term F j satisfying F j (t) Hs (Ω) C j e γt (1 + t) α with γ > σ 0 and α > 0. A resolvent estimate for the linearised equation then yields that U i j also satisfies this inequality. 16 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Now we examine w = u ε u ap. Choosing N so that Nσ 0 > u ap L ([0,T0 ε τ] Ω) + 1/2, a quick energy estimate on the equation satisfied by w yields w(t) L 2 (Ω) C N+1(1 + t) 1/4 (g ε (t)) N+1. 17 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Now we examine w = u ε u ap. Choosing N so that Nσ 0 > u ap L ([0,T0 ε τ] Ω) + 1/2, a quick energy estimate on the equation satisfied by w yields w(t) L 2 (Ω) C N+1(1 + t) 1/4 (g ε (t)) N+1. As a result, u ε (t) u sh ( εt) N+1 L (1 + t) 1/4 C 2 (Ω A (t)) 1g ε (t) C j (g ε (t)) j j=2 17 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Now we examine w = u ε u ap. Choosing N so that Nσ 0 > u ap L ([0,T0 ε τ] Ω) + 1/2, a quick energy estimate on the equation satisfied by w yields w(t) L 2 (Ω) C N+1(1 + t) 1/4 (g ε (t)) N+1. As a result, u ε (t) u sh ( εt) N+1 L (1 + t) 1/4 C 2 (Ω A (t)) 1g ε (t) C j (g ε (t)) j (1 + t) 1/4 (C 1/2)g ε (t) for t T ε 0 τ, with τ > τ independent of ε. j=2 17 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Now we examine w = u ε u ap. Choosing N so that Nσ 0 > u ap L ([0,T0 ε τ] Ω) + 1/2, a quick energy estimate on the equation satisfied by w yields w(t) L 2 (Ω) C N+1(1 + t) 1/4 (g ε (t)) N+1. As a result, u ε (t) u sh ( εt) N+1 L (1 + t) 1/4 C 2 (Ω A (t)) 1g ε (t) C j (g ε (t)) j (1 + t) 1/4 δ 0 for t = T ε 0 T, with T > τ chosen independently of ε. j=2 18 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Comparison with the Dirichlet boundary condition. WKB expansion based on the Prandtl ansatz (3): O(1) boundary layers. The study time T ε 0 is shorter. One misses instability in L, ε 1/4 δ 0 replaces δ 0 in the final equality, but we do get instability of the vorticity. 19 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Robustness of the method. This linear instability implies nonlinear instability argument requires: a most unstable eigenmode for the linearised system, a resolvent estimate for the linearised system, and good energy estimates for the equation solved by the difference between the exact and approximate solutions. 10 Rousset & Tzvetkov, multiple references 11 arxiv 1502.05647 20 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Robustness of the method. This linear instability implies nonlinear instability argument requires: a most unstable eigenmode for the linearised system, a resolvent estimate for the linearised system, and good energy estimates for the equation solved by the difference between the exact and approximate solutions. Thus it has also been used to obtain transverse instability of solitary waves in numerous settings: NLS, KP-I, water-waves 10, Euler-Korteweg 11,... 10 Rousset & Tzvetkov, multiple references 11 arxiv 1502.05647 20 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

1 Stability and instability in the 2D incompressible model 2 21 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

3D isentropic, compressible Navier-Stokes equation on the half space Ω = {X = (x, y, z) R 2 ]0, + [}: { t ρ + div (ρu) = 0 ρ t u + ρu u = div Σ + ρf, (9) modelling the motion of a viscous, non-heat-conducting fluid with typically supersonic speeds. 22 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

3D isentropic, compressible Navier-Stokes equation on the half space Ω = {X = (x, y, z) R 2 ]0, + [}: { t ρ + div (ρu) = 0 ρ t u + ρu u = div Σ + ρf, (9) modelling the motion of a viscous, non-heat-conducting fluid with typically supersonic speeds. Σ is the stress tensor: Σ = εµ(ρ)( + t )u + ελ 0 (ρ)div ui 3 P(ρ)I 3. (10) The viscosity coefficients µ > 0 and µ + λ 0 > 0 can depend smoothly on the density. The pressure P(ρ) is given by a barotropic law: P(ρ) = kρ γ with k > 0 and γ > 1. 22 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

F is a smooth force term. We solve for t R with the fluid at rest for negative times: Navier boundary condition: (ρ, u, F ) = (1, 0, 0) for t < 0. u 3 z=0 = 0 and (ε 1 Σ n + au) τ z=0 = 0, where v τ = (v 1, v 2 ) is the tangential part of a 3D vector v, and a R. 23 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

F is a smooth force term. We solve for t R with the fluid at rest for negative times: Navier boundary condition: (ρ, u, F ) = (1, 0, 0) for t < 0. u 3 z=0 = 0 and [µ(ρ) z u τ ] z=0 = 2au τ z=0, (11) where v τ = (v 1, v 2 ) is the tangential part of a 3D vector v, and a R. 24 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Global existence is assured for weak solutions with constant viscosity coefficients, and only under certain conditions for variable coefficients 12. Strong solutions are only local in time 13. In the constant coefficient case, convergence of weak solutions to a strong solution of the Euler equations has been obtained by Sueur. 12 Bresch, Desjardins & Gérard-Varet, J. Math. Pures Appl., 2007 13 Solonnikov, Zap. Naučn. Sem. LOMI Steklov., 1976 25 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Method: a priori estimates in conormal Sobolev spaces, similar to works on the incompressible Navier-Stokes system with Navier BCs (fixed boundary and free surface 14 ). This avoids constructing approximate solutions. 14 Masmoudi & Rousset, Arch. Ration. Mech. Anal., 2012 & preprint 26 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Method: a priori estimates in conormal Sobolev spaces, similar to works on the incompressible Navier-Stokes system with Navier BCs (fixed boundary and free surface 14 ). This avoids constructing approximate solutions. Conormal regularity is necessary to get uniform estimates for boundary layers. Indeed, we expect the ansatz (ρ, u)(t, y, z) (r, v)(t, y, z) + ε(r, V )(t, y, z/ ε), with (r, v) solving the isentropic Euler equation: Lipschitz control at best. 14 Masmoudi & Rousset, Arch. Ration. Mech. Anal., 2012 & preprint 26 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Consider this set of derivatives that are tangent to Ω: Z 0 = t, Z 1,2 = y1,2, Z 3 = φ(z) z = z 1 + z z, 27 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Consider this set of derivatives that are tangent to Ω: Z 0 = t, Z 1,2 = y1,2, Z 3 = φ(z) z = z 1 + z z, For m N and p {2, }, we introduce the notations f (t) 2 m,p = (Z α f )(t) 2 L p (Ω), for T 0. α N 4, α m and f m,p,t = sup f (t) m,p t [0,T ] 27 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The solutions U = (ρ 1, u) we consider will satisfy T E m (T, U) = U 2 m,2,t + zu τ 2 m 1,2,T + z (ρ, u 3 )(s) 2 m 1,2 0 ds T + z u τ 2 1,,T + ( z ρ, tz ρ)(s) 2 1, 0 ds < +. 28 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The solutions U = (ρ 1, u) we consider will satisfy T E m (T, U) = U 2 m,2,t + zu τ 2 m 1,2,T + z (ρ, u 3 )(s) 2 m 1,2 0 ds T + z u τ 2 1,,T + ( z ρ, tz ρ)(s) 2 1, The force term will satisfy, for every T 0, N m (T, F ) = and F (t, x) = 0 for t < 0. 0 ds < +. sup ( F (t) 2 m,2 + zf (t) 2 m 1,2 + F (t) 2 2, ) < + t [ T,T ] 28 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Theorem 2.1 (submitted, 2014 15 ) Under the conditions stated previously, for ε 0 > 0 and m 7, there exists T > 0 such that, for every 0 < ε < ε 0, there exists a unique solution U to the isentropic Navier-Stokes system (9) with Navier boundary condition (11) in the class of functions satisfying E m (T, U) < +. Moreover, there is uniformly no vacuum on this time interval. The family of solutions (U ε ) 0<ε<ε0 constructed above converges in L 2 ([0, T ] Ω) and in L ([0, T ] Ω) towards the unique solution of the isentropic Euler equation ((9) with ε = 0 and the non-penetration boundary condition) in the same class. 15 arxiv 1410.2811 29 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Remarks. The results extend to less regular forces or different initial conditions, providing the compatibility conditions yield uniform bounds on N m (0, U). 16 Wang & Williams, Ann. Inst. Fourier, 2012 30 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Remarks. The results extend to less regular forces or different initial conditions, providing the compatibility conditions yield uniform bounds on N m (0, U). Wang and Williams 16 proved similar results in the constant coefficient case by constructing approximate solutions, and needs prior existence of the solution to the Euler system. 16 Wang & Williams, Ann. Inst. Fourier, 2012 30 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Idea of the proof: a priori estimates + compactness argument. 31 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Idea of the proof: a priori estimates + compactness argument. Theorem 2.2 Under the conditions of Theorem 2.1, for ε 0 > 0, there exist T > 0 and a positive, increasing function Q : R + R + such that, for every 0 < ε < ε 0 and t T, t E m (t, U) + z u 3 2 1,,t + ε u(s) 2 m,2 + 2 u τ (s) 2 m,2 0 ds (t + ε)q(e m (t, U) + N m (t, F )). 31 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Idea of the proof: a priori estimates + compactness argument. Theorem 2.2 Under the conditions of Theorem 2.1, for ε 0 > 0, there exist T > 0 and a positive, increasing function Q : R + R + such that, for every 0 < ε < ε 0 and t T, t E m (t, U) + z u 3 2 1,,t + ε u(s) 2 m,2 + 2 u τ (s) 2 m,2 0 ds (t + ε)q(e m (t, U) + N m (t, F )). Throughout, ρ is assumed to be uniformly bounded from below - this is checked by a bootstrap argument using the fact that, along the flow of u, t ρ = ρdiv u, and div u is bounded. 31 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The energy estimates on U use the symmetrisable hyperbolic-parabolic structure of the equation. The system on U α = Z α U is (DA 0 )(ρ) t U α + 3 (DA j )(U) xj U α ε(0, div (µ u) α +(µ+λ 0 ) div u α )) j=1 = (0, (ρf + εσ( U)) α ) C α, with the DA j (U) symmetric matrices, σ( U) the remainder of the order-two terms and C α containing commutators between Z α and the operators involved in the equation (e.g. [Z α, DA j xj ]U). Basic inequalities are used in the energy estimate. 32 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

We now need estimates on z U: these are obtained component by component. 33 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

We now need estimates on z U: these are obtained component by component. The norms of z u 3 are easily controlled by using the equation of conservation of mass: ρ z u 3 = t ρ u ρ ρ( x u 1 + y u 2 ) 33 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

To deal with z u τ, we use a modified vorticity W = (curl u) τ 2a µ(ρ) u τ, so that we have ρ t W + ρu W εµ W = H and W = 0 on the boundary. 34 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

To deal with z u τ, we use a modified vorticity W = (curl u) τ 2a µ(ρ) u τ, so that we have ρ t W + ρu W εµ W = H and W = 0 on the boundary. The commutator estimates are standard, except for [Z α, ρu 3 z ]W, in which a Hardy-type inequality is used to avoid the appearance of 2 u. 34 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

To deal with z u τ, we use a modified vorticity W = (curl u) τ 2a µ(ρ) u τ, so that we have ρ t W + ρu W εµ W = H and W = 0 on the boundary. The commutator estimates are standard, except for [Z α, ρu 3 z ]W, in which a Hardy-type inequality is used to avoid the appearance of 2 u. L bounds are obtained using a variant of the maximum principle, on W in the constant coefficient case, on (curl u) τ otherwise. 34 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The mass equation is used to control R = z ρ: lε( t R + u R) + γp(ρ)r = h, with l = 2µ + λ 0. The conormal energy estimates are L 2 in time. 35 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

The mass equation is used to control R = z ρ: with l = 2µ + λ 0. lε( t R + u R) + γp(ρ)r = h, The conormal energy estimates are L 2 in time. By following the flow of u, we reduce the equation to the ODE lε t R + γpr = h, and the Duhamel formula yields the L estimates. 35 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

1 Stability and instability in the 2D incompressible model 2 36 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

A recent preprint obtains a similar result 17 using the natural energy structure of the Navier-Stokes system to get estimates on (u, P), using H 1 co bounds on P to control ρ 1,. 17 Wang, Xin & Yong, arxiv 1501.01718 37 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

A recent preprint obtains a similar result 17 using the natural energy structure of the Navier-Stokes system to get estimates on (u, P), using H 1 co bounds on P to control ρ 1,. They thus exhibit weaker boundary layers on ρ compared to u, and, on the domain T (0, 1), they obtain the rates of convergence: U ε (t) U E (t) L 2 Cε 2 and U ε (t) U E (t) L Cε 2/5. 17 Wang, Xin & Yong, arxiv 1501.01718 37 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Try the same approach as ours on the full Navier-Stokes equation: t ρ + div (ρu) = 0 ρ t u + ρu u = div Σ + ρf ρ t θ + ρu θ + Pdiv u = εκ θ + Σ : u + ρu F. (12) Internal energy proportional to the temperature θ Neumann boundary condition for θ Ideal gas pressure law P(ρ, θ) = kρθ 38 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Disaster strikes! The component-by-component examination of the normal derivatives has worked for the incompressible and isentropic models, but fails in this case! The problem: the terms involving div u in the equations on ρ and θ. The field z V = ( z ρ, z u 3, z θ) cannot be decoupled. 39 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Disaster strikes! The component-by-component examination of the normal derivatives has worked for the incompressible and isentropic models, but fails in this case! The problem: the terms involving div u in the equations on ρ and θ. The field z V = ( z ρ, z u 3, z θ) cannot be decoupled. Work in progress: we can get Hco m 1 bounds: Green function study? W 1, co estimates on z V, but we still need 39 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions

Thank you for your attention. 40 / 40 Matthew Paddick Boundary layers for Navier-Stokes equations with slip boundary conditions