A finite difference numerical method, based on the VOF approach for trackin
g interface distortions, is presented. It is capable of accurately simulati
ng the fluid flow of multiple immiscible fluids for metallurgical applicati
ons. This volume tracking method is based on piecewise linear reconstructio
ns of interfaces, density distributions based on a shifted grid approach, a
nd a fully kernel-based CSF method for surface force modelling. Second orde
r temporal and spatial accuracy are achieved using improved Euler time-step
ping enhancement of a two-step projection algorithm, supported by a multigr
id-preconditioned GMRES solver that enabled large density ratios (1.30 000)
between the fluids and fine scale flow phenomena to be resolved. The code
was used to simulate the rise of an air bubble in water and in liquid pig i
ron and was able to capture the time dependent oscillation of the bubble. T
he bubble velocity varied with the instantaneous shape of the bubble. The a
veraged terminal velocity of the gas bubble in water was in good agreement
with published experimental data. Splash formation from a top submerged gas
injection lance was simulated to illustrate the capability of the code in
resolving the break up and fragmentation of liquid drops for possible use i
n the study of bath smelting processes.