Reliable and accurate short-term prediction of wind speed at hub height is very important to optimize the integration of wind energy into existing electrical systems. To this end, a coupled model based on the Weather Research Forecasting (WRF) model and Open Source Field Operation and Manipulation (OpenFOAM) Computational Fluid Dynamics (CFD) model is proposed to improve the forecast of the wind fields over complex terrain regions. The proposed model has been validated with the quality-controlled observations of 15 turbine sites in a target wind farm in Japan. The numerical results show that the coupled model provides more precise forecasts compared to the WRF alone forecasts, with the overall improvements of 26%, 22% and 4% in mean error (ME), root mean square error (RMSE) and correlation coefficient (CC), respectively. As the first step to explore further improvement of the coupled system, the polynomial chaos expansion (PCE) approach is adopted to quantitatively evaluate the effects of several parameters in the coupled model. The statistics from the uncertainty quantification results show that the uncertainty in the inflow boundary conditions to the CFD model affects more dominantly the hub-height wind prediction in comparison with other parameters in the turbulence model, which suggests an effective approach to parameterize and assimilate the coupling interface of the model.