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
* GROUP 1. Run identifiers and other preliminaries.
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
TEXT(ON-LINE WIND TURBINES:144
NY=22;NZ=40;NX=1
* GROUP 4. y-direction grid specification.
NY=25
YVLAST=4.700E+02
YFRAC(1)=4.255E-02;YFRAC(2)=6.383E-02
YFRAC(3)=8.511E-02;YFRAC(4)=1.064E-01
YFRAC(5)=1.277E-01;YFRAC(6)=1.383E-01
YFRAC(7)=1.596E-01;YFRAC(8)=1.809E-01
YFRAC(9)=2.234E-01;YFRAC(10)=2.319E-01
YFRAC(11)=2.404E-01;YFRAC(12)=2.489E-01
YFRAC(13)=2.574E-01;YFRAC(14)=2.660E-01
YFRAC(15)=2.819E-01;YFRAC(16)=2.979E-01
YFRAC(17)=3.617E-01;YFRAC(18)=4.415E-01
YFRAC(19)=5.213E-01;YFRAC(20)=6.011E-01
YFRAC(21)=6.809E-01;YFRAC(22)=7.606E-01
YFRAC(23)=8.404E-01;YFRAC(24)=9.202E-01
YFRAC(25)=1.000
* GROUP 5. z-direction grid specification.
NZ=36
ZWLAST=1.300E+03
ZFRAC(1)=3.077E-02;ZFRAC(2)=6.154E-02
ZFRAC(3)=8.462E-02;ZFRAC(4)=1.154E-01
ZFRAC(5)=1.538E-01;ZFRAC(6)=1.923E-01
ZFRAC(7)=2.308E-01;ZFRAC(8)=2.692E-01
ZFRAC(9)=3.077E-01;ZFRAC(10)=3.462E-01
ZFRAC(11)=3.846E-01;ZFRAC(12)=4.231E-01
ZFRAC(13)=4.615E-01;ZFRAC(14)=4.846E-01
ZFRAC(15)=5.077E-01;ZFRAC(16)=5.154E-01
ZFRAC(17)=5.231E-01;ZFRAC(18)=5.385E-01
ZFRAC(19)=5.538E-01;ZFRAC(20)=5.692E-01
ZFRAC(21)=5.846E-01;ZFRAC(22)=6.000E-01
ZFRAC(23)=6.154E-01;ZFRAC(24)=6.308E-01
ZFRAC(25)=6.462E-01;ZFRAC(26)=6.615E-01
ZFRAC(27)=6.769E-01;ZFRAC(28)=6.923E-01
ZFRAC(29)=7.000E-01;ZFRAC(30)=7.308E-01
ZFRAC(31)=7.692E-01;ZFRAC(32)=8.077E-01
ZFRAC(33)=8.462E-01;ZFRAC(34)=8.846E-01
ZFRAC(35)=9.231E-01;ZFRAC(36)=1.000
* GROUP 7. Variables (including porosities) named,
* stored & solved.
SOLVE(W1,V1,P1)
SOLUTN(P1,Y,Y,Y,N,N,N)
STORE(ENUT,PRPS)
TURMOD(KEMODL)
* GROUP 9. Properties of the medium (or media).
RHO1=1.2200E+00
ENUL=1.4650E-05
* GROUP 11. Initialization of fields of variables,
* porosities, etc.
* BLOCKAGE boundary condition, name B1
CONPOR(0,CELL,1,1,-1,-1,-3,-3)
* BLOCKAGE boundary condition, name B2
CONPOR(0,CELL,1,1,-1,-2,-4,-4)
* BLOCKAGE boundary condition, name B21
CONPOR(0,CELL,1,1,-1,-3,-5,-5)
* BLOCKAGE boundary condition, name B3
CONPOR(0,CELL,1,1,-1,-4,-6,-13)
* BLOCKAGE boundary condition, name B4
CONPOR(0,CELL,1,1,-1,-6,-14,-14)
* BLOCKAGE boundary condition, name B5
CONPOR(0,CELL,1,1,-1,-7,-15,-17)
* BLOCKAGE boundary condition, name B6
CONPOR(0,CELL,1,1,-1,-6,-18,-19)
* BLOCKAGE boundary condition, name B7
CONPOR(0,CELL,1,1,-1,-5,-20,-30)
* BLOCKAGE boundary condition, name B8
CONPOR(0,CELL,1,1,-1,-4,-31,-35)
* BLOCKAGE boundary condition, name B9
CONPOR(0,CELL,1,1,-1,-3,-36,-36)
*
FIINIT(W1)=4.0
FIINIT(KE)= 1.0
FIINIT(EP)= 3.0E-04
* GROUP 13. Boundary conditions and special sources
* INLET boundary condition, name INL
PATCH(INL,LOW,1,NX,1,NY,1,1,1,1)
COVAL(INL,W1,0.,4.)
COVAL(INL,P1,FIXFLU,rho1*4.)
COVAL(INL,KE,ONLYMS, 1.0)
COVAL(INL,EP,ONLYMS, .0003)
* Propeller momentum sinks
PATCH(WAKE1,PHASEM,1,1,10,16,16,16,1,1)
COVAL(WAKE1,W1,-0.6,0.0)
PATCH(WAKE2,PHASEM,1,1,9,14,29,29,1,1)
COVAL(WAKE2,W1,-0.6,0.0)
* Free North boundary
PATCH(FREE,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(FREE,P1,fixp,0.0)
COVAL(FREE,W1,0.,SAME)
COVAL(FREE,V1,0.,SAME)
COVAL(FREE,KE,ONLYMS,SAME)
COVAL(FREE,EP,ONLYMS,SAME)
* OUTLET boundary condition, name OUT
PATCH(OUT,HIGH,1,1,4,NY,NZ,NZ,1,1)
COVAL(OUT,W1,0.,SAME)
COVAL(OUT,V1,0.,SAME)
COVAL(OUT,P1,1.E3,0.)
COVAL(OUT,KE,ONLYMS,SAME)
COVAL(OUT,EP,ONLYMS,SAME)
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 15. Termination criteria for sweeps and
* outer iterations.
LSWEEP=250
RELAX(P1,LINRLX,0.4)
real(fac);fac=100.
REAL(SCALEL,SCALEU);SCALEL= 1.000E-01;SCALEU= 3.00
RELAX(W1,FALSDT,fac*SCALEL/SCALEU)
RELAX(V1,FALSDT,fac*SCALEL/SCALEU)
RELAX(KE,FALSDT, fac*SCALEL/SCALEU)
RELAX(EP,FALSDT, fac*SCALEL/SCALEU)
* Assign cell-indicies of spot-point monitoring location
IZMON=23;IYMON=12
* GROUP 23. Variable-by-variable field printout and plot
* and/or tabulation of spot-values and residuals.
ITABL=3;NPLT=1
NAMSAT=MOSG
tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=1.960E+00; EX(V1)=3.132E-01; EX(W1)=2.452E+00
EX(KE)=1.561E+00; EX(EP)=8.638E-03; EX(VPOR)=8.278E-01
EX(PRPS)=3.493E+01; EX(ENUT)=4.731E+01
LIBREF=102
STOP