Page 334 - Introduction to Computational Fluid Dynamics
P. 334
P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
313
APPENDIX C. 2D CARTESIAN CODE
Subroutines PR1D This routine is used to print out 1D variables. May 25, 2005 11:59
Function STAN This function routine is called from subroutine BOUND to im-
plement a wall-function boundary condition for temperature and mass-fraction
variables in a turbulent flow. The routine evaluates the Stanton number based on
specification of PF, which the USER may change if required. For mass-fraction
variables, the function is called from routine OMEGA.
Functions FINTW, FINTE, FINTS, and FINTN These function routines evalu-
ate variable at cell faces, w, e, s, and n, respectively, using linear interpolation.
COMMON BLOCK COM2D.FOR
C*******************************
C THIS IS COM2D.FOR
C ******************************
C IT AND JT CHANGE WITH THE PROBLEM
PARAMETER(IT=37,JT=37,MXGR=37)
PARAMETER(GREAT=1.0E+20,SMALL=1.0E-20,PI=3.1415926)
LOGICAL TURBUL,STEADY,UNSTDY,FTRAN,CONMAS,AXISYMM,BSOR
LOGICAL UPWIND,HYBRID,POWER,SLVE,IREAD,IWRITE
LOGICAL GRCELL,GRNODE
C
COMMON/CFLOW/TURBUL,STEADY,UNSTDY,FTRAN,CONMAS,AXISYMM
COMMON/SCHEME/UPWIND,HYBRID,POWER,SLVE(7),IREAD,IWRITE
COMMON/CGRID/IN,JN,INM,JNM,IPREF,JPREF,CORP,NPERIOD,MFREQ
1 ,GRCELL,GRNODE
COMMON/CONTR1/CC,MXIT,CCTM,MXSTEP,DELT,STIME,TTIME
COMMON/CONTR2/RP(9),NSWEEP(7),NITER,RSDU(7),FDIF(7),RNORM(7)
COMMON/CPROP/DENSIT,VISCOS,PR(7),PRT(7),RHO(IT,JT)
1 ,SPHEAT,SPH(IT,JT),VIS(IT,JT),VIST(IT,JT),GAMMA
COMMON/CTURB/CD1,CD2,CMU,ELOG,CAPPA
COMMON/CTAG/NTAG(IT,JT),NTAGW(IT,JT),NTAGE(IT,JT)
1 ,NTAGS(IT,JT),NTAGN(IT,JT)
COMMON/COFV/AW(IT,JT),AE(IT,JT),AS(IT,JT),AN(IT,JT)
COMMON/CSOR/SU(IT,JT),SP(IT,JT),BSOR(9)
COMMON/CAP/AP1(IT,JT),AP(IT,JT),APU(IT,JT),APV(IT,JT)
COMMON/CVAR/U(IT,JT),V(IT,JT),W(IT,JT),P(IT,JT)
1 ,E(IT,JT),D(IT,JT),T(IT,JT),PP(IT,JT),PSM(IT,JT)
COMMON/CVAO/UO(IT,JT),VO(IT,JT),WO(IT,JT),PO(IT,JT)
1 ,EO(IT,JT),DO(IT,JT),TO(IT,JT),RHOO(IT,JT)
COMMON/CORD/X(IT),Y(JT),XC(IT),YC(JT),R(JT),RC(JT)
1 ,DXMI(IT),DYMI(JT),DXP(IT),DYP(JT),VOL(IT,JT)
COMMON/CHEAT/QW(IT,JT),TAUW(IT,JT),AMW(IT,JT)