GROUP 1. Run title and other preliminaries
TEXT(LAM-BRE_2D Y-Z CONJ CHANNEL TEST   :T217
TITLE
  DISPLAY
  This case provides a test of the Lam-Bremhorst low-Reynolds number
  k-e model for 2d steady developing flow in a channel. As a test of
  symmetry both halves of the channel are considered in the y-z
  plane. In addition the temperature equation is solved via the
  TEM1 variable with conduction through outer slabs of fireclay
  brick on either side of the channel. At each wall a uniformly-
  distributed heat flux is applied to the fluid. The case provides
  a test of the correct representation of blockages via the
  "conjugate-heat-transfer" facility. The results of this
  calculation may be compared with calculations of the same case
  ( excepting conduction in the solid ) without blockages and also
  with blockages represented by the "porosity" facility. A
  relatively coarse grid is used so as to permit convergence in
  about 400 sweeps. The problem has been set up to perform 30
  sweeps only, so as to reduce execution times in performing
  routine testing procedures.
  ENDDIS
 
REAL(HGHT,REY,WIN,DHYDR,REYH,FRIC,US,TIN,HIN);TIN=300.
HGHT=0.1; REY=1.E5;WIN=10.0;DHYDR=2.*HGHT; REYH=2.*REY
FRIC=1./(1.82*LOG10(REYH)-1.64)**2
US=WIN*(FRIC/8.)**0.5;ENULA=WIN*HGHT/REY;HIN=1.E3*TIN
    GROUP 3. X-DIRECTION GRID SPECIFICATION
    GROUP 4. Y-DIRECTION GRID SPECIFICATION
INTEGER(NYG,NYS,NYGF,NYGL,NYG2,NYGH,JJM,JJJ)
REAL(LENGTH,DELT1,DELY,KFAC,AA,BB)
NYS=4;NYGF=NYS+1
  ** define first dely from wall and the grid-expansion
     factor Kfac which defines a constant ratio of lengths of
     two adjacent cells.
DELT1=0.5*ENULA/US;KFAC=1.8;DELY=DELT1/(0.5*HGHT)
  ** calculate NYG from dely & Kfac
AA=(0.5/DELY)*(KFAC-1.0)+1.0;AA=LOG(AA)/LOG(KFAC)+1.0001
NYG2=AA;NYG=2*NYG2;NYGH=NYS+NYG2;NYGL=NYS+NYG;NY=2*NYS+NYG
  ** define uniform grid initially
IREGY=1;GRDPWR(Y,NY,YVLAST,1.0)
  ** define grid in the bottom solid
AA=NYS
DO JJ=1,NYS
+ YFRAC(JJ)=JJ*HGHT/NYS
ENDDO
  ** compute expanding grid from south boundary over one
     half of the channel width
ARRAY(ARR1,REAL,NYG)
ARR1(1)=DELY
DO JJ=2,NYG2
+ JJM=JJ-1
+ DELY=KFAC*DELY
+ ARR1(JJ)=ARR1(JJM)+DELY
+ BB=ARR1(JJ)
ENDDO
ARR1(NYG2)=0.5
  ** create symmetrical grid in the second half of the channel
JJJ=0
DO JJ=NYG-1,NYG2+1,-1
+ JJJ=JJJ+1
+ ARR1(JJ)=1.-ARR1(JJJ)
ENDDO
ARR1(NYG)=1.0
JJJ=0
DO JJ=NYGF,NYGL
+ JJJ=JJJ+1
+ YFRAC(JJ)=ARR1(JJJ)*HGHT+YFRAC(NYS)
ENDDO
JJJ=0
DO JJ=NYGL+1,NY
+ JJJ=JJJ+1
+ YFRAC(JJ)=YFRAC(NYGL)+JJJ*HGHT/NYS
ENDDO
    GROUP 5. Z-direction grid specification
LENGTH=5.*HGHT;NZ=10;GRDPWR(Z,NZ,LENGTH,1.0)
    GROUP 7. Variables stored, solved & named
NAME(49)=WDIS;NAME(50)=ENUT;NAME(40)=BLOK
NAME(39)=TEM1;NAME(38)=PRPS;STORE(WDIS)
SOLVE(P1,V1,W1,C1,TEM1);SOLUTN(P1,Y,Y,Y,P,P,P)
  ** activate whole-field solution for TEM1 & C1
SOLUTN(TEM1,Y,Y,Y,P,P,P);SOLUTN(C1,Y,Y,Y,P,P,P)
STORE(PRPS,BLOK);STOVAR(TEM1,IVARBK);ISOLBK=1
TURMOD(KEMODL-LOWRE);STORE(ENUT);KELIN=1
    GROUP 8. TERMS (IN DIFFERENTIAL EQUATIONS) & DEVICES
    GROUP 9. PROPERTIES OF THE MEDIUM (OR MEDIA)
;ENUL=FILE;PRNDTL(TEM1)=CONDFILE;
    GROUP 11. INITIALIZATION OF VARIABLE OR POROSITY FIELDS
FIINIT(W1)=WIN;FIINIT(TEM1)=TIN;FIINIT(P1)=1.E-3
REAL(TKEIN,EPIN);TKEIN=2.*US*US
EPIN=0.1643*TKEIN**1.5/(0.1*0.5*HGHT)
FIINIT(KE)=TKEIN;FIINIT(EP)=EPIN
    material 34 is air
INIADD=F;FIINIT(PRPS)=34.0;FIINIT(BLOK)=1.0
  ** remove -ves from patch limits as now unable to overwrite
     wall-function covals; define them separately
    material 152 is fireclay brick
CONPOR(SOLIDS,-1.0,CELL,1,NX,1,NYS,1,NZ)
INIT(SOLIDS,PRPS,0.0,152.0);COVAL(SOLIDS,BLOK,0.0,2.0)
 
CONPOR(SOLIDN,-1.0,CELL,1,NX,(NYS+NYG+1),NY,1,NZ)
INIT(SOLIDN,PRPS,0.0,152.0);COVAL(SOLIDN,BLOK,0.0,3.0)
    GROUP 13. BOUNDARY CONDITIONS AND SPECIAL SOURCES
  ** air inflow boundary
PATCH(INL,LOW,1,1,NYS+1,NYS+NYG,1,1,1,1)
COVAL(INL,P1,FIXFLU,WIN);COVAL(INL,TEM1,ONLYMS,TIN)
COVAL(INL,W1,ONLYMS,WIN);COVAL(INL,KE,ONLYMS,TKEIN)
COVAL(INL,EP,ONLYMS,EPIN);COVAL(INL,C1,ONLYMS,1.0)
  ** outflow boundary
PATCH(OUTL,HIGH,1,1,NYS+1,NYS+NYG,NZ,NZ,1,1)
COVAL(OUTL,P1,1.E3,0.0);COVAL(OUTL,TEM1,ONLYMS,SAME)
COVAL(OUTL,W1,ONLYMS,0.0);COVAL(OUTL,KE,ONLYMS,0.0)
COVAL(OUTL,EP,ONLYMS,0.0);COVAL(OUTL,C1,ONLYMS,SAME)
  ** cannot use solids minus feature as patches are
     set up in EARTH & so cannot overwrite covals
  ** north wall
PATCH(SOLID-NW,NWALL,1,NX,NYGL,NYGL,1,NZ,1,1)
COVAL(SOLID-NW,W1,LOGLAW,SAME);COVAL(SOLID-NW,KE,1.0,ZERO)
  ** south wall
PATCH(SOLID-SW,SWALL,1,NX,NYGF,NYGF,1,NZ,1,1)
COVAL(SOLID-SW,W1,LOGLAW,SAME);COVAL(SOLID-SW,KE,1.0,ZERO)
  ** the "conjugate-heat-transfer" facility defaults to a
     finite diffusive link at the boundary, so zero the
     unwanted links at the fluid-solid interface
PATCH(GP12DFNT,CELL,1,NX,NYGL,NYGL,1,NZ,1,1)
COVAL(GP12DFNT,KE,0.0,0.0);COVAL(GP12DFNT,EP,0.0,0.0)
COVAL(GP12DFNT,C1,0.0,0.0)
PATCH(GP12DFNB,CELL,1,NX,NYS,NYS,1,NZ,1,1)
COVAL(GP12DFNB,KE,0.0,0.0);COVAL(GP12DFNB,EP,0.0,0.0)
COVAL(GP12DFNB,C1,0.0,0.0)
  ** the "conjugate-heat-transfer" facility uses an old treatment
     for lateral momentum convection for momentum control volumes
     adjacent to solids; zero links for consistency with porosity
     treatment.
PATCH(GP12DFN1,CELL,1,NX,NYGL-1,NYGL-1,1,NZ,1,1)
COVAL(GP12DFN1,V1,0.0,0.0)
PATCH(GP12DFN2,CELL,1,NX,NYS,NYS,1,NZ,1,1)
COVAL(GP12DFN2,V1,0.0,0.0)
PATCH(GP12CNN1,CELL,1,NX,NYGL-1,NYGL-1,1,NZ,1,1)
COVAL(GP12CNN1,V1,0.0,0.0)
PATCH(GP12CNN2,CELL,1,NX,NYS,NYS,1,NZ,1,1)
COVAL(GP12CNN2,V1,0.0,0.0)
  ** applied heat flux from outer slab to the fluid
PATCH(HEATN,NORTH,1,NX,NYGL,NYGL,1,NZ,1,1)
COVAL(HEATN,TEM1,FIXFLU,0.5E3/LENGTH)
PATCH(HEATS,SOUTH,1,NX,NYS+1,NYS+1,1,NZ,1,1)
COVAL(HEATS,TEM1,FIXFLU,0.5E3/LENGTH)
    GROUP 15. TERMINATION OF SWEEPS
RESREF(P1)=1.E-12*WIN*HGHT
RESREF(W1)=RESREF(P1)*WIN; RESREF(V1)=RESREF(W1)
RESREF(TEM1)=RESREF(P1)*HIN; RESREF(C1)=RESREF(P1)
RESREF(KE)=RESREF(P1)*TKEIN; RESREF(EP)=RESREF(P1)*EPIN
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,1.0); REAL(DTF);DTF=5.*ZWLAST/WIN/NZ
RELAX(V1,FALSDT,DTF); RELAX(W1,FALSDT,DTF)
RELAX(KE,FALSDT,DTF); RELAX(EP,FALSDT,DTF)
    GROUP 18. Limits on variables or increments to them
VARMIN(C1)=1.E-10
    GROUP 21. Print-out of variables
OUTPUT(BLOK,N,N,N,N,N,N);OUTPUT(PRPS,N,N,N,N,N,N)
    GROUP 22. Monitor print-out
IZMON=NZ/2;IYMON=NYGL;YPLS=T
    GROUP 24. DUMPS FOR RESTARTS
LSWEEP=30;TSTSWP=-1;NPRINT=LSWEEP;NYPRIN=1;NZPRIN=1
  ** Define fluid & solid properties via Q1
  ** mat no. rho enul cp kond expan
  **   1        gas
  **   2        solid
CSG10='q1'
  MATFLG=T;NMAT=2
  34 1.0 1.E-5 1.E3 0.01 0
  152 2000 1.0  960 1.0 0