Film boiling has various industrial applications specially in heat exchangers. Studying this phenomena on complex geometries and investigating heat transfer coefficient is designed by many industries. A front tracking method is used for capturing the interface. In this method, marker points are assumed to represent the interface. Elements are created by connecting the marker point together. The points and elements are connected to next and pervious ones. In fact, there are two type of grid. One grid is concerned with interface and the other one is relevant to computational domain. The gradient of indicator function can be determined by interpolating the front points to computational grid and with solving a poisson’s equation throughout the domain. Theromophysical properties can be obtained by indicator function. Therefore, using projection method, Navier-Stocks and energy equations are discretized and solved in case of phase change. So the temperature, pressure and velocity fields are found in the whole domain. In next step, the velocity of front is interpolated from velocity on computational grid. The new positions of the front is located based on the current position, time step and new velocity field and the new indicator function can be computed from the new position of the interface. Film boiling is simulated on one cylinder or more. The effects of Grashof number, Jacoob number on a cylinder, effects of separation distance, relative angels and diameter for two cylinders and effects of different configuration of multicylinders in staggered and normal array and effects of number of rows of cylinders is evaluated. The nusselt number is computed on every cylinder. Key words : Film boiling, Front tracking method, Complex geometries, Heat transfer