Gaussian beams for high frequency waves: some recent developments
|
|
- Kelley Gallagher
- 5 years ago
- Views:
Transcription
1 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 Waves Sanya, China. July, 2.
2 Outline High frequency waves and geometrical optics 2 Single-scale Gaussian beams for Schrödinger Beam ingredients Beam width Gaussian wavepacket transforms Numerical results 3 Multiscale Gaussian beams for wave equations 4 Eulerian Gaussian beams 5 Conclusion
3 The parabolic wave equation (Schrödinger) ı U t + H(, ı )U =, R d, t >, U t= = U (). H(, p) is a real principal symbol: H(, p) = V() + p 2 2. is a small positive parameter; V is a smooth potential. U is a compactly supported L 2 -function, presumably highly oscillatory. real(u)
4 An eample real(u) Figure: Schrödinger under a potential hill. = T = π. Solution at
5 The problem The difficulties When the initial condition contains oscillations of a small wavelength, the Schrödinger equation propagates these oscillations in space and time. Resolving such small oscillations by direct numerical methods requires an enormous computational grid and is very costly in practice. Methods based on geometrical optics are sought as an alternative for capturing such highly oscillatory phenomena.
6 Geometrical optics The traditional geometrical-optics method starts from the WKBJ ansatz consisting of a real phase and a real amplitude, U(, t) A(, t)e ı τ(,t). The WKBJ ansatz breaks down at the caustics, in the neighborhood of which the phase function is multivalued and the amplitude function blows up.
7 High frequency waves and geometrical optics Single-scale Gaussian beams for Schrödinger Multiscale Gaussian beams for wave e t t Geometrical optics for Schrödinger (Cont.)
8 t Gaussian beams Gaussian beam (Babich 72, Hörmander 7, Cerveny-Popov-Psencik 82, Popov 82, Ralston 83, Tanushev-Q-Ralston 7) constructs a global comple phase and a global comple amplitude that satisfy the eikonal equation and the transport equation approimately near a specified ray path. Away from the ray path, the quadratic imaginary part of the phase provides a rapidly decaying Gaussian profile. This gives rise to a single Gaussian-like asymptotic solution, which is accurate near the underlying ray path. Since the equation is linear, a superposition yields a global asymptotic solution
9 Beam ingredients Eikonal and transport equations ı U t + H(, ı )U =, R d, t >, () U(, ) = U () (2) Asymptotic solutions: A(, t)e ı τ(,t), so that the equation () and its associated initial conditions (2) are satisfied approimately with a small error when is small. Consider the leading orders in powers of : τ t (, t) + H(, τ (, t)) =, (3) A t (, t) + H p A + A 2 (trace(h ppτ + H p )) =. (4)
10 Beam ingredients A single Gaussian beam () A single Gaussian beam is concentrated near a ray path which is the -projection of a certain bicharacteristic. The beam ingredients are computed by solving: ẋ = d = H p ((t), p(t)), dt t= =, ṗ = dp = H ((t), p(t)), dt p t= = p, dτ dt = p(t) H p ((t), p(t)) H((t), p(t)), τ t= = τ, dm(t) + HpM T + MH p + MH pp M + H =, M t= = M = ıǫi dt da dt + A 2 (trace(h ppm + H p )) =, A t= = A, τ(t) = τ((t), t), p(t) = τ ((t), t), and M(t) = τ ((t), t).
11 Beam ingredients A single Gaussian beam (2) We define two global, smooth approimate functions for the phase and amplitude along the ray path γ: τ(, t) τ(t) + p(t) ( (t)) + 2 ( (t))t M(t)( (t)) A(, t) A((t), t) = A(t), which are accurate near the ray path γ = {((t), t) : t }.
12 Beam ingredients Taylor epansion for beams τ(t)=τ(t,(t)) p(t)=τ (t,(t)) γ={(t,(t))} ((t),p(t),m(t),a(t)) M(t)=τ (t,(t)).5 τ(t,)=τ(t,(t))+τ (t,(t)) ( (t))+.5τ (t,(t)) ( (t)) 2 A(t,)=A(t)=A(t,(t)) t
13 t Beam ingredients A single Gaussian beam (3) We have a single-beam asymptotic solution concentrated on γ = {((t), t) : t }, Φ(, t) = A(, t) ep(ı τ(, t) ). (5) Because Im(τ(, t)) = 2 ( (t))t Im(M(t))( (t)), Φ(, t) has a Gaussian profile concentrated on γ: ( ep ) 2 ( (t))t Im(M(t))( (t))
14 Beam width Beam width and long-term beam propagation Beam width and reinitialization A simple analysis reveals that the width of a Gaussian beam in the case of the Schrödinger equation is related to the derivatives of the underlying potentials, and a beam width may grow eponentially during the evolution process. So the beam loses its localized significance, leading to deteriorating accuracy in the Taylor epansion for the phase function and high cost in the beam summation.
15 Beam width Beam width and Hessian of potential A single Gaussian beam centered at a ray trajectory {((t), t) : t } propagates according to a quadratic approimation to potential centered at = (t) for t : V() = V((t))+V ((t))( (t))+ 2 ( (t))t V ((t))( (t) The properties of V ((t)) control how the beam width grows. Assume that the quadratic potential function is given at : V() = V( ) + V ( )( ) + 2 ( ) t V ( )( ). The Riccati equation for the Hessian of the phase function can be written as: Ṁ = M 2 V ( ), M t= = ıǫi, (6) where ǫ is a given positive number, and I is the identity.
16 Beam width Hessian of potential: V ( ). Solving the resulting Riccati equation for M yields M(t) = tǫ2 + ıǫ + t 2 ǫ 2. Thus the beam decays like ep ( 2 ) ( Im(M(t))( (t))2 = ep ) ǫ 2 + t 2 ǫ2( (t))2 (+t The beam width is roughly proportional to 2 ǫ 2 ) ǫ, which grows almost linearly as time t increases.
17 Beam width Hessian of potential: V ( ) >. Solving the Riccati equation for M yields (ω = V ( )) M(t) = ω(ǫ2 ω 2 ) sin ωt cos ωt + ıǫω 2 ω 2 cos 2 ωt + ǫ 2 sin 2. ωt Thus the beam decays like ( ) ǫω 2 ep 2 (ω 2 cos 2 ωt + ǫ 2 sin 2 ( (t))2. ωt) The beam width is roughly proportional to (ω 2 cos 2 ωt + ǫ 2 sin 2 ωt) ǫω 2, which is bounded from above by ǫ if ǫ ω or by ǫ if ω 2 ǫ ω. In particular, when ǫ = ω, the beam width does not change.
18 Beam width Hessian of potential: V ( ) < Solving the Riccati equation for M yields (ω = V ( )) M(t) = ǫ2 ω( + )( e 4ωt ) + ı4ǫe 2ωt ω 2 ( + e 2ωt ) 2 + ǫ2 ( e ω 2ωt ). 2 2 The beam decays very slowly like ( ep 4ǫe 2ωt ( 2 [( + e 2ωt ) 2 + ǫ2 ( e ω 2ωt ) 2 (t))2 ] 2 The beam width is roughly proportional to e ωt [ω 2 ( + e 2ωt ) 2 + ǫ 2 ( e 2ωt ) 2 ] 4ǫω 2, which grows eponentially as t increases. ).
19 Beam width Hessian of potential: V ( ) < (Cont.) This implies that the beam loses its localized significance, leading to deteriorating accuracy in the Taylor epansion for the phase function and high cost in beam summation. In practice, along a ray trajectory the Hessian of a potential function may change from positive definite to negative definite, so that the beam width may grow eponentially unepectedly. In the multi-dimensional case, the growth of the beam width depends on the Hessian of the potential in more complicated ways. However, some general features will be the same as in the one-dimensional case.
20 Beam width Long time propagation by reinitialization One natural solution to this issue is to monitor the widths of the Gaussian beams and reinitialize the Gaussian beam representation before any one of the beams becomes too wide. Although the reinitialization idea is rather straightforward, it is rather difficult to combine with eisting methods of beam initialization for reasons related to representation and efficiency. The reinitialization idea fits perfectly with the fast Gaussian wavepacket transform based beam algorithm (Q-Ying (JCP )).
21 Beam width Long time propagation by reinitialization: cont. As an eample, we take the case V < to demonstrate the benefit of reinitialization. Suppose that V ( ) = 4π 2, ω = V ( ) = 2π, and the final time T = 2.. Without reinitialization the beam width at the final time T is roughly e ωt e 4π = π 2π If we choose t such that e ω t 2π, and we divide [, T] into subintervals, then the beam width in each T t subinterval is roughly. In this particular eample, t.25, and = 8. T t This simple case analysis demonstrates that it is critical to carry out reinitialization.
22 Beam width Beam ingredients The functions, (t), p(t), τ(t), M(t), A(t), τ(, t), A(, t), and Φ(, t), are uniquely determined by the initial data, p, τ, M, and A. We denote these initial data collectively by a tuple α = (, p, τ, M, A ). The solutions are denoted, respectively, by α (t), p α (t), M α (t), A α (t), τ α (, t), A α (, t), and Φ α (, t).
23 Beam width Beam summation For a given tuple α = (, p, τ, M, A ), the Gaussian beams Φ α (, t) have a simple Gaussian envelope. For a general initial condition U(, ), one needs to find a set I of tuples such that at time t = U(, ) α I Φ α (, ). Once this initial decomposition is given, the linearity of the equation gives the Gaussian beam solution U(, t) α I Φ α (, t).
24 Beam width How to generate a beam decomposition Related works Beam decompositions for wave equations: Cerveny 82, Hill (Geop 9), Leung-Q-Burridge (Geop 7), Tanushev-Q-Ralston (SIAM MMS 7). FBI-transform based beam decompositions for Schrödinger equations: Leung-Q (JCP 9, ). For Schrödinger equations the most efficient beam decomposition method is based on single-scale Gaussian wavepacket transforms: Q-Ying (JCP ).
25 Gaussian wavepacket transforms Motivation At time t = each Gaussian beam takes the form ( ( ı A ep p ( ) + )) 2 ( ) T M ( ), where p = O() and the Hessian M is purely imaginary and of order O(). This is a modulated Gaussian function that oscillates at a wavelength of order O( ) and has an effective support of width O( /2 ) in space. The initialization step is equivalent to decomposing U(, ) into a linear combination of such Gaussian functions. Gaussian wavepacket transforms.
26 Gaussian wavepacket transforms Partitioning the Fourier domain Start by partitioning the Fourier domain R d into d-dimensional boes of size W in each dimension. Denote these boes by B i where i = (i, i 2,...,i d ) is a multiinde with integer components and the center of each B i by ξ i = (ξ i,, ξ i,2,...,ξ i,d ) R d.
27 Gaussian wavepacket transforms Illustrations
28 Gaussian wavepacket transforms Gaussian windows in the Fourier domain Each B i is associated with a smooth window function g i (ξ) that is compactly supported in a bo centered at ξ i with width L = 2W (i.e., d s= [ξ i,s W, ξ i,s + W]). g i (ξ) is also required to approimate a Gaussian profile with σ = W/2. g i (ξ) e ξ ξi 2 σ
29 Gaussian wavepacket transforms A partition of unity Choose g i (ξ) to satisfy three admissible conditions (Q-Ying JCP ). For each B i define the conjugate filter h i (ξ): h i (ξ) = g i (ξ) i (g i(ξ)) 2. By construction, the products of g i (ξ) and h i (ξ) form a partition of unity: i g i(ξ)h i (ξ) =.
30 Gaussian wavepacket transforms Two sets of wavepacket functions We introduce two sets of functions {φ i,k ()}, {ψ i,k ()}, which are defined in the Fourier domain by ˆφ i,k (ξ) = k ξ 2πı e L d/2 L g i (ξ), k Z d, ˆψ i,k (ξ) = k ξ 2πı e L d/2 L h i (ξ), k Z d. Taking the inverse Fourier transforms gives their definitions in the spatial domain: φ i,k () = L d/2 R d e 2πı( k L ) ξ g i (ξ)dξ, k Z d ψ i,k () = L d/2 R d e 2πı( k L ) ξ h i (ξ)dξ, k Z d
31 Gaussian wavepacket transforms Two sets of wavepacket functions (Cont.) The definitions of g i (ξ) and φ i,k () imply that φ i,k () ( ) d π L σ e 2πı( k L ) ξ i e σ2 π 2 k L 2. φ i,k () is approimately a Gaussian function that is spatially centered at k/l, oscillates at frequency ξ i, and has an O(σ) = O( 2) effective width in the Fourier domain and an O(/σ) = O( 2) effective width in the spatial domain. The functions {φ i,k ()} fit eactly into the profile of a Gaussian beam with: ξ i. Qualitatively, ψ i,k () is also a wavepacket with slightly larger support in compared to φ i,k ().
32 Gaussian wavepacket transforms Two sets of wavepacket functions (Cont.) Figure: Typical profiles of φ i,k () and ψ i,k ().
33 Gaussian wavepacket transforms Two frames {φ i,k ()} and {ψ i,k ()} are two frames of L 2 (R d ). Lemma (Q-Ying (JCP )) There eists constants C and C 2 such that for any f L 2 (R d ) C f 2 2 i,k φ i,k, f 2 C 2 f 2 2, C f 2 2 i,k ψ i,k, f 2 C 2 f 2 2.
34 Gaussian wavepacket transforms Dual frames {φ i,k ()} and {ψ i,k ()} are dual frames. Lemma (Q-Ying (JCP )) For any f L 2 (R d ), f() = ψ i,k, f φ i,k (). i,k Lemma 2.2 offers a way to decompose any function f L 2 (R d ) into a sum of Gaussian-like functions.
35 Gaussian wavepacket transforms Forward and inverse Gaussian wavepacket transforms Given f, the forward Gaussian wavepacket transform computes the coefficients {c i,k } defined by c i,k = ψ i,k, f = ˆψ i,k,ˆf (7) where, is the usual L 2 (R d ) inner product and ˆf denotes the Fourier transform of f. By Lemma 2. each coefficient c i,k is at most of order O() and many coefficients are negligible. Given a set of coefficients {c i,k }, the inverse Gaussian wavepacket transform synthesizes a function u() defined by u() = c i,k φ i,k (). (8) i,k
36 Numerical results -D eample: potential hill () Initial condition: ( U(, t = ) = ep 25(.5) 2) ep ( ıτ () ), τ () = ln (ep(5(.5)) + ep( 5(.5))). 5 In this case, we take V() = cos(2π( +.5)), = 256π, η =., N = 24 and T = 2..
37 Numerical results potential hill (2): no reinitialization.5.5 real u real u.5.5 (a) (b) real u.2 real u (c) (d)
38 Numerical results potential hill (3): reinitialization.5.5 real u real u.5.5 (a) (b) real u real u (c) (d)
39 Numerical results potential hill (4): = 2 4 π, T =., N = 27, and Rein = real u real(u) (a) (b) real (u).5.5 position density.5.5 (c) (d)
40 Numerical results 2-D eample V(, y) = +.5 sin(2π) cos(2πy), and the initial condition is given by ( U(t = ) = ep 25(.5) 2) ( ) ıτ () ep, τ () = ln (ep(5(.5)) + ep( 5(.5))). 5 We use N N = = 52π. η =.. We compute the solution up to the final time T =.56 without reinitialization and with reinitialization.
41 y y Numerical results 2-D eample: without reinitialization () (a) (b) real u real u (c) (d)
42 Numerical results 2-D eample: without reinitialization (2) real u real u (e) (f) real u real u (g) y (h) y
43 y y Numerical results 2-D eample: with reinitialization () (a) (b) real u real u (c) (d)
44 Numerical results 2-D eample: with reinitialization (2) real u real u (e) (f) real u real u (g) y (h) y
45 z z y u u u Numerical results 3-D eample: with reinitialization T =.56 (a) y (b) z (c) (d) y (e) (f) y
46 Multiscale Gaussian beams for wave equations U tt V 2 () U =, R d, t >, U t= = f (), U t t= = f 2 (). V is smooth, positive, and bounded away from zero. f () and f 2 () are compactly supported L 2 -functions, presumably highly oscillatory.
47 Multiscale Gaussian-beam setup Beam ingredients can be computed according to the eikonal and transport equations. Two polarized wave modes for the wave equation. For wave equations, the most efficient beam decomposition method is based on fast multiscale Gaussian wavepacket transforms: Q-Ying (SIAM MMS ). This is inspired by the parabolic scaling principle (Smith 98, Candes-Demanet 5).
48 Partitioning the Fourier domain Start by partitioning the Fourier domain R d into Cartesian coronae {C l } for l : C = [ 4, 4] d, C l = {ξ = (ξ, ξ 2,...,ξ d ) : ma s d ξ s [4 l, 4 l ]}, l 2. ξ C l implies that ξ = O(4 l ). Each corona C l is further partitioned into boes B l,i = d [2 l i s, 2 l (i s + )], s= where the integer multiinde i = (i, i 2,...,i d ) ranges over all possible choices that satisfy B l,i C l. All boes in a fied C l have the same length W l = 2 l in each dimension and the center of the bo B l,i is denoted by ξ l,i = (ξ l,i,, ξ l,i,2,...,ξ l,i,d ).
49 Illustrations
50 Gaussian windows in the Fourier domain To each bo B l,i, we associate a smooth function g l,i (ξ), which is compactly supported in a bo centered at ξ l,i with size L l = 2W l in each dimension. g l,i (ξ) is also required to approimate a Gaussian profile with σ l = W l /2. «2 g l,i (ξ) e ξ ξl,i σ l
51 Multiscale Gaussian wavepacket transforms Different from Single-scale Gaussian wavepacket transforms (Q-Ying JCP), we introduce angular resolution as well in Multiscale Gaussian wavepacket transforms (Q-Ying MMS). Define a partition of unity based on these associated Gaussian windows. We introduce two sets of functions {φ l,i,k ()}, {ψ l,i,k ()}, ˆφ l,i,k (ξ) = L d/2 l ˆψ l,i,k (ξ) = L d/2 l k ξ 2πı L e l g l,i (ξ), k Z d, k ξ 2πı L e l h l,i (ξ), k Z d. These define dual frames and consist of multiscale Forward and inverse Gaussian wavepacket transforms.
52 2-D eample V(, y) =. +.5 sin(2π) cos(2πy), and the initial conditions are given by U t= = 2. sin(256π) cos(256πy), U t t= = 2. cos(256π) sin(256πy). (9) We use N N = to discretize [, ] [, ]. η =.5. We compute the solution up to the final time T =.6 without reinitialization and with reinitialization.
53 High frequency waves and geometrical optics Single-scale Gaussian beams for Schrödinger Multiscale Gaussian beams for wave e 2-D eample: without reinitialization () y y
54 2-D eample: without reinitialization (2) u u (c) y (d) y u u 2 (e) y (f) y
55 High frequency waves and geometrical optics Single-scale Gaussian beams for Schrödinger Multiscale Gaussian beams for wave e 2-D eample: with reinitialization () y y
56 2-D eample: with reinitialization (2) u.5.5 u (c) y (d) y u 2 u..2.3 (e) y (f) y
57 Eulerian Gaussian beams: Methodology Based on phase-space PDE formulation of geometrical optics, we can reformulate the ODEs of beam ingredients into PDEs in the full phase space by using the material derivative along ray trajectories (Leung-Q-Burridge Geop7, Leung-Q JCP9, JCP). The initial beam decomposition is achieved by using the Fourier-Biros-Iagnitzer (FBI) transform (Leung-Q JCP9,).
58 Conclusions Two main computational aspects of the Gaussian beam methods for high frequency waves: beam initialization for general initial condition and long time propagation. For the beam initialization problem, we proposed fast (multiscale) Gaussian wavepacket transforms and developed based on them a new efficient algorithm for beam initialization for general initial data. Long time propagation was addressed by reinitializing the beam representation using this algorithm when the beam widths get too wide.
59 Conclusions: cont. Smith 98, Candes-Demanet 5 represented the solution operator of the wave equation in the curvelet frame. However, numerical tests show that these representations often have a significant constant factor, which limits their practical application. We can also view our Gaussian beam method as a way to represent the solution operator of the Schrödinger equation or the wave equation. Since the Gaussian beams (indeed by continuous parameters) are much more fleible and general than the Gaussian wavepackets (indeed by the discrete (i, k) parameters), what we gain in the Gaussian beam method is a representation of the solution operator that is one-to-one and more efficient.
60 Acknowledgement Many thanks to the organizers for the invitation.
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 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 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 informationFAST HUYGENS SWEEPING METHODS FOR SCHRÖDINGER EQUATIONS IN THE SEMI-CLASSICAL REGIME
FAST HUYGENS SWEEPING METHODS FOR SCHRÖDINGER EQUATIONS IN THE SEMI-CLASSICAL REGIME SHINGYU LEUNG, JIANLIANG QIAN, AND SUSANA SERNA Abstract. We propose fast Huygens sweeping methods for Schrödinger equations
More informationc 2014 International Press Vol. 21, No. 1, pp , March
METHODS AND APPLICATIONS OF ANALYSIS. c International Press Vol., No., pp., March FAST HUYGENS SWEEPING METHODS FOR SCHRÖDINGER EQUATIONS IN THE SEMI-CLASSICAL REGIME SHINGYU LEUNG, JIANLIANG QIAN, AND
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 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 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 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 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 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 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 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 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 informationY m = y n e 2πi(m 1)(n 1)/N (1) Y m e 2πi(m 1)(n 1)/N (2) m=1
The Discrete Fourier Transform (Bretherton notes): 1 Definition Let y n, n = 1,..., N be a sequence of N possibly comple values. The discrete Fourier transform () of this sequence is the sequence Y m,
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
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 informationc. Better work with components of slowness vector s (or wave + k z 2 = k 2 = (ωs) 2 = ω 2 /c 2. k=(kx,k z )
.50 Introduction to seismology /3/05 sophie michelet Today s class: ) Eikonal equation (basis of ray theory) ) Boundary conditions (Stein s book.3.0) 3) Snell s law Some remarks on what we discussed last
More informationReproducing formulas associated with symbols
Reproducing formulas associated with symbols Filippo De Mari Ernesto De Vito Università di Genova, Italy Modern Methods of Time-Frequency Analysis II Workshop on Applied Coorbit space theory September
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 5 Fourier Analysis of Finite Difference Methods In this lecture, we determine the stability of PDE discretizations using Fourier analysis. First, we begin with Fourier analysis of PDE s, and then
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationFourier transform of tempered distributions
Fourier transform of tempered distributions 1 Test functions and distributions As we have seen before, many functions are not classical in the sense that they cannot be evaluated at any point. For eample,
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
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 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 informationM445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
More informationA Padé approximation to the scalar wavefield extrapolator for inhomogeneous media
A Padé approimation A Padé approimation to the scalar wavefield etrapolator for inhomogeneous media Yanpeng Mi, Zhengsheng Yao, and Gary F. Margrave ABSTRACT A seismic wavefield at depth z can be obtained
More informationFourier Reconstruction from Non-Uniform Spectral Data
School of Electrical, Computer and Energy Engineering, Arizona State University aditya.v@asu.edu With Profs. Anne Gelb, Doug Cochran and Rosemary Renaut Research supported in part by National Science Foundation
More informationThe Quantum Theory of Atoms and Molecules
The Quantum Theory of Atoms and Molecules The postulates of quantum mechanics Dr Grant Ritchie The postulates.. 1. Associated with any particle moving in a conservative field of force is a wave function,
More informationSignals and Spectra - Review
Signals and Spectra - Review SIGNALS DETERMINISTIC No uncertainty w.r.t. the value of a signal at any time Modeled by mathematical epressions RANDOM some degree of uncertainty before the signal occurs
More informationBabich s Expansion and the Fast Huygens Sweeping Method for the Helmholtz Wave Equation at High Frequencies
Babich s Epansion and the Fast Huygens Sweeping Method for the Helmholt Wave Equation at High Frequencies Wangtao Lu Jianliang Qian Robert Burridge May 4, 25 Abstract In some applications, it is reasonable
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 informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı & Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering August 2 nd, 2010 Outline PART II Pseudodifferential
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Packet The Group Velocity Phys 375 Overview and Motivation: Last time we looked at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
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 informationClassical field theory 2012 (NS-364B) Feynman propagator
Classical field theory 212 (NS-364B Feynman propagator 1. Introduction States in quantum mechanics in Schrödinger picture evolve as ( Ψt = Û(t,t Ψt, Û(t,t = T exp ı t dt Ĥ(t, (1 t where Û(t,t denotes the
More informationSuperposition of electromagnetic waves
Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationGaussian beam diffraction in inhomogeneous media of cylindrical symmetry
Optica Applicata, Vol. XL, No. 3, 00 Gaussian beam diffraction in inhomogeneous media of cylindrical symmetry PAWEŁ BERCZYŃSKI, YURI A. KRAVTSOV, 3, GRZEGORZ ŻEGLIŃSKI 4 Institute of Physics, West Pomeranian
More informationOften, in this class, we will analyze a closed-loop feedback control system, and end up with an equation of the form
ME 32, Spring 25, UC Berkeley, A. Packard 55 7 Review of SLODEs Throughout this section, if y denotes a function (of time, say), then y [k or y (k) denotes the k th derivative of the function y, y [k =
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationModelling in Biology
Modelling in Biology Dr Guy-Bart Stan Department of Bioengineering 17th October 2017 Dr Guy-Bart Stan (Dept. of Bioeng.) Modelling in Biology 17th October 2017 1 / 77 1 Introduction 2 Linear models of
More informationModelling and Specifying Dispersive Laser Cavities
Modelling and Specifying Dispersive Laser Cavities C. Sean Bohun (UOIT), Yuri Cher (Toronto), Linda J. Cummings (NJIT), Mehran Ebrahimi (UOIT), Peter Howell (Oxford), Laurent Monasse (CERMICS-ENPC), Judith
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More information221A Lecture Notes Path Integral
A Lecture Notes Path Integral Feynman s Path Integral Formulation Feynman s formulation of quantum mechanics using the so-called path integral is arguably the most elegant. It can be stated in a single
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 informationEmpirical Wavelet Transform
Jérôme Gilles Department of Mathematics, UCLA jegilles@math.ucla.edu Adaptive Data Analysis and Sparsity Workshop January 31th, 013 Outline Introduction - EMD 1D Empirical Wavelets Definition Experiments
More informationSuperposition of waves (review) PHYS 258 Fourier optics SJSU Spring 2010 Eradat
Superposition of waves (review PHYS 58 Fourier optics SJSU Spring Eradat Superposition of waves Superposition of waves is the common conceptual basis for some optical phenomena such as: Polarization Interference
More informationMathematical Tripos Part III Michaelmas 2017 Distribution Theory & Applications, Example sheet 1 (answers) Dr A.C.L. Ashton
Mathematical Tripos Part III Michaelmas 7 Distribution Theory & Applications, Eample sheet (answers) Dr A.C.L. Ashton Comments and corrections to acla@damtp.cam.ac.uk.. Construct a non-zero element of
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationGaussian beams in inhomogeneous anisotropic layered structures
Gaussian beams in inhomogeneous anisotropic layered structures Vlastislav Červený 1 ) and Ivan Pšenčík 2 ) 1) Charles University, Faculty of Mathematics and Physics, Department of Geophysics, Ke Karlovu
More informationDispersive numerical schemes for Schrödinger equations
Dispersive numerical schemes for Schrödinger equations Enrique Zuazua joint work with L. Ignat zuazua@bcamath.org Basque Center for Applied Mathematics (BCAM), Bilbao, Basque Country, Spain IMA Workshop:
More informationNonlinear Oscillations and Chaos
CHAPTER 4 Nonlinear Oscillations and Chaos 4-. l l = l + d s d d l l = l + d m θ m (a) (b) (c) The unetended length of each spring is, as shown in (a). In order to attach the mass m, each spring must be
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 informationContinuous Wave Data Analysis: Fully Coherent Methods
Continuous Wave Data Analysis: Fully Coherent Methods John T. Whelan School of Gravitational Waves, Warsaw, 3 July 5 Contents Signal Model. GWs from rotating neutron star............ Exercise: JKS decomposition............
More informationEE16B - Spring 17 - Lecture 11B Notes 1
EE6B - Spring 7 - Lecture B Notes Murat Arcak 6 April 207 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Interpolation with Basis Functions Recall that
More informationKrylov Subspace Spectral Methods for the Time-Dependent Schrödinger Equation with Non-Smooth Potentials
Krylov Subspace Spectral Methods for the Time-Dependent Schrödinger Equation with Non-Smooth Potentials James V. Lambers Abstract This paper presents modifications of Krylov Subspace Spectral (KSS) Methods,
More information= k, (2) p = h λ. x o = f1/2 o a. +vt (4)
Traveling Functions, Traveling Waves, and the Uncertainty Principle R.M. Suter Department of Physics, Carnegie Mellon University Experimental observations have indicated that all quanta have a wave-like
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science
Time: March 10, 006, -3:30pm MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.097 (UG) Fundamentals of Photonics 6.974 (G) Quantum Electronics Spring 006
More informationAn inverse scattering problem in random media
An inverse scattering problem in random media Pedro Caro Joint work with: Tapio Helin & Matti Lassas Computational and Analytic Problems in Spectral Theory June 8, 2016 Outline Introduction and motivation
More informationMULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER
More informationLECTURE 3: Quantization and QFT
LECTURE 3: Quantization and QFT Robert Oeckl IQG-FAU & CCM-UNAM IQG FAU Erlangen-Nürnberg 14 November 2013 Outline 1 Classical field theory 2 Schrödinger-Feynman quantization 3 Klein-Gordon Theory Classical
More informationLecture 16 February 25, 2016
MTH 262/CME 372: pplied Fourier nalysis and Winter 2016 Elements of Modern Signal Processing Lecture 16 February 25, 2016 Prof. Emmanuel Candes Scribe: Carlos. Sing-Long, Edited by E. Bates 1 Outline genda:
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 informationPart 2 Introduction to Microlocal Analysis
Part 2 Introduction to Microlocal Analysis Birsen Yazıcı& Venky Krishnan Rensselaer Polytechnic Institute Electrical, Computer and Systems Engineering March 15 th, 2010 Outline PART II Pseudodifferential(ψDOs)
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 informationTHE PARAXIAL WAVE EQUATION GAUSSIAN BEAMS IN UNIFORM MEDIA:
THE PARAXIAL WAVE EQUATION GAUSSIAN BEAMS IN UNIFORM MEDIA: In point-to-point communication, we may think of the electromagnetic field as propagating in a kind of "searchlight" mode -- i.e. a beam of finite
More informationMultiresolution Analysis
Multiresolution Analysis DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Frames Short-time Fourier transform
More informationWAVELETS, SHEARLETS AND GEOMETRIC FRAMES: PART II
WAVELETS, SHEARLETS AND GEOMETRIC FRAMES: PART II Philipp Grohs 1 and Axel Obermeier 2 October 22, 2014 1 ETH Zürich 2 ETH Zürich, supported by SNF grant 146356 OUTLINE 3. Curvelets, shearlets and parabolic
More informationLecture 5 (Sep. 20, 2017)
Lecture 5 8.321 Quantum Theory I, Fall 2017 22 Lecture 5 (Sep. 20, 2017) 5.1 The Position Operator In the last class, we talked about operators with a continuous spectrum. A prime eample is the position
More informationFourier Analysis Fourier Series C H A P T E R 1 1
C H A P T E R Fourier Analysis 474 This chapter on Fourier analysis covers three broad areas: Fourier series in Secs...4, more general orthonormal series called Sturm iouville epansions in Secs..5 and.6
More informationOn the Structure of Anisotropic Frames
On the Structure of Anisotropic Frames P. Grohs ETH Zurich, Seminar for Applied Mathematics ESI Modern Methods of Time-Frequency Analysis Motivation P. Grohs ESI Modern Methods of Time-Frequency Analysis
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 informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
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 informationQUANTUM MECHANICS I PHYS 516. Solutions to Problem Set # 5
QUANTUM MECHANICS I PHYS 56 Solutions to Problem Set # 5. Crossed E and B fields: A hydrogen atom in the N 2 level is subject to crossed electric and magnetic fields. Choose your coordinate axes to make
More informationEM waves: energy, resonators. Scalar wave equation Maxwell equations to the EM wave equation A simple linear resonator Energy in EM waves 3D waves
EM waves: energy, resonators Scalar wave equation Maxwell equations to the EM wave equation A simple linear resonator Energy in EM waves 3D waves Simple scalar wave equation 2 nd order PDE 2 z 2 ψ (z,t)
More informationChapter 17. Finite Volume Method The partial differential equation
Chapter 7 Finite Volume Method. This chapter focusses on introducing finite volume method for the solution of partial differential equations. These methods have gained wide-spread acceptance in recent
More informationCS711008Z Algorithm Design and Analysis
CS711008Z Algorithm Design and Analysis Lecture 5 FFT and Divide and Conquer Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 56 Outline DFT: evaluate a polynomial
More information13.1 Ion Acoustic Soliton and Shock Wave
13 Nonlinear Waves In linear theory, the wave amplitude is assumed to be sufficiently small to ignore contributions of terms of second order and higher (ie, nonlinear terms) in wave amplitude In such a
More information2. Higher-order Linear ODE s
2. Higher-order Linear ODE s 2A. Second-order Linear ODE s: General Properties 2A-1. On the right below is an abbreviated form of the ODE on the left: (*) y + p()y + q()y = r() Ly = r() ; where L is the
More informationThe method of stationary phase
Chapter The method of stationary phase In this chapter we introduce a very useful analytical tool which will allow us to find asymptotic expansions for integrals that cannot, in many cases, be calculated
More informationModal Decomposition and the Time-Domain Response of Linear Systems 1
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.151 Advanced System Dynamics and Control Modal Decomposition and the Time-Domain Response of Linear Systems 1 In a previous handout
More informationThe Hilbert Transform
The Hilbert Transform David Hilbert 1 ABSTRACT: In this presentation, the basic theoretical background of the Hilbert Transform is introduced. Using this transform, normal real-valued time domain functions
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationPendulum with a vibrating base
Pendulum with a vibrating base Chennai, March 14, 2006 1 1 Fast perturbations- rapidly oscillating perturbations We consider perturbations of a Hamiltonian perturbed by rapid oscillations. Later we apply
More informationAmmonia molecule, from Chapter 9 of the Feynman Lectures, Vol 3. Example of a 2-state system, with a small energy difference between the symmetric
Ammonia molecule, from Chapter 9 of the Feynman Lectures, Vol 3. Eample of a -state system, with a small energy difference between the symmetric and antisymmetric combinations of states and. This energy
More informationMultiscale Frame-based Kernels for Image Registration
Multiscale Frame-based Kernels for Image Registration Ming Zhen, Tan National University of Singapore 22 July, 16 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image
More informationReverberation Impulse Response Analysis
MUS424: Signal Processing Techniques for Digital Audio Effects Handout #18 Jonathan Abel, David Berners April 27, 24 Lecture #9: April 27, 24 Lecture Notes 7 Reverberation Impulse Response Analysis 1 CCRMA
More informationThe Evolution of Large-Amplitude Internal Gravity Wavepackets
The Evolution of Large-Amplitude Internal Gravity Wavepackets Sutherland, Bruce R. and Brown, Geoffrey L. University of Alberta Environmental and Industrial Fluid Dynamics Laboratory Edmonton, Alberta,
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationConstruction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling
Construction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling Rujie Yin Department of Mathematics Duke University USA Email: rujie.yin@duke.edu arxiv:1602.04882v1 [math.fa] 16 Feb
More informationTheoretical Photochemistry WiSe 2016/17
Theoretical Photochemistry WiSe 2016/17 Lecture 8 Irene Burghardt burghardt@chemie.uni-frankfurt.de) http://www.theochem.uni-frankfurt.de/teaching/ Theoretical Photochemistry 1 Topics 1. Photophysical
More informationMultiple scale methods
Multiple scale methods G. Pedersen MEK 3100/4100, Spring 2006 March 13, 2006 1 Background Many physical problems involve more than one temporal or spatial scale. One important example is the boundary layer
More information