TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2017, Date 01/09/17
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
The case considered is the simulation of the constant-stress-
layer region of a 2d turbulent neutral atmospheric boundary
layer using the Wilcox 1988 k-w model and fully-rough logarithmic
wall functions. In this region the velocity and turbulence
variables are given by:
u/utau=ln(z/z0)/K ; k=utau^2/sqrt(cmucd)
w=utau/(sqrt(cmucd)*K*z); enut=K*z*utau
where utau=friction velocity, K=0.41, cmucd=0.09 and z0 is the
roughness height. A calculation is made using the following
wind inlet-profile data set by the wind object: Wind reference
speed: Ur=5m/s at a reference height Zr=10m, wind direction
from the west, and terrain roughness height z0=0.03 m. The
friction velocity is given by utau=Ur*K/ln(zr/z0)=0.353m/s
and so the shear_stress/density=0.125 m^2/s^2. The wind object
is used which sets the wind inlet profiles, a fixed-pressure
condition at the downstream boundary, and a constant shear-
stress condition at the sky boundary via an entrainment and
diffusive boundary condition. The height(z) of the solution
domain is 100m, and the domain length is 200m(x). The wind enters
from the west boundary, and leaves via the east boundary. The
flow is 2d, although the simulation is 3d. The calculation
demonstrates that PHOENICS preserves the wind inlet profiles
and returns what is an essentially uniform pressure field.
ENDDIS
save1begin
The inlet velocity and turbulence log-law profiles aren't
preserved exactly because a spike in the k profile is produced
at the second cell from the wall, which originates from the
standard discretisation practices employed for the turbulence
production rate Pk in the turbulence transport equations and the
vertical diffusion terms in the axial momentum equation (see
library case T306 for more details).
save1end
************************************************************
IRUNN = 1 ;LIBREF = 703
************************************************************
Group 1. Run Title
TEXT(T703:2DZY ABL WIND FLOW- Wilcox 1988 k-w)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,30,5,30)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(137)=EVDS ;NAME(138)=KDIS
NAME(139)=UDIS ;NAME(140)=EVIN
NAME(141)=OMIN ;NAME(142)=EPIN
NAME(143)=KEIN ;NAME(144)=UIN
NAME(145)=VABS ;NAME(146)=STRS
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1,U1,V1,W1)
* Stored variables list
STORE(ENUT,EL1,DEN1,STRS,VABS,UIN,KEIN,EPIN)
STORE(OMIN,EVIN,UDIS,KDIS,EVDS)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
TURMOD(KWMODL)
************************************************************
Echo save-block settings for Group 7
save7begin
REAL(VKC,XLEN,REFH,zdz0)
REAL(USTAR,USTAR2,Z0,UIN)
REAL(KCNS,USTAR3,EPCON)
UIN =5.0 ! reference velocity
REFH=10. ! reference height
Z0=0.03 ! roughness height
VKC=0.41 ! von karman's constant
ZDZ0=REFH/Z0
USTAR=UIN*VKC/LOG(zdz0) ! friction velocity
USTAR
USTAR2=USTAR*USTAR;KCNS=USTAR2/0.3;USTAR3=USTAR2*USTAR
KCNS
EPCON=USTAR3/VKC
*** store fully-developed boundary-layer values
(stored of UIN is (:ustar:*loge(ZG/:Z0:)/:vkc:)!zslstr)
(stored of KEIN is :kcns:!zslstr)
(stored of EPIN is (:epcon:/ZG)!zslstr)
(stored of EVIN is :vkc:*:ustar:*ZG!zslstr)
(stored of OMIN is :epcon:/(ZG*0.09*:kcns:)!zslstr)
** discrepancy with log-law values
(stored of UDIS is 100.*(U1/UIN-1.0) with IF(IX.LT.NX))
(stored of KDIS is 100.*(KE/KEIN-1.0))
(stored of EVDS is 100.*(ENUT/EVIN-1.0))
save7end
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
DVO1DT =3.41E-03
PRT(KE)=2. ;PRT(OMEG)=2.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0.
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
PATCH(KWSOURCE, PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KWSOURCE, KE, GRND4, GRND4)
COVAL(KWSOURCE, OMEG, GRND4, GRND4)
BUOYD =1.189
BUOYE =20.
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(KE ,FALSDT,1. )
RELAX(OMEG,FALSDT,1. )
************************************************************
Group 18. Limits
VARMAX(OMEG)=1.0E+10 ;VARMIN(OMEG)=1.0E-10
************************************************************
Group 19. EARTH Calls To GROUND Station
GENK = T
PARSOL = F
IENUTA = 10
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MAXINC,KE,R,0.025)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=2.405E-04 ;EX(U1)=5.768
EX(V1)=2.113E-04 ;EX(W1)=1.426E-04
EX(KE)=0.4166 ;EX(EP)=0.0118
EX(VABS)=5.768 ;EX(STRS)=3.855E-03
EX(OMEG)=0.3141 ;EX(DEN1)=1.189
EX(EL1)=15.64 ;EX(ENUT)=5.521
EX(EVDS)=1.031E+00
EX(KDIS)=1.247E+00
EX(UDIS)=2.055E-01
EX(EVIN)=5.531E+00
EX(OMIN)=3.055E-01
EX(EPIN)=1.141E-02
EX(KEIN)=4.151E-01
EX(UIN )=5.767E+00
************************************************************
Group 21. Print-out of Variables
OUTPUT(EP ,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 23 ;IYMON = 1 ;IZMON = 5
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.,-1.,0.)
GVIEW(UP,0.,0.,1.)
GVIEW(VDIS,109.981827)
GVIEW(CENTRE,100.,25.,50.)
> DOM, SIZE, 2.000000E+02, 5.000000E+01, 1.000000E+02
> DOM, MONIT, 1.500000E+02, 5.000000E+00, 7.402118E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, RSET_X_1, 30, 1.000000E+00,G
> GRID, RSET_Y_1, 5, 1.000000E+00
> GRID, RSET_Z_1, 30, 1.050000E+00,G
> DOM, INI_AMB, YES
> DOM, INI_BUOY, YES
> OBJ, NAME, WIND1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, windgr
> OBJ, ROTATION24, 7
> OBJ, VISIBLE, NO
> OBJ, TYPE, WIND
> OBJ, WIREFRAME, YES
> OBJ, PRESSURE, 1.01325E+05
> OBJ, COEFFICIENT, 1.0E+04
> OBJ, VELOCITY, 5.
> OBJ, WIND_DIR, 270.
> OBJ, AXIS_DIR, 0.
> OBJ, PROFILE, Logarithmic
> OBJ, REF_HEIGHT, 10.
> OBJ, RGHNS_HEIGHT,0.03
> OBJ, UP-DIR, Z
> OBJ, SKY, YES
> OBJ, GROUND, YES
> OBJ, GROUND-TEMP, ADIABATIC
STOP