this paper is to propose a simple and robust flux splitting finite volume solver to compute an incompressible multi-fluid model based on preconditioning type artificial compressibility [1,2]. To implement unsteady flow calculations, a dual time stepping
