Cp = 0.7737 * 4.0 * 1.0/3.0 * math.pow((1 - 1.0/3.0), 2)
Cp = np.array([Cp, Cp])
# Define site measurements
windDirection = 270
airDensity = 1.1716
# Define turbine locations
dist = 3.0
turbineX = np.array([0.0, dist*rotorDiameter[0]])
turbineY = np.array([0.0, 0.0])
res = 200
Ueff = np.zeros(res)
y = np.linspace(-1.5*rotorDiameter[0], 1.5*rotorDiameter[0], res)
theta = np.zeros_like(y)
for i in np.arange(0, res):
turbineY[1] = y[i]
_, temp = jensen_topHat(Vinf, rotorDiameter, axialInd, turbineX, turbineY, k, windDirection, Cp, airDensity)
print temp
Ueff[i] = temp[1]
theta[i] = y[i]/dist
print y[i], theta[i]
plt.figure()
plt.plot(theta, Ueff/Vinf)
plt.show()
# Display retunrs
from jensen_topHat import jensen_topHat, powerCalc import numpy as np import math # Define turbine locations turbineX = np.array([1164.7, 947.2, 1682.4, 1464.9, 1982.6, 2200.1]) turbineY = np.array([1024.7, 1335.3, 1387.2, 1697.8, 2060.3, 1749.7]) # Define turbine characteristics Cp = 0.7737 * 4.0 * 1.0/3.0 * math.pow((1 - 1.0/3.0),2) Cp = np.array([Cp, Cp, Cp, Cp, Cp, Cp]) axialInd = np.array([1.0/3.0, 1.0/3.0, 1.0/3.0, 1.0/3.0, 1.0/3.0, 1.0/3.0]) rotorDiameter = np.array([126.4, 126.4, 126.4, 126.4, 126.4, 126.4]) # Define turbine measurements Vinf = 8.0 k = 0.1 # Define site measurements windDirection = 240 airDensity = 1.1716 wt_power, wt_velocity = jensen_topHat(Vinf, rotorDiameter, axialInd, turbineX, turbineY, k, windDirection, Cp, airDensity) print "effective windspeeds (m/s): ", wt_velocity print "Power of each turbine: ", wt_power print "Wind farm total power: ", np.sum(wt_power)
