TALK=T;RUN(1,1)
PHOTON USE
p
GR OU Y 2 X 1 1 Z 3 3 COL 2
GR OU Y 1 X 1 1 Z 3 3 COL 2
GR OU Y 3 X 1 1 Z 4 4 COL 3
GR OU Y 1 X 1 1 Z 4 4 COL 3
GR OU Y 4 X 1 1 Z 5 5 COL 4
GR OU Y 1 X 1 1 Z 5 5 COL 4
GR OU Y 5 X 1 1 Z 6 13 COL 5
GR OU Y 1 X 1 1 Z 6 13 COL 5
GR OU Y 7 X 1 1 Z 14 14 COL 6
GR OU Y 1 X 1 1 Z 14 14 COL 6
GR OU Y 8 X 1 1 Z 15 17 COL 7
GR OU Y 1 X 1 1 Z 15 17 COL 7
GR OU Y 7 X 1 1 Z 18 19 COL 8
GR OU Y 1 X 1 1 Z 18 19 COL 8
GR OU Y 6 X 1 1 Z 20 30 COL 9
GR OU Y 1 X 1 1 Z 20 30 COL 9
GR OU Y 5 X 1 1 Z 31 35 COL 10
GR OU Y 1 X 1 1 Z 31 35 COL 10
GR OU Y 4 X 1 1 Z 36 36 COL 11
GR OU Y 1 X 1 1 Z 36 36 COL 11
GR OU X 2 Y 10 16 Z 16 16 COL 5
GR OU X 1 Y 10 16 Z 16 16 COL 5
GR OU Y 17 X 1 1 Z 16 16 COL 5
GR OU Y 10 X 1 1 Z 16 16 COL 5
GR OU X 2 Y 9 14 Z 29 29 COL 6
GR OU X 1 Y 9 14 Z 29 29 COL 6
GR OU Y 15 X 1 1 Z 29 29 COL 6
GR OU Y 9 X 1 1 Z 29 29 COL 6
GR OU X 1; GR OU X M; GR OU Y 1; GR OU Y M
red
msg ON-LINE WIND TURBINES
msg Geometry
msg Press Enter to continue
pause
gr OFF;red
msg ON-LINE WIND TURBINES
MSG Grid
GR X 1
msg Press Enter to continue
pause
gr OFF;red
gr ou x 1
msg ON-LINE WIND TURBINES
msg Velocity field
VE X 1 SH
msg Press Enter to continue
pause
ve OFF;red
msg ON-LINE WIND TURBINES
msg Streamwise velocity contours
con w1 x 1 sh;int 10
msg Press Enter to continue
pause
con OFF;red
msg ON-LINE WIND TURBINES
msg Turbulence energy contours
con ke x 1 fil;0.01
msg Press E to end
ENDUSE
DISPLAY
This run introduces the influence of the stand of widely-spaced
trees on the wake fields generated by two on-line wind turbines
located on the complex terrain of Valaam island of Ladoga lake
in Nothern Russia. The standard K-E model is modified to take
into account the plant airflow interactions resulting in
momentum sinks. The corresponding GROUND codings are made by way
of PLANT option. The main details can be found in paper by
S.Zhubrin et al."Practical evalution of turbulence models
applied to on-line wind turbines on a complex terrain",
presented at International Seminar of AirFlows on Complex
Terrain, University of Karlsrue, February 1995
ENDDIS
PLANTBEGIN
* Plant airflow interactions
- form darg
PATCH(INTER,PHASEM,1,1,5,9,7,12,1,1)
CO=0.5*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,V1,GRND,0.0)
CO=0.5*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,W1,GRND,0.0)
- wake production
VAL=0.5*0.6*10.5*(V1*V1+W1*W1)**1.5
COVAL(INTER,KE,fixflu,GRND)
- turbulence energy loss
CO=2.*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,KE,GRND,0.0)
- dissipation production
VAL=0.75*EP/KE*0.6*10.5*(V1*V1+W1*W1)**1.5
COVAL(INTER,EP,fixflu,GRND)
- dissipation loss
CO=3.*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,EP,GRND,0.0)
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(ON-LINE WIND TURBINES:144 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 102
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =1.
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 25
YVLAST =470.
YFRAC(1)=0.04255 ;YFRAC(2)=0.06383
YFRAC(3)=0.08511 ;YFRAC(4)=0.1064
YFRAC(5)=0.1277 ;YFRAC(6)=0.1383
YFRAC(7)=0.1596 ;YFRAC(8)=0.1809
YFRAC(9)=0.2234 ;YFRAC(10)=0.2319
YFRAC(11)=0.2404 ;YFRAC(12)=0.2489
YFRAC(13)=0.2574 ;YFRAC(14)=0.266
YFRAC(15)=0.2819 ;YFRAC(16)=0.2979
YFRAC(17)=0.3617 ;YFRAC(18)=0.4415
YFRAC(19)=0.5213 ;YFRAC(20)=0.6011
YFRAC(21)=0.6809 ;YFRAC(22)=0.7606
YFRAC(23)=0.8404 ;YFRAC(24)=0.9202
YFRAC(25)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 36
ZWLAST =1300.
ZFRAC(1)=0.03077 ;ZFRAC(2)=0.06154
ZFRAC(3)=0.08462 ;ZFRAC(4)=0.1154
ZFRAC(5)=0.1538 ;ZFRAC(6)=0.1923
ZFRAC(7)=0.2308 ;ZFRAC(8)=0.2692
ZFRAC(9)=0.3077 ;ZFRAC(10)=0.3462
ZFRAC(11)=0.3846 ;ZFRAC(12)=0.4231
ZFRAC(13)=0.4615 ;ZFRAC(14)=0.4846
ZFRAC(15)=0.5077 ;ZFRAC(16)=0.5154
ZFRAC(17)=0.5231 ;ZFRAC(18)=0.5385
ZFRAC(19)=0.5538 ;ZFRAC(20)=0.5692
ZFRAC(21)=0.5846 ;ZFRAC(22)=0.6
ZFRAC(23)=0.6154 ;ZFRAC(24)=0.6308
ZFRAC(25)=0.6462 ;ZFRAC(26)=0.6615
ZFRAC(27)=0.6769 ;ZFRAC(28)=0.6923
ZFRAC(29)=0.7 ;ZFRAC(30)=0.7308
ZFRAC(31)=0.7692 ;ZFRAC(32)=0.8077
ZFRAC(33)=0.8462 ;ZFRAC(34)=0.8846
ZFRAC(35)=0.9231 ;ZFRAC(36)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(149)=PRPS
NAME(150)=ENUT
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
VIST = 150
PRPS = 149
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1.22 ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.465E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=4. ;FIINIT(KE)=1.
FIINIT(EP)=3.0E-04 ;FIINIT(PRPS)=-1.
FIINIT(ENUT)=1.0E-10
PATCH(CMP1 ,INIVAL, 1, 1, 1, 1, 3, 3, 1, 1)
INIT(CMP1 ,PRPS,0. ,198. )
PATCH(CMP2 ,INIVAL, 1, 1, 1, 2, 4, 4, 1, 1)
INIT(CMP2 ,PRPS,0. ,198. )
PATCH(CMP3 ,INIVAL, 1, 1, 1, 3, 5, 5, 1, 1)
INIT(CMP3 ,PRPS,0. ,198. )
PATCH(CMP4 ,INIVAL, 1, 1, 1, 4, 6, 13, 1, 1)
INIT(CMP4 ,PRPS,0. ,198. )
PATCH(CMP5 ,INIVAL, 1, 1, 1, 6, 14, 14, 1, 1)
INIT(CMP5 ,PRPS,0. ,198. )
PATCH(CMP6 ,INIVAL, 1, 1, 1, 7, 15, 17, 1, 1)
INIT(CMP6 ,PRPS,0. ,198. )
PATCH(CMP7 ,INIVAL, 1, 1, 1, 6, 18, 19, 1, 1)
INIT(CMP7 ,PRPS,0. ,198. )
PATCH(CMP8 ,INIVAL, 1, 1, 1, 5, 20, 30, 1, 1)
INIT(CMP8 ,PRPS,0. ,198. )
PATCH(CMP9 ,INIVAL, 1, 1, 1, 4, 31, 35, 1, 1)
INIT(CMP9 ,PRPS,0. ,198. )
PATCH(CMP10 ,INIVAL, 1, 1, 1, 3, 36, 36, 1, 1)
INIT(CMP10 ,PRPS,0. ,198. )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(INL ,LOW , 1, 1, 1, 25, 1, 1, 1, 1)
COVAL(INL ,P1 , FIXFLU ,4.88 )
COVAL(INL ,W1 ,0. ,4. )
COVAL(INL ,KE ,0. ,1. )
COVAL(INL ,EP ,0. ,3.0E-04 )
PATCH(WAKE1 ,PHASEM, 1, 1, 10, 16, 16, 16, 1, 1)
COVAL(WAKE1 ,W1 ,-0.6 ,0. )
PATCH(WAKE2 ,PHASEM, 1, 1, 9, 14, 29, 29, 1, 1)
COVAL(WAKE2 ,W1 ,-0.6 ,0. )
PATCH(FREE ,NORTH , 1, 1, 25, 25, 1, 36, 1, 1)
COVAL(FREE ,P1 ,1. ,0. )
COVAL(FREE ,V1 ,0. , SAME )
COVAL(FREE ,W1 ,0. , SAME )
COVAL(FREE ,KE ,0. , SAME )
COVAL(FREE ,EP ,0. , SAME )
PATCH(OUT ,HIGH , 1, 1, 4, 25, 36, 36, 1, 1)
COVAL(OUT ,P1 ,1000. ,0. )
COVAL(OUT ,V1 ,0. , SAME )
COVAL(OUT ,W1 ,0. , SAME )
COVAL(OUT ,KE ,0. , SAME )
COVAL(OUT ,EP ,0. , SAME )
PATCH(INTER ,PHASEM, 1, 1, 5, 9, 7, 12, 1, 1)
COVAL(INTER ,V1 , GRND ,0. )
COVAL(INTER ,W1 , GRND ,0. )
COVAL(INTER ,KE , GRND ,0. )
COVAL(INTER ,EP , GRND ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 250 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.4)
RELAX(V1,FALSDT,3.333333)
RELAX(W1,FALSDT,3.333333)
RELAX(KE,FALSDT,3.333333)
RELAX(EP,FALSDT,3.333333)
RELAX(PRPS,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
KELIN = 0
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=1.96 ;EX(V1)=0.3132
EX(W1)=2.452 ;EX(KE)=1.561
EX(EP)=8.638E-03 ;EX(PRPS)=34.93
EX(ENUT)=47.310001
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 12 ;IZMON = 23
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP