2009年度计算流体力学考试题-Read.docVIP

! 一维非稳态导热通用程序(不变部分) ! This is a general purpose program to solve 1-D diffusion. ! problem in the form of: ! ρcdt/dz=1/a(x)d/dx(a(x)Γdt/dx)+s !******************2003.7 revised************************ !...................Define Variables......................... MODULE VARIABLES INTEGER,PARAMETER::L1=130 REAL,DIMENSION(L1):: X,XF,XM,XP,R,RF,AP REAL,DIMENSION(L1)::AE,AW,CN,T,TA REAL,DIMENSION(L1)::TG,GM,RC INTEGER:: K=1,KM=1,KP=1,OM=1 INTEGER:: JB,JE,KE,KI,KF,KN,KR,KT,LS,MD,M1,M2,NF REAL:: AEC,AI,BE,BI,DF,DS,DT,EP,EX,PW,TU,TM,XE,XI REAL:: A1,A2,T2,TC,SC,SP,RO,TE,DN,LM END !..................Main Program............... PROGRAM MAIN USE VARIABLES IMPLICIT NONE INTEGER I OPEN(1,FILE=q.dat) OPEN(2,FILE=temp.dat) NF=1 KN=1 TU=0 150 KT=1 CALL Speci !First to specify the problem CALL Grid !Set up grid points 200 CALL Difsor !Specify the diff-coeff and source term 220 CALL InterOutput !Output intermediate results CALL Coeff !Set up coefficients of discretization equation CALL TDMA !Solve the algebraic equation by TDMA IF(LS.EQ.2.OR.LS.EQ.4) THEN IF(DF.GT.EP) THEN DO I=1,M1 TA(I)=TA(I)+OM*(T(I)-TA(I)) END DO DF=0 KT=KT+1 GOTO 200 END IF END IF CALL GPRINT IF(LS.EQ.3.OR.LS.EQ.4) THEN IF(TU.LT.TM) THEN DO I=1,M1 TG(I)=T(I) END DO KT=1 IF(LS.EQ.3) THEN GOTO 220 ELSE GOTO 200 END IF END IF END IF !.............special results print out, if not, just leave it open.......... CALL SPRINT IF(NF.NE.KM) THEN NF=NF+1 GOTO 150 END IF CLOSE(2) CLOSE(1) END !...................Subroutine................. SUBROUTINE SETUP USE VARIABLES REAL,DIMENSION(L1)::P,Q ENTRY COEFF !................coefficients of boundary points.......... IF(KI.LE.1) THEN AP(1)=1 AE(1)=0 AW(1)=0 CN(1)=AI ELSE AE(1)=GM(1)/XM(2) AP(1)=AE(1)+BI AW(1)=0