# Application of GO Management in Bistatic RCS Computation Using the Vector Parabolic Equation Method

Authors: X. J. Zhong; T. J. Cui; J. F. Zhang; W. M. Yu
Journal: International Journal of Antennas and Propagation
 (2013)

---

## Abstract

The parabolic equation (PE) method is a good choice in solving large-scale problems, but the resultant matrix is usually ill conditioned. In this letter, we introduce the geometric optics (GO) management in the calculation of bistatic radar cross sections using three-dimensional vector PE method. This method manages the object surface by GO, and hence the ill-conditioned problem can be avoided. Examples are given using the presented method, original method, and the method of moments. Results show the validity and stability of the presented method.


---

## Body

## 1. Introduction

The parabolic equation (PE) method is extremely powerful and had been successfully used in solving long-range forward propagation problems [1, 2]. By considering only the energy propagating along a preferred direction, the parabolic partial differential equation is used to approximate the elliptic wave equation, which can be solved by using the simple marching techniques. This method could also be used to calculate the radar cross section (RCS) of electrically large objects [3–13]. Since the paraxial direction could be chosen independently from the incident direction, the bistatic RCS of an object in a given angular sector could be computed by choosing the appropriate paraxial direction. However, the PE method usually suffers from the ill condition of the resultant matrix.In this letter, a combination of PE and the geometric optics (GO) method is proposed to compute the bistatic RCS which avoids solving the ill-conditioned equations introduced by the field components coupling at the surface points. The proposed method also makes the bistatic RCS computation more robust than the original method. Examples are given to show the advantage of the proposed method.


## 2. Vector Parabolic Equation Method

In this letter,e-iωt time dependence of fields is assumed. Working with Cartesian coordinates (x,y,z), the three-dimensional (3D) scalar wave equation for the scattering electric fields outside the objects could be given as
(1)∇2Exs+k2Exs=0,∇2Eys+k2Eys=0,∇2Ezs+k2Ezs=0.
Introduce the reduced functions associated with the field components:
(2)ums(x,y,z)=exp(-ikx)Ems(x,y,z),
where m may be x, y, or z. The definition of ums is linked to the choice of the paraxial direction x. Denoting ums by u, assuming the objects embedded in a vacuum, and substituting (2) into (1), u could be found satisfying the following equation:
(3)∂2u∂x2+2ik∂u∂x+∂2u∂y2+∂2u∂z2=0.The outgoing parabolic equation foru is then written as
(4)∂u∂x=-ik(1-Q)u
in which the pseudodifferential operator Q is defined by
(5)Q=1k2∂2∂y2+1k2∂2∂z2+1.Using the first-order Taylor expansions of the square root and the exponential, the well-known standard parabolic equation (SPE) could be obtained:(6)∂u∂x=i2k(∂2u∂y2+∂2u∂z2).
Denote the discretization sizes along the x, y, and z directions as Δx, Δy, Δz, and let Δy=Δz. To solve the partial differential equation (6), the marching technique could be used to get the solution at the range x+Δx from the fields at range x. Under the assumption of homogeneous background and letting the index of refraction n equal to 1, (6) is then discretized as
(7)ui,jm+1-ui,jmΔx=i2k(ui-1,jm+1+ui+1,jm+1-2ui,jm+1Δy2+ui,j-1m+1+ui,j+1m+1-2ui,jm+1Δz2).
The above formula could also be given as below
(8)aui-1,jm+1+bui,j-1m+1+cui,jm+1+dui,j+1m+1+eui+1,jm+1=ui,jm
in which
(9)a=b=d=e=-iΔx2kΔz2,c=1+i2ΔxΔy2k.Using the iterative solver such as BICGSTAB to solve the equations formed by (8), the fields at range x+Δx could be obtained from the fields at range x.


## 3. GO Management at Surface Points

In the vector parabolic equation (VPE) method, to consider the field components coupling on the object surfaces, it needs to solve equations at the object surface points [3]:
(10)nxEy(P)-nyEz(P)=0,nxEz(P)-nzEx(P)=0,nyEx(P)-nxEy(P)=0,∇·E→(P)s=0,
where Ex, Ey, Ez are the electric field components at the surface point P, E→(P) is the total electric field at P, and n^=x^nx+y^ny+z^nz is the outward unit normal vector at P.Now, the GO management is introduced to avoid solving the above equations. Leti^  be the incident unit vector and r^ be the reflection unit vector. Then the total incident electric field could be decomposed as
(11)EiTE(P)=E→(P)·ϕ^iTEEiTM(P)=E→(P)·ϕ^iTM
in which
(12)ϕ^iTE=n^×i^|n^×i^|,ϕ^iTM=ϕ^iTE×i^|ϕ^iTE×i^|.Let Rv and Rh be the reflection coefficients of the TE and TM fields on the object surface, respectively, and DF be the divergence factor at the point. The reflected TE and TM fields could be obtained by(13)ErTE(P)=EiTE(P)·Rv·DF,ErTM(P)=EiTM(P)·Rh·DF.Furthermore, the total reflected field could be given as(14)E→r(P)=ErTE(P)·ϕ^rTE+ErTM(P)·ϕ^rTM,
where
(15)ϕ^rTE=n^×r^|n^×r^|,ϕ^rTM=ϕ^rTE×r^|ϕ^rTE×r^|.So the major steps of using GO management are as follows: firstly, acquire the incident field on the object surface according to the incident angle; secondly, decompose the incident field using formula (11); thirdly, calculate the total reflected field using formulas (13) and (14).Once the total reflected field is obtained, the initial scattering field at each direction is determined.


## 4. RCS Calculations

The RCS calculation could be processed using the marching techniques with formula (8). At each range x where the objects are involved, the GO management should be taken after the marching process is finished so that the field components coupling effects on the object surfaces could be considered. Since the marching technique is used, the computation complexity of the proposed method is mainly in proportion to the transverse size of the objects. Hence the method could solve the scattering problems of very large objects on single personal computers.Examples are given to verify the presented method. The first example is a perfectly electrical conducting (PEC) cuboid with thickness of 0.1 m, height of 20 m, and length of 20 m. The thickness direction is alongx. The incident signal is E→=eik→·r→z^ with the incident angle of (π/2,0), k→=-x^, operating at the frequency of 300 MHz.Figure1 shows the RCS calculated by PE methods, in which PE denotes the result of PE with GO management, PE-O denotes the original PE method [3], and ϕ is the receiving angle with the unit of degree. As comparison, the full-wave method of moments is also used to calculate the problem, in which the multilevel fast multipole method (MLFMM) [14, 15] has been adopted to accelerate calculations. The CPU (Xeon E5420@2.5 GHz) times used by PE methods are both about ten seconds (the iteration number of bicgstab codes for PE is 6 and for PE-O is 8; the relative residual is 1e-6), while the CPU time for MLFMM is about 4180 seconds. Though much fast as they are, the results calculated by the PE methods have excellent agreements to those from MLFMM, especially in the paraxial direction.Figure 1
RCS of a PEC cuboid calculated by PEs and MLFMA.A PEC cylinder with the radius of 2.0 m, height of 3.0 m, and the axis along thez direction is also calculated using the same incident signal. Figure 2 gives the results from PE and MLFMM. Different from the first example, the original PE method failed to complete the computation due to bad condition numbers. From Figure 2, it is shown that the RCS calculated by PE (with iteration number of 25, the relative residual is 1e-6) coincides very well with that from MLFMM. The difference between PE and MLFMM results is little larger on the directions deviate from the paraxial direction; this is caused by the shortcoming of the PE method, and the larger angle deviated from the paraxial direction, the less accurate the results would be. This problem could be alleviated through changing the paraxial direction. We also observe that the PE method could be more accurate when the objects become larger, and an example of electrical large object is shown later.Figure 2
RCS of a PEC cylinder calculated by PE and MLFMA.To further validate the proposed method, a complicated PEC object with the cross section was shown in Figure3 and the thickness of 0.3 m is also considered under the same incident signal. Figure 4 illustrates the results from PE and MLFMM, which have also excellent agreements.Figure 3
The cross section of a PEC object.Figure 4
RCS of the PEC object shown in Figure3 calculated from PE and MLFMM.Finally, a PEC square plate with the length of 100 m is considered to test the computation capability when an electrically large object is involved, as shown in Figure5. Under the illumination of the same incident field, the computed RCS by VPE and MLFMM are illustrated in Figure 6, from which we observe that the computational results have excellent agreements. Although the object is electrically large, it takes less than two minutes to compute the results using the VPE method on a PC.Figure 5
PEC square plate with length of 100 m.Figure 6
RCS of PEC square calculated by PE and MLFMM.


## 5. Conclusion

The RCS computation is accelerated by using the VPE formulas together with the GO management at the object surfaces. Since it does not need to solve equations involved at the objects’ surfaces, the computation is very efficient and robust and could handle more complex objects than original PE method. Numerical examples show the validity and capability of the proposed method.



---

*Source: 560970-2013-07-25.xml*