Тепло и массообмен в пограничных слоях - Патанкар С.
Скачать (прямая ссылка):
IF(KEX.EQ. 1)CALL MASS(XU,XD,AME)
CALL WALL
26 CALL OUTPUT CALL COEFF Подстановка скоростей иа свободных границах
IiF(KEX.EQ.2) U(NP3) = SQRT(U (NP3) *U (NP3) —2.* (XD—XU) * DPDX/RHO (NP3))
IF(KIN.EQ.2) U (1) =SQRT (U (II) *U ((11) —2.* (XD-XU) *DPDX/RHO (1)) CALL SOLVE (AU,BU,CU,U,NP3)
Подстановка скоростей на линии симметрии IF(KIN.NE.3)
Переход к 71 U (1) =U(2)
IF (KRAD.EQ.O) U (1) = ,75*U (2) + ,25*U(3)
71 IF(KEX.EQ.3)U(NP3)=.75*U(NP2)+.25*U(NP1)
72 Uiil=U(l)*U(l)
U13=U(1)*U<3)
U33 = U(3)*U(3)
SQ = 84.*U11 —12 *U13 + 9.*U33
U(2) = (16.*U1 1—4.*U13 + U33) /(2.* (U(l) +U(3)) +SQRT(SQ)) Y(2)=Y(3)*(U(2)+U(3)-2.*U(1))*.5/(U(2)+U(3)+U(1))
Переход к 74
73 IF (KRAD.NE.O)
Переход к 89
U (2) = (4.*U (1)—U (3) )/3.
Y (2) =0.
Переход к 74 89 U(2) = U(1)
Y (2) = Y (3) /3.
74 Переход к (75.70.77) КЕХ
75 U(NP2)=U(NP1)/(1.+2.*BETA)
Y (NP2) = Y(NP3) — (Y (NP3)—Y (NP1)) * BETA/(2. +BET A)
Переход к 78
76 Uil 1 = U (NP1) *U(NP1)
U13 = U(NP1) *U(NP3)
U33=U(NP3) *U (NP3)
SQ=84.*U33—12.*UI13+9.*U11
U('NP2) = (16.*U33—4.*U13 + U11)/(2.* (U(NP1) +U(NP3)) +SQRT (SQ))
Y(NP2) = Y(NP3) — (Y(NP3) —Y(NP1)) * (U(NP2) +U(NP1)—2,*U(NP3)) *.5/ 1 (U (NP2) +U(NP1) +U(NP3)
Переход к 78
77 U(NP2) = (4.*U(NP3)— U(NPl))/3. ,
Y(NP2) =Y(NTP3)
78 Продолжение IF(NEQ.EQ.l) Переход к 45 DO 08 J= 1,NPH
Начальные профили других зависимых переменных READ (5,43) F (J, 1), (F( J, 1), 1 = 3,NP1) ,F (J,NP3)
Вычисление соответствующих величин скольжения GAMA(J) = .143 Переход к (81,82,83),KIN
81 F(J,2)=F(Jjl) + (F(J,3)^F(J,l))*(l.+BETA— GAMA(J))/(1.+
+BETA+GAMA (J))
Переход к 84
82 G=(U(2)+U(3)-8.*U(1))/(5.*(U(2)+U(3))+8.*U(1))
GF= (1.—PREF(J)) /(l. + PREF (J))
GF=(G + GF)/(,1. + G*GF)
F (J,2) =F(J,3)*GF+1.—GF)*F(J,1)
Перехот к 84
83 F (J,2) = F(J,1)
IF(KRAD.EQ.0)F(J,2) = (4.*F(J,1)—F (J,3)) /3.
84 Переход к (85,86,87) КЕХ
85 F(J,NP2) =F(J,NP3) + (F(J,NP1)—F(J,NP3)) * (l. + BETA—GAMA(J))/(1.+
1В ET A +'G AM A (J))
Переход к 88
86 G= (U(NP2) +U(NP1)—8.*U(NP3))/(5.* (U(NP2) +U(NP1)) +8.*U(NP3)) GF= (1.—PREF(J) )/(l.+PREF(J))
GF=(G + GF)/(1+G*GF)
F(J,NP2) = F(J,NP 1) *GF+ (1.—GF) *F(J,NP3)
Переход к 88
87 F (J,NP2) = (4,*F(J,NP3)—F (J,NP 1)) /3.
88 Продолжение
45 Продолжение CALL DENSTY
С Вычисление радиусов
CALL RAD (XU,R (1) ,CSALFA)
IF(CSALFA.EQ.O.OR.KRAD.EQ.O) Переход к 27 DO 28 1 =2,NP3
28 R(l)=R(l)+Y(il)*CSALFA
72 Продолжение
IF(NEQ.EQ.l) Переход к 30 DO 45 J=l, NPH DO 46 '1 = 2,NP2 AU(I)=A(J,1)
BU(I)=B(J,I)
46 CU(1)=C(J,I)
DO 47 1 = 1,NP3
47 SC (1) =F (J,l)
CALL SOLVE (AU,BU,CU,SC,NP3)
DO 48 I = 1.NP3
48 F(J,1) =SC(1)
IF (KASE.EQ.2) Переход к 81 С Подстановка F-величин на стенке
IF(KI'N.EQ.1.AND.INDI(J).EQ.2)F(J,1) = ((1. + BETA2 + GAMA(J))*F 1 (J,2) —(1. + ВЕТА—GAMA(J)) *F(J,3)) *'5/GAMA(J) IF(KEX.EQ.1.AND.INDE(J).EQ.2)F(J,NP3) = ((1. + BETA + GAMA(J))*F 1 (J,NP2) —,('1. + BETA—GAMA (J)) *F (FJ.NPl ))* ,5/GAMA (J)
С Подстановка F-величин на линии симметрии
81 IF(KIN.NE.3) Переход к 82 F(J,1)=F(J,2)
IF(KRAD.EQ.0)F(J,1) =.75*F(J,2) + .25*F(J,3)
82 IF(KEX.EQ.3)F(J,NP3) =.75*F(J,NP2) +.25*F(J,NP1)
45 Продолжение 30 XP = XU XU = XD
PEI ='PEI + DX* (R(l) *AMI—R(NP3) *AME)
С Предельное условие IF (XU.LT.XL)
Переход к 15 Переход к 16 Окончание
BEGIN Подпрограмма BEGIN
COMMON /GEN/PEI,AMI,AME,DPDX,PREF(2),PR(2),P (2),DEN,AMU,
1XU,XD,XP,XL,DX,INTG,CSALFA
1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD
1/B/BETA,GAMA (2) ,TAUI,TAUE,AJ1 (2) ,AJE (2) ,INDI (2) ,INDE (2)
1/V/U (43) ,F (2,43) ,R (43.RHO (43) ,OM (43) ,Y(43)
С Характеристики задачи
READ (5,42)KRAD,NEQ,KEX,KIN,N
42 Формат 411, 12 KASE=2
IF(KIN.EQ.1.0R.KEX.EQ.1)KASE = 1 XU = 0.
NPH = NEQ—1 NP1 =N+il NP2 = N + 2 NP3=N + 3 С Начальный профиль окорости
READ (5,43) Y(1) ,U(1),Y (1) ,U(1) ,1 = 3,NP 1) ,Y (NP3),U (NP3)
43 Формат (6F10.0)
С Вычисление скоростей скольжения и расстояний ВЕТА = .143.
Переход к (71,72,73) KIN 71 U(2)=J(3)/(I.42.*BETA)
Y(2) = Y(3) *ВЕТА/(2. + BETA)
Переход к 74 Переход к 29
27 DO 30 I = 2.NP3
30 R (1) = R (1)
29 Продолжение
С Вычисление (о-величин
0м(1)=0.
ОМ (2) =0.
•DO 49 I=.3,NP2.
49 OM(l)=OM(I— l)+.5*(PHO(I)*U(I)*R(I)+RHO(I— 1)*
1U(I—1)*R(I— 1))*(Y(I) —Y(I— 1))
PEI = OM(NP2)
DO 59 I = 3,NP1 59 OM(I)=OM(I)/PEI OM(NP2) =>1.
