SET INITIAL TRV

CALL TURVS(REL,DELTA,IMAX,UE,DUEDXU,TAU,TRV,DTRV,АКР, 1APZ,ATR,PC0N)

DO 16 I = 1ДИАХ DTRV(I) 0.

16 CONTINUE

COMPUTE CC,AA AND EF

CALL FEPAR(IMAX)

DO 18 I = 1,IMAX DO 17 J = 1,5

17 ABC(J,I) = 0.

18 CONTINUE

BEGINNING OF DOWNSTREAM LOOP

ХСИ = XO

X = XO AIM = IMAP DU = l./AIM KH = 2

DO 27 N = 1,NHAX KST = KH

COMPUTE PRESSURE GRADIENT PARAMETERS

DO 19 KA = KST.NPG К = KA

IF(X .GT. XUE(K))GOTO 19 KH = K-l

IF(K+1 .GT. NPG)K = NPG-1 DUM = XUE(K+1) - XUE(K)

XI = (XUE(K-l)-XUE(K))/DUM X4 = (X - XUE(K))/DUM

X5 = (X + DX - XUE(K))/DUM

DUM = UEE(K-1)-UEE(K)-XI*XIMUEE(K + 1)-UEE(K))

Al = DUM/XI/<1.-XI)

A2 = UEE(K+1)-UEE(K)-A1

UE = UEE(K) + A1*X4 + A2*X4*X4

UEP = UEE(K) + A1*X5 + A2*X5*X5

DUM = UEX(K-1)-UEX(K)-XI*XIMUEX(K+1)-UEX(K))

Al = DUM/XI/(1.-XI)

A2 = UEX(K+1)-UEX(K)-A1

DUEDXU = (UEX(K)+A1*X4+A2*X4*X4)/UE

DUEDXP = (UEX(K) + A1*X5 + A2*X5*X5)/UE

GOTO 20

19 CONTINUE

CALCULATE TRIDIAGONAL TERMS

20 CALL ABCDdMAX.DX,DUEDXU,DUEDXP,UE,UEP,BETA)

CALL BANFAC(ABC.IMAX.1)

CALL BANSOL{D,DTHETA,ABC,IMAX,1)

Рис. 15.13 (продолжение).

Рис. 15.13 (окончание).

122 с

123 с RAT IS THE TYPICAL FRACTIONAL CHANGE IN THETA

124 с

125 RAT = ABS(TAU(IREF)*DTHETA(IREF))

126 IF(0.5*DX .LT. DXMDGOTO 21

127 IF(RAT .GT. RATCH)DX = 0.5*DX

128 UEP = 0.5*(UE+UEP)

129 DUEDXP = 0.5MDUEDXU+DUEDXP)

130 IF(RAT .GT. RATCH)GOTO 20

131 21 X = X + DX

132 IF(1.5*DX .GT. DXMA)GOTO 22

133 IF(RAT .LT. 0.1*RATCH)DX = 1.5*DX

134 22 CONTINUE

135 с

136 с EVALUATE NEW TAU ARRAY

137 с

138 DO 23 I=1,IMAX

139 23 TAUdl = l./(l./TAU(I)+DTHETA(I))

140 CF = 2.*TAU(1)/SQRT(REL)

141 IF(TAU(1) .LT. 0.)GOTO 25

142 с

143 с CALCULATE DISPLACEMENT AND MOMENTUM THICKNESS

144 с

145 DELTA = 0.

146 THKMOM=0.

147 IH = IMAX/2

148 DO 24 lA = 1,IH

149 I = 2*IA - 1

150 AI = I - 1

151 UI DU*AI

152 DELTA = DELTA + 2./TAU(I) + 4./TAU(I+l)

153 THKMOM = THKMOM + 2.*UI/TAU(I) + 4.*(UI+DU)/TAU(I+l)

154 24 CONTINUE

155 DELTA = (DELTA+l./TAU(IMAX)-l./TAU(l))*DU/3.

156 THKMOM = (THKM0M+l./TAU(IMAX))*DU/3.

157 SHPRTR = DELTA/THKMOM

158 DELTA = DELTA/UE/SQRT(REL)

159 THKMOM = THKMOM/UE/SQRT(REL)

160 RTH = REL*UEP*THKMOM

161 с

162 с CALCULATE EDDY VISCOSITY

163 с

164 CALL TURVS(REL,DELTA,IMAX,UEP,DUEDXP,TAU,TRV,DTRV,АКР,

165 lAPZ,ATR,PCON)

166 с

167 IF(X .LT. XCH)GOTO 27

168 XCH = XCH + DXCH

169 25 WRITE(6,26)N,DX,X,TAU{1),CF,DELTA,THKMOM,SHPRTR,RTH,TRV(4)

170 26 FORMATC N=M3, DX=,F5.3, X=\F4.2, TAU (1) = , F5. 3,

171 1* CF=.F7.5, D-TH=,F6.4, M-TH=,F6.4. SH=,F5.3,

172 2* RTH=,E10.3, TRV(4)=,F5.3)

173 с

174 С TEST FOR SEPARATION AND X .GT. XMAX

175 С

176 IF(TAU(1) .LT. 0.)GOTO 28

177 IF{X+0.001 .GE. XMAX)GOTO 28

178 27 CONTINUE

179 28 CONTINUE

180 STOP

181 END

Таблица 15.3. Параметры, используемые в программе DOROD

следующим (внутренним) формулам:

СС .,=, СС = . СС,; =-. (15.67)

\ - и] 1 - Ам

EF , =-2/3, £f , = -0.5 + 0.5 ;--

(15.68)