“…Our starting point is a parallel, adaptive, multigrid, local quasi-Newton coupled solver for two-dimensional steady Navier -Stokes flows on staggered grids called pamg [3][4][5]. Its extension to multiphase flows, pamg-multiphase, requires a careful treatment because the differential operators are far more complex.…”