Page 346 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
May 25, 2005
0 521 85326 5
0521853265appc
325
APPENDIX C. 2D CARTESIAN CODE
1 ,3.0,3.3,3.6,3.9,4.2,4.5,5.0,5.5,6.5,7.0,7.5,8.0,8.5,9.0 11:59
1 ,9.7,10.5,11.5,12.5,14.0,16.0,18.0,21.0,25.0,29.0,33.0
1 ,38.0,43.0,48.0,53.0,58.0/
END
Natural Convection Evaporation – Chapter 5
The USER file that follows shows implementation of the mass transfer equation in
subroutine OMEGA. In this subroutine, first coefficients of the discretised equation
(AE, AW, AN, and AS) are evaluated through CALL COEF(0,SC,0.9), where Pr t =
0.9 is inserted though not required in actual calculations because the flow is laminar.
Then, since there is no source term (case of inert mass-transfer), no update of Su
3
and Sp is made. Now, boundary conditions are given where the mass transfer flux
AMW(I,1)atthesouthwallisevaluated.Then,theequationissolvedthroughCALL
SOLVE(0,RPO,RSU), where RPO is the underrelaxation factor. In the ADSORB
subroutine, the source term in the v equation is added to account for buoyancy.
Density is taken to be constant.

C **************************************
C THIS IS USER FILE NATURAL CONVECTION MASS TRANSFER
C **************************************
PROGRAM MAIN
INCLUDE ’COM2D.FOR’
OPEN(UNIT=6,FILE=’OO’)
WRITE(*,*)’-------- output is in OO file --------------’
C **** INITIAL DATA
INM=IN-1
JNM=JN-1
SC=0.614
OBR=50
OWT=50
C
WRITE(*,*)’GIVE ----- MXIT,IREAD,GRM ’
READ(*,*)MXIT,IREAD,GRM
CALL MAINPR
STOP
END
C **************************************
SUBROUTINE TITLE
INCLUDE ’COM2D.FOR’
C **************************************

3 In a combustion problem, source terms must be calculated for each = ω j of interest.