Fokas Methods applied to a Boundary Valued Problem for Conjugate - - PowerPoint PPT Presentation

fokas methods applied to a boundary valued problem for
SMART_READER_LITE
LIVE PREVIEW

Fokas Methods applied to a Boundary Valued Problem for Conjugate - - PowerPoint PPT Presentation

Fokas Methods applied to a Boundary Valued Problem for Conjugate Conductivity Equations Conference in memory of Gennadi Henkin Joint work with Slah Chaabi and Franck Wielonsky (Universit e Aix-Marseille) Moscow, the 15th of september, 2016


slide-1
SLIDE 1

Fokas Methods applied to a Boundary Valued Problem for Conjugate Conductivity Equations

Conference in memory of Gennadi Henkin

Joint work with Slah Chaabi and Franck Wielonsky (Universit´ e Aix-Marseille)

Moscow, the 15th of september, 2016

slide-2
SLIDE 2

Different motivations.

◮ Ω is a smooth domain in C.

slide-3
SLIDE 3

Different motivations.

◮ Ω is a smooth domain in C. ◮ L an elliptic differential operator in Ω

L = a(x, y) ∂2 ∂x2 + b(x, y) ∂2 ∂x∂y + c(x, y) ∂2 ∂y2 + +d(x, y) ∂ ∂x + e(x, y) ∂ ∂y + f (x, y)

slide-4
SLIDE 4

Different motivations.

◮ Ω is a smooth domain in C. ◮ L an elliptic differential operator in Ω

L = a(x, y) ∂2 ∂x2 + b(x, y) ∂2 ∂x∂y + c(x, y) ∂2 ∂y2 + +d(x, y) ∂ ∂x + e(x, y) ∂ ∂y + f (x, y)

◮ Lax-Milgram theorem gives the existence and unicity of the

solution to Lu = 0 if u is given on ∂Ω (and continuous).

slide-5
SLIDE 5

Different motivations.

◮ Ω is a smooth domain in C. ◮ L an elliptic differential operator in Ω

L = a(x, y) ∂2 ∂x2 + b(x, y) ∂2 ∂x∂y + c(x, y) ∂2 ∂y2 + +d(x, y) ∂ ∂x + e(x, y) ∂ ∂y + f (x, y)

◮ Lax-Milgram theorem gives the existence and unicity of the

solution to Lu = 0 if u is given on ∂Ω (and continuous).

◮ Construction of fundamental solutions gives a representation

formula in terms of u and the normal derivative ∂

nu on the

boundary ∂Ω : we need to much informations.

slide-6
SLIDE 6

Different motivations.

◮ Ω is a smooth domain in C. ◮ L an elliptic differential operator in Ω

L = a(x, y) ∂2 ∂x2 + b(x, y) ∂2 ∂x∂y + c(x, y) ∂2 ∂y2 + +d(x, y) ∂ ∂x + e(x, y) ∂ ∂y + f (x, y)

◮ Lax-Milgram theorem gives the existence and unicity of the

solution to Lu = 0 if u is given on ∂Ω (and continuous).

◮ Construction of fundamental solutions gives a representation

formula in terms of u and the normal derivative ∂

nu on the

boundary ∂Ω : we need to much informations.

◮ Is it possible to obtain a relation between u and ∂ nu on ∂Ω

without computing u in Ω ?

slide-7
SLIDE 7

One other motivation : What is a Tokamak ?

slide-8
SLIDE 8

Mathematical Formulation

slide-9
SLIDE 9

Mathematical Formulation

◮ Poloidal field ψ and its normal derivative ∂ nψ are known on

∂D(a, r), ψ = C on the boundary ∂P of the plasma and div 1 x ∇ψ

  • = 0

in D(a, r) \ P

slide-10
SLIDE 10

Mathematical Formulation

◮ Poloidal field ψ and its normal derivative ∂ nψ are known on

∂D(a, r), ψ = C on the boundary ∂P of the plasma and div 1 x ∇ψ

  • = 0

in D(a, r) \ P

◮ Generalized axisymetrical potential for α ∈ R :

div(xα∇u) = 0 ⇔ ∆u + α x ∂u ∂x = 0.

slide-11
SLIDE 11

The Fokas method Fokas (1997) : a simple example.

qt − qxx = 0 sur x ≥ 0, t ≥ 0 (H)

slide-12
SLIDE 12

The Fokas method Fokas (1997) : a simple example.

qt − qxx = 0 sur x ≥ 0, t ≥ 0 (H)

◮ Formulation in terms of Lax pair : splitting of this equation

into two ordinary differential equations, compatible if and only if q is a solution to (H) µx − ikµ = q µt − µxx = 0

slide-13
SLIDE 13

The Fokas method Fokas (1997) : a simple example.

qt − qxx = 0 sur x ≥ 0, t ≥ 0 (H)

◮ Formulation in terms of Lax pair : splitting of this equation

into two ordinary differential equations, compatible if and only if q is a solution to (H) µx − ikµ = q µt − µxx = 0

  • µx − ikµ = q

µt + k2µ = qx + ikq

slide-14
SLIDE 14

The Fokas method Fokas (1997) : a simple example.

qt − qxx = 0 sur x ≥ 0, t ≥ 0 (H)

◮ Formulation in terms of Lax pair : splitting of this equation

into two ordinary differential equations, compatible if and only if q is a solution to (H) µx − ikµ = q µt − µxx = 0

  • µx − ikµ = q

µt + k2µ = qx + ikq

  

  • µe−ikx+k2t

x

= qe−ikx+k2t

  • µe−ikx+k2t

t

= (qx + ikq)e−ikx+k2t

slide-15
SLIDE 15

The Fokas method Fokas (1997) : a simple example.

qt − qxx = 0 sur x ≥ 0, t ≥ 0 (H)

◮ Formulation in terms of Lax pair : splitting of this equation

into two ordinary differential equations, compatible if and only if q is a solution to (H) µx − ikµ = q µt − µxx = 0

  • µx − ikµ = q

µt + k2µ = qx + ikq

  

  • µe−ikx+k2t

x

= qe−ikx+k2t

  • µe−ikx+k2t

t

= (qx + ikq)e−ikx+k2t

d

  • µe−ikx+k2t

= e−ikx+k2t (qdx + (qx + ikq)dt)

slide-16
SLIDE 16

◮ ν = e−ikx+k2t (qdx + (qx + ikq)dt) is closed

slide-17
SLIDE 17

◮ ν = e−ikx+k2t (qdx + (qx + ikq)dt) is closed ◮ q(x, t) and qx(x, t) decay to 0 at ∞ sufficiently fast.

slide-18
SLIDE 18

◮ ν = e−ikx+k2t (qdx + (qx + ikq)dt) is closed ◮ q(x, t) and qx(x, t) decay to 0 at ∞ sufficiently fast. ◮ arg k ∈

4 , 7π 4

x=0

e−ikxq(x, 0)dx = ∞

t=0

ek2t(qx(0, t) + ikq(0, t))dt (GR)

slide-19
SLIDE 19

◮ ν = e−ikx+k2t (qdx + (qx + ikq)dt) is closed ◮ q(x, t) and qx(x, t) decay to 0 at ∞ sufficiently fast. ◮ arg k ∈

4 , 7π 4

x=0

e−ikxq(x, 0)dx = ∞

t=0

ek2t(qx(0, t) + ikq(0, t))dt (GR)

◮ Let us assume first that q(x, 0), qx(0, t), q(0, t) are known.

Reconstruction of q in R+ × R+ is equivalent to the reconstruction of µ since d(µe−ikx+k2t) = ν.

slide-20
SLIDE 20

◮ ν = e−ikx+k2t (qdx + (qx + ikq)dt) is closed ◮ q(x, t) and qx(x, t) decay to 0 at ∞ sufficiently fast. ◮ arg k ∈

4 , 7π 4

x=0

e−ikxq(x, 0)dx = ∞

t=0

ek2t(qx(0, t) + ikq(0, t))dt (GR)

◮ Let us assume first that q(x, 0), qx(0, t), q(0, t) are known.

Reconstruction of q in R+ × R+ is equivalent to the reconstruction of µ since d(µe−ikx+k2t) = ν.

◮ Integration of ν between points of the boundary and (x, t).

slide-21
SLIDE 21

µ1(x, t, k) = eikx−k2t

  • γ1

ν

slide-22
SLIDE 22

µ1(x, t, k) = eikx−k2t

  • γ1

ν µ2(x, t, k) = eikx−k2t

  • γ2

ν

slide-23
SLIDE 23

µ46(x, t, k) = eikx−k2t

  • γ1

ν µ5(x, t, k) = eikx−k2t

  • γ2

ν

slide-24
SLIDE 24

µ46(x, t, k) = eikx−k2t

  • γ46

ν µ5(x, t, k) = eikx−k2t

  • γ5

ν µ123(x, t, k) = eikx−k2t

  • γ123

ν

slide-25
SLIDE 25

Riemann-Hilbert Problems

slide-26
SLIDE 26

Riemann-Hilbert Problems

◮ ◮ µ(x, t, k) = O

1 k

slide-27
SLIDE 27

Riemann-Hilbert Problems

◮ ◮ µ(x, t, k) = O

1 k

  • ◮ µ+ − µ− = φ
slide-28
SLIDE 28

Riemann-Hilbert Problems

◮ ◮ µ(x, t, k) = O

1 k

  • ◮ µ+ − µ− = φ

◮ µ =

1 2πi

  • L

φ(k′)dk′ k′ − k (Plemelj Formula)

slide-29
SLIDE 29

Lax Pairs and closed differential form for the GASP equation.

◮ Lα(u) = ∆u + α

x ∂u ∂x = 0 in H = {z, Re z > 0}.

slide-30
SLIDE 30

Lax Pairs and closed differential form for the GASP equation.

◮ Lα(u) = ∆u + α

x ∂u ∂x = 0 in H = {z, Re z > 0}.

◮ uz ¯ z +

α 2(z + ¯ z)(uz + u¯

z) = 0.

(GASP)

slide-31
SLIDE 31

Lax Pairs and closed differential form for the GASP equation.

◮ Lα(u) = ∆u + α

x ∂u ∂x = 0 in H = {z, Re z > 0}.

◮ uz ¯ z +

α 2(z + ¯ z)(uz + u¯

z) = 0.

(GASP)

   φz(z, k) = (k + ¯ z)α/2(k − z)α/2−1uz(z), φ¯

z(z, k) = (k + ¯

z)α/2−1(k − z)α/2u¯

z(z)

slide-32
SLIDE 32

Lax Pairs and closed differential form for the GASP equation.

◮ Lα(u) = ∆u + α

x ∂u ∂x = 0 in H = {z, Re z > 0}.

◮ uz ¯ z +

α 2(z + ¯ z)(uz + u¯

z) = 0.

(GASP)

   φz(z, k) = (k + ¯ z)α/2(k − z)α/2−1uz(z), φ¯

z(z, k) = (k + ¯

z)α/2−1(k − z)α/2u¯

z(z) ◮ W (z, k) =

  • (k − z)(k + ¯

z) α/2−1 (k + ¯ z)uz(z)dz + (k − z)u¯

z(z)d ¯

z

  • Note that, when α ∈ 2N∗, the differential form has no

singularity in Ω and k may be any complex number. Otherwise, for α ∈ R \ 2N∗, W (z, k) has a pole or a branching point in k or −¯ k if one of this point is in Ω.

slide-33
SLIDE 33

Lax Pairs and closed differential form for the GASP equation.

◮ Lα(u) = ∆u + α

x ∂u ∂x = 0 in H = {z, Re z > 0}.

◮ uz ¯ z +

α 2(z + ¯ z)(uz + u¯

z) = 0.

(GASP)

   φz(z, k) = (k + ¯ z)α/2(k − z)α/2−1uz(z), φ¯

z(z, k) = (k + ¯

z)α/2−1(k − z)α/2u¯

z(z) ◮ W (z, k) =

  • (k − z)(k + ¯

z) α/2−1 (k + ¯ z)uz(z)dz + (k − z)u¯

z(z)d ¯

z

  • Note that, when α ∈ 2N∗, the differential form has no

singularity in Ω and k may be any complex number. Otherwise, for α ∈ R \ 2N∗, W (z, k) has a pole or a branching point in k or −¯ k if one of this point is in Ω.

◮ Lα(u) = 0 ⇔ L2−α(xα−1u) = 0.

slide-34
SLIDE 34

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

slide-35
SLIDE 35

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

φ(z, k) =

  • W (z′, k)
slide-36
SLIDE 36

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

φ(z, k) =

  • W (z′, k)

◮ poles of order m in −zr, −z, z, zr if m > 0.

slide-37
SLIDE 37

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

φ(z, k) =

  • W (z′, k)

◮ poles of order m in −zr, −z, z, zr if m > 0. ◮ φ(z, −k) = −φ(z, k)

slide-38
SLIDE 38

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

φ(z, k) =

  • W (z′, k)

◮ poles of order m in −zr, −z, z, zr if m > 0. ◮ φ(z, −k) = −φ(z, k) ◮ Jump on (z, zr) ∪ (−zr, −z) equal to J(k) = −

  • C

W (z′, k)

slide-39
SLIDE 39

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

φ(z, k) =

  • W (z′, k)

◮ poles of order m in −zr, −z, z, zr if m > 0. ◮ φ(z, −k) = −φ(z, k) ◮ Jump on (z, zr) ∪ (−zr, −z) equal to J(k) = −

  • C

W (z′, k)

◮ J has no singularity in z and zr.

slide-40
SLIDE 40

α = −2m, m ∈ N.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

[(k − z)(k + z)]m+1

φ(z, k) =

  • W (z′, k)

◮ poles of order m in −zr, −z, z, zr if m > 0. ◮ φ(z, −k) = −φ(z, k) ◮ Jump on (z, zr) ∪ (−zr, −z) equal to J(k) = −

  • C

W (z′, k)

◮ J has no singularity in z and zr. ◮ φ ∼k→+∞ u(z)−u(zr) k2m−1

slide-41
SLIDE 41

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)
slide-42
SLIDE 42

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)

φ(z, k) ∼k→+∞ u(z) − u(zr) k

slide-43
SLIDE 43

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)

φ(z, k) ∼k→+∞ u(z) − u(zr) k

φ regular in z and −z, polar singularities in zr et −zr.

slide-44
SLIDE 44

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)

φ(z, k) ∼k→+∞ u(z) − u(zr) k

φ regular in z and −z, polar singularities in zr et −zr.

φzr,−zr (z, k) polar part in this point

slide-45
SLIDE 45

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)

φ(z, k) ∼k→+∞ u(z) − u(zr) k

φ regular in z and −z, polar singularities in zr et −zr.

φzr,−zr (z, k) polar part in this point

φ − φzr,−zr analytic outside (z, zr) ∪ (−zr, −z) and has jump equal to J(z, k)

slide-46
SLIDE 46

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)

φ(z, k) ∼k→+∞ u(z) − u(zr) k

φ regular in z and −z, polar singularities in zr et −zr.

φzr,−zr (z, k) polar part in this point

φ − φzr,−zr analytic outside (z, zr) ∪ (−zr, −z) and has jump equal to J(z, k)

◮ Plemelj formula

  • φ(z, k) −

φzr,−zr (z, k) = 1 2πi

  • (−zr,−z)∪(z,zr)
  • J(z, k′)dk′

k′ − k

slide-47
SLIDE 47

◮ Renormalization :

φ(z, k) = ((k − z)(k + z))mφ(z, k)

  • J(z, k) = ((k − z)(k + z))mJ(z, k)

φ(z, k) ∼k→+∞ u(z) − u(zr) k

φ regular in z and −z, polar singularities in zr et −zr.

φzr,−zr (z, k) polar part in this point

φ − φzr,−zr analytic outside (z, zr) ∪ (−zr, −z) and has jump equal to J(z, k)

◮ Plemelj formula

  • φ(z, k) −

φzr,−zr (z, k) = 1 2πi

  • (−zr,−z)∪(z,zr)
  • J(z, k′)dk′

k′ − k

u(z) − u(zr) = 2 Re ar − 1 πIm

  • (z,zr)
  • J(z, k′)dk′
slide-48
SLIDE 48

Computation of the residue ar of φ in zr

◮ uzdz = 1 2(ut + iun)ds,

uzdz = 1

2(ut − iun)ds.

slide-49
SLIDE 49

Computation of the residue ar of φ in zr

◮ uzdz = 1 2(ut + iun)ds,

uzdz = 1

2(ut − iun)ds. ◮ φ(z, k) =

  • W (z′, k) =

[(k − iy′)ut + ix′un)]ds [(k − z′)(k + z′)]m+1

slide-50
SLIDE 50

Computation of the residue ar of φ in zr

◮ uzdz = 1 2(ut + iun)ds,

uzdz = 1

2(ut − iun)ds. ◮ φ(z, k) =

  • W (z′, k) =

[(k − iy′)ut + ix′un)]ds [(k − z′)(k + z′)]m+1

◮ φ(z, k) =

  • (k − z′)−m−1w(z′, k)dz′ where

w(z′, k) = (k + z′)−m−1((k − iy′)ut + ix′un)τ −1(z′)dz′ τ(z′) is the unit tangent vector to C in z′.

slide-51
SLIDE 51

Computation of the residue ar of φ in zr

◮ uzdz = 1 2(ut + iun)ds,

uzdz = 1

2(ut − iun)ds. ◮ φ(z, k) =

  • W (z′, k) =

[(k − iy′)ut + ix′un)]ds [(k − z′)(k + z′)]m+1

◮ φ(z, k) =

  • (k − z′)−m−1w(z′, k)dz′ where

w(z′, k) = (k + z′)−m−1((k − iy′)ut + ix′un)τ −1(z′)dz′ τ(z′) is the unit tangent vector to C in z′.

◮ m integrations by parts give

φzr,−zr , then ar.

slide-52
SLIDE 52

Let u be a solution to the equation ∆u + αx−1∂xu = 0, α = −2m, m ∈ N, in the domain D with smooth tangential derivatives ut and normal derivatives un on the boundary C. u(z) = − 1 πIm

  • (z,zr)
  • (k − z)(k + ¯

z) mJ(z, k)dk + 2Re ar + u(zr), (1) where ar can be explicitly computed in terms of the tangential derivative along C of ut and un, up to the order m − 1, in zr. Function J(z, k) is given by J(z, k) = −

  • C

W (z′, k), where W (z, k) is the differential form W (z, k) =

  • (k − z)(k + ¯

z) −m−1 ((k + ¯ z)uz(z)dz + (k − z)u¯

z(z)d ¯

z) (2) =

  • (k − z)(k + ¯

z) −m−1 ((k − iy)ut(z) + ixun(z)) ds, (3) with z = x + iy and ds the unit length element on C.

slide-53
SLIDE 53

The case α = −1.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

  • (k − z)(k + z)
slide-54
SLIDE 54

The case α = −1.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

  • (k − z)(k + z)

◮ Riemann surface Sz : two copies of C, Sz,1 et Sz,2 glued

together along the branching cut [−z, z].

slide-55
SLIDE 55

The case α = −1.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

  • (k − z)(k + z)

◮ Riemann surface Sz : two copies of C, Sz,1 et Sz,2 glued

together along the branching cut [−z, z].

slide-56
SLIDE 56

The case α = −1.

◮ W (z, k) = (k + z)uzdz + (k − z)uzdz

  • (k − z)(k + z)

◮ Riemann surface Sz : two copies of C, Sz,1 et Sz,2 glued

together along the branching cut [−z, z].

◮ ◮ λ1(z, k) ∼ k when k → ∞1 on the sheet above Sz,1

λ2(z, k) ∼ −k when k → ∞2 on the sheet below Sz,2

slide-57
SLIDE 57

Integration of the form.

◮ φ(z, k) =

(k + z′)uzdz′ + (k − z′)uzdz′ λ(z′, k)

slide-58
SLIDE 58

Integration of the form.

◮ φ(z, k) =

(k + z′)uzdz′ + (k − z′)uzdz′ λ(z′, k)

◮ Sheet 1

slide-59
SLIDE 59

Integration of the form.

◮ φ(z, k) =

(k + z′)uzdz′ + (k − z′)uzdz′ λ(z′, k)

◮ Sheet 1 ◮ Sheet 2

slide-60
SLIDE 60

Integration of the form.

◮ φ(z, k) =

(k + z′)uzdz′ + (k − z′)uzdz′ λ(z′, k)

◮ Sheet 1 ◮ Sheet 2 ◮ φ(z, ∞1) = −φ(z, ∞2).

slide-61
SLIDE 61

slide-62
SLIDE 62

◮ ◮

φ(z, k) = 1 4iπ

  • Cup∪−Cup

J(z, k′) λ(z, k) λ1(z, k′) + 1

  • dk′

k′ − k + 1 4iπ

  • Clow∪−Clow

J(z, k′) λ(z, k) λ2(z, k′) + 1

  • dk′

k′ − k , (4)

slide-63
SLIDE 63

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

slide-64
SLIDE 64

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

slide-65
SLIDE 65

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

◮ In full generality, no

slide-66
SLIDE 66

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

◮ In full generality, no ◮ Example : Ω = D, W = uz k−z dz + u¯

z

k+z dz. ∆u = 0.

slide-67
SLIDE 67

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

◮ In full generality, no ◮ Example : Ω = D, W = uz k−z dz + u¯

z

k+z dz. ∆u = 0. ◮ uzdz = 1 2(ut + iun)ds, u¯ z = 1 2(ut − iun)ds.

slide-68
SLIDE 68

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

◮ In full generality, no ◮ Example : Ω = D, W = uz k−z dz + u¯

z

k+z dz. ∆u = 0. ◮ uzdz = 1 2(ut + iun)ds, u¯ z = 1 2(ut − iun)ds. ◮

  • ∂D

un(z)

  • 1

k − z − 1 k + z

  • ds(z) =

i

  • ∂D

ut(z)

  • 1

k − z + 1 k + z

  • ds(z)
slide-69
SLIDE 69

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

◮ In full generality, no ◮ Example : Ω = D, W = uz k−z dz + u¯

z

k+z dz. ∆u = 0. ◮ uzdz = 1 2(ut + iun)ds, u¯ z = 1 2(ut − iun)ds. ◮

  • ∂D

un(z)

  • 1

k − z − 1 k + z

  • ds(z) =

i

  • ∂D

ut(z)

  • 1

k − z + 1 k + z

  • ds(z)

◮ ut ≡ 0 et

1 k − z =

  • n=1

zn−1 kn , 1 ¯ z + k =

  • n=1

(−z)n−1 kn

slide-70
SLIDE 70

The notion of good Lax pair.

  • ∂Ω

W (ut, un, k)ds = 0 for all k ∈ K ⊂ C.

◮ Is this relation useful ? In other words, if u = 0 ⇒ un = 0 on

∂Ω ? If u is known on ∂Ω, is it possible to obtain un on ∂Ω ?

◮ In full generality, no ◮ Example : Ω = D, W = uz k−z dz + u¯

z

k+z dz. ∆u = 0. ◮ uzdz = 1 2(ut + iun)ds, u¯ z = 1 2(ut − iun)ds. ◮

  • ∂D

un(z)

  • 1

k − z − 1 k + z

  • ds(z) =

i

  • ∂D

ut(z)

  • 1

k − z + 1 k + z

  • ds(z)

◮ ut ≡ 0 et

1 k − z =

  • n=1

zn−1 kn , 1 ¯ z + k =

  • n=1

(−z)n−1 kn

un(z) =

  • n=−∞

bnzn,

+∞

  • n=1

bn−1 + (−1)nbn−1 kn = 0.

slide-71
SLIDE 71

Good Lax pair for α = −2(m − 1).

◮ ut ≡ 0.

slide-72
SLIDE 72

Good Lax pair for α = −2(m − 1).

◮ ut ≡ 0. ◮

  • ∂D(a,1)

xun(z)ds(z) [(k − z)(k + z)]m = 0, ∀k ∈ C \ (Da ∪ D−a)

slide-73
SLIDE 73

Good Lax pair for α = −2(m − 1).

◮ ut ≡ 0. ◮

  • ∂D(a,1)

xun(z)ds(z) [(k − z)(k + z)]m = 0, ∀k ∈ C \ (Da ∪ D−a)

◮ Da → D. z − a = 1 z−a.

slide-74
SLIDE 74

Good Lax pair for α = −2(m − 1).

◮ ut ≡ 0. ◮

  • ∂D(a,1)

xun(z)ds(z) [(k − z)(k + z)]m = 0, ∀k ∈ C \ (Da ∪ D−a)

◮ Da → D. z − a = 1 z−a. ◮

  • T

zm−1f (z)dz (z − (k − a))m z +

1 k+a

m = 0, f (z) = (x+a)un(z+a).

slide-75
SLIDE 75

Good Lax pair for α = −2(m − 1).

◮ ut ≡ 0. ◮

  • ∂D(a,1)

xun(z)ds(z) [(k − z)(k + z)]m = 0, ∀k ∈ C \ (Da ∪ D−a)

◮ Da → D. z − a = 1 z−a. ◮

  • T

zm−1f (z)dz (z − (k − a))m z +

1 k+a

m = 0, f (z) = (x+a)un(z+a).

◮ f is real valued, f (z) = g(z) + g(1/z), g ∈ H(D),

g(1/z) ∈ H(C \ D).

slide-76
SLIDE 76

  • T

zm−1(g(z) + g(1/z))dz (z − (k − a))m z +

1 k+a

m = 0.

slide-77
SLIDE 77

  • T

zm−1(g(z) + g(1/z))dz (z − (k − a))m z +

1 k+a

m = 0.

−1 k + a m zm−1g(z) (z − (k − a))m (m−1) −1 k + a

  • +
  • 1

k − a m zm−1¯ g(z) (z + (k + a))m (m−1) 1 k − a

  • = 0.
slide-78
SLIDE 78

  • T

zm−1(g(z) + g(1/z))dz (z − (k − a))m z +

1 k+a

m = 0.

−1 k + a m zm−1g(z) (z − (k − a))m (m−1) −1 k + a

  • +
  • 1

k − a m zm−1¯ g(z) (z + (k + a))m (m−1) 1 k − a

  • = 0.

Φ(z) =

m−1

  • p=0

(−1)m−1+pzm−1−p (2m − p − 2)! (m − p − 1)!p!× ×(zm−1g(z))(p)(z2 + 2az + 1)p,

slide-79
SLIDE 79

ϕ(z) = − z 1 + 2az ϕ

1 k + a

  • =

1 k − a.

slide-80
SLIDE 80

ϕ(z) = − z 1 + 2az ϕ

1 k + a

  • =

1 k − a.

Φ(z) + ¯ Φ(ϕ(z)) = 0, z ∈ A.

slide-81
SLIDE 81

ϕ(z) = − z 1 + 2az ϕ

1 k + a

  • =

1 k − a.

Φ(z) + ¯ Φ(ϕ(z)) = 0, z ∈ A.

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) αp = (−1)m−1+p (2m − p − 2)! (m − p − 1)!p! h(z) = zm−1g(z).

slide-82
SLIDE 82

ϕ(z) = − z 1 + 2az ϕ

1 k + a

  • =

1 k − a.

Φ(z) + ¯ Φ(ϕ(z)) = 0, z ∈ A.

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) αp = (−1)m−1+p (2m − p − 2)! (m − p − 1)!p! h(z) = zm−1g(z).

slide-83
SLIDE 83

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) = P2m−2

slide-84
SLIDE 84

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) = P2m−2

◮ h(z) = +∞ n=0 an(z − z1)n.

slide-85
SLIDE 85

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) = P2m−2

◮ h(z) = +∞ n=0 an(z − z1)n. ◮ m−1 ℓ=0 an−ℓβn ℓ = 0, n ≥ 2m − 1.

slide-86
SLIDE 86

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) = P2m−2

◮ h(z) = +∞ n=0 an(z − z1)n. ◮ m−1 ℓ=0 an−ℓβn ℓ = 0, n ≥ 2m − 1. ◮ For N ∈ {0, . . . , m − 2} and n = 2m − 1 + N,

βn

N+1 = · · · = βn m−1 = 0

and βn

0 = 0, . . . , βn N = 0.

slide-87
SLIDE 87

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) = P2m−2

◮ h(z) = +∞ n=0 an(z − z1)n. ◮ m−1 ℓ=0 an−ℓβn ℓ = 0, n ≥ 2m − 1. ◮ For N ∈ {0, . . . , m − 2} and n = 2m − 1 + N,

βn

N+1 = · · · = βn m−1 = 0

and βn

0 = 0, . . . , βn N = 0. ◮ h is a polynomial of degree less than 2m − 2.

slide-88
SLIDE 88

Φ(z) =

m−1

  • p=0

αpzm−1−p(z2 + 2az + 1)ph(p)(z) = P2m−2

◮ h(z) = +∞ n=0 an(z − z1)n. ◮ m−1 ℓ=0 an−ℓβn ℓ = 0, n ≥ 2m − 1. ◮ For N ∈ {0, . . . , m − 2} and n = 2m − 1 + N,

βn

N+1 = · · · = βn m−1 = 0

and βn

0 = 0, . . . , βn N = 0. ◮ h is a polynomial of degree less than 2m − 2. ◮ g is a polynomial de degree less than m − 1.

slide-89
SLIDE 89

  • T

zm−1f (z) (z − (k − a))m(z + (k + a)−1)m = 0, ∀k ∈ C\{Da∪−Da}, where f (z) = g(z) + g 1 z

  • .
slide-90
SLIDE 90

  • T

zm−1f (z) (z − (k − a))m(z + (k + a)−1)m = 0, ∀k ∈ C\{Da∪−Da}, where f (z) = g(z) + g 1 z

  • .

◮ We put ξ = −(k + a)−1, we get

  • T

zm−1f (z) (ξz + 2aξ + 1)m (z − ξ)m = 0 ∀ξ ∈ D \ D

2a 4a2 − 1, 1 4a2 − 1

  • .
slide-91
SLIDE 91

∂m−1 ∂zm−1

  • zm−1f (z)

(ξz + 2aξ + 1)m

  • |z=ξ

= 0.

slide-92
SLIDE 92

∂m−1 ∂zm−1

  • zm−1f (z)

(ξz + 2aξ + 1)m

  • |z=ξ

= 0.

◮ Let us denote F(z) = zm−1f (z) ∈ C2m−2[z]. m−1

  • k=0

(−1)k(m + k − 1)! k!(m − 1 − k)! ξk (ξ2 + 2aξ + 1)m+k F (m−1−k)(ξ) = 0.

slide-93
SLIDE 93

∂m−1 ∂zm−1

  • zm−1f (z)

(ξz + 2aξ + 1)m

  • |z=ξ

= 0.

◮ Let us denote F(z) = zm−1f (z) ∈ C2m−2[z]. m−1

  • k=0

(−1)k(m + k − 1)! k!(m − 1 − k)! ξk (ξ2 + 2aξ + 1)m+k F (m−1−k)(ξ) = 0.

◮ We can extend F to C.

slide-94
SLIDE 94

∂m−1 ∂zm−1

  • zm−1f (z)

(ξz + 2aξ + 1)m

  • |z=ξ

= 0.

◮ Let us denote F(z) = zm−1f (z) ∈ C2m−2[z]. m−1

  • k=0

(−1)k(m + k − 1)! k!(m − 1 − k)! ξk (ξ2 + 2aξ + 1)m+k F (m−1−k)(ξ) = 0.

◮ We can extend F to C. ◮ Denotes z1 = −a +

√ a2 − 1 and z2 = −a − √ a2 − 1 ∀k ∈ {0, 1, . . . , m − 1}, F (k)(z1) = 0, F (k)(z2) = 0.

slide-95
SLIDE 95

∂m−1 ∂zm−1

  • zm−1f (z)

(ξz + 2aξ + 1)m

  • |z=ξ

= 0.

◮ Let us denote F(z) = zm−1f (z) ∈ C2m−2[z]. m−1

  • k=0

(−1)k(m + k − 1)! k!(m − 1 − k)! ξk (ξ2 + 2aξ + 1)m+k F (m−1−k)(ξ) = 0.

◮ We can extend F to C. ◮ Denotes z1 = −a +

√ a2 − 1 and z2 = −a − √ a2 − 1 ∀k ∈ {0, 1, . . . , m − 1}, F (k)(z1) = 0, F (k)(z2) = 0.

◮ F ≡ 0.

slide-96
SLIDE 96

∂m−1 ∂zm−1

  • zm−1f (z)

(ξz + 2aξ + 1)m

  • |z=ξ

= 0.

◮ Let us denote F(z) = zm−1f (z) ∈ C2m−2[z]. m−1

  • k=0

(−1)k(m + k − 1)! k!(m − 1 − k)! ξk (ξ2 + 2aξ + 1)m+k F (m−1−k)(ξ) = 0.

◮ We can extend F to C. ◮ Denotes z1 = −a +

√ a2 − 1 and z2 = −a − √ a2 − 1 ∀k ∈ {0, 1, . . . , m − 1}, F (k)(z1) = 0, F (k)(z2) = 0.

◮ F ≡ 0. ◮ Thank you very much for your attention !