Advances in Mathematical Physics
Volume 2017 (2017), Article ID 3138720, 8 pages
https://doi.org/10.1155/2017/3138720

Research Article

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

Dong-yang Sun,1,2 Wei-min Bao,¹ and Xiao-ping Li¹

¹School of Aerospace Science and Technology, Xidian University, Xi’an 710071, China
²Northwest Institute of Nuclear Technology, Xi’an 710024, China

Correspondence should be addressed to Dong-yang Sun; sundongyang2001@163.com

Received 18 July 2017; Accepted 29 August 2017; Published 9 October 2017

Academic Editor: Eugen Radu

Copyright © 2017 Dong-yang Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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. Figure 1 illustrates the geometry of this canonical problem, where the electric field has the horizontal polarization with the strict Snell transmission coefficient as shown in the following [8, 9]:

Figure 1: Plane wave obliquely strikes a half-space interface.

The material permittivity is and . The material permeability is and . Here is the permittivity of vacuum space, is the permeability of vacuum space, is angle frequency, and are relative permittivity, and are electrical conductivity, is incident angle, and is refraction angle. Using the relation [10], the transmission coefficient is represented as follows:

where is the Laplace Transform variable.

Equation (2) can be transformed to (3) through some cumbersome transformations.

where , , , , , , , , and .

The inverse Laplace Transform of (2) is

where is deformation Bessel function and is Dirac function.

The derivation for the transmission coefficient is cumbersome, whose details are recorded in Appendix A 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]:

where is the attenuation wave vector, whose scalar value is determined by (6); is propagation distance vector, whose scalar value is determined by (7):

where is vertical depth.

In (6), , , ; using the power series expansion of binomial expression (see (8)), the can be represented as (9) where high order terms are ignored:

The magnitude of the is proportional to and 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 . Snell’s Law of refraction is

The refraction angle is determined by the incident angle, the medium parameters, and the transmission wave frequency (or ):

Combining with the results of upper derivation and using the series expansion, the attenuation coefficient is represented as follows:

where and .

The increases the difficulty of inverse Laplace Transform for (12). The typical items are

The inverse Laplace Transforms of the typical items listed in (13) are

The derivation details are recorded in Appendix B 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

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

The analytic expression of (16) is too cumbersome to simplify. In order to obtain the final transmission electric field , 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 parameter in attenuation coefficient has a large value from its expression. The exponential function attenuates rapidly if is not too smaller than such as seawater [11]; attenuates to zero within one nanosecond as shown in Figure 2; therefore and attenuate to zero within one nanosecond .

Figure 2: 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 , dry soil , freshwater , and so on. Figure 2 illustrates the duration of about seawater and wet soil for contrast.

Secondly, the item within the sum from n equals zero to in (12) is unharmonious; when is greater than one, the attenuation coefficient may be vertiginous. There is a similar restriction in (4), but the parameter is less than the parameter 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:

Table 1: The effectual depth about attenuation coefficient.

Medium
Freshwater
Wet soil
Dry soil

Effectual depth (m) when = 0°
47.8
16.8
1062

Effectual depth (m) when = 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 to layered medium can be approximated by a double exponential function in time-domain [12]:

where , = 50 (kV/m), α = 4 × 10⁷ (1/s), and β = 6 × 10⁸ (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 = 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 equal . 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 about . 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 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 Figure 8 when 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 Figure 9 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.

Appendix

A. The Derivation Details of

The reflection coefficient in frequency domain is as follows [8]:

The material permittivity and permeability are

The reflection coefficient is represented as follows:

Using the relation , the reflection coefficient is

The transmission coefficient is

Try to normalize the coefficient of under the radical

Define three parameters to simplify the expression:

where , , and . Define and ; then

Use to multiply the denominator and numerator of equation (A.8):

Simplify the expression as follows:

where , , , and .

Define , whose first derivative is

There are some known Laplace Transform pairs in mathematics handbook [13]:

Homologous transformation relations used in this paper are

Combined with the expressions of upper derivation, the final transmission coefficient in time-domain is

B. The Derivation Details of and

There are some known Laplace Transform pairs in mathematics handbook [13]:

There are two approaches to obtain from . The first approach is relatively simple:

The second approach is to use partial integration.

The approach to obtaining is to use partial integration. In partial integral calculation, the derivative of one-order deformation Bessel function is used:

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

C. Christopoulos, “Modeling and simulation for EMC-Part I,” IEEE Electromagnetic Compatibility Magazine, vol. 4, no. 1, pp. 47–56, 2015.

H. P. Neff and D. A. Reed, “Plane Wave (EMP) Incidence on a Finitely Conducting Plane Earth with the Magnetic Field Intensity Parallel to the Earths Surface , Theoretical Note 351,” Air Force Weapons Lab, 1984.

B.-Y. Sun and H. Zhou, “HEMP environment over two-layer lossy soil,” Nuclear Electronics and Detection Technology, vol. 33, no. 5, pp. 572–575, 2013.

D. G. Dudley, T. M. Papazoglou, and R. C. White, “On the interaction of a transient electromagnetic plane wave and a lossy half-space,” Journal of Applied Physics, vol. 45, no. 3, pp. 1171–1175, 1974.

J. A. Klassen, Time-domain Analysis of One-dimensional Electromagnetic Scattering by Lossy Media, Rep. FEL-90-A211, TNO Phys. Electron. Lab., The Hague, 1990.

P. R. Barnes and F. M. Tesche, “On the Direct Calculation of a Transient Plane Wave Reflected from a Finitely Conducting Half Space,” IEEE Transactions on Electromagnetic Compatibility, vol. 33, no. 2, pp. 90–96, 1991.

S. He and S. Ström, “Time-Domain Propagating Modes in a Finitely Conducting Half-Space and Calculation of the Transient Reflection,” IEEE Transactions on Electromagnetic Compatibility, vol. 37, no. 2, pp. 277–282, 1995.

K. Zhang and D. Li, Electromagnetic Theory for Microwaves and Optoelectronics, Publishing House of Electronics Industr, 1998.

B. S. Guru and H. R. Hiziroglu, Electromagnetic Field Theory Fundamentals, China Machine Press, 2nd edition, 2005.

C. R. Paul, Analysis of Multi-conductor Transmission Lines, China electric power press, 2nd edition, 2013.

D. B. Ge and B. Wei, Electromagnetic Wave Theory, Science Publishing Company, Beijing, China, 2011.

Y.-Z. Xie, Z.-J. Wang, Q.-S. Wang, H. Zhou, and B.-Y. Sun, “High altitude nuclear electromagnetic pulse waveform standards: A review,” Qiangjiguang Yu Lizishu/High Power Laser and Particle Beams, vol. 15, no. 8, pp. 781–787, 2003.

W. Hackbusch, H. R. Schwarz, E. Zeidler, and E. Zeidler, Teubner-Taschenbuch der Mathematik, Science Publishing Company, 2003.