two dimensional water waves in holomorphic coordinateslmichel/confgilles/slides/tataru.pdfform...

30
Two dimensional water waves in holomorphic coordinates Daniel Tataru University of California, Berkeley Nice, June 2014 This is joint work with Mihaela Ifrim and, in part, with John Hunter D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 1 / 30

Upload: ngothuy

Post on 10-Aug-2019

213 views

Category:

Documents


0 download

TRANSCRIPT

Two dimensional water waves in holomorphiccoordinates

Daniel Tataru

University of California, Berkeley

Nice, June 2014

This is joint work with Mihaela Ifrim and, in part, with John Hunter

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 1 / 30

Two dimensional water-waves

The setting:

inviscid incompressible fluid flow (Euler equations)

infinite bottom

free boundary

irrotational

with gravity or surface tension

periodic or nonperiodic setting

Question:

Study long time solutions (i.e. the stability of the trivial steady state)in two settings:

lifespan bounds for small data

global solutions for small localized data

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 2 / 30

The standard formulation

Fluid domain: Ω(t), free boundary Γ(t).Velocity field u, pressure p, gravity g, surface tension σ.Euler equations in Ω(t):

ut + u · ∇u = ∇p− gjdiv u = 0

curl u = 0

u(0, x) = u0(x)

Boundary conditions on Γt:∂t + u · ∇ is tangent to

⋃Γt (kinematic)

p = −2σH on Γt (dynamic)

H = mean curvature of the boundary.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 3 / 30

Irrotational flows

Velocity potential

u = ∇φ, ∆φ = 0 in Ωt

Dynamic boundary condition:

φt +1

2|∇φ|2 + y = 0 on Γt

Equations reduced to the boundary in Eulerian formulation:∂tη −G(η)φ = 0

∂tφ+ gη − σH(η) +1

2|∇φ|2 − 1

2

(∇η∇ψ +G(η)ψ)2

1 + |∇η|2= 0.

H(η) = div

(∇η

1 + |∇η|2

), G(η) = Dirichlet to Neuman operator

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 4 / 30

Overview:

Features:

Hamiltonian structure

dispersive character

few resonant interactions (g/t)

symmetries: translation, scaling (g/t)

gauge freedom = choice of coordinates (e.g. Euler or Lagrange)

Difficulties:

quasilinear evolution

degenerate higher order hyperbolic evolution

nonlocal equation

quadratic interactions

many nonresonant interactions (g+t)

solitons (g+t)

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 5 / 30

Overview:Earlier work:

Local well-posedness in Sobolev spaces: Nalimov ’74, Yoshihara,Wu, Lannes, Coutand-Shkoller, Lindblad, Shatah-Zeng ...Alazard-Burq-Zuily.

Long time solutions for small localized data (g): Wu ’09 (almostglobal) Ionescu-Pusateri ’13, (global), Alazard-Delort ’13 (global)

Global solutions for small localized data in 3-d:Germain-Masmoudi-Shatah ’10(g) , ’12 (t), Wu ’10 (g).

Our goals:

To use the formulation of the equations in holomorphiccoordinates (Nalimov ’74) in order to provide a simpler approachto the local problem

To introduce a modified energy method, improving on the normalform method, which yields an easier route to long time solutions

To develop a more robust and simpler way, based on wave-packets,in order to obtain global solutions.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 6 / 30

Holomorphic coordinates:Holomorphic coordinates:

Z : =z ≤ 0 → Ωt, α− iβ → Z(α+ iβ)

Boundary condition at infinity:

Z(α)− α→ 0 (nonperiodic) Z(α)− α periodic (periodic)

Free boundary parametrization:

Z : R→ Ωt, α→ Z(α)

Perturbation of steady state:

W = Z − α

Holomorphic velocity potential (v = ∇φ):

Q = φ+ iψ

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 7 / 30

Water wave equations in holomorphic coords.

P - Projection onto negative wavenumbers

Equations for (W = Z − α,Q):Wt + F (1 +Wα) = 0,

Qt + FQα + P

[|Qα|2

J

]− igW+ iσP

[Wαα

J1/2(1 +Wα)− Wαα

J1/2(1 + Wα)

]= 0.

where

F = P

[Qα − Qα

J

], J = |1 +Wα|2.

Conserved energy (Hamiltonian):

E(W,Q) =

∫=(QQα) +

1

2g(|W |2 −<(W 2Wα)

)+

1

4σ(J

12 − 1−<Wα) dα

Symmetries:

Translations in α and t.

Scaling (W (t, α), Q(t, α))→ (λ−2W (λt, λ2α), λ−3Q(λt, λ2α))

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 8 / 30

Physical parameters:

R =Qα

1 +Wα- velocity field on the free boundary

b = 2<P[

R

1 + Wα

]- advection coeff. (high freq. velocity limit)

a = 2=P [RRα] > 0 - normal derivative of the pressure = g + a

Alternate quasilinear system for diagonal variables (W = Wα, R):Wt + P [bWα] + P

[(1 + W)Rα

1 + W

]= G

Rt + P [bRα] + iP

[g + a

1 + W

]+ iσ

1

1 + WP

[Wα

J1/2(1 + W)

= K,

where

G = (1 + W)P

[Rα

1 + W+

RWα

(1 + W)2

]+ [P,W]

(Rα

1 + W+

RWα

(1 + W)2

)K = − P

[RRα

]+ iσP

[Wα

J1/2(1 + W)

represent perturbative terms in the equation.D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 9 / 30

Cubic lifespan bounds

Theorem

Consider the two dimensional differentiated water wave equation withinitial data in (W0, R0) ∈ H1 of size ε,

‖(W0, R0)‖H0∩H1 . ε (gravity waves)

respectively

‖(W0, R0)‖H0∩H2 . ε (capillary waves)

Then the solutions have a lifespan of at least

Tε ≈ ε−2

Proof idea: quasilinear modified energy method

Bounds for these and higher norms propagate on same timescale.

The result applies both in the periodic and nonperiodic setting.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 10 / 30

Normal forms and long time existenceQuestion: Find improved lifespan estimates for small data solutions.

(i) Equations with quadratic nonlinearities:

d

dtE(u) . ‖u‖E(u)

For data ‖u(0)‖ = ε 1 this leads by Gronwall to a lifespan Tε ≈ ε−1

(ii) Equations with cubic nonlinearities:

d

dtE(u) . ‖u‖2E(u)

For data ‖u(0)‖ = ε 1 this leads by Gronwall to a lifespan Tε ≈ ε−2

(iii) Normal form method (Shatah ’85): transform an equation with aquadratic nonlinearity into one with a cubic one via a normal formtransformation,

u→ v = u+B(u, u) + higher

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 11 / 30

Dispersive character (gravity waves)Linearization around zero:

∂tw + rα = 0

∂tr − igw = 0

Energy:‖(w, r)‖2H0 = g‖w‖2L2 + ‖r‖2

H12

Dispersion relation:τ = ±

√g|ξ|, ξ ≤ 0

Group velocity of waves:

v = ±√g

2√|ξ|

Quasilinear model: (∂t + b∂α)w + rα = 0

(∂t + b∂α)r − i(g + a)w = 0

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 12 / 30

Dispersive character (gravity waves)Linearization around zero:

∂tw + rα = 0

∂tr − igw = 0

Energy:‖(w, r)‖2H0 = g‖w‖2L2 + ‖r‖2

H12

Dispersion relation:τ = ±

√g|ξ|, ξ ≤ 0

Group velocity of waves:

v = ±√g

2√|ξ|

Bilinear resonant interactions:

|ξ|12 ± |η|

12 ± |ξ + η|

12 = 0 =⇒ ξη(ξ + η) = 0

No resonances away from frequency zero !D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 13 / 30

Normal forms for gravity water waves

Existence of a normal form transformation is related to the absence ofresonant bilinear interactions. For 2-d gravity waves in holomorphiccoordinates, such a normal form transformation exists and is given by

W = W − 2M<WWα,

Q = Q− 2M<WQα,

The normal variables solve an equation of the form∂tW + Qα = cubic and higher

∂tQ− iW = cubic and higher

However, the cubic and higher nonlinearities also contain higherderivatives, so one cannot close the energy estimates. This is related tothe fact that the normal form transformation is not invertible, andfurther to the fact that the water wave equation is quasilinear, ratherthan semilinear.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 14 / 30

Dispersive character (capillary waves)Linearization around zero:

∂tw + rα = 0

∂tr + iσwαα = 0

Energy:‖(w, r)‖2H1 = σ‖wα‖2L2 + ‖r‖2

H12

Dispersion relation:

τ = ±√σ|ξ|

32 , ξ ≤ 0

Group velocity of waves:

v = ±3

2

√σ|ξ|

Quasilinear model: (∂t + b∂α)w + rα = 0

(∂t + b∂α)r − iσJ−12wαα = 0

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 15 / 30

Dispersive character (capillary waves)Linearization around zero:

∂tw + rα = 0

∂tr + iσwαα = 0

Energy:‖(w, r)‖2H1 = σ‖wα‖2L2 + ‖r‖2

H12

Dispersion relation:

τ = ±√σ|ξ|

32 , ξ ≤ 0

Group velocity of waves:

v = ±3

2

√σ|ξ|

Bilinear resonant interactions:

|ξ|32 ± |η|

32 ± |ξ + η|

32 = 0 =⇒ ξη(ξ + η) = 0

No resonances away from frequency zero !D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 16 / 30

Normal forms for capillary water waves

Existence of a normal form transformation is related to the absence ofresonant bilinear interactions. For 2-d capillary waves in holomorphiccoordinates, such a normal form transformation exists and is given by

W = W + L1(W,W ) + L0(Q,Q)

Q = Q+ L1(W,Q)

where L0, L1 stand for bilinear paraproduct type operators of order0, 1. The normal variables solve an equation of the form

∂tW + Qα = cubic and higher

∂tQ− iW = cubic and higher

as before. Again, the cubic and higher nonlinearities also containhigher derivatives, so one cannot close the energy estimates.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 17 / 30

The modified energy methodIdea: Modify the energy rather than the equation in order to get cubicenergy estimates.

Step 1: Construct a cubic normal form energy

EnNF (W,Q) = (quadratic+ cubic)(‖W (n)‖2L2 + ‖Q(n)‖2H

12)

Thend

dtEnNF (W,Q) = quartic+ higher

Here higher derivatives arise on the right, making it impossible to close.

Step 2: Switch EnNF (W,Q) to diagonal variables EnNF (W, R).

Step 3: To account for the fact that the equation is quasilinear,replace the leading order terms in EnNF (W, R) with their naturalquasilinear counterparts to obtain a good cubic quasilinear energyEn(W, R).D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 18 / 30

Energy estimates for gravity waves

We use the control norms:

A = ‖W‖L∞ + ‖W(1 + W)−1‖L∞ + ‖D12R‖L∞ , (scale invariant)

B = ‖D12W‖BMO + ‖Rα‖BMO (controlled by ‖(W, R)‖H1)

and construct energy functionals En(W, R) with these properties:

(i) Energy equivalence:

En(W, R) = (1 +O(A))‖(W, R)‖2Hn

(ii) Cubic scale invariant energy estimate:

d

dtEn(W, R) .A ABE

n(W, R)

Similar H0 bound for linearized eqn (w, r) = (w, q −Rw).

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 19 / 30

Energy estimates for capillary waves

We use the control norms:

A = ‖W‖L∞ + ‖W(1 + W)−1‖L∞ + ‖D−12R‖L∞ , (scale invariant)

B = ‖D32W‖BMO + ‖Rα‖BMO (controlled by ‖(W, R)‖H2)

and construct energy functionals En(W, R) with these properties:

(i) Energy equivalence:

En(W, R) = (1 +O(A))‖(W, R)‖2Hn

(ii) Cubic scale invariant energy estimate:

d

dtEn(W, R) .A ABE

n(W, R)

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 20 / 30

Water waves on the real line: heuristics #1Dispersive decay for the linear equation

Wt +Qα = 0

Qt − iW = 0

with smooth localized data of size ε:

|W |+ |Q| . εt−12

This shows that one would expect our control norms to decay like

A,B ≈ εt−12

Henced

dtE . ε2t−1E

which provides uniform bounds up to an exponential time

Tε = ecε−2

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 21 / 30

Water waves on the real line: heuristics #2Wt + Qα = cubic(W , Q) + higher

Qt − iW = cubic(W , Q) + higher

Based on the dispersion relation and the group velocity of waves

τ = ±√|ξ|, v = ± 1

2√|ξ|

for small data localized near 0, one expects asymptotics

(W , Q) ≈ γ(t, α/t)t−12 e

it2

(t

2α, 1

)with a better function γ. Substituting in the above equations, we getthe asymptotic evolution

γ(t, v) = ct−1γ(t, v)|γ(t, v)|2 +O(εt−1−), γ(1, α) = O(ε)

The global boundedness of gamma follows.D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 22 / 30

New tool: scaling symmetryScaling symmetry for gravity waves:

(W (t, α), Q(t, α))→ (λ−2W (λt, λ2α), λ−3Q(λt, λ2α))

Generator:

S(W,Q) = ((S − 2)W, (S − 3)Q), S = t∂t + 2α∂α

S(W,Q) solves the linearized equation, good energy estimates.

Scaling symmetry for capillary waves:

(W (t, α), Q(t, α))→ (λ−1W (λ32 t, λα), λ−

12Q(λ

32 t, λα)).

Generator:

S(W,Q) = ((S − 1)W, (S − 1

2)Q), S =

3

2t∂t + α∂α.

S(W,Q) solves the linearized equation, good energy estimates.D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 23 / 30

Global gravity waves for small localized dataWeighted function space WH for localized data:

‖(W,Q)(t)‖2WH = ‖(W,Q)‖2H0 + ‖(W, R)‖2H6 + ‖AS(W,Q)‖2H0

Theorem

Assume that the initial data for the water wave equation satisfies:

‖(W,Q)(0)‖WH . ε

Then the solution exists globally in time, with energy bounds

‖(W,Q)(t)‖WH . εtCε2

and pointwise decay

A(t) +B(t) .ε√t

Original result by Ionescu-Pusateri and Alazard-Zuily (with longerproofs and with more regularity)

.D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 24 / 30

Global capillary waves for small localized data

Weighted function space WH for localized data:

‖(W,Q)(t)‖2WH = ‖(W,Q)‖2Hδ∩H10 + ‖S(W,Q)‖2Hδ∩H1

Theorem

Assume that the initial data for the water wave equation satisfies:

‖(W,Q)(0)‖WH . ε

Then the solution exists globally in time, with energy bounds

‖(W,Q)(t)‖WH . εtCε2

and pointwise decay

C(t) := ‖(D12W,Q)‖W 4,∞ .

ε√t

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 25 / 30

Bootstrap argument for gravity waves:1. Energy estimates:

‖(W,Q)(t)‖WH . e∫ T0 AB(s)ds‖(W,Q)(0)‖WH . εtCε

2

based on the cubic energy functionals

2. Initial pointwise bounds:

(A+B)(t) . t−12ω(t, α/t)‖(W,Q)(t)‖WH, ω(t, v) = t−

118 +(v+v−1)−

12

akin to the Klainerman-Sobolev inequalities, only harderuses the normal form variablesbetter decay for v away from 1, reduces problem to t−δ < v < tδ.

3. Final pointwise bounds:

A(t) +B(t) . εt−12

Uses the method of testing with wave-packets to derive anasymptotic equation.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 26 / 30

Bootstrap argument for capillary waves

All is similar except for the energy estimates, which are done in twoways, both for the full equation and for the linearization:

1. High frequency energy estimates (e.g. in H10)

Uses the modified energy metod to obtain a cubic energyfunctional.

2. Low frequency energy estimates (i.e. in Hδ)Uses the normal form metod to obtain a cubic energy functional.

Not all cubic terms in the normal form equation are good, but thebad ones are nonresonant and can be eliminated with a furthercubic normal form correction.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 27 / 30

Asymptotic equations in NLS context

A. Hayashi-Naumkin, refined by Kato-Pusateri; derive an asymptoticequation for the Fourier transform of the solutions,

d

dtu(t, ξ) = λit−1 u(t, ξ)|u(t, ξ)|+OL∞(t−1−ε).

B. Lindblad-Soffer; derive an asymptotic equation in the physicalspace along rays,

(t∂t + x∂x)u(t, x) = λitu(t, x)|u(t, x)|2 +OL∞(t−ε).

C. Deift-Zhou used the inverse scattering method to obtain longrange asymptotics

D. Ifrim-T. wave packet testing: test the NLS solution with anapproximate wave packet type linear wave.

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 28 / 30

Testing by wave packets

Wave packets:

q = vχ

(v(α− vt)√

t

)eit2

4α , w = −iqt

The γ function:γ(t, v) = 〈(W , Q), (w,q)〉H0

Good approximation (also for derivatives):

(W , Q) = γ(t, α/t)t−12 e

it2

(t

2α, 1

)+O(t−

58 )

Asymptotic equation for γ:

γ(t, v) = ct−1γ(t, v)|γ(t, v)|2 +O(εt−1−), γ(1, α) = O(ε)

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 29 / 30

Happy birthday Gilles !

D. Tataru (UC Berkeley) 2-d water waves Nice, June 2014 30 / 30