In this paper, we study electromagnetic (EM) wave scattering problem by many small impedance bodies. A numerical method for solving this problem is presented. The problem is solved under the physical assumptions ka << 1 , where a is the characteristic size of the bodies and k is the wave number. This problem is solved asymptotically and numerical experiments are provided to illustrate the idea of the method. Error estimate for the asymptotic solution is also discussed.
Electromagnetic scattering is the effect caused by EM waves such as light or radio waves hitting an object. The waves will then be scattered and the scattered field contains useful information about the object, see [
In [
smaller than the distance d between neighboring bodies, d = O ( a 2 − κ 3 ) where
κ ∈ [ 0,1 ) , and this distance d is much less than the wave length λ , k a ≪ 1 where k is the wave number. The distribution of these small bodies is assumed to follow this law
N ( Δ ) = 1 a 2 − κ ∫ Δ N ( x ) d x [ 1 + o ( 1 ) ] , a → 0 , (1)
in which Δ is an arbitrary open subset of the domain Ω that contains all the small bodies, N ( Δ ) is the number of the small bodies in Δ , and N is the distribution function of the bodies
N ( x ) ≥ 0, N ( x ) ∈ C ( Ω ) . (2)
The boundary impedance of the bodies is of the form ζ = h a − κ , where h is a continuous function such that Im h ≥ 0 , and κ ∈ [ 0,1 ) . The function h and constant κ can be chosen as desired.
To make the paper self-contained, the theory of EM wave scattering by one and many small impedance bodies is given in Sections 2 and 3. In Section 4, a numerical method for solving the EM scattering problem is presented. The solution of this problem is computed asymptotically and error analysis of the asymptotic solutions is also provided.
Let D be a bounded domain of one small body, a be its radius, and S be its smooth boundary, S ∈ C 1, γ , γ ∈ ( 0,1 ] . Assume that the dielectric permittivity ϵ and magnetic permittivity μ are constants. Let E and H denote the electric field and magnetic field, respectively. E 0 is the incident field and v E is the scattered field. The electromagnetic wave scattering by one small impedance body problem can be stated as follows
∇ × E = i ω μ H , in D ′ : = ℝ 3 \ D , (3)
∇ × H = − i ω ϵ E , in D ′ , (4)
[ N , [ E , N ] ] = ζ [ N , H ] , Re ζ ≥ 0, (5)
E = E 0 + v E , (6)
E 0 = E e i k α ⋅ x , E ⋅ α = 0, α ∈ S 2 , (7)
∂ v E ∂ r − i k v E = o ( 1 r ) , r : = | x | → ∞ , (8)
where ω > 0 is the frequency, k = 2 π / λ = ω ϵ μ is the wave number, k a ≪ 1 , λ is the wave length, ζ is the boundary impedance of the body, and α is a unit vector that indicates the direction of the incident wave E 0 . This incident wave satisfies the relation ∇ ⋅ E 0 = 0 . The scattered field v E satisfies the radiation condition (8). Here, N is the outward pointing unit normal to the surface S.
It is known from [
E ( x ) = E 0 ( x ) + ∇ × ∫ S g ( x , t ) J ( t ) d t , g ( x , t ) : = e i k | x − t | 4 π | x − t | , (9)
where E 0 is the incident plane wave defined in (7) and J is an unknown pseudovector. J is assumed to be tangential to S and can be found from the impedance boundary condition (5). Here E is a vector in ℝ 3 and ∇ × E is a pseudovector.
Once we have E, from (3) H can be found by the formula
H = ∇ × E i ω μ , (10)
The asymptotic formula of E when the radius a of the body D tends to zero is
E ( x ) = E 0 ( x ) + [ ∇ x g ( x , x 1 ) , Q ] , (11)
where | x − x 1 | ≫ a , and the point x 1 is an arbitrary point inside the small body D, see [
Q = ∫ S J ( t ) d t . (12)
The analytical formula for Q is derived in [
Theorem 1. One has
Q = − ζ | S | i ω μ τ ∇ × E 0 (13)
where
τ : = I 3 − b , b = ( b j m ) : = 1 | S | ∫ S N j ( s ) N m ( s ) d s , (14)
and | S | is the surface area of S.
Here, 1 ≤ j , m ≤ 3 correspond to x , y , and z coordinates in ℝ 3 , I 3 is a 3 × 3 identity matrix, and N j , 1 ≤ j ≤ 3 is the j-th component of the outer unit normal vector to the surface S.
Now, consider a domain Ω containing M small bodies D m , 1 ≤ m ≤ M , and S m are their corresponding smooth boundaries. Let D : = ∪ m = 1 M D m ⊂ Ω and D' be the complement of D in ℝ 3 . We assume that S = ∪ m = 1 M S m ∈ C 1 , γ , γ ∈ ( 0 , 1 ] . We also assume that the dielectric permittivity ϵ and magnetic permittivity μ are constants. Let E and H denote the electric field and magnetic field, respectively. E 0 is the incident field and v is the scattered field. The electromagnetic wave scattering by many small impedance bodies problem involves solving the following system
∇ × E = i ω μ H , in D ′ : = ℝ 3 \ D , D : = ∪ m = 1 M D m , (15)
∇ × H = − i ω ϵ E , in D ′ , (16)
[ N , [ E , N ] ] = ζ m [ N , H ] , on S m , 1 ≤ m ≤ M , (17)
E = E 0 + v , (18)
E 0 = E e i k α ⋅ x , E ⋅ α = 0, α ∈ S 2 . (19)
where v satisfies the radiation condition (8), ω > 0 is the frequency, k = 2 π / λ
is the wave number, k a ≪ 1 , a : = 1 2 m a x m diam D m , α is a unit vector that
indicates the direction of the incident wave E 0 , and ζ m is the boundary impedance of the body D m . These ζ m ’s are given by the following formula
ζ m = h ( x m ) a κ , 0 ≤ κ < 1 , x m ∈ D m , 1 ≤ m ≤ M , (20)
where h ( x ) is a continuous function in a bounded domain Ω ,
Re h ( x ) ≥ 0, Im ε = σ ω ≥ 0, (21)
ε = ε 0 , μ = μ 0 in Ω ′ : = ℝ 3 \ Ω . (22)
The distribution of small bodies D m , 1 ≤ m ≤ M , in Ω satisfies the following assumption
N ( Δ ) = 1 a 2 − κ ∫ Δ N ( x ) d x [ 1 + o ( 1 ) ] , a → 0, (23)
where N ( Δ ) is the number of small bodies in Δ , Δ is an arbitrary open subset of Ω ,
N ( x ) ≥ 0, N ( x ) ∈ C ( Ω ) (24)
and κ ∈ [ 0,1 ) is the parameter from (20).
From (15) and (16) we have
∇ × ∇ × E = k 2 E , k 2 = ω 2 ϵ μ , (25)
if μ = const . Once we have E, then H can be found from this relation
H = ∇ × E i ω μ . (26)
From (26) and (25), one can get (16). Thus, we need to find only E which satisfies the boundary condition (17). It was proved in [
E ( x ) = E 0 ( x ) + ∑ m = 1 M ∇ × ∫ S m g ( x , t ) J m ( t ) d t . (27)
where
Q m : = ∫ S m J m ( t ) d t . (28)
When a → 0 , the asymptotic solution for the electric field is given by
E ( x ) = E 0 ( x ) + ∑ m = 1 M [ ∇ g ( x , x m ) , Q m ] . (29)
Therefore, instead of finding J m ( t ) ,1 ≤ m ≤ M , we can just find Q m . The analytic formula for Q m is derived in [
Q m = − ζ m | S m | i ω μ τ m ∇ × E e m , 1 ≤ m ≤ M . (30)
The effective field acting on the m-th body is defined as
E e ( x m ) = E 0 ( x m ) + ∑ j ≠ m M [ ∇ g ( x , x j ) , Q j ] | x = x m , (31)
and E e m : = E e ( x m ) , x m is a point in D m . When a → 0 , the effective field E e ( x ) is asymptotically equal to the field E ( x ) in (29) as proved in [
From (20), (30), and (31), one gets
E e m = E 0 m − c a 2 − κ i ω μ ∑ j ≠ m M [ ∇ g ( x , x j ) | x = x m , τ ∇ × E e j ] h j , (32)
where c is a positive constant depending on the shape of the body S m , | S m | = c a 2 , τ m = τ : = I 3 − b , and
b = ( b j n ) : = 1 | S m | ∫ S m N j ( s ) N n ( s ) d s . (33)
Here, 1 ≤ j , n ≤ 3 correspond to x , y , and z coordinates in ℝ 3 , I 3 is a 3 × 3 identity matrix, and N j , 1 ≤ j ≤ 3 , is the j-th component of the outer unit normal vector to the surface S m .
Our goal is to find E e m in (32). Take curl of (32), set x = x j , and let A m : = ∇ × E e m , we have
A j = A 0 j − c a 2 − κ i ω μ ∑ m ≠ j M k 2 g ( x j , x m ) τ A m h m + ( τ A m ⋅ ∇ x ) ∇ g ( x , x m ) | x = x j h m , (34)
where 1 ≤ j ≤ M , see [
This linear system can be solved directly using Gaussian elimination method. Then E can be computed as follows
E e m = E 0 m − c a 2 − κ i ω μ ∑ j ≠ m M [ ∇ g ( x , x j ) | x = x m , τ A j ] ] h j . (35)
Matrix τ is of size 3x3 and can be computed as follows
τ : = I 3 − b = 2 3 I 3 , b = ( b j m ) : = 1 | S | ∫ S N j ( s ) N m ( s ) d s . (36)
Let A i : = ( X i , Y i , Z i ) then one can rewrite (34) as
F x ( i ) = X i + ∑ j ≠ i M a i j X j + ∑ j ≠ i M b i j Y j + ∑ j ≠ i M c i j Z j , (37)
F y ( i ) = Y i + ∑ j ≠ i M a ′ i j X j + ∑ j ≠ i M b ′ i j Y j + ∑ j ≠ i M c ′ i j Z j , (38)
F z ( i ) = Z i + ∑ j ≠ i M a ″ i j X j + ∑ j ≠ i M b ″ i j Y j + ∑ j ≠ i M c ″ i j Z j , (39)
where by the subscripts x , y , z the corresponding coordinates are denoted, e.g. F ( i ) = ( F x , F y , F z ) ( i ) , F ( i ) : = A 0 i , and
a i j : = [ k 2 g ( i , j ) + ∂ x ∇ g ( i , j ) x ] D j , (40)
b i j : = ∂ y ∇ g ( i , j ) x D j , (41)
c i j : = ∂ z ∇ g ( i , j ) x D j , (42)
a ′ i j : = ∂ x ∇ g ( i , j ) y D j , (43)
b ′ i j : = [ k 2 g ( i , j ) + ∂ y ∇ g ( i , j ) y ] D j , (44)
c ′ i j : = ∂ z ∇ g ( i , j ) y D j , (45)
a ″ i j : = ∂ x ∇ g ( i , j ) z D j , (46)
b ″ i j : = ∂ y ∇ g ( i , j ) z D j , (47)
c i j : = [ k 2 g ( i , j ) + ∂ z ∇ g ( i , j ) z ] D j , (48)
in which D j : = 2 3 c a 2 − κ i ω μ h j .
The error of the method presented in Section 4 can be estimated as follows. From the solution E of the electromagnetic scattering problem by many small bodies given in (27)
E ( x ) = E 0 ( x ) + ∑ m = 1 M ∇ × ∫ S m g ( x , t ) J m ( t ) d t , (49)
we can rewrite it as
E ( x ) = E 0 ( x ) + ∑ m = 1 M [ ∇ g ( x , x m ) , Q m ] + ∑ m = 1 M ∇ × ∫ S m [ g ( x , t ) − g ( x , x m ) ] J m ( t ) d t . (50)
Comparing this with the asymptotic formula for E when a → 0 given in (29)
E ( x ) = E 0 ( x ) + ∑ m = 1 M [ ∇ g ( x , x m ) , Q m ] , (51)
we have the error of this asymptotic formula is
Error = | ∑ m = 1 M ∇ × ∫ S m [ g ( x , t ) − g ( x , x m ) ] J m ( t ) d t | (52)
~ 1 4 π ( a k 2 d + a k d 2 + a d 3 ) ∑ m = 1 M | Q m | , (53)
where d = m i n m | x − x m | and
Q m : = ∫ S m J m ( t ) d t ≃ − ζ m | S m | i ω μ τ m ∇ × E e m = − ζ m | S m | i ω μ τ A m , a → 0, (54)
because
| ∇ [ g ( x , t ) − g ( x , x m ) ] | = O ( a k 2 d + a k d 2 + a d 3 ) , a = m a x m | t − x m | . (55)
Thus, to reduce the error, one needs to reduce the quantity k a .
To illustrate the idea of the method described in 4, consider a domain Ω as a unit cube placed in the first octant such that the origin is one of its vertex. This domain Ω contains M small bodies. The small bodies are particles which are distributed uniformly in the unit cube. The following physical parameters are used to solve the problem
・ Speed of wave, c = 3.0 E + 10 cm / sec .
・ Frequency, ω = 1.0 E + 14 Hz .
・ Wave number, k = 2.094395 e + 04 cm − 1 .
・ Constant κ = 0.9 .
・ Volume of the domain Ω that contains all the particles, | Ω | = 1 cm 3 .
・ Direction of plane wave, α = ( 1 , 0 , 0 ) .
・ Vector E = ( 0,1,0 ) .
・ Function N ( x ) = M a 2 − κ / | Ω | where M is the total number of particles and a is the radius of one particle.
・ Function h ( x ) = 1 .
・ Function μ ( x ) = 1 .
・ The distance between two neighboring particles, d = 1 / ( b − 1 ) cm, where b is the number of particles on a side of the cube.
・ Vector A 0 : A 0 m : = A 0 ( x m ) = ∇ × E 0 ( x ) | x = x m = ∇ × E e i k α ⋅ x | x = x m .
The radius a of the particles is chosen variously so that it satisfies and dissatisfies the assumption k a ≪ 1 . The numerical solution to EM wave scattering problem by many small impedance bodies is computed for M = 125 and 1000 particles.
M = 125 , d = 2.50 E − 01 | |||||
---|---|---|---|---|---|
a | 1.00E−04 | 5.00E−05 | 1.00E−05 | 5.00E−06 | 1.00E−06 |
Norm of E | 1.12E+01 | 1.12E+01 | 1.12E+01 | 1.12E+01 | 1.12E+01 |
Error of E | 7.07E−08 | 1.65E−08 | 5.61E−10 | 1.31E−10 | 4.46E−12 |
M = 1000 , d = 1.25 E − 01 | |||||
---|---|---|---|---|---|
a | 1.00E−04 | 5.00E−05 | 1.00E−05 | 5.00E−06 | 1.00E−06 |
Norm of E | 3.16E+01 | 3.16E+01 | 3.16E+01 | 3.16E+01 | 3.16E+01 |
Error of E | 1.13E−06 | 2.63E−07 | 8.96E−09 | 2.09E−09 | 7.11E−11 |
In this paper, we present a numerical method for solving EM wave scattering problem by many small impedance bodies. For illustration, numerical experiments are also provided. The solution to the EM wave scattering problem can be computed numerically and asymptotically using the described method, and the result is highly accurate if the assumption k a ≪ 1 is satisfied.
Tran, N.T. (2017) Numerical Method for Solving Electromagnetic Scattering Problem by Many Small Impedance Bodies. American Journal of Computational Mathematics, 7, 435-443. https://doi.org/10.4236/ajcm.2017.74031