Multi-d shock waves and surface waves

Size: px
Start display at page:

Download "Multi-d shock waves and surface waves"

Transcription

1 Multi-d shock waves and surface waves S. Benzoni-Gavage University of Lyon (Université Claude Bernard Lyon 1 / Institut Camille Jordan) HYP2008 conference, June 11, / 26

2 Outline Theory Eamples 1 Theory Eamples / 26

3 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. Φ= / 26

4 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. Φ= 0 1 Basic assumption: hyperbolicity in t-direction, i.e. for all u U R n, the matri A 0 (u) := df 0 (u) is nonsingular, and for all ν R d, the matri A 0 (u) 1 A j (u) ν j only has real semisimple eigenvalues. 3 / 26

5 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. Φ= 0 1 A 0 (u, ν) := A 0 (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) 3 / 26

6 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. Φ= 0 1 A 0 (u, ν) := A 0 (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) Shock is noncharacteristic iff both matrices A 0 (u ±, Φ) are nonsingular along Φ = 0. 3 / 26

7 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. Φ= 0 1 A 0 (u, ν) := A 0 (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) Shock is classical (or Laian) iff dime u (A 0 (u, Φ)) + dime s (A 0 (u +, Φ)) = n + 1, dime u (A 0 (u +, Φ) + dime s (A 0 (u, Φ)) = n 1. 3 / 26

8 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. Φ= 0 1 A 0 (u, ν) := A 0 (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) Shock is undercompressive iff dime u (A 0 (u, Φ)) + dime s (A 0 (u +, Φ)) = n + 1 p, dime u (A 0 (u +, Φ) + dime s (A 0 (u, Φ)) = n 1 + p. 3 / 26

9 Theory Eamples General equations for a shock wave t t f 0 (u) + j f j (u) = 0 n, Φ(t, ) 0, Φ= 0 d [f 0 (u)] t Φ + [f j (u)] j Φ = 0 n, Φ(t, ) = 0. [g 0 (u)] t Φ + [g j (u)] j Φ = 0 p, Φ(t, ) = 0. 1 A 0 (u, ν) := A 0 (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) Shock is undercompressive iff dime u (A 0 (u, Φ)) + dime s (A 0 (u +, Φ)) = n + 1 p, dime u (A 0 (u +, Φ) + dime s (A 0 (u, Φ)) = n 1 + p. 3 / 26

10 Linear analysis Theory Eamples Ref. planar shock t = σ t d d 1 [Lopatinskĭı 70], [Kreiss 70], [Blokhin 82], [Majda 83], [Freistühler 98]. 4 / 26

11 Linear analysis Theory Eamples Ref. planar shock t Change frame = σ = 0. d 1 [Lopatinskĭı 70], [Kreiss 70], [Blokhin 82], [Majda 83], [Freistühler 98]. 4 / 26

12 Linear analysis Theory Eamples Ref. planar shock t d Change frame = σ = 0. Change coordinates (t, ) (t, y) := (t, 1,..., d 1, Φ(t, )), Φ(t, ) = d χ(t, 1,..., d 1 ). 1 [Lopatinskĭı 70], [Kreiss 70], [Blokhin 82], [Majda 83], [Freistühler 98]. 4 / 26

13 Linear analysis Theory Eamples Ref. planar shock t d Change frame = σ = 0. Change coordinates (t, ) (t, y) := (t, 1,..., d 1, Φ(t, )), Φ(t, ) = d χ(t, 1,..., d 1 ). Linearize eqns about (u, χ) = (u, 0), u := u ±, y d 0. 1 [Lopatinskĭı 70], [Kreiss 70], [Blokhin 82], [Majda 83], [Freistühler 98]. 4 / 26

14 Normal modes analysis Theory Eamples Transmission problem: { A0 (u) t u + A j (u) j u = 0 n, y d 0, [F 0 (u)] t χ + [F j (u)] j χ = [df d (u) u], y d = 0. 5 / 26

15 Normal modes analysis Theory Eamples Transmission problem: { A0 (u) t u + A j (u) j u = 0 n, y d 0, [F 0 (u)] t χ + [F j (u)] j χ = [df d (u) u], y d = 0. Fourier-Laplace transform (t, ˇy) (τ, ˇη) shooting ODE problem. 5 / 26

16 Normal modes analysis Theory Eamples Transmission problem: { A0 (u) t u + A j (u) j u = 0 n, y d 0, [F 0 (u)] t χ + [F j (u)] j χ = [df d (u) u], y d = 0. Fourier-Laplace transform (t, ˇy) (τ, ˇη) shooting ODE problem. A d (u, ν) := A d (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) 5 / 26

17 Normal modes analysis Theory Eamples Transmission problem: { A0 (u) t u + A j (u) j u = 0 n, y d 0, [F 0 (u)] t χ + [F j (u)] j χ = [df d (u) u], y d = 0. Fourier-Laplace transform (t, ˇy) (τ, ˇη) shooting ODE problem. A d (u, ν) := A d (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) Normal modes: χ = X e τt+iη j y j, u = U(y d ) e τt+iη j y j with U L 2 (R), Re (τ) > 0, U(0+) E u (A d (u, τ, i ˇη)) and U(0 ) E s (A d (u, τ, i ˇη)). 5 / 26

18 Normal modes analysis Theory Eamples Transmission problem: { A0 (u) t u + A j (u) j u = 0 n, y d 0, [F 0 (u)] t χ + [F j (u)] j χ = [df d (u) u], y d = 0. Fourier-Laplace transform (t, ˇy) (τ, ˇη) shooting ODE problem. A d (u, ν) := A d (u) 1 (A 0 (u)ν 0 + A j (u) ν j ) Neutral modes of finite energy, or surface waves: χ = X e iη 0t+iη j y j, u = U(y d ) e iη 0t+iη j y j with still U L 2 (R), U(0+) E u (A d (u, iη 0, i ˇη)) and U(0 ) E s (A d (u, iη 0, i ˇη)). 5 / 26

19 Surface waves Theory Eamples t d 1 u d 1 6 / 26

20 Surface waves Theory Eamples t Isotropic elasticity tt u = λ u + (λ + µ) divu, 2 > 0, 2 u u 2 = 0, 2 = 0, d µ 1 u 1 + (2λ + µ) 2 u 2 = 0, 2 = 0. 1 u d [Rayleigh1885] (see also [Serre 06]) 1 6 / 26

21 Surface waves Theory Eamples t Isotropic elasticity tt u = λ u + (λ + µ) divu, 2 > 0, 2 u u 2 = 0, 2 = 0, d µ 1 u 1 + (2λ + µ) 2 u 2 = 0, 2 = 0. 1 u d For λ > 0, λ + µ > 0, Rayleigh waves, or Surface Acoustic Waves, of speed less than λ. [Rayleigh1885] (see also [Serre 06]) 1 6 / 26

22 Surface waves Theory Eamples t Classical shocks in gas dynamics [Bethe 42], [D yakov 54], [Iordanskĭı 57], [Kontorovič 58], [Erpenbeck 62], [Majda 83], [Blokhin 82]. d 1 u d 1 [Menikoff Plohr 89], [SBG Serre 07]. [Jenssen-Lyng 04], 6 / 26

23 Surface waves Theory Eamples t Classical shocks in gas dynamics [Bethe 42], [D yakov 54], [Iordanskĭı 57], [Kontorovič 58], [Erpenbeck 62], [Majda 83], [Blokhin 82]. 1 u d d There eist neutral modes iff 1 M < k 1 + M 2 (r 1), where M = Mach number behind the shock, r = v p /v b with v p,b = volume past/behind the shock, k = 2 + M 2 (v b v p) T p s. 1 [Menikoff Plohr 89], [SBG Serre 07]. [Jenssen-Lyng 04], 6 / 26

24 Surface waves Theory Eamples t Phase boundaries [SBG 98-99], [SBG Freistühler 04] d 1 u d 1 6 / 26

25 Surface waves Theory Eamples t Phase boundaries [SBG 98-99], [SBG Freistühler 04] 1 u d For nondissipative subsonic phase boundaries there eist surface waves, of speed less than u b u p. d 1 6 / 26

26 Outline 1 Theory Eamples / 26

27 Constant coefficients linear problem Interior operator L(u) := A 0 (u) t + A j (u) j Boundary operator B(u) := [F 0 (u)] t + [F j (u)] j [df d (u) ] 8 / 26

28 Constant coefficients linear problem Interior operator L(u) := A 0 (u) t + A j (u) j Boundary operator B(u) := [F 0 (u)] t + [F j (u)] j [df d (u) ] Maimal a priori estimate γ e γt u 2 L 2 + e γt u yd =0 2 L 2 + e γt χ 2 H 1 γ 1 γ e γt L(u)u 2 L 2 + e γt B(u)(χ, u) 2 L 2 8 / 26

29 Constant coefficients linear problem Interior operator L(u) := A 0 (u) t + A j (u) j Boundary operator B(u) := [F 0 (u)] t + [F j (u)] j [df d (u) ] Maimal a priori estimate γ e γt u 2 L 2 + e γt u yd =0 2 L 2 + e γt χ 2 H 1 γ 1 γ e γt L(u)u 2 L 2 + e γt B(u)(χ, u) 2 L 2 OK under uniform Kreiss-Lopatinskĭı condition, i.e. without neutral modes. (Proof based on Kreiss symmetrizers technique.) 8 / 26

30 Constant coefficients linear problem Interior operator L(u) := A 0 (u) t + A j (u) j Boundary operator B(u) := [F 0 (u)] t + [F j (u)] j [df d (u) ] A priori estimate with loss of derivatives γ e γt u 2 L 2 + e γt u yd =0 2 L 2 + e γt χ 2 H 1 γ 1 γ 3 e γt L(u)u 2 L 2 (R + ;H 1 γ ) + 1 γ 2 e γt B(u)(χ, u) 2 H 1 γ 9 / 26

31 Constant coefficients linear problem Interior operator L(u) := A 0 (u) t + A j (u) j Boundary operator B(u) := [F 0 (u)] t + [F j (u)] j [df d (u) ] A priori estimate with loss of derivatives γ e γt u 2 L 2 + e γt u yd =0 2 L 2 + e γt χ 2 H 1 γ 1 γ 3 e γt L(u)u 2 L 2 (R + ;H 1 γ ) + 1 γ 2 e γt B(u)(χ, u) 2 H 1 γ Takes into account neutral modes. (Proof still based on Kreiss symmetrizers technique [Coulombel 02], [Sablé-Tougeron 88].) 9 / 26

32 Fully nonlinear problem Local-in-time eistence of smooth solutions 10 / 26

33 Fully nonlinear problem Local-in-time eistence of smooth solutions under uniform Kreiss Lopatinskĭı condition [Majda 83], [Blokhin 82], [Métivier et al ], 10 / 26

34 Fully nonlinear problem Local-in-time eistence of smooth solutions under uniform Kreiss Lopatinskĭı condition [Majda 83], [Blokhin 82], [Métivier et al ], under mere Kreiss Lopatinskĭı condition [Coulombel Secchi 08]: with neutral modes and characteristic modes ; application to subsonic phase boundaries and compressible 2d-vorte sheets. (Proof using Nash Moser iteration scheme.) 10 / 26

35 Outline 1 Theory Eamples / 26

36 Fully nonlinear problem { A0 (u) t u + A j (u) j u + A d (u, χ) d u = 0 n, y d 0, [F 0 (u)] t χ + [F j (u)] j χ = [F d (u)], y d = 0. A d (u, χ) := A d (u) A 0 (u) t χ A j (u) j χ 12 / 26

37 Fully nonlinear problem { A0 (u) t u + A j (u) j u + A d (u, χ) d u = 0 n, y d 0, J χ + h(u) = 0 n+p, y d = 0. A d (u, χ) := A d (u) A 0 (u) t χ A j (u) j χ 1... J = / 26

38 Fully nonlinear problem { A0 (u) t u + A j (u) j u + A d (u, χ) d u = 0 n, y d 0, J χ + h(u) = 0 n+p, y d = 0. Asymptotic epansion u = u + ε u(η 0 t + η j y j, y d, εt) + ε 2 ü(η 0 t + η j y j, y d, εt) + h.o.t. χ = ε χ(η 0 t + η j y j, y d, εt) + ε 2 χ(η 0 t + η j y j, εt) + h.o.t. [SBG-Rosini 08], [Hunter 89], [Parker 88]. 12 / 26

39 Approimate problems ξ := η 0 t + η j y j, z := y d, τ := εt. First order { (η0 A 0 (u) + η j A j (u)) ξ u + A d (u) z u = 0 n, z 0, Jη ξ χ + dh(u) u = 0 n+p, z = 0, 13 / 26

40 Approimate problems ξ := η 0 t + η j y j, z := y d, τ := εt. First order { (η0 A 0 (u) + η j A j (u)) ξ u + A d (u) z u = 0 n, z 0, Jη ξ χ + dh(u) u = 0 n+p, z = 0, Second order { (η0 A 0 (u) + η j A j (u)) ξ ü + A d (u) z ü = Ṁ, z 0, Jη ξ χ + dh(u) ü = Ġ, z = 0, 13 / 26

41 Approimate problems ξ := η 0 t + η j y j, z := y d, τ := εt. First order { (η0 A 0 (u) + η j A j (u)) ξ u + A d (u) z u = 0 n, z 0, Jη ξ χ + dh(u) u = 0 n+p, z = 0, Second order { (η0 A 0 (u) + η j A j (u)) ξ ü + A d (u) z ü = Ṁ, z 0, Jη ξ χ + dh(u) ü = Ġ, z = 0, Ṁ := A 0 (u) τ u + (η 0 da 0 (u) + η j da j (u)) u ξ u + da d (u) u z u ( ξ χ)(η 0 A 0 (u) + η j A j (u)) z u 13 / 26

42 Approimate problems ξ := η 0 t + η j y j, z := y d, τ := εt. First order { (η0 A 0 (u) + η j A j (u)) ξ u + A d (u) z u = 0 n, z 0, Jη ξ χ + dh(u) u = 0 n+p, z = 0, Second order { (η0 A 0 (u) + η j A j (u)) ξ ü + A d (u) z ü = Ṁ, z 0, Jη ξ χ + dh(u) ü = Ġ, z = 0, Ġ := ( τ χ)e d2 h(u) ( u, u). 13 / 26

43 Transformation of approimate problems Fourier transform ξ k, First order { ik (η0 A 0 (u) + η j A j (u)) u + A d (u) z u = 0 n, z 0, ik Jη χ + dh(u) u = 0 n+p, z = 0, Second order { ik (η0 A 0 (u) + η j A j (u)) ü + A d (u) z ü = Ṁ, z 0, ik Jη χ + dh(u) ü = Ġ, z = / 26

44 Transformation of approimate problems Fourier transform ξ k, elimination of χ, χ. First order { ik (η0 A 0 (u) + η j A j (u)) u + A d (u) z u = 0 n, z 0, C(u; η) u = 0 n+p 1, z = 0, Second order { ik (η0 A 0 (u) + η j A j (u)) ü + A d (u) z ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = / 26

45 Transformation of approimate problems Fourier transform ξ k, elimination of χ, χ. First order { L(u; kη) u = 0n, z 0, C(u; η) u = 0 n+p 1, z = 0, Second order { L(u; kη) ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = 0. L(u; η) := iη 0 A 0 (u) + iη j A j (u) + A d (u) z. 14 / 26

46 Resolution of approimate problems First order level Eistence of linear surface wave = L 2 (dz) solution u 1 of { L(u; η) u1 = 0 n, z 0, C(u; η) u 1 = 0 n+p 1 z = / 26

47 Resolution of approimate problems First order level Eistence of linear surface wave = L 2 (dz) solution u 1 of { L(u; η) u1 = 0 n, z 0, C(u; η) u 1 = 0 n+p 1 z = 0. Homogeneity = other square integrable solutions of first order system of the form u(k, z, τ) = W (k, τ) u 1 (kz), k > / 26

48 Resolution of approimate problems First order level Eistence of linear surface wave = L 2 (dz) solution u 1 of { L(u; η) u1 = 0 n, z 0, C(u; η) u 1 = 0 n+p 1 z = 0. Homogeneity = other square integrable solutions of first order system of the form u(k, z, τ) = W (k, τ) u 1 (kz), k > 0. Amplitude function: F 1 (W ) =: w(ξ, τ). 15 / 26

49 Resolution of approimate problems First order level Eistence of linear surface wave = L 2 (dz) solution u 1 of { L(u; η) u1 = 0 n, z 0, C(u; η) u 1 = 0 n+p 1 z = 0. Homogeneity = other square integrable solutions of first order system of the form u(k, z, τ) = W (k, τ) u 1 (kz), k > 0. Amplitude function: F 1 (W ) =: w(ξ, τ). Second order level Solvability condition by means of an adjoint problem. 15 / 26

50 Solvability of second order problem ( ) { L(u; kη) ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = 0, 16 / 26

51 Solvability of second order problem ( ) { L(u; kη) ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = 0, L(u; η) := iη 0 A 0 (u) + iη j A j (u) + A d (u) z, 16 / 26

52 Solvability of second order problem ( ) { L(u; kη) ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = 0, L(u; η) := iη 0 A 0 (u) + iη j A j (u) + A d (u) z, C(u; η)u = C + u(0+) C u(0 ), 16 / 26

53 Solvability of second order problem ( ) { L(u; kη) ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = 0, L(u; η) := iη 0 A 0 (u) + iη j A j (u) + A d (u) z, C(u; η)u = C + u(0+) C u(0 ), ( ) ( Ad (u ) 0 D = 0 A d (u + ) D+ ) N + P ( C C + ). 16 / 26

54 Solvability of second order problem ( ) { L(u; kη) ü = Ṁ, z 0, C(u; η)ü = T (η)ġ, z = 0, L(u; η) := iη 0 A 0 (u) + iη j A j (u) + A d (u) z, C(u; η)u = C + u(0+) C u(0 ), ( ) ( ) Ad (u ) 0 D = 0 A d (u + ) D+ N + P ( C C + ). There eists a L 2 (dz) solution ü of ( ) iff v Ṁ dz + (v(0 ) v(0+) )PT Ġ = 0, with v solution of { L(u; kη) v = 0 n, z 0, D + v(0+) D v(0 ) = 0 n p+1, z = / 26

55 Resulting amplitude equation Nonlocal generalisation of Burgers equation: τ w + ξ Q[w] = 0, F (Q[w])(k) = + Λ(k l, l)ŵ(k l)ŵ(l)dl. with piecewise smooth kernel Λ, homogeneous degree / 26

56 Resulting amplitude equation Nonlocal generalisation of Burgers equation: F (Q[w])(k) = τ w + ξ Q[w] = 0, + Λ(k l, l)ŵ(k l)ŵ(l)dl. with piecewise smooth kernel Λ, homogeneous degree 0. Recover classical inviscid Burgers equation if Λ 1/2 (arises in case of neutral modes of infinite energy [Artola-Majda 87]). 17 / 26

57 Nonlocal Burgers equations In Fourier variables: τ ŵ + ik + Λ(k l, l)ŵ(k l, τ)ŵ(l, τ)dl = 0. Eistence of smooth solutions? Well-posedness? 18 / 26

58 Nonlocal Burgers equations In Fourier variables: τ ŵ + ik + Λ(k l, l)ŵ(k l, τ)ŵ(l, τ)dl = 0. Eistence of smooth solutions? Well-posedness? Properties of Λ: 18 / 26

59 Nonlocal Burgers equations In Fourier variables: τ ŵ + ik + Λ(k l, l)ŵ(k l, τ)ŵ(l, τ)dl = 0. Eistence of smooth solutions? Well-posedness? Properties of Λ: Λ(k, l) = Λ(l, k) 18 / 26

60 Nonlocal Burgers equations In Fourier variables: τ ŵ + ik + Λ(k l, l)ŵ(k l, τ)ŵ(l, τ)dl = 0. Eistence of smooth solutions? Well-posedness? Properties of Λ: Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) 18 / 26

61 Nonlocal Burgers equations In Fourier variables: τ ŵ + ik + Λ(k l, l)ŵ(k l, τ)ŵ(l, τ)dl = 0. Eistence of smooth solutions? Well-posedness? Properties of Λ: Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) Λ(1, 0 ) = Λ(1, 0+) 18 / 26

62 Nonlocal Burgers equations In Fourier variables: τ ŵ + ik + Λ(k l, l)ŵ(k l, τ)ŵ(l, τ)dl = 0. Eistence of smooth solutions? Well-posedness? Properties of Λ: Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) Λ(k + ξ, ξ) = Λ(k, ξ) 18 / 26

63 Hamiltonian nonlocal Burgers equations Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) Λ(k + ξ, ξ) = Λ(k, ξ) = Hamiltonian structure : H[w] := 1 3 τ w + δh[w] = 0, Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dk dl. 19 / 26

64 Hamiltonian nonlocal Burgers equations Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) Λ(k + ξ, ξ) = Λ(k, ξ) = Hamiltonian structure : H[w] := 1 3 τ w + δh[w] = 0, Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dk dl. = Local eistence of smooth periodic solutions [Hunter 06] (also see [Alì Hunter Parker 02]). 19 / 26

65 Stable nonlocal Burgers equations Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) Λ(1, 0 ) = Λ(1, 0+) = a priori estimates, 20 / 26

66 Stable nonlocal Burgers equations Λ(k, l) = Λ(l, k) Λ( k, l) = Λ(k, l) Λ(1, 0 ) = Λ(1, 0+) = a priori estimates, and eventually local H 2 well-posedness [SBG 08]. 20 / 26

67 A priori estimates Local Burgers: d dτ ( n ξ w)2 ξ w L ( n ξ w)2. 21 / 26

68 A priori estimates Local Burgers: d dτ ( n ξ w)2 ξ w L ( n ξ w)2. Nonlocal Burgers: d ( ξ n dτ w)2 F ( ξ w) L 1 ( n ξ w)2. 21 / 26

69 A priori estimates L 2 estimate (n = 0) : d w 2 dξ = d ŵ 2 dk = dτ dτ ( ) 2 Re i k Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk 22 / 26

70 A priori estimates L 2 estimate (n = 0) : d w 2 dξ = d ŵ 2 dk = dτ dτ ( ) 2 Re i k Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk 2 Λ L ŵ 2 L 2 by Fubini and Cauchy-Schwarz! kŵ(k) dk 22 / 26

71 A priori estimates H 1 estimate (n = 1) : d ( ξ w) 2 dξ = d k 2 ŵ 2 dk = dτ dτ ( ) 2 Re i k 3 Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk 23 / 26

72 A priori estimates H 1 estimate (n = 1) : d ( ξ w) 2 dξ = d k 2 ŵ 2 dk = dτ dτ ( ) 2 Re i k 3 Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk = ( 4 Re ) i k 2 (k l) Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk 23 / 26

73 A priori estimates H 1 estimate (n = 1) : d ( ξ w) 2 dξ = d k 2 ŵ 2 dk = dτ dτ ( ) 2 Re i k 3 Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk = ( 4 Re ) i k 2 (k l) Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk. k l... Λ L ξ w L 1 ξ w 2 L / 26

74 A priori estimates H 1 estimate (cont.) ( ) Re i k 2 (k l) Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk = k > l i i k > l k > l k 2 (k l) Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk k 2 (k l) Λ(l k, l)ŵ(l k)ŵ( l)ŵ(k) dldk. 24 / 26

75 A priori estimates H 1 estimate (cont.) ( ) Re i k 2 (k l) Λ(k l, l)ŵ(k l)ŵ(l)ŵ( k) dldk k > l i k > l k(k l) ((k l) Λ(l k, l) kλ(k, l)) ŵ(l k)ŵ( l)ŵ(k) dldk. after change of variables (k, l) (k l, l) in first integral. = 24 / 26

76 Local well-posedness Theorem ([SBG 08]) If Λ is smooth outside the lines k = 0, l = 0, and k + l = 0, homogeneous degree zero, preserves real-valued functions, and satisfies the stability condition Λ(1, 0 ) = Λ( 1, 0 ), then for all w 0 H 2 (R) there eists T > 0 and a unique solution w C (0, T ; H 2 (R)) C 1 (0, T ; H 1 (R)) such that w(0) = w 0 of the nonlocal Burgers equation of kernel Λ, and the mapping is continuous. H 2 (R) C (0, T ; H 2 (R)) w 0 w 25 / 26

77 Blow-up criterion The solution w can be etended beyond T provided that T 0 F ( ξw) L 1 is finite. 26 / 26

78 Blow-up criterion The solution w can be etended beyond T provided that T 0 F ( ξw) L 1 is finite. Applications 26 / 26

79 Blow-up criterion The solution w can be etended beyond T provided that T 0 F ( ξw) L 1 is finite. Applications elasticity: Λ(k + ξ, ξ) = Λ(k, ξ), 26 / 26

80 Blow-up criterion The solution w can be etended beyond T provided that T 0 F ( ξw) L 1 is finite. Applications elasticity: Λ(k + ξ, ξ) = Λ(k, ξ), phase boundaries: Λ(1, 0 ) Λ(1, 0+)! 26 / 26

Weakly nonlinear surface waves and subsonic phase boundaries

Weakly nonlinear surface waves and subsonic phase boundaries Weakly nonlinear surface waves and subsonic phase boundaries S. Benzoni-Gavage & M. D. Rosini May 14, 28 Abstract The aim of this work is twofold. In a first, abstract part, it is shown how to derive an

More information

Nonlinear stability of compressible vortex sheets in two space dimensions

Nonlinear stability of compressible vortex sheets in two space dimensions of compressible vortex sheets in two space dimensions J.-F. Coulombel (Lille) P. Secchi (Brescia) CNRS, and Team SIMPAF of INRIA Futurs Evolution Equations 2006, Mons, August 29th Plan 1 2 3 Related problems

More information

2D compressible vortex sheets. Paolo Secchi

2D compressible vortex sheets. Paolo Secchi 2D compressible vortex sheets Paolo Secchi Department of Mathematics Brescia University Joint work with J.F. Coulombel EVEQ 2008, International Summer School on Evolution Equations, Prague, Czech Republic,

More information

Piecewise Smooth Solutions to the Burgers-Hilbert Equation

Piecewise Smooth Solutions to the Burgers-Hilbert Equation Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang

More information

Hyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan

Hyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global

More information

The Existence of Current-Vortex Sheets in Ideal Compressible Magnetohydrodynamics

The Existence of Current-Vortex Sheets in Ideal Compressible Magnetohydrodynamics ARMA manuscript No. (will be inserted by the editor) The Existence of Current-Vortex Sheets in Ideal Compressible Magnetohydrodynamics Yuri Trahinin Abstract We prove the local-in-time existence of solutions

More information

Local well-posedness of nonlocal Burgers equations

Local well-posedness of nonlocal Burgers equations Local well-poseness of nonlocal Burgers equations Sylvie Benzoni-Gavage To cite this version: Sylvie Benzoni-Gavage. Local well-poseness of nonlocal Burgers equations. Differential Integral Equations,

More information

RECTILINEAR VORTEX SHEETS OF INVISCID LIQUID-GAS TWO-PHASE FLOW: LINEAR STABILITY

RECTILINEAR VORTEX SHEETS OF INVISCID LIQUID-GAS TWO-PHASE FLOW: LINEAR STABILITY RECTILINEAR VORTEX SHEETS OF INVISCID LIQUID-GAS TWO-PHASE FLOW: LINEAR STABILITY LIZHI RUAN, DEHUA WANG, SHANGKUN WENG, AND CHANGJIANG ZHU Abstract. The vortex sheet solutions are considered for the inviscid

More information

REGULARITY OF WEAKLY WELL POSED HYPERBOLIC MIXED PROBLEMS WITH CHARACTERISTIC BOUNDARY

REGULARITY OF WEAKLY WELL POSED HYPERBOLIC MIXED PROBLEMS WITH CHARACTERISTIC BOUNDARY REGULARITY OF WEAKLY WELL POSED HYPERBOLIC MIXED PROBLEMS WITH CHARACTERISTIC BOUNDARY ALESSANDRO MORANDO AND PAOLO SECCHI Abstract. We study the mixed initial-boundary value problem for a linear hyperbolic

More information

Decay profiles of a linear artificial viscosity system

Decay profiles of a linear artificial viscosity system Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter

More information

arxiv: v1 [math.ap] 15 Apr 2016

arxiv: v1 [math.ap] 15 Apr 2016 arxiv:604.04546v [math.ap] 5 Apr 06 Geometric optics for surface waves in nonlinear elasticity Jean-François Coulombel & Mark Williams December 5, 07 CNRS and Université de Nantes, Laboratoire de mathématiques

More information

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows

Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,

More information

Viscous Boundary Value Problems for Symmetric Systems with Variable Multiplicities

Viscous Boundary Value Problems for Symmetric Systems with Variable Multiplicities Viscous Boundary Value Problems for Symmetric Systems with Variable Multiplicities Olivier Gues, Guy Métivier, Mark Williams, Kevin Zumbrun June 19, 2006 Abstract Extending investigations of Métivier&Zumbrun

More information

On the Dependence of Euler Equations on Physical Parameters

On the Dependence of Euler Equations on Physical Parameters On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,

More information

Corrugation instability and bifurcation of slow viscous MHD shocks in a duct. Kevin Zumbrun

Corrugation instability and bifurcation of slow viscous MHD shocks in a duct. Kevin Zumbrun Corrugation instability and bifurcation of slow viscous MHD shocks in a duct Kevin Collaborators: Blake Barker Rafael Monteiro Department of Mathematics Indiana University Sponsored by NSF Grants no. DMS-0300487

More information

Hyperbolic Conservation Laws Past and Future

Hyperbolic Conservation Laws Past and Future Hyperbolic Conservation Laws Past and Future Barbara Lee Keyfitz Fields Institute and University of Houston bkeyfitz@fields.utoronto.ca Research supported by the US Department of Energy, National Science

More information

Instability of Finite Difference Schemes for Hyperbolic Conservation Laws

Instability of Finite Difference Schemes for Hyperbolic Conservation Laws Instability of Finite Difference Schemes for Hyperbolic Conservation Laws Alberto Bressan ( ), Paolo Baiti ( ) and Helge Kristian Jenssen ( ) ( ) Department of Mathematics, Penn State University, University

More information

Generic Singularities of Solutions to some Nonlinear Wave Equations

Generic Singularities of Solutions to some Nonlinear Wave Equations Generic Singularities of Solutions to some Nonlinear Wave Equations Alberto Bressan Deartment of Mathematics, Penn State University (Oberwolfach, June 2016) Alberto Bressan (Penn State) generic singularities

More information

Rates of Convergence to Self-Similar Solutions of Burgers Equation

Rates of Convergence to Self-Similar Solutions of Burgers Equation Rates of Convergence to Self-Similar Solutions of Burgers Equation by Joel Miller Andrew Bernoff, Advisor Advisor: Committee Member: May 2 Department of Mathematics Abstract Rates of Convergence to Self-Similar

More information

Math 575-Lecture 26. KdV equation. Derivation of KdV

Math 575-Lecture 26. KdV equation. Derivation of KdV Math 575-Lecture 26 KdV equation We look at the KdV equations and the so-called integrable systems. The KdV equation can be written as u t + 3 2 uu x + 1 6 u xxx = 0. The constants 3/2 and 1/6 are not

More information

Dirichlet s principle and well posedness of steady state solutions in peridynamics

Dirichlet s principle and well posedness of steady state solutions in peridynamics Dirichlet s principle and well posedness of steady state solutions in peridynamics Petronela Radu Work supported by NSF - DMS award 0908435 January 19, 2011 The steady state peridynamic model Consider

More information

Wave breaking in the short-pulse equation

Wave breaking in the short-pulse equation Dynamics of PDE Vol.6 No.4 29-3 29 Wave breaking in the short-pulse equation Yue Liu Dmitry Pelinovsky and Anton Sakovich Communicated by Y. Charles Li received May 27 29. Abstract. Sufficient conditions

More information

On critical behaviour in systems of Hamiltonian partial differential equations

On critical behaviour in systems of Hamiltonian partial differential equations On critical behaviour in systems of Hamiltonian partial differential equations arxiv:1311.7166v1 [math-ph] 7 Nov 13 B. Dubrovin, T. Grava, C. Klein, A. Moro Abstract We study the critical behaviour of

More information

On the weakly nonlinear Kelvin-Helmholtz instability of tangential discontinuities in MHD

On the weakly nonlinear Kelvin-Helmholtz instability of tangential discontinuities in MHD On the weakly nonlinear Kelvin-Helmholtz instability of tangential discontinuities in MHD John K. Hunter Department of Mathematics University of California at Davis Davis, California 9566, U.S.A. and J.

More information

A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations

A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction

More information

Statistical Geometry Processing Winter Semester 2011/2012

Statistical Geometry Processing Winter Semester 2011/2012 Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian

More information

Chapter 6 Controllability and Obervability

Chapter 6 Controllability and Obervability Chapter 6 Controllability and Obervability Controllability: whether or not the state-space equation can be controlled from input. Observability: whether or not the initial state can be observed from output.

More information

Given the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w.

Given the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w. Selected problems from the tetbook J. Neustupa, S. Kračmar: Sbírka příkladů z Matematiky I Problems in Mathematics I I. LINEAR ALGEBRA I.. Vectors, vector spaces Given the vectors u, v, w and real numbers

More information

Second-order optimality conditions for state-constrained optimal control problems

Second-order optimality conditions for state-constrained optimal control problems Second-order optimality conditions for state-constrained optimal control problems J. Frédéric Bonnans INRIA-Saclay & CMAP, Ecole Polytechnique, France Joint work with Audrey Hermant, INRIA-Saclay & CMAP

More information

M 340L CS Homework Set 1

M 340L CS Homework Set 1 M 340L CS Homework Set 1 Solve each system in Problems 1 6 by using elementary row operations on the equations or on the augmented matri. Follow the systematic elimination procedure described in Lay, Section

More information

Geometric optics expansions for linear hyperbolic boundary value problems and optimality of energy estimates for surface waves.

Geometric optics expansions for linear hyperbolic boundary value problems and optimality of energy estimates for surface waves. Geometric optics expansions for linear hyperbolic boundary value problems and optimality of energy estimates for surface waves. Antoine Benoit To cite this version: Antoine Benoit. Geometric optics expansions

More information

Shooting methods for numerical solutions of control problems constrained. by linear and nonlinear hyperbolic partial differential equations

Shooting methods for numerical solutions of control problems constrained. by linear and nonlinear hyperbolic partial differential equations Shooting methods for numerical solutions of control problems constrained by linear and nonlinear hyperbolic partial differential equations by Sung-Dae Yang A dissertation submitted to the graduate faculty

More information

2 Introduction T(,t) = = L Figure..: Geometry of the prismatical rod of Eample... T = u(,t) = L T Figure..2: Geometry of the taut string of Eample..2.

2 Introduction T(,t) = = L Figure..: Geometry of the prismatical rod of Eample... T = u(,t) = L T Figure..2: Geometry of the taut string of Eample..2. Chapter Introduction. Problems Leading to Partial Dierential Equations Many problems in science and engineering are modeled and analyzed using systems of partial dierential equations. Realistic models

More information

SELF-SIMILAR SOLUTIONS FOR THE 2-D BURGERS SYSTEM IN INFINITE SUBSONIC CHANNELS

SELF-SIMILAR SOLUTIONS FOR THE 2-D BURGERS SYSTEM IN INFINITE SUBSONIC CHANNELS Bull. Korean Math. oc. 47 010, No. 1, pp. 9 37 DOI 10.4134/BKM.010.47.1.09 ELF-IMILAR OLUTION FOR THE -D BURGER YTEM IN INFINITE UBONIC CHANNEL Kyungwoo ong Abstract. We establish the existence of weak

More information

Applied Mathematics Course Contents Matrices, linear systems Gauss-Jordan elimination, Determinants, Cramer's rule, inverse matrix.

Applied Mathematics Course Contents Matrices, linear systems Gauss-Jordan elimination, Determinants, Cramer's rule, inverse matrix. Applied Mathematics Course Contents Matrices, linear systems Gauss-Jordan elimination, Determinants, Cramer's rule, inverse matri. Vector spaces, inner product and norm, linear transformations. Matri eigenvalue

More information

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

More information

The Scalar Conservation Law

The Scalar Conservation Law The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan

More information

arxiv: v1 [math.ap] 28 May 2009

arxiv: v1 [math.ap] 28 May 2009 Wave breaking in the short-pulse equation arxiv:95.4668v [math.ap] 8 May 9 Yue Liu, Dmitry Pelinovsky, and Anton Sakovich Department of Mathematics, University of Teas at Arlington, Arlington, TX, 769,

More information

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Analytical Solution Techniques J. Kevorkian University of Washington Wadsworth & Brooks/Cole Advanced Books & Software Pacific Grove, California C H A P T E R 1 The Diffusion

More information

Theory and applications of time-frequency analysis

Theory and applications of time-frequency analysis Theory and applications of time-frequency analysis Ville Turunen (ville.turunen@aalto.fi) Abstract: When and how often something happens in a signal? By properly quantizing these questions, we obtain the

More information

Rayleigh Waves in Linear Elasticity as a Propagation of Singularities Phenomenon Michael E. Taylor

Rayleigh Waves in Linear Elasticity as a Propagation of Singularities Phenomenon Michael E. Taylor Rayleigh Waves in Linear Elasticity as a Propagation of Singularities Phenomenon Michael E. Taylor Abstract. We examine the propagation of surface waves known as Rayleigh waves from the perspective of

More information

Propagation and quenching in a reactive Burgers-Boussinesq system.

Propagation and quenching in a reactive Burgers-Boussinesq system. Propagation and quenching in a reactive system., University of New Meico, Los Alamos National Laboratory, University of Chicago October 28 Between Math and Astrophysics Motivation The model Reaction Compression

More information

Semiclassical Limit for the Schrödinger-Poisson Equation in a Crystal

Semiclassical Limit for the Schrödinger-Poisson Equation in a Crystal Semiclassical Limit for the Schrödinger-Poisson Equation in a Crystal PHILIPPE BECHOUCHE Université de Nice Sophia-Antipolis Universidad de Granada NORBERT J. MAUSER Universität Wien Courant Institute

More information

First-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012

First-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012 First-order overdetermined systems for elliptic problems John Strain Mathematics Department UC Berkeley July 2012 1 OVERVIEW Convert elliptic problems to first-order overdetermined form Control error via

More information

Lecture 3: Asymptotic Methods for the Reduced Wave Equation

Lecture 3: Asymptotic Methods for the Reduced Wave Equation Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed

More information

Elliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations

Elliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations 5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have

More information

A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem

A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract

More information

Stability and Shoaling in the Serre Equations. John D. Carter. March 23, Joint work with Rodrigo Cienfuegos.

Stability and Shoaling in the Serre Equations. John D. Carter. March 23, Joint work with Rodrigo Cienfuegos. March 23, 2009 Joint work with Rodrigo Cienfuegos. Outline The Serre equations I. Derivation II. Properties III. Solutions IV. Solution stability V. Wave shoaling Derivation of the Serre Equations Derivation

More information

Modulation and water waves

Modulation and water waves Summer School on Water Waves................Newton Institute August 2014 Part 1 Thomas J. Bridges, University of Surrey Modulation of uniform flows and KdV h 0 u 0 00000000000000000000000000000 11111111111111111111111111111

More information

Decay rates for partially dissipative hyperbolic systems

Decay rates for partially dissipative hyperbolic systems Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods

More information

THERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING

THERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING THERMOACOUSTIC TOMOGRAPHY ARISING IN BRAIN IMAGING PLAMEN STEFANOV AND GUNTHER UHLMANN Abstract. We study the mathematical model of thermoacoustic and photoacoustic tomography when the sound speed has

More information

The Gaussian free field, Gibbs measures and NLS on planar domains

The Gaussian free field, Gibbs measures and NLS on planar domains The Gaussian free field, Gibbs measures and on planar domains N. Burq, joint with L. Thomann (Nantes) and N. Tzvetkov (Cergy) Université Paris Sud, Laboratoire de Mathématiques d Orsay, CNRS UMR 8628 LAGA,

More information

c. Better work with components of slowness vector s (or wave + k z 2 = k 2 = (ωs) 2 = ω 2 /c 2. k=(kx,k z )

c. Better work with components of slowness vector s (or wave + k z 2 = k 2 = (ωs) 2 = ω 2 /c 2. k=(kx,k z ) .50 Introduction to seismology /3/05 sophie michelet Today s class: ) Eikonal equation (basis of ray theory) ) Boundary conditions (Stein s book.3.0) 3) Snell s law Some remarks on what we discussed last

More information

Hyperbolic Systems of Conservation Laws. I - Basic Concepts

Hyperbolic Systems of Conservation Laws. I - Basic Concepts Hyperbolic Systems of Conservation Laws I - Basic Concepts Alberto Bressan Mathematics Department, Penn State University Alberto Bressan (Penn State) Hyperbolic Systems of Conservation Laws 1 / 27 The

More information

A special choice of the kernel amounts to coupling (1) with an elliptic equation: (3) t u + div x (F (u) r) = 0, r + r = u.

A special choice of the kernel amounts to coupling (1) with an elliptic equation: (3) t u + div x (F (u) r) = 0, r + r = u. Waves in Systems of Conservation Laws: Large vs Small Scales Denis Serre Abstract. Several results about the stability of waves in systems of conservation/balance laws, disseminated in the litterature,

More information

INTRODUCTION TO PDEs

INTRODUCTION TO PDEs INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial

More information

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

Stability of Small Viscosity Noncharacteristic Boundary Layers. Guy Métivier MAB Université de Bordeaux I Talence Cedex, France Stability of Small Viscosity Noncharacteristic Boundary Layers Guy Métivier MAB Université de Bordeaux I 33405 Talence Cedex, France June 12, 2003 Contents 1 An Example 13 1.1 The equation...........................

More information

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that

More information

Stability of an abstract wave equation with delay and a Kelvin Voigt damping

Stability of an abstract wave equation with delay and a Kelvin Voigt damping Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability

More information

The Navier Stokes Equations for Incompressible Flows: Solution Properties at Potential Blow Up Times

The Navier Stokes Equations for Incompressible Flows: Solution Properties at Potential Blow Up Times The Navier Stokes Equations for Incompressible Flows: Solution Properties at Potential Blow Up Times Jens Lorenz Department of Mathematics and Statistics, UNM, Albuquerque, NM 873 Paulo Zingano Dept. De

More information

Nonlinear Wave Propagation in 1D Random Media

Nonlinear Wave Propagation in 1D Random Media Nonlinear Wave Propagation in 1D Random Media J. B. Thoo, Yuba College UCI Computational and Applied Math Seminar Winter 22 Typeset by FoilTEX Background O Doherty and Anstey (Geophysical Prospecting,

More information

Some Collision solutions of the rectilinear periodically forced Kepler problem

Some Collision solutions of the rectilinear periodically forced Kepler problem Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University

More information

First Order Equations

First Order Equations 10 1 Linear and Semilinear Equations Chapter First Order Equations Contents 1 Linear and Semilinear Equations 9 Quasilinear Equations 19 3 Wave Equation 6 4 Sstems of Equations 31 1 Linear and Semilinear

More information

MICROLOCAL ANALYSIS METHODS

MICROLOCAL ANALYSIS METHODS MICROLOCAL ANALYSIS METHODS PLAMEN STEFANOV One of the fundamental ideas of classical analysis is a thorough study of functions near a point, i.e., locally. Microlocal analysis, loosely speaking, is analysis

More information

MATH 220 solution to homework 5

MATH 220 solution to homework 5 MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly

More information

New developments for increased performance of the SBP-SAT finite difference technique

New developments for increased performance of the SBP-SAT finite difference technique New developments for increased performance of the SBP-SAT finite difference technique Jan Nordström 1 and Peter Eliasson 2 1 Department of Mathematics, Computational Mathematics, University of Linköping,

More information

Microlocal Analysis : a short introduction

Microlocal Analysis : a short introduction Microlocal Analysis : a short introduction Plamen Stefanov Purdue University Mini Course, Fields Institute, 2012 Plamen Stefanov (Purdue University ) Microlocal Analysis : a short introduction 1 / 25 Introduction

More information

AN ASYMPTOTIC PRESERVING SCHEME BASED ON A NEW FORMULATION FOR NLS IN THE SEMICLASSICAL LIMIT

AN ASYMPTOTIC PRESERVING SCHEME BASED ON A NEW FORMULATION FOR NLS IN THE SEMICLASSICAL LIMIT AN ASYMPTOTIC PRESERVING SCHEME BASED ON A NEW FORMULATION FOR NLS IN THE SEMICLASSICAL LIMIT CHRISTOPHE BESSE, RÉMI CARLES, AND FLORIAN MÉHATS Abstract. We consider the semiclassical limit for the nonlinear

More information

Computation of entropic measure-valued solutions for Euler equations

Computation of entropic measure-valued solutions for Euler equations Computation of entropic measure-valued solutions for Euler equations Eitan Tadmor Center for Scientific Computation and Mathematical Modeling (CSCAMM) Department of Mathematics, Institute for Physical

More information

Chapter 9. Derivatives. Josef Leydold Mathematical Methods WS 2018/19 9 Derivatives 1 / 51. f x. (x 0, f (x 0 ))

Chapter 9. Derivatives. Josef Leydold Mathematical Methods WS 2018/19 9 Derivatives 1 / 51. f x. (x 0, f (x 0 )) Chapter 9 Derivatives Josef Leydold Mathematical Methods WS 208/9 9 Derivatives / 5 Difference Quotient Let f : R R be some function. The the ratio f = f ( 0 + ) f ( 0 ) = f ( 0) 0 is called difference

More information

THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS

THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS Journal of Hyperbolic Differential Equations Vol., No. 4 005 909 917 c World Scientific Publishing Company THE ELLIPTICITY PRINCIPLE FOR SELF-SIMILAR POTENTIAL FLOWS VOLKER ELLING, and TAI-PING LIU, Dept.

More information

Nonlinear Evolution Equations 1

Nonlinear Evolution Equations 1 Nonlinear Evolution Equations 1 John K. Hunter October 1996 1 c John K. Hunter 1996 Contents 1 Introduction 1 1.1 Evolution equations....................... 1 1.2 Blow up..............................

More information

Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary

Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary Weimin Sheng (Joint with Li-Xia Yuan) Zhejiang University IMS, NUS, 8-12 Dec 2014 1 / 50 Outline 1 Prescribing

More information

Gaussian beams for high frequency waves: some recent developments

Gaussian beams for high frequency waves: some recent developments Gaussian beams for high frequency waves: some recent developments Jianliang Qian Department of Mathematics, Michigan State University, Michigan International Symposium on Geophysical Imaging with Localized

More information

R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant

R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant CONTENTS I. Introductory Remarks S1. General Information about the Variety of Solutions.

More information

Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities

Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Margaret Beck Joint work with Anna Ghazaryan, University of Kansas and Björn Sandstede, Brown University September

More information

Model reduction for linear systems by balancing

Model reduction for linear systems by balancing Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,

More information

Transformed Companion Matrices as a Theoretical Tool in the Numerical Analysis of Differential Equations

Transformed Companion Matrices as a Theoretical Tool in the Numerical Analysis of Differential Equations Transformed Companion Matrices as a Theoretical Tool in the Numerical Analysis of Differential Equations Winfried Auzinger Vienna University of Technology (work in progress, with O. Koch and G. Schranz-Kirlinger)

More information

THE INSTABILITY OF BOURGAIN-WANG SOLUTIONS FOR THE L 2 CRITICAL NLS FRANK MERLE, PIERRE RAPHAËL, AND JEREMIE SZEFTEL

THE INSTABILITY OF BOURGAIN-WANG SOLUTIONS FOR THE L 2 CRITICAL NLS FRANK MERLE, PIERRE RAPHAËL, AND JEREMIE SZEFTEL THE INSTABILITY OF BOURGAIN-WANG SOLUTIONS FOR THE L CRITICAL NLS FRANK MERLE, PIERRE RAPHAËL, AND JEREMIE SZEFTEL Abstract. We consider the two dimensional L critical nonlinear Schrödinger equation i

More information

The Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems

The Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems The Fattorini Criterion for the Stabilizability of Parabolic Systems and its Application to MHD flow and fluid-rigid body interaction systems Mehdi Badra Laboratoire LMA, université de Pau et des Pays

More information

Dark and gray spatial optical solitons in Kerr-type nonlocal media

Dark and gray spatial optical solitons in Kerr-type nonlocal media Dark and gray spatial optical solitons in Kerr-type nonlocal media Shigen Ouyang and Qi Guo Laboratory of Photonic Information Technology, South China Normal University, Guangzhou, 510631, P. R. China

More information

Waves and characteristics: Overview 5-1

Waves and characteristics: Overview 5-1 Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables

More information

Bayesian inverse problems with Laplacian noise

Bayesian inverse problems with Laplacian noise Bayesian inverse problems with Laplacian noise Remo Kretschmann Faculty of Mathematics, University of Duisburg-Essen Applied Inverse Problems 2017, M27 Hangzhou, 1 June 2017 1 / 33 Outline 1 Inverse heat

More information

Friedrich symmetric systems

Friedrich symmetric systems viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary

More information

Modeling of Variable Lamé s Modulii for a FGM Generalized Thermoelastic Half Space

Modeling of Variable Lamé s Modulii for a FGM Generalized Thermoelastic Half Space 75 Modeling of Variable Lamé s Modulii for a FGM Generalized Thermoelastic Half Space Abstract In this work we consider a problem in the contet of the generalized theory of thermoelasticity for a half

More information

= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer.

= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer. Problem Sheet,. i) Draw the graphs for [] and {}. ii) Show that for α R, α+ α [t] dt = α and α+ α {t} dt =. Hint Split these integrals at the integer which must lie in any interval of length, such as [α,

More information

Quasi-neutral limit of the non-isentropic Euler Poisson system

Quasi-neutral limit of the non-isentropic Euler Poisson system Proceedings of the Royal Society of Edinburgh, 136A, 1013 1026, 2006 Quasi-neutral limit of the non-isentropic Euler Poisson system Yue-Jun Peng Laboratoire de Mathématiques, CNRS UMR 6620, Université

More information

ON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM

ON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM Indian J. Pure Appl. Math., 44(3): 343-354, June 2013 c Indian National Science Academy ON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM S. Mesloub and T. Hayat Mathematics Department, College

More information

STABILITY OF LARGE-AMPLITUDE SHOCK WAVES OF COMPRESSIBLE NAVIER STOKES EQUATIONS

STABILITY OF LARGE-AMPLITUDE SHOCK WAVES OF COMPRESSIBLE NAVIER STOKES EQUATIONS STABILITY OF LARGE-AMPLITUDE SHOCK WAVES OF COMPRESSIBLE NAVIER STOKES EQUATIONS Kevin Zumbrun Mathematics Department, Indiana University Bloomington, IN 4745-43, USA email: kzumbrun@indiana.edu September

More information

Well-Posedness and Adiabatic Limit for Quantum Zakharov System

Well-Posedness and Adiabatic Limit for Quantum Zakharov System Well-Posedness and Adiabatic Limit for Quantum Zakharov System Yung-Fu Fang (joint work with Tsai-Jung Chen, Jun-Ichi Segata, Hsi-Wei Shih, Kuan-Hsiang Wang, Tsung-fang Wu) Department of Mathematics National

More information

Coercivity of high-frequency scattering problems

Coercivity of high-frequency scattering problems Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.

More information

Inexact Alternating Direction Method of Multipliers for Separable Convex Optimization

Inexact Alternating Direction Method of Multipliers for Separable Convex Optimization Inexact Alternating Direction Method of Multipliers for Separable Convex Optimization Hongchao Zhang hozhang@math.lsu.edu Department of Mathematics Center for Computation and Technology Louisiana State

More information

Nonlinear eigenvalue problems for higher order model equations

Nonlinear eigenvalue problems for higher order model equations Nonlinear eigenvalue problems for higher order model equations L.A. Peletier Mathematical Institute, Leiden University PB 951, 3 RA Leiden, The Netherlands e-mail: peletier@math.leidenuniv.nl 1 Introduction

More information

Conical Shock Waves for Isentropic Euler System

Conical Shock Waves for Isentropic Euler System Conical Shock Waves for Isentropic Euler System Shuxing Chen Institute of Mathematical Research, Fudan University, Shanghai, China E-mail: sxchen@public8.sta.net.cn Dening Li Department of Mathematics,

More information

Lecture 2 Supplementary Notes: Derivation of the Phase Equation

Lecture 2 Supplementary Notes: Derivation of the Phase Equation Lecture 2 Supplementary Notes: Derivation of the Phase Equation Michael Cross, 25 Derivation from Amplitude Equation Near threshold the phase reduces to the phase of the complex amplitude, and the phase

More information

The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations

The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan

More information

The Interior Transmission Eigenvalue Problem for Maxwell s Equations

The Interior Transmission Eigenvalue Problem for Maxwell s Equations The Interior Transmission Eigenvalue Problem for Maxwell s Equations Andreas Kirsch MSRI 2010 epartment of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research

More information

The De Giorgi-Nash-Moser Estimates

The De Giorgi-Nash-Moser Estimates The De Giorgi-Nash-Moser Estimates We are going to discuss the the equation Lu D i (a ij (x)d j u) = 0 in B 4 R n. (1) The a ij, with i, j {1,..., n}, are functions on the ball B 4. Here and in the following

More information

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory

Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Gennady El 1, Roger Grimshaw 1 and Noel Smyth 2 1 Loughborough University, UK, 2 University of Edinburgh, UK

More information

Natural States and Symmetry Properties of. Two-Dimensional Ciarlet-Mooney-Rivlin. Nonlinear Constitutive Models

Natural States and Symmetry Properties of. Two-Dimensional Ciarlet-Mooney-Rivlin. Nonlinear Constitutive Models Natural States and Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Nonlinear Constitutive Models Alexei Cheviakov, Department of Mathematics and Statistics, Univ. Saskatchewan, Canada Jean-François

More information