In this paper we present our results on numerical study of vapor bubbles growing in quiescent superheated liquid, as effect of liquid evaporation at the interface. Height Function interface reconstruction algorithm is coupled with an evaporation model based on continuum field representation of source terms. The flow solver is a finitevolume CFD code. Interface is tracked within a Volume-Of-Fluid framework. Continuum-Surface-Force method accounts for surface tension effects. Vapor bubble heat-transfer-controlled growth is simulated for three different working fluids: water, HFE-7100 and R134a. Accuracy of interface reconstruction algorithm is of maximum importance. Unbalance between pressure gradients and surface tension forces at interface leads to the growth of an unphysical velocity field which switches original only diffusive heat transfer mechanism to combined diffusive-convective one. Height Function algorithm reduces the magnitude of this unreal velocity field. Standard test cases are considered to assess the performances of implemented version, through comparison with the widely used Youngs algorithm.