GROUP 1. Run title and other preliminaries
 
  DISPLAY
    This is a continuation of library case number 486
    and forms the transient part of the calculation
    in which a dye is 'injected' into the flow field.
    Since it requires the PHIDA file from library ref
    486, it must be run directly after no 486.
  ENDDIS
 
TEXT(Transient Dye Injection Into Steady Flow
TITLE
    GROUP 2. Transience; time-step specification
STEADY=F;TFIRST=0.0;SUBGRD(T,1,1,0.001,1.5)
    GROUP 3. X-direction grid specification
CARTES=F;GRDPWR(X,2,0.01,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,5,0.2,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,5,1.0,1.0)
    GROUP 7. Variables stored, solved & named
STORE(P1,U1,V1,W1);SOLUTN(C1,Y,Y,Y,N,N,Y)
    GROUP 8. Terms (in differential equations) & devices
TERMS(C1,Y,Y,N,Y,Y,N)
    GROUP 9. Properties of the medium (or media)
ENUT=1.E-5 * 100
    GROUP 13. Boundary conditions and special sources
REAL(VELIN);VELIN=20.0
    Inlet 1
PATCH(INLET,HIGH,1,NX,1,2,1,1,1,1)
COVAL(INLET,P1,FIXFLU,VELIN);COVAL(INLET,U1,ONLYMS,0.0)
COVAL(INLET,V1,ONLYMS,0.0);COVAL(INLET,W1,ONLYMS,VELIN)
COVAL(INLET,C1,ONLYMS,1.0)
    Outlet Boundary
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,FIXP,0.0)
    GROUP 15. Termination of sweeps
FSWEEP=1;LSWEEP=30
    GROUP 16. Termination of iterations
REAL(MDOT);MDOT=RHO1*VELIN;RESREF(P1)=MDOT*1.E-12
RESREF(C1)=MDOT*1.E-12;RESREF(U1)=VELIN*MDOT*1.E-12
RESREF(V1)=VELIN*MDOT*1.E-12;RESREF(W1)=VELIN*MDOT*1.E-12
LITER(P1)=25;LITER(C1)=-30
    GROUP 17. Under-relaxation devices
REAL(RELAXN);RELAXN=ZWLAST/VELIN/NZ;RELAX(P1,LINRLX,1.0)
RELAX(U1,FALSDT,RELAXN*5.);RELAX(V1,FALSDT,RELAXN*5.)
RELAX(W1,FALSDT,RELAXN*5.)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 21. Print-out of variables
OUTPUT(C1,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
IXMON=1;IYMON=3;IZMON=3;NPRMON=1;TSTSWP=1
    GROUP 23. Field print-out and plot control
ITABL=3;IPLTL=LSWEEP;IPLTF=FSWEEP;NPRINT=LSWEEP
IZPRF=1;IZPRL=NZ;NUMCLS=8;NPLT=1
RESTRT(ALL);SAVE=T
    GROUP 24. Dumps for restarts