This paper provides a comprehensive overview of the numerical methods used to model the solar wind, integrating fluid, kinetic, and hybrid perspectives. Beginning with the foundations of solar wind theory and the development of magnetohydrodynamic (MHD) models, we discuss how fluid-based formulations enable the simulation of global structures such as coronal mass ejections, shocks, and large-scale variations in the heliospheric magnetic field. To address processes that fall outside the scope of MHD, we examine kinetic modeling based on the Vlasov–Maxwell equations, emphasizing its capability to reproduce non-Maxwellian particle distributions, wave–particle interactions, temperature anisotropies, and collisionless heating. Hybrid approaches that merge MHD with kinetic techniques are highlighted as essential tools for capturing the multi-scale nature of the solar wind, particularly in regions where macroscopic flows couple to microphysical dynamics. The paper further reviews major numerical strategies used in solar wind simulations, comparing explicit and implicit time integration, adaptive mesh refinement, Particle-in-Cell (PIC) methods, and semi-Lagrangian approaches. Key stability considerations—including boundary-condition selection, the Courant–Friedrichs–Lewy (CFL) constraint, appropriate spatial and velocity-space resolution, and the targeted use of artificial diffusion—are discussed in relation to their impact on accuracy and robustness. Example simulations demonstrate the ability of advanced models to reproduce observed proton and electron temperature profiles from the Sun out to 1AU. Overall, numerical modeling plays a central role in interpreting solar wind observations and predicting space-weather conditions, and ongoing advances in computational methods continue to strengthen our understanding of heliospheric plasma dynamics.
Somayeh Taran (Mon,) studied this question.