GROUP 1. Run title and other preliminaries
TEXT(Pipe; Outlet Restriction Added    
TITLE
  DISPLAY
  In this case an area restriction at the outlet has
  been added as well.
 
                     hot pipe wall
               /////////////////////////////////
               --------------------------------|
           --->          - - - - -             |
  uniform flow          -    ||     -  -       |
  at inlet --->   - -  -     ||obstacle  - - ->  exit
             __________ _ ___||____________axis____
 
  ENDDIS
 
    GROUP 2. Transience; time-step specification
  ** Activate unsteady terms in the solved equations
STEADY=F
IREGT=1; GRDPWR(T,5,0.5,1.0)
 
    GROUP 3. X-direction grid specification
CARTES=F
 
   *  Write objects
RSET(D,PIPE,1.,0.01,0.5,1,0)
RSET(B,PIPEW,0,0.01,0,1.,0,0.5,11,0)
RSET(B,OBSTC,0,0,0.25,1.,0.006,0,12,0)
RSET(B,OUTLT,0,0,0.5,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,4,0.004,1.0)
IREGY=3; GRDPWR(Y,4,0.004,1.0)
 
    GROUP 5. Z-direction grid specification
NREGZ=2
IREGZ=1; GRDPWR(Z,5,0.25,1.0)
IREGZ=2; GRDPWR(Z,5,0.25,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)
SOLUTN(P1,P,P,Y,P,P,P)
SOLUTN(H1,P,P,Y,P,P,P)
SOLUTN(C1,P,P,Y,P,P,P)
SOLUTN(C2,P,P,Y,P,P,P)
SOLUTN(C3,P,P,Y,P,P,P)
NAME(C1)=A; NAME(C2)=B; NAME(C3)=C
 
    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 13. Boundary conditions and special sources
 
    1. Wall--- at north boundary
 
WALL (PIPEW,NORTH,#1,#1,#NREGY,#NREGY,#1,#NREGZ,#1,#NREGT)
COVAL(PIPEW,W1,1.0,0.0);COVAL(PIPEW,H1,1./PRNDTL(H1),1.0)
COVAL(PIPEW,A,1.0,1.0);COVAL(PIPEW,B,10.0,1.0)
COVAL(PIPEW,C,0.1,1.0)
 
    2. Inlet--- uniform flow
 
INLET(UNIFORM,LOW,#1,#1,#1,#NREGY,#1,#1,#1,#NREGT)
VALUE(UNIFORM,P1,0.1);VALUE(UNIFORM,W1,0.1)
VALUE(UNIFORM,H1,0.0);VALUE(UNIFORM,A,0.0)
VALUE(UNIFORM,B,0.0);VALUE(UNIFORM,C,0.0)
 
    3. Outlet--- fix pressure at restricted exit
 
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)
 
    4. Obstacle--- at the half length of the pipe
 
PATCH(OBSTC,HIGH,#1,#1,#1,#2,#1,#1,#1,#NREGT)
COVAL(OBSTC,W1,FIXVAL,0.0)
 
    GROUP 14. Downstream pressure for PARAB=.TRUE.
    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
  ** Print out flow field values only from IZ=16 to IZ=20
NYPRIN=2;NZPRIN=1;NTPRIN=1
PATCH(TIMEPLOT,PROFIL,1,1,1,NY,NZ-1,NZ-1,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,20,1,20,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)