Film boiling is a heat transfer mechanism that might appear in different processes, such as cryogenics, metallurgical, and nuclear reactors during abnormal operating conditions, as happens during a loss of coolant accident. In this research, film boiling around a finite vertical cylinder was studied by means of computational fluid dynamics (CFD) simulations, considering six cases that include three different levels of surface temperature and two different shapes of the cylinder ends: flat and hemispherical ends. Volume of fluid (VOF) method for the treatment of multiphase flow was used, and a user-defined function was programed to consider the exchange of mass and energy between the phases. The simulations were performed with a vertical cylinder of 32 mm in diameter and 32 or 64 mm in high for flat ends or hemispherical ends, respectively, placed in a two-dimensional axisymmetric domain of 0.125 m × 0.25 m. Results obtained for the heat flux show a periodic fluctuating behavior in time as a consequence of periodical variations in the thickness of the vapor film around the cylinder. A wavy liquid–vapor interface is observed as is reported in the experimental works. The simulations results are compared with the experimental values reported in literature as well as with values obtained from correlations. The results show that the computational code used captures reasonably well the physics involved in the film boiling, being obtained that average heat flux to the case of hemispherical ends is 15.6% higher than for the case of flat ends, versus 15.2% showing experiments and 1.6% calculated combining correlations for the individual surfaces. It shows that use of correlations in this way is not appropriate in film boiling because it does not take into account the interactions between the different surfaces.