Abstract Wave observations are as yet an important technique to resolve fine structures in the plasma that cannot be identified by different instruments. In this way, we proposed the generation of electrostatic solitary waves during the interaction between the solar wind particles and Venus’s atmosphere at high altitude. The plasma system is treated as a multicomponent unmagnetized plasma consists of background electrons, positive ions (H+ and O+), and streaming solar wind electrons and protons. The dispersion relation is derived for linear waves and the stability/instability of electrostatic wavepackets is investigated. The dependence of the instability growth rate on the ion beam speed has been analyzed. Stability analysis reveals the occurrence of an imaginary frequency part in three regions. Using the reductive perturbation theory, the set of fluid equations is reduced to the Korteweg-de Vries (KdV) equation to describe the evolution of small but finite amplitude soliton waves. The model predicts the propagation of electrostatic soliton waves with an electric field of 1.2 mV/m and a time duration of 0.4 ms. The output of the fast Fourier transform of the electric field pulse is a broadband in the frequency range of ~ 3.2 to 199.5 kHz. It is proposed that the model can be a decent possibility for clarifying the observation of a wide variety of plasma oscillations by Galileo flyby of Venus.