As the lander approaches the lunar surface, the engine plumes impinge on the lunar regolith and entrain lunar dust from the surface. This plume–surface interaction and the resulting dispersion of lunar dust form a multi-physics, multi-scale problem, which becomes even more complex under multi-engine conditions. This study employed the direct simulation Monte Carlo method to simulate the plume–surface interaction flow field of a four-engine lunar lander at various landing altitudes and lunar surface angles. Flow characteristics were analyzed, and the impact of the plume and backflow on the lander was assessed. Subsequently, lunar dust simulation was conducted using the plume field as a basis. The study determined the spatial distribution of particles with different diameters at various landing altitudes and surface angles, as well as their impact velocities on the lander. Furthermore, taking into account the variations in the lander's altitude and attitude, a dynamic simulation of lunar dust during the landing process was conducted. This process resulted in the dynamic distribution of lunar dust during landing, laying the groundwork for real-time simulation of lunar dust distribution and reliable visualization during landing simulations. These findings are valuable for assessing and mitigating the hazards posed by lunar dust.