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.