A three dimensional coupled VOF and Level set (VOSET) method with and without phase change on general curvilinear grids. (21st September 2020)