An analytic methodology was developed for free surface motions between liquid metal coolant and cover gas in order to calculate the phenomena of gas entrainment in hot pool surface through IHX EMP and reactor core. The methodology was setup by applying the first order VOF convection model to CFX4 general purpose fluid dynamics analysis code. The methodology was validated by applying it to an experimental apparatus designed for free surface motions of KALIMER reactor. The distributions of free surface calculated by the present methodology were almost coincident with the experimental data. The developed methodology was applied to the KALIMER reactor of full power operating condition. The shapes of the free surface were nearly uniform. From the results, it was found that the altitude of the free surface from the IHX inlet nozzle of KALIMER reactor is high enough not to affect to free surface motions of generating gas bubbles from the turbulent shear flows such as hydraulic jump and water falls.