An aerodynamical model for studying three-dimensional flow fields about wind turbine rotors is presented. The developed algorithm combines a three-dimensional Navier-Stokes solver with a so-called actuator line technique in which the loading is distributed along lines representing the blade forces. The loading is determined iteratively using a blade-element approach and tabulated airfoil data. Computations are carried out for a 500 kW Nordtank wind turbine equipped with three LM19.1 blades. The computed power production is found to be in good agreement with measurements. The computations give detailed information about basic features of wind turbine wakes, including distributions of interference factors and vortex structures. The model serves in particular to analyze and verify the validity of the basic assumptions employed in the simple engineering models.