“…While the vast majority of related documented studies rely on explicit and segregated methods (such as projection methods), we describe here a fully coupled method for dealing with velocity-pressure coupling, to take advantage of its robustness, principally for the simulation of two-phase flows at large density and viscosity ratios, which has been little studied in the literature. Furthermore, it has been demonstrated that the fully-coupled solver hence preserving the consistency of the discretized system with the continuous equations at large contrast ratio in most cases, which allows to simplify the implementation and to potentially overcome the momentum conserving method contributions ( Sussman [32], Raessi [33], Le chenadec [34], Desjardins [35], Ménard [36], ...). To characterize the interface location, a PLIC VOF approach is used in the present work.…”