Simulation and analysis of two phase flow that crosses over tube bundles is crucial in safety analysis and design of kettle reboilers and steam generators. The geometry complexity of the tube bundle flow field increases the difficulty of the conventional numerical analysis. One of the methods to reduce the numerical calculations cost, is to use the porous media theory instead of the complete tube bundle modeling. Drag and tube bundle resistance force equations have been used in the porous media analysis. Based on available experimental results, two tube bundle arrangements have been considered. Due to existence of symmetric geometry and uniform energy source over the tube bundle, the two dimensional symmetric models has been used as well. It was observed that the predicted pressure drop in this research has acceptable adaptation with the experimental results. Meanwhile, by considering different outlet boundary conditions, calculated void fraction is compared to the experimental results and showed better accuracy than similar CFD research. It was observed that the enhancement of the tube bundle thermal power increases the void fraction in the heating area of the reboiler.