The paper introduces a numerical simulation method for Synthetic Aperture Radar (SAR) imaging of submerged body wakes by integrating hydrodynamics, electromagnetic scattering, and SAR imaging simulation. This work is helpful for better understanding SAR images of submerged body wakes. Among these, the hydrodynamic model consists of two sets of ocean dynamics closely related to SAR imaging, namely the wake of the submerged body and wind waves. For the wake, we simulated it using computational fluid dynamics (CFD) numerical methods. Furthermore, we compared and computed the electromagnetic scattering characteristics of wakes under various navigation parameters and sea surface conditions. Following that, based on the operational principles and imaging theory of synthetic aperture radar (SAR), we established the SAR raw echo signal of the wake. Employing a Range-Doppler (RD) algorithm, we generated simulated SAR images of the wake. The results indicate that utilizing Computational Fluid Dynamics (CFD) numerical methods enables the simulation of wake characteristics generated by the motion of a submerged body with different velocities. The backscattering features of wakes are closely associated with the relative orientation between the wake and the radar line of sight. Under specific wind speeds, the wake gets masked within the sea surface background, resulting in less discernible characteristics of the wake in SAR images. This suggests that at lower speeds of submerged body or under specific wind conditions, the detectability of the wake in SAR images significantly diminishes.

The airborne Synthetic Aperture Radar (SAR) has the advantages of high-precision real-time observation of wave height variations and portability in the high frequency band, such as the Ku band. In view of the Four Fast Fourier Transform (4-FFT) algorithm, combined with a Gaussian operator, a Laplacian of Gaussian (LoG) Phase Unwrapping (PU) expression was derived. Then, an Adaptive LoG (ALoG) algorithm was proposed based on adaptive variance, further optimizing the algorithm through iteration. Building the models of Kelvin wake on the sea surface and height to phase, the interferometric phase of wave height can be simulated. These PU algorithms were qualitatively and quantitatively evaluated. The Principal Component Analysis (PCA) scores of the ALoG iteration (ALoGI) algorithm are the best under the tested noise levels of the simulation. Through a simulation experiment, it has been proven that the superiority of the ALoGI algorithm in high spatial resolution inversion for the sea-ship surface height of the Kelvin wake, with good stability and noise resistance.

Microwave camera provides 3-D high-resolution radar images at video frame rates, enabling the capture of dynamic target features. Multiple-input–multiple-output (MIMO) array-based 3-D radar imaging system requires fewer antennas, which effectively reduces hardware costs. Due to the limited computational resources of the miniaturized MIMO microwave camera, real-time processing of a large amount of 3-D echo data requires an imaging algorithm that has both real-time performance and large output spatial extent. This paper presents the limited output spatial extent and spatial aliasing in existing MIMO wavenumber domain algorithms through theoretical derivation and simulation. To suppress aliasing while expanding the output spatial extent, an optimization approach for the wavenumber domain algorithms is proposed. The improved wavenumber domain algorithms divide the target area into multiple sub-blocks, and a broader range of imaging results is obtained through independent imaging of the sub-blocks and a spatial aliasing suppression filter. Simulation results show that the improved wavenumber domain algorithms effectively suppress the aliasing energy of each sub-block while maintaining the advantage of low time complexity. Expansion of output spatial extent in existing MIMO wavenumber domain algorithms is achieved.