TEXT(YX BFC UNIFORM FLOW THROUGH A BOX
TITLE
  DISPLAY
  This problem concerns uniform flow through a skewed mesh in a
  rectangular domain. A scalar step profile is convected through
  the domain using the MINMOD scheme for the velocities, and
  van Leer harmonic for the scalar. This case provides a basic
  test of the higher-order schemes on non-orthogonal meshes.
  ENDDIS
  PHOTON USE
  P
 
 
 
   0.20443E+04 0.15633E+04 CR
  gr z 1;vec z 1 sh
  msg Velocity vectors & mesh
  msg press  to continue
  pause
  cl;GR ou Z 1;CON C1 Z 1 FI;.01
  msg van Leer harmonic scheme for scalar concentration
  msg press  to continue
  pause
  ENDUSE
 
REAL(VTOT,YLEN,DTF);VTOT=1.0;YLEN=1.0
   GROUP 6.  Body-fitting and other grid distortions.
BFC=T
GSET(P,P01,0.0000,0.0000,0.0000)
GSET(P,P02,1.0000,0.0000,0.0000)
GSET(P,P03,1.0000,1.0000,0.0000)
GSET(P,P04,0.0000,1.0000,0.0000)
GSET(L,L01,P01,P02,20,1.0)
GSET(L,L02,P02,P03,20,1.0)
GSET(L,L03,P03,P04,20,1.0)
GSET(L,L04,P04,P01,20,1.0)
GSET(F,F01,P01,-,P02,-,P03,-,P04,-)
GSET(D,20,20,1,1.0000,1.0000,0.1000)
GSET(M,F01,+I+J,1,1,1,TRANS)
GSET(P,P01,0.0000,0.0000,0.0000)
GSET(P,P02,1.0000,0.0000,0.0000)
GSET(P,P03,1.0000,0.3000,0.0000)
GSET(P,P04,1.0000,0.6000,0.0000)
GSET(P,P05,1.0000,1.0000,0.0000)
GSET(P,P06,0.0000,1.0000,0.0000)
GSET(P,P07,0.0000,0.6000,0.0000)
GSET(P,P08,0.0000,0.3000,0.0000)
GSET(L,L01,P01,P02,20,-1.2)
GSET(L,L02,P02,P03,7,1.0)
GSET(L,L03,P03,P04,5,1.0)
GSET(L,L04,P04,P05,8,1.0)
GSET(L,L05,P05,P06,20,-1.2)
GSET(L,L06,P06,P07,8,1.0)
GSET(L,L07,P07,P08,5,1.0)
GSET(L,L08,P08,P01,7,1.0)
GSET(L,L09,P08,P03,20,-1.2)
GSET(L,L10,P07,P04,20,1.2)
GSET(F,F01,P01,-,P02,-,P03,-,P08,-)
GSET(F,F02,P08,-,P03,-,P04,-,P07,-)
GSET(F,F03,P07,-,P04,-,P05,-,P06,-)
GSET(D,20,20,1,1.0000,1.0000,0.1000)
GSET(M,F01,+I+J,1,1,1,TRANS)
GSET(M,F02,+I+J,1,8,1,TRANS)
GSET(M,F03,+I+J,1,13,1,TRANS)
GSET(C,K:NZ+1:,F,K1,1,NX,1,NY,+,0.0,0.0,0.1,INC,1.0)
NONORT=T
   GROUP 7.  Variables named, stored & solved.
SOLVE(P1,U1,V1);SOLUTN(P1,Y,Y,Y,N,N,N);SOLVE(C1)
   GROUP 8.  Terms (in differential equations) and devices.
TERMS(C1,N,Y,N,N,P,P)
RHO1=1.161;ENUL=1.5890E-05;ENUT=0.
SCHEME(MINMOD,U1,V1);SCHEME(VANLH,C1)
   GROUP 10. Interphase-transfer processes and properties.
   GROUP 11. Initialization of fields of variables, etc.
FIINIT(UCRT)=0.0;FIINIT(VCRT)=VTOT
FIINIT(U1)=0.0;FIINIT(V1)=VTOT
   GROUP 12. Convection and diffusion adjustments
   GROUP 13. Boundary conditions and special sources
   * INLET boundary condition, name INB
INLET(BFCINA,SOUTH,1,10,1,1,1,#1,#1,#NREGT)
VALUE(BFCINA,U1,GRND1);VALUE(BFCINA,UCRT,0.0)
VALUE(BFCINA,V1,GRND1);VALUE(BFCINA,VCRT,VTOT)
VALUE(BFCINA,P1,GRND1);VALUE(BFCINA,C1,1)
BFCA=RHO1
   * INLET boundary condition, name INB
INLET(BFCINB,SOUTH,11,NX,1,1,#1,#1,#1,#NREGT)
VALUE(BFCINB,U1,GRND1);VALUE(BFCINB,UCRT,0.0)
VALUE(BFCINB,V1,GRND1);VALUE(BFCINB,VCRT,VTOT)
VALUE(BFCINB,P1,GRND1);VALUE(BFCINB,C1,0)
   * OUTLET boundary condition, name OUTT
OUTLET(OUTT,NORTH,#1,#1,NY,NY,#1,#1,#1,#NREGT)
VALUE(OUTT,U1,SAME);VALUE(OUTT,V1,SAME)
VALUE(OUTT,P1,0);VALUE(OUTT,C1,SAME)
   GROUP 14. Downstream pressure (for free parabolic flow).
   GROUP 15. Termination criteria for sweeps
LSWEEP=120
   GROUP 16. Termination criteria for inner iterations.
   GROUP 17. Under-relaxation and related devices.
DTF=YLEN/VTOT/NY
RELAX(P1,LINRLX,0.3);RELAX(U1,FALSDT,DTF)
RELAX(V1,FALSDT,DTF);RELAX(C1,FALSDT,DTF)
   GROUP 22. Location of spot-value.
IXMON=11;IYMON=11;TSTSWP=-1
   GROUP 23. Variable-by-variable field printout and residuals.
ITABL=3;NPLT=10