A numerical formulation for the Conservative Level-Set Method of Olsson and Kreiss is implemented using the finite volume approach. Both convective term of the interface advection equation and compression term of the reinitialization equation are discretized by the first order upwind plus flux limiter scheme. Time integration is solved by the TVD Third Order Runge Kutta method, whereas gradients are computed by the linear Least-Squares approach. The method is tested computing moving interfaces with external velocity fields on unstructured meshes. Additionally, Level-Set method was coupled to a variable density Navier-Stokes solver for incompressible flow TERMOFLUIDS. Physical properties are assumed to vary smoothly in a narrow band around the interface to avoid numerical instabilities. Pressure velocity coupling is solved using a two step projection method where in the first step the convection, diffusion and source terms are computed explicitily using the Adams-Bashforth time integration. In the second step the velocity is corrected by adding the gradient of the pressure variable. Numerical solution of pressure Poisson equation is calculated with an efficient preconditioned conjugate gradient solver. Mass conservation and convergence properties of the model are investigated for cases with large density and viscosity ratios. For the purpose of validation, numerical results are compared with numerical and experimental data for standard benchmark problems.

Peer Reviewed