L(190)
GROUP 1. Run titlee and other preliminaries
TEXT(Boundary Layer K-L Turbulen Model
TITLE
The flow of case 190 is considered; but here the Prandtl one-
equation model is used, the mixing-length Lm being
specified according to Escudier's formula, ie: Lm=k*y for
y/d<0.09/k, and Lm=0.09*d for y/d>0.09. Here k is the von
Karman's constant and y is the normal distance from the wall.
GROUP 7. Variables stored, solved & named
TURMOD(KLMODL);GENK=F
GROUP 9. Properties of the medium (or media)
** Select Prandtl-Kolmogorov Formula : Group 9/Sect. 5 of GREX3
EL1=MIXLENBL;EL1A=0.0;EL1B=AK
GROUP 11. Initialization of variable or porosity fields
FIINIT(KE)=WFREE**2*0.01
GROUP 13. Boundary conditions and special sources
** South Wall Boundary
COVAL(WFUN,KE,LOGLAW,LOGLAW)
** North Free Boundary
COVAL(FREE,KE,ONLYMS,1.E-10)
** Inlet Boundary
COVAL(PROF,KE,ONLYMS,GRND3)
GROUP 21. Print-out of variables
OUTPUT(KE,Y,Y,Y,Y,Y,Y)
GROUP 23. Field print-out and plot control
PLOT(IZEQNZ,KE,0.,0.0)