**  CHEMKIN - 2D Thermal Diffusion Demo. XY-plane
 
  DISPLAY
    A simple case to demonstrate the effect of thermal
    diffusion on the fields of H2 (in C1) and N2 (in C2).
    The when the temperature gradients are positive, the
    thermal diffusion acts against the bulk mass-flux and
    acts to establish a negative gradient in H2 in the
    cross-flow direction. When the temperature gradient
    is reversed, the thermal diffusion acts with the bulk
    mass-flux causing the overall values of H2 to be smaller
    and a positive gradient to be established in the cross-
    flow direction.
  ENDDISPLAY
 
  PHOTON USE
  p
  phi
 
 
 
 
  view z
  norm
  gr ou z 1
  msg Contours of the temperature
  cont tem1 z 1 fil;.001
  msg
  msg Press  to continue
  pause
  gr ou z 1
  msg Contours of the light species (H2)
  cont h2 z 1 fil;.001
  msg
  msg Press  to continue
  pause
  msg
  msg Press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(Thermal Diffusion Test-Case: XY-plane
TITLE
    GROUP 3. X-direction grid specification
NX=5
GRDPWR(X,NX,1.,1.)
    GROUP 4. Y-direction grid specification
NY=5
GRDPWR(Y,NY,1.,1.)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,C1,C2)
STORE(TEM1,ELH)
    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=T;INIFLD=T
FIINIT(C1)=0.5; FIINIT(C2)=0.5
REAL(TAT1,TATNX);TAT1=0.;TATNX=TAT1-100.
PATCH(WHOLEX,LINVLX,1,NX,1,NY,1,NZ,1,LSTEP)
INIT(WHOLEX,TEM1,TATNX-TAT1,TAT1)
PATCH(WHOLEY,LINVLY,1,NX,1,NY,1,NZ,1,LSTEP)
INIT(WHOLEY,TEM1,TATNX-TAT1,TAT1)
    GROUP 13. Boundary conditions and special sources
INLET(IN1,SOUTH,1,NX,1,1,1,NZ,1,LSTEP)
VALUE(IN1,P1,0.5); VALUE(IN1,V1,0.5)
OUTLET(OUT1,NORTH,1,NX,NY,NY,1,NZ,1,LSTEP)
PATCH(H2SO,CELL,1,NX,NY,NY,1,NZ,1,LSTEP)
COVAL(H2SO,C1,FIXFLU,0.002)
    GROUP 15. Termination of sweeps
LSWEEP=80
    GROUP 16. Termination of iterations
LITER(C1)=20
    GROUP 19. Data communicated by satellite to GROUND
LSG61=T
ENULA=GRND9
ENULB=GRND9
CHSOC=1.
CSG4='tdif'
    GROUP 21. Print-out of variables
NXPRIN=1;NYPRIN=1;NZPRIN=1
    GROUP 22. Spot-value print-out
IXMON=NX/2+1;IYMON=NY/2+1
    GROUP 24. Dumps for restarts
TSTSWP=-1