In this paper, the characteristics of Gaussian beam propagation through turbulent plasma sheath are studied. According to the generalized Huygens-Fresnel principle, the random phase screen is generated by power spectrum inversion method based on the fast Fourier transform. The random phase screen is used to simulate turbulence effect. The thickness of the plasma sheath is of about centimeter order of magnitude. Compared with the single fast Fourier transform algorithm, the double fast Fourier transform algorithm is not prone to under-sampling and can obtain good image results, even if the diffraction distance is 1 mm. Therefore, double fast Fourier transform algorithm is used for investigating the beam propagation between two phase screens. The turbulence effect of the plasma sheath surrounding a hypersonic vehicle is simulated by the multi-random phase screens. When the flight altitude is 45 km and the flight speed is 18 Mach, the intensity of refractive index fluctuation variance ranges from ${10}^{ - 11} $ to $ {10}^{ - 14} $ indicated by analyzing the plasma flow field around the hypersonic vehicle. The characteristics of the Gaussian beam propagation through the turbulent plasma are numerically simulated. The results show that the refractive index fluctuation, wavelength and transmission distance are important factors affecting the Gaussian beam quality. The larger the refractive index variance, the more severe the spot dispersion and the more obvious the light intensity fluctuation. As wavelength is longer, the ability of the Gaussian beam to resist turbulence becomes stronger and the dispersion of the light spot and the intensity fluctuation are smaller. The beam distortion and the spot dispersion become more severe as the transmission distance is longer.