Computational High Frequency Wave Propagation
|
|
- Margery McBride
- 5 years ago
- Views:
Transcription
1 Computational High Frequency Wave Propagation Olof Runborg CSC, KTH Isaac Newton Institute Cambridge, February 2007 Olof Runborg (KTH) High-Frequency Waves INI, / 52
2 Outline 1 Introduction, background 2 Geometrical optics (GO) 3 Formulations of GO and numerical methods 1 Eikonal equations 2 Ray tracing 3 Wavefront methods 4 Kinetic formulation 5 Phase space methods 4 Other models Olof Runborg (KTH) High-Frequency Waves INI, / 52
3 Wave equation Cauchy problem for scalar wave equation u tt c(x) 2 u = 0, (t, x) R + R d, u(0, x) = A(x)e iωφ(x ), u t (0, x) = ωb(x)e iωφ(x ), where c(x) (variable) speed of propagation. Helmholtz equation (u = v exp(iωt)) v + n(x) 2 ω 2 v = f (x), +radiation condition, where n(x) := 1/c(x) index of refraction, ω angular frequency. Scattering problem Similar versions for elastic wave equation, Maxwell equations, Schrödinger equation,.... Olof Runborg (KTH) High-Frequency Waves INI, / 52
4 Computational challenge Simulation at high frequencies a major challenge! I.e., high relative to size of computational domain in time and space, ω 1 T, ω c = k 1 X. (When 0 t T and 0 x X.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
5 Computational challenge Simulation at high frequencies a major challenge! I.e., high relative to size of computational domain in time and space, ω 1 T, ω c = k 1 X. (When 0 t T and 0 x X.) High frequency short wave length highly oscillatory solutions many gridpoints. Olof Runborg (KTH) High-Frequency Waves INI, / 52
6 Computational challenge Simulation at high frequencies a major challenge! I.e., high relative to size of computational domain in time and space, ω 1 T, ω c = k 1 X. (When 0 t T and 0 x X.) High frequency short wave length highly oscillatory solutions many gridpoints. Direct numerical solution resolves wavelength: #gridpoints ω d at least. Olof Runborg (KTH) High-Frequency Waves INI, / 52
7 Computational challenge Simulation at high frequencies a major challenge! I.e., high relative to size of computational domain in time and space, ω 1 T, ω c = k 1 X. (When 0 t T and 0 x X.) High frequency short wave length highly oscillatory solutions many gridpoints. Direct numerical solution resolves wavelength: #gridpoints ω d at least. Often unrealistic approach for applications in e.g. optics, electromagnetics, geophysics, acoustics,... Olof Runborg (KTH) High-Frequency Waves INI, / 52
8 Computational challenge cont. Challenge met by Efficient direct numerical methods (cost increases with ω for fixed accuracy) High frequency approximations (accuracy increases with ω for fixed cost) Prescribed tolerance methods for certain problems, [Bruno, Chandler-Wilde,... ] (cost independent of ω for fixed accuracy) Olof Runborg (KTH) High-Frequency Waves INI, / 52
9 Geometrical optics Helmholtz equation u + ω 2 n(x) 2 u = 0. Write solution on the form u(x) = A(x, ω)e iωφ(x). Olof Runborg (KTH) High-Frequency Waves INI, / 52
10 Geometrical optics Helmholtz equation u + ω 2 n(x) 2 u = 0. Write solution on the form u(x) = A(x, ω)e iωφ(x). (a) Amplitude A(x) (b) Phase φ(x) Olof Runborg (KTH) High-Frequency Waves INI, / 52
11 Geometrical optics Amplitude A and phase φ vary on a much coarser scale than u. (And varies little with ω.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
12 Geometrical optics Amplitude A and phase φ vary on a much coarser scale than u. (And varies little with ω.) Geometrical optics approximation considers A and φ as ω. Olof Runborg (KTH) High-Frequency Waves INI, / 52
13 Geometrical optics Amplitude A and phase φ vary on a much coarser scale than u. (And varies little with ω.) Geometrical optics approximation considers A and φ as ω. Good accuracy for large ω. Computational cost ω-independent. u(x) = A(x)e iωφ(x) + O(1/ω). Olof Runborg (KTH) High-Frequency Waves INI, / 52
14 Geometrical optics Amplitude A and phase φ vary on a much coarser scale than u. (And varies little with ω.) Geometrical optics approximation considers A and φ as ω. Good accuracy for large ω. Computational cost ω-independent. u(x) = A(x)e iωφ(x) + O(1/ω). Waves propagate as rays, c.f. visible light. Not all wave effects captured correctly (diffraction, caustics). Olof Runborg (KTH) High-Frequency Waves INI, / 52
15 Geometrical optics Amplitude A and phase φ vary on a much coarser scale than u. (And varies little with ω.) Geometrical optics approximation considers A and φ as ω. Good accuracy for large ω. Computational cost ω-independent. u(x) = A(x)e iωφ(x) + O(1/ω). Waves propagate as rays, c.f. visible light. Not all wave effects captured correctly (diffraction, caustics). More generally, multiple crossing waves: u(x) = N A n (x)e iωφn(x) + O(1/ω). n=1 Several amplitude and phase functions. Olof Runborg (KTH) High-Frequency Waves INI, / 52
16 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). Olof Runborg (KTH) High-Frequency Waves INI, / 52
17 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Olof Runborg (KTH) High-Frequency Waves INI, / 52
18 Limitations of geometrical optics Diffracted waves Olof Runborg (KTH) High-Frequency Waves INI, / 52
19 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Olof Runborg (KTH) High-Frequency Waves INI, / 52
20 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Creeping rays Olof Runborg (KTH) High-Frequency Waves INI, / 52
21 Limitations of geometrical optics Creeping rays Olof Runborg (KTH) High-Frequency Waves INI, / 52
22 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Creeping rays Olof Runborg (KTH) High-Frequency Waves INI, / 52
23 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Creeping rays Disappear in GO limit, but is of size 1/ω α, 0 < α < 1. (Hence 1/ω.) GTD used in numerical simulations by e.g. adding diffracted waves at corners. Olof Runborg (KTH) High-Frequency Waves INI, / 52
24 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Creeping rays Disappear in GO limit, but is of size 1/ω α, 0 < α < 1. (Hence 1/ω.) GTD used in numerical simulations by e.g. adding diffracted waves at corners. Caustics is another problem: Amplitude A infinite, but should be A ω α, 0 < α < 1. Olof Runborg (KTH) High-Frequency Waves INI, / 52
25 Limitations of geometrical optics Caustics Concentration of rays. A when ω. Olof Runborg (KTH) High-Frequency Waves INI, / 52
26 Limitations of geometrical optics GO approximation breaks down when the wavelength variations in c(x). By the geometrical theory of diffraction (GTD) ([Keller, 59]) errors can be much larger than O(1/ω) because of boundaries: Diffracted waves Creeping rays Disappear in GO limit, but is of size 1/ω α, 0 < α < 1. (Hence 1/ω.) GTD used in numerical simulations by e.g. adding diffracted waves at corners. Caustics is another problem: Amplitude A infinite, but should be A ω α, 0 < α < 1. Olof Runborg (KTH) High-Frequency Waves INI, / 52
27 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
28 Numerical and modeling issues Goal: Find phase (traveltime) and amplitude in a domain, on a regular grid. Issues: Multiple arrivals, crossing rays, superposition multivaluedness. Eulerian/Lagrangian model. Fixed/moving grids. PDE/ODE model. Complexity. Olof Runborg (KTH) High-Frequency Waves INI, / 52
29 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
30 Eikonal- and transport equations Derivation Make the WKB ansatz u(x) = e iωφ(x) k=0 a k (x)(iω) k. Olof Runborg (KTH) High-Frequency Waves INI, / 52
31 Eikonal- and transport equations Derivation Make the WKB ansatz u(x) = e iωφ(x) k=0 a k (x)(iω) k. Enter this into Helmholtz equation u + ω 2 n 2 u = 0, ω 2 (n 2 φ 2 )u + iω(2 φ a 0 + φa 0 )e iωφ(x) + O(ω 0 ) = 0. Olof Runborg (KTH) High-Frequency Waves INI, / 52
32 Eikonal- and transport equations Derivation Make the WKB ansatz u(x) = e iωφ(x) k=0 a k (x)(iω) k. Enter this into Helmholtz equation u + ω 2 n 2 u = 0, ω 2 (n 2 φ 2 )u + iω(2 φ a 0 + φa 0 )e iωφ(x) + O(ω 0 ) = 0. Equate coefficients of powers of ω to zero: Olof Runborg (KTH) High-Frequency Waves INI, / 52
33 Eikonal- and transport equations Derivation Make the WKB ansatz u(x) = e iωφ(x) k=0 a k (x)(iω) k. Enter this into Helmholtz equation u + ω 2 n 2 u = 0, ω 2 (n 2 φ 2 )u + iω(2 φ a 0 + φa 0 )e iωφ(x) + O(ω 0 ) = 0. Equate coefficients of powers of ω to zero: For ω 2 we get the eikonal equation φ = n(x), Olof Runborg (KTH) High-Frequency Waves INI, / 52
34 Eikonal- and transport equations Derivation Make the WKB ansatz u(x) = e iωφ(x) k=0 a k (x)(iω) k. Enter this into Helmholtz equation u + ω 2 n 2 u = 0, ω 2 (n 2 φ 2 )u + iω(2 φ a 0 + φa 0 )e iωφ(x) + O(ω 0 ) = 0. Equate coefficients of powers of ω to zero: For ω 2 we get the eikonal equation φ = n(x), For ω 1 we get the transport equation 2 φ a 0 + a 0 φ = 0. Olof Runborg (KTH) High-Frequency Waves INI, / 52
35 Eikonal- and transport equations Derivation Make the WKB ansatz u(x) = e iωφ(x) k=0 a k (x)(iω) k. Enter this into Helmholtz equation u + ω 2 n 2 u = 0, ω 2 (n 2 φ 2 )u + iω(2 φ a 0 + φa 0 )e iωφ(x) + O(ω 0 ) = 0. Equate coefficients of powers of ω to zero: For ω 2 we get the eikonal equation φ = n(x), For ω 1 we get the transport equation 2 φ a 0 + a 0 φ = 0. Discard rest of terms in series (ω ). Let A = a 0. Olof Runborg (KTH) High-Frequency Waves INI, / 52
36 Eikonal equation Hamilton Jacobi type nonlinear PDE. Can be solved efficiently on fixed Eulerian grids. Olof Runborg (KTH) High-Frequency Waves INI, / 52
37 Eikonal equation Hamilton Jacobi type nonlinear PDE. Can be solved efficiently on fixed Eulerian grids. Stationary version. φ = n(x). Fast marching [Sethian,Tsitsiklis] or fast sweeping methods [Zhao, Tsai, et al]. Olof Runborg (KTH) High-Frequency Waves INI, / 52
38 Eikonal equation Hamilton Jacobi type nonlinear PDE. Can be solved efficiently on fixed Eulerian grids. Stationary version. φ = n(x). Fast marching [Sethian,Tsitsiklis] or fast sweeping methods [Zhao, Tsai, et al]. Time-dependent version. Wave equation plus ansatz u(t, x) A(t, x)e iωϕ(t,x) give ϕ t + n(x) 1 ϕ = 0. Upwind, high-resolution (ENO, WENO) finite difference methods [Osher, Shu, et al] Olof Runborg (KTH) High-Frequency Waves INI, / 52
39 Eikonal equation Hamilton Jacobi type nonlinear PDE. Can be solved efficiently on fixed Eulerian grids. Stationary version. φ = n(x). Fast marching [Sethian,Tsitsiklis] or fast sweeping methods [Zhao, Tsai, et al]. Time-dependent version. Wave equation plus ansatz u(t, x) A(t, x)e iωϕ(t,x) give ϕ t + n(x) 1 ϕ = 0. Upwind, high-resolution (ENO, WENO) finite difference methods [Osher, Shu, et al] (Note, if IC and BC match, ϕ = φ t.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
40 Eikonal equation Properties Eikonal equation: φ = n(x), Olof Runborg (KTH) High-Frequency Waves INI, / 52
41 Eikonal equation Properties Eikonal equation: φ = n(x), Ansatz only treats one wave. In general crossing waves u(x) A 1 (x)e iωφ 1(x) + A 2 (x)e iωφ 2(x) + Olof Runborg (KTH) High-Frequency Waves INI, / 52
42 Eikonal equation Properties Eikonal equation: φ = n(x), Ansatz only treats one wave. In general crossing waves u(x) A 1 (x)e iωφ 1(x) + A 2 (x)e iωφ 2(x) + Nonlinear equation, no superposition principle Olof Runborg (KTH) High-Frequency Waves INI, / 52
43 Eikonal equation Properties Eikonal equation: φ = n(x), Ansatz only treats one wave. In general crossing waves u(x) A 1 (x)e iωφ1(x) + A 2 (x)e iωφ2(x) + Nonlinear equation, no superposition principle Viscosity solution, kinks Olof Runborg (KTH) High-Frequency Waves INI, / 52
44 Eikonal equation Properties Eikonal equation: φ = n(x), Ansatz only treats one wave. In general crossing waves u(x) A 1 (x)e iωφ1(x) + A 2 (x)e iωφ2(x) + Nonlinear equation, no superposition principle Viscosity solution, kinks First arrival property: φ visc (x) = min n φ n (x) Olof Runborg (KTH) High-Frequency Waves INI, / 52
45 Eikonal equation First arrival property (a) Correct solution (b) Eikonal equation Note: φ traveltime of the wave φ(x) =constant (level sets) represent wave fronts Olof Runborg (KTH) High-Frequency Waves INI, / 52
46 Eikonal equation Example Olof Runborg (KTH) High-Frequency Waves INI, / 52
47 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
48 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
49 Ray tracing Rays are the (bi)characteristics (x(t), p(t)) of the eikonal equation, given by ODEs dx dt = c(x) 2 p, dp dt = c(x) c(x), p(t) is local ray direction, "slowness" vector. Olof Runborg (KTH) High-Frequency Waves INI, / 52
50 Ray tracing Rays are the (bi)characteristics (x(t), p(t)) of the eikonal equation, given by ODEs dx dt = c(x) 2 p, dp dt = c(x) c(x), p(t) is local ray direction, "slowness" vector. Ray tracing method of characteristics for eikonal equation. Olof Runborg (KTH) High-Frequency Waves INI, / 52
51 Ray tracing Rays are the (bi)characteristics (x(t), p(t)) of the eikonal equation, given by ODEs dx dt = c(x) 2 p, dp dt = c(x) c(x), p(t) is local ray direction, "slowness" vector. Ray tracing method of characteristics for eikonal equation. If valid at t = 0, then for all t > 0: φ(x(t)) = p(t), (local ray direction wavefronts) Olof Runborg (KTH) High-Frequency Waves INI, / 52
52 Ray tracing Rays are the (bi)characteristics (x(t), p(t)) of the eikonal equation, given by ODEs dx dt = c(x) 2 p, dp dt = c(x) c(x), p(t) is local ray direction, "slowness" vector. Ray tracing method of characteristics for eikonal equation. If valid at t = 0, then for all t > 0: φ(x(t)) = p(t), (local ray direction wavefronts) φ(x(t)) = t, (phase traveltime) Olof Runborg (KTH) High-Frequency Waves INI, / 52
53 Ray tracing Rays are the (bi)characteristics (x(t), p(t)) of the eikonal equation, given by ODEs dx dt = c(x) 2 p, dp dt = c(x) c(x), p(t) is local ray direction, "slowness" vector. Ray tracing method of characteristics for eikonal equation. If valid at t = 0, then for all t > 0: φ(x(t)) = p(t), (local ray direction wavefronts) φ(x(t)) = t, (phase traveltime) p(t) = 1/c(x(t)), (can reduce to p S d 1, in 2D: p θ) Olof Runborg (KTH) High-Frequency Waves INI, / 52
54 Ray tracing Rays are the (bi)characteristics (x(t), p(t)) of the eikonal equation, given by ODEs dx dt = c(x) 2 p, dp dt = c(x) c(x), p(t) is local ray direction, "slowness" vector. Ray tracing method of characteristics for eikonal equation. If valid at t = 0, then for all t > 0: φ(x(t)) = p(t), (local ray direction wavefronts) φ(x(t)) = t, (phase traveltime) p(t) = 1/c(x(t)), (can reduce to p S d 1, in 2D: p θ) There are also ODEs for the amplitude along rays. Olof Runborg (KTH) High-Frequency Waves INI, / 52
55 Ray tracing Numerics Ray equations (with p = 1/c) Methods: dx dt = c(x) 2 p, dp dt = c(x) c(x), Olof Runborg (KTH) High-Frequency Waves INI, / 52
56 Ray tracing Numerics Ray equations (with p = 1/c) dx dt = c(x) 2 p, dp dt = c(x) c(x), Methods: Standard numerical ODE methods, e.g. Runge Kutta. Olof Runborg (KTH) High-Frequency Waves INI, / 52
57 Ray tracing Numerics Ray equations (with p = 1/c) dx dt = c(x) 2 p, dp dt = c(x) c(x), Methods: Standard numerical ODE methods, e.g. Runge Kutta. If c(x) is piecewise constant rays piecewise straight lines refracted/reflected at interfaces by Snell s law. (Geometrical problem, shooting and bouncing rays (SBR).) Olof Runborg (KTH) High-Frequency Waves INI, / 52
58 Ray tracing Numerics Ray equations (with p = 1/c) dx dt = c(x) 2 p, dp dt = c(x) c(x), Methods: Standard numerical ODE methods, e.g. Runge Kutta. If c(x) is piecewise constant rays piecewise straight lines refracted/reflected at interfaces by Snell s law. (Geometrical problem, shooting and bouncing rays (SBR).) Properties: Rays can cross. (No problem with multivaluedness.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
59 Ray tracing Numerics Ray equations (with p = 1/c) dx dt = c(x) 2 p, dp dt = c(x) c(x), Methods: Standard numerical ODE methods, e.g. Runge Kutta. If c(x) is piecewise constant rays piecewise straight lines refracted/reflected at interfaces by Snell s law. (Geometrical problem, shooting and bouncing rays (SBR).) Properties: Rays can cross. (No problem with multivaluedness.) Solution only given along rays. (Hard to interpolate traveltime to regular grid.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
60 Ray tracing Numerics Ray equations (with p = 1/c) dx dt = c(x) 2 p, dp dt = c(x) c(x), Methods: Standard numerical ODE methods, e.g. Runge Kutta. If c(x) is piecewise constant rays piecewise straight lines refracted/reflected at interfaces by Snell s law. (Geometrical problem, shooting and bouncing rays (SBR).) Properties: Rays can cross. (No problem with multivaluedness.) Solution only given along rays. (Hard to interpolate traveltime to regular grid.) Diverging rays. (Hard to cover full domain.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
61 Example: Ray tracing solution Olof Runborg (KTH) High-Frequency Waves INI, / 52
62 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
63 Wavefronts Wavefront given by φ =constant. Suppose initial wavefront is γ(α) and φ(γ) = 0. Let x(t, α), p(t, α) be the rays emitted orthogonal from this curve: x(0, α) = γ(α), p(0, α) = γ (α) c γ (α). Wavefronts also given by: x(t = constant, α). More general definition. Olof Runborg (KTH) High-Frequency Waves INI, / 52
64 Phase space Wavefront in xy-space typically non-smooth. Caustics appear. Olof Runborg (KTH) High-Frequency Waves INI, / 52
65 Phase space Wavefront in xy-space typically non-smooth. Caustics appear. Introduce phase space (x, p), where p S d 1 is local ray direction Olof Runborg (KTH) High-Frequency Waves INI, / 52
66 Phase space Wavefront in xy-space typically non-smooth. Caustics appear. Introduce phase space (x, p), where p S d 1 is local ray direction Observation: Wavefront is a smooth curve in phase space. 2D problems: 1D curve in 3D phase space (x, y, θ). 3D problems: 2D surface in 5D phase space (x, y, z, θ, α). Olof Runborg (KTH) High-Frequency Waves INI, / 52
67 Phase space Wavefront in xy-space typically non-smooth. Caustics appear. Introduce phase space (x, p), where p S d 1 is local ray direction Observation: Wavefront is a smooth curve in phase space. 2D problems: 1D curve in 3D phase space (x, y, θ). 3D problems: 2D surface in 5D phase space (x, y, z, θ, α). Wavefront in phase space sweeps out a smooth surface the Lagrangian submanifold. Olof Runborg (KTH) High-Frequency Waves INI, / 52
68 Wavefront tracking Directly solve for wavefront given by x(t = const, α). Suppose γ(α) is the initial wavefront, x(0, α) = γ(α). Follow ensemble of rays x(t, α) = c 2 p, t x(0, α) = γ(α), p(t, α) = c t c, p(0, α) = γ (α) c γ (α). Note: In principle we do not need to track p. Moving front in normal direction a possibility x t = c x α x α (since 0 = αφ(x(t, α)) = x α φ = x α p) But not good numerically since wavefront non-smooth! Olof Runborg (KTH) High-Frequency Waves INI, / 52
69 Wavefront construction [Vinje, Iversen, Gjöystdal, Lambaré,... ] Solve for x(t, α) and p(t, α). Discretize in α and trace rays for α 1, α 2, α 3,... where α j = j α. x(t n, α j ) Olof Runborg (KTH) High-Frequency Waves INI, / 52
70 Wavefront construction [Vinje, Iversen, Gjöystdal, Lambaré,... ] Solve for x(t, α) and p(t, α). Discretize in α and trace rays for α 1, α 2, α 3,... where α j = j α. Insert new rays adaptively by interpolation when front resolution deteriorates. E.g.: If x(t n, α j+1 ) x(t n, α j ) tol then insert new ray at α j+1/2. x(t n, α j ) Olof Runborg (KTH) High-Frequency Waves INI, / 52
71 Wavefront construction [Vinje, Iversen, Gjöystdal, Lambaré,... ] Solve for x(t, α) and p(t, α). Discretize in α and trace rays for α1, α2, α3,... where αj = j α. Insert new rays adaptively by interpolation when front resolution deteriorates. E.g.: If x(tn, αj+1 ) x(tn, αj ) tol then insert new ray at αj+1/2. Interpolate traveltime/phase/amplitude onto regular grid. x(tn, αj ) Olof Runborg (KTH) High-Frequency Waves INI, / 52
72 Example: Wavefront tracking solution Multiple arrivals ok. Lagrangian method. Interpolation can be complicated. Olof Runborg (KTH) High-Frequency Waves INI, / 52
73 Wavefront tracking, Eulerian versions Level set methods [Osher, Tsai, Cheng, Liu, Jin, Qian,... ] Level set methods: Represent wavefront implicitly as the zero level set of a signed distance function φ(x). γ = {x R d : φ(x) = 0}. Olof Runborg (KTH) High-Frequency Waves INI, / 52
74 Level set methods Wavefront in phase space higher co-dimension: Represent as intersection of zero level sets of several functions. In 2D, γ = {x R 3 : φ 1 (x) = 0, φ 2 (x) = 0}, φ j = φ j (x, y, θ). Level set functions satisfy linear hyperbolic PDEs: t φ j + c(x, y) cos θ x φ j + c(x, y) sin θ y φ j + (c x sin θ c y cos θ) θ φ j = 0. Use local level set method to reduce complexity. Olof Runborg (KTH) High-Frequency Waves INI, / 52
75 Wavefront tracking, Eulerian versions Segment projection method [Engquist, Tornberg, OR] Wavefront γ is given as a union of curve segments γ j. The segments are chosen such that they can be represented by a function of one coordinate variable: f i (t, x) and g j (t, y) From ray equations simple PDEs can be derived for f i and g j. Ex: t f i + u x f i = v. Connectivity of segments is also maintained. (a) Interface (b) x-segments (c) y-segments Olof Runborg (KTH) High-Frequency Waves INI, / 52
76 Segment projection method Example Plane wave entering from the left. (a) Local Ray Directions (b) Wave Fronts Olof Runborg (KTH) High-Frequency Waves INI, / 52
77 Segment projection method Example, segments Olof Runborg (KTH) High-Frequency Waves INI, / 52
78 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
79 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
80 Kinetic formulation View rays as trajectories of particles. Let f (t, x, p) be the particle (photon) density in phase space. Olof Runborg (KTH) High-Frequency Waves INI, / 52
81 Kinetic formulation View rays as trajectories of particles. Let f (t, x, p) be the particle (photon) density in phase space. Bicharacteristic equations Liouville equation f t + c 2 p x f c c pf = 0. Olof Runborg (KTH) High-Frequency Waves INI, / 52
82 Kinetic formulation View rays as trajectories of particles. Let f (t, x, p) be the particle (photon) density in phase space. Bicharacteristic equations Liouville equation f t + c 2 p x f c c pf = 0. A ray at x in direction p 0 (t, x) and amplitude A(t, x) represented as f (t, x, p) = A 2 (t, x)δ(p p 0 (t, x)). Olof Runborg (KTH) High-Frequency Waves INI, / 52
83 Kinetic formulation View rays as trajectories of particles. Let f (t, x, p) be the particle (photon) density in phase space. Bicharacteristic equations Liouville equation f t + c 2 p x f c c pf = 0. A ray at x in direction p 0 (t, x) and amplitude A(t, x) represented as f (t, x, p) = A 2 (t, x)δ(p p 0 (t, x)). Since p = n(x), the density f supported on sphere p = n(x). Olof Runborg (KTH) High-Frequency Waves INI, / 52
84 Kinetic formulation Relation to Wigner theory Liouville equation can also be derived directly from wave eq. through e.g. Wigner measures [Tartar, Lions, Paul, Gerard, Mauser, Markowich, Poupaud,... ] Olof Runborg (KTH) High-Frequency Waves INI, / 52
85 Kinetic formulation Relation to Wigner theory Liouville equation can also be derived directly from wave eq. through e.g. Wigner measures [Tartar, Lions, Paul, Gerard, Mauser, Markowich, Poupaud,... ] Let f be limit of the Wigner tranform of Helmholtz solution f = lim ω F y pu(x + y/2ω)u(x y/2ω). Then f satisfies Liouville eq. f t + c 2 p x f c c pf = 0. Olof Runborg (KTH) High-Frequency Waves INI, / 52
86 Kinetic formulation Relation to Wigner theory Liouville equation can also be derived directly from wave eq. through e.g. Wigner measures [Tartar, Lions, Paul, Gerard, Mauser, Markowich, Poupaud,... ] Let f be limit of the Wigner tranform of Helmholtz solution f = lim ω F y pu(x + y/2ω)u(x y/2ω). Then f satisfies Liouville eq. f t + c 2 p x f c c pf = 0. Relationship to wave equation solution: u = Ae iωφ f = A 2 δ (p φ). Olof Runborg (KTH) High-Frequency Waves INI, / 52
87 Kinetic formulation Relation to Wigner theory Liouville equation can also be derived directly from wave eq. through e.g. Wigner measures [Tartar, Lions, Paul, Gerard, Mauser, Markowich, Poupaud,... ] Let f be limit of the Wigner tranform of Helmholtz solution f = lim ω F y pu(x + y/2ω)u(x y/2ω). Then f satisfies Liouville eq. f t + c 2 p x f c c pf = 0. Relationship to wave equation solution: u = Ae iωφ f = A 2 δ (p φ). lim F y pa(x + y/2ω)a(x y/2ω)e iω(φ(x+y/2ω) φ(x y/2ω)) ω = F y p A 2 (x)e iy φ = A 2 δ (p φ). Olof Runborg (KTH) High-Frequency Waves INI, / 52
88 Kinetic formulation Relation to Wigner theory Liouville equation can also be derived directly from wave eq. through e.g. Wigner measures [Tartar, Lions, Paul, Gerard, Mauser, Markowich, Poupaud,... ] Let f be limit of the Wigner tranform of Helmholtz solution f = lim ω F y pu(x + y/2ω)u(x y/2ω). Then f satisfies Liouville eq. f t + c 2 p x f c c pf = 0. Relationship to wave equation solution: u = Ae iωφ f = A 2 δ (p φ). Note: Loss of phase information. Olof Runborg (KTH) High-Frequency Waves INI, / 52
89 Moment equations Derived from Liouville equation in phase space + closure assumption for a system of equations representing the moments. (C.f. hydrodynamic limit from Boltzmann eq.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
90 Moment equations Derived from Liouville equation in phase space + closure assumption for a system of equations representing the moments. (C.f. hydrodynamic limit from Boltzmann eq.) Full equation for f expensive to solve numerically. 6 independent variables in 3D. Olof Runborg (KTH) High-Frequency Waves INI, / 52
91 Moment equations Derived from Liouville equation in phase space + closure assumption for a system of equations representing the moments. (C.f. hydrodynamic limit from Boltzmann eq.) Full equation for f expensive to solve numerically. 6 independent variables in 3D. Moment eq is a PDE description in the small (t, x)-space. Fixed Eulerian grids can be used. Olof Runborg (KTH) High-Frequency Waves INI, / 52
92 Moment equations Derived from Liouville equation in phase space + closure assumption for a system of equations representing the moments. (C.f. hydrodynamic limit from Boltzmann eq.) Full equation for f expensive to solve numerically. 6 independent variables in 3D. Moment eq is a PDE description in the small (t, x)-space. Fixed Eulerian grids can be used. Arbitrary good superposition. N crossing waves allowed. (But larger N means a larger system of PDEs must be solved.) [Brenier, Corrias, Engquist, OR] (wave equation), [Gosse, Jin, Li, Markowich, Sparber] (Schrödinger) Olof Runborg (KTH) High-Frequency Waves INI, / 52
93 Moment equations Derivations, homogenous case (c 1) Starting point is f t + p x f = 0. Let p = (p 1, p 2 ). Define the moments, m ij = p1 i pj R 2 2 fdp. From p1 i pj R 2 2 (f t + p x f )dp = 0, we get the infinite (valid i, j 0) system of moment equations (m ij ) t + (m i+1,j ) x + (m i,j+1 ) y = 0. Olof Runborg (KTH) High-Frequency Waves INI, / 52
94 Moment equations Derivations, homogenous case, cont. Make the closure assumption N f (x, p, t) = A 2 k δ(p p k ), k=1 ( ) cos p k = θk. sin θ k The moments take the form m ij = N A 2 k cosi θ k sin j θ k. k=1 Corresponds to a maximum of N waves at each point. Choose equations for moments m 2k 1,0 and m 0,2k 1, k = 1,..., N. Gives closed system of 2N equations with 2N unknowns (the A k s and θ k s). Olof Runborg (KTH) High-Frequency Waves INI, / 52
95 Moment equations Examples Ex. N = 1 ( u1 u 2 ) t + u 1 2 u 2 1 +u2 2 u 1 u 2 u 2 1 +u2 2 x + u 1 u 2 u 2 1 +u2 2 u 2 2 u 2 1 +u2 2 y = 0. where u 1 = m 10 = A 2 cos θ and u 2 = m 01 = A 2 sin θ. For N 2, F 0 (u) t + F 1 (u) x + F 2 (u) y = 0. where F 0 (u), F 1 (u) and F 2 (u) are complicated non-linear functions. PDE = weakly hyperbolic system of conservation laws, (with source terms when c varies) Flux functions in conservation law can be difficult to evaluate. Olof Runborg (KTH) High-Frequency Waves INI, / 52
96 Moment equations Wedge example (a) N = 1 (b) N = 2 Olof Runborg (KTH) High-Frequency Waves INI, / 52
97 Geometrical optics models and numerical methods u + ω 2 n(x) 2 u = 0 Rays dx dt = c 2 p, dp dt Ray tracing = c c Wavefront methods Kinetic f t + c 2 p x f 1 c c pf = 0 Moment methods, Phase space methods Eikonal φ = n(x) Hamilton Jacobi methods Olof Runborg (KTH) High-Frequency Waves INI, / 52
98 Phase Space Methods Eikonal solvers, wavefront tracking, etc. essentially constructed to solve the problem for one single initial data, e.g. a point source or a plane wave. Solves for a "sheet" in phase space (the Lagrangian submanifold). Some applications demand solution for many different sets of initial data (e.g. inverse problems). Phase space methods: Get solution in the whole phase space. (Corresponds to solving for all possible initial data.) Olof Runborg (KTH) High-Frequency Waves INI, / 52
99 Phase Space Methods, cont. Equation to solve typically of the type time-independent Liouville equation/whole level set equation. (But other interpretation of unknown.) Computational cost: Suppose N discretization points in each dimension. In 2D eikonal solvers/wavefront construction typically cost O(N 2 ) for one initial data. Phase space methods cost O(N 3 log N) for all initial data. Combining advantages of using a fixed grid (cf eikonal solvers) and capturing multiple arrivals (cf ray methods) relatively easy. Examples: Phase Flow Method [Candès, Ying], Fast phase space method [Fomel, Sethian]. Olof Runborg (KTH) High-Frequency Waves INI, / 52
100 Physical optics Integral formulation of scattering problem u s (x) = G( x x ) u(x ) n dx, x R d \ Ω, Ω Approximate u n by geometrical optics solution. E.g. if u inc = exp(iωα x) is a plane wave, Ω convex, then { u n (u inc + us GO ) 2iωα ˆn(x)e iωα x, x illuminated, = n 0, x in shadow. Cost of evaluating solution still depends on ω. Olof Runborg (KTH) High-Frequency Waves INI, / 52
101 Physical optics Exact physical optics For convex Ω, make ansatz u n = iωa(x, ω)eiωα x then A(x, ω) smooth, uniformly in ω, except at shadow boundaries. Discretize and solve A(x, ω) at cost independent of ω. [Bruno, Chandler Wilde,... ] Olof Runborg (KTH) High-Frequency Waves INI, / 52
102 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Olof Runborg (KTH) High-Frequency Waves INI, / 52
103 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Error normally O(1/ω). Boundaries and caustic cause bigger errors. GTD reduces them. Olof Runborg (KTH) High-Frequency Waves INI, / 52
104 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Error normally O(1/ω). Boundaries and caustic cause bigger errors. GTD reduces them. Trade-off point in ω increases with increasing computer power and is problem and accuracy dependent. Olof Runborg (KTH) High-Frequency Waves INI, / 52
105 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Error normally O(1/ω). Boundaries and caustic cause bigger errors. GTD reduces them. Trade-off point in ω increases with increasing computer power and is problem and accuracy dependent. Challenges for geometrical optics methods include capturing multiple arrivals on fixed grids at reasonable complexity. Olof Runborg (KTH) High-Frequency Waves INI, / 52
106 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Error normally O(1/ω). Boundaries and caustic cause bigger errors. GTD reduces them. Trade-off point in ω increases with increasing computer power and is problem and accuracy dependent. Challenges for geometrical optics methods include capturing multiple arrivals on fixed grids at reasonable complexity. Phase space methods useful when many similar problems need to be solved. Olof Runborg (KTH) High-Frequency Waves INI, / 52
107 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Error normally O(1/ω). Boundaries and caustic cause bigger errors. GTD reduces them. Trade-off point in ω increases with increasing computer power and is problem and accuracy dependent. Challenges for geometrical optics methods include capturing multiple arrivals on fixed grids at reasonable complexity. Phase space methods useful when many similar problems need to be solved. Hybrid methods coupling elements of asymptotic models and direct methods a possibility. Olof Runborg (KTH) High-Frequency Waves INI, / 52
108 Concluding remarks Complexity of simulation based on high-frequency approximations are not ω-dependent. Error normally O(1/ω). Boundaries and caustic cause bigger errors. GTD reduces them. Trade-off point in ω increases with increasing computer power and is problem and accuracy dependent. Challenges for geometrical optics methods include capturing multiple arrivals on fixed grids at reasonable complexity. Phase space methods useful when many similar problems need to be solved. Hybrid methods coupling elements of asymptotic models and direct methods a possibility. Methods generalize to Maxwell, elastic wave eq, etc. Olof Runborg (KTH) High-Frequency Waves INI, / 52
109 References B. Engquist and O. Runborg. Computational high-frequency wave propagation. Acta Numerica, 12: , J.-D. Benamou. An introduction to Eulerian geometrical optics ( ). J. Sci. Comput., 19(1-3):63 93, O. Runborg. Mathematical models and numerical methods for high frequency waves. To appear in Commun. Comp. Phys., Olof Runborg (KTH) High-Frequency Waves INI, / 52
A 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 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 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 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 informationGaussian 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 informationGaussian Beam Approximations of High Frequency Waves
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
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 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 informationA LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS
A LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS SHI JIN AND STANLEY J. OSHER Abstract. We develop a level set method for the
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 informationRecent computational methods for high frequency waves in heterogeneous media ABSTRACT. Shi Jin. 1. Introduction
Recent computational methods for high frequency waves in heterogeneous media Shi Jin ABSTRACT In this note, we review our recent results on the Eulerian computation of high frequency waves in heterogeneous
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 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 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 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 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 informationSemiclassical computational methods for quantum dynamics with bandcrossings. Shi Jin University of Wisconsin-Madison
Semiclassical computational methods for quantum dynamics with bandcrossings and uncertainty Shi Jin University of Wisconsin-Madison collaborators Nicolas Crouseilles, Rennes Mohammed Lemou, Rennes Liu
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationNumerical methods for hyperbolic systems with singular coefficients: well-balanced scheme, Hamiltonian preservation, and beyond
Proceedings of Symposia in Applied Mathematics Numerical methods for hyperbolic systems with singular coefficients: well-balanced scheme, Hamiltonian preservation, and beyond Shi Jin Abstract. This paper
More informationDevelopment of level set methods for computing the semiclassical limit of Schrödinger equations with potentials
Retrospective Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 28 Development of level set methods for computing the semiclassical limit of Schrödinger equations with
More informationA Numerical Approach To The Steepest Descent Method
A Numerical Approach To The Steepest Descent Method K.U. Leuven, Division of Numerical and Applied Mathematics Isaac Newton Institute, Cambridge, Feb 15, 27 Outline 1 One-dimensional oscillatory integrals
More informationA Fast Method for Solving the Helmholtz Equation Based on Wave Splitting JELENA POPOVIC
A Fast Method for Solving the Helmholtz Equation Based on Wave Splitting JELENA POPOVIC Licenciate Thesis Stockholm, Sweden 2009 TRITA-CSC-A 2009:11 ISSN-1653-5723 ISRN KTH/CSC/A- -09/11-SE ISBN 978-91-7415-370-5
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationS. Marmorat (POEMS - INRIA Rocquencourt) F. Collino (CERFACS) J.-D. Benamou (POEMS - INRIA Rocquencourt)
Numerical MicroLocal Analysis (NMLA) S. Marmorat (POEMS - INRIA Rocquencourt) F. Collino (CERFACS) J.-D. Benamou (POEMS - INRIA Rocquencourt) Given frequency domain wave data, the proposed algorithm gives
More informationFast Hybrid Algorithms for High Frequency Scattering
Fast Hybrid Algorithms for High Frequency Scattering Björn Engquist, Khoa Tran, and Lexing Ying Citation: AIP Conference Proceedings 1106, 3 (2009); doi: 10.1063/1.3117111 View online: http://dx.doi.org/10.1063/1.3117111
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 informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationSeismology and Seismic Imaging
Seismology and Seismic Imaging 4. Ray theory N. Rawlinson Research School of Earth Sciences, ANU Seismology lecture course p.1/23 The ray approximation Here, we consider the problem of how body waves (P
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
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 informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More information1. Introduction. In this paper, we construct and study a numerical scheme for the Liouville equation in d-dimension:
SIAM J. NUMER. ANAL. Vol. 44 No. 5 pp. 8 88 c 6 Society for Industrial and Applied Mathematics A HAMILTONIAN-PRESERVING SCHEME FOR THE LIOUVILLE EQUATION OF GEOMETRICAL OPTICS WITH PARTIAL TRANSMISSIONS
More informationSeismic imaging and optimal transport
Seismic imaging and optimal transport Bjorn Engquist In collaboration with Brittany Froese, Sergey Fomel and Yunan Yang Brenier60, Calculus of Variations and Optimal Transportation, Paris, January 10-13,
More informationIntroduction to Multiscale Modeling
Introduction to Multiscale Modeling Bjorn Engquist University of Texas at Austin Summer School on Multiscale Modeling and Analysis, University of Texas at Austin, July 21 - August 8, 2008 Outline I Overview
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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationRadiative transfer equations with varying refractive index: a mathematical perspective
Radiative transfer equations with varying refractive index: a mathematical perspective Guillaume Bal November 23, 2005 Abstract This paper reviews established mathematical techniques to model the energy
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 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 informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationLecture 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 informationBasic principles of the seismic method
Chapter 2 Basic principles of the seismic method In this chapter we introduce the basic notion of seismic waves. In the earth, seismic waves can propagate as longitudinal (P) or as shear (S) waves. For
More informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
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 informationLecture 15 Notes: 07 / 26. The photoelectric effect and the particle nature of light
Lecture 15 Notes: 07 / 26 The photoelectric effect and the particle nature of light When diffraction of light was discovered, it was assumed that light was purely a wave phenomenon, since waves, but not
More informationWave propagation in discrete heterogeneous media
www.bcamath.org Wave propagation in discrete heterogeneous media Aurora MARICA marica@bcamath.org BCAM - Basque Center for Applied Mathematics Derio, Basque Country, Spain Summer school & workshop: PDEs,
More informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
More informationA STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS *
Journal of Computational Mathematics Vol.8, No.4, 010, 55 568. http://www.global-sci.org/jcm doi:10.408/jcm.1003-m0016 A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
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 informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationModule 2: First-Order Partial Differential Equations
Module 2: First-Order Partial Differential Equations The mathematical formulations of many problems in science and engineering reduce to study of first-order PDEs. For instance, the study of first-order
More informationA New Finite-Difference Method for General Long-Range Rotorcraft Acoustics: Initial Comparisons with Intermediate- Range Data
A New Finite-Difference Method for General Long-Range Rotorcraft Acoustics: Initial Comparisons with Intermediate- Range Data Subhashini Chitta, John Steinhoff (B.H. Goethert Professor, UTSI), Andrew Wilson
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 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 informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationIntroduction to PML in time domain
Introduction to PML in time domain Alexander Thomann Introduction to PML in time domain - Alexander Thomann p.1 Overview 1 Introduction 2 PML in one dimension Classical absorbing layers One-dimensional
More informationLecture notes 5: Diffraction
Lecture notes 5: Diffraction Let us now consider how light reacts to being confined to a given aperture. The resolution of an aperture is restricted due to the wave nature of light: as light passes through
More informationMathematical and computational methods for semiclassical Schrödinger equations
Acta Numerica (2012), pp. 1 89 c Cambridge University Press, 2012 DOI: 10.1017/S0962492902 Printed in the United Kingdom Mathematical and computational methods for semiclassical Schrödinger equations Shi
More informationOptimal transport for Seismic Imaging
Optimal transport for Seismic Imaging Bjorn Engquist In collaboration with Brittany Froese and Yunan Yang ICERM Workshop - Recent Advances in Seismic Modeling and Inversion: From Analysis to Applications,
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
More informationAnalysis of a Fast Method for Solving the High Frequency Helmholtz Equation in One Dimension
BIT manuscript No. will be inserted by the editor Analysis of a Fast Method for Solving the High Frequency Helmholtz Equation in One Dimension Jelena Popovic Olof Runborg Received: date / Accepted: date
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
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 informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationDiffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch)
Diffraction by Edges András Vasy (with Richard Melrose and Jared Wunsch) Cambridge, July 2006 Consider the wave equation Pu = 0, Pu = D 2 t u gu, on manifolds with corners M; here g 0 the Laplacian, D
More informationEulerian calculations of wave breaking and multi-valued solutions in a traveling wave tube
Eulerian calculations of wave breaking and multi-valued solutions in a traveling wave tube J.G. Wöhlbier, S. Jin, 2 and S. Sengele 3 Los Alamos National Laboratory, ISR-6, Radio Frequency and Accelerator
More informationLecture 11: Polarized Light. Fundamentals of Polarized Light. Descriptions of Polarized Light. Scattering Polarization. Zeeman Effect.
Lecture 11: Polarized Light Outline 1 Fundamentals of Polarized Light 2 Descriptions of Polarized Light 3 Scattering Polarization 4 Zeeman Effect 5 Hanle Effect Fundamentals of Polarized Light Electromagnetic
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationWaves, Polarization, and Coherence
Waves, Polarization, and Coherence Lectures 5 Biophotonics Jae Gwan Kim jaekim@gist.ac.kr, X 2220 School of Information and Communication Engineering Gwangju Institute of Sciences and Technology Outline
More informationConceptual Physics. Luis A. Anchordoqui. Department of Physics and Astronomy Lehman College, City University of New York. Lesson VI October 3, 2017
Conceptual Physics Luis A. Anchordoqui Department of Physics and Astronomy Lehman College, City University of New York Lesson VI October 3, 2017 https://arxiv.org/abs/1711.07445 L. A. Anchordoqui (CUNY)
More informationLecture 16. Wednesday, May 25, φ t + V φ = 0 (1) In upwinding, we choose the discretization of the spatial derivative based on the sign of u.
Lecture 16 Wednesday, May 25, 2005 Consider the linear advection equation φ t + V φ = 0 (1) or in 1D φ t + uφ x = 0. In upwinding, we choose the discretization of the spatial derivative based on the sign
More informationNanoscale shift of the intensity distribution of dipole radiation
Shu et al. Vol. 26, No. 2/ February 2009/ J. Opt. Soc. Am. A 395 Nanoscale shift of the intensity distribution of dipole radiation Jie Shu, Xin Li, and Henk F. Arnoldus* Department of Physics and Astronomy,
More informationExtended isochron rays in prestack depth (map) migration
Extended isochron rays in prestack depth (map) migration A.A. Duchkov and M.V. de Hoop Purdue University, 150 N.University st., West Lafayette, IN, 47907 e-mail: aduchkov@purdue.edu (December 15, 2008)
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
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 informationLecture 19 Optical MEMS (1)
EEL6935 Advanced MEMS (Spring 5) Instructor: Dr. Huikai Xie Lecture 19 Optical MEMS (1) Agenda: Optics Review EEL6935 Advanced MEMS 5 H. Xie 3/8/5 1 Optics Review Nature of Light Reflection and Refraction
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 informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationLight.notebook May 03, 2016
Unit 4 Light LIGHT.1 Describe the ray model of light. 16.1 LIGHT.2 Predict the effect of distance on light s illuminance. 16.1 LIGHT.3 Explain polarization and the Doppler effect. 16.2 LIGHT.4 Describe
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 informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics Diffraction I Basic Physics M.P. Vaughan Diffraction Electromagnetic waves Geometric wavefront The Principle of Linear Superposition Diffraction regimes Single
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 informationNumerical dispersion and Linearized Saint-Venant Equations
Numerical dispersion and Linearized Saint-Venant Equations M. Ersoy Basque Center for Applied Mathematics 11 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 The Saint-Venant equations
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
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 informationElectricity & Optics
Physics 24100 Electricity & Optics Lecture 26 Chapter 33 sec. 1-4 Fall 2017 Semester Professor Koltick Interference of Light Interference phenomena are a consequence of the wave-like nature of light Electric
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationThe Fast Sweeping Method. Hongkai Zhao. Department of Mathematics University of California, Irvine
The Fast Sweeping Method Hongkai Zhao Department of Mathematics University of California, Irvine Highlights of fast sweeping mehtod Simple (A Gauss-Seidel type of iterative method). Optimal complexity.
More informationLecture 11: Introduction to diffraction of light
Lecture 11: Introduction to diffraction of light Diffraction of waves in everyday life and applications Diffraction in everyday life Diffraction in applications Spectroscopy: physics, chemistry, medicine,
More informationHOMOGENIZATION OF METRIC HAMILTON-JACOBI EQUATIONS
HOMOGENIZATION OF METRIC HAMILTON-JACOBI EQUATIONS ADAM M. OBERMAN, RYO TAKEI, AND ALEXANDER VLADIMIRSKY Abstract. In this work we provide a novel approach to homogenization for a class of static Hamilton-Jacobi
More informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.306 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Monday March 12, 2018. Turn it in (by 3PM) at the Math.
More information