**  CHEMKIN - 1D Thermal Diffusion Demo. X-dir'n
 
 
    DISPLAY
 
    This simple case demonstrates the effects of thermal
    diffusion on a single light species, H2, diffusing in
    a heavy carrier species, N2. 
    
    A source of H2 is placed
    near the outlet, and N2 flows in at the opposite end.
    
    Thermal (Soret) diffusion is the diffusion of species 
    due to temperature gradients, and the model used here 
    makes the light species diffuse more rapidly towards
    colder surfaces.
 
    ENDDISPLAY
 
    GROUP 1. Run title and other preliminaries
TEXT(C206: 1DX Thermal Diffusion Test
TITLE
    GROUP 3. X-direction grid specification
GRDPWR(X,10,1.,1.)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,C1,C2);STORE(TEM1)
    GROUP 9. Properties of the medium (or media)
PRNDTL(C1)=-GRND9;PRNDTL(C2)=-GRND9
    GROUP 11. Initialization of variable or porosity fields
TEMP0=300.
INIADD=F;INIFLD=T
FIINIT(C1)=0.5; FIINIT(C2)=0.5
  ** Specify negative temperature gradient
REAL(TAT1,TATNY);TAT1=0.;TATNY=TAT1-100.
PATCH(WHOLE,LINVLX,1,NX,1,NY,1,NZ,1,LSTEP)
INIT(WHOLE,TEM1,TATNY-TAT1,TAT1)
INIT(WHOLE,C1,0.2,0.4)
    GROUP 13. Boundary conditions and special sources
  ** Injection of H2 in the outlet cell
PATCH(H2SO,CELL,NX,NX,1,NY,1,NZ,1,LSTEP)
COVAL(H2SO,C1,FIXFLU,1.E-1)
  ** Inflow of pure N2 at west boundary
INLET(IN1,WEST,1,1,1,NY,1,NZ,1,LSTEP)
VALUE(IN1,P1,1.); VALUE(IN1,U1,1.)
  ** Outflow of N2 and H2
OUTLET(OUT1,EAST,NX,NX,1,NY,1,NZ,1,LSTEP)
    GROUP 15. Termination of sweeps
RELAX(C1,FALSDT,10.0)
VARMIN(C1)=1.E-10;VARMIN(C2)=1.E-10
RESREF(C1)=RESREF(C1)/100
LSWEEP=50
    GROUP 19. Data communicated by satellite to GROUND
  Activate the CHEMKIN Interface
LSG61=T
  Set the CHEMKIN input file names
CSG4='tdif'
  Set the reference pressure in atmospheres
CHSOC=1.
  ** Activate "mixture-averaged" species-diffusion model
ENULA=GRND9
  ** Activate species thermal-diffusion model 
ENULB=GRND9
    GROUP 21. Print-out of variables
NXPRIN=1;NYPRIN=1;NZPRIN=1
    GROUP 22. Spot-value print-out
IXMON=5;ITABL=3
    GROUP 24. Dumps for restarts
    Indicate that the temperature should not be constant
   CHKIBEGIN
   CONTEM F
   CHKIEND
TSTSWP=-1
DISTIL=T
EX(P1)=1.000E-03;EX(U1)=1.000E+00;EX(C1)=3.254E-02
EX(C2)=9.675E-01;EX(TEM1)=4.500E+01
STOP