Abstract

Hybrid scheme combining shooting and bouncing ray with semi-deterministic facet model is proposed to analyze composite scattering from ship-ocean scene in this study. This model can deal with complex electromagnetic interaction between ship and sea surface. Thus, scattering properties of composite ship-ocean scenes with influence of various parameters (such as incident angle and wind speed) can be studied and analyzed efficiently. Studying such properties is of significance for target detection and high-resolution radar imaging in sea environments. Accuracy and performance of this method are validated and evaluated by comparing with multilevel fast multipole method of FEKO for electrically small objects. All simulation results indicate that the proposed method is suitable for providing preliminary radar cross section prediction of electrically large composite model.

1. Introduction

Electromagnetic scattering of electrically large targets above rough sea surface is common military scene in target detection and tracking [13]. In recent years, extensive endeavors have been devoted to simulation of such complex composite scattering. Method of Moment (MoM) [4] is widely used numerical algorithm, which is adequately developed into fast algorithms, such as Generalized Forward-Backward Method (GFBM) [5, 6] and multilevel fast multipole algorithm (MLFMA) [7, 8]. However, such methods require considerable computational time when dealing with electrically large targets. Asymptotic methods such as Iterative Physical Optics (IPO) which contains a large number of iterations to ensure the accuracy also face poor computing efficiency. Traditional four-path model [9] uses complex reflection coefficient to evaluate coupling interaction of ship-ocean. However, actual sea surface is usually assumed to be two-scale rough surface model [10]: capillary waves modulated by gravity wave in large scale. Thus, scattering field is not only concentrated on specular direction. Besides, four-path model cannot consider mutual shade between ship and ocean. Shooting and Bouncing Ray (SBR) [11, 12] seems more flexible and suitable for such issue. With ray tracing capacity of SBR method, multiple interactions between ship and ocean can be determined. Nevertheless, SBR method requires more refined meshes involving capillary waves on the sea surface. This requirement enormously increases computing burden. For large-scope sea surfaces, semi-deterministic facet scattering model (SDFM) [13, 14] was developed in relevant publications. In this model, elementary radar echoes from respective facets are computed by a semideterministic scheme considering both specular and diffuse configurations.

This paper aims to establish feasible scheme that conforms to sea surface two-scale model by combining SDFM and SBR method. With SDFM considering Bragg scattering from capillary waves, sea surface can be meshed into larger triangle patches, which adapt well to SBR method for ray tracing between the ship and ocean. Additionally, based on improved backward SBR tracing process and kd-tree acceleration treatment mentioned in [15], computation time rapidly decreases. In the proposed composite scattering model, SBR method is used to calculate scattering from multiple interactions of ship-ocean and scattering from ship, whereas scattering from sea surface is calculated by SDFM. Hybrid SBR-SDFM method should be time-saving and highly efficient to handle radar cross section calculation of electrically large ship-ocean composite model.

2. SBR-SDFM Method for Composite Scattering

2.1. Composite Scattering Model of Ship-Ocean

Figure 1 shows simplified composite scattering model of ship above rough sea surface. Under framework of the SBR-SDFM method, composite scattering model is divided into four major parts.

(1) Single Scattering from Ship. This part is calculated equivalently by physical optics. During the computation, not only self-shade of ship but also shade to ship by sea surface should be taken into consideration. As shown in Figure 1, scattering calculation does not involve part of ship submerged in sea surface.

(2) Multiple Scattering from Ship. This part is calculated by SBR. Ship model usually contains numerous dihedral structures; these structures can generate large number of multiple reflections, which develop into major components of backscattered field. Hence, multiple scattering of ship should be considered.

(3) Single Scattering from Sea Surface. SDFM method is adopted to calculate scattering contribution of sea surface. In the theory of SDFM, sea wave is envisaged as a superposition of gravity wave configuration approximated by planar facets and capillary wave configuration added on the planar facets at resonant wavenumber. This assumption better conforms to the sea surface two-scale model than PO or SBR. Then, with Bragg scattering from capillary waves being considered separately, sea surface can be meshed into larger triangle patches, which can adapt well to the SBR method for ray tracing between ship and ocean. Thus, computing time is comparatively reduced. Furthermore, visibility of facets on the sea surface should be checked to eliminate shade from ship and other facets on the sea surface.

(4) Coupling Scattering between Ship and Sea Surface. This part will be discussed in Section 2.2 in detail.

Then, scattered field of composite model can be developed by creating superposition of obtained scattered fields: where , , , and are scattering field of four parts mentioned above, respectively.

2.2. Coupling Scattering between Ship and Ocean

SBR method is applied to describe coupling of ship and sea surface (see Figure 2). There are four main ray paths included: ship-sea, sea-ship, ship-sea-ship, and sea-ship-sea. First, the ship and the sea surface are meshed into triangle patches, which are numbered and stored into kd-trees, respectively. Taking ray path ship-sea as an example, given the incident wave direction, patches on the ship model illuminated directly by the incident wave will be identified (i.e., triangle T1) and stored into computer memory as an array. Then the ray starts from illuminated patches on the ship and propagates along the reflected direction, implementing intersection tests with previously built kd-tree of sea surface to determine illuminated patches on the sea surface. Once one patch is found to be illuminated (i.e., triangle S2), backward ray tracing will be performed on its neighbors (i.e., triangles S1 and S3; their IDs are always adjacent to illuminated patches) to check their visibility. Finally, calculation will be performed on far field scattering of all illuminated patches on the sea surface.

In this study, scattering field from sea surface is computed based on SDFM. Bragg scattering field from each facet can be obtained by Fuks’ model, with fair approximation that Bragg wave is a single frequency sinusoidal wave to simplify computation, and the polarization factor is modulated by gravity wave. The Bragg scattering field can be expressed as follows:where and represent the combined vectors of incident and scattered wave number, respectively; is the distance from the observation to the facet centers; is the electromagnetic wave number; and are the polarization of the incident and scattered wave; denotes the scattering amplitude that can be expressed aswhere is the relative permittivity of the sea water. is defined by , and is the electromagnetic wave number. is the polarization factor.

in (3) represents the fluctuations of the capillary on the facet, which is envisaged to be a sinusoidal wave that can result in the Bragg resonance and it can be written as where , are the wavenumber and circle frequency of the Bragg resonance and is the coordinate of the respective point on the facet. is the amplitude of the sinusoidal short wave, is the sea spectrum, and is the area of each facet.

We have performed simulation for an individual facet on the sea surface (see Figure 3) to verify the correctness of the model. Consider a sinusoidal wave expressed by , where is the size of the facet and is the spatial wavelength satisfying the Bragg resonance. The dielectric constant of the sea surface is and the incident wave frequency is set to 2.0 GHz. The incident direction is set by , . Figures 4 and 5 show numerical results of SDFM compared with those simulated by MLFMM in commercial software FEKO. Results of SDFM are consistent with those of MLFMM. Table 1 presents the computational time of results in Figure 4. The computer configuration is Intel(R) Xeon(R) CPU (2.1 GHz, 32.0 GB).

As for the ray path sea-ship, instead of patches on the ship, patches on the sea surface illuminated directly by the incident wave will be determined first, and the ray tracing process will be implemented from sea surface to ship. Finally, scattering field of illuminated patches on the ship will be calculated by PO. Scattering field approximated by PO integral can be written as follows:where ; is the current density at on the patch.

The high order coupling interaction of sea-ship-sea or ship-sea-ship is similar to the second-order interaction described above. The difference is that high order coupling interaction requires one more ray tracing process.

2.3. Multiple Scattering of Ship

There are some simulation results to verify the necessity of taking multiple scattering of the ship into account. Figure 6 displays dimensions of the ship. Figure 7 shows monostatic scattering results of the ship; results include both first-order scattering and second-order scattering.

Results present visible difference. When multiple scattering of the ship is considered, the backscattering field will be significantly enhanced.

3. Results and Discussion

3.1. Validation for the Composite Scattering Model
3.1.1. Cube on Flat Plate

The composite model of a flat plate () with a cube located over at height of 0.03 m is simulated at frequency of 10 GHz for HH-polarization (see Figure 8). Figure 9 shows bistatic RCS of the model when side length of the cube is 0.03 m, and the incident angles are set to , . The monostatic RCS is shown in Figure 10. The results of SBR method agree well with results based on MLFMM. Furthermore, higher accuracy is obtained compared with the four-path model.

3.1.2. Sphere on Rough Sea Surface

The diameter of the sphere is 0.3 m and the sphere is above the sea surface by 0.15 m (see Figure 11). The sea area is  m2, and the wind speed is  m/s. The incident angle is , whereas the frequency of the incident wave is 4 GHz. Both the sphere and sea surface are assumed to be perfect electronic conductor. Figure 12 shows the bistatic RCS results of the model for HH-polarization. The simulation results calculated by our method are in good agreement with those by FEKO-MLFMM.

3.2. Composite Scattering of the Ship-Ocean

In the following, SBR-SDFM method is applied to the RCS prediction of the battle cruiser above the sea surface. As presented in Figure 13 (dimensions of the ship are the same as those in Figure 6), large size of model causes almost impossible calculation and analysis of scattering problem using numerical algorithm at high frequencies.

In our simulation, the incident wave frequency is fixed at 10 GHz. The permittivity of the sea surface is . First, bistatic RCS for HH-polarization is calculated with the incident direction set at , , and the scattering direction is , . The wind speed is selected by  m/s. Figure 14 illustrates bistatic RCS results. For this scenario, flank of the ship and the sea surface form strong coupling, which can significantly increase the backscattered field intensity. It can be seen from the results that when scattering angle varies from to , scattering from both ship and coupling interaction preponderate in the total scattered field because of dihedral reflectors on the ship and dihedral reflectors formed by flank of the ship and the sea surface. Nevertheless, in the specular direction, scattering from the sea surface dominates the total scattering field.

As is known, the wind speed is a principal factor that influences the composite scattering field. The simulation results of bistatic RCS for the HH-polarization in different wind speeds 5 m/s, 7 m/s, and 10 m/s are shown in Figure 15. The incident wave slants at , , and the scattering direction is , . As the wind speed increases, the fluctuation of sea surface becomes stronger, which can attenuate the scattering intensity in specular direction but enhance that in diffuse direction. From the comparison results in Figure 15, one can see that scattering field intensity in most diffuse directions increases and the scattering field intensity in specular direction decreases, as the wind speed increases.

Figure 16 indicates the monostatic RCS results of the ship-ocean model at different wind speed. The incident and scattering angles are set by , . As the wind speed decreases, specular scattering from sea surface is enhanced. This phenomenon consequently amplifies the coupling interaction between ship and ocean, further increasing the composite scattering.

4. Conclusion

This paper proposes suitable and efficient hybrid SBR-SDFM method for calculation of composite scattering of electrically large ship above sea surface. SBR and SDFM are, respectively, applied to the ship and sea surface. For the key problem coupling interaction between the ship and the sea surface, SBR is used to provide accurate simulation. Results for the validation of the method indicate good agreement with exact numerical method of FEKO.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Science Foundation for Distinguished Young Scholars of China (Grant no. 61225002), the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (Grant no. 61621005), the National Natural Science Foundation of China (Grant no. 61431010), and the National Science Foundation of China (Grant no. 61501361).