GROUP 1. Run title and other preliminaries TEXT(Obstacle,Outlet Restriction And Cart TITLE DISPLAY A cartesian coordinate system is used to encapsulate a rectangular duct,whose topology is as in the previous case. ENDDIS GROUP 2. Transience; time-step specification STEADY=F IREGT=1; GRDPWR(T,5,0.5,1.0) GROUP 3. X-direction grid specification * Write objects RSET(D,DUCT,1.,0.01,0.2,1,0) RSET(B,DUCTWW,0,0.01,0,1.,0,0.2,11,0) RSET(B,OBSTC,0,0,0.1,1.,0.008,0,12,0) RSET(B,OUTLT,0,0,0.2,1.,0.002,0,5,0) GROUP 4. Y-direction grid specification NREGY=3 IREGY=1; GRDPWR(Y,2,0.002,1.0) IREGY=2; GRDPWR(Y,6,0.006,1.0) IREGY=3; GRDPWR(Y,2,0.002,1.0) GROUP 5. Z-direction grid specification NREGZ=2 IREGZ=1; GRDPWR(Z,5,0.1,1.0) IREGZ=2; GRDPWR(Z,5,0.1,1.0) GROUP 7. Variables stored, solved & named ** Solve for three extra variables treated as temperatures with different Prandtl numbers. SOLVE(P1,V1,W1,H1,C1,C2,C3) NAME(C1)=A; NAME(C2)=B; NAME(C3)=C GROUP 8. Terms (in differential equations) & devices TERMS(H1,N,Y,Y,Y,Y,Y) GROUP 9. Properties of the medium (or media) ENUL=1.E-5;PRNDTL(H1)=0.7 PRNDTL(A)=1.0;PRNDTL(B)=0.1;PRNDTL(C)=10.0 GROUP 11. Initialization of variable or porosity fields FIINIT(W1)=0.1 GROUP 13. Boundary conditions and special sources ** Wall WALL (DUCTW,NORTH,#1,#1,#NREGY,#NREGY,#1,#NREGZ,#1,#NREGT) COVAL(DUCTW,W1,1.0,0.0);COVAL(DUCTW,H1,1./PRNDTL(H1),1.0) COVAL(DUCTW,A,1.0,1.0);COVAL(DUCTW,B,10.0,1.0) COVAL(DUCTW,C,0.1,1.0) ** Inlet INLET(INLET,LOW,#1,#1,#1,#NREGY,#1,#1,#1,#NREGT) VALUE(INLET,P1,0.1);VALUE(INLET,W1,0.1) VALUE(INLET,H1,0.0);VALUE(INLET,A,0.0) VALUE(INLET,B,0.0);VALUE(INLET,C,0.0) ** Outlet PATCH(OUTLT,HIGH,#1,#1,#1,#1,#NREGZ,#NREGZ,#1,#NREGT) COVAL(OUTLT,P1,FIXVAL,0.0);COVAL(OUTLT,V1,ONLYMS,0.0) COVAL(OUTLT,W1,ONLYMS,0.0);COVAL(OUTLT,A,ONLYMS,0.0) COVAL(OUTLT,B,ONLYMS,0.0);COVAL(OUTLT,C,ONLYMS,0.0) ** Obstacle PATCH(OBSTC,HIGH,#1,#1,#1,#2,#1,#1,#1,#NREGT) COVAL(OBSTC,W1,FIXVAL,0.0) GROUP 15. Termination of sweeps LSWEEP=20 RESREF(P1)=1.E-6;RESREF(V1)=1.E-6;RESREF(H1)=1.E-6 RESREF(W1)=1.E-6 RESREF(A)=1.E-6;RESREF(B)=1.E-6;RESREF(C)=1.E-6 SPEDAT(SET,GXMONI,TRANSIENT,L,F) GROUP 22. Spot-value print-out NPRMON=20;IYMON=5 GROUP 22. Spot-value print-out IYMON=19;IZMON=10;IPLTL=30 GROUP 23. Field print-out and plot control PATCH(TIMEPLOT,PROFIL,1,1,3*NY/4,3*NY/4,NZ/2,NZ/2,1,LSTEP) PATCH(TIMEPLOT,PROFIL,1,1,3*NY/4,3*NY/4,NZ/2,NZ/2,1,LSTEP) PLOT(TIMEPLOT,W1,0.0,0.2);PLOT(TIMEPLOT,A,0.0,1.0) PLOT(TIMEPLOT,B,0.0,1.0);PLOT(TIMEPLOT,C,0.0,1.0) PATCH(TIMEPLOT,CONTUR,1,1,1,NY,1,NZ,1,LSTEP) PLOT(TIMEPLOT,W1,0.0,10);PLOT(TIMEPLOT,H1,0.0,10) PLOT(TIMEPLOT,A,0.0,10);PLOT(TIMEPLOT,B,0.0,10) PLOT(TIMEPLOT,C,0.0,10)