Abstract: The isothermal dispersed bubbly flow in a vertical tube was simulated by Eulerian two-fluid method with two different turbulence models (RANS and LES). With proper inter-phase interaction models, the time averaged results of flow filed agree well with experiment results. The wall peak bubble distribution is captured and two-phase velocity fields are accurate. Compared with RANS SST turbulence model, the LES with WALE sub-grid model can get better agreement, and it also gives transient and detail information of the flow field.