This study uses the large eddy simulation model to analyze the heliostat's mean and peak aerodynamic loads due to the atmospheric boundary layer. A synthetic method based on the Fourier technique called consistent discrete random flow generation (CDRFG), is proposed. With the CDRFG technique, key flow parameters, including mean velocity profile, turbulent intensities, integral length scales, and turbulent spectra generated in wind tunnels, can be replicated while also satisfying the divergence-free condition. A three-facet heliostat with an elevation angle of a = 45° and the rear aligned to the inflow was analyzed. The heliostat behaves like a lifting surface in this orientation, accentuating the aerodynamic effect. The methodology proposed in this study can accurately reproduce flow statistics and predict the peak loads. Compared to experimental data, differences of 6.92% for drag, 12.59% for lift, and 6.72% for overturning were observed. Furthermore, the simulation reveals the generation of wingtip vortices on the sides of the heliostat, which contribute to the aerodynamic load. Overall, this technique has been demonstrated to be effective in replicating the atmospheric boundary layer and predicting the aerodynamic coefficients of heliostats.