10-30-2008 12:13 AM
10-30-2008 11:38 AM
10-30-2008 12:07 PM - edited 10-30-2008 12:07 PM
I think I have an idea of where the problem is. Unfortunately, I need to do some of my own work, so I may not be able to work through a solution.
The formulas all use a single data point except for the integration which calculates the angular velocity Wm. Your VI is working with arrays of wind speed. Those arrays have many points initialized to values which do not represent the actual data, but are still used by the integration and other parts of the program.
I think if you can convert the Wm calculation to work with incremental data, then the model should work. Integration simply looks at all the data since the start of recording. In the physical world the angular speed is an indication of past forces applied to the turbine (Velocity = integral of acceleration over time). But the velocity at time t can be calculated from the velocity at t-1 and the force (torque) at t.
I hope to be able to do something with this later today or tomorrow.
Lynn
10-30-2008 01:22 PM - edited 10-30-2008 01:24 PM
Thanks for your kidness,
I understand what you say, you are right. I'm trying to make a discrete representation of the model, for programming it into a Matlab Script block. maybe works.
Thank you
10-30-2008 03:43 PM - edited 10-30-2008 03:44 PM
I had no problems seeing the formulas, but others such as Mathan might.
Don't put your images on a 3rd party server.
Attach them to your message then post it. Then edit your message to add the image link that points to the location on the NI server. Then everyone should be able to see it.![]()
For security reasons, many users may have firewalls that block the images coming from 3rd party sources.
10-30-2008 03:53 PM
OK. I think I have something which is close to what you want.
All the calculations are done on the current value of the wind speed. This what a real system has: the current wind and the angular velocity which was imparted by the previous wind.
Using [i] to designate the values calculated from the element i of the array of wind speeds Vm, I have the following equations:
Lambda[i] = R*Wm[i-1]/V[i]. Note that this is not the Lambda(sub)i from the equations you posted. Lambda[i] is the Tip Speed Ratio calculated for interval i.
Cp[i] = function of Lambda[i] as given in your post above.
Pm[i] = 0.5*A*rho*Cp[i]*(Vm[i])^3
Tm[i] = Pm[i]/Wm{i-1]. If Wm[i-1] = 0 then a finite value is substituted. See discussion below.
Wm[i] = Wm[i-1] + Tm[i]/It. This is the incremental change in the velocity. The "integration" of the acceleration term (Tm/It) is based on di = 1.
When I implement this I get results which are very similar to your Matlab results except for lambda.
Comment of the Torque calculation: The given formula for torque is valid, however it does not represent the physics of the production of the torque. The torque is not produced by the power generated, it is produced by the change of momentum of the air passing the turbine blades. If the formula for the torque was based on the wind and some geometric factors, the infinities produced by dividing by zero could be avoided. I did a brief search but could not find any straightforward formula.
Lynn
10-30-2008 04:29 PM - edited 10-30-2008 04:30 PM
Amazing!!, Thank you very much,
10-30-2008 04:48 PM
I run my Matlab simulation again, and the result of lambda is the same, now there is no objection to your work.
thanks a lot.
10-31-2008 09:19 AM
Glad I could be of assistance.
Sometimes someone with a different way of looking at things can offer a suggestion which would not have occurred to the person who was originally involved with the problem.
Lynn
10-31-2008 10:59 AM
Yes, your solution is wonderful.
Now I have to keep working, the torque signal must be sent to a frequency inverter which controls an induction machine, so to have a tangible form of torque produced by the turbine.