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