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