# Effects of Reentry Plasma Sheath on Mutual-Coupling Property of Array Antenna

Authors: B. W. Bai; X. P. Li; Y. M. Liu; J. Xu
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

A plasma sheath enveloping a reentry vehicle would cause the failure of on-board antennas, which is an important effect that contributes to the “blackout” problem. The method of replacing the on-board single antenna with the array antennas and using beamforming technology has been proposed to mitigate “blackout” problem by many other researchers. Because the plasma sheath is a reflective medium, plasma will alter the mutual coupling between array elements and degrade the beamforming performance of array antenna. In this paper, the effects of the plasma sheath on the mutual coupling properties between adjacent array elements are studied utilizing the algorithm of finite integration technique. Results show that mutual coupling coefficients of array elements are deteriorating more seriously with the decrease of collision frequency. Moreover, when electron density and collision frequency are both large, plasma sheath improves the mutual coupling property of array elements; this conclusion suggests that replacing the on-board single antenna with the array antennas and using beamforming technology can be adopted to mitigate the blackout problem in this condition.


---

## Body

## 1. Introduction

A spacecraft reentering the Earth’s atmosphere at a hypersonic speed is enveloped by plasma fluid because of the shock wave heating of ambient gas and ablating of heat shield, causing air molecules and heat shield material dissociated and ionized [1]. Plasma fluid is called reentry plasma sheath, and it contains so many free electrons and ions that it attenuates radio frequency signals greatly through the reflection and absorption effects. The primary problem [2] associated with the reentry plasma sheath is the “communications blackout” phenomenon during which performance of on-board electromagnetic systems can be severely degraded or even completely invalid. The problem is of great importance, because the phase of losing tracking, telemetry, and command (TT&C) data and GPS navigation signals up to several minutes is often the most crucial part of the vehicle’s flight. Thus, the problem has been studied since the Gemini and Apollo space programs and attracts more attentions with the increasing space activity.One of many significant effects that contribute to the “blackout” problem is the failure of on-board antenna which is affected greatly by the reentry plasma sheath [3]. Performance degradations and failures of on-board antennas covered by the reentry plasma sheath have been studied by many authors for several decades. Swift et al. [4–6] analysed characteristics of S-band telemetry slot and waveguide antennas radiating into the nonuniform reentry plasma sheath. Galejs [7–9] analysed and measured admittance deviations of plasma-covered slot antenna using the variational method. Based on the impedance sheet notion, Fante [10] described a simple technique for the admittance and radiation pattern calculations of slot antenna covered with thin plasma slab. Hodara [11] analysed the radiation characteristics deviations of slot antennas covered with a magnetized plasma sheath. Qian et al. [12] analysed the performance degradations of C-band microstrip patch antenna covered by nonuniform plasma sheath using PLRC-FDTD algorithm. Recently, Bai [13–15] studies the admittance and radiation degradations of GPS and Beidou navigation patch antenna covered by reentry plasma sheath.It is well known that the degradation of on-board antenna performance would lower the RF signal strength. However, by using beamforming technology, array antennas can provide a certain amount of gain to the communication link between the spacecraft and ground stations, and the received or transmitted RF signal strength would be increased. Thus, replacing the on-board single antenna with the array antennas can be adopted to mitigate the blackout problem. The process of this mitigation method is shown as follows. First, the directions of received or transmitted wave are calculated according to the navigation results and the flight attitude information of the spacecraft. Then, using beamforming technology, on-board array antenna is adjusted to point to the calculated direction. Consequently, the antenna gain and signal strength is increased, and a certain amount of signal gain is obtained to the communication link between the spacecraft and ground station.However, the reentry plasma sheath is a reflective medium and would change the input impedance of array element, alter the mutual coupling between array elements, and degrade the beamforming performance of array antenna [16, 17]. As a result, it is necessary to figure out this issue before this “blackout” mitigating method comes into use. Therefore, the mutual coupling property between the array elements covered by reentry plasma sheath is studied in this paper.The remainder of this paper is organized as follows. Section2 outlines an approximate stratified model of the nonuniform reentry plasma sheath for the EM simulation and utilizes the finite integration technique to analyse the interaction between the reentry plasma sheath and the array antenna. In Section 3, taking S-band telemetry four elements array antenna for example, self-coupling coefficient property of array element and the mutual-coupling coefficient characteristic between the adjacent and diagonal array elements covered by the reentry plasma sheath are calculated and discussed systematically. The conclusions and comments for future study are provided in Section 4.


## 2. Simulation Model

### 2.1. Nonuniform Plasma Sheath Stratified Modeling

In practice, the reentry plasma sheath is spatial nonuniform and variation profiles of electron density can be characterized by a bi-Gaussian function shown below [18]:(1)
Ne
z
=
N
e
peak
exp
⁡
[
-
a
1
z
-
z
0
2
]
(
0
≤
z
<
z
0
)
N
e
peak
exp
⁡
[
-
a
2
z
-
z
0
2
]
z
≥
z
0
.Generally, the nonuniform plasma sheath can be modeled approximately by several adjacent uniform thin plasma slabs according to variation profiles of electron density. A stratified modeling method is described in [15]. The schematic diagram of the stratified model is the dotted line shown in Figure 1. Assuming that the plasma sheath is divided into N layers, the complex relative permittivity of the nth plasma slab can be determined by(2)
ε
~
n
=
1
-
ω
p
,
n
2
ω
2
+
ν
2
-
j
·
ν
ω
·
ω
p
,
n
2
ω
2
+
ν
2
·
ε
0
,where ω
p
,
n is the plasma frequency in the nth layer and ν is the collision frequency. The electromagnetic characteristic of the nonuniform reentry plasma sheath can be established by combination of the plasma frequency ω
p
,
n and the collision frequency ν of each uniform plasma slab.Figure 1
Variation profile of an electron density radially outward from vehicle surface.


### 2.2. Plasma-Covered Array Antenna Simulation Technique

A square patch array antenna of four elements is designed. The array antenna is a simple 2-by-2 array, which is used to provide a representative result that reflective plasma degrades beamforming performance of array antenna in the presence of mutual coupling. The basic structure and dimensions of array element with an off-centre point feeding are shown in Figure2(a). The layout of the array antenna on a plane is shown in Figure 2(b), and the adjacent elements are spaced at a distance of half a wavelength. The array antenna operating frequency is designed to be 2.3 GHz which is the operating frequency of S-band TT&C system. The relative dielectric constant and loss tangent of the substrate are 3 and 0.0013, respectively. The EM simulation model for array antenna covered with nonuniform plasma is shown in Figure 3. The interaction between the reentry plasma sheath and the array antenna is analysed using the algorithm of finite integration technique [19].Figure 2
Basic structure and dimensions of antenna element (a) and layout of the array antenna (b).
(a)
(b)Figure 3
Plasma-covered array antenna simulation models.


## 2.1. Nonuniform Plasma Sheath Stratified Modeling

In practice, the reentry plasma sheath is spatial nonuniform and variation profiles of electron density can be characterized by a bi-Gaussian function shown below [18]:(1)
Ne
z
=
N
e
peak
exp
⁡
[
-
a
1
z
-
z
0
2
]
(
0
≤
z
<
z
0
)
N
e
peak
exp
⁡
[
-
a
2
z
-
z
0
2
]
z
≥
z
0
.Generally, the nonuniform plasma sheath can be modeled approximately by several adjacent uniform thin plasma slabs according to variation profiles of electron density. A stratified modeling method is described in [15]. The schematic diagram of the stratified model is the dotted line shown in Figure 1. Assuming that the plasma sheath is divided into N layers, the complex relative permittivity of the nth plasma slab can be determined by(2)
ε
~
n
=
1
-
ω
p
,
n
2
ω
2
+
ν
2
-
j
·
ν
ω
·
ω
p
,
n
2
ω
2
+
ν
2
·
ε
0
,where ω
p
,
n is the plasma frequency in the nth layer and ν is the collision frequency. The electromagnetic characteristic of the nonuniform reentry plasma sheath can be established by combination of the plasma frequency ω
p
,
n and the collision frequency ν of each uniform plasma slab.Figure 1
Variation profile of an electron density radially outward from vehicle surface.


## 2.2. Plasma-Covered Array Antenna Simulation Technique

A square patch array antenna of four elements is designed. The array antenna is a simple 2-by-2 array, which is used to provide a representative result that reflective plasma degrades beamforming performance of array antenna in the presence of mutual coupling. The basic structure and dimensions of array element with an off-centre point feeding are shown in Figure2(a). The layout of the array antenna on a plane is shown in Figure 2(b), and the adjacent elements are spaced at a distance of half a wavelength. The array antenna operating frequency is designed to be 2.3 GHz which is the operating frequency of S-band TT&C system. The relative dielectric constant and loss tangent of the substrate are 3 and 0.0013, respectively. The EM simulation model for array antenna covered with nonuniform plasma is shown in Figure 3. The interaction between the reentry plasma sheath and the array antenna is analysed using the algorithm of finite integration technique [19].Figure 2
Basic structure and dimensions of antenna element (a) and layout of the array antenna (b).
(a)
(b)Figure 3
Plasma-covered array antenna simulation models.


## 3. Results and Discussion

In this section, reflection coefficient property of an array element and mutual coupling characteristic between adjacent and diagonal array elements covered by reentry plasma sheath of different parameters are studied and calculated.In the calculation, with respect to the typical reentry plasma generated in the RAM-C project [20], the bi-Gaussian plasma density distribution parameters range is confined as N
e
peak~1016/m3~1018/m3, ν~0.1 GHz~20 GHz and the plasma sheath thickness is set as 10 cm. In order to simplify the multivariable problem, the growth index a
1 and the fall-off index a
2 are both selected to be 0.5, and the position of N
e
peak is fixed at the centre of the reentry plasma sheath. The simulation results are shown below.


### 3.1. Validity Analysis

To check the validity of this simulation method, we carry out a comparison between results using the simulation method and results using XFDTD simulation tool. In this validation, we give some typical comparison results, and the plasma parametersNe and ν are set as 1016/m3 and 0.1 GHz, 1017/m3 and 1 GHz, and 1018/m3 and 5 GHz, respectively. The material for plasma is determined by the Debye-Drude dispersion relation model, and the layered plasma model is also applied to model the nonuniform plasma sheath In the XFDTD simulation. The S
21 parameter which is mutual coupling coefficient from the number 1 array element to the number 2 array element is presented. As shown in Figure 4, the solid lines are the results from this paper simulation method, and the scatters are the results from the XFDTD simulation. It is shown that the results from different simulation methods are coincident well with each other. Thus, the simulation method in this paper is effective and corrective.Figure 4
Comparisons between the results from the analysis method and the XFDTD simulation method.


### 3.2. Self-Coupling Coefficient of Array Element

Since the plasma is a reflective medium, the beamforming of array antenna would be affected by the wave reflected back into its array element. We need to calculate the self-coupling coefficient of array element.Because the structure of every array element is the same, theS
11 parameter which is the self-coupling coefficient of number 1 array element is presented as a demonstration. Setting N
e
peak as 0, 1016/m3, 1017/m3, and 1018/m3 for each calculation, self-coupling coefficients for collision frequencies 0.1 GHz, 1 GHz, 5 GHz, and 20 GHz are shown in Figures 5(a)~5(d), respectively.Figure 5
(a) Self-coupling coefficient (ν
=
0.1 GHz) of array element. (b) Self-coupling coefficient (ν
=
1 GHz) of array element. (c) Self-coupling coefficient (ν
=
5 GHz) of array element. (d) Self-coupling coefficient (ν
=
20 GHz) of array element.
(a)
(b)
(c)
(d)The results indicate that self-coupling coefficients of the array element are deteriorating more seriously with the increase of electron density. For instance, when the collision frequency is 1 GHz, the S
11 values for N
e
peak 1016/m3, 1017/m3, and 1018/m3 are −24.5 dB, −10.2 dB, and −2.6 dB, respectively.Besides, the lower the collision frequency, the more serious the deterioration of self-coupling coefficient. For instance, when the N
e
peak is 1017/m3, the S
11 values for the collision frequencies 0.1 GHz, 1 GHz, and 5 GHz are −6.5 dB, −10.2 dB, and −19 dB, respectively.Moreover, for a large electron density (e.g., 1018/m3), the resonant frequency of the array element is adjusted to a higher frequency.


### 3.3. Mutual Coupling Coefficient of Adjacent Array Elements

In the former calculation, the results reveal that the plasma sheath reflects EM waves back into the array antenna. Thus, the mutual coupling characteristic of adjacent elements caused by the plasma sheath is studied in this calculation.TheS
21 parameter which is mutual coupling coefficient of the number 1 array element to the number 2 array element is presented as a demonstration for mutual coupling properties of adjacent elements. By setting plasma parameters the same as above, the mutual coupling coefficients are shown in Figures 6(a)~6(d).Figure 6
(a) Mutual coupling coefficient (ν
=
0.1 GHz) of adjacent elements. (b) Mutual coupling coefficient (ν
=
1 GHz) of adjacent elements. (c) Mutual coupling coefficient (ν
=
5 GHz) of adjacent elements. (d) Mutual coupling coefficient (ν
=
20 GHz) of adjacent elements.
(a)
(b)
(c)
(d)The results indicate that the magnitude of mutual coupling coefficients is deteriorating more seriously with the decrease of collision frequency. In particular, for the collision frequency lower than the antenna operating frequency (Figures5(a) and 5(b)), the mutual coupling coefficients are increased when the N
e
peak are 1016/m3 and 1017/m3, whereas the mutual coupling coefficients are decreased when the N
e
peak is 1018/m3. For the collision frequency higher than the antenna frequency (Figures 5(c) and 5(d)), the mutual coupling coefficients change a little when the N
e
peak is 1016/m3 and 1017/m3, whereas the mutual coupling coefficients are decreased when the N
e
peak is 1018/m3.The results also indicate that the phase variation relative to initial phase of mutual coupling coefficients turn bigger with the increase ofN
e
peak and decrease of collision frequency.


### 3.4. Mutual Coupling Coefficient of Diagonal Array Element

TheS
31 parameter which is mutual coupling coefficient of the number array element to the number 3 array element is presented as a demonstration for mutual coupling properties of diagonal adjacent elements. By setting plasma parameters the same as above, the mutual coupling coefficients are shown in Figures 7(a)~7(d), respectively.Figure 7
(a) Mutual coupling coefficient (ν
=
0.1 GHz) of diagonal elements. (b) Mutual coupling coefficient (ν
=
1 GHz) of diagonal elements. (c) Mutual coupling coefficient (ν
=
5 GHz) of diagonal elements. (d) Mutual coupling coefficient (ν
=
20 GHz) of diagonal elements.
(a)
(b)
(c)
(d)The results indicate that when theN
e
peak is 1018/m3 and the collision frequency is 0.1 GHz, the plasma sheath deteriorates the magnitude of mutual coupling coefficients of nonadjacent elements. As for the other plasma parameters, the mutual coupling coefficients are decreased. We can conclude that the plasma sheath improves the mutual coupling property of diagonal array elements.


## 3.1. Validity Analysis

To check the validity of this simulation method, we carry out a comparison between results using the simulation method and results using XFDTD simulation tool. In this validation, we give some typical comparison results, and the plasma parametersNe and ν are set as 1016/m3 and 0.1 GHz, 1017/m3 and 1 GHz, and 1018/m3 and 5 GHz, respectively. The material for plasma is determined by the Debye-Drude dispersion relation model, and the layered plasma model is also applied to model the nonuniform plasma sheath In the XFDTD simulation. The S
21 parameter which is mutual coupling coefficient from the number 1 array element to the number 2 array element is presented. As shown in Figure 4, the solid lines are the results from this paper simulation method, and the scatters are the results from the XFDTD simulation. It is shown that the results from different simulation methods are coincident well with each other. Thus, the simulation method in this paper is effective and corrective.Figure 4
Comparisons between the results from the analysis method and the XFDTD simulation method.


## 3.2. Self-Coupling Coefficient of Array Element

Since the plasma is a reflective medium, the beamforming of array antenna would be affected by the wave reflected back into its array element. We need to calculate the self-coupling coefficient of array element.Because the structure of every array element is the same, theS
11 parameter which is the self-coupling coefficient of number 1 array element is presented as a demonstration. Setting N
e
peak as 0, 1016/m3, 1017/m3, and 1018/m3 for each calculation, self-coupling coefficients for collision frequencies 0.1 GHz, 1 GHz, 5 GHz, and 20 GHz are shown in Figures 5(a)~5(d), respectively.Figure 5
(a) Self-coupling coefficient (ν
=
0.1 GHz) of array element. (b) Self-coupling coefficient (ν
=
1 GHz) of array element. (c) Self-coupling coefficient (ν
=
5 GHz) of array element. (d) Self-coupling coefficient (ν
=
20 GHz) of array element.
(a)
(b)
(c)
(d)The results indicate that self-coupling coefficients of the array element are deteriorating more seriously with the increase of electron density. For instance, when the collision frequency is 1 GHz, the S
11 values for N
e
peak 1016/m3, 1017/m3, and 1018/m3 are −24.5 dB, −10.2 dB, and −2.6 dB, respectively.Besides, the lower the collision frequency, the more serious the deterioration of self-coupling coefficient. For instance, when the N
e
peak is 1017/m3, the S
11 values for the collision frequencies 0.1 GHz, 1 GHz, and 5 GHz are −6.5 dB, −10.2 dB, and −19 dB, respectively.Moreover, for a large electron density (e.g., 1018/m3), the resonant frequency of the array element is adjusted to a higher frequency.


## 3.3. Mutual Coupling Coefficient of Adjacent Array Elements

In the former calculation, the results reveal that the plasma sheath reflects EM waves back into the array antenna. Thus, the mutual coupling characteristic of adjacent elements caused by the plasma sheath is studied in this calculation.TheS
21 parameter which is mutual coupling coefficient of the number 1 array element to the number 2 array element is presented as a demonstration for mutual coupling properties of adjacent elements. By setting plasma parameters the same as above, the mutual coupling coefficients are shown in Figures 6(a)~6(d).Figure 6
(a) Mutual coupling coefficient (ν
=
0.1 GHz) of adjacent elements. (b) Mutual coupling coefficient (ν
=
1 GHz) of adjacent elements. (c) Mutual coupling coefficient (ν
=
5 GHz) of adjacent elements. (d) Mutual coupling coefficient (ν
=
20 GHz) of adjacent elements.
(a)
(b)
(c)
(d)The results indicate that the magnitude of mutual coupling coefficients is deteriorating more seriously with the decrease of collision frequency. In particular, for the collision frequency lower than the antenna operating frequency (Figures5(a) and 5(b)), the mutual coupling coefficients are increased when the N
e
peak are 1016/m3 and 1017/m3, whereas the mutual coupling coefficients are decreased when the N
e
peak is 1018/m3. For the collision frequency higher than the antenna frequency (Figures 5(c) and 5(d)), the mutual coupling coefficients change a little when the N
e
peak is 1016/m3 and 1017/m3, whereas the mutual coupling coefficients are decreased when the N
e
peak is 1018/m3.The results also indicate that the phase variation relative to initial phase of mutual coupling coefficients turn bigger with the increase ofN
e
peak and decrease of collision frequency.


## 3.4. Mutual Coupling Coefficient of Diagonal Array Element

TheS
31 parameter which is mutual coupling coefficient of the number array element to the number 3 array element is presented as a demonstration for mutual coupling properties of diagonal adjacent elements. By setting plasma parameters the same as above, the mutual coupling coefficients are shown in Figures 7(a)~7(d), respectively.Figure 7
(a) Mutual coupling coefficient (ν
=
0.1 GHz) of diagonal elements. (b) Mutual coupling coefficient (ν
=
1 GHz) of diagonal elements. (c) Mutual coupling coefficient (ν
=
5 GHz) of diagonal elements. (d) Mutual coupling coefficient (ν
=
20 GHz) of diagonal elements.
(a)
(b)
(c)
(d)The results indicate that when theN
e
peak is 1018/m3 and the collision frequency is 0.1 GHz, the plasma sheath deteriorates the magnitude of mutual coupling coefficients of nonadjacent elements. As for the other plasma parameters, the mutual coupling coefficients are decreased. We can conclude that the plasma sheath improves the mutual coupling property of diagonal array elements.


## 4. Conclusions

This paper presents an approximate layered EM simulation model for the nonuniform reentry plasma sheath which obeys bi-Gaussian distributions from NASA’s report and utilizes the algorithm of finite integration technique to analyze the effects of the plasma sheath on the reflection coefficient property of an array element and mutual coupling property between array elements. The results indicate that reflection coefficients of an array element are deteriorating more seriously with increase of the electron density and decrease of the collision frequency. The resonant frequency of the array element is adjusted to a higher frequency by a larger electron density. The mutual coupling coefficients of array elements are deteriorating more seriously with the decrease of collision frequency. Moreover, when electron density and collision frequency are both large, plasma sheath improves the mutual coupling property of array elements, this conclusion suggests that replacing the on-board single antenna with the array antennas and using beamforming technology can be adopted to mitigate the blackout problem in this condition.



---

*Source: 542392-2015-04-07.xml*