Coastal vegetation plays a key role in wave attenuation on sloping beaches. However, accurately predicting the vegetation-related drag coefficient (CD) under breaking waves on slopes remains a challenge. In this study, the Reynolds-averaged Navier–Stokes equations were modified by incorporating a vegetation resistance term to consider the vegetation-induced momentum loss. The stabilized k-ω shear stress transport turbulence model was improved by introducing the vegetation-induced turbulence generation and dissipation terms. The numerical model was validated by two laboratory experiments. Based on experimental wave height data, the CD was determined via a calibration method. The numerical results were found to be highly sensitive to the CD. The calibrated CD achieved excellent agreement with experimental data (RMSE < 0.2). In addition, the calibrated CD showed a high correlation (R2 = 0.85) with the Iribarren number ξ. Finally, a CD prediction formula expressed by the combination of Reynolds number and ξ was proposed via regression analysis, and its performance was evaluated by two additional cases. The limitations and applicability of the CD prediction formula were also discussed. The proposed numerical model and CD prediction formula have the potential to be useful for the numerical simulation of wave interactions with vegetation on slopes.