# Parallel Computation of RCS of Electrically Large Platform with Coatings Modeled with NURBS Surfaces

Authors: Ying Yan; Wei-Xin Zhao; Xun-Wang Zhao; Yu Zhang; Chang-Hong Liang
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

The significance of Radar Cross Section (RCS) in the military applications makes its prediction an important problem. This paper uses large-scale parallel Physical Optics (PO) to realize the fast computation of RCS to electrically large targets, which are modeled by Non-Uniform Rational B-Spline (NURBS) surfaces and coated with dielectric materials. Some numerical examples are presented to validate this paper’s method. In addition, 1024 CPUs are used in Shanghai Supercomputer Center (SSC) to perform the simulation of a model with the maximum electrical size 1966.7λ for the first time in China. From which, it can be found that this paper’s method can greatly speed the calculation and is capable of solving the real-life problem of RCS prediction.


---

## Body

## 1. Introduction

Since the invention of radar, many applications have been found for its use in military, which include the detection and the tracking of ships, aircrafts, and ground vehicles. Modern radar systems make it virtually impossible for vehicles to approach targets without being detected at great distance. For this single reason, the military has put great emphasis on the reduction of RCS of military vehicles. The word “stealth” has been used to describe military vehicles specially designed with a low RCS.PO technique is a powerful tool to analyze the target with electrically large dimensions, especially to predict its RCS. However, many researches of PO technique focus on the problem of target with PEC surface, not only the triangular surface modeling technique [1–6] but also the NURBS modeling technique [7–10], but the study of PO theory to calculate the body with dielectric coatings is not comprehensive. As we know, the platform usually has dielectric coatings on its surface with different aims. Specially, in the military applications, coating on the surface has become a more and more significant method to reduce the vehicle’s RCS. Therefore, this paper studies the PO theory to calculate RCS of target with dielectric coatings. Moreover, there also have been some research papers that describe the PO technique to compute RCS of body with coatings [11–14], but most of them are serial significant code. When the model simulated has very large electrical dimensions, the simulation time will be very long, and the computational efficiency will be pretty low as well. To make the PO technique for model with coatings more practical, this paper also proposes a large-scale parallel scheme, which is based on MPI (Message Passing Interface) [15] library. And also a test in SSC [16] using 1024 CPUs is performed.In addition, as a popular modeling technique, NURBS [7–10] is used to achieve the object of accurate modeling.In conclusion, this paper realizes the rapid prediction of RCS to electrically large platform coated with dielectric materials by using more than 1000 CPUs. In Section2, the basic theory and computational platforms are presented; Section 3 presents some numerical examples to validate the accuracy and effectiveness of this paper’s approach; finally, Section 4 makes the conclusion.


## 2. Basic Theory

### 2.1. Initializing of the NURBS Model

As we know, NURBS is a very powerful tool for geometric design due to several reasons including fast and numerically stable algorithms and unified mathematical basis for representing complex objects [17]. Using it to model target geometry is good for target with curved and smooth surfaces so as to guarantee the accuracy of the model and make it more realistic.However, although target modeling is easy with NURBS, applying some numerical tools directly on NURBS is not so straightforward. Since the surface point function and derivative computation is simpler for Rational Bezier surfaces, NURBS surfaces are usually converted to Rational Bezier surfaces by Cox-de Boor Algorithm [18], which can be directly realized by the commercial modeling software, such as Rhinoceros. Therefore, in this paper, Rational Bezier surfaces are used for surface modeling.Before the simulation, we should firstly generate a NURBS model using the commercial software, such as Rhinoceros. In this case, NURBS model is flexible for arbitrary target. “.obj” is a file format to store the NURBS model and the file of this format is simple to handle in the computation. Therefore, this paper uses this format to store the data of input model.The initialing of the NURBS model can be summarized into two steps:(1) the NURBS surfaces composing the model should be converted to Rational Bezier surfaces, and this step can be directly realized by Rhinoceros; (2) the converted model is imported, and the code in this paper can automatically mesh each Bezier surface to a lot of small planar patches, which means that it is not necessary to use any commercial software to do the operation of mesh. In addition, the mesh size of two directions U and V, of the Bezier surface can be adjusted to acquire the best effect, which is shown in Figure 1.Illustrating of Rational Bezier surfaces with different mesh sizes.
(a)
Model of a Rational Bezier surface
(b)
Both directions ofU and V have 15 segments
(c)
DirectionU has 30 segments and V has 50 segmentsAs illustrated in Figure1, it can be found that the mesh size of each direction will affect the accuracy in modeling so that the computational precision will be influenced as well. Therefore, the choice of the mesh sizes in each direction is important.


### 2.2. RCS Calculation of Target with Dielectric Coatings

Consider the situation depicted in Figure2. An arbitrary scattered body is placed at the origin, with the observation point being at coordinates (x,y,z). For RCS calculation, the observation point is considered to be in the far field of the target. Therefore, the vectors R→ and r→ are approximately parallel.Figure 2
Far field scattering from an arbitrary body.The formula to calculate the scattered electric field is [19]
(1)E→(r→)=1jwε⋅4π∬s′[3-k2R2+j3kRR5e-jkRR→×(R→×J→s(r→′))+2J→s(r→′)1+jkRR3e-jkR]dS′,
where R=|r→-r→′|. Using the approximate formulas of far field can obtain the expression of scattered integral equation [19]
(2)|kR|≫1,R=|r→-r→′|≈r-r^⋅r→′,|r→|≫|r→′|,(3)E→s(r,θ,ϕ)=Eθθ^+Eϕϕ^=-jkZ04π·e-jkRR∬AJ→s(r→′)ejkgdS,r^·r→′=g.Generally, the magnetic field intensity of the incident wave can be obtained from the electric field intensity, thenJ→s=2n^×H→i is used to obtain the induced surface current, which can also be wrote in the form of J→s=(Jxx^+Jyy^+Jzz^)ejkh, where h=r^i·r→=xui+yvi+zwi.Usually, after the model has been meshed into a lot of small patches, the global coordinate system is not consistent with the patch’s local coordinate system. But when the global system is transformed to the local(x′′,y′′,z′′), each patch will be on its local x′′oy′′ plane and with z^′′ being the unit normal vector, and this is convenient for the computation. So we have performed this coordinate transformation and then obtained the induced current J→s in the local coordinate system:(4)J→s=(Jx′′x^+Jy′′y^)ejkh,h=r^i⋅r→=xui+yvi+zwi.In (4), the components of surface-induced current J→x′′,J→y′′ are [14, 20, 21]
(5)J→x′′=-Eiθ′′cosϕ′′Z0Γ//-Eiϕ′′cosθ′′sinϕ′′Z0Γ⊥,J→y′′=-Eiθ′′sin⁡ϕ′′Z0Γ//-Eiϕcos⁡θ′′cos⁡ϕ′′Z0Γ⊥,Γ//=-Z0cosθ′′2Rs+Z0cosθ′′,Γ⊥=-Z02Rscosθ′′+Z0,
where Rs is the surface impedance. When Rs=0, the surface is a perfect electric conductor (PEC).It is necessary to indicate that the reflection coefficientsΓ∥,Γ⊥in (5) can only be applied to the model with impedance surface. In the situation of coated model, Γ is obtained by using the theory of wave matrix [14]
(6)[c1b1]=∏n=1N1τn[ejϕnΓne-jϕnΓnejϕne-jϕn][cN+1bN+1]=[A11A12A21A22][cN+1bN+1].
In (6), b1 and c1 are the reflected wave and the incident wave in the first layer of coatings; bN,cN are the reflected wave and the incident wave in the last layer. ϕn=βntn=(2π/λn)tn, with λn being the wavelength in layer n and tn  being the thickness of layer n. Finally, the overall reflection coefficient is Γ=b1/c1=A21/A11.SubstitutingΓ into (5) can obtain the induced surface current in local coordinate system. Obviously, (6) has to be used twice to compute TE and TM situations, respectively.When the surface current is obtained, the total scattered field can be calculated by substitutingJx′′, Jy′′ into (3). And also, it is necessary to mention that the amplitude of the current is approximate constant over the surface. For the reason that this paper’s method automatically mesh the NURBS surface into some small patches, and when the mesh size is small enough, each patch will be planar:
(7)E→S(r→)=-jkZ04π⋅e-jkRR(Jx′′x^+Jy′′y^)∬s′ejk(g+h)dS′.Finally, Ludwig integral [22, 23] formula is used to calculate Ic=∬s′ejk(g+h)dS′. Substituting the integral result into (7) can get the total scattered field in local coordinate system and then transform it back to global coordinate system and RCS is obtained by the formula σ=limR→∞4πR2|E→S|2/|E→I|2.


### 2.3. Large-Scale Parallel Scheme

As we know, the model has to be meshed into a lot of elements before the computation. When the model is of electrically large dimensions, the number of meshed elements will be large. In this case, the computational quantity will be tremendous when Bistatic RCS with solid angles (both the azimuth angleϕ and the polarizing angle θ are changing simultaneously) is computed. If the serial code is still used, it will take a long time to perform the simulation, which is not practical in the engineering. Therefore, the large-scale parallel technique is employed, the scheme of which is described as follows.At a certain frequency, the parallel scheme is divided into two steps:(1) since the direction of incident wave is fixed, the sheltering procedure should be performed once, the computation of induced surface current as well. With this idea in mind, different patches are assigned to different processes periodically, and then each process conducts the sheltering, respectively. This is beneficial to obtain a good load balance. The scheme of assignment is illustrated in Figure 3, which takes 9 patches and 3 processes as an example. Patches 1, 4, and 7 are assigned to process 0, patches 2, 5, 8 to process 1, and patches 3, 6, 9 to process 2. When the sheltering procedure is finished, each process will continually calculate the induced surface current, respectively. Finally, all the results of different processes are gathered by the main process 0; (2) in scattered direction, all the azimuth angles ϕ and the polarizing angles θ are assigned to all processes periodically, and then they start to calculate the PO integral, respectively, as illustrated in Figure 4. Finally, the results of all the processes are collected by the main process 0 [19].Figure 3
Assign patches to different processes periodically.Figure 4
Assign scattered angles to different processes periodically.From what has been presented above, the large-scale parallel scheme is summarized as follows:(1)
MPI initializing;(2)
assign different patches to all the processes periodically, and then each process performs the sheltering procedure and computes the induced surface current, respectively; the main process collects all the processes’ results of sheltering and induced surface current at last;(3)
assign the scattered angles to all the processes periodically, and then each process computes the PO integral and RCS, respectively. Finally, the main process gathers all the processes’ results;(4)
MPI terminates.


### 2.4. Description of the Computational Platforms

To illustrate the versatility of the solver in this paper, two representative computation platforms have been chosen.(1)
Personal Computer: Quad core Intel i5 processor (2.67 GHz) with 4 GB RAM and 500 GB of hard disk.(2)
Shanghai Supercomputer Center (SSC) [16]: The 37 nodes from Magic-cube Machine with a total of 592 AMD CPU cores (1.9 GHz per CPU and 4 cores on each CPU): 16 CPU cores on each node and 4 GB RAM per core, and a total amount of RAM approximately equal to 2.3 TB. No hard disk storage is available for computation. InfiniBand is used for the network’s interconnection.


## 2.1. Initializing of the NURBS Model

As we know, NURBS is a very powerful tool for geometric design due to several reasons including fast and numerically stable algorithms and unified mathematical basis for representing complex objects [17]. Using it to model target geometry is good for target with curved and smooth surfaces so as to guarantee the accuracy of the model and make it more realistic.However, although target modeling is easy with NURBS, applying some numerical tools directly on NURBS is not so straightforward. Since the surface point function and derivative computation is simpler for Rational Bezier surfaces, NURBS surfaces are usually converted to Rational Bezier surfaces by Cox-de Boor Algorithm [18], which can be directly realized by the commercial modeling software, such as Rhinoceros. Therefore, in this paper, Rational Bezier surfaces are used for surface modeling.Before the simulation, we should firstly generate a NURBS model using the commercial software, such as Rhinoceros. In this case, NURBS model is flexible for arbitrary target. “.obj” is a file format to store the NURBS model and the file of this format is simple to handle in the computation. Therefore, this paper uses this format to store the data of input model.The initialing of the NURBS model can be summarized into two steps:(1) the NURBS surfaces composing the model should be converted to Rational Bezier surfaces, and this step can be directly realized by Rhinoceros; (2) the converted model is imported, and the code in this paper can automatically mesh each Bezier surface to a lot of small planar patches, which means that it is not necessary to use any commercial software to do the operation of mesh. In addition, the mesh size of two directions U and V, of the Bezier surface can be adjusted to acquire the best effect, which is shown in Figure 1.Illustrating of Rational Bezier surfaces with different mesh sizes.
(a)
Model of a Rational Bezier surface
(b)
Both directions ofU and V have 15 segments
(c)
DirectionU has 30 segments and V has 50 segmentsAs illustrated in Figure1, it can be found that the mesh size of each direction will affect the accuracy in modeling so that the computational precision will be influenced as well. Therefore, the choice of the mesh sizes in each direction is important.


## 2.2. RCS Calculation of Target with Dielectric Coatings

Consider the situation depicted in Figure2. An arbitrary scattered body is placed at the origin, with the observation point being at coordinates (x,y,z). For RCS calculation, the observation point is considered to be in the far field of the target. Therefore, the vectors R→ and r→ are approximately parallel.Figure 2
Far field scattering from an arbitrary body.The formula to calculate the scattered electric field is [19]
(1)E→(r→)=1jwε⋅4π∬s′[3-k2R2+j3kRR5e-jkRR→×(R→×J→s(r→′))+2J→s(r→′)1+jkRR3e-jkR]dS′,
where R=|r→-r→′|. Using the approximate formulas of far field can obtain the expression of scattered integral equation [19]
(2)|kR|≫1,R=|r→-r→′|≈r-r^⋅r→′,|r→|≫|r→′|,(3)E→s(r,θ,ϕ)=Eθθ^+Eϕϕ^=-jkZ04π·e-jkRR∬AJ→s(r→′)ejkgdS,r^·r→′=g.Generally, the magnetic field intensity of the incident wave can be obtained from the electric field intensity, thenJ→s=2n^×H→i is used to obtain the induced surface current, which can also be wrote in the form of J→s=(Jxx^+Jyy^+Jzz^)ejkh, where h=r^i·r→=xui+yvi+zwi.Usually, after the model has been meshed into a lot of small patches, the global coordinate system is not consistent with the patch’s local coordinate system. But when the global system is transformed to the local(x′′,y′′,z′′), each patch will be on its local x′′oy′′ plane and with z^′′ being the unit normal vector, and this is convenient for the computation. So we have performed this coordinate transformation and then obtained the induced current J→s in the local coordinate system:(4)J→s=(Jx′′x^+Jy′′y^)ejkh,h=r^i⋅r→=xui+yvi+zwi.In (4), the components of surface-induced current J→x′′,J→y′′ are [14, 20, 21]
(5)J→x′′=-Eiθ′′cosϕ′′Z0Γ//-Eiϕ′′cosθ′′sinϕ′′Z0Γ⊥,J→y′′=-Eiθ′′sin⁡ϕ′′Z0Γ//-Eiϕcos⁡θ′′cos⁡ϕ′′Z0Γ⊥,Γ//=-Z0cosθ′′2Rs+Z0cosθ′′,Γ⊥=-Z02Rscosθ′′+Z0,
where Rs is the surface impedance. When Rs=0, the surface is a perfect electric conductor (PEC).It is necessary to indicate that the reflection coefficientsΓ∥,Γ⊥in (5) can only be applied to the model with impedance surface. In the situation of coated model, Γ is obtained by using the theory of wave matrix [14]
(6)[c1b1]=∏n=1N1τn[ejϕnΓne-jϕnΓnejϕne-jϕn][cN+1bN+1]=[A11A12A21A22][cN+1bN+1].
In (6), b1 and c1 are the reflected wave and the incident wave in the first layer of coatings; bN,cN are the reflected wave and the incident wave in the last layer. ϕn=βntn=(2π/λn)tn, with λn being the wavelength in layer n and tn  being the thickness of layer n. Finally, the overall reflection coefficient is Γ=b1/c1=A21/A11.SubstitutingΓ into (5) can obtain the induced surface current in local coordinate system. Obviously, (6) has to be used twice to compute TE and TM situations, respectively.When the surface current is obtained, the total scattered field can be calculated by substitutingJx′′, Jy′′ into (3). And also, it is necessary to mention that the amplitude of the current is approximate constant over the surface. For the reason that this paper’s method automatically mesh the NURBS surface into some small patches, and when the mesh size is small enough, each patch will be planar:
(7)E→S(r→)=-jkZ04π⋅e-jkRR(Jx′′x^+Jy′′y^)∬s′ejk(g+h)dS′.Finally, Ludwig integral [22, 23] formula is used to calculate Ic=∬s′ejk(g+h)dS′. Substituting the integral result into (7) can get the total scattered field in local coordinate system and then transform it back to global coordinate system and RCS is obtained by the formula σ=limR→∞4πR2|E→S|2/|E→I|2.


## 2.3. Large-Scale Parallel Scheme

As we know, the model has to be meshed into a lot of elements before the computation. When the model is of electrically large dimensions, the number of meshed elements will be large. In this case, the computational quantity will be tremendous when Bistatic RCS with solid angles (both the azimuth angleϕ and the polarizing angle θ are changing simultaneously) is computed. If the serial code is still used, it will take a long time to perform the simulation, which is not practical in the engineering. Therefore, the large-scale parallel technique is employed, the scheme of which is described as follows.At a certain frequency, the parallel scheme is divided into two steps:(1) since the direction of incident wave is fixed, the sheltering procedure should be performed once, the computation of induced surface current as well. With this idea in mind, different patches are assigned to different processes periodically, and then each process conducts the sheltering, respectively. This is beneficial to obtain a good load balance. The scheme of assignment is illustrated in Figure 3, which takes 9 patches and 3 processes as an example. Patches 1, 4, and 7 are assigned to process 0, patches 2, 5, 8 to process 1, and patches 3, 6, 9 to process 2. When the sheltering procedure is finished, each process will continually calculate the induced surface current, respectively. Finally, all the results of different processes are gathered by the main process 0; (2) in scattered direction, all the azimuth angles ϕ and the polarizing angles θ are assigned to all processes periodically, and then they start to calculate the PO integral, respectively, as illustrated in Figure 4. Finally, the results of all the processes are collected by the main process 0 [19].Figure 3
Assign patches to different processes periodically.Figure 4
Assign scattered angles to different processes periodically.From what has been presented above, the large-scale parallel scheme is summarized as follows:(1)
MPI initializing;(2)
assign different patches to all the processes periodically, and then each process performs the sheltering procedure and computes the induced surface current, respectively; the main process collects all the processes’ results of sheltering and induced surface current at last;(3)
assign the scattered angles to all the processes periodically, and then each process computes the PO integral and RCS, respectively. Finally, the main process gathers all the processes’ results;(4)
MPI terminates.


## 2.4. Description of the Computational Platforms

To illustrate the versatility of the solver in this paper, two representative computation platforms have been chosen.(1)
Personal Computer: Quad core Intel i5 processor (2.67 GHz) with 4 GB RAM and 500 GB of hard disk.(2)
Shanghai Supercomputer Center (SSC) [16]: The 37 nodes from Magic-cube Machine with a total of 592 AMD CPU cores (1.9 GHz per CPU and 4 cores on each CPU): 16 CPU cores on each node and 4 GB RAM per core, and a total amount of RAM approximately equal to 2.3 TB. No hard disk storage is available for computation. InfiniBand is used for the network’s interconnection.


## 3. Results and Discussion

### 3.1. Results Validation

In this section, the accuracy of this paper’s method is validated through three different examples: Section3.1.1 a missile’s head whose surfaces are coated with dielectric materials; Section 3.1.2 an aircraft whose surfaces are coated with dielectric materials; Section 3.1.3 a Rational Bezier surface coated with dielectric materials.


#### 3.1.1. The Model of a Missile’s Head

As illustrated in Figure5, the model simulated is the head of a missile with dielectric materials coated on its surfaces. The model consists of 4 Rational Bezier surfaces and the geometrical dimensions are 8.1 m × 10.1 m × 10.1 m. Its surfaces are coated with a layer of dielectric material, whose parameters are presented in Table 1. The simulation frequency is 300 MHz.Table 1
Information of the material coating.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
Thicknesst (mm)
16.0
0.275
1.0
5Figure 5
RationalBezier surfaces separated from a missile model.In this example, Bistatic RCS is computed, and the results obtained are compared with those from the commercial software FEKO [24] to validate the accuracy of this paper’s method. The compared algorithms are chosen to be PO and RWG MoM [25] in FEKO.The incident wave is of VV (θθ) polarization and along the direction of θ=0∘,ϕ=0∘ and the direction of the scattered wave is θ=0∘~360∘,ϕ=0∘.Figure6 presents the results’ comparison. In Figure 6(a), PO algorithm in FEKO is chosen to be the compared technique to validate this paper’s method. It can found from the figure that this paper’s results are in good agreement with those of FEKO, except some tiny differences near θ=110∘,250∘; in Figure 6(b), we choose a more exact solution—RWG MoM to compare with this paper’s method. From the figure, we can find that the results’ comparison among polarizing angles θ=20°~120° and θ=230°~300° exists some differences; this is reasonable because of the approximation introduced in the PO current. In addition, the results agree well among other angle ranges.Comparison of this paper’s results and the results of FEKO.
(a)
This paper’s method versus PO algorithm in FEKO
(b)
This paper’s method versus RWG MoM algorithm in FEKO


#### 3.1.2. The Model of an Aircraft

In this example, Bistatic RCS of an aircraft coated with dielectric materials is computed. As illustrated in Figure7, the model’s geometrical dimensions are 11.6 m × 7.0 m × 2.92 m, and it is necessary to mention that the model has no structures such as engine and air inlets (i.e., it is not a real aircraft).Figure 7
Model of an aircraft.(
1
) Accuracy Validation
In this part, the model is coated with one layer of coating, which is the first kind of material presented in Table2.
The simulation frequency is 1 GHz. The incident wave is of VV (θθ) polarization and along the direction of the aircraft’s nose. The direction of the scattered wave is ϕ=0∘~360∘,θ=90∘. This simulation is performed on the first kind of computational platform described in Section 2.4, and the results obtained are compared with those of FEKO. Since this model’s electrical dimensions are big and it will take long time to perform the simulation by using RWG MoM, PO algorithm in FEKO is chosen to be the validation method.
Figure8 illustrates the comparison of Bistatic RCS on XOY plane, which indicates that the results of this paper are in good agreement with the results of FEKO, which shows that this paper’s method is correct.
Moreover, Table3 compares the memory required in the simulation, which indicates that this paper’s method needs less memory than that of FEKO.Table 2
Information about the coating.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
Thicknesst (mm)
1st material
2.5
0.00045
1.0
5
2nd material
16
0.275
1.0
6.25Table 3
Comparison of memory required.
Memory required (MB)
This paper’s method
5.3
PO algorithm in FEKO
15.3Figure 8
This paper’s results versus the results of FEKO.(
2
) RCS Reduction by Using Coatings
In this part, two situations are simulated to compare the variation of RCS:(1) PEC model and (2) model coated with two layers of dielectric materials.
The model simulated is the same as that in part 1 of this example. Two layers of coating, which are the two kinds of materials presented in Table2, are coated on the aircraft. The frequencies range from 0.1 GHz to 10 GHz. The polarization of the incident wave is VV (θθ), and the direction of backward scattering is along the nose of the aircraft, namely, θ=90∘,ϕ=0∘. This simulation is performed on the first kind of computational platform described in Section 2.4.
Figure9 indicates that two layers of coating will significantly reduce the model’s RCS. To be specific, there is a reduction of about 10 dB at approximate 3 GHz, and at the range of 7 Hz to 10 Hz, RCS is reduced a lot. This example shows the great influence of coating on reducing the target’s RCS.Figure 9
Aircraft with PEC surface versus Aircraft coated with two layers of materials.


#### 3.1.3. The Model of a Rational Bezier Surface

In this example, a Rational Bezier surface’s Monostatic RCS is computed. The results obtained are compared with those from Higher-order MoM [26] in HOBBIES [27, 28] to specify in what situation this paper’s method is valid.The model’s geometrical dimensions are 2 m× 1.1 m × 2 m. It is coated with one layer of material and the parameters of which are presented in Table 1. In Figures 10(a) and 10(b), the simulation frequencies are 1 GHz and 5 GHz, respectively. The incident wave is of the VV (θθ) polarization, and the direction of the backward scattering is ϕ=270°, θ=0°~180°. Figure 10(c) is to verify the accuracy of this paper’s method at various frequencies which range from 300 MHz to 5 GHz, and the backward scattering is along ϕ=270°, θ=90°.This paper’s method versus Higher-order MoM in HOBBIES.
(a)
1 GHz
(b)
5 GHz
(c)
300 MHz~5 GHzFrom the results illustrated in Figure10, we can find that both the results of 1 GHz and 5 GHz do not agree well at the angle ranges of θ=0°~30° and θ=150°~180°, which is because there is a structure of edge in the model among these angles, and the diffraction effect that contributes to the scattered field is not considered so that the results near the edge will not be accurate. Moreover, it can also be found that the results of 5 GHz agree better than that of 1 GHz at the angle range of θ=30∘~150∘, which indicates that the results will be more accurate when the frequency is higher. Meanwhile, we can also find the better agreement at high frequencies from Figure 10(c).From what has been described above, this paper’s method is accurate when the incident direction is not near the structure of edge of the model and the simulated frequency is high.


### 3.2. Parallel Speed-Up Ratio

In this example, a formation of 7 aircrafts’ Bistatic RCS is calculated by using the large-scale parallel technique, which is based on MPI. The model is presented in Figure11. Each aircraft has the geometrical dimensions of 19.2 m × 12.2 m × 4.3 m and contains the structure of engine. The model is coated with one layer of material, and the parameters of which are presented in Table 1.Figure 11
A formation of 7 aircrafts.As the simulation frequency is 5.0 GHz, the model is automatically meshed into 1094870 elements by this solver. The incident wave of VV (θθ) polarization is along the direction of θ=90∘, ϕ=270∘, and the scattered direction is θ=0∘~180∘, ϕ=0∘~360∘. The second kind of computational platform described in Section 2.4 is used to perform this test.Table4 indicates the time comparison of using different numbers of CPU: it takes 164052 seconds by using 1 CPU, namely, it needs more than 45 hours to perform the simulation. However, only 2670 seconds is required by using 64 CPUs, which means that the computational efficiency is greatly improved. Meanwhile, the illustrating of parallel speed-up ratio in Figure 12 is linear. Namely, with more CPUs, the computation time will be reduced linearly and also the parallel speed-up ratio will be increased linearly. This is consistent with the theory that parallel PO has a linear parallel speed-up ratio [17].Table 4
Time comparison of using different numbers of CPU.
Number of CPU
1
4
8
16
32
64
Wall Clock Time (Second)
164052
41054
20459
10475
5293
2670Figure 12
Parallel speed-up ratio.Finally, the Bistatic RCS are also illustrated in Figure13.Bistatic RCS.
(a)
XOY plane
(b)
YOZ plane


### 3.3. Large-Scale Parallel Computation by Using More Than 1000 CPUs

In this example, an electrically large model of a ship on the sea is simulated by using 1024 CPUs. The whole model consists of two parts: a ship and the sea, which are illustrated in Figure14. The ship is a PEC model and the sea is of dielectric material whose parameters are presented in Table 5. The geometrical dimensions of the ship are 30 m × 2.2 m × 5.2 m, and those of the sea are 59 m × 59 m. As the simulation frequency is 10 GHz, the ship and the sea have the electrical dimensions of 1000 λ× 73.3 λ× 173.3 λ and 1966.7 λ× 1966.7 λ, respectively. In addition, the elements of the ship generated are 740708, and those of the sea are 8809388. Therefore, the whole model (ship + sea) has the total elements of 9550096.Table 5
Dielectric material of the sea.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
81.5
0.0
1.0Figure 14
Model of a ship on the sea.Monostatic RCS is calculated. The incident wave is of HH (ϕϕ) polarization, and the backward scattering direction is θ=0∘~90∘, ϕ=0∘~360∘ (i.e., the space above the sea).From the description above, we can find that it is unrealistic to solve this problem by RWG MoM algorithm because of the large quantity of elements. Even if PO technique is used, the serial code will also encounter so large computational quantity that requires very long time for simulation. Thus, we employ large-scale parallel technique described in Section2.3 and use 1024 CPUs in SSC to solve this problem. 3D and 2D results are presented in Figures 15 and 16, respectively. Finally, the computation information is listed in Table 6.Table 6
Information about the simulation.
Project’s name
Number of CPU
Wall Clock Time (Second)
A ship on the sea
1024
1366Figure 15
3D Monostatic RCS (unit : dBsm).2D Monostatic RCS.
(a)
θ
=
9
0
∘, ϕ=0∘~360∘
(b)
θ
=
-
9
0
∘
~
9
0
∘, ϕ=0∘
(c)
θ
=
-
9
0
∘
~
9
0
∘, ϕ=90∘From the results presented above, we can find that the simulation time is significantly reduced by using 1024 CPUs; it only takes less than half an hour to solve this problem. In addition, it is necessary to indicate that this is the first time in China employing more than 1000 CPUs simultaneously to solve such a large problem by using large-scale parallel PO technique.


## 3.1. Results Validation

In this section, the accuracy of this paper’s method is validated through three different examples: Section3.1.1 a missile’s head whose surfaces are coated with dielectric materials; Section 3.1.2 an aircraft whose surfaces are coated with dielectric materials; Section 3.1.3 a Rational Bezier surface coated with dielectric materials.


### 3.1.1. The Model of a Missile’s Head

As illustrated in Figure5, the model simulated is the head of a missile with dielectric materials coated on its surfaces. The model consists of 4 Rational Bezier surfaces and the geometrical dimensions are 8.1 m × 10.1 m × 10.1 m. Its surfaces are coated with a layer of dielectric material, whose parameters are presented in Table 1. The simulation frequency is 300 MHz.Table 1
Information of the material coating.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
Thicknesst (mm)
16.0
0.275
1.0
5Figure 5
RationalBezier surfaces separated from a missile model.In this example, Bistatic RCS is computed, and the results obtained are compared with those from the commercial software FEKO [24] to validate the accuracy of this paper’s method. The compared algorithms are chosen to be PO and RWG MoM [25] in FEKO.The incident wave is of VV (θθ) polarization and along the direction of θ=0∘,ϕ=0∘ and the direction of the scattered wave is θ=0∘~360∘,ϕ=0∘.Figure6 presents the results’ comparison. In Figure 6(a), PO algorithm in FEKO is chosen to be the compared technique to validate this paper’s method. It can found from the figure that this paper’s results are in good agreement with those of FEKO, except some tiny differences near θ=110∘,250∘; in Figure 6(b), we choose a more exact solution—RWG MoM to compare with this paper’s method. From the figure, we can find that the results’ comparison among polarizing angles θ=20°~120° and θ=230°~300° exists some differences; this is reasonable because of the approximation introduced in the PO current. In addition, the results agree well among other angle ranges.Comparison of this paper’s results and the results of FEKO.
(a)
This paper’s method versus PO algorithm in FEKO
(b)
This paper’s method versus RWG MoM algorithm in FEKO


### 3.1.2. The Model of an Aircraft

In this example, Bistatic RCS of an aircraft coated with dielectric materials is computed. As illustrated in Figure7, the model’s geometrical dimensions are 11.6 m × 7.0 m × 2.92 m, and it is necessary to mention that the model has no structures such as engine and air inlets (i.e., it is not a real aircraft).Figure 7
Model of an aircraft.(
1
) Accuracy Validation
In this part, the model is coated with one layer of coating, which is the first kind of material presented in Table2.
The simulation frequency is 1 GHz. The incident wave is of VV (θθ) polarization and along the direction of the aircraft’s nose. The direction of the scattered wave is ϕ=0∘~360∘,θ=90∘. This simulation is performed on the first kind of computational platform described in Section 2.4, and the results obtained are compared with those of FEKO. Since this model’s electrical dimensions are big and it will take long time to perform the simulation by using RWG MoM, PO algorithm in FEKO is chosen to be the validation method.
Figure8 illustrates the comparison of Bistatic RCS on XOY plane, which indicates that the results of this paper are in good agreement with the results of FEKO, which shows that this paper’s method is correct.
Moreover, Table3 compares the memory required in the simulation, which indicates that this paper’s method needs less memory than that of FEKO.Table 2
Information about the coating.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
Thicknesst (mm)
1st material
2.5
0.00045
1.0
5
2nd material
16
0.275
1.0
6.25Table 3
Comparison of memory required.
Memory required (MB)
This paper’s method
5.3
PO algorithm in FEKO
15.3Figure 8
This paper’s results versus the results of FEKO.(
2
) RCS Reduction by Using Coatings
In this part, two situations are simulated to compare the variation of RCS:(1) PEC model and (2) model coated with two layers of dielectric materials.
The model simulated is the same as that in part 1 of this example. Two layers of coating, which are the two kinds of materials presented in Table2, are coated on the aircraft. The frequencies range from 0.1 GHz to 10 GHz. The polarization of the incident wave is VV (θθ), and the direction of backward scattering is along the nose of the aircraft, namely, θ=90∘,ϕ=0∘. This simulation is performed on the first kind of computational platform described in Section 2.4.
Figure9 indicates that two layers of coating will significantly reduce the model’s RCS. To be specific, there is a reduction of about 10 dB at approximate 3 GHz, and at the range of 7 Hz to 10 Hz, RCS is reduced a lot. This example shows the great influence of coating on reducing the target’s RCS.Figure 9
Aircraft with PEC surface versus Aircraft coated with two layers of materials.


### 3.1.3. The Model of a Rational Bezier Surface

In this example, a Rational Bezier surface’s Monostatic RCS is computed. The results obtained are compared with those from Higher-order MoM [26] in HOBBIES [27, 28] to specify in what situation this paper’s method is valid.The model’s geometrical dimensions are 2 m× 1.1 m × 2 m. It is coated with one layer of material and the parameters of which are presented in Table 1. In Figures 10(a) and 10(b), the simulation frequencies are 1 GHz and 5 GHz, respectively. The incident wave is of the VV (θθ) polarization, and the direction of the backward scattering is ϕ=270°, θ=0°~180°. Figure 10(c) is to verify the accuracy of this paper’s method at various frequencies which range from 300 MHz to 5 GHz, and the backward scattering is along ϕ=270°, θ=90°.This paper’s method versus Higher-order MoM in HOBBIES.
(a)
1 GHz
(b)
5 GHz
(c)
300 MHz~5 GHzFrom the results illustrated in Figure10, we can find that both the results of 1 GHz and 5 GHz do not agree well at the angle ranges of θ=0°~30° and θ=150°~180°, which is because there is a structure of edge in the model among these angles, and the diffraction effect that contributes to the scattered field is not considered so that the results near the edge will not be accurate. Moreover, it can also be found that the results of 5 GHz agree better than that of 1 GHz at the angle range of θ=30∘~150∘, which indicates that the results will be more accurate when the frequency is higher. Meanwhile, we can also find the better agreement at high frequencies from Figure 10(c).From what has been described above, this paper’s method is accurate when the incident direction is not near the structure of edge of the model and the simulated frequency is high.


## 3.1.1. The Model of a Missile’s Head

As illustrated in Figure5, the model simulated is the head of a missile with dielectric materials coated on its surfaces. The model consists of 4 Rational Bezier surfaces and the geometrical dimensions are 8.1 m × 10.1 m × 10.1 m. Its surfaces are coated with a layer of dielectric material, whose parameters are presented in Table 1. The simulation frequency is 300 MHz.Table 1
Information of the material coating.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
Thicknesst (mm)
16.0
0.275
1.0
5Figure 5
RationalBezier surfaces separated from a missile model.In this example, Bistatic RCS is computed, and the results obtained are compared with those from the commercial software FEKO [24] to validate the accuracy of this paper’s method. The compared algorithms are chosen to be PO and RWG MoM [25] in FEKO.The incident wave is of VV (θθ) polarization and along the direction of θ=0∘,ϕ=0∘ and the direction of the scattered wave is θ=0∘~360∘,ϕ=0∘.Figure6 presents the results’ comparison. In Figure 6(a), PO algorithm in FEKO is chosen to be the compared technique to validate this paper’s method. It can found from the figure that this paper’s results are in good agreement with those of FEKO, except some tiny differences near θ=110∘,250∘; in Figure 6(b), we choose a more exact solution—RWG MoM to compare with this paper’s method. From the figure, we can find that the results’ comparison among polarizing angles θ=20°~120° and θ=230°~300° exists some differences; this is reasonable because of the approximation introduced in the PO current. In addition, the results agree well among other angle ranges.Comparison of this paper’s results and the results of FEKO.
(a)
This paper’s method versus PO algorithm in FEKO
(b)
This paper’s method versus RWG MoM algorithm in FEKO


## 3.1.2. The Model of an Aircraft

In this example, Bistatic RCS of an aircraft coated with dielectric materials is computed. As illustrated in Figure7, the model’s geometrical dimensions are 11.6 m × 7.0 m × 2.92 m, and it is necessary to mention that the model has no structures such as engine and air inlets (i.e., it is not a real aircraft).Figure 7
Model of an aircraft.(
1
) Accuracy Validation
In this part, the model is coated with one layer of coating, which is the first kind of material presented in Table2.
The simulation frequency is 1 GHz. The incident wave is of VV (θθ) polarization and along the direction of the aircraft’s nose. The direction of the scattered wave is ϕ=0∘~360∘,θ=90∘. This simulation is performed on the first kind of computational platform described in Section 2.4, and the results obtained are compared with those of FEKO. Since this model’s electrical dimensions are big and it will take long time to perform the simulation by using RWG MoM, PO algorithm in FEKO is chosen to be the validation method.
Figure8 illustrates the comparison of Bistatic RCS on XOY plane, which indicates that the results of this paper are in good agreement with the results of FEKO, which shows that this paper’s method is correct.
Moreover, Table3 compares the memory required in the simulation, which indicates that this paper’s method needs less memory than that of FEKO.Table 2
Information about the coating.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
Thicknesst (mm)
1st material
2.5
0.00045
1.0
5
2nd material
16
0.275
1.0
6.25Table 3
Comparison of memory required.
Memory required (MB)
This paper’s method
5.3
PO algorithm in FEKO
15.3Figure 8
This paper’s results versus the results of FEKO.(
2
) RCS Reduction by Using Coatings
In this part, two situations are simulated to compare the variation of RCS:(1) PEC model and (2) model coated with two layers of dielectric materials.
The model simulated is the same as that in part 1 of this example. Two layers of coating, which are the two kinds of materials presented in Table2, are coated on the aircraft. The frequencies range from 0.1 GHz to 10 GHz. The polarization of the incident wave is VV (θθ), and the direction of backward scattering is along the nose of the aircraft, namely, θ=90∘,ϕ=0∘. This simulation is performed on the first kind of computational platform described in Section 2.4.
Figure9 indicates that two layers of coating will significantly reduce the model’s RCS. To be specific, there is a reduction of about 10 dB at approximate 3 GHz, and at the range of 7 Hz to 10 Hz, RCS is reduced a lot. This example shows the great influence of coating on reducing the target’s RCS.Figure 9
Aircraft with PEC surface versus Aircraft coated with two layers of materials.


## 3.1.3. The Model of a Rational Bezier Surface

In this example, a Rational Bezier surface’s Monostatic RCS is computed. The results obtained are compared with those from Higher-order MoM [26] in HOBBIES [27, 28] to specify in what situation this paper’s method is valid.The model’s geometrical dimensions are 2 m× 1.1 m × 2 m. It is coated with one layer of material and the parameters of which are presented in Table 1. In Figures 10(a) and 10(b), the simulation frequencies are 1 GHz and 5 GHz, respectively. The incident wave is of the VV (θθ) polarization, and the direction of the backward scattering is ϕ=270°, θ=0°~180°. Figure 10(c) is to verify the accuracy of this paper’s method at various frequencies which range from 300 MHz to 5 GHz, and the backward scattering is along ϕ=270°, θ=90°.This paper’s method versus Higher-order MoM in HOBBIES.
(a)
1 GHz
(b)
5 GHz
(c)
300 MHz~5 GHzFrom the results illustrated in Figure10, we can find that both the results of 1 GHz and 5 GHz do not agree well at the angle ranges of θ=0°~30° and θ=150°~180°, which is because there is a structure of edge in the model among these angles, and the diffraction effect that contributes to the scattered field is not considered so that the results near the edge will not be accurate. Moreover, it can also be found that the results of 5 GHz agree better than that of 1 GHz at the angle range of θ=30∘~150∘, which indicates that the results will be more accurate when the frequency is higher. Meanwhile, we can also find the better agreement at high frequencies from Figure 10(c).From what has been described above, this paper’s method is accurate when the incident direction is not near the structure of edge of the model and the simulated frequency is high.


## 3.2. Parallel Speed-Up Ratio

In this example, a formation of 7 aircrafts’ Bistatic RCS is calculated by using the large-scale parallel technique, which is based on MPI. The model is presented in Figure11. Each aircraft has the geometrical dimensions of 19.2 m × 12.2 m × 4.3 m and contains the structure of engine. The model is coated with one layer of material, and the parameters of which are presented in Table 1.Figure 11
A formation of 7 aircrafts.As the simulation frequency is 5.0 GHz, the model is automatically meshed into 1094870 elements by this solver. The incident wave of VV (θθ) polarization is along the direction of θ=90∘, ϕ=270∘, and the scattered direction is θ=0∘~180∘, ϕ=0∘~360∘. The second kind of computational platform described in Section 2.4 is used to perform this test.Table4 indicates the time comparison of using different numbers of CPU: it takes 164052 seconds by using 1 CPU, namely, it needs more than 45 hours to perform the simulation. However, only 2670 seconds is required by using 64 CPUs, which means that the computational efficiency is greatly improved. Meanwhile, the illustrating of parallel speed-up ratio in Figure 12 is linear. Namely, with more CPUs, the computation time will be reduced linearly and also the parallel speed-up ratio will be increased linearly. This is consistent with the theory that parallel PO has a linear parallel speed-up ratio [17].Table 4
Time comparison of using different numbers of CPU.
Number of CPU
1
4
8
16
32
64
Wall Clock Time (Second)
164052
41054
20459
10475
5293
2670Figure 12
Parallel speed-up ratio.Finally, the Bistatic RCS are also illustrated in Figure13.Bistatic RCS.
(a)
XOY plane
(b)
YOZ plane


## 3.3. Large-Scale Parallel Computation by Using More Than 1000 CPUs

In this example, an electrically large model of a ship on the sea is simulated by using 1024 CPUs. The whole model consists of two parts: a ship and the sea, which are illustrated in Figure14. The ship is a PEC model and the sea is of dielectric material whose parameters are presented in Table 5. The geometrical dimensions of the ship are 30 m × 2.2 m × 5.2 m, and those of the sea are 59 m × 59 m. As the simulation frequency is 10 GHz, the ship and the sea have the electrical dimensions of 1000 λ× 73.3 λ× 173.3 λ and 1966.7 λ× 1966.7 λ, respectively. In addition, the elements of the ship generated are 740708, and those of the sea are 8809388. Therefore, the whole model (ship + sea) has the total elements of 9550096.Table 5
Dielectric material of the sea.
Relative Electric Permittivityεr
Loss Tangenttanδ
Relative magnetic Permeabilityμr
81.5
0.0
1.0Figure 14
Model of a ship on the sea.Monostatic RCS is calculated. The incident wave is of HH (ϕϕ) polarization, and the backward scattering direction is θ=0∘~90∘, ϕ=0∘~360∘ (i.e., the space above the sea).From the description above, we can find that it is unrealistic to solve this problem by RWG MoM algorithm because of the large quantity of elements. Even if PO technique is used, the serial code will also encounter so large computational quantity that requires very long time for simulation. Thus, we employ large-scale parallel technique described in Section2.3 and use 1024 CPUs in SSC to solve this problem. 3D and 2D results are presented in Figures 15 and 16, respectively. Finally, the computation information is listed in Table 6.Table 6
Information about the simulation.
Project’s name
Number of CPU
Wall Clock Time (Second)
A ship on the sea
1024
1366Figure 15
3D Monostatic RCS (unit : dBsm).2D Monostatic RCS.
(a)
θ
=
9
0
∘, ϕ=0∘~360∘
(b)
θ
=
-
9
0
∘
~
9
0
∘, ϕ=0∘
(c)
θ
=
-
9
0
∘
~
9
0
∘, ϕ=90∘From the results presented above, we can find that the simulation time is significantly reduced by using 1024 CPUs; it only takes less than half an hour to solve this problem. In addition, it is necessary to indicate that this is the first time in China employing more than 1000 CPUs simultaneously to solve such a large problem by using large-scale parallel PO technique.


## 4. Conclusion

In this paper, a parallel method to predict RCS of electrically large NURBS model with dielectric coatings is presented. A real-life model with extra big electrical dimensions is simulated by using 1024 CPUs, which demonstrates that this paper’s method is capable of solving the practical challenging problem.



---

*Source: 348052-2012-06-25.xml*