Esempio n. 1
0
    else:
        q = (solidity-soln1)/(soln2-soln1)

    param = np.array([men1,sdv1,rat1,wdt1,spr1,scl1,men2,sdv2,rat2,wdt2,spr2,scl2,men3,sdv3,rat3,wdt3,spr3,scl3,men4,sdv4,rat4,wdt4,spr4,scl4,p,q])

elif cfd_data == 'vort2':
    loc1,loc2,loc3,spr1,spr2,skw1,skw2,scl1,scl2,scl3 = vorticity2(tsr,solidity)
    param = np.array([loc1,loc2,loc3,spr1,spr2,skw1,skw2,scl1,scl2,scl3])

    # print param
    # import time
    #
    # time.sleep(10)

elif cfd_data == 'velo2':
    spr1,pow1,pow2,spr2,skw,scl1,scl2,scl3 = velocity2(tsr,solidity)
    param = np.array([spr1,pow1,pow2,spr2,skw,scl1,scl2,scl3])



vel = velocity_field(xt,yt,xt + x0,yt + y0,velf,dia,tsr,solidity,cfd_data,param)

print '\nNormalized velocity at (',x0,',',y0,') from the turbine =',vel,'\n' # output velocity (normalized by free stream wind speed)

## Plotting
fs = 25 # 18# font size for plots

# Option to plot velocity profiles
vel_slice = True
vel_slice = False # comment this out if desired on
        sdv1 = np.zeros_like(tsr)
        sdv2 = np.zeros_like(tsr)
        sdv3 = np.zeros_like(tsr)
        sdv4 = np.zeros_like(tsr)
        rat1 = np.zeros_like(tsr)
        wdt1 = np.zeros_like(tsr)
        spr1 = np.zeros_like(tsr)
        spr2 = np.zeros_like(tsr)
        spr3 = np.zeros_like(tsr)
        spr4 = np.zeros_like(tsr)
        scl1 = np.zeros_like(tsr)
        scl2 = np.zeros_like(tsr)
        scl3 = np.zeros_like(tsr)

        for i in range(np.size(tsr)):
            men,sdv,rat,wdt,spr,scl = velocity2(tsr[i],solidity[k])
            men1[i] = men[0]
            sdv1[i] = sdv[0]
            sdv2[i] = sdv[1]
            sdv3[i] = sdv[2]
            sdv4[i] = sdv[3]
            rat1[i] = rat[0]
            wdt1[i] = wdt[0]
            spr1[i] = spr[0]
            spr2[i] = spr[1]
            spr3[i] = spr[2]
            spr4[i] = spr[3]
            scl1[i] = scl[0]
            scl2[i] = scl[1]
            scl3[i] = scl[2]