Chapter 5 Fast Multipole Methods
|
|
- Myron Willis
- 6 years ago
- Views:
Transcription
1 Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast computation of matrix-vector products in case of densely populated matrices of special structure. We consider the FMM in case of acoustic scattering. Solutions of the Helmholtz equation are given in terms of near-field and far-field expansions. Lemma 5.1 Expansion of an exterior solution Let z 1 lr 3 and B ρ1 (z 1 ) := {x lr 3 x z 1 < ρ 1 }. An outgoing solution u of the Helmholtz equation in lr 3 \ B ρ1 (z 1 ), given by means of the acoustic single layer potential u(x) = G(x, y) ϕ(y) dσ(y), (5.1) admits the representation u(x) = B ρ1 (z 1 ) α ν ( x z 1 x z 1 ) h(1) ν (k x z 1 ), (5.2) where α ν are homogeneous harmonic polynomials of degree ν given by α ν (ˆx) = ik +ν m= ν Yν m (ˆx) B ρ1 (z 1 ) j ν (k y z 1 ) Yν m ( y z 1 ) ϕ(y)dσ(y). (5.3) y z 1 The series converges uniformly on compact subsets of lr 3 \ B ρ1 (z 1 ). Proof: we obtain Using the representation G(x, y) = ik u(x) = = ik B ρ1 (z 1 ) +ν = ik m= ν B ρ1 (z 1 ) 2ν + 1 +ν m= ν P ν (ˆx ŷ) h (1) ν (k x ) j ν (k y ) = Yν m (ˆx) Yν m (ŷ) h (1) ν (k x ) j ν (k y ), G(x, y) ϕ(y) dσ(y) = h (1) ν (k x z 1 ) Yν m ( x z 1 x z 1 ) j ν (k y z 1 ) Yν m ( y z 1 ) ϕ(y) dσ(y) = y z 1
2 2 Ronald H.W. Hoppe = α ν ( x z 1 x z 1 ) h(1) ν (k x z 1 ). Lemma 5.2 Expansion of an interior solution Let z 2 lr 3 and B ρ2 (z i ) := {x lr 3 x z 2 < ρ 2 }. A solution v of the Helmholtz equation in B ρ2 (z 2 ), given by means of the acoustic single layer potential (5.1) with B ρ1 (z 1 ) replaced by B ρ2 (z 2 ), admits the representation v(x) = β ν ( x z 1 x z 1 ) j ν(k x z 2 ), (5.4) where β ν are homogeneous harmonic polynomials of degree ν given by β ν (ˆx) = ik +ν m= ν Yν m (ˆx) B ρ2 (z 2 ) h (1) ν (k y z 1 ) Yν m ( y z 1 ) ϕ(y)dσ(y).(5.5) y z 1 The series converges uniformly in B ρ2 (z 2 ). Proof: The assertion can be proved in the same way as that of Lemma 5.1. The behavior of an outgoing solution u in the far field can be described by the far-field representative. Lemma 5.3 Far-field representative Let u be an exterior solution of the Helmholtz equation as given by (5.2). Then, there holds uniformly for ẑ, where lim u(z 1 + t ẑ) ikt exp( ikt) = Ψ(ẑ) (5.6) t Ψ(ẑ) = ( i) ν α ν (ẑ) (5.7) is called the far-field representative. Proof: We refer to Thms. 2.5,2.15 in D. Colton and R. Kress; Inverse Acoustic and Electromagnetic Scattering Theory, Springer, Berlin-Heidelberg- New York, Likewise, the behavior of an interior solution v in the near field is given by means of the near-field representative.
3 Computational Electromagnetics; Chapter 1 3 Lemma 5.4 Near-field representative Let v be an interior solution of the Helmholtz equation as given by (5.4). Then, there holds v(x) = 1 Γ(ẑ) exp(ik(x z 2 ) ẑ) dσ(ẑ), (5.8) where Γ(ẑ) = is called the near-field representative. ( i) ν β ν (ẑ) (5.9) Proof: We have Γ(ẑ)exp(ik(x z 2 ) ẑ)dσ(ẑ) = ( i) ν β ν (ẑ)exp(ik(x z 2 ) ẑ)dσ(ẑ). b b Using the Jacobi-Anger expansion of exp(ik(x z 2 ) ẑ) as given by exp(ik(x z 2 ) ẑ) = we obtain = Since l=0 i l (2l + 1)j l (k x z 2 )P l ( x z 2 x z 2 ẑ), ( i) ν β ν (ẑ) i l (2l + 1)j l (k x z 2 )P l ( x z 2 ẑ) dσ(ẑ) = l=0 x z 2 ( i) ν i l x z 2 (2l + 1)j l (k x z 2 ) β ν (ẑ)p l ( ẑ) dσ(ẑ). l=0 x z 2 x z 2 β ν (ŷ)p l ( x z 2 ẑ) dσ(ẑ) = 2ν + 1 β ν( x z 2 x z 2 )δ νl, in the previous equality we can take advantage of which gives the assertion. β ν (ẑ) P l ( x z 2 ẑ) dσ(ẑ) = v(x), x z 2
4 4 Ronald H.W. Hoppe 5.2 Far-to-far, local-to-local, and far-to-local translations A powerful tool in the calculus of near- and far-field representatives is given by the translation operator: Lemma 5.5 Translation of near- and far-field representatives (i) Let Ψ 1 and Ψ 2 be far-field representatives of u with respect to z 1 and z 2. Then there holds Ψ 2 (ẑ) = exp(ik(z 2 z 1 ) ẑ) Ψ 1 (ẑ). (5.10) (ii) Likewise, if Γ 1 and Γ 2 are near-field representatives of v with respect to z 1 and z 2, we have Γ 2 (ẑ) = exp(ik(z 2 z 1 ) ẑ) Γ 1 (ẑ). (5.11) Proof: we refer to V. Rokhlin; Diagonal forms of translation operators for the Helmholtz equation in three dimensions, Appl. Comput. Harmonic Analysis 93, 1 82, Corollary 5.6 Special far-field representatives (i) The far-field representative of a point charge is given by exp(ik x y ) ) = ikh (1) 0 (k x y x y Ψ(ẑ) = ik exp(ik(z 1 y) ẑ). (ii) Let r. Then, the far-field representative of is given by exp(ik x y ) grad y r = x y exp(ik x y ) 1 (y x) r = (ik ) x y x y y x Ψ(ẑ) = k 2 exp(ik(z 1 y) ẑ) ẑ r. (iii) Assume that an outgoing solution u τ (x) = k(x, y) ϕ(y) dσ(y) τ is generated by τ Γ, ϕ C(Γ), and k(x, y) = α G(x, y) + β G(x,y). Then, n(y) its far-field representative is given by Ψ τ (ẑ) = ˆk(ẑ, y) ϕ(y) dσ(y), τ
5 Computational Electromagnetics; Chapter 1 5 where ˆk(ẑ, y) = 1 exp(ik(z 1 y) ẑ) (αik + βk 2 ẑ n(y)). Proof: Lemma 5.3 readily gives the far-field representative of ikh (1) 0 ( x y ). Assertion (i) is then a direct consequence of the transformation rule (5.10). The proof of (ii) follows in the same way observing grad y exp(ik x y ) x y Finally, taking into account that r = ik 2 h (1) 1 (k x y ) (x y) r x y Ψ τ (ẑ) = lim (u τ (z 1 + tẑ) exp( kt) ikt) = t = ϕ(y) ( lim k(z 1 + tẑ) exp( ikt)) dσ(y), t the proof of (iii) follows from (i) and (ii). τ The previous results allow us to perform local-to-local and far-to-far translations. Another important tool is the far-to-local translation which enables us to a far-field representative into a local near-field representative. Theorem 5.7 The far-to-local translation operator Let u be an outgoing solution of the Helmholtz equation in lr 3 \ B ρ1 (z 1 ) with the associated far-field representative Ψ. Moreover, assume B ρ2 (z 2 ) lr 3 with δ := z 2 z 1 > ρ 1 + ρ 2. In B ρ2 (z 2 ) we consider the near-field representative Γ n (ẑ) := ( i m (2m + 1) h (1) m (k z 2 z 1 ) P m ( z 2 z 1 z 2 z 1 ẑ)) }{{} =: µ n (z 2 z 1,ẑ). Ψ(ẑ). (5.12) Then, the function u n defined according to u n (x) := 1 Γ n (ẑ) exp(ik(x z 2 ) ẑ) dσ(ẑ) (5.13) converges to u uniformly in B ρ2 (z 2 ) u n (x) u(x) 0 (n ). Proof: First, we consider the special case u(x) = G(x, y) and y B ρ1 (z 1 ). In view of Corollary 5.6 (i) the far-field representative of u is given by Ψ(ẑ) = i k exp(ik(z 1 y) ẑ).
6 6 Ronald H.W. Hoppe Setting d := z 2 z 1, for x B ρ2 (z 2 ) we obtain u n (x) = ik Γ n (ẑ) exp(ikẑ (x z 2 )) dσ(ẑ) = = ik i m (2m + 1)h (1) m (k d ) P m (ˆd ẑ)exp(ik(x z 2 + z 1 y) ẑ)dσ(ẑ). The Jacobi-Anger expansion of exp(ikw ẑ), w := (x z 2 ) + (z 1 y) leads to u n (x) = ik i m (2m + 1)h (1) m (k d ) i l (2l + 1) j l (k w ) P l (ŵ ẑ) P m (ˆd ẑ)dσ(ẑ). l=0 Taking advantage of the orthogonality of the Legendre polynomials we finally get P l (ŵ ẑ) P m (ˆd ẑ) dσ(ẑ) = 2m + 1 δ lm P m (ŵ ˆd), u n (x) = ik ( 1) m (2m + 1) h (1) m (k d ) j m (k w ) P m (ŵ bd) ˆ = = ik (2m + 1) h (1) m (k d ) j m (k x z 1 ) P m ( ŵ ˆd) = = 4 π G n (d, w). Since d + w = x y, we know that u n (x) = G n (x, ) G(x, ) (n ) uniformly for x B ρ2 (z 2 ). The uniform convergence of grad G n to grad G follows by Theorem 2.10 in D. Colton and R. Kress; Inverse Acoustic and Electromagnetic Scattering Theory, Springer, Berlin-Heidelberg-New York, Let us finally consider u (x) = k(x, y) ϕ(y) dσ(y), where Γ, ϕ C(Γ) and the kernel k is given by k(x, y) = α G(x, y) + β G(x, y) n(y).
7 Computational Electromagnetics; Chapter 1 7 According to Corollary 5.8 (iii), we obtain the far-field representative Ψ (ẑ) = ˆk(ẑ, y) ϕ(y) dσ(y). The uniform convergence of G n to G and grad G n to grad G implies, that for n = [k(ρ 1 + ρ 2 )] + l: u,n (x) u (x) C z 1 z 2 ( ρ 1 + ρ 2 z 1 z 2 )l 1 ( α + n z 1 z 2 β ) ϕ,. Remark 5.10: We know from Chapter 3 that n has to be chosen sufficiently large, i.e., n > k(ρ 1 + ρ 2 ). In practice, the integrals over the unit sphere will be approximated by means of an appropriate quadrature formula: We denote by λ i, 0 i N, the zeroes of the (N +1)-st Legendre polynomial P N+1 and by a j := 2 (1 λ 2 j)(p N+1(λ j )) 2, 0 j N the weighting factors of the associated Gauss-Legendre quadrature formula. Setting θ j := arccos(λ j ), 0 j N, ψ l := 2πl 2N + 1, 0 l 2N, we define ˆξ ν = ˆξ (j,l) = cos(ψ l ) sin(θ j ) sin(ψ l ) cos(θ j ) cos(θ j ), w ν = w (j,l ) = 2π 2N + 1 a j. Then, the quadrature formula f(ˆx) dσ(ˆx) = w ν f(ˆξ) + R N (f) = ν Λ N ˆQ N (f) + R N (f), (5.14) where Λ N := {0, 1,..., N} {0, 1,..., 2N}, is exact for all polynomials up to the degree 2N. The following result shows that the convergence is not effected by the application of the quadrature formula ˆQ N. Theorem 5.10 Far-to-local translation by quadrature Under the assumptions of Theorem 5.8 let ũ n (x) be the approximation of u(x) = G(x, y) obtained by quadrature, i.e., ũ n (x) = ik ˆQ N (Γ n ( ) exp(ik(x z 2 ) ).
8 8 Ronald H.W. Hoppe For N = n + 1, we also have uniformly for x B ρ2 (z 2 ). Proof: and obtain: ũ n (x) u(x) (n ) As in the proof of Theorem 5.8 we set d := z 2 z 1, w : (x z 2 ) + (z 1 y) ũ n (x) u n (x) = k n i m (2m + 1) h (1) m (k d ) i l (2l + 1) j l (k w ) l=0 ( P l (ŵ ẑ)p m ( bd ˆ ẑ)dσ(ẑ) ˆQ N (P l (ŵ )P m (ˆd ))). Taking account of the facts that ˆQ N is exact for all polynomials up to degree 2N and has positive weights and that P l (t) 1, t 1, l 0, we get ũ n (x) u n (x) k (2m + 1) h (1) m (k d ) (2l + 1) j l (k w ). l=2n m In view of the asymptotics (2l + 1) j l (k w ) = we conclude that the infinite series (k w ) l (2l 1) (1 + O(1 l )) = 1 2 ( e 2 l=2n m can be estimated by the first term (2l + 1) j l (k w ) k w ) l (1 + O( 1 l l )),. (4N 2m + 1) j 2N m (k w ), provided 2N m is sufficiently large. We thus obtain ũ n (x) u n (x) C (2m + 1) h (1) m (k d ) (4N 2m + 1) j 2N m (k w ). Moreover, since h (1) m (k d ) (m ), the sum over the spherical Hankel functions can be estimated by the last term in the sum, i.e., ũ n (x) u n (x) C(2n + 1) h (1) n (k d ) (4N 2n + 1) j 2N n (k w ).
9 Computational Electromagnetics; Chapter 1 9 We are interested that the error ũ n (x) u ( x) is of the same magnitude as u n (x) u ( x). Therefore, we must have j n+1 (k w ) (4N 2n + 1) j 2N n (k w ), which can be achieved for N = n N 25 n Figure 5.1: Dependence of the relative error on n Figure 5.1 displays the dependence of the relative error ũ n (x) u(x) u(x) for n = 1, 2,..., 25, ρ 1 = ρ 2 = 1/4, z 2 z 1 = 1, and k =. We see that n > k(ρ 1 + ρ 2 ) is required and that the method breaks down for n > N.
Chapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationHomogenization of the Transmission Eigenvalue Problem for a Periodic Media
Homogenization of the Transmission Eigenvalue Problem for a Periodic Media Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
More informationFastMMLib: a generic Fast Multipole Method library. Eric DARRIGRAND. IRMAR Université de Rennes 1
FastMMLib: a generic Fast Multipole Method library Eric DARRIGRAND joint work with Yvon LAFRANCHE IRMAR Université de Rennes 1 Outline Introduction to FMM Historically SVD principle A first example of
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationSome negative results on the use of Helmholtz integral equations for rough-surface scattering
In: Mathematical Methods in Scattering Theory and Biomedical Technology (ed. G. Dassios, D. I. Fotiadis, K. Kiriaki and C. V. Massalas), Pitman Research Notes in Mathematics 390, Addison Wesley Longman,
More informationLecture 2: The Fast Multipole Method
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 2: The Fast Multipole Method Gunnar Martinsson The University of Colorado at Boulder Research support by: Recall:
More informationRepresentation of sound fields for audio recording and reproduction
Representation of sound fields for audio recording and reproduction F. M. Fazi a, M. Noisternig b and O. Warusfel b a University of Southampton, Highfield, SO171BJ Southampton, UK b Institut de Recherche
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationThe Interior Transmission Eigenvalue Problem for Maxwell s Equations
The Interior Transmission Eigenvalue Problem for Maxwell s Equations Andreas Kirsch MSRI 2010 epartment of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research
More informationLet b be the distance of closest approach between the trajectory of the center of the moving ball and the center of the stationary one.
Scattering Classical model As a model for the classical approach to collision, consider the case of a billiard ball colliding with a stationary one. The scattering direction quite clearly depends rather
More informationProblem Set #4: 4.1,4.7,4.9 (Due Monday, March 25th)
Chapter 4 Multipoles, Dielectrics Problem Set #4: 4.,4.7,4.9 (Due Monday, March 5th 4. Multipole expansion Consider a localized distribution of charges described by ρ(x contained entirely in a sphere of
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationHarmonic Functions and Brownian motion
Harmonic Functions and Brownian motion Steven P. Lalley April 25, 211 1 Dynkin s Formula Denote by W t = (W 1 t, W 2 t,..., W d t ) a standard d dimensional Wiener process on (Ω, F, P ), and let F = (F
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationPHYS 502 Lecture 8: Legendre Functions. Dr. Vasileios Lempesis
PHYS 502 Lecture 8: Legendre Functions Dr. Vasileios Lempesis Introduction Legendre functions or Legendre polynomials are the solutions of Legendre s differential equation that appear when we separate
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationFRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES
Progress In Electromagnetics Research, PIER 5, 3 38, 000 FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES Q. A. Naqvi and A. A. Rizvi Communications Lab. Department of Electronics Quaidi-i-Azam University
More informationIntroduction to the Boundary Element Method
Introduction to the Boundary Element Method Salim Meddahi University of Oviedo, Spain University of Trento, Trento April 27 - May 15, 2015 1 Syllabus The Laplace problem Potential theory: the classical
More informationOrthogonality Sampling for Object Visualization
Orthogonality ampling for Object Visualization Roland Potthast October 31, 2007 Abstract The goal of this paper is to propose a new sampling algorithm denoted as orthogonality sampling for the detection
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationA coupled BEM and FEM for the interior transmission problem
A coupled BEM and FEM for the interior transmission problem George C. Hsiao, Liwei Xu, Fengshan Liu, Jiguang Sun Abstract The interior transmission problem (ITP) is a boundary value problem arising in
More informationFast algorithms for integral formulations of steady-state radiative transfer equation
Fast algorithms for integral formulations of steady-state radiative transfer equation Yuwei Fan, Jing An, Lexing Ying February 11, 218 Abstract We investigate integral formulations and fast algorithms
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationAppendix A Vector Analysis
Appendix A Vector Analysis A.1 Orthogonal Coordinate Systems A.1.1 Cartesian (Rectangular Coordinate System The unit vectors are denoted by x, ŷ, ẑ in the Cartesian system. By convention, ( x, ŷ, ẑ triplet
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationFactorization method in inverse
Title: Name: Affil./Addr.: Factorization method in inverse scattering Armin Lechleiter University of Bremen Zentrum für Technomathematik Bibliothekstr. 1 28359 Bremen Germany Phone: +49 (421) 218-63891
More informationHigh-Order Corrected Trapezoidal Quadrature Rules for Functions with a Logarithmic Singularity in 2-D
High-Order Corrected Trapezoidal Quadrature Rules for Functions with a Logarithmic Singularity in 2-D Juan C. Aguilar 1 Yu Chen 2 April 24, 2002 Abstract In this report we construct correction coefficients
More informationUniqueness in determining refractive indices by formally determined far-field data
Applicable Analysis, 2015 Vol. 94, No. 6, 1259 1269, http://dx.doi.org/10.1080/00036811.2014.924215 Uniqueness in determining refractive indices by formally determined far-field data Guanghui Hu a, Jingzhi
More informationInverse Obstacle Scattering
, Göttingen AIP 2011, Pre-Conference Workshop Texas A&M University, May 2011 Scattering theory Scattering theory is concerned with the effects that obstacles and inhomogenities have on the propagation
More informationSOLUTION OF POISSON S EQUATION. Contents
SOLUTION OF POISSON S EQUATION CRISTIAN E. GUTIÉRREZ OCTOBER 5, 2013 Contents 1. Differentiation under the integral sign 1 2. The Newtonian potential is C 1 2 3. The Newtonian potential from the 3rd Green
More informationPerron method for the Dirichlet problem.
Introduzione alle equazioni alle derivate parziali, Laurea Magistrale in Matematica Perron method for the Dirichlet problem. We approach the question of existence of solution u C (Ω) C(Ω) of the Dirichlet
More informationarxiv: v3 [math.ap] 4 Jan 2017
Recovery of an embedded obstacle and its surrounding medium by formally-determined scattering data Hongyu Liu 1 and Xiaodong Liu arxiv:1610.05836v3 [math.ap] 4 Jan 017 1 Department of Mathematics, Hong
More informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationDiscreteness of Transmission Eigenvalues via Upper Triangular Compact Operators
Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research
More informationEXISTENCE OF GUIDED MODES ON PERIODIC SLABS
SUBMITTED FOR: PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 8 EXISTENCE OF GUIDED MODES ON PERIODIC SLABS Stephen
More informationFast Multipole Methods: Fundamentals & Applications. Ramani Duraiswami Nail A. Gumerov
Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov Week 1. Introduction. What are multipole methods and what is this course about. Problems from physics, mathematics,
More informationINFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station,
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationGreen s Functions, Boundary Integral Equations and Rotational Symmetry
Green s Functions, Boundary Integral Equations and Rotational Symmetry...or, How to Construct a Fast Solver for Stokes Equation Saibal De Advisor: Shravan Veerapaneni University of Michigan, Ann Arbor
More informationEfficient Evaluation of the Rokhlin Translator in Multilevel Fast Multipole Algorithm
Ilari Hänninen and Jukka Sarvas, Efficient evaluation of the Rokhlin translator in multilevel fast multipole algorithm, Helsinki University of Technology, Electromagnetics Laboratory Report Series, Report
More informationarxiv: v1 [math.ap] 21 Dec 2018
Uniqueness to Inverse Acoustic and Electromagnetic Scattering From Locally Perturbed Rough Surfaces Yu Zhao, Guanghui Hu, Baoqiang Yan arxiv:1812.09009v1 [math.ap] 21 Dec 2018 Abstract In this paper, we
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationModified Bessel functions : Iα, Kα
Modified Bessel functions : Iα, Kα The Bessel functions are valid even for complex arguments x, and an important special case is that of a purely imaginary argument. In this case, the solutions to the
More informationInverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing
Inverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing Isaac Harris Texas A & M University College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: F. Cakoni, H.
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationAngular Momentum - set 1
Angular Momentum - set PH0 - QM II August 6, 07 First of all, let us practise evaluating commutators. Consider these as warm up problems. Problem : Show the following commutation relations ˆx, ˆL x = 0,
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationTopological properties of Z p and Q p and Euclidean models
Topological properties of Z p and Q p and Euclidean models Samuel Trautwein, Esther Röder, Giorgio Barozzi November 3, 20 Topology of Q p vs Topology of R Both R and Q p are normed fields and complete
More informationThe Factorization Method for Maxwell s Equations
The Factorization Method for Maxwell s Equations Andreas Kirsch University of Karlsruhe Department of Mathematics 76128 Karlsruhe Germany December 15, 2004 Abstract: The factorization method can be applied
More informationNotes on Transmission Eigenvalues
Notes on Transmission Eigenvalues Cédric Bellis December 28, 2011 Contents 1 Scattering by inhomogeneous medium 1 2 Inverse scattering via the linear sampling method 2 2.1 Relationship with the solution
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationLaplace equation. In this chapter we consider Laplace equation in d-dimensions given by. + u x2 x u xd x d. u x1 x 1
Chapter 6 Laplace equation In this chapter we consider Laplace equation in d-dimensions given by u x1 x 1 + u x2 x 2 + + u xd x d =. (6.1) We study Laplace equation in d = 2 throughout this chapter (excepting
More informationELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES. Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia
ELECTROMAGNETIC SCATTERING FROM PERTURBED SURFACES Katharine Ott Advisor: Irina Mitrea Department of Mathematics University of Virginia Abstract This paper is concerned with the study of scattering of
More informationElectrostatic Imaging via Conformal Mapping. R. Kress. joint work with I. Akduman, Istanbul and H. Haddar, Paris
Electrostatic Imaging via Conformal Mapping R. Kress Göttingen joint work with I. Akduman, Istanbul and H. Haddar, Paris Or: A new solution method for inverse boundary value problems for the Laplace equation
More informationWavelet Frames on the Sphere for Sparse Representations in High Angular Resolution Diusion Imaging. Chen Weiqiang
Wavelet Frames on the Sphere for Sparse Representations in High Angular Resolution Diusion Imaging Chen Weiqiang Overview 1. Introduction to High Angular Resolution Diusion Imaging (HARDI). 2. Wavelets
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationAngular Momentum. Classically the orbital angular momentum with respect to a fixed origin is. L = r p. = yp z. L x. zp y L y. = zp x. xpz L z.
Angular momentum is an important concept in quantum theory, necessary for analyzing motion in 3D as well as intrinsic properties such as spin Classically the orbital angular momentum with respect to a
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationMechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations
Appl. Math. Mech. -Engl. Ed., 32(12), 1505 1514 (2011) DOI 10.1007/s10483-011-1519-7 c Shanghai University and Springer-Verlag Berlin Heidelberg 2011 Applied Mathematics and Mechanics (English Edition)
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More information44 CHAPTER 3. CAVITY SCATTERING
44 CHAPTER 3. CAVITY SCATTERING For the TE polarization case, the incident wave has the electric field parallel to the x 3 -axis, which is also the infinite axis of the aperture. Since both the incident
More informationMax-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig Hierarchical matrices and the High-Frequency Fast Multipole Method for the Helmholtz Equation with Decay revised version: March 014)
More informationScattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
More informationFinite and Boundary Element Methods in Acoustics
Finite and Boundary Element Methods in Acoustics W. Kreuzer, Z. Chen, H. Waubke Austrian Academy of Sciences, Acoustics Research Institute ARI meets NuHAG Kreuzer, Chen, Waubke (ARI) FEM-BEM-FMM ARI meets
More information(1) u i = g(x i, x j )q j, i = 1, 2,..., N, where g(x, y) is the interaction potential of electrostatics in the plane. 0 x = y.
Encyclopedia entry on Fast Multipole Methods. Per-Gunnar Martinsson, University of Colorado at Boulder, August 2012 Short definition. The Fast Multipole Method (FMM) is an algorithm for rapidly evaluating
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More informationAn adaptive fast multipole boundary element method for the Helmholtz equation
An adaptive fast multipole boundary element method for the Helmholtz equation Vincenzo Mallardo 1, Claudio Alessandri 1, Ferri M.H. Aliabadi 2 1 Department of Architecture, University of Ferrara, Italy
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 3, pp. 1002 1026 c 2006 Society for Industrial and Applied Mathematics ACOUSTIC SCATTERING BY MILDLY ROUGH UNBOUNDED SURFACES IN THREE DIMENSIONS SIMON N. CHANDLER-WILDE,
More informationThe Factorization Method for a Class of Inverse Elliptic Problems
1 The Factorization Method for a Class of Inverse Elliptic Problems Andreas Kirsch Mathematisches Institut II Universität Karlsruhe (TH), Germany email: kirsch@math.uni-karlsruhe.de Version of June 20,
More informationMatrix construction: Singular integral contributions
Matrix construction: Singular integral contributions Seminar Boundary Element Methods for Wave Scattering Sophie Haug ETH Zurich November 2010 Outline 1 General concepts in singular integral computation
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationAkira Ishimaru, Sermsak Jaruwatanadilok and Yasuo Kuga
INSTITUTE OF PHYSICS PUBLISHING Waves Random Media 4 (4) 499 5 WAVES IN RANDOMMEDIA PII: S959-774(4)789- Multiple scattering effects on the radar cross section (RCS) of objects in a random medium including
More informationElectromagnetism Answers to Problem Set 9 Spring 2006
Electromagnetism 70006 Answers to Problem et 9 pring 006. Jackson Prob. 5.: Reformulate the Biot-avart law in terms of the solid angle subtended at the point of observation by the current-carrying circuit.
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationThe Mollifier Theorem
The Mollifier Theorem Definition of the Mollifier The function Tx Kexp 1 if x 1 2 1 x 0 if x 1, x R n where the constant K is chosen such that R n Txdx 1, is a test function on Rn. Note that Tx vanishes,
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationON THE MATHEMATICAL BASIS OF THE LINEAR SAMPLING METHOD
Georgian Mathematical Journal Volume 10 (2003), Number 3, 411 425 ON THE MATHEMATICAL BASIS OF THE LINEAR SAMPLING METHOD FIORALBA CAKONI AND DAVID COLTON Dedicated to the memory of Professor Victor Kupradze
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationA FAST DIRECT SOLVER FOR THE INTEGRAL EQUATIONS OF SCATTERING THEORY ON PLANAR CURVES WITH CORNERS
A FAST DIRECT SOLVER FOR THE INTEGRAL EQUATIONS OF SCATTERING THEORY ON PLANAR CURVES WITH CORNERS JAMES BREMER Abstract. We describe an approach to the numerical solution of the integral equations of
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationUniqueness in an inverse acoustic obstacle scattering problem for both sound-hard and sound-soft polyhedral scatterers
INSTITUTE OF PHYSICS PUBLISHING Inverse Problems 22 (2006) 515 524 INVERSE PROBLEMS doi:10.1088/0266-5611/22/2/008 Uniqueness in an inverse acoustic obstacle scattering problem for both sound-hard and
More informationExplicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems
z Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems Johan Helsing Lund University Talk at Integral equation methods: fast algorithms and applications, Banff,
More information