Hydrofoil shapes are used for the marine turbine blades to capture kinetic energy from water currents effectively. Predicting transitions is a critical concern when studying the hydrofoil boundary layer. This paper analyzed the transitional behavior of the boundary layer in the National Advisory Committee of Aeronautics (NACA) hydrofoil, NACA0009, with a blunt trailing edge using the Algebraic Local-Correlation-based Transition Modeling (Algebraic LCTM) model. First, through sensitivity analysis, the effects of the maximum y+ (the dimensionless distance y to the wall), grid expansion ratio, number of normal and streamlined grids, and timescale on transition prediction were studied. The results indicate that finer y+ value and appropriate grid expansion ratios can improve the accuracy of transition prediction, while the influence of timescale on the prediction results is relatively small within the range of Courant number theory values. Second, further analysis was conducted on the transition prediction performance under different Reynolds numbers. It was found that the model predictions were consistent with experimental values at low Reynolds numbers, but the predicted transition position was advanced at high Reynolds numbers, mainly because of the significant disparity in eddy viscosity coefficients within the free flow field. In the study of leading-edge roughness bands' impact on boundary layer transition for hydrofoil, the introduction of roughness significantly expedited the transition process. The Algebraic LCTM model outperformed the gamma (γ) transition model, reducing prediction errors by 5–40% for boundary layer parameters and maintaining errors between 0.005 and 4% for wake vortex shedding frequency, as opposed to the γ model's 0–23%. The results of this study provide a theoretical basis for hydrofoil design.