This thesis is concerned with the application of adaptive local quasi-Newton coupled multigrid (ALQNMG) solvers to the numerical simulation of viscous incompressible fluids, using the multi-fluid model. The ALQNMG ...