A simple method to compute standard two-fluid models

We describe in this paper, a tool to compute approximate solutions of standard two-fluid models with an equilibrium pressure assumption. The basic approach takes its grounds in the two-fluid two-pressure formalism, and takes advantage of the relaxation techniques. The method may be used to compute either the single- or the two-pressure model, depending on the size of mesh, which is used. It is also shown on the basis of a simple numerical experiment that the local equilibrium assumption may lead to a blow-up of the numerical solution on fine meshes, even if one accounts for drag stabilizing effects.