Disclosed is a technique for simulating wind interaction with wind turbines. A turbine blade is divided into radial sections. The effect that each of these radial sections has on the velocities in Eulerian computational cells they overlap is determined. The effect is determined using Lagrangian techniques such that the calculations need not include wind components in the radial direction. A force on each radial section of turbine blade is determined. This force depends on the axial and azimuthal components of the fluid flow in the computational cell and the geometric properties of the turbine blade. The force on the turbine blade is fed back to effect the fluid flow in the computational cell for the next time step.