9512.net
甜梦文库
当前位置:首页 >> >>

The spectrum of the integral operator of electromagnetic scattering


The spectrum of the integral operator of electromagnetic scattering
Jussi Rahola Center for Scienti c Computing P.O. Box 405 FIN-02101 ESPOO FINLAND Jussi.Rahola@csc. The scattering of electromagnetic radiation is governed by an integral equation for the electric eld. The problem can be cast as a surface or a volume integral equation. Here we are interested in the volume integral equation because the potential applications involve scattering from inhomogeneous and anisotropic materials. We consider only the simplest discretization of the problem: the scatterer is divided into cubic cells and the electric eld is assumed constant inside each cell. Two example scatterers are given in Figures 1 and 2. The computational cells are shown as small spheres. The computational problem is the solution of a system of linear equations where the coe cient matrix

Introduction

is full and complex symmetric. We have applied several iterative methods in this setting and QMR was chosen as the iterative method for the production code. The matrix-vector product can be computed with two special methods: the fast Fourier transform and the fast multipole method. Figure 1: An example of a sphere where a number of particles have been removed. Figure 2: A cluster of particles generated by the di usion-limited aggregation (DLA) process.

The scattering integral equation for the electric eld E(r) is Z E(r) = Einc (r)+k 3 (m(r )2 ?1)G(r; r ) E(r ) d3 r ;
0 0 0 0

Scattering Integral Equation

where Einc is the incoming eld, G is the dyadic Green's function and

V

rr G(r; r ) = 1 + 2 g (jr ? r j) k
0 0

ikr e g(r) = 4 kr : Here m is the complex refractive index and k is the wave number. The dyadic rr contains the second-order partial derivatives, so it is essentially the Hessian matrix. The integral equation is a strongly singular one. We de ne the integral operator as Z (KE)(r) = k3 (m(r )2 ? 1)G(r; r ) E(r ) d3 r :
0 0 0 0

Now the integral equation can be written as (1 ? K)E = Einc

V

The spectrum of the integral operator K consists of those complex numbers z for which the operator (z ? K) does not have an inverse that is bounded and de ned everywhere. 1 K) which The operator (z ?K) is the same as z (1? z corresponds to a scattering problem with the same particle but with refractive index m de ned by
0

Spectrum of the integral operator

m 2 ? 1 = (m2 ? 1)=z:
0

The problem is now to nd those refractive indices m for which the scattering problem does not have an inverse that is bounded and de ned everywhere. Then the points of the spectrum are given by
0

z = (m2 ? 1)=(m 2 ? 1):
0

If the scatterer is a homogeneous sphere, the scattering problem has an analytical solution, the so called Lorenz-Mie solution.

When we computed the eigenvalues of the coe cient matrix arising from the simple discretization of a homogeneous sphere, we found that most of the eigenvalues line on a line. Figure 3 gives the numerically computed eigenvalues together with some analytically computed points of the spectrum. The points of the spectrum align nicely with the computed eigenvalues which lie outside the line. We were unable to indentify the line with the analytical spectrum. However, the line corresponds to a scattering problem with a purely imaginary refractive index which physically corresponds to non-oscillating decaying radiation. Figure 4 gives the numerical eigenvalues arising from the discretization of the DLA particle.

Numerical eigenvalues vs. spectrum

0.2

0.1

0

?0.1

?0.2

?0.3

?0.4

?0.5 0.5

1

1.5

2

2.5

3

Figure 3: Numerically computed eigenvalues (dots) and analytically computed points of the spectrum (rings) for the sphere

0.04 0.02 0 -0.02 -0.04 Imag -0.06 -0.08 -0.1 -0.12 -0.14 -0.16 0.7

0.8

0.9

1 Real

1.1

1.2

1.3

Figure 4: Numerically computed eigenvalues for the DLA particle.

The convergence of iterative methods is determined by the distribution of eigenvalues of the coe cient matrix. Let us consider the xed point problem

Convergence of Iterative Methods

x = Kx + b: (1) For any Krylov-subspace method, the residual rk = b ? (1 ? K)xk can decrease with the rate krk k = ( (K))k kr 0 k; (2) where the optimal reduction factor (K) is given by 1 : (3) (K) = ( ) = j (1) j Here is the conformal map that maps the unbounded component of the complement of spectrum of K, (K), to the outside of the closed unit disk (see Figure 5). Here we assume that the spectrum is simply connected. Figure 6 shows the ideal convergence rate together with the convergence of QMR and GMRES for a spherical scatterer.
K

φ

α
d1

0

d2

1

Figure 5: The conformal map maps the outside of the line to the outside of the unit circle.

10

1

10

0

10

-1

10 Residual

-2

10

-3

10

-4

10

-5

Ideal 10
-6

QMR GMRES

10

-7

2

3

4

5

6

7 Iterations

8

9

10

11

12

Figure 6: The convergence of QMR and GMRES together with the ideal convergence rate for a spherical scatterer.

References
1] K. Lumme and J. Rahola. Light scattering by porous dust particles in the discrete-dipole approximation. Astrophys. J. 425, 653{667, 1994. 2] O. Nevanlinna. Convergence of Iterations for Linear Equations. Birkhauser, Basel, 1993. 3] J. Rahola. Solution of dense systems of linear equations in electromagnetic scattering calculations. Licenciate's thesis, Helsinki University of Technology, 1994. 4] J. Rahola. Diagonal forms of the translation operators in the fast multipole algorithm for scattering problems. Helsinki University of Technology, Institute of Mathematics, Research Reports A 349, 1995. 5] J. Rahola. Solution of dense systems of linear equations in the discrete-dipole approximation. SIAM J. Sci. Comp. 17, No. 1, 1996. In press.



更多相关文章:
...the integral operator of electromagnetic scattering.pdf
The spectrum of the integral operator of electromagnetic scattering_专业资料。The scattering of electromagnetic radiation is governed by an integral equation for...
Iterative Solution of Dense Linear Systems.pdf
For the surface integral equation of electromagnetic scattering we have applied...The points in the spectrum of the integral operator corresponding to o -...
Note on the asymptotic approximation of a double integral ....pdf
approximation of a double integral with an angular spectrum represen_专业资料...These double integrals typically appear in electromagnetic scattering problems. ...
Spectral analysis of X-ray pulsars with the INTEGRAL ....pdf
Spectral analysis of X-ray pulsars with the INTEGRAL observatory_专业资料。...le, a resonance cyclotron scattering feature. 4. RESULTS As it was ...
Model Spectra of Rotation Powered Pulsars in the INTEGRAL ....pdf
Model Spectra of Rotation Powered Pulsars in the INTEGRAL Range_专业资料。...for the case when scatterings considerably decrease the energy of the ...
Absolutely continuous spectrum of perturbed Stark operators_....pdf
Moreover, if the integral (6) converges for a.e. ; the absolu 5 tely continuous part of the spectral measure of the operator Hq lls the whole ...
On the discrete spectrum of non-analytic matrix-valued ....pdf
cient conditions for the spectrum of matrixvalued Friedrichs model to be ?...(y, x) and the integral operator with the kernel K2 (x, y) in the ...
Spectrum of the Multi-Reggeon Compound States in Multi-Colour....pdf
Spectrum of the Multi-Reggeon Compound States in...The high-energy asymptotics of the scattering ...The lowest integral may be expressed in terms of...
...of V0332+53 over the 20042005 outburst with Integral.pdf
over the 20042005 outburst with Integral_专业资料... and a cyclotron resonant scattering feature (CRSF...of the spectrum over the outburst, we used the...
Spectrum of CHL Dyons (I)_图文.ppt
Spectrum of CHL Dyons (I)_经济学_高等教育_教育...Here NL is the number operator of 24 leftmoving...integral lattice in 28 dimensional Lorentzian space...
Trace estimates and invariance of the essential spectrum.pdf
Trace estimates and invariance of the essential spectrum_专业资料。We provide...p . (4) 4 hence V g(F ) has the integral kernel V (x)kg (x ?...
Excitation spectrum of d-wave Fermi surface deformation.pdf
Excitation spectrum of d-wave Fermi surface ...(l ≤ 2) neighbor hopping integral, and we de... The dFSD is generated by forward scattering ...
INTEGRAL observations of TeV plerions_图文.pdf
INTEGRAL observations of TeV plerions_专业资料。...Combining spectral data of the ASCA satellite with...scattering of accelerated electrons on the Cosmic ...
...Transient Black Hole Candidates observed in the INTEGRAL ....pdf
Monitoring of Transient Black Hole Candidates observed in the INTEGRAL survey_...The spectrum of this peak, that we show in the bottom panel ?gure, ...
...Zeros of the Associated Legendre Functions with Integral ....pdf
Associated Legendre Functions with Integral Degree_...(inverse) scattering theory of the one-dimensional.... . , ?κn represesent the spectral values ...
Effect of muon-nuclear inelastic scattering on high-energy ....pdf
scattering on nuclei, as well as the integral energy spectra of atmospheric...[4] for the electromagnetic strucγ ture function F2 is applied and the ...
On the multiwavelength spectrum of the microquasar 1E 1740.7-....pdf
On the multiwavelength spectrum of the microquasar...In the 10100 keV band, RXTE and INTEGRAL ...IC scattering and gamma-ray absorption in the ...
...of the Iterated Logarithm for Some Stochastic Integrals_....pdf
A. R Let Lt be a continuous version of the stochastic integral 0t < ... a(A) is the in mum of the spectrum of the self-adjoint operator ?...
INTEGRAL Cross-calibration Status.pdf
Below we present the comparison of INTEGRAL Crab spectra with the canonical ...The scattering of the ISGRI ?ux is suppressed for OSA 3.0 in comparison ...
Light-cone formulation and spin spectrum of non-critical ....pdf
-Schwarz sector to a tachyon free subspace of the fermion parity operator....direct integral of Hilbert spaces over the spectrum of momentum operators dp...
更多相关标签:

All rights reserved Powered by 甜梦文库 9512.net

copyright ©right 2010-2021。
甜梦文库内容来自网络,如有侵犯请联系客服。zhit325@126.com|网站地图