Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Computational high frequency waves in heterogeneous media Shi Jin University of Wisconsin-Madison High frequency waves • • Fig. 1. The electromagnetic spectrum, which encompasses the visible region of light, extends from gamma rays with wave lengths of one hundredth of a nanometer to radio waves with wave lengths of one meter or greater. • High frequency waves: wave length/domain of computation <<1 Difficulty of high frequecy wave computation • Consider the example of visible lights in this lecture room: wave length: » 10-6 m computation domain » m 1d computation: 106 » 107 2d computation: 1012 » 1014 3d computation: 1018 » 1021 do not forget time! Time steps: 106 » 107 An Example: Linear Schrodinger Equation Semiclassical limit of the linear schrodinger equation If one can find the asymptotic (semiclassical) limit as -> 0 then one can just solve the limiting equation numerically. This allows coarse (underresolved) computation Semiclassical limit in the phase space Wigner Transform A convenient tool to study the semiclassical limit (Lions-Paul; Gerard, Markowich, Mauser, Poupaud) Moments of the Wigner function The connection between W and is established through the moments The semiclassical limit Other high frequency waves • The Wigner method is generic: linear wave—geometrical optics elastic waves electromagnetic waves, etc • We always end up at Liouville equation or coupled Liouville systems Ryzhik-Papanicoulou-Keller • Most recent numerical methods for high frequency waves are based on Liouville equations Liouville equation ft + r H¢ rx f - rx H ¢ r f = 0 f(t, x, ) is the density distribution of a classical particle at position x, time t, with momentum H=H(x, ) is the Hamiltonian The bicharacterisitcs of this equation is a Hamiltonian system: dx/dt = r H d/dt = -rx H Classical mechanics: H=1/2 ||2+V(x) (=> Newton’s second law) Geometrical optics: H = c(x) || The Liouville equation can be solved by method of characteristics if H is smooth Discontinuous Hamiltonians • H=1/2||2+V(x): V(x) is discontinuous- potential barrier, • H=c(x)||: c(x) is discontinuous-different index of refraction • quantum tunneling effect, semiconductor devise modeling, plasmas, geometric optics, interfaces between different materials, etc. • The semiclassical limit breaks down at barrier/interface Analytic issues ft + r H¢ rx f - rx H ¢ r f = 0 • The PDE does not make sense for discontinuous H. What is a weak solution? dx/dt = r H d/dt = -rx H • How to define a solution of systems of ODEs when the RHS is discontinuous or/and measure-valued? Numerical issues • for H=1/2||2+V(x) • since V’(x)= 1 at a discontinuity of V, one can smooth out V then Dv_i=O(1/x), thus t=O( x ) poor resoultion (for complete transmission) wrong solution (for partial transmission) Mathematical and Numerical Approaches Q: what happens before we take the high frequency limit? Snell-Decartes Law of refraction • When a plane wave hits the interface, the angles of incident and transmitted waves satisfy (n=c0/c) Our Idea: An interface condition • We introduce an interface condition for f that connects (the good) Liouville equations on both sides of the interface. - f(x+, +)=Tf(x-, )+R f(x+, -+) for +>0 H(x+, +)=H(x-,-) R: reflection rate T: transmission rate R+T=1 • • • • T, R defined from the original “microscopic” problems This gives a mathematically well-posed problem that is physically relavant We can show the interface condition is equivalent to Snell’s law in geometrical optics A new method of characteristics (bifurcate at interfaces) Solution to Hamiltonian System with discontinuous Hamiltonians T • • R Particles cross over or be reflected by the corresponding transmission or reflection coefficients (probability) Based on this definition we have also developed particle methods (both deterministic and Monte Carlo) methods Key idea in Hamiltonian-preserving schemes • consider a standard finite difference approximation V: piecewise linear approximation—allow good CFL fI,j+1/2, f-i+1/2,j ---- upwind discretization f+i+1/2, j ---- incorporating the interface condition (motivated by Berthame & Simeoni: kinetic scheme for shallow water equations with bottom topography) Scheme I (finite difference formulation) • If at xi+1/2 V is continuous, then f+i+1/2,j= f-i+1/2,j; • Otherwise, For j>0, f+i+1/2,j = f(x+i+1/2, +) = T f-(x-i+1/2, -) +R f(x+i+1/2, -+) The transmitted and reflected fluxes 1) if the particle is transmitted 2) If the particle is reflected: The New CFL condition • Note the discrete derivative of V is defined only on continuous points of V, thus t=O( x, ) Positivity, stability, l1-convergence • for first order scheme (forward Euler in time + upwind in space), under the “good” CFL condition if fn >0, then fn+1 > 0; k fn+1kl1 (x, ) · k fnkl1 (x, ) n k f k1 · C k f0k1 (except for measure-valued initial data) l1-convergence Curved interface Geometrical optics • The same idea has also been extended to geometric optics H = c(x) || with partial transmission and reflection We build in Snell’s Law into the flux References: J-Wen, semiclassical limit of Schrodiger, Comm Math Sci ’05 J-Wen, geometrical optics, JCP 06, SINUM Quantum barrier A semiclassical approach for thin barriers (with Kyle Novak, SIAM MMS, JCP) • Barrier width in the order of De Broglie length, separated by order one distance • Solve a time-independent Schrodinger equation for the local barrier/well to determine the scattering data • Solve the classical liouville equation elsewhere, using the scattering data at the interface • Ben Abdallah-Degond-Gamba (1d quantum-classical coupling) • Ben Abdallah-Tang (1d, stationary computation) A step potential ( V(x)=1/2 H(x) ) Resonant tunnelling Circular barrier (Schrodinger with =1/400) Circular barrier (semiclassical model) Circular barrier (classical model) Diffraction grating: Semiclassical Semiclasical vs Schrodinger (=1/800) Computation of diffraction (with Dongsheng Yin) Transmissions, reflections and diffractions (Type A interface) Type B interface Hamiltonian preserving+Geometric Theory of Diffraction • We uncorporate Keller’s GTD theory into the interface condition: A type B interface Another type B interface A type A interface Half plane Other applications/developments • Elastic waves (with X. Liao, JHDE) • High frequency waves in random media (with X. Liao and X. Yang, SISC) • Fast phase-flow particle method for the Liouville equations (with H. Wu and Z. Huang) Computational cost (=10-6) • Full simulation of original problem for x » t » O()=O(10-6) Dimension total cost 2d, O(1018) 3d O(1024) • Liouville based solver for diffraction x » t » O(1/3) = O(10-2) Dimension 2d, 3d total cost O(1010) O(1014) Can be much less with local mesh refinement Summary • Unique, physically relevant solution to linear transport equation with discontinuous and measure-valued solution • Probability solution to Hamiltonian system with discontinuous Hamiltonians • Finite-difference, finite-volume methods with interface condition built into the numerical flux • Particle (both deterministic and Monte-Carlo) methods • Able to compute (partial) transmission, reflection, and diffraction (by building geometrical theory of diffraction into the interface condition) for many high frequency waves (geometrical optics, semiclassical limit of Schrodinger, elastic wave, thin quantum barrier, etc.): Liouville equation + interface conditions