TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2017, Date 20/01/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 standard k-f model and fully-rough logarithmic
wall functions. In this region the velocity and turbulence
variables are given by:
u/u*=ln(z/z0)/K ; k=u**2/sqrt(cmucd)
f=u*/(sqrt(cmucd)*K*z); enut=K*z*u*
where u*=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 u*=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 a uniform pressure field.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 312
************************************************************
Group 1. Run Title
TEXT(NEUTRAL ABL WIND FLOW - K-F MODEL: T312 )
************************************************************
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(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)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
TURMOD(KOMODL)
************************************************************
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
No PATCHes used for this Group
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(VABS)=5.768
EX(STRS)=3.855E-03 ;EX(OMEG)=0.3141
EX(DEN1)=1.189 ;EX(EL1)=15.64
EX(ENUT)=5.521
************************************************************
Group 21. Print-out of Variables
************************************************************
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(DEPTH,3.0E+04)
GVIEW(VDIS,110.552322)
GVIEW(CENTRE,100.,0.5,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