GROUP 1. Run title and other preliminaries TEXT(RSTM_2D PARABOLIC PLANE WALL JET :T604 TITLE DISPLAY The flow considered is a heated non-buoyant plane submerged jet directed along a plane adiabatic wall into a stagnant environment. For the wall jet, experiments indicate non- coincidence of the lateral poistions of zero shear stress and maximum velocity. This feature can be predicted by a Reynolds- stress transport model (RSTM), but not with an eddy-viscosity model which forces the shear stress to vanish at the location of the velocity maximum. However, for practical purposes this failing does not adversely affect other aspects of the computed flow field. ENDDIS The turbulence is simulated with RSTM, which may be one of three variants, namely: the IPM pressure-strain model (IRSMHM=0); the model coefficients of Gibson & Younis plus IPM (IRSMHM=1); the QIM pressure-strain model (IRSMHM=2). Modelled tranpsort equations are solved for the turbulent heat fluxes. The calculations are started at the jet discharge, and the parabolic marching integration is carried out until both the mean flow and turbulence profiles become self similar. Fifty lateral radial grid cells are employed across the jet, together with a forward step size of about 8% (DZW1) of the local jet width. The y-extent of the grid increases linearly with downstream distance so as to accommodate the spread of the jet. The calculation is terminated about 60 slot widths from the jet discharge. The grid distribution has not been optimised, particularly in the wall-layer region where the grid is relatively coarse and the near-wall grid is too close to the wall. The main results of the calculations are compared with those of experiment in the table below. plane wall jet IPM IPM/GY QIM Data dyw/dz .091 .077 .092 .076 dyt/dz .104 .090 .107 - vw,max/wm**2 .018 .015 .018 .015 The table indicates the values of the half-width spreading rates of the velocity and temperature fields, and the maximum values of the normalised values of the cross-stream turbulent shearing stress and turbulent heat flux. No grid sensitivity studies have been performed, but these results are in close agreement with those reported in the literature. IRSMSM=2;IRSMHM=2;CARTES=T ** Jet-Discharge values REAL(REYNO,TJET,WJET,TKEIN,EPSIN,HSLOT,GMIXL,DYLDZ,DTF) INTEGER(NYS) REYNO=1.25E4;TJET=1.0;HSLOT=.0125;WJET=5.;TKEIN=0.01*WJET*WJET GMIXL=0.09*HSLOT;EPSIN=.1643*TKEIN**1.5/GMIXL GROUP 3-5. Grid specification PARAB=T;NY=50;NYS=44;NREGY=2 IREGY=1;GRDPWR(Y,NYS,HSLOT,1.4) IREGY=2;GRDPWR(Y,6,0.25*HSLOT,1.1) DYLDZ=0.09*2.5;NZ=140;DZW1=0.08;DTF=0.004 DTF=0.001 AZYV=1.0;ZWADD=YVLAST/DYLDZ;AZDZ=PROPY;IPARAB=1 GROUP 7. Variables stored, solved & named SOLVE(P1,V1,W1,H1) PATCH(WAL1,SWALL,1,1,1,1,1,NZ,1,1) TURMOD(REYSTRS,DTF,WAL1) GROUP 8. Terms (in differential equations) & devices DIFCUT=0.0 GROUP 9. Properties of the medium (or media) ENUL=WJET*HSLOT/REYNO GROUP 11. Initialization of variable or porosity fields FIINIT(KE)=TKEIN;FIINIT(EP)=EPSIN FIINIT(U2RS)=2.*TKEIN/3. FIINIT(V2RS)=FIINIT(U2RS);FIINIT(W2RS)=FIINIT(U2RS) FIINIT(VWRS)=0.3*TKEIN;FIINIT(H1)=TJET GROUP 13. Boundary conditions and special sources 1. Outer Boundary -- free stream PATCH(HIGHY,NORTH,1,1,NY,NY,1,NZ,1,1) COVAL(HIGHY,P1,1.0E4,0.0);COVAL(HIGHY,H1,ONLYMS,0.0) COVAL(HIGHY,W1,ONLYMS,0.0);COVAL(HIGHY,V1,ONLYMS,0.0) 2. Inlet Boundary-- uniform flow INLET(UNIFORM,LOW,1,1,1,NYS,1,1,1,1) VALUE(UNIFORM,P1,WJET);VALUE(UNIFORM,W1,WJET) VALUE(UNIFORM,H1,TJET) VALUE(UNIFORM,KE,TKEIN);VALUE(UNIFORM,EP,EPSIN) VALUE(UNIFORM,W2RS,2.*TKEIN/3.);VALUE(UNIFORM,V2RS,2.*TKEIN/3.) VALUE(UNIFORM,U2RS,2.*TKEIN/3.) GROUP 16. Termination of iterations LITHYD=30 GROUP 18. Limits on variables or increments to them VARMIN(W1)=1.E-10;VARMIN(EP)=1.E-8*RHO1*HSLOT*WJET*EPSIN VARMIN(U2RS)=1.E-10;VARMIN(V2RS)=1.E-10;VARMIN(W2RS)=1.E-10 VARMIN(H1)=1.E-10 GROUP 22. Monitor print-out IZMON=NZ/2;IYMON=NY/2;ITABL=3;NPLT=5;IPLTL=LITHYD TSTSWP=-1;NYPRIN=1;NZPRIN=NZ GROUP 23. Field print-out and plot control ORSIZ=0.4;PATCH(IZEQNZ,PROFIL,1,1,1,NY,NZ,NZ,1,1) PLOT(IZEQNZ,W1,0.0,0.0) GROUP 24. Dumps for restarts RESREF(P1)=1.E-8*HSLOT*WJET; RESREF(W1)=1.E-8*RHO1*HSLOT*WJET*WJET RESREF(H1)=1.E-8*RHO1*HSLOT*TJET*WJET; RESREF(V1)=RESREF(W1) RG(1)=HSLOT;LG(1)=T;IG(1)=14