Analysis of a 1D incompressible two-fluid model including artificial diffusion

Abstract

This article examines a 1D incompressible two-fluid model including artificial tensor diffusion. The aim is to obtain a formulation that provides convergent numerical solutions for all flow conditions within the stratified and the stratified wavy flow regime. With appropriate simplifications, the two-fluid model reduces to one momentum balance, one mass conservation and two algebraic equations. It has previously been established that a formulation that is well posed in possessing exclusively real characteristics can be obtained by including an axial diffusion term in the momentum balance. In this article, however, we demonstrate that this is not sufficient to obtain a system suitable for numerical simulations. Although the unbounded growth rates of the standard two-fluid model are eliminated, linear stability theory predicts that infinitesimal wavelengths still experience finite growth. This entails that grid refinement always will result in new unstable wavelengths being resolved. On the other hand, if artificial axial diffusion is added to both the mass and the momentum equations as suggested here, a cut-off wavelength is established below which all wavelengths are stable. Thus, a numerically converging model is formed, which retains the long-wavelength properties of the standard two-fluid model. The conclusions of the mathematical analysis are substantiated by numerical simulations of ID gravity waves.