Case 4   Cross-flow heat exchanger
 
    GROUP 1. Run title
TEXT(Case 4 Steady Crossflow          
TITLE
    User-defined variables
    T1IN = inlet temperature of shell fluid
    T2IN = inlet temperature of tube fluid
    FLO1 = mass-flow rate of shell fluid
    FLO2 = mass-flow rate of tube fluid
    COEF1 = heat-transfer coefficient on shell-fluid side
    COEF2 = heat-transfer coefficient on tube-fluid side
    COEF12 = overall heat-transfer coefficient from fluid 1
           to fluid 2
 
REAL(T1IN,T2IN,FLO1,FLO2,COEF1,COEF2,COEF12)
T1IN=1.0;T2IN=0.0;FLO1=0.1;FLO2=0.1;COEF1=0.5;COEF2=0.5
COEF12=1.0/(1.0/COEF1+1.0/COEF2)
 
    GROUP 3. X-direction grid specification
 
    The heat exchanger is a rectangular box, 1m high,
    1m wide and 1m long.
 
GRDPWR(X,10,1.0,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,10,1.0,1.0)
 
    GROUP 7. Variables stored, solved , named
SOLVE(H1,H2)
NAME(H1)=1STH;NAME(H2)=2NDH
 
STORE(EPOR,NPOR)
 
    GROUP 8. Terms (in differential equations)  devices
    Cut out the built-in sources, convection and diffusion terms
TERMS(1STH,N,N,N,Y,Y,Y);TERMS(2NDH,N,N,N,Y,N,N)
 
    GROUP 11. Initialization of variable or porosity fields
FIINIT(1STH)=T1IN;FIINIT(2NDH)=T2IN
FIINIT(EPOR)=0.5;FIINIT(NPOR)=0.5
 
    GROUP 13. Boundary conditions and special sources
 
    West boundary; shell fluid inlet ; 10 cells in west wall
PATCH(INLET1,CELL,1,1,1,NY,1,1,1,1000)
COVAL(INLET1,1STH,FLO1/10.0,T1IN)
 
    South boundary; tube fluid inlet; 10 cells in south wall
PATCH(INLET2,CELL,1,NX,1,1,1,1,1,1000)
COVAL(INLET2,2NDH,FLO2/10.0,T2IN)
 
    Flow of shell fluid
PATCH(NEW1,CELL,2,NX,1,NY,1,1,1,1000)
COVAL(NEW1,1STH,FLO1/10.0,LOCNE)
 
    Flow of tube fluid
PATCH(NES1,CELL,1,NX,2,NY,1,1,1,1000)
COVAL(NES1,2NDH,FLO2/10.0,LOCNE)
 
    Heat-exchange with tube fluid, throughout the exchanger.
PATCH(NEPLUS,VOLUME,1,NX,1,NY,1,NZ,1,1000)
COVAL(NEPLUS,1STH,COEF12,LOCNE)
 
    Heat-exchange with shell fluid, throughout the exchanger.
PATCH(NEMINUS,VOLUME,1,NX,1,NY,1,NZ,1,1000)
COVAL(NEMINUS,2NDH,COEF12,LOCNE)
 
    GROUP 14. Downstream pressure for PARAB=.TRUE.
    GROUP 15. Termination of sweeps
RESREF(1STH)=1.E-6*FLO1;RESREF(2NDH)=RESREF(1STH)
LSWEEP=50
 
    GROUP 21. Print-out of variables
    Print-out of porosities is suppressed.
OUTPUT(EPOR,N,N,N,N,N,N);OUTPUT(NPOR,N,N,N,N,N,N)
OUTPUT(1STH,N,N,Y,Y,Y,Y);OUTPUT(2NDH,N,N,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
IXMON=NX/2;IYMON=NY/2;TSTSWP=5
    GROUP 23. Field print-out and plot control
NXPRIN=NX/5;IPLTL=LSWEEP;IPROF=1;ORSIZ=0.4;XZPR=T;NPLT=1
 
PATCH(TABMAP,CONTUR,1,NX,1,NY,1,1,1,1000)
PLOT(TABMAP,1STH,0.0,20.0);PLOT(TABMAP,2NDH,0.0,20.0)
 
    GROUP 24. Dumps for restarts