If you're still working on this project, send me a PM. I was tasked with designing a horizontal axis wind turbine for my sophomore college project. We were required to write a simulation in Matlab and I found my code to be accurate within ~35% (hopefully most of this discrepancy comes from error within the motor efficiency).
The model we used was basically conservation of energy. The air speed was calculated as it was approaching the blade, and then after it left the blade. Because we were limited to using only foam-board as the blades, the values of lift and drag coefficients were approximated as rectangular prisms. The lift force was integrated (in the code we just broke each blade into ~1000 small rectangles using for loops instead of using numerical integration) along each blade and converted to torque. The total torque was summed and there we have our power curve! Now to optimize it we just adjust the load of the turbine so that it reaches its peak power.
The hardest part was definitely calculating 'a' and 'a prime', two values that determine the wind speeds (I'm not sure if these are just random variables assigned by our professors to help us understand it more easily or if they're actually used in aerodynamics). I don't remember the equations off the top of my head, but in order to calculate the wind velocities there was two pretty complicated equations that were required to converge. I spent hours on end trying to find the requirements for convergence, but I never did. I ended up just tossing out the divergent solutions and the final product wasn't too bad.
Let me know if you'd like any more specific details, or if you want a copy of the code we used (you must have Matlab to view it).