Inverse synthetic aperture radar range profile compensation of plasma-sheathenveloped reentry object

Yaocong XIE(謝曜聰),Xiaoping LI(李小平),Fangfang SHEN(沈方芳),Bowen BAI(白博文),Lei SHI(石磊)and Xuyang CHEN(陳旭陽)

School of Aerospace Science and Technology,Xidian University,Xi’an 710071,People’s Republic of China

Abstract The scattering points in a plasma sheath characterized with coupled velocities can cause pulse compression mismatching,which results in displacement and energy diffusion in the onedimension range profile.To solve this problem,we deduce the echo model of the plasma-sheathenveloped reentry object.By estimating the coupled velocities,we propose a compensation method to correct the defocus of an inverse synthetic aperture radar(ISAR)image in range dimension to improve the quality of the ISAR images.The simulation results suggest that the echoes from different regions of the surface of the reentry object have various coupling velocities,and the higher the coupled velocity,the more serious the displacement and energy diffusion in the range dimension.Our proposed method can correct the range dimension aberration.Two measurement metrics were used to evaluate the improvement of the compensation method.

Keywords:plasma sheath,inverse synthetic aperture radar imaging,defocus,coupled velocity,reentry object

1.Introduction

Inverse synthetic aperture radar(ISAR)imaging is an effective technique by which images of moving targets in the range-Doppler domain can be obtained[1,2].Two-dimensional high-resolution images of the reentry object can be obtained by ISAR imaging for subsequent operations of target recognition,which is of great significance to the safe recovery of spacecraft.When a reentry spacecraft enters into the Earth’s near-space atmosphere,the effect of the generated aerodynamic heat can form a plasma sheath covering the surface of the reentry object[3–5],absorbing,attenuating and reflecting electromagnetic(EM)waves[6–8].

For the ISAR imaging of high-speed targets,research on the compensation method of the moving-in-pulse duration(MPD)has attracted much attention.Fenget alproposed an effective algorithm based on Radon transform for the range profile of a hypersonic target,through which the range profile aberration and ISAR image blurring can be compensated[9].Tianet alstudied the integrated cubic phase function and estimated the second-order phase coefficient to compensate for the high-speed movement on ISAR images[10].However,these studies mainly focused on the influence of the target’s high velocity on ISAR defocusing and ignored the problem of the generated plasma sheath enveloping the target during the flight.

In addition to the above problem,most of the related studies focused on the influence of the plasma sheath imposed on radar detection.Yuanet alanalyzed the relationship between the plasma parameters and the reflection that is a combination of the power reflected from the plasma interface.They derived that the reflections of the incident wave power depended on the thickness of plasma for different incident wave frequencies and various plasma parameters[11].Dinget aldiscovered that the velocity of the plasma sheath had an effect on the one-dimensional range profile,which caused a false target and even detection failures[12].In their approach,the reentry object was treated as a single scattering point rather than multiple scattering points,which is a prerequisite for high-resolution ISAR imaging.Zhanget alproposed a wideband radar scattering center model of a plasma-sheathcovered reentry target,in which the reflection characteristics and velocity of the plasma sheath are considered[13].Despite the fact that their research proposed a basic model,they did not conduct a profound analysis of the problems caused by the plasma sheath.Xuet alpresented the MPD model for hypersonic targets and evaluated the performances of pulse compression under the stop-and-go and MPD models by theoretical analysis and simulations[14].However,the effect of plasma was not considered in their study.

Unfortunately,few studies concerning the ISAR defocusing of the reentry object can be found so far.This work analyzes the reflection characteristics of the plasma sheath to obtain the radar echo model of the plasma-sheath-enveloped reentry object.It is found that in the areas with high electron density,the EM wave cannot penetrate through the plasma sheath,and the radar scattering points are distributed at a certain depth in the plasma sheath.These scattering points are affected by the flow field of the plasma sheath,exhibiting coupled velocities.The coupled velocities will further introduce linear frequency modulation(LFM)in the pulse compression results,causing displacement and energy diffusion in the one-dimension range profile.Because the coupled velocities of scattering points differ in the plasma sheath,certain regions of the ISAR image will be distorted and defocused to some extent in the range dimension.By estimating the coupled velocity,we propose a compensation method to correct the range dimension aberration.The simulation results and two measurement metrics have verified the effectiveness of our proposed method.

The remainder of this paper is organized as follows.Section 2 introduces the reflection coefficient of the plasma sheath and the reflection depths with different electron densities.Section 3 deduces the echo model of the plasma sheath and discusses the formation mechanism of displacement and energy diffusion in the range dimension.Section 4 proposes the compensation method.Section 5 simulates the results of ISAR imaging using the compensation method,and evaluates the improvement effect.Section 6 summarizes the conclusions of this paper.

2.Reflection coefficient and reflection depth of plasma sheath

The plasma is an ionized medium that can absorb,attenuate and reflect EM waves.Functioning as the most important parameter of plasma,the reflection characteristics of plasma primarily correlate with the electron density[15].Currently,the most commonly used method to analyze the parameters of the plasma sheath is numerical simulation of computational fluid dynamics(CFD)[16,17].In this work,we adopt the simulation data from the CFD of RAM-C II to establish the blunt cone model[18,19].In the study of the RAM-C II project,the altitude of the reentry object ranges from 20 to 80 km.At high altitude,the air density and temperature are lower,so the plasma has lower electron density,and can be penetrated by the EM wave.At low altitude,the air density and temperature are higher,so the plasma has a higher electron density that is difficult for the EM wave to penetrate through.In this study,we choose the reentry object with different velocities at the altitude of 50 km,which is an intermediate transition state.The reentry object at 50 km altitude has the plasma characteristics of both low electron density and high electron density when the velocity is changed.For the application scenario of analyzing the displacement and energy diffusion caused by the coupled velocity of the plasma sheath,50 km is a typical value.

Figure 1(a)shows the electron density distributions of the plasma sheath enveloping the surface of a blunt cone reentry object with 15 Mach velocity and 50 km altitude.It can be seen that the electron density varies by several orders of magnitude at different positions and depths of the plasma sheath,ranging from 1 × 1017to 1 × 1020m-3.Furthermore,along the surface of the reentry object,the closer to the stagnation point,the higher the electron density.Along the vertical direction of the surface,a higher electron density appears closer to the surface.According to the research of the RAM-C II project,a non-uniform plasma sheath can be generally stratified into a model of evenly distributed plasma with multiple layers,and the multiple-layered model of plasma has been adopted in several related studies[12,13,20].Figure 1(b)is a schematic diagram showing the uniform plasma sheath withNlayers at the local part of the spacecraft surface.

Figure 1.Electron density distribution(a)and stratification model(b)of the plasma sheath.

Figure 2.Schematic diagram of transmission line method.

Figure 3.Amplitude of reflection coefficient with different depths where the peak electron density varies from 1 × 1017to 1 × 10 20 m-3.

Figure 4.Relative velocity distribution and schematic diagram of the coupled velocities of scattering points.

Figure 5.Flowchart of the compensation method.

The propagation characteristics of EM waves in the plasma sheath are similar to those of EM waves in a lossy microwave transmission line,and can be written as a cascade of impedance of different characteristic waves.The reflection coefficient of the stratified plasma sheath can be calculated by the transmission line method(TLM)[21,22].A schematic diagram of the TLM is presented in figure 2.With respect to thenth layer plasma,we have

whereknandZnare the propagation constant and the effective impedance of the plasma,fis the frequency of the incident wave,μ0is the vacuum permeability,andneis the complex permittivity of thenth layer.necan be expressed as

whereeis the unit charge,meis the electron mass,0eis the vacuum permittivity,Veis the collision frequency correlating with the temperature and pressure of the plasma,andNe,nis the electron density of thenth layer.

Supposing thatnθis the incident angle of thenth layer,the transmission matrix of thenth layer plasma sheath can be expressed as

When EM waves penetrate through the plasma sheath,the reflection comes from both the plasma sheath and the reentry object surface.For ease of calculation,we presuppose that the surface is made of metal.Then the reflection coefficient of the plasma-sheath-enveloped surface can be obtained using the following equation:

When the frequency of the incident EM waves is lower than the plasma frequency,the EM waves will be attenuated and cannot penetrate through the plasma[13].Considering that the electron density distribution of the plasma has a double-Gaussian layered structure on the surface of the highspeed target[20],the outer layer electron densities of the plasma are low,so the EM wave can penetrate through them.As the incident depth increases,the inner layer electron densities become higher.When the plasma frequency of thenth layer is greater than the incident wave frequency,the incident depth of the EM wave cannot increase at the position of thenth layer,and the position of thenth layer is the incident depth of the EM wave.When the thickness of thenth layer is low,the actual incident depth can be deeper.

Through analyzing the electron density distribution of the plasma sheath on the surface of a blunt cone reentry object and calculating the reflection coefficient,we conclude that the radar wave cannot penetrate through the high-electron-density area of the plasma sheath and can be totally reflected at a certain depth.For areas with low electron density,the radar wave can penetrate through the plasma sheath and be reflected by the surface of the reentry object.

3.Echo model of plasma sheath and mechanism of ISAR defocus

Under the influence of the viscous resistance of the air and reentry object surface,the flow field of plasma has relative flow with the surface.The flow field at different depths has different flow velocities,and the velocity exhibits the characteristics of spatial three-dimensional gradient distribution,where the different colors shown in figure 4(a)denote the relative velocity between the plasma sheath and the blunt cone reentry object with 15 Mach velocity and 50 km altitude.The closer to the surface of the reentry object,the lower the relative velocity.

Supposing the relative velocity of thekth scattering point at a certain depth in the plasma sheath isu k,and the moving velocity of the reentry object isv,then the moving velocity of thekth scattering point can be written as

After performing inverse Fourier transform,the onedimensional range profile of the plasma sheath echoes can be expressed as

4.Compensation method in range dimension

The coupled velocity of the scattering points is the main cause of pulse compression mismatching.Because of the coupled velocity,the linear term is changed and the quadratic term in the pulse compression result cannot be eliminated in equation(17)when processing the matched filtering,which further causes displacement and energy diffusion in the range dimension of the ISAR image.Therefore,it is necessary to reconstruct a modified matched filter to correct the linear term and to eliminate the quadratic term.The modified matched filter can be expressed as

Then,by conducting pulse compression to equation(15)using the modified matched filter,we obtain

After performing inverse Fourier transform,the onedimensional range profile of the plasma sheath echoes processed can be expressed as

Step 4.Perform FRFT with respect to the pulse compression data to obtain the result in the frequency domain,from which the rotation angle and the frequency-domain data corresponding to the maximum peak are extracted.Obtain the coupling velocity of the strongest scattering point by its linear frequency modulation.The one-dimensional range profile of the strongest scattering point in the time domain can be filtered by inverse FRFT transform.In this way,every scattering point can be filtered,and the compensated onedimensional range profile can be obtained respectively.Then translational motion compensation(TMC)is used to realize envelope alignment and phase focus[24],and we can obtain the focused ISAR image by using the range-Doppler(RD)imaging algorithm[25,26].

The flowchart of the compensation method is shown in figure 5.

5.Simulation

In order to evaluate the improvement effect of our proposed method,we simulate the displacement and energy diffusion in the one-dimensional range profile using both the traditionally matched filtering and the modified matched filtering methods.Using the CFD’s flow-field simulation data of the RAM-C II,we simulate the defocus phenomenon and the results of the compensation method in the one-dimensional range profile and in the two-dimensional ISAR image,respectively.The effectiveness of the proposed method is verified by the peak signal-to-noise ratio(PSNR)and structural simulation index method(SSIM).

5.1.Simulation results of displacement and energy diffusion at different velocities

In general,the velocity of the reentry object ranges from 10 to 30 Mach.Considering the relative velocity between the plasma sheath and the reentry object,the coupled velocity of the plasma sheath is less than or equal to that of the reentry object.Supposing that the altitude of the object is 50 km,the carrier frequency is 10 GHz,and the velocity of reentry object is 25 Mach,the displacement and energy diffusion with different bandwidths and pulse widths can be calculated by equations(19)and(21).The results are shown in figure 6.The increase in bandwidth will weaken the displacement in range dimension,and the increase in pulse width will exacerbate the displacement and energy diffusion.

In order to simulate the effect of our proposed method,the pulse width is set at 100 μs and the bandwidth is set at 1 GHz,which are the typical parameters for an ISAR system.We calculate the displacement and energy diffusion at different coupled velocities using both the traditional and the modified pulse compression processing methods.The simulation results are shown in figure 7.We select three coupled velocity estimation errors(CVEEs),which are 1%,5% and 10%,respectively,to calculate the improvement under different conditions.The error of estimated echo delay has been considered.Figure 7(a)shows the displacement caused by the coupled velocity.The higher the coupled velocity,the greater the displacement.The proposed method significantly reduces the displacement in the one-dimensional range profile,outperforming the effect of the traditional method.The increase in CVEE slightly enlarges the displacement.In figure 7(b),the proposed method almost eliminates the effect of energy diffusion,whereas the result of the traditional method exhibits apparent energy diffusion.The deterioration caused by CVEE does not exceed the resolution cell.

Figure 6.Displacement and energy diffusion with different bandwidths and pulse widths.

Figure 7.Displacement and energy diffusion at different coupled velocities using the traditional and the proposed method.

The displacement and energy diffusion of the scattering point in the range dimension caused by the plasma sheath can be calculated using equations(19)and(21).We can obtain the range resolution by the bandwidth.Combining with the analysis results in figure 7,the correction error caused by CVEE can be obtained,which are converted to the number of range cells.For a 15 Mach target,the displacement number is one range cell.For a 20 Mach target,the displacement is two range cells,and for a 25 Mach target,it is six range cells.The energy diffusion is less than one range cell for these three different velocities.

5.2.Simulation results of the compensation method for onedimensional range profile

We select three scattering points with typical coupled velocities to conduct the simulation and to verify the effectiveness of our proposed compensation method.The distance between the scattering points in range dimension is 2.5 m.The modulation type of the radar signal is LFM.The pulse width is 100 μs,the bandwidth is 1 GHz,and the signal-to-noise ratio(SNR)of the echo is ?15 dB.Figure 8 shows the onedimensional range profiles of three scattering points,which are obtained using the traditional method and the compensation method,respectively.The black line in figure 8 represents the reference one-dimensional range profile of the three scattering points without considering the plasma.The blue line represents the one-dimensional range profile of the three scattering points in plasma.The coupled velocities of the three scattering points are set as 5,15 and 25 Mach,respectively,which are within the velocity range of the reentry object[5].Figure 8(a)demonstrates the result obtained using the traditional method,which is calculated using equation(17).The three scattering points have the problem of displacement and energy diffusion in the one-dimensional range profile,to which higher coupled velocity implies greater displacement and diffusion.Figures 8(b)–(d)represent the results using the compensation method with the three estimated velocities.It can be observed that each estimated velocity can compensate for the scattering points of the corresponding velocity to the reference position after eliminating its energy diffusion.Figure 8(e)shows the filtered and merged results of the focused scattering points,in which the displacement and energy diffusion of these scattering points are corrected.Figures 8(e)and(f)show the focused results when the SNR is ?20 dB and ?25 dB.It can be found that our compensation method performs well in low-SNR scenarios.

5.3.Simulation results of the compensation method for ISAR imaging

In this subsection,we analyze the simulation results of ISAR imaging,which are obtained using the compensation method.The pulse repetition frequency is 1 kHz and the SNR of the echo is ?20 dB.We select nine regions on the surface of the blunt cone model from the CFD’s flow-field simulation data from RAM-C II,which is shown in figure 9.These nine regions include the stagnation point,the middle part and the tail of the reentry object.

We simulate the flight state of the reentry object in three typical scenarios,in which the velocities of the object are set as 15,20 and 25 Mach,respectively.The flight altitude is set as 50 km and the observing angle of the radar is 15°.The amplitudeAkand the phase shiftkφcaused by the plasma can be calculated by equations(6)and(7).The peak electron densityNe,the coupled velocityand other parameters originate from the CFD flow-field simulation data of RAM-C II[18,19].According to the amplitudes of the reflection coefficient at different depths,the depth of the scattering point can be obtained by finding the position of maximum amplitude along the thickness direction.The plasma flow-field data covering the nine regions of the surface are listed in tables 1–3,respectively.The fluctuation range of electron density is usually several orders of magnitude,and the variation range is much greater than the collision frequency.Although both the electron density and collision frequency are closely correlated with the plasma characteristics,the characteristic difference caused by the wider variation range of electron density is much greater than that caused by the collision frequency.The use of fixed collision frequency does not affect the final analysis results.For ease of calculation,we adopt a fixed collision frequency of 1 GHz as the parameter of plasma.

As can be seen from the data in the tables,when the velocities of the reentry object are 15 and 20 Mach,the electron density at the tail of the object is lower,and the radar wave can penetrate through the plasma sheath.The scattering points of these regions are distributed on the surface of the object,and their coupled velocities are identical to the velocity of the reentry object,as shown by the 4th to 7th regions in table 1 and by the 5th to 6th regions in table 2.

Figure 10(a)shows the normal ISAR image of the reentry object as the reference image.Figures 10(b)–(d)are the ISAR images of the plasma-sheath-enveloped reentry object with the velocities of 15,20 and 25 Mach,respectively,which are obtained using the traditional processing methods.The radar wave penetrates through the plasma sheath at the tail of the object with velocities of 15 and 20 Mach.Therefore,the most serious displacement occurs at the tail of the object due to the maximum coupled velocity,as shown in figures 10(b)and(c).Meanwhile,for the regions where EM waves cannot penetrate,the displacement of scattering points is weak,because their coupled velocities are lower than the velocity of the reentry object.This suggests that the higher the coupled velocity,the more serious the displacement.

Figure 8.One-dimensional range profiles of three scattering points corrected using the compensation method.

Figure 9.Schematic diagram of the reentry object.

Table 1.Parameters of plasma sheath covering the nine regions with velocity of 15 Mach.

Table 2.Parameters of plasma sheath covering the nine regions with velocity of 20 Mach.

Table 3.Parameters of plasma sheath covering the nine regions with velocity of 25 Mach.

We use FRFT to obtain the evaluated chirp rateμ?kto calculate the estimated velocity ?vk.The results are shown in table 4.

Figures 11(a)–(c)show the results of ISAR images processed by the compensation method with the velocities of 15,20 and 25 Mach,respectively.The problems of displacement and energy diffusion of the scattering points are effectively corrected.The compensation of displacement corrects the image distortion,and the energy focusing improves the SNR.The ISAR images of the reentry object at the three velocities are well focused.

Figure 10.The obtained ISAR images of the reentry object by traditional method.

Figure 11.ISAR images of reentry object processed by the compensation method.

In order to evaluate the improvement on ISAR imaging,the PSNR and the SSIM[27]are used to represent the improvement effect of the proposed method.The values of PSNR and the SSIM results are listed in table 5,suggesting that the quality of ISAR imaging is significantly improved.

Table 4.Results of parameter estimation by FRFT.

6.Conclusions

Table 5.Evaluation results of ISAR image using traditional method and proposed method.

With respect to plasma-sheath-enveloped reentry objects,this paper analyzes the scattering characteristics and the coupled velocities of the scattering points in the plasma sheath,from which the radar echo model of the plasma sheath is obtained.The coupled velocity will introduce LFM into pulse compression results,causing the displacement and energy diffusion,which further results in ISAR image defocus in the range dimension.By estimating the LFM parameters using FRFT,we can calculate the coupled velocity of the scattering points inversely.Based on the estimated coupled velocities,we propose a compensation method by reconstructing a modified matched filter to correct the displacement and energy diffusion in the onedimensional range profile.The simulation results have validated the effectiveness of our proposed method,laying the theoretical foundation for the ISAR imaging of plasmasheath-covered reentry objects.

Although this work analyzes the influence of the plasma sheath in the range dimension of ISAR images,the influence of the azimuth dimension is not considered.This will be the focus of our future research.

Acknowledgments

This work was supported by National Natural Science Foundation of China(No.61971330).

ORCID iDs