† Corresponding author. E-mail:
Project supported by the National Basic Research Program of China (Grant No. 2014CB340203), the National Natural Science Foundation of China (Grant Nos. 61431010 and 61501350), and the Natural Science Foundation of Shaanxi Province, China (Grant Nos. 2018JM6016 and 2016JM1001).
The propagation characteristics of oblique incidence terahertz (THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain (SO-FDTD) method combined with the phase matching condition. The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab. Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions, which makes it a potential candidate to overcome the ionization blackout problem.
When vehicles fly at hypersonic velocities within the atmosphere, they become enveloped in a plasma sheath that prevents radio communication between the aerospace and the ground control center, which is the well-known blackout problem.[1–4] In order to mitigate or eliminate the reentry blackout problem, a number of techniques have been proposed, such as high power, high frequencies, low frequencies, lasers, magnetic fields, and so on.[3,5–7] The terahertz wave (THz), commonly referred to as T-rays, sub-millimeter, or far infrared, refers to the electromagnetic radiation at the frequency from 0.1 THz to 10 THz, has attracted much attention recently partially because it is one of the potential candidates in the high frequencies communication methods. The techniques for the generation of THz wave were developed rapidly in recent years[8–11] and promising applications of THz waves can be found in various fields.[12] As the basis of communication between the hypersonic vehicles and the ground control center, the interaction between plasma and THz wave (absorption, reflection, and transmission) is a key issue to be analyzed. Over the past few years, the propagation characteristics of the THz wave in plasma have been investigated by different methodologies, such as the analytical method,[13–17] experimental method,[18] scattering matrix method (SMM),[19,20] propagation matrix method (PMM),[21] finite-difference time-domain (FDTD) method,[22–28] Wentzel–Kramer–Brillouin (WKB) method,[29,30] and so on.
The FDTD method has been widely used to solve the interaction problems between the THz wave and plasma because of its capability in simulating the propagation characteristics of a wideband THz wave through inhomogeneous dispersive plasma. Wang et al.[22,27,28] investigated the propagation properties of the THz wave in a time-varying dusty plasma slab using the auxiliary differential equation finite-difference time-domain (ADE-FDTD) method. Zheng et al.[23] studied the propagation characteristics of the THz wave in unmagnetized plasma by the Z transform finite-difference time-domain (ZT-FDTD) method. Theoretical results were compared with experimental measurements in Ref. [20], where good agreements were achieved. The broadband propagation characteristics of the THz wave in an anisotropic magnetized plasma were investigated by Tian et al.[24] using the JE convolution FDTD (JEC-FDTD) method. Gao et al.[25] investigated the characteristics of reflection, transmission, and absorption of THz waves in a sandwich type microplasma structure by the ZT-FDTD method.
Although studies of propagation characteristics of the THz wave in plasma media can be found in several works, most of the previous research focused on the normally incident wave cases, i.e., the incident plane wave propagates at an angle of zero with respect to the normal of the media interfaces, limited attention was paid to the oblique incidence cases. In practice, the angle of incidence changes over a wide range during flight, which indicates that the oblique incidence wave cases require more consideration. Recently, using the impedance matching method, Rahmani et al.[31] investigated the transmission, reflection, and absorption of s-polarized wave obliquely incident on the plasma slab which has a bell-like electron density distribution. The SMM was adopted by Yuan et al.[32] to analyze the transmission of obliquely incident THz waves in the plasma sheaths covering a blunt coned vehicle, in which the plasma sheaths were generated based on a hypersonic fluid model. Cao et al.[33] studied the propagation characteristics of obliquely incident THz waves though dusty plasma using the PMM, where the electron density distribution of the dusty plasma was assumed to be parabolic. Based on the twice reflection model, Tian et al.[34] investigated the obliquely incident wave through a nonuniform magnetized plasma slab, in which both the electron density and the collision frequency across the plasma were assumed to be in a Gaussian profile. To solve the problem under study, the non-uniform plasma was divided into a series of layers with uniform electron density.[34] Nevertheless, this dividing method fails when one encounters plasma with a high gradient electron density distribution. To overcome this difficulty, the FDTD is a good choice.
The THz wave obliquely incident on a plasma sheath, which is non-uniform in half space, is a two-dimensional problem for the FDTD method. Thus it will cost a lot of memory and computation time in the simulation. To save computational memory and time, in this paper, the FDTD method combined with the phase matching condition is utilized to solve the propagation problem of obliquely incident THz waves through an inhomogeneous plasma. With the help of phase matching condition, the problem is simplified to a pseudo one-dimensional problem, which can greatly reduce the memory required and improve the efficiency. To describe the distribution of the electron density in an inhomogeneous plasma, the Gaussian profile is introduced which corresponds to the flow field around the hypersonic aircraft at a high altitude of space.[35]
The rest of this paper is structured as follows. In Section
A geometric sketch of the problem under study is shown in Fig.
As shown in Fig.
The fields in plasma can be described by Maxwell’s equations in the frequency domain
For the TEz mode, the propagation equation in the xoy plane can be expressed as
For simplicity, two new auxiliary variables
The permittivity of unmagnetized plasma is given as[36]
Substituting Eq. (
With the transition relationship from the frequency domain to time domain, jω → ∂/∂t, we can obtain the expressive forms of Eqs. (
If we assume that the function in the time domain has the form
The expression of
For the TM case, the transmission coefficient and absorption coefficient can be calculated in the same way by substituting electric field Ez for magnetic field Hz.
Numerical simulations are implemented based on the theoretical treatments presented in Section
Before we start to analyze the propagation characteristics of the THz wave through an inhomogeneous plasma, the validity of the method adopted in this paper is verified by comparing the results obtained by our method with those obtained by an analytical method.[35] The propagation characteristics of THz waves with different incident angles in a homogeneous plasma slab are taken into consideration.
The FDTD problem space consists of 474 cells, and the homogeneous plasma occupies 71–404 cells. The cell size is Δy = 1.5 × 10−5 m, and the thickness of the plasma layer is 5~mm. Both ends of the cell space are treated by uniaxial perfectly matched layer absorbing boundary condition to eliminate unwanted reflections. The plasma parameters are
Comparisons of the transmission coefficients of the THz wave in the plasma slab calculated by our method and those obtained using the analytical solution are displayed in Fig.
The plasma density stemmed from the reentry hypersonic flight in the stagnation region could reach 1021 m−3 and the plasma frequency corresponding to this electron density can reach 0.28 THz.[24] In the following simulations, the Gaussian electron density profile is chosen to stand for the variation of electron density, which corresponds to the non-uniform plasma flow field around the hypersonic aircraft at a high altitude of space[35]
In order to investigate the influence of the incident angle on the propagation characteristics of a THz wave in plasma, the incident angle of the THz wave is varied in steps of 15° from 0° to 90° in the simulations. The transmission and absorption coefficients for TE wave are calculated and shown in Fig.
Comparisons of transmission and absorption coefficients for plasma slab of different thickness are displayed in Fig.
Figure
The variations of power transmission and absorption coefficients for THz wave obliquely incident on a non-uniform plasma slab are analyzed by the SO-FDTD method combined with the phase matching condition. The factors that may affect the propagation characteristics of THz wave through plasma, including the incident angle of the incident wave, the thickness of the plasma, and the polarization of the wave, are studied. The results reveal that the incident angle can affect the propagation properties greatly, while the polarization of the incident wave has little influence on the propagation properties for the cases under study. When the THz wave is employed to the communication in the ionization blackout situation, the incident angle should be considered, and the frequency corresponding to the maximum absorption should be avoided. This study will be useful to provide a reference for the exploration of potential techniques to overcome ionization blackout problem as well as to provide insights to the interactions between THz wave and plasma.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] |