FFarmTools/FFarmInputSetup.py at master · kshaler/FFarmTools · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
import os, sys
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
from stochasticTurbulenceTools_mod import stochasticTurbulence
from FFarmCaseCreation import FFarmCaseCreation, WriteFFarmFile

plt.rc("font",size=14)

ParamsFile = './FFarmParams.txt'

Pars=pd.read_csv(ParamsFile, delimiter='\t',header=None,skiprows=1,index_col=0,nrows=5,names=['Var','Val'])
xlocs_tmp=pd.read_csv(ParamsFile, delimiter='\t',header=None,skiprows=6,index_col=0,nrows=1)
ylocs_tmp=pd.read_csv(ParamsFile, delimiter='\t',header=None,skiprows=7,index_col=0,nrows=1)
zlocs_tmp=pd.read_csv(ParamsFile, delimiter='\t',header=None,skiprows=8,index_col=0,nrows=1)

xlocs=[xlocs_tmp[1][0], xlocs_tmp[2][0], xlocs_tmp[3][0], xlocs_tmp[4][0], xlocs_tmp[5][0]]
ylocs=[ylocs_tmp[1][0], ylocs_tmp[2][0], ylocs_tmp[3][0], ylocs_tmp[4][0], ylocs_tmp[5][0]]
zlocs=[zlocs_tmp[1][0], zlocs_tmp[2][0], zlocs_tmp[3][0], zlocs_tmp[4][0], zlocs_tmp[5][0]]

D=Pars['Val']['D']
HubHt=Pars['Val']['HubHt']
high_extent_X=Pars['Val']['high_extent_X']
high_extent_Y=Pars['Val']['high_extent_Y']
nTurbs=int(Pars['Val']['nTurbs'])

Case = FFarmCaseCreation(nTurbs,'Neutral_TS1')
tpath = '../turbineModel/Test18' ## do not include .fst extention
Case.Turb(D,HubHt,tpath,cmax=3)

## If ReadTS=True, TurbSim parameters will be read directly from the TurbSim.inp or TurbSim.bts file
ReadTS = True

## if ReadINP = True, TurbSim parameters will be read directly from the TurbSim.inp file. Otherwise, TurbSim.bts file will be used.
ReadINP = False

if ReadTS == True:
    TSpath = '.'
    if ReadINP == True:
        TSFile = os.path.join(TSpath+os.sep+Case.prefix+'.inp')

    else:
        TSFile = os.path.join(TSpath+os.sep+Case.prefix+'.bts')
        IFdata = stochasticTurbulence(Case.D,prefix=Case.prefix)
        IFdata.readBTS(TSpath, HubHt)

        IFdata.kHub = IFdata.z2k(HubHt)

        Vhub = IFdata.u[:,IFdata.jHub,IFdata.kHub].mean()
else:
    print('TurbSim parameters must be entered directly.')

Case.turbLocs(xlocs,ylocs,zlocs)
Case.discretization(Vhub)
Case.highResDomain(IFdata,Vhub,Case.dX_High_desired,high_extent_X,high_extent_Y)
Case.lowResDomain(Case.dX_Low_desired,Vhub)

FFpath = './'
OrigFile = os.path.join(FFpath+Case.prefix+'.fstf')
NewFile = os.path.join(FFpath+Case.prefix+'_mod.fstf')
WriteFFarmFile(OrigFile, NewFile, Case, NewFile=False)


fig = plt.figure(figsize=(13.5,10))
ax  = fig.add_subplot(111,aspect="equal")

xmax_low = Case.X0_Low+Case.dX_Low*Case.nX_Low
ymax_low = Case.Y0_Low+Case.dY_Low*Case.nY_Low

# low-res box
ax.plot([Case.X0_Low,xmax_low,xmax_low,Case.X0_Low,Case.X0_Low],
        [Case.Y0_Low,Case.Y0_Low,ymax_low,ymax_low,Case.Y0_Low],'--k',lw=2,label='Low')

# high-res boxes
for wt in range(Case.nTurbs):
    xmax_high = Case.X0_High[wt]+Case.dX_High*Case.nX_High
    ymax_high = Case.Y0_High[wt]+Case.dY_High*Case.nY_High
    ax.plot([Case.X0_High[wt],xmax_high,xmax_high,Case.X0_High[wt],Case.X0_High[wt]],
            [Case.Y0_High[wt],Case.Y0_High[wt],ymax_high,ymax_high,Case.Y0_High[wt]],
            '-',
            label="HighT{0}".format(wt+1))
    ax.plot(Case.x[wt],Case.y[wt],'x',ms=8,mew=2,label="WT{0}".format(wt+1))

plt.legend(bbox_to_anchor=(1.05,1.015),frameon=False)

ax.set_xlabel("x-location [m]")
ax.set_ylabel("y-location [m]")

fig.tight_layout
fig.savefig('FFarmLayout.pdf',bbox_to_inches='tight',dpi=500)