Gaussian Beam Approximations of High Frequency Waves
|
|
- Deborah Page
- 5 years ago
- Views:
Transcription
1 Gaussian Beam Approximations of High Frequency Waves Olof Runborg Department of Mathematics, KTH Joint with Hailiang Liu, Iowa, James Ralston, UCLA, Nick Tanushev, Z-Terra Corp. ESF Exploratory Workshop IPP Garching, October 2013 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
2 High frequency waves Cauchy problem for scalar wave equation u tt c(x) 2 u = 0, (t, x) R + R n, u(0, x) = A(x)e iφ(x)/ε, u t (0, x) = 1 ε B(x)eiφ(x)/ε, where c(x) (variable) smooth speed of propagation. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
3 High frequency waves Cauchy problem for scalar wave equation u tt c(x) 2 u = 0, (t, x) R + R n, u(0, x) = A(x)e iφ(x)/ε, u t (0, x) = 1 ε B(x)eiφ(x)/ε, where c(x) (variable) smooth speed of propagation. High frequency short wave length highly oscillatory solutions many gridpoints. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
4 High frequency waves Cauchy problem for scalar wave equation u tt c(x) 2 u = 0, (t, x) R + R n, u(0, x) = A(x)e iφ(x)/ε, u t (0, x) = 1 ε B(x)eiφ(x)/ε, where c(x) (variable) smooth speed of propagation. High frequency short wave length highly oscillatory solutions many gridpoints. Direct numerical solution resolves wavelength: #gridpoints ε n at least cost ε n 1 at least Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
5 High frequency waves Cauchy problem for scalar wave equation u tt c(x) 2 u = 0, (t, x) R + R n, u(0, x) = A(x)e iφ(x)/ε, u t (0, x) = 1 ε B(x)eiφ(x)/ε, where c(x) (variable) smooth speed of propagation. High frequency short wave length highly oscillatory solutions many gridpoints. Direct numerical solution resolves wavelength: #gridpoints ε n at least cost ε n 1 at least Often unrealistic approach for applications in e.g. optics, electromagnetics, geophysics, acoustics, quantum mechanics,... Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
6 High frequency material Cauchy problem for scalar wave equation u tt c ε (x) u = 0, (t, x) R + R n, u(0, x) = A(x), u t (0, x) = B(x), where c ε (x) R d d has variations on length scale ε. The functions A(x) and B(x) are smooth (and independent of ε). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
7 High frequency material Cauchy problem for scalar wave equation u tt c ε (x) u = 0, (t, x) R + R n, u(0, x) = A(x), u t (0, x) = B(x), where c ε (x) R d d has variations on length scale ε. The functions A(x) and B(x) are smooth (and independent of ε). Direct numerical solution resolves wavelength: #gridpoints ε n at least cost ε n 1 at least Prohibitively expensive when ε 1, particularly in higher dimensions. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
8 High frequency material Cauchy problem for scalar wave equation u tt c ε (x) u = 0, (t, x) R + R n, u(0, x) = A(x), u t (0, x) = B(x), where c ε (x) R d d has variations on length scale ε. The functions A(x) and B(x) are smooth (and independent of ε). Direct numerical solution resolves wavelength: #gridpoints ε n at least cost ε n 1 at least Prohibitively expensive when ε 1, particularly in higher dimensions. Our approach: Heterogeneous Multiscale Method (HMM) [E,Engquist,2001]. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
9 High frequency material Cauchy problem for scalar wave equation u tt c ε (x) u = 0, (t, x) R + R n, u(0, x) = A(x), u t (0, x) = B(x), where c ε (x) R d d has variations on length scale ε. The functions A(x) and B(x) are smooth (and independent of ε). Direct numerical solution resolves wavelength: #gridpoints ε n at least cost ε n 1 at least Prohibitively expensive when ε 1, particularly in higher dimensions. Our approach: Heterogeneous Multiscale Method (HMM) [E,Engquist,2001]. Solve small micro problems (localized in time and space) to probe effective dynamics, which is approximated on coarse grid ( x ε). Method cost (essentially) independent of ε. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
10 Geometrical optics Wave equation u tt c(x) 2 u = 0. Write solution on the form u(t, x) = a(t, x, ε)e iφ(t,x)/ε. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
11 Geometrical optics Wave equation u tt c(x) 2 u = 0. Write solution on the form u(t, x) = a(t, x, ε)e iφ(t,x)/ε. (a) Amplitude a(x) (b) Phase φ(x) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
12 Geometrical optics a, φ vary on a much coarser scale than u. (And varies little with ε.) Geometrical optics approximation considers a and φ as ε 0. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
13 Geometrical optics a, φ vary on a much coarser scale than u. (And varies little with ε.) Geometrical optics approximation considers a and φ as ε 0. Phase and amplitude satisfy eikonal and transport equations φ 2 t c(y) 2 φ 2 = 0, φ a a t + c φ + c2 φ φ tt a = 0. 2c φ Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
14 Geometrical optics a, φ vary on a much coarser scale than u. (And varies little with ε.) Geometrical optics approximation considers a and φ as ε 0. Phase and amplitude satisfy eikonal and transport equations φ 2 t c(y) 2 φ 2 = 0, φ a a t + c φ + c2 φ φ tt a = 0. 2c φ Ray tracing: x(t), p(t) bicharacteristics of the eikonal equation, dx dt = c(x) 2 p, dp dt = c(x) c(x), φ(t, x(t)) = φ(0, x(0)). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
15 Geometrical optics a, φ vary on a much coarser scale than u. (And varies little with ε.) Geometrical optics approximation considers a and φ as ε 0. Phase and amplitude satisfy eikonal and transport equations φ 2 t c(y) 2 φ 2 = 0, φ a a t + c φ + c2 φ φ tt a = 0. 2c φ Ray tracing: x(t), p(t) bicharacteristics of the eikonal equation, dx dt = c(x) 2 p, dp dt = c(x) c(x), φ(t, x(t)) = φ(0, x(0)). Good accuracy for small ε. Computational cost ε-independent. u(t, x) = a(t, x)e iφ(t,x)/ε + O(ε). (#DOF and cost independent of ε) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
16 Geometrical optics The ansatz u(t, x) = a(t, x)e iφ(t,x)/ε, generally breaks down in finite time if valid at t = 0. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
17 Geometrical optics The ansatz u(t, x) = a(t, x)e iφ(t,x)/ε, generally breaks down in finite time if valid at t = 0. Refraction of waves gives rise to multiple crossing waves N u(t, x) = a n (t, x)e iφn(t,x)/ε n=1 Several amplitude and phase functions. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
18 Geometrical optics The ansatz u(t, x) = a(t, x)e iφ(t,x)/ε, generally breaks down in finite time if valid at t = 0. Refraction of waves gives rise to multiple crossing waves N u(t, x) = a n (t, x)e iφn(t,x)/ε n=1 Several amplitude and phase functions. Caustics appear at points of transition = concentration of rays. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
19 Geometrical optics The ansatz u(t, x) = a(t, x)e iφ(t,x)/ε, generally breaks down in finite time if valid at t = 0. Refraction of waves gives rise to multiple crossing waves N u(t, x) = a n (t, x)e iφn(t,x)/ε n=1 Several amplitude and phase functions. Caustics appear at points of transition = concentration of rays. Geometrical optics predicts infinite amplitude at caustics. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
20 Geometrical optics The ansatz u(t, x) = a(t, x)e iφ(t,x)/ε, generally breaks down in finite time if valid at t = 0. Refraction of waves gives rise to multiple crossing waves u(t, x) = N a n (t, x)e iφn(t,x)/ε n=1 Several amplitude and phase functions. Caustics appear at points of transition = concentration of rays. Geometrical optics predicts infinite amplitude at caustics. Handling multiphase solutions tricky for numerical methods with fixed grids. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
21 Caustics Concentration of rays. GO amplitude a(t, y) but should be a(t, y) ε α, 0 < α < 1. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
22 Gaussian beams Approximate, localized, solutions to the wave equation/schrodinger with a Gaussian profile (width ε). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
23 Gaussian beams Approximate, localized, solutions to the wave equation/schrodinger with a Gaussian profile (width ε). Studied in e.g. Geophysics [Cerveny, Popov, Babich, Psencik, Klimes, Kravtsov,... ], Quantum Mechanics,[Heller, Hagedorn, Herman, Kluk, Kay,... ], Plasma Physics, [Pereverzev, Peeters, Maj,... ], Mathematics [Ralston, Hörmander,... ] Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
24 Gaussian beams Approximate, localized, solutions to the wave equation/schrodinger with a Gaussian profile (width ε). Studied in e.g. Geophysics [Cerveny, Popov, Babich, Psencik, Klimes, Kravtsov,... ], Quantum Mechanics,[Heller, Hagedorn, Herman, Kluk, Kay,... ], Plasma Physics, [Pereverzev, Peeters, Maj,... ], Mathematics [Ralston, Hörmander,... ] No breakdown at caustics. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
25 Gaussian beams Gaussian beams are of the same form as geometrical optics solutions, v(t, y) = A(t, y)e iφ(t,y)/ε, centered around a geometrical optics ray x(t), A(t, y) = a(t, y x(t)), Φ(t, y) = φ(t, y x(t)). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
26 Gaussian beams Gaussian beams are of the same form as geometrical optics solutions, v(t, y) = A(t, y)e iφ(t,y)/ε, centered around a geometrical optics ray x(t), A(t, y) = a(t, y x(t)), Φ(t, y) = φ(t, y x(t)). The phase Φ will now have a positive imaginary part away from the ray x(t). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
27 Gaussian beams Gaussian beams are of the same form as geometrical optics solutions, v(t, y) = A(t, y)e iφ(t,y)/ε, centered around a geometrical optics ray x(t), A(t, y) = a(t, y x(t)), Φ(t, y) = φ(t, y x(t)). The phase Φ will now have a positive imaginary part away from the ray x(t). Imaginary part of φ y 2 v(t, y) e y x(t) 2 /ε, Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
28 Gaussian beams Gaussian beams are of the same form as geometrical optics solutions, v(t, y) = A(t, y)e iφ(t,y)/ε, centered around a geometrical optics ray x(t), A(t, y) = a(t, y x(t)), Φ(t, y) = φ(t, y x(t)). The phase Φ will now have a positive imaginary part away from the ray x(t). Imaginary part of φ y 2 v(t, y) e y x(t) 2 /ε, Gaussian with width ε Localized around x(t). (Moves along the space time ray.) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
29 Gaussian beams Gaussian beams are of the same form as geometrical optics solutions, v(t, y) = A(t, y)e iφ(t,y)/ε, centered around a geometrical optics ray x(t), A(t, y) = a(t, y x(t)), Φ(t, y) = φ(t, y x(t)). The phase Φ will now have a positive imaginary part away from the ray x(t). Imaginary part of φ y 2 v(t, y) e y x(t) 2 /ε, Gaussian with width ε Localized around x(t). (Moves along the space time ray.) Phase Φ(t, y) and amplitude A(t, y) approximated by polynomials locally around x(t) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
30 Gaussian beams Gaussian beams are of the same form as geometrical optics solutions, v(t, y) = A(t, y)e iφ(t,y)/ε, centered around a geometrical optics ray x(t), A(t, y) = a(t, y x(t)), Φ(t, y) = φ(t, y x(t)). The phase Φ will now have a positive imaginary part away from the ray x(t). Imaginary part of φ y 2 v(t, y) e y x(t) 2 /ε, Gaussian with width ε Localized around x(t). (Moves along the space time ray.) Phase Φ(t, y) and amplitude A(t, y) approximated by polynomials locally around x(t) Φ(t, y) and A(t, y) solve eikonal and transport equation only upto O( y x m ). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
31 First order beams The simplest ("first order") Gaussian beams are of the form v(t, y) = a 0 (t)e iφ(t,y)/ε, Φ(t, y) = φ(t, y x(t)), where φ(t, y) = φ 0 (t) + y p(t) y M(t)y. i.e. A(t, y) approximated to 0th order, and Φ(t, y) to 2nd order. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
32 First order beams The simplest ("first order") Gaussian beams are of the form v(t, y) = a 0 (t)e iφ(t,y)/ε, Φ(t, y) = φ(t, y x(t)), where φ(t, y) = φ 0 (t) + y p(t) y M(t)y. i.e. A(t, y) approximated to 0th order, and Φ(t, y) to 2nd order. We require that Φ(t, y) solves eikonal to order O( y x 3 ) and A(t, y) solves transport equation to order O( y x ). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
33 First order beams Let us thus require that Φ A A t + c Φ Φ 2 t c(y) 2 Φ 2 = O( y x(t) 3 ), + c2 Φ Φ tt A = O( y x(t) ), 2c Φ Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
34 First order beams Let us thus require that Φ A A t + c Φ We obtain ODEs for φ 0, x, p, M, a 0. ẋ(t) = c(x) 2 p, φ 0 (t) = 0, Φ 2 t c(y) 2 Φ 2 = O( y x(t) 3 ), + c2 Φ Φ tt A = O( y x(t) ), 2c Φ ṗ(t) = c(x)/c(x), Ṁ(t) = D MB B T M MCM, ȧ 0 (t) = a ( ) 0 c(x)p c(x) c(x) 3 p Mp + c(x) 2 Tr[M], 2 where B, C, D are matrix functions involving x, p and c(x). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
35 First order beams Let us thus require that Φ A A t + c Φ We obtain ODEs for φ 0, x, p, M, a 0. ẋ(t) = c(x) 2 p, φ 0 (t) = 0, Φ 2 t c(y) 2 Φ 2 = O( y x(t) 3 ), + c2 Φ Φ tt A = O( y x(t) ), 2c Φ ṗ(t) = c(x)/c(x), Ṁ(t) = D MB B T M MCM, ȧ 0 (t) = a ( ) 0 c(x)p c(x) c(x) 3 p Mp + c(x) 2 Tr[M], 2 where B, C, D are matrix functions involving x, p and c(x). ODEs easy to solve numerically. Beams easy to evaluate: v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
36 First order beams Let us thus require that Φ A A t + c Φ We obtain ODEs for φ 0, x, p, M, a 0. ẋ(t) = c(x) 2 p, φ 0 (t) = 0, Φ 2 t c(y) 2 Φ 2 = O( y x(t) 3 ), + c2 Φ Φ tt A = O( y x(t) ), 2c Φ ṗ(t) = c(x)/c(x), Ṁ(t) = D MB B T M MCM, ȧ 0 (t) = a ( ) 0 c(x)p c(x) c(x) 3 p Mp + c(x) 2 Tr[M], 2 Asymptotic order of accuracy is ( ) 1 v tt c(y) 2 v = O. ε Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
37 Higher order beams More generally, we can construct higher order beams. Let where, for order K beams, v(t, y) = a(t, y x(t))e iφ(t,y x(t))/ε, Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
38 Higher order beams More generally, we can construct higher order beams. Let where, for order K beams, v(t, y) = a(t, y x(t))e iφ(t,y x(t))/ε, The phase is a Taylor polynomial of order K + 1, φ(t, y) = φ 0 (t) + y p(t) + y 1 K +1 2 M(t)y + β =3 1 β! φ β(t)y β. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
39 Higher order beams More generally, we can construct higher order beams. Let where, for order K beams, v(t, y) = a(t, y x(t))e iφ(t,y x(t))/ε, The phase is a Taylor polynomial of order K + 1, φ(t, y) = φ 0 (t) + y p(t) + y 1 K +1 2 M(t)y + A is now a finite WKB expansion, a(t, y) = K /2 1 j=0 ε j a j (t, y) β =3 1 β! φ β(t)y β. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
40 Higher order beams More generally, we can construct higher order beams. Let where, for order K beams, v(t, y) = a(t, y x(t))e iφ(t,y x(t))/ε, The phase is a Taylor polynomial of order K + 1, φ(t, y) = φ 0 (t) + y p(t) + y 1 K +1 2 M(t)y + A is now a finite WKB expansion, a(t, y) = K /2 1 j=0 ε j a j (t, y) β =3 1 β! φ β(t)y β. Each amplitude term a j is a Taylor polynomial to order K 2j 1 a j (t, y) = K 2j 1 β =0 1 β! a j,β(t)y β Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
41 Higher order beams We now require that Φ(t, y) = φ(t, y x)) solves eikonal equation to order y x K +2 a j (t, y x) solve higher order transport equations to order y x K 2j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
42 Higher order beams We now require that Φ(t, y) = φ(t, y x)) solves eikonal equation to order y x K +2 a j (t, y x) solve higher order transport equations to order y x K 2j Again, this gives ODEs for all Taylor coefficients, ẋ(t) = c(x) 2 p, φ0 (t) = 0, ṗ(t) = c(x)/c(x), Ṁ(t) = D MB B T M MCM, ȧ j,β (t) =..., φ β (t) =..., Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
43 Higher order beams We now require that Φ(t, y) = φ(t, y x)) solves eikonal equation to order y x K +2 a j (t, y x) solve higher order transport equations to order y x K 2j Again, this gives ODEs for all Taylor coefficients, ẋ(t) = c(x) 2 p, φ0 (t) = 0, ṗ(t) = c(x)/c(x), Ṁ(t) = D MB B T M MCM, ȧ j,β (t) =..., φ β (t) =..., Asymptotic order of accuracy is ( v tt c(y) 2 v = O ε K /2 1). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
44 Gaussian beams Properties v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
45 Gaussian beams Properties v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y Φ(t, x(t)) = φ(t, 0) = φ 0 (t) is real valued Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
46 Gaussian beams Properties v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y Φ(t, x(t)) = φ(t, 0) = φ 0 (t) is real valued If M(0) is symmetric and IM(0) is positive definite then this is true for M(t) (which exists) for all t > 0. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
47 Gaussian beams Properties v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y Φ(t, x(t)) = φ(t, 0) = φ 0 (t) is real valued If M(0) is symmetric and IM(0) is positive definite then this is true for M(t) (which exists) for all t > 0. a 0 (t) exists everywhere (no blow-up at caustics) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
48 Gaussian beams Properties v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y Φ(t, x(t)) = φ(t, 0) = φ 0 (t) is real valued If M(0) is symmetric and IM(0) is positive definite then this is true for M(t) (which exists) for all t > 0. a 0 (t) exists everywhere (no blow-up at caustics) Shape of beam remains Gaussian Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
49 Gaussian beams Properties v(t, y) = a 0 (t)e iφ(t,y x(t))/ε, φ(t, y) = φ 0 (t) + y p(t) y M(t)y Φ(t, x(t)) = φ(t, 0) = φ 0 (t) is real valued If M(0) is symmetric and IM(0) is positive definite then this is true for M(t) (which exists) for all t > 0. a 0 (t) exists everywhere (no blow-up at caustics) Shape of beam remains Gaussian For high order need cutoff in a neighborhood of central ray to avoid spurious growth. v(t, y) = a(t, y x(t))e iφ(t,y x(t))/ε ϱ(y x(t)) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
50 Superpositions of Gaussian beams To approximate more general solutions, use superpositions of beams. Let v(t, y; z) be a beam starting from the point y = z and define u GB (t, y) = ε n 2 v(t, y; z)dz K 0 (n dimension, K 0 compact set) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
51 Superpositions of Gaussian beams To approximate more general solutions, use superpositions of beams. Let v(t, y; z) be a beam starting from the point y = z and define u GB (t, y) = ε n 2 a 0 (t; z)e i ε [φ 0(t;z)+(y x(t;z)) p(t;z)+ K 0 (n dimension, K 0 compact set) 1 2 (y x(t;z) M(t;z)(y x(t;z))] Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
52 Superpositions of Gaussian beams To approximate more general solutions, use superpositions of beams. Let v(t, y; z) be a beam starting from the point y = z and define u GB (t, y) = ε n 2 a 0 (t; z)e i ε [φ 0(t;z)+(y x(t;z)) p(t;z)+ K 0 (n dimension, K 0 compact set) 1 2 (y x(t;z) M(t;z)(y x(t;z))] By linearity of the wave equation equation a sum of solutions is also a solution. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
53 Superpositions of Gaussian beams To approximate more general solutions, use superpositions of beams. Let v(t, y; z) be a beam starting from the point y = z and define u GB (t, y) = ε n 2 a 0 (t; z)e i ε [φ 0(t;z)+(y x(t;z)) p(t;z)+ K 0 (n dimension, K 0 compact set) 1 2 (y x(t;z) M(t;z)(y x(t;z))] By linearity of the wave equation equation a sum of solutions is also a solution. u GB (t, y) is an asymptotic solution with initial data u GB (0, y). Sufficient to describe e.g. WKB data: K -th order beams s.t. A(y)e iφ(y)/ε u GB (0, ) = O(ε K /2 ), E [Tanushev, 2007]. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
54 Superpositions of Gaussian beams To approximate more general solutions, use superpositions of beams. Let v(t, y; z) be a beam starting from the point y = z and define u GB (t, y) = ε n 2 a 0 (t; z)e i ε [φ 0(t;z)+(y x(t;z)) p(t;z)+ K 0 (n dimension, K 0 compact set) 1 2 (y x(t;z) M(t;z)(y x(t;z))] By linearity of the wave equation equation a sum of solutions is also a solution. u GB (t, y) is an asymptotic solution with initial data u GB (0, y). Sufficient to describe e.g. WKB data: K -th order beams s.t. A(y)e iφ(y)/ε u GB (0, ) = O(ε K /2 ), E [Tanushev, 2007]. Prefactor normalizes beams appropriately, u GB E = O(1). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
55 Superpositions of Gaussian beams More general phase space superposition: Let v(t, y; z, p) be a beam starting from the point y = z with momentum p and define u GB (t, y) = ε n K 0 v(t, y; z, p)dzdp. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
56 Superpositions of Gaussian beams More general phase space superposition: Let v(t, y; z, p) be a beam starting from the point y = z with momentum p and define u GB (t, y) = ε n K 0 v(t, y; z, p)dzdp. u GB (t, y) is an asymptotic solution with initial data u GB (0, y) = ε n K 0 v(0, y; z, p)dzdp. Can describe more general data. (C.f. FBI transform.) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
57 Numerical methods Approximate superposition integral by sum (trapezoidal rule) u GB (t, y) = ε n 2 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. K 0 j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
58 Numerical methods Approximate superposition integral by sum (trapezoidal rule) u GB (t, y) = ε n 2 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. K 0 Lagrangian methods Solve ODEs z j with standard methods. Similar to ray tracing but with all the additional Taylor coefficients computed along the rays (M, a j,β, φ β,... ) [Hill, Klimes,... ] j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
59 Numerical methods Approximate superposition integral by sum (trapezoidal rule) u GB (t, y) = ε n 2 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. K 0 Lagrangian methods Solve ODEs z j with standard methods. Similar to ray tracing but with all the additional Taylor coefficients computed along the rays (M, a j,β, φ β,... ) [Hill, Klimes,... ] Eulerian methods obtain parameters from solving PDEs on fixed grids [Leung, Qian, Burridge,07]], [Jin, Wu, Yang,08], [Jin, Wu, Yang, Huang, 09], [Leung, Qian,09], [Qian,Ying,10],... j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
60 Numerical methods Approximate superposition integral by sum (trapezoidal rule) u GB (t, y) = ε n 2 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. K 0 Lagrangian methods Solve ODEs z j with standard methods. Similar to ray tracing but with all the additional Taylor coefficients computed along the rays (M, a j,β, φ β,... ) [Hill, Klimes,... ] Eulerian methods obtain parameters from solving PDEs on fixed grids [Leung, Qian, Burridge,07]], [Jin, Wu, Yang,08], [Jin, Wu, Yang, Huang, 09], [Leung, Qian,09], [Qian,Ying,10],... Wavefront methods solve for parameters on a wave front [Motamed, OR,09] j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
61 Numerical methods u GB (t, y) = ε n 2 Numerical issues K 0 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. Cost number of beams since each beam is O(1). For accuracy need z ε width of beams. cost O(ε n/2 ) C.f. direct solution of wave equations, at least O(ε (n+1) ) j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
62 Numerical methods u GB (t, y) = ε n 2 Numerical issues K 0 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. Cost number of beams since each beam is O(1). For accuracy need z ε width of beams. cost O(ε n/2 ) C.f. direct solution of wave equations, at least O(ε (n+1) ) For phase space superposition would get O(ε n ) but can often be improved (e.g. support in p ε for WKB data) j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
63 Numerical methods u GB (t, y) = ε n 2 Numerical issues K 0 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. Cost number of beams since each beam is O(1). For accuracy need z ε width of beams. cost O(ε n/2 ) C.f. direct solution of wave equations, at least O(ε (n+1) ) For phase space superposition would get O(ε n ) but can often be improved (e.g. support in p ε for WKB data) Spreading of beams Wide beams large Taylor approximation errors j Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
64 Numerical methods u GB (t, y) = ε n 2 Numerical issues K 0 v(t, y; z)dz ε n 2 v(t, y; z j ) z n. Cost number of beams since each beam is O(1). For accuracy need z ε width of beams. cost O(ε n/2 ) C.f. direct solution of wave equations, at least O(ε (n+1) ) For phase space superposition would get O(ε n ) but can often be improved (e.g. support in p ε for WKB data) Spreading of beams Wide beams large Taylor approximation errors Initial data approximation Many degrees of freedom. Can have huge impact on accuracy at later times. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36 j
65 Approximation errors Let := tt + c(y) 2. Suppose u is exact solution of wave equation and ũ is the Gaussian beam approximation u = 0, ũ = O(ε K /2 1 ). What is the norm error in ũ, i.e. u ũ? Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
66 Approximation errors Let := tt + c(y) 2. Suppose u is exact solution of wave equation and ũ is the Gaussian beam approximation u = 0, ũ = O(ε K /2 1 ). What is the norm error in ũ, i.e. u ũ? Use ε-scaled energy norm u 2 E := ε2 2 This is O(1) for WKB type initial data, R n u t 2 c(y) 2 + u 2 dy. u(0, x) = A(x)e iφ(x)/ε u(0, ) E = O(1). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
67 Approximation errors Use well-posedness (stability) estimate for wave equation solutions w: w(t, ) E w(0, ) E + εc(t ) sup w(t, ) L 2, 0 t T. t [0,T ] Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
68 Approximation errors Use well-posedness (stability) estimate for wave equation solutions w: w(t, ) E w(0, ) E + εc(t ) sup w(t, ) L 2, 0 t T. t [0,T ] Since, u = 0 and by linearity Hence, assuming u(0, x) = ũ(0, x), [ũ u] = ũ. ũ(t, ) u(t, ) E εc(t ) sup ũ(t, ) L 2, 0 t T. t [0,T ] Error in ũ how well it satisfies equation, plus one order in ε Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
69 Approximation Errors Single Gaussian Beam By earlier construction and ũ GB (t, x) = O(ε K /2 1 ). ũ(t, ) u(t, ) E εc(t ) sup ũ(t, ) L 2, 0 t T. t [0,T ] Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
70 Approximation Errors Single Gaussian Beam By earlier construction and ũ GB (t, x) = O(ε K /2 1 ). ũ(t, ) u(t, ) E εc(t ) sup ũ(t, ) L 2, 0 t T. t [0,T ] Hence, after appropriate normalization (so ũ E = 1) and assuming zero initial data error, [Ralston, 82] Estimate is sharp. sup ũ(t, ) u(t, ) E O(ε K /2 ). t [0,T ] Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
71 Approximation Errors Superpositions Norm estimates of u u GB only rather recently derived [Swart, Rousse, Liu, Ralston, Tanushev, Bougacha, Alexandre, Lu, Yang,... ] Need to check how well u GB satisfies equation u GB := tt u GB + c(y) 2 u GB, u GB (t, y) = ε n 2 K 0 v(t, y; z)dz. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
72 Approximation Errors Superpositions Norm estimates of u u GB only rather recently derived [Swart, Rousse, Liu, Ralston, Tanushev, Bougacha, Alexandre, Lu, Yang,... ] Need to check how well u GB satisfies equation u GB := tt u GB + c(y) 2 u GB, u GB (t, y) = ε n 2 By linearity, u GB = ε n/2 K 0 v(t, y; z)dz. K 0 v(t, y; z)dz. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
73 Approximation Errors Superpositions Norm estimates of u u GB only rather recently derived [Swart, Rousse, Liu, Ralston, Tanushev, Bougacha, Alexandre, Lu, Yang,... ] Need to check how well u GB satisfies equation u GB := tt u GB + c(y) 2 u GB, u GB (t, y) = ε n 2 By linearity and simple estimate, u GB ε n/2 K 0 v(t, y; z) dz. K 0 v(t, y; z)dz. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
74 Approximation Errors Superpositions Norm estimates of u u GB only rather recently derived [Swart, Rousse, Liu, Ralston, Tanushev, Bougacha, Alexandre, Lu, Yang,... ] Need to check how well u GB satisfies equation u GB := tt u GB + c(y) 2 u GB, u GB (t, y) = ε n 2 By linearity and simple estimate, u GB ε n/2 K 0 v(t, y; z) dz. After scaling and Cauchy Schwarz, K 0 v(t, y; z)dz. u GB (y) 2 2 Cε n K 0 v(t, y; z) 2 2 dz Cε n ε K 2+n/2. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
75 Approximation Errors Superpositions Norm estimates of u u GB only rather recently derived [Swart, Rousse, Liu, Ralston, Tanushev, Bougacha, Alexandre, Lu, Yang,... ] Need to check how well u GB satisfies equation u GB := tt u GB + c(y) 2 u GB, u GB (t, y) = ε n 2 By linearity and simple estimate, u GB ε n/2 K 0 v(t, y; z) dz. After scaling and Cauchy Schwarz, K 0 v(t, y; z)dz. u GB (y) 2 2 Cε n K 0 v(t, y; z) 2 2 dz Cε n ε K 2+n/2. Gives error estimate for u GB u u GB E εc u GB 2 Cε K /2 n/4. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
76 Approximation errors superpositions Basic estimate u(t, ) u GB (t, ) E O(ε K /2 n/4 ). is not sharp. E.g. it does not predict convergence for first order beams in 2D. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
77 Approximation errors superpositions Basic estimate u(t, ) u GB (t, ) E O(ε K /2 n/4 ). is not sharp. E.g. it does not predict convergence for first order beams in 2D. Problem: The step u GB ε n/2 K 0 v(t, y; z) dz ignores cancellations between neighbouring beams. Very bad except at caustics where beams interfere constructively. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
78 Approximation errors superpositions Basic estimate u(t, ) u GB (t, ) E O(ε K /2 n/4 ). is not sharp. E.g. it does not predict convergence for first order beams in 2D. Problem: The step u GB ε n/2 K 0 v(t, y; z) dz ignores cancellations between neighbouring beams. Very bad except at caustics where beams interfere constructively. Gives a dependence on dimension n in estimate. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
79 Error estimate Theorem (Liu, Tanushev, O.R.,2010) For the wave equation, u(t, ) u GB (t, ) E O(ε K /2 ). For the Schrödinger equation, u(t, ) u GB (t, ) L 2 O(ε K /2 ). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
80 Error estimate Theorem (Liu, Tanushev, O.R.,2010) For the wave equation, u(t, ) u GB (t, ) E O(ε K /2 ). For the Schrödinger equation, u(t, ) u GB (t, ) L 2 O(ε K /2 ). Superposition in physical space. Initial data approximated on a submanifold of phase space (WKB data). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
81 Error estimate Theorem (Liu, Tanushev, O.R.,2010) For the wave equation, u(t, ) u GB (t, ) E O(ε K /2 ). For the Schrödinger equation, u(t, ) u GB (t, ) L 2 O(ε K /2 ). Superposition in physical space. Initial data approximated on a submanifold of phase space (WKB data). Convergence of all beams independent of dimension and presence of caustics. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
82 Error estimate Theorem (Liu, Tanushev, O.R.,2010) For the wave equation, u(t, ) u GB (t, ) E O(ε K /2 ). For the Schrödinger equation, u(t, ) u GB (t, ) L 2 O(ε K /2 ). Superposition in physical space. Initial data approximated on a submanifold of phase space (WKB data). Convergence of all beams independent of dimension and presence of caustics. Result also for general scalar, strictly hyperbolic m-th order PDEs. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
83 Error estimate Theorem (Liu, Tanushev, O.R.,2010) For the wave equation, u(t, ) u GB (t, ) E O(ε K /2 ). For the Schrödinger equation, u(t, ) u GB (t, ) L 2 O(ε K /2 ). Superposition in physical space. Initial data approximated on a submanifold of phase space (WKB data). Convergence of all beams independent of dimension and presence of caustics. Result also for general scalar, strictly hyperbolic m-th order PDEs. Cf. [Bougacha, Akian, Alexandre, 2009], [Rousse, Swart, 2009] and [Lu, Yang, 2011] for other settings. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
84 Gaussian beams for Helmholtz For time-harmonic waves consider Helmholtz equation u + (iαε 1 + ε 2 )n 2 u = g, x R d. where n(x) = 1/c(x), α=damping and g supported on a co-dimension one manifold. (Ex. g = g 0 (x 2 )δ(x 1 )/ε.) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
85 Gaussian beams for Helmholtz For time-harmonic waves consider Helmholtz equation u + (iαε 1 + ε 2 )n 2 u = g, x R d. where n(x) = 1/c(x), α=damping and g supported on a co-dimension one manifold. (Ex. g = g 0 (x 2 )δ(x 1 )/ε.) "Blobs" "Fat rays" localized around geometrical optics ray Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
86 Gaussian beams for Helmholtz For time-harmonic waves consider Helmholtz equation u + (iαε 1 + ε 2 )n 2 u = g, x R d. where n(x) = 1/c(x), α=damping and g supported on a co-dimension one manifold. (Ex. g = g 0 (x 2 )δ(x 1 )/ε.) "Blobs" "Fat rays" localized around geometrical optics ray To leading order gaussian transversely to ray Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
87 Gaussian beams for Helmholtz Same ansatz, v = a(s, y x(s))e iφ(s,y x(s))/ε, centered around a geometrical optics ray x(s) but s not time. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
88 Gaussian beams for Helmholtz Same ansatz, v = a(s, y x(s))e iφ(s,y x(s))/ε, centered around a geometrical optics ray x(s) but s not time. First order beams are of the form φ = φ 0 (s) + y p(s) y M(s)y, a = a 0(s), i.e. a approximated to 0th order, and φ to 2nd order. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
89 Gaussian beams for Helmholtz Same ansatz, v = a(s, y x(s))e iφ(s,y x(s))/ε, centered around a geometrical optics ray x(s) but s not time. First order beams are of the form φ = φ 0 (s) + y p(s) y M(s)y, a = a 0(s), i.e. a approximated to 0th order, and φ to 2nd order. Similar ODEs for a 0, x, p, M, φ 0 as in the time-dependent case. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
90 Gaussian beams for Helmholtz Same ansatz, v = a(s, y x(s))e iφ(s,y x(s))/ε, centered around a geometrical optics ray x(s) but s not time. First order beams are of the form φ = φ 0 (s) + y p(s) y M(s)y, a = a 0(s), i.e. a approximated to 0th order, and φ to 2nd order. Similar ODEs for a 0, x, p, M, φ 0 as in the time-dependent case. Similar properties as in time-dependent case: Phase φ evaluated on ray = φ 0 (s) is real valued If M(0) is symmetric and IM(0) is positive definite then this is true for M(s) (which exists) for all s > 0. a 0 (s) exists everywhere (no blow-up at caustics) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
91 Gaussian beams for Helmholtz Extension off ray v(y) = a(s, y x(s))e iφ(s,y x(s))/ε, v(y) y x(s) η x(s*) How to evaluate "(s, y x(s))" in expression for beam? Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
92 Gaussian beams for Helmholtz Extension off ray v(y) = a(s, y x(s ))e iφ(s,y x(s ))/ε, v(y) y x(s) s = s (y) η x(s*) How to evaluate "(s, y x(s))" in expression for beam? No distinguished "time" variable Extend beam by Taylor expansion transversely to ray: Let s = s (y) such that x(s ) is closest point on ray to y. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
93 Gaussian beams for Helmholtz Extension off ray v(y) = a(s, y x(s ))e iφ(s,y x(s ))/ε v(y) y x(s) ϱ(y x(s )), s = s (y) η x(s*) How to evaluate "(s, y x(s))" in expression for beam? No distinguished "time" variable Extend beam by Taylor expansion transversely to ray: Let s = s (y) such that x(s ) is closest point on ray to y. Only well-defined close enough to ray Cutoff ϱ(y) needed also for first order beams (size η) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
94 Gaussian beams for Helmholtz Source Helmholtz with source on Σ = {y : ρ(y) = 0}. u + (iαε 1 + ε 2 )n 2 u = 1 ε g(y)δ(ρ(y)). Σ u (x) x(s) η x 0 u + (x) ρ < 0 ρ > 0 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
95 Gaussian beams for Helmholtz Source Helmholtz with source on Σ = {y : ρ(y) = 0}. u + (iαε 1 + ε 2 )n 2 u = 1 ε g(y)δ(ρ(y)). Σ u (x) x(s) η x 0 u + (x) ρ < 0 ρ > 0 Beams shoot out orthogonally in each direction from Σ Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
96 Gaussian beams for Helmholtz Source Helmholtz with source on Σ = {y : ρ(y) = 0}. u + (iαε 1 + ε 2 )n 2 u = 1 ε g(y)δ(ρ(y)). Σ u (x) x(s) η x 0 u + (x) ρ < 0 ρ > 0 Beams shoot out orthogonally in each direction from Σ Gives beams v ± (y), with v + (y) = 0 when ρ(y) < 0 etc. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
97 Gaussian beams for Helmholtz Source Helmholtz with source on Σ = {y : ρ(y) = 0}. Lu =: u + (iαε 1 + ε 2 )n 2 u = 1 ε g(y)δ(ρ(y)). Σ u (x) x(s) η x 0 u + (x) ρ < 0 ρ > 0 Beams shoot out orthogonally in each direction from Σ Gives beams v ± (y), with v + (y) = 0 when ρ(y) < 0 etc. Note that v + = v on Σ, but φ + = φ so that L(v + + v ) δ(ρ(y))+ smooth part. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
98 Gaussian beams for Helmholtz Superposition Σ Lu =: u + (iαε 1 + ε 2 )n 2 u u (x) x(s) = 1 ε g(y)δ(ρ(y)). η x 0 u + (x) ρ < 0 ρ > 0 Let v ± (y; z) be the beams starting from z Σ and define superposition u GB (y) = ε n 1 2 [v + (y; z) + v (y; z)]da z (1) Σ Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
99 Gaussian beams for Helmholtz Superposition Σ Lu =: u + (iαε 1 + ε 2 )n 2 u u (x) x(s) = 1 ε g(y)δ(ρ(y)). η x 0 u + (x) ρ < 0 ρ > 0 Let v ± (y; z) be the beams starting from z Σ and define superposition u GB (y) = ε n 1 2 [v + (y; z) + v (y; z)]da z (1) Choose initial data for beam v ± (z; z) such that Σ with g g. Lu GB (y) 1 ε g(y)δ(ρ(y)) + f GB Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
100 Error estimate K th order beams, Helmholtz case Theorem (Liu, Ralston, Tanushev, O.R., 2013) Assume Smooth, compactly supported source g(x) Index of refraction n(x) smooth and constant for x > R No trapped rays: L s.t. x(±l) > 2R if x(0) < R, p(0) = n(x(0)) No initial data error g = g. Then with C independent of ε and α, u u GB L2 ( x <R) Cε K /2, Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
101 Error estimate K th order beams, Helmholtz case Theorem (Liu, Ralston, Tanushev, O.R., 2013) Assume Smooth, compactly supported source g(x) Index of refraction n(x) smooth and constant for x > R No trapped rays: L s.t. x(±l) > 2R if x(0) < R, p(0) = n(x(0)) No initial data error g = g. Then with C independent of ε and α, u u GB L2 ( x <R) Cε K /2, Superposition in physical space. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
102 Error estimate K th order beams, Helmholtz case Theorem (Liu, Ralston, Tanushev, O.R., 2013) Assume Smooth, compactly supported source g(x) Index of refraction n(x) smooth and constant for x > R No trapped rays: L s.t. x(±l) > 2R if x(0) < R, p(0) = n(x(0)) No initial data error g = g. Then with C independent of ε and α, u u GB L2 ( x <R) Cε K /2, Superposition in physical space. Convergence of all beams independent of dimension and presence of caustics. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
103 Sketch of proof, wave equation Use energy estimate u GB (t, ) u(t, ) E u GB (0, ) u(0, ) E + Cε sup u GB (t, ) L 2, t [0,T ] Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
104 Sketch of proof, wave equation Use energy estimate u GB (t, ) u(t, ) E u GB (0, ) u(0, ) E + Cε sup u GB (t, ) L 2, t [0,T ] The residual is of the form u GB (t, y) = ε K /2 q J j=1 ε r j T ε j [f j ](t, y) + O(ε ), where r j 0, J finite and f j L 2 (all independent of ε). T ε j : L 2 L 2 belongs to a class of oscillatory integral operators. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
105 Sketch of proof, wave equation Use energy estimate u GB (t, ) u(t, ) E u GB (0, ) u(0, ) E + Cε sup u GB (t, ) L 2, t [0,T ] The residual is of the form u GB (t, y) = ε K /2 q J j=1 ε r j T ε j [f j ](t, y) + O(ε ), where r j 0, J finite and f j L 2 (all independent of ε). T ε j : L 2 L 2 belongs to a class of oscillatory integral operators. Together we get (if initial data exact) u GB (t, ) u(t, ) E C(T )ε K /2 J j=1 ε r j T ε j L 2 f j L 2 + O(ε ) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
106 Sketch of proof, cont. We have u GB (t, ) u(t, ) E C(T )ε K /2 J j=1 T ε j L 2 + O(ε ) where, in its simplest form, T ε [w](t, y) := ε n+ α 2 w(z)(y x(t; z)) α e iφ(t,y x(t;z);z)/ε dz, K 0 for some multi-index α, Gaussian beam phase φ and geometrical optics rays x(t; z) with x(0; z) = z. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
107 Sketch of proof, cont. We have u GB (t, ) u(t, ) E C(T )ε K /2 J j=1 T ε j L 2 + O(ε ) where, in its simplest form, T ε [w](t, y) := ε n+ α 2 w(z)(y x(t; z)) α e iφ(t,y x(t;z);z)/ε dz, K 0 for some multi-index α, Gaussian beam phase φ and geometrical optics rays x(t; z) with x(0; z) = z. Result follows if we prove that T ε is bounded in L 2 independent of ε, T ε L 2 C. This is the key estimate of our proof. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
108 Sketch of proof, cont. Estimate of T ε L 2, where T ε [w](t, y) := ε n+ α 2 w(z)(y x(t; z)) α e iφ(t,y x(t;z);z)/ε dz. K 0 Main difficulty: no globally invertible map x(0; z) = z x(t; z) because of caustics. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
109 Sketch of proof, cont. Estimate of T ε L 2, where T ε [w](t, y) := ε n+ α 2 w(z)(y x(t; z)) α e iφ(t,y x(t;z);z)/ε dz. K 0 Main difficulty: no globally invertible map x(0; z) = z x(t; z) because of caustics. Mapping (x(0; z), p(0; z)) (x(t; z), p(t; z) is however globally invertible and smooth. Gives the "non-squeezing" property, c 1 z z p(t; z) p(t; z ) + x(t; z) x(t; z ) c 2 z z. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
110 Sketch of proof, cont. Estimate of T ε L 2, where T ε [w](t, y) := ε n+ α 2 w(z)(y x(t; z)) α e iφ(t,y x(t;z);z)/ε dz. K 0 Main difficulty: no globally invertible map x(0; z) = z x(t; z) because of caustics. Mapping (x(0; z), p(0; z)) (x(t; z), p(t; z) is however globally invertible and smooth. Gives the "non-squeezing" property, c 1 z z p(t; z) p(t; z ) + x(t; z) x(t; z ) c 2 z z. Allows us to use stationary phase arguments close to caustics, and carefully control cancellations of oscillations there (similar to [Swart,Rousse], [Bougacha, Akian, Alexandre]). Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
111 Approximation errors Remarks The estimate u(t, ) u GB (t, ) E O(ε K /2 ) is sharp for individual beams (relative error). But for superpositions? Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
112 Approximation errors Remarks The estimate u(t, ) u GB (t, ) E O(ε K /2 ) is sharp for individual beams (relative error). But for superpositions? Predicts convergence rate of first order beam to be only O( ε). These beams are based on same high frequency approximation as geometrical optics which has O(ε) accuracy. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
113 Approximation errors Remarks The estimate u(t, ) u GB (t, ) E O(ε K /2 ) is sharp for individual beams (relative error). But for superpositions? Predicts convergence rate of first order beam to be only O( ε). These beams are based on same high frequency approximation as geometrical optics which has O(ε) accuracy. Numerical experiments suggests a better rate for odd order beams Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
114 Approximation errors Remarks The estimate u(t, ) u GB (t, ) E O(ε K /2 ) is sharp for individual beams (relative error). But for superpositions? Predicts convergence rate of first order beam to be only O( ε). These beams are based on same high frequency approximation as geometrical optics which has O(ε) accuracy. Numerical experiments suggests a better rate for odd order beams For the Helmholtz case we have proved [Motamed, OR] that u(x) u GB(x) O(ε K /2 ) for the Taylor expansion part of the error away from caustics. This gives O(ε) for first order beams. Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
115 Approximation errors Remarks The estimate u(t, ) u GB (t, ) E O(ε K /2 ) is sharp for individual beams (relative error). But for superpositions? Predicts convergence rate of first order beam to be only O( ε). These beams are based on same high frequency approximation as geometrical optics which has O(ε) accuracy. Numerical experiments suggests a better rate for odd order beams For the Helmholtz case we have proved [Motamed, OR] that u(x) u GB(x) O(ε K /2 ) for the Taylor expansion part of the error away from caustics. This gives O(ε) for first order beams. More error cancellations coming in for odd order beams? ( no gain in using even order beams) Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
116 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
117 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
118 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
119 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
120 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
121 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
122 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
123 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
124 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
125 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
126 Numerical examples Cusp caustic Consider the test case where Φ(0, y) = y 1 + y 2 2, A(0, y) = e 10 y 2. Cusp caustic at t = 0.5 Two fold caustics at t > 0.5 Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
127 Numerical examples Cusp caustic, convergence Energy norm Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
128 Numerical examples Cusp caustic, convergence Max norm Olof Runborg (KTH) Gaussian Beam Approximation IPP Garching, / 36
Gaussian Beam Approximations
Gaussian Beam Approximations Olof Runborg CSC, KTH Joint with Hailiang Liu, Iowa, and Nick Tanushev, Austin Technischen Universität München München, January 2011 Olof Runborg (KTH) Gaussian Beam Approximations
More informationERROR ESTIMATES FOR GAUSSIAN BEAM SUPERPOSITIONS August 6, 2010
ERROR ESTIMATES FOR GAUSSIAN BEAM SUPERPOSITIONS August 6, 00 HAILIANG LIU, OLOF RUNBORG, AND NICOLAY M. TANUSHEV Abstract. Gaussian beams are asymptotically valid high frequency solutions to hyperbolic
More informationComputational High Frequency Wave Propagation
Computational High Frequency Wave Propagation Olof Runborg CSC, KTH Isaac Newton Institute Cambridge, February 2007 Olof Runborg (KTH) High-Frequency Waves INI, 2007 1 / 52 Outline 1 Introduction, background
More informationGaussian 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 informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationGAUSSIAN BEAM METHODS FOR STRICTLY HYPERBOLIC SYSTEMS
GAUSSIAN BEAM METHODS FOR STRICTLY HYPERBOLIC SYSTEMS HAILIANG LIU AND MAKSYM PRYPOROV Abstract. In this work we construct Gaussian beam approximations to solutions of the strictly hyperbolic system with
More informationFROZEN GAUSSIAN APPROXIMATION FOR HIGH FREQUENCY WAVE PROPAGATION
FROZEN GAUSSIAN APPROXIMATION FOR HIGH FREQUENCY WAVE PROPAGATION JIANFENG LU AND XU YANG Abstract. We propose the frozen Gaussian approimation for computation of high frequency wave propagation. This
More informationGaussian beam methods for the Schrödinger equation with discontinuous potentials
Noname manuscript No. (will be inserted by the editor) Gaussian beam methods for the Schrödinger equation with discontinuous potentials Shi Jin Dongming Wei Dongsheng Yin Received: date / Accepted: date
More informationSUMMARY INTRODUCTION THEORY
Seismic modeling using the frozen Gaussian approximation Xu Yang, University of California Santa Barbara, Jianfeng Lu, Duke University, and Sergey Fomel, University of Texas at Austin SUMMARY We adopt
More informationGAUSSIAN BEAM METHODS FOR THE SCHRÖDINGER EQUATION WITH DISCONTINUOUS POTENTIALS
GAUSSIAN BEAM METHODS FOR THE SCHRÖDINGER EQUATION WITH DISCONTINUOUS POTENTIALS SHI JIN, DONGMING WEI, AND DONGSHENG YIN Abstract. We propose Eulerian and Lagrangian Gaussian beam methods for the Schrödinger
More informationA Survey of Computational High Frequency Wave Propagation I
A Survey of Computational High Frequency Wave Propagation I Bjorn Engquist University of Texas at Austin CSCAMM Workshop on High Frequency Wave Propagation, University of Maryland, September 19-22, 2005
More informationComputational Methods of Physical Problems with Mixed Scales
Computational Methods of Physical Problems with Mixed Scales Shi Jin University of Wisconsin-Madison IPAM Tutorial, March 10-12, 2008 Scales and Physical Laws Figure from E & Engquist, AMS Notice Connections
More informationStochastic Regularity of a Quadratic Observable of High Frequency Waves
Research in the Mathematical Sciences manuscript No. (will be inserted by the editor) Stochastic Regularity of a Quadratic Observable of High Frequency Waves G. Malenová M. Motamed O. Runborg Received:
More informationGaussian Beam Decomposition of High Frequency Wave Fields
Gaussian Beam Decomposition of High Frequency Wave Fields Nicolay M. Tanushev, Björn Engquist, Richard Tsai Department of Mathematics The University of Texas at Austin Austin, TX 787 April, 9 Abstract
More informationIn a uniform 3D medium, we have seen that the acoustic Green s function (propagator) is
Chapter Geometrical optics The material in this chapter is not needed for SAR or CT, but it is foundational for seismic imaging. For simplicity, in this chapter we study the variable-wave speed wave equation
More informationEulerian semiclassical computational methods in quantum dynamics. Shi Jin Department of Mathematics University of Wisconsin-Madison
Eulerian semiclassical computational methods in quantum dynamics Shi Jin Department of Mathematics University of Wisconsin-Madison outline Classical mechanics: Hamiltonian system, discontinuous Hamiltonian,
More informationSemi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime
Commun. Comput. Phys. doi:.48/cicp.99.63s Vol. 9, No. 3, pp. 668-687 March Semi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime Shi Jin,, Hao Wu and
More informationSolving Wave Equation Using Asymptotic Approaches
Solving Wave Equation Using Asymptotic Approaches Thesis by: Mi Yu June 4, 2018 Abstract Unlike ordinary differential equations, partial differential equations have more than one independent variables,
More informationEULERIAN GAUSSIAN BEAM METHOD FOR HIGH FREQUENCY WAVE PROPAGATION IN THE REDUCED MOMENTUM SPACE
EULERIAN GAUSSIAN BEAM METHOD FOR HIGH FREQUENCY WAVE PROPAGATION IN THE REDUCED MOMENTUM SPACE HAO WU AND XU YANG Abstract. Eulerian Gaussian beam method is an efficient wa to compute high frequenc wave
More informationc 2003 International Press
COMM. MATH. SCI. Vol., No. 3, pp. 593 6 c 3 International Press COMPUTATIONAL HIGH-FREQUENCY WAVE PROPAGATION USING THE LEVEL SET METHOD, WITH APPLICATIONS TO THE SEMI-CLASSICAL LIMIT OF SCHRÖDINGER EQUATIONS
More informationLocal smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping
Local smoothing and Strichartz estimates for manifolds with degenerate hyperbolic trapping H. Christianson partly joint work with J. Wunsch (Northwestern) Department of Mathematics University of North
More informationSemi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime
Semi-Eulerian and High Order Gaussian Beam Methods for the Schrödinger Equation in the Semiclassical Regime Shi Jin, Hao Wu, and Xu Yang October 9, 29 Abstract A novel Eulerian Gaussian beam method was
More informationA LEVEL SET FRAMEWORK FOR TRACKING MULTI-VALUED SOLUTIONS OF NONLINEAR FIRST-ORDER EQUATIONS
A LEVEL SET FRAMEWORK FOR TRACKING MULTI-VALUED SOLUTIONS OF NONLINEAR FIRST-ORDER EQUATIONS HAILIANG LIU, LI-TIEN CHENG, AND STANLEY OSHER Abstract. We introduce a level set method for the computation
More information1. Introduction We are interested in the Gaussian beam methods for the numerical approximation of the Schrödinger equation:
COMMUN. MATH. SCI. Vol. 6, No. 4, pp. 995 2 c 28 International Press GAUSSIAN BEAM METHODS FOR THE SCHRÖDINGER EQUATION IN THE SEMI-CLASSICAL REGIME: LAGRANGIAN AND EULERIAN FORMULATIONS SHI JIN, HAO WU,
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationA LEVEL SET FRAMEWORK FOR CAPTURING MULTI-VALUED SOLUTIONS OF NONLINEAR FIRST-ORDER EQUATIONS
A LEVEL SET FRAMEWORK FOR CAPTURING MULTI-VALUED SOLUTIONS OF NONLINEAR FIRST-ORDER EQUATIONS HAILIANG LIU, LI-TIEN CHENG, AND STANLEY OSHER Abstract. We introduce a level set method for the computation
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationFull-wave Electromagnetic Field Simulations in the Lower Hybrid Range of Frequencies
Full-wave Electromagnetic Field Simulations in the Lower Hybrid Range of Frequencies P.T. Bonoli, J.C. Wright, M. Porkolab, PSFC, MIT M. Brambilla, IPP, Garching, Germany E. D Azevedo, ORNL, Oak Ridge,
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More informationMicrolocal 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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationCHAPTER 6 Quantum Mechanics II
CHAPTER 6 Quantum Mechanics II 6.1 6.2 6.3 6.4 6.5 6.6 6.7 The Schrödinger Wave Equation Expectation Values Infinite Square-Well Potential Finite Square-Well Potential Three-Dimensional Infinite-Potential
More informationMATH 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 informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationGAUGE-INVARIANT FROZEN GAUSSIAN APPROXIMATION METHOD FOR THE SCHRÖDINGER EQUATION WITH PERIODIC POTENTIALS
GAUGE-INVARIANT FROZEN GAUSSIAN APPROXIMATION METHOD FOR THE SCHRÖDINGER EQUATION WITH PERIODIC POTENTIALS RICARDO DELGADILLO, JIANFENG LU, AND XU YANG Abstract. We develop a gauge-invariant frozen Gaussian
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
More informationNonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation
Nonlinear Solitary Wave Basis Functions for Long Time Solution of Wave Equation John Steinhoff Subhashini Chitta Nonlinear Waves --Theory and Applications Beijing - June, 2008 Objective Simulate short
More informationEvolution of semiclassical Wigner function (the higher dimensio
Evolution of semiclassical Wigner function (the higher dimensional case) Workshop on Fast Computations in Phase Space, WPI-Vienna, November 2008 Dept. Appl. Math., Univ. Crete & IACM-FORTH 1 2 3 4 5 6
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationODE Final exam - Solutions
ODE Final exam - Solutions May 3, 018 1 Computational questions (30 For all the following ODE s with given initial condition, find the expression of the solution as a function of the time variable t You
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationThe X-ray transform for a non-abelian connection in two dimensions
The X-ray transform for a non-abelian connection in two dimensions David Finch Department of Mathematics Oregon State University Corvallis, OR, 97331, USA Gunther Uhlmann Department of Mathematics University
More informationCan There Be a General Theory of Fourier Integral Operators?
Can There Be a General Theory of Fourier Integral Operators? Allan Greenleaf University of Rochester Conference on Inverse Problems in Honor of Gunther Uhlmann UC, Irvine June 21, 2012 How I started working
More informationCHAPTER 6 Quantum Mechanics II
CHAPTER 6 Quantum Mechanics II 6.1 The Schrödinger Wave Equation 6.2 Expectation Values 6.3 Infinite Square-Well Potential 6.4 Finite Square-Well Potential 6.5 Three-Dimensional Infinite-Potential Well
More information1. Introduction We are interested in developing the Gaussian beam method for the Dirac equation in the semi-classical regime
GAUSSIAN BEAM METHODS FOR THE DIRAC EQUATION IN THE SEMI-CLASSICAL REGIME HAO WU, ZHONGYI HUANG, SHI JIN, AND DONGSHENG YIN Abstract. The Dirac equation is an important model in relativistic quantum mechanics.
More informationInverse problems for hyperbolic PDEs
Inverse problems for hyperbolic PDEs Lauri Oksanen University College London Example: inverse problem for the wave equation Let c be a smooth function on Ω R n and consider the wave equation t 2 u c 2
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationQuasi-neutral limit for Euler-Poisson system in the presence of plasma sheaths
in the presence of plasma sheaths Department of Mathematical Sciences Ulsan National Institute of Science and Technology (UNIST) joint work with Masahiro Suzuki (Nagoya) and Chang-Yeol Jung (Ulsan) The
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationThe wave function for a particle of energy E moving in a constant potential V
Chapter 37 WKB quantization The wave function for a particle of energy E moving in a constant potential V is ψ = Ae i pq (37.1) with a constant amplitude A, and constant wavelength λ = 2π/k, k = p/, and
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationFast Huygens sweeping methods for Helmholtz equations in inhomogeneous media in the high frequency regime
Fast Huygens sweeping methods for Helmholtz equations in inhomogeneous media in the high frequency regime Songting Luo Jianliang Qian Robert Burridge December 2, 23 Abstract In some applications, it is
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationz x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.
Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These
More informationComputation of the Semiclassical Limit of the Schrödinger Equation with Phase Shift by a Level Set Method
Computation of the Semiclassical Limit of the Schrödinger Equation with Phase Shift by a Level Set Method Shi Jin and Xu Yang Abstract In this paper, we show how the level set method, developed in [8,
More information4. Supplementary Notes on Time and Space Evolution of a Neutrino Beam
Lecture Notes for Quantum Physics II & III 8.05 & 8.059 Academic Year 1996/1997 4. Supplementary Notes on Time and Space Evolution of a Neutrino Beam c D. Stelitano 1996 As an example of a two-state system
More informationIntroduction to Wavelet Based Numerical Homogenization
Introduction to Wavelet Based Numerical Homogenization Olof Runborg Numerical Analysis, School of Computer Science and Communication, KTH RTG Summer School on Multiscale Modeling and Analysis University
More information1. Introduction. We are interested in developing efficient numerical methods for computing the Schrödinger equation on unbounded domain
COPUTATION OF THE SCHRÖDINGER EQUATION IN THE SEICLASSICAL REGIE ON UNBOUNDED DOAIN XU YANG AND JIWEI ZHANG Abstract. The study of this paper is two-fold: On the one hand, we generalize the high-order
More informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationWeak Feller Property and Invariant Measures
Weak Feller Property and Invariant Measures Martin Ondreját, J. Seidler, Z. Brzeźniak Institute of Information Theory and Automation Academy of Sciences Prague September 11, 2012 Outline 2010: Stochastic
More informationTheory 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 informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationStrichartz Estimates in Domains
Department of Mathematics Johns Hopkins University April 15, 2010 Wave equation on Riemannian manifold (M, g) Cauchy problem: 2 t u(t, x) gu(t, x) =0 u(0, x) =f (x), t u(0, x) =g(x) Strichartz estimates:
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationTopics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality
Topics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality Po-Lam Yung The Chinese University of Hong Kong Introduction While multiplier operators are very useful in studying
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationMICROLOCAL 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 informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationLECTURE NOTES ON GEOMETRIC OPTICS
LECTURE NOTES ON GEOMETRIC OPTICS PLAMEN STEFANOV We want to solve the wave equation 1. The wave equation (1.1) ( 2 t c 2 g )u = 0, u t=0 = f 1, u t t=0 = f 2, in R t R n x, where g is a Riemannian metric
More informationSoliton-like Solutions to NLS on Compact Manifolds
Soliton-like Solutions to NLS on Compact Manifolds H. Christianson joint works with J. Marzuola (UNC), and with P. Albin (Jussieu and UIUC), J. Marzuola, and L. Thomann (Nantes) Department of Mathematics
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationQuantum Theory. Thornton and Rex, Ch. 6
Quantum Theory Thornton and Rex, Ch. 6 Matter can behave like waves. 1) What is the wave equation? 2) How do we interpret the wave function y(x,t)? Light Waves Plane wave: y(x,t) = A cos(kx-wt) wave (w,k)
More informationGaussian beams and the propagation of singularities
Gaussian beams and the propagation of singularities Dominic Dold, Karen Habermann, Ellen Powell March 17, 2014 Contents 1 Introduction 3 2 Geometrical optics 5 2.1 Construction of the approximation...............................
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationPARTIAL DIFFERENTIAL EQUATIONS
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: P.A. Markowich bgmt2008 http://www.damtp.cam.ac.uk/group/apde/teaching/pde_partii.html In addition to the sets of lecture notes written by previous lecturers ([1],[2])
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationWebster s horn model on Bernoulli flow
Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular
More informationRecovering point sources in unknown environment with differential data
Recovering point sources in unknown environment with differential data Yimin Zhong University of Texas at Austin April 29, 2013 Introduction Introduction The field u satisfies u(x) + k 2 (1 + n(x))u(x)
More informationCLASSIFICATION 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 information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationGaussian integrals and Feynman diagrams. February 28
Gaussian integrals and Feynman diagrams February 28 Introduction A mathematician is one to whom the equality e x2 2 dx = 2π is as obvious as that twice two makes four is to you. Lord W.T. Kelvin to his
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationDispersive Equations and Hyperbolic Orbits
Dispersive Equations and Hyperbolic Orbits H. Christianson Department of Mathematics University of California, Berkeley 4/16/07 The Johns Hopkins University Outline 1 Introduction 3 Applications 2 Main
More informationt L(q, q)dt (11.1) ,..., S ) + S q 1
Chapter 11 WKB and the path integral In this chapter we discuss two reformulations of the Schrödinger equations that can be used to study the transition from quantum mechanics to classical mechanics. They
More information