# Analytic Calculation of Transmission Field in Homogeneously Layered Mediums Excited by EMP

Authors: Dong-yang Sun; Wei-min Bao; Xiao-ping Li
Journal: Advances in Mathematical Physics
 (2017)

---

## Abstract

This paper presents an analytic derivation for the time-domain transmission across layered mediums. The transmission coefficient and attenuation coefficient are obtained in the time-domain from general electromagnetic theory. The transmission electric field can be obtained within a few seconds by convolving the coefficients with incident EMP. The results are accordant with the FDTD method, and this approach can deal with the multilayer mediums problem. The limitations of this approach are discussed in this paper.


---

## Body

## 1. Introduction

Numerical simulation has become an indispensable tool for EMC studies, and these numerical methods can potentially simulate many physical problems [1]. Meanwhile, the analytic method has not been used widely, which is restricted by complex structure in practical problem. But analytic expressions for some ideal models can reveal the physical essential of these problems. And they are the origin of some numerical methods. It is the contention of this paper that the analytic derivation can offer an efficient method for some ideal physical problems in the time-domain.Some installations such as subway or shelter are built under the ground. The depth increases the computational costs of numerical simulation about these installations excited by EMP (electromagnetic pulse) or other electromagnetic excitations over the ground. There are two steps for solving this kind of problem: first, establishing the time-domain transmission field across the ground; second, calculating the response of the installation excited by the transmission field using numerical simulation method. This paper addresses the first step.The layered medium is a canonical structure in the electromagnetic field theory. The transmission in the frequency domain is determined by Snell’s law of refraction. Several authors obtain the time-domain transmission field by performing a numerical FFT [2, 3]. Some improved time-domain methods are presented for reflected field, which do not refer to the transmission [4–7].In this paper, the analytic derivation for the transmission coefficient and attenuation coefficient are presented, and the transmission electric field across layered medium is obtained by convolving the coefficients with incident EMP. The transmission electric field in monolayer medium is accordant with the FDTD method. The transmission electric field in dry soil and freshwater are compared, whose peculiarity is conforming to physical laws. This approach is also suitable for multilayer problem, and the limitations of this approach are discussed.


## 2. The Transmission Coefficient and Attenuation Coefficient in the Time-Domain

### 2.1. Transmission Coefficient in the Time-Domain

A plane wave obliquely striking a half-space interface from the optically thinner medium to optically denser medium in the frequency domain is considered. Figure1 illustrates the geometry of this canonical problem, where the electric field has the horizontal polarization with the strict Snell transmission coefficient TE(ω) as shown in the following [8, 9]:(1)TEω=2cos⁡θicos⁡θi+ε2μ1/ε1μ21-ε1μ1/ε2μ2sin2⁡θi.Figure 1
Plane wave obliquely strikes a half-space interface.The material permittivity isε1=ε0εr1-jσ1/ω and ε2=ε0εr2-jσ2/ω. The material permeability is μ1=μ0 and μ2=μ0. Here ε0 is the permittivity of vacuum space, μ0 is the permeability of vacuum space, ω is angle frequency, εr1 and εr2 are relative permittivity, σ1 and σ2 are electrical conductivity, θi is incident angle, and θt is refraction angle. Using the relation s↔jω [10], the transmission coefficient is represented as follows:(2)TEs=2cos⁡θicos⁡θi+ε0εr2+σ2/s/ε0εr1+σ1/s-sin2⁡θi,where s is the Laplace Transform variable.Equation (2) can be transformed to (3) through some cumbersome transformations.(3)TEs-γ=a+bcd∑n=1∞Sn,where γ=γ1+γ2/2, γ1=σ1/ε0εr1, γ2=(σ2-σ1sin2⁡θi)/(ε0εr2-ε0εr1sin2⁡θi), a=2χ1/(1+χ1), b=-4χ1/1+χ12, c=(γ2-γ1)/2, d=c((χ1-1)/(χ1+1)), χ1=cos⁡θiεr1/(εr2-εr1sin2⁡θi), and S=d/(s+s2-c2).The inverse Laplace Transform of (2) is(4)TEt=aδte-γt+bcde-γtt∑n=1∞ndcnInct,where Inct=∑k=0∞(1/k!Γk+n+1)ct/2n+2k is deformation Bessel function and δ(t) is Dirac function.The derivation for the transmission coefficient is cumbersome, whose details are recorded in AppendixA at the end of the paper. The transmission coefficient of vertical polarization wave is not listed, which can be obtained using similar approach.


### 2.2. Attenuation Coefficient in the Time-Domain

The attenuation coefficient of electric wave propagation in lossy medium is calculated by an exponential function as shown in the following [11]:(5)KEω=exp⁡-kI·r,where kI is the attenuation wave vector, whose scalar value is determined by (6); r is propagation distance vector, whose scalar value is determined by (7):(6)kIω=ω2με21+σωε2-1,(7)r=zcos⁡θt,where z is vertical depth.In (6), μ=μ2=μ0, ε=ε2=ε0εr2-jσ2/ω, s↔jω; using the power series expansion of binomial expression (see (8)), the kI(ω) can be represented as (9) where high order terms are ignored:(8)1-a=1-12a-18a2-116a3-⋯,(9)kIs=σ22sμ0sε0εr2+σ2=σ22μ0ε0εr2+σ2/s.The magnitude of thekIs is proportional to s and s represents the frequency. This relationship means that high frequency wave attenuates more rapidly than low frequency wave in lossy medium.Similarly, (7) needs to be represented by the Laplace Transform variable s. Snell’s Law of refraction is(10)sin⁡θtsin⁡θi=μ1ε1μ2ε2=ε0εr1+σ1/sε0εr2+σ2/s.The refraction angleθt is determined by the incident angle, the medium parameters, and the transmission wave frequency s (or ω):(11)cos⁡θt=1-sin2⁡θisε0εr1+σ1sε0εr2+σ2.Combining with the results of upper derivation and using the series expansion, the attenuation coefficient is represented as follows:(12)KEs=exp⁡-xss+y=∑n=0∞-1nxnn!ss+yn/2,where x=zσ2μ0/2ε0εr2-ε0εr1sin2⁡θi and y=(σ2-sin2⁡θiσ1)/(ε0εr2-ε0εr1sin2⁡θi).Then/2 increases the difficulty of inverse Laplace Transform for (12). The typical items are(13)κ1s=ss+y,κ2s=ss+y.The inverse Laplace Transforms of the typical items listed in (13) are(14)κ1t=δt+12ye-y/2tI1y2t-I0y2t,κ2t=δt-ye-yt.The derivation details are recorded in AppendixB at the end of the paper. The inverse Laplace Transforms of other items in (12) can be obtained by several convolutions. Therefore the attenuation coefficient expression in time-domain is(15)KEt=δt-xκ1t+x22!κ2t-x33!κ1t∗κ2t+x44!κ2t∗κ2t-x55!κ1t∗κ2t∗κ2t+⋯.In order to obtain the transmission field in the time-domain, the transmission coefficient and the attenuation coefficient are represented in the time-domain. The transmission electric field across lossy medium is(16)ETt=Eit∗TEt∗KEt.The analytic expression of (16) is too cumbersome to simplify. In order to obtain the final transmission electric field ETt, a numerical convolution is chosen to manipulate (16) by MATLAB software; meanwhile the high order terms in (15) are ignored for simplification.


### 2.3. The Limitations of Coefficients

Firstly, the parametery in attenuation coefficient has a large value from its expression. The exponential function e-yt attenuates rapidly if σ2 is not too smaller than εr2 such as seawater (εr2=81,σ2=4) [11]; ye-yt attenuates to zero within one nanosecond as shown in Figure 2; therefore κ1(t) and κ2(t) attenuate to zero within one nanosecond (I1(t)<I0(t)).Figure 2
e
-
y
t about seawater and wet soil.The attenuation coefficient of seawater is ineffectual for some transient electromagnetic wave; for example, the duration of EMP is tens of nanoseconds; that is, it is ineffectual to good conductor. Otherwise the attenuation coefficient is effectual to some other usual media like wet soil(εr2=10,σ2=10-3), dry soil (εr2=4,σ2=10-5), freshwater (εr2=81,σ2=10-3), and so on. Figure 2 illustrates the duration of e-yt about seawater and wet soil for contrast.Secondly, the itemxn within the sum fromn equals zero to +∞ in (12) is unharmonious; when x is greater than one, the attenuation coefficient may be vertiginous. There is a similar restriction in (4), but the parameter d is less than the parameter c from their expressions; the transmission coefficient is effectual naturally. By the restriction, the effectual depth of attenuation coefficient is restricted by (17) and is listed in Table 1 for some usual mediums:(17)z<2ε0εr2-ε0εr1sin2⁡θiσ2μ0.Table 1
The effectual depth about attenuation coefficient.
Medium
Freshwater
Wet soil
Dry soil
Effectual depth (m) whenθi = 0°
47.8
16.8
1062
Effectual depth (m) whenθi = 90°
47.5
15.9
920


## 2.1. Transmission Coefficient in the Time-Domain

A plane wave obliquely striking a half-space interface from the optically thinner medium to optically denser medium in the frequency domain is considered. Figure1 illustrates the geometry of this canonical problem, where the electric field has the horizontal polarization with the strict Snell transmission coefficient TE(ω) as shown in the following [8, 9]:(1)TEω=2cos⁡θicos⁡θi+ε2μ1/ε1μ21-ε1μ1/ε2μ2sin2⁡θi.Figure 1
Plane wave obliquely strikes a half-space interface.The material permittivity isε1=ε0εr1-jσ1/ω and ε2=ε0εr2-jσ2/ω. The material permeability is μ1=μ0 and μ2=μ0. Here ε0 is the permittivity of vacuum space, μ0 is the permeability of vacuum space, ω is angle frequency, εr1 and εr2 are relative permittivity, σ1 and σ2 are electrical conductivity, θi is incident angle, and θt is refraction angle. Using the relation s↔jω [10], the transmission coefficient is represented as follows:(2)TEs=2cos⁡θicos⁡θi+ε0εr2+σ2/s/ε0εr1+σ1/s-sin2⁡θi,where s is the Laplace Transform variable.Equation (2) can be transformed to (3) through some cumbersome transformations.(3)TEs-γ=a+bcd∑n=1∞Sn,where γ=γ1+γ2/2, γ1=σ1/ε0εr1, γ2=(σ2-σ1sin2⁡θi)/(ε0εr2-ε0εr1sin2⁡θi), a=2χ1/(1+χ1), b=-4χ1/1+χ12, c=(γ2-γ1)/2, d=c((χ1-1)/(χ1+1)), χ1=cos⁡θiεr1/(εr2-εr1sin2⁡θi), and S=d/(s+s2-c2).The inverse Laplace Transform of (2) is(4)TEt=aδte-γt+bcde-γtt∑n=1∞ndcnInct,where Inct=∑k=0∞(1/k!Γk+n+1)ct/2n+2k is deformation Bessel function and δ(t) is Dirac function.The derivation for the transmission coefficient is cumbersome, whose details are recorded in AppendixA at the end of the paper. The transmission coefficient of vertical polarization wave is not listed, which can be obtained using similar approach.


## 2.2. Attenuation Coefficient in the Time-Domain

The attenuation coefficient of electric wave propagation in lossy medium is calculated by an exponential function as shown in the following [11]:(5)KEω=exp⁡-kI·r,where kI is the attenuation wave vector, whose scalar value is determined by (6); r is propagation distance vector, whose scalar value is determined by (7):(6)kIω=ω2με21+σωε2-1,(7)r=zcos⁡θt,where z is vertical depth.In (6), μ=μ2=μ0, ε=ε2=ε0εr2-jσ2/ω, s↔jω; using the power series expansion of binomial expression (see (8)), the kI(ω) can be represented as (9) where high order terms are ignored:(8)1-a=1-12a-18a2-116a3-⋯,(9)kIs=σ22sμ0sε0εr2+σ2=σ22μ0ε0εr2+σ2/s.The magnitude of thekIs is proportional to s and s represents the frequency. This relationship means that high frequency wave attenuates more rapidly than low frequency wave in lossy medium.Similarly, (7) needs to be represented by the Laplace Transform variable s. Snell’s Law of refraction is(10)sin⁡θtsin⁡θi=μ1ε1μ2ε2=ε0εr1+σ1/sε0εr2+σ2/s.The refraction angleθt is determined by the incident angle, the medium parameters, and the transmission wave frequency s (or ω):(11)cos⁡θt=1-sin2⁡θisε0εr1+σ1sε0εr2+σ2.Combining with the results of upper derivation and using the series expansion, the attenuation coefficient is represented as follows:(12)KEs=exp⁡-xss+y=∑n=0∞-1nxnn!ss+yn/2,where x=zσ2μ0/2ε0εr2-ε0εr1sin2⁡θi and y=(σ2-sin2⁡θiσ1)/(ε0εr2-ε0εr1sin2⁡θi).Then/2 increases the difficulty of inverse Laplace Transform for (12). The typical items are(13)κ1s=ss+y,κ2s=ss+y.The inverse Laplace Transforms of the typical items listed in (13) are(14)κ1t=δt+12ye-y/2tI1y2t-I0y2t,κ2t=δt-ye-yt.The derivation details are recorded in AppendixB at the end of the paper. The inverse Laplace Transforms of other items in (12) can be obtained by several convolutions. Therefore the attenuation coefficient expression in time-domain is(15)KEt=δt-xκ1t+x22!κ2t-x33!κ1t∗κ2t+x44!κ2t∗κ2t-x55!κ1t∗κ2t∗κ2t+⋯.In order to obtain the transmission field in the time-domain, the transmission coefficient and the attenuation coefficient are represented in the time-domain. The transmission electric field across lossy medium is(16)ETt=Eit∗TEt∗KEt.The analytic expression of (16) is too cumbersome to simplify. In order to obtain the final transmission electric field ETt, a numerical convolution is chosen to manipulate (16) by MATLAB software; meanwhile the high order terms in (15) are ignored for simplification.


## 2.3. The Limitations of Coefficients

Firstly, the parametery in attenuation coefficient has a large value from its expression. The exponential function e-yt attenuates rapidly if σ2 is not too smaller than εr2 such as seawater (εr2=81,σ2=4) [11]; ye-yt attenuates to zero within one nanosecond as shown in Figure 2; therefore κ1(t) and κ2(t) attenuate to zero within one nanosecond (I1(t)<I0(t)).Figure 2
e
-
y
t about seawater and wet soil.The attenuation coefficient of seawater is ineffectual for some transient electromagnetic wave; for example, the duration of EMP is tens of nanoseconds; that is, it is ineffectual to good conductor. Otherwise the attenuation coefficient is effectual to some other usual media like wet soil(εr2=10,σ2=10-3), dry soil (εr2=4,σ2=10-5), freshwater (εr2=81,σ2=10-3), and so on. Figure 2 illustrates the duration of e-yt about seawater and wet soil for contrast.Secondly, the itemxn within the sum fromn equals zero to +∞ in (12) is unharmonious; when x is greater than one, the attenuation coefficient may be vertiginous. There is a similar restriction in (4), but the parameter d is less than the parameter c from their expressions; the transmission coefficient is effectual naturally. By the restriction, the effectual depth of attenuation coefficient is restricted by (17) and is listed in Table 1 for some usual mediums:(17)z<2ε0εr2-ε0εr1sin2⁡θiσ2μ0.Table 1
The effectual depth about attenuation coefficient.
Medium
Freshwater
Wet soil
Dry soil
Effectual depth (m) whenθi = 0°
47.8
16.8
1062
Effectual depth (m) whenθi = 90°
47.5
15.9
920


## 3. The Transmission Electric Field Excited by EMP

### 3.1. Verification by the FDTD Method

The incident EMP from air(εr1=1,σ1=0) to layered medium can be approximated by a double exponential function in time-domain [12]:(18)Eit=k0E0e-αt-e-βt,where k0=1.3, E0 = 50 (kV/m), α = 4 × 107 (1/s), and β = 6 × 108 (1/s); Figure 3 illustrates the waveform of the EMP.Figure 3
Incident EMP (IEC61000-2-9).Convolving the coefficients with the incident EMP, the transmission electric fields are obtained whenθi = 0 in freshwater and dry soil. Use the FDTD method to calculate the transmission electric fields excited by the same EMP. Two kinds of results are compared, shown as in Figure 4.Figure 4
Analytical result is compared with the FDTD method. Black solid line refers to transmission electric field in dry soil by FDTD. Red dashed line refers to transmission electric field in dry soil by analytic. Blue solid line refers to transmission electric field in freshwater by FDTD. Pink dashed line refers to transmission electric field in freshwater by analytic.The analytic calculation result is accordant with the FDTD method. The analytic calculation is finished within a few seconds by MATLAB software performing the convolution operation.


### 3.2. The Transmission Electric Field in Lossy Medium

Considering the attenuation of the lossy medium, convolving the two coefficients with the incident EMP, the transmission electric field is obtained in different depth whenθi equal π/6. The results in freshwater and dry soil are relatively illustrated in Figures 5 and 6.Figure 5
Transmission electric fields in freshwater.Figure 6
Transmission electric fields in dry soil.Without considering the attenuation coefficient, the transmission electric field does not change aboutz. The value of transmission electric field in freshwater is less and attenuates more rapidly than the value in dry soil if considering the attenuation coefficient. Change the incident angle from 0° to 90° by 10°; the peak values of transmission electric field in dry soil and freshwater are relatively illustrated in Figure 7 when z equal 10 m.Figure 7
Peak values with different incident angles.The peak values decrease as the attenuation coefficient as long as the incident angle is increasing.


### 3.3. The Transmission Electric Field in Multilayer Mediums

The transient plane wave strikes multilayer mediums always from one optically thinner medium to optically denser medium; the transmission coefficient obtained above is effectual. Consider the geometry of the multilayer mediums as shown in Figure8 when θi equal 0°; the transmission electric fields in three mediums are illustrated in Figure 9.Figure 8
Normal incidence at multilayer mediums interfaces.Figure 9
Transmission electric fields in multiple mediums.The results do not refer to attenuation and the incident angle equals 0°. So the results in Figure9 may be the maximal transmission electric field. It is apparent that if considering the attenuation or increasing the incident angle, a lesser transmission electric field ought to be gained.There is a question here: how can these incident angles be determined potentially in different interfaces when the first incident angle is not 0°?


## 3.1. Verification by the FDTD Method

The incident EMP from air(εr1=1,σ1=0) to layered medium can be approximated by a double exponential function in time-domain [12]:(18)Eit=k0E0e-αt-e-βt,where k0=1.3, E0 = 50 (kV/m), α = 4 × 107 (1/s), and β = 6 × 108 (1/s); Figure 3 illustrates the waveform of the EMP.Figure 3
Incident EMP (IEC61000-2-9).Convolving the coefficients with the incident EMP, the transmission electric fields are obtained whenθi = 0 in freshwater and dry soil. Use the FDTD method to calculate the transmission electric fields excited by the same EMP. Two kinds of results are compared, shown as in Figure 4.Figure 4
Analytical result is compared with the FDTD method. Black solid line refers to transmission electric field in dry soil by FDTD. Red dashed line refers to transmission electric field in dry soil by analytic. Blue solid line refers to transmission electric field in freshwater by FDTD. Pink dashed line refers to transmission electric field in freshwater by analytic.The analytic calculation result is accordant with the FDTD method. The analytic calculation is finished within a few seconds by MATLAB software performing the convolution operation.


## 3.2. The Transmission Electric Field in Lossy Medium

Considering the attenuation of the lossy medium, convolving the two coefficients with the incident EMP, the transmission electric field is obtained in different depth whenθi equal π/6. The results in freshwater and dry soil are relatively illustrated in Figures 5 and 6.Figure 5
Transmission electric fields in freshwater.Figure 6
Transmission electric fields in dry soil.Without considering the attenuation coefficient, the transmission electric field does not change aboutz. The value of transmission electric field in freshwater is less and attenuates more rapidly than the value in dry soil if considering the attenuation coefficient. Change the incident angle from 0° to 90° by 10°; the peak values of transmission electric field in dry soil and freshwater are relatively illustrated in Figure 7 when z equal 10 m.Figure 7
Peak values with different incident angles.The peak values decrease as the attenuation coefficient as long as the incident angle is increasing.


## 3.3. The Transmission Electric Field in Multilayer Mediums

The transient plane wave strikes multilayer mediums always from one optically thinner medium to optically denser medium; the transmission coefficient obtained above is effectual. Consider the geometry of the multilayer mediums as shown in Figure8 when θi equal 0°; the transmission electric fields in three mediums are illustrated in Figure 9.Figure 8
Normal incidence at multilayer mediums interfaces.Figure 9
Transmission electric fields in multiple mediums.The results do not refer to attenuation and the incident angle equals 0°. So the results in Figure9 may be the maximal transmission electric field. It is apparent that if considering the attenuation or increasing the incident angle, a lesser transmission electric field ought to be gained.There is a question here: how can these incident angles be determined potentially in different interfaces when the first incident angle is not 0°?


## 4. Conclusion

For the purpose of increasing the computational efficiency in the transient electromagnetics, the transmission issue in layered mediums is addressed with analytic method. The transmission coefficient and attenuation coefficient are precise in frequency domain, performing inverse Laplace Transform, and convolving them with transient incident wave; the transmission field can be obtained. There are some complicated steps in performing inverse Laplace Transform; the attenuation coefficient in time-domain has a restricted depth and is ineffectual for good conductor such as seawater. The truncation of the infinite sum before performing convolution degrades the accuracy.Several examples of this approach are presented. The transmission electric field calculated by analytic method in monolayer medium is accordant with the FDTD method. The transmission electric field in freshwater is less than dry soil and attenuates more rapidly. The peak values of electric field and the attenuation coefficient decrease as long as the incident angle is increasing as expected. The maximal transmissions in multilayer mediums are calculated finally.



---

*Source: 3138720-2017-10-09.xml*