Page 341 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
11:59
May 25, 2005
CB908/Date
0521853265appc
320
STOP 0 521 85326 5 APPENDIX C. 2D CARTESIAN CODE
END
C **************************************
SUBROUTINE TITLE
INCLUDE ’COM2D.FOR’
C **************************************
WRITE(6,*)’*******************************’
WRITE(6,*)’ PROGRAM TO CALCULATE SUDDEN EXPANSION D2/D1=2’
WRITE(6,*)’ KRALL AND SPARROW ’
WRITE(6,*)’*******************************’
RETURN
END
C **************************************
SUBROUTINE INIT
INCLUDE ’COM2D.FOR’
C **************************************
C INITIAL GUESS
RATIO=10*0.563E-3/VISCOS
UBAR=U2
DO 10 I=1,IN
QW(I,JN)=1.0
DO 10 J=1,JNM
U(I,J)=UBAR
U(1,J)=(D2/D1)**2*UBAR
IF(J.GE.JB1)U(1,J)=0.0
IF(TURBUL)THEN
E(I,J)=0.1*0.1*U(I,J)**2
D(I,J)=CMU*DENSIT*E(I,J)**2/VISCOS/RATIO
ENDIF
10 CONTINUE
RETURN
END
C **************************************
SUBROUTINE BSPEC
INCLUDE ’COM2D.FOR’
C **************************************
C **** PROVIDE BOUNDARY & BLOCKED REGIONS
CHARACTER*10 BLOCK,WEST,EAST,SOUTH,NORTH
CHARACTER*10 INFLOW,EXIT1,SYMM,EXIT2,WALLT,WALLQ,PERIOD
DATA BLOCK,WEST,EAST,SOUTH,NORTH
1 /’BLOCK’,’WEST’,’EAST’,’SOUTH’,’NORTH’/