A weakly nonlinear (WN) theoretical model for the two-dimensional incompressible Kelvin–Helmholtz instability (KHI) is proposed. Its solution form is the complete expansion in real space. The transition from linear to nonlinear growth is analytically studied via third-order solutions of plane KHI initiated by a single-mode surface perturbation. The difference between the WN growth of the Rayleigh–Taylor instability (RTI) in plane geometry and the WN growth of the KHI in plane geometry is discussed. It is found that there are resonance solutions in the higher harmonics of KHI but not in RTI. The vertex of spikes and bubbles is deflected because of the shearing effect in KHI compared with that in RTI. The spike velocity increases with the Atwood number before a particular time and then inversely after that time. There is no such reversal in RTI. However, the bubble velocity with the Atwood number of KHI has the same evolution rule as that of RTI. In addition to the influence of initial perturbation on the nonlinear saturation amplitude, the influence of Atwood number on the nonlinear saturation amplitude of the fundamental mode is obtained. Finally, compared with the numerical simulations, at the same initial conditions, the analytical solutions up to the third harmonics are well consistent with the simulated datum in the linear and weakly nonlinear stages for a widely varied Atwood number.
Read full abstract