In the present research, numerical simulation of the characteristic chart and steady-state Wakefield flow around a marine propeller is conducted. Solutions were performed using the open-source OpenFOAM software and the steady incompressible simple-Foam solver. The gradients were calculated using the linear Gauss algorithm, and the pressure equation was solved with the multi-grid method. In this research, characteristic chart simulation of the propeller was carried out for the entire operational conditions and the effect of using Realizable-k-ε and k-ε-v^2-f turbulence models on the results was investigated. The results were found to be in good agreement in all conditions except for the near bollard region. In this region, the propeller inlet angle of attack severely increased, and the two equation model predicted the thrust coefficient with 24% error, while implementing the four equation model significantly developed the results and decreased the error to 5%. The wake region parameters were also investigated in the numerical simulations at different longitudinal and radial cross sections behind the propeller which showed good agreement compared with the available experimental data. Wake region investigation showed that the flow behavior in downstream cross sections is similar to the corresponding upstream section with smaller variation ranges and for the swirling flow behind the propeller, the maximum and minimum angular position of the wake components rotates. The obtained results also show that the wake axial velocity component deviation is extremely large at the blade tip.