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