Rod bundle flow is an important research field related to reactor cooling in nuclear power plants. Owing to the rapid development of computerized performance assessments, interest in coolant flow analysis using computational fluid dynamics has garnered research interest. Rod bundle flow research data compared with experimental results under various conditions are thus needed. To address this, a boiling model verification study was conducted with reference to experiments. This study adopts the Reynolds-averaged Navier–Stokes equation, a practical analysis method compared to direct numerical simulation and large eddy simulation, including turbulence modeling, to predict the flow of coolant inside a rod bundle. This study also investigates void behavior in low-pressure subcooled flow boiling using a Eulerian approach (two-fluid model). The rod bundle has a length of 0.59 m and a hydraulic diameter of approximately 14.01 mm. At the cross-section at a height of 0.58 m, near the exit, numerical results were compared with the experimental values of the volume fraction of vapor and interfacial area concentration. The simulation results showed good agreement with the experimental data for six different initial conditions with constant density.