# Solution Methods

PISO scheme with neighborhood correction factor of 1 is used for pressure-velocity coupling. Spatial discretization schemes are selected suitable for the Volume of fluid method. The gradient is calculated using least squares cell based method. This algorithm is the most accurate algorithm available with competent computation time. Pressure and momentum equations are solved using the PRESTO! And QUICK algorithms respectively. Geo-Reconstruct scheme is used for the volume fraction calculations. This scheme calculates the volume fraction of the phases in each cell and regenerates it using a straight line with slope in each cell. The scheme is very useful in problems where the interface among phases is to be traced accurately. A very smooth interface is generated using Geo-Reconstruct scheme. For the transient formulation, non-iterative time advancement was selected which means all the fluid property equations were solved once per time step.

# Mesh

An Eulerian approach is used in the simulation. In the Eulerian approach, rather than following all the particles in a fluid motion, fixed number of points in space are selected and properties of fluid at a particular time step are recorded. In the Eulerian approach, the fluid properties are a function of space and time. The rectangular domain was divided into a grid mesh with each computational cell of equal volume. The mesh size plays an important part in the computation time of the simulation. Considering the accuracy of solution and computation time, a mesh was generated such that each control volume was square with side of 25 microns.