# Geometry

To study the oscillation of a levitated zirconium drop in helium atmosphere, a 2D axisymmetric model is created as shown in Fig. 2. A rectangular domain (4 mm x 2 mm) which includes both zirconium and helium was created. On the left edge of the rectangular, a semicircle with diameter of 2.5 mm was created to define the zirconium drop.

# Numerical Solver

The pressure-based solver is used for the numerical simulation. The reason for this is the use of volume of fluid model for multiphase simulations discussed later. Also, a pressure-based solver is used as the fluid can be assumed to be incompressible. The pressure based solver satisfies mass conservation by solving the pressure equation. The transient time solver is implemented and simulation is carried out until the desired time-length of data is obtained to analyze and extract frequency from it.

# Multiphase Modelling

Volume of fluid model consisting of two phases is used for the simulations. Volume of fluid method is used for problems where the multiple phases do not mix and form a defined boundary at all time. Volume of fluid method solves single momentum equation for all the phases. A new variable whose value is between 0 and 1 is introduced which traces the volume fraction of all the phases in each control volume. The properties of the individual phases at a particular time step are function of volume fraction variable and global property at given time step. Volume fraction parameters were solved using explicit scheme.