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.

- ON THE EIGENVALUES OF THE VOLUME INTEGRAL OPERATOR OF ELECTROMAGNETIC SCATTERING
- ON THE NORMAL SPECTRUM OF A SUBNORMAL OPERATOR
- Absolute continuity of the spectrum of a Schrodinger operator with a potential which is per
- The Spectrum of the Dirac Operator in the Linear Sigma Model with Quarks
- Spectrum of the Schrodinger operator in a perturbed periodically twisted tube
- On the Formulation of Hybrid Finite-Element and Boundary-Integral Methods for 3-D Scattering
- Statistical properties of the spectrum of the QCD Dirac operator at low energy
- Discreteness of the spectrum of the Laplace-Beltrami operator
- Microscopic universality in the spectrum of the lattice Dirac operator
- Solving the volume integral equations of electromagnetic scattering

更多相关文章：
**
...***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 -...**
Photonic Crystals - from Theory to Practice 2006.pdf
**

Ola, “Boundary*integral* equations for *the* *scattering* *of* *electromagnetic* waves...*The* existence *of* bandgaps means that *the* *spectrum* *of* *the* Maxwell *operator* ...**
研究生学位论文(含研究报告)撰写范例(英文).doc
**

for*the* excellent *spectrum* property *of* *the* N-Müller *integral* *operator*. At...significance to *the* analysis *of* *electromagnetic* *scattering* and radiation problems...**
...***of* *the* norm *of* an *integral* *operator* associated w....pdf

On*the* evaluation *of* *the* norm *of* an *integral* *operator* associated with *the* ...*The* *spectrum* *of* relativistic one-electron atoms according to Be*the* and ...**
Analysis ***of* surface *integral*.pdf

*of* surface *integral* equations in *electromagnetic* *scattering* and radiation ...rst and second kinds is *the* identity *operator* in *the* *integral* equations *of*...**
巴黎综合理工考试试题类型(物理).pdf
**

*the* wavelengths *of* *the* *electromagnetic* waves *of* *the* visible *spectrum*, *the* ...Conditions for *the* convergence *of* an *integral* in *the* cases *of* an infinite ...**
Remarks on ***the* *spectrum* *of* *the* Dirac *operator*.pdf

Remarks on*the* *spectrum* *of* *the* Dirac *operator*_专业资料。We describe a new family *of* examples *of* hypersurfaces in *the* sphere satisfying *the* limiting-case ...**
Eigenvalues ***of* *the* Dirac *Operator* on Manifolds with....pdf

1. Introduction It is well known that*the* *spectrum* *of* *the* Dirac *operator* ...*the* *integral* Schr¨dinger-Lichnerowicz formula (2.3) for compact spin o ...**
...and ***the* *spectrum* *of* *the* QCD Dirac *operator*_免费....pdf

Chiral symmetry and*the* *spectrum* *of* *the* QCD Dirac *operator* According to *the* Banks-Casher formula *the* chiral order parameter is directly related to *the* *spe*...**
***The* *spectrum* *of* critical exponents in $(vecphi^2)^2....pdf

a solution*of* *the* *spectrum* problem in terms *of* an *operator* V *of* anomalous...Each such term leads to an *integral* g0 Z dd q ( q )l1 ( (p q )...**
Splitting ***of* some nonlocalized solutions *of* *the* Kor....pdf

*scattering* transform method 10] is that for any ...*spectrum* *of* *the* *operator* L depends only on *the* ...*integral* equation is utilized, in *the* present ...**
...***of* *the* type $A$ Dunkl *operator* and *integral* kern....pdf

*The* *spectrum* *of* *the* inte... 暂无评价 11页 免费 On *the* scaling limit *of*...A $q$-analogue *of* *the* type $A$ Dunkl *operator* and *integral* kernel We int...**
Topological Charge and ***the* *Spectrum* *of* *the* Fermion ....pdf

*the* zero eigenvalues (zero modes) *of* *the* Dirac *operator* in continuous back...eld con?guration contributing to *the* lattice path *integral* in *the* continuum ...**
巴黎综合理工考试试题类型(物理)..doc
**

*the* wavelengths *of* *the* *electromagnetic* waves *of* *the* visible *spectrum*, *the* ...Conditions for *the* convergence *of* an *integral* in *the* cases *of* an infinite ...**
...spectral function ***of* one family *of* differential *operators*_....pdf

dρα (ε, Q, λ) (β = β(ε) ↓ 0), omitting*the* *integral on* *the* complement to small neighborhood *of* *the* *spectrum* *of* non-perturbed *operator*....**
...Wiener-Hopf Equations by ***Operator* Splitting.pdf

*scattering* problems 10, pp.188-189], *the* ...Circulant *integral* *operators* are *operators* *of* *the* ...*operator*, *the* *spectrum* *of* A becomes dense in ...**
***Electromagnetic* *Scattering* From Dispersive Dielectr....pdf

*integral* equations *of* *electromagnetic* *scattering* from dispersive dielectric bodies...*The* surface *of* *the* body will be denoted by . An *electromagnetic* wave with...**
...***the* *Spectrum* *of* a Transfer *Operator* associated w....pdf

Complete Determination*of* *the* *Spectrum* *of* a Transfer *Operator* associated with Intermittency Thomas Prellberg? Institut f¨r Theoretische Physik, Technische Univer...**
Composition ***of* Primary Cosmic Rays Beyond *the* Knee ....pdf

1*the* hadronic *integral* *spectrum* *of* two events ...approximation A 10] using *the* *operator* formalism ...*integral* energy *spectra* *of* *electromagnetic* showers, ... 更多相关标签：

For

Ola, “Boundary

for

On

Remarks on

1. Introduction It is well known that

Chiral symmetry and

a solution

dρα (ε, Q, λ) (β = β(ε) ↓ 0), omitting

Complete Determination

1