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