Close-packed contra-rotating vertical-axis turbines have potential advantages in wind and hydrokinetic power generation. This paper describes the development of a numerical model of a vertical axis turbine with a torque-controlled system using an actuator line model (ALM). The developed model, coupled with the open-source OpenFOAM computational fluid dynamics (CFD) code, is used to examine the characteristics of turbulent flow behind a single two-bladed vertical-axis turbine (VAT). The flow field containing the turbine is simulated by solving the unsteady Reynolds-averaged Navier-Stokes (URANS) equations with a k - ω shear stress transport (SST) turbulence model. The numerical model is validated against experimental measurements from a two-bladed H-type wind turbine. Turbine loading is predicted, and the vorticity distribution is investigated in the vicinity of the turbine. Satisfactory overall agreement is obtained between numerical predictions and measured data on thrust coefficients. The model captures important three-dimensional flow features that contribute to wake recovery behind a vertical-axis turbine, which will be useful for future studies of close-packed rotors with a large number of blades.