In the realm of marine science and engineering, hydrofoils play a pivotal role in the efficiency and performance of marine turbines and water-jet pumps. In this investigation, the boundary layer characteristics of an NACA0009 hydrofoil with a blunt trailing edge are focused on. The effectiveness of both the two-equation gamma theta (γ-Re_θt) transition model and the one-equation intermittency (γ) transition model in forecasting boundary layer behavior is evaluated. When considering natural transition, these two models outperform the shear stress transport two-equation (SST k-ω) turbulence model, notably enhancing the accuracy of predicting boundary layer flow distribution for chord-length Reynolds numbers (Re_L) below 1.6 × 10⁶. However, as Re_L increases, both transition models deviate from experimental values, particularly when Re_L is greater than 2 × 10⁶. The results indicate that the laminar separation bubble (LSB) is sensitive to changes in angles of attack (AOA) and Re_L, with its formation observed at AOA greater than 2°. The dimensions of the LSB, including the initiation and reattachment points, are found to contract as Re_L increases while maintaining a constant AOA. Conversely, an increase in AOA at similar Re_L values leads to a reduced size of the LSB. The findings are essential for the design and performance optimization of water-jet pumps, particularly in predicting and flow separation and transition phenomena.