Flow boiling heat transfer is routinely encountered in nuclear reactors, steam engines and other engineering applications. Although several researchers have carried out different numerical and experimental investigations on flow boiling, the underlying physics of the interfacial interaction is still a complex phenomenon to understand in detail. Hence, the numerical simulation and optimisation regarding the adoption of engineering flow boiling parameters have been conducted in this study using the Modified Eulerian-Eulerian Model (MEEM) and Wall Heat Balance Model (WHBM). To predict interpenetrating flow fields and to provide detailed relevant information on the flow behaviour, this study considered a uniform axial heating profile for a cylindrical flow channel. The Raynolds Average Navier Stokes (RANS) equation with an appropriate turbulence model are used to predict the effect of turbulence on the mean flow field, while the MEEM multiphase sub-models are employed to predict the temperature distribution along the wall, the average void fraction, tracking of the single bubble detachment diameter, heat balance at the wall, effect of surface roughness on heat transfer, the effect of aspect ratio, and the critical heat flux. The results obtained from this study are compared with the selected numerical investigation and experimental data presented in the open literature. The present study shows a better approximate prediction (with minimal uncertainties) of both the subcooled boiling heat transfer and the saturated boiling heat transfer. In summary, this study agrees with extant theories and experimental predictions. Thus, it has provided more profound insights into flow boiling heat transfer particularly for flow in a vertical pipe.