One-dimensional two-fluid model for wavy flow beyond the Kelvin-Helmholtz instability: Limit cycles and chaos

A 1D TFM numerical simulation of near horizontal stratified two-phase flow is performed where the TFM, including surface tension and viscous stresses, is simplified to a two-equation model using the fixed-flux approximation. As the angle of inclination of the channel increases so does the driving body force, so the flow becomes KH unstable, and waves grow and develop nonlinearities. It is shown that these waves grow until they reach a limit cycle due to viscous dissipation at wave fronts. Upon further inclination of the channel, chaos is observed. The appearance of chaos in a 1D TFM implies a nonlinear process that transfers energy intermittently from long wavelengths where energy is produced to short wavelengths where energy is dissipated by viscosity, so that an averaged energy equilibrium in frequency space is attained. This is comparable to the well-known turbulent stability mechanism of the multi-dimensional Navier-Stokes equations, i.e., chaos implies Lyapunov stability, but in this case it is strictly a two-phase phenomenon.