The present study is focused on modeling of dynamic stall behavior of a pitching airfoil. The deep stall regime is in particular considered. A model is proposed, which has a low implementation and computational complexity but yet is able to deal with different types of dynamic stall conditions, including those characterized by multiple vortex shedding at the airfoil leading edge. The proposed model is appraised against an extensive data set of experimental (α,CL) curves for NACA0012. The results of an existing widely used model, having comparable complexity, are also shown for comparison. The proposed model is able to well reproduce not only the classic curves of deep dynamic stall but also the curves characterized by lift oscillations at high angles of attack due to the shedding of multiple vortices. Furthermore, the model appears to be robust to variations of its parameters from the optimal values and of the airfoil geometry. Finally, the model is successfully implemented in a commercial CFD software and applied to the simulation of a vertical axis wind turbine within the actuator cylinder approach. The accuracy of the prediction of the turbine power coefficient in the whole rotation cycle is very good for the optimal working condition of the turbine, for which the model parameters were calibrated. Fairly good accuracy is also obtained in significantly different working conditions without any further calibration.