The current study is focused on investigating a Jacobian-Free Newton-Krylov (JFNK) method to obtain a fully-implicit solution for two phase flows. In the JFNK formulation, the Jacobian matrix is not directly determined ...