GROUP 1. Run title and other preliminaries
TEXT(Temporal Decay Of Turbu; K-E Model
TITLE
  DISPLAY
  The case considered is the temporal decay of turbulence in a
  turbulent parcel of fluid. This case provides a good basic test
  of a two-equation turbulence model because the turbulence
  transport equations may be reduced to  two ordinary differential
  equations to which an analytical solution exists. The solution
  is
 
  k=(k0**a + a*e0*(t0-t)/k0**(1-a)) ** (1/a)
 
  where k0 is the k value at t0 (the value of the constant a that
  fits the experimental data is 1.0-1.92 ie. -0.92) and e0 is the
  dissipation rate at t0.
  ENDDIS
REAL(TLEN,TKEIN,EPIN)
LSTEP=20;TLEN=0.05
 
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,LSTEP,TLEN,1.0)
 
    GROUP 7. Variables stored, solved & named
TURMOD(KEMODL);STORE(ENUT)
 
    GROUP 9. Properties of the medium (or media)
RHO1=1.0;ENUL=1.E-5
 
    GROUP 11. Initialization of variable or porosity fields
TKEIN=1.0;FIINIT(KE)=TKEIN;EPIN=7.5;FIINIT(EP)=EPIN
 
    GROUP 16. Termination of iterations
LSWEEP=4;RESREF(KE)=1.E-5;RESREF(EP)=1.E-5
 
    GROUP 21. Print-out of variables
OUTPUT(KE,Y,Y,Y,N,Y,Y);OUTPUT(EP,Y,Y,Y,N,Y,Y)
 
    GROUP 22. Monitor print-out
NPLT=1;ITABL=2;IPLTL=40
SPEDAT(SET,GXMONI,TRANSIENT,L,F)  
    GROUP 23. Field print-out and plot control
PATCH(TIMEPLOT,PROFIL,1,1,1,1,1,1,1,LSTEP)
PLOT(TIMEPLOT,KE,0.2,1.0);PLOT(TIMEPLOT,EP,0.,10.)
NXPRIN=1;NPRINT=10;NTPRIN=10