Transverse stability of periodic waves in water-wave models
Mariana Haragus
Institut FEMTO-ST and LMB Universit´ e Bourgogne Franche-Comt´ e, France
Transverse stability of periodic waves in water-wave models Mariana - - PowerPoint PPT Presentation
Transverse stability of periodic waves in water-wave models Mariana Haragus Institut FEMTO-ST and LMB Universit e Bourgogne Franche-Comt e, France ICERM, April 28, 2017 Water-wave problem gravity/gravity-capillary waves
Institut FEMTO-ST and LMB Universit´ e Bourgogne Franche-Comt´ e, France
three-dimensional inviscid fluid layer constant density gravity/gravity and surface tension irrotational flow
depth at rest h
velocity potential φ; free surface
h + ηmean curvature K =
1+η2
x+η2 z
1+η2
x+η2 z
parameters ρ,
g, σ, hcharacteristic length
hcharacteristic velocity
inverse square of the Froude number
Weber number
2
2
variable domain (free surface) nonlinear boundary conditions
symmetries, Hamiltonian structures many particular solutions
1 3
1 3
longitudinal co-periodic perturbations transverse periodic perturbations
a periodic wave
U∗( x) is an equilibrium z x1 the system
U z = DU t + F( U) is reversible/Hamiltonian;2 the linear operator L =
F ′( U∗) possesses a pair of3 the operators D and L are closed in X with D(L) ⊂ D(D);
1 For any λ ∈ R sufficiently small, the linearized system
U z = DU t + L U2
U∗ is transversely linearly unstable.boundary conditions
xu + u + 1
yu
xu + u + 1
yu
boundary conditions
boundary conditions
boundary conditions
space of symmetric functions (
x → − x)e (0, 2π) × L2 e(0, 2π) × H1
ε + L1 ε
L0
ε
η ω φ ξ = ω β −εk2
a βηxx + (1 + ǫ)η − kaφx |y=1
ξ −εk2
a φxx − φyy
, L1
ε
η ω φ ξ = g1 g2 G1 G2 g1 = (1 + εk2
a η2 ax )1/2
β
1 1 + εηa 1 yφay ξ dy
ω β g2 = 1
a φax φx −
φay φy (1 + εηa)2 + εφ2
ay η
(1 + εηa)3 − ε3k2
a y2η2 ax φay φy
(1 + εηa)2 − ε3k2
a y2ηax φ2 ay ηx
(1 + εηa)2 + ε3k2
a y2η2 ax φ2 ay η
(1 + εηa)3 +
a yφay φx + εk2 a yφax φy −
2ε2k2
a y2ηax φay φy
1 + εηa − ε2k2
a y2φ2 ay ηx
1 + εηa + ε3k2
a y2ηax φ2 ay η
(1 + εηa)2
+ εk2
a βηxx − εk2 a β
(1 + ε3k2
a η2 ax )3/2
G1 = − εηaξ 1 + εηa + (1 + ε3k2
a η2 ax )1/2
β(1 + εηa)
1 1 + εηa 1 yφay ξ dy
G2 =
(1 + εηa) + εφaη (1 + ηa)2
− ε2k2
a [ηaφx + φax η − yφay ηx − yηax φy ]x
+ ε2k2
a
ε2y2η2
ax φay η
(1 + εηa)2 − εy2η2
ax φy
1 + εηa − 2εy2ηax φay ηx 1 + εηa
scaling λ = εℓ,
decomposition φ( x,
y) = φ1( x) + φ2( x, y)ω = β (1 + ǫ3η⋆2
x )1/2 (η† + ikη) −
1 1 + ǫη⋆ 1 yΦ⋆
y ξdy,
ξ = (1 + ǫη⋆)(Φ† + ikΦ) − ǫyΦ⋆
y (η† + ikη)
(1 + ǫ) ǫ2 η − 1 ǫ2 Φx |y=1 − 1 ǫ βηxx − ikβ(hǫ
1 + ikη) = hǫ 2
− 1 ǫ Φxx − 1 ǫ2 Φyy − ik(Hǫ
1 + ikΦ) = Hǫ 2 ,
hǫ
2 = ω† − gǫ 2 ,
Hǫ
2 = ξ† − Gǫ 2
hǫ
1 =
ω β − ikη = − 1 β(1 + ǫη⋆) 1 yΦ⋆
y [−ǫyΦ⋆ y (ikη + η†) + (1 + ǫη⋆)(ikΦ + Φ†)]dy
+
(1 + ǫ3η⋆2
x )1/2 − 1
η† (1 + ǫ3η⋆2
x )1/2 ,
Hǫ
1 = ξ − ikΦ
= (1 + ǫη⋆)Φ† + ikǫη⋆Φ − ǫyΦ⋆
y (η† + ikη).
Bǫ(η, Φ) = −ǫηx + Bǫ
0 + Bǫ 1 ,
Bǫ = ǫη⋆Φy 1 + ǫη⋆ + ǫΦ⋆
y η
(1 + ǫη⋆)2
, Bǫ
1
= ǫ2η⋆
x Φx + ǫ2Φ⋆ x ηx +
ǫ4η⋆2
x Φ⋆ y η
(1 + ǫη⋆)2 − ǫ3η⋆2
x Φy
1 + ǫη⋆ − −ˆ Φyy + q2 ˆ Φ = ǫ2( ˆ Hǫ
2 + ik ˆ
Hǫ
1 ),
0 < y < 1 ˆ Φy = 0, y = 0 ˆ Φy − ǫµ2 ˆ Φ 1 + ǫ + βq2 = − ǫ3iµ(ˆ hǫ
2 + ikβˆ
hǫ
1)
1 + ǫ + βq2 + ˆ Bǫ
0 + ˆ
Bǫ
1 , y = 1
G(y, ζ) = cosh qy cosh q (1 + ǫ + βq2) cosh q(1 − ζ) + (ǫµ2/q) q2 − (1 + ǫ + βq2)q tanh q − ǫ cosh qζ cosh q (1 + ǫ + βq2) cosh q(1 − y) + (ǫµ2/q) q2 − (1 + ǫ + βq2)q tanh q −
ˆ Φ1 = 1 + ǫ ǫ2(k2(1 + ǫ) + µ2 + (β − 1/3)µ4) × 1 ǫ2( ˆ ξ† − iµ ˆ Gǫ
2,2 + ik ˆ
Hǫ
1 )dζ − ǫq2
1 ˆ pǫ
2 dζ
− ǫ3iµ(ˆ hǫ
2 + ikβˆ
hǫ
1 )
1 + ǫ + βq2 + ǫ2µ2 ˆ pǫ
2 |ζ=1
1 + ǫ + βq2
ˆ Φ2 = − 1 G1( ˆ ξ† − iµ ˆ Gǫ
2,2 + ik ˆ
Hǫ
1 )dζ −
1 G1ζ ˆ Gǫ
2,1dζ +
1 (ǫk2 + µ2)G1 ˆ pǫ
2 dζ + ǫˆ
pǫ
2
− G1|ζ=1
ǫiµ(ˆ hǫ
2 + ikβˆ
hǫ
1)
1 + ǫ + βq2 + µ2 ˆ pǫ
2 |ζ=1
1 + ǫ + βq2
ˆ Φ = − 1 Gǫ2( ˆ ξ† − ıµ ˆ Gǫ
2,2 + ık ˆ
Hǫ
1 )dζ −
1 Gζǫ2 ˆ Gǫ
2,1dζ +
ǫ3ıµG|ζ=1(ˆ hǫ
2 + ıkβˆ
hǫ
1 )
1 + ǫ + βq2 + 1 ǫq2G ˆ pǫ
2 dζ + ǫˆ
pǫ
2 −
ǫ2µ2 1 1 ǫq2G ˆ pǫ
2 dζ + ǫˆ
pǫ
2 −
ǫ2µ2G|ζ=1 ˆ pǫ
2 |ζ=1
1 + ǫ + βq2 = 1 Gǫˆ pǫ
2ζζdζ − ǫG|ζ=1 ˆ
pǫ
2ζ|ζ=1 =
1 Gǫ2( ˆ Gǫ
2,0)ζζdζ − G|ζ=1 ˆ
Bǫ
0 ,
ˆ Φ1 + ˆ Φ2 = − 1 Gǫ2( ˆ ξ† − ıµ ˆ Gǫ
2,2 + ık ˆ
Hǫ
1 )dζ −
1 Gζǫ2 ˆ Gǫ
2,1dζ
+ ǫ3ıµG|ζ=1(ˆ hǫ
2 + ıkβˆ
hǫ
1)
1 + ǫ + βq2 + 1 ǫq2G ˆ pǫ
2 dζ + ǫˆ
pǫ
2 −
ǫ2µ2G|ζ=1 ˆ pǫ
2 |ζ=1
1 + ǫ + βq2 ,
ıµˆ hǫ
2
= ıµˆ ω† + ıµF
1 ǫ2 1
x Φx −
Φ⋆
y Φy
(1 + ǫη⋆)2 + ǫΦ⋆2
y η
(1 + ǫη⋆)3 − ǫ3y2η⋆2
x Φ⋆ y Φy
(1 + ǫη⋆)2 − ǫ3y2η⋆
x Φ⋆2 y ηx
(1 + ǫη⋆)2 + ǫ4y (1 + µ2 ǫ F 1
y Φx + yΦ⋆ x Φy −
2ǫy2η⋆
x Φy Φ⋆ y
1 + ǫη⋆ − ǫy2Φ⋆2
y ηx
1 + ǫη⋆ + ǫ2y2η⋆
x Φ⋆2 y η
(1 + ǫη⋆)2
βµ2 ǫ F
F−1
hǫ
2
1 + ǫ + βq2
1 + ǫ + βq2 F[(Φ⋆
1x Φ1x )x ]
1 + ǫ + βq2 F 1 yΦ⋆
x Φ2y dy
1 1 + ǫ + βq2 F 1
2x Φ1x + Φ⋆ x Φ2x −
Φ⋆
y Φy
ǫ(1 + ǫη⋆)2 + Φ⋆2
y η
(1 + ǫη⋆)3 − ǫ2y2η⋆2
x Φ⋆ y Φy
(1 + ǫη⋆)2 − ǫ2y2η⋆
x Φ⋆2 y ηx
(1 + ǫη⋆)2 + ǫ3y2η⋆
x Φ⋆2 y η
(1 + ǫη⋆)3
ıµ 1 + ǫ + βq2 F 1
y Φx −
2ǫy2η⋆
x Φ⋆ y Φy
1 + ǫη⋆ − ǫy2Φ⋆2
y ηx
1 + ǫη⋆ + ǫ2y2η⋆
x Φ⋆2 y η
(1 + ǫη⋆)2
βıµ 1 + ǫ + βq2 F
(1 + ǫ3η⋆2
x )3/2 − ηx
+ F−1
ω† 1 + ǫ + βq2
−F−1
1 + ǫ + βq2 F[(Φ⋆
1x Φ1x )x ]
1 + ǫ + βq2 F 1 yΦ⋆
x Φ2y dy
F−1
1 + ǫ + βq2 F 1 yΦ⋆
x Φ2y dy
1 + ǫ + βq2 F
1x Φ2|y=1 −
1 Φ⋆
1x Φ2dy +
1 yΦ⋆
2x Φ2y
=
1 + ǫ + βq2 µ1/2F[Φ⋆
1x Φ2|y=1] −
1 1 + ǫ + βq2 1 (Φ⋆
1x Φ2)xdy
+ µ 1 + ǫ + βq2 1 yΦ⋆
2x Φ2y dy
= ǫ−1/4(L(Φ2))x + (L(Φ2, Φ2x , ǫ1/2Φ2y ))x , F−1
hǫ
2
1 + ǫ + βq2
1 + ǫ + βq2 F[(Φ⋆
1x Φ1x )x ]
+ ǫ−1/2(L(ǫΦx , ǫ2Φ2y , ǫ4η, ǫ3ηx )x + (L(ǫΦ1x, Φ2, Φ2x , Φ2y , ǫ2η, ǫ4ηx )x + H. F−1
hǫ
1
1 + ǫ + βq2
(L(Φ2, ǫ2η))x + ǫ2k2(L(Φ1))x + H, F−1
pǫ
2 |ζ=1
1 + ǫ + βq2
F−1
1 ˆ pǫ
2 dζ
k2L(ǫΦ2, ǫ2η) + (L(Φ2, Φ2x , ǫη, ǫηx ))x 1 (ξ† − (Gǫ
2,2)x + ıkHǫ 1 )dζ = (η⋆Φ1x )x + (Φ⋆ 1x η)x + (L(Φ2x, Φ2y , ǫη, ǫηx ))x +
(β − 1/3)Φ1xxxx − Φ1xx + k2(1 + ǫ)Φ1 = (η⋆Φ1x )x + (Φ⋆
1x η)x + F−1
1 + ǫ + βq2 F[(Φ⋆
1x Φ1x )x ]
η = F−1
Φ1 1 + ǫ + βq2
Bε,ℓ small relatively bounded perturbation of B0,ℓ
Bε,ℓ small relatively bounded perturbation of B0,ℓ
spectrum of ∂
xA∂ x is known (KP-I):perturbation arguments . . . . . .
longitudinal co-periodic perturbations transverse periodic perturbations
solutions of the Kawahara equation
depend analytically upon (a, c) ∈ (−a0, a0) × (−c0, c0) ka,c = k0(c) + c
k0(c) = 1+√1+4c
2
1/2 , k(a, c) =
n≥1
k2n(c)a2n pa,c(z) = ac cos(z) + c
m,n
(n, m ≥ 0, n + m ≥ 2, n − m = ±1) explicit Taylor expansions for
if −ω2 is an isolated eigenvalue then
u∗ is transverselyif −ω2 belongs to the essential spectrum
scaling:
z = k a, x,1 the linear operator Λ∂
z − B a, acting in L2(0, 2π) has a2 the periodic wave ϕ
a, is transversely linearly unstablethe operator Λ∂
z − B a, is realperturbation argument: the negative eigenvalue of −B
a,0 is a triple eigenvalue all other eigenvalues are negative
0 is a triple eigenvalue all other eigenvalues are negative
σ1(B
a, ) ⊂ V , V neighborhood of 0σ2(B
a, ) ⊂ {ν ∈ C ; Re ν < − m}0 is a triple eigenvalue all other eigenvalues are negative
0 is a triple eigenvalue all other eigenvalues are negative
a = 0use symmetries and show that 0 is a double eigenvalue third eigenvalue: compute an expansion for small
a, . . .fully localized/bounded perturbations
J skew-adjoint operator L self-adjoint operator
J skew-adjoint operator L self-adjoint operator
relies upon the existence of a conserved higher-order energy
speed
> 12π-periodic, even profile φ
satisfying the KdV equation v ′′( x) + v( x) + 3 v 2( x) = 0known explicitly!
2π-periodic coefficients in
xAnsatz
w( x, y, t) = eλ t+ ip y W ( x),2D bounded perturbations: space
C b(R) and p ∈ R.continuous spectrum . . .
skew-adjoint operator J (γ) = (∂
x + iγ)self-adjoint operator
find commuting operators M
, p(γ)show that suitable linear combination of M
, p(γ) andresulting operator satisfies the commutativity relation cannot obtain positive operators . . .
resulting operator satisfies the commutativity relation cannot obtain positive operators . . .
decompose:
MKdV is obtained from a higher order conserved functional:
compute MKP directly from the commutativity relation:
there exist constants
b such that the operatorsthe commutativity relation holds the general counting criterion implies that the spectra of
there exist constants
b such that the operatorsthe commutativity relation holds the general counting criterion implies that the spectra of