Тепло и массообмен в пограничных слоях - Патанкар С.
Скачать (прямая ссылка):
Переход к 43
U21 =ABS (.5* (U(2) +U(3)—U(l))
IF(U21.LT.DIF)
Переход к 47
YIP = SQRT(DIF/U21)*(Y(2)+Y(3))*.5 Переход к 44
47 J=2
48 J = J+1
YJI = U(J)—U(l)
IF(ABS (UJI) .GE.DIF)
Переход к 49 Переход к 48
49 Al = l.
IF (UJl.LT.O.)Al =—1.
YIP = Y(J—1) + (Y(J)—Y(J—1))*(U(I)+A1*DIF—U(J—1))/(U(J)— -U(J-l))
Переход к 44
43 YIP = 0.
С Поиск вблизи Е-границы
44 IF(KEX.NE.2) Переход к 45
U21 = ABS (.5* (U (NP1) + U (UP2)) —U (NP3))
IF(U21 .LT.DIF) Переход к 50
YEM = SQRT (DIF/U21) * (.5* (Y(NP1) +Y(NP2))—,(NP3)) + Y(NP3) Переход к 46
50 J = NP2
51 J=J—1
UJI = U(J)— U(NP3)
IF (ABS (UJI) .GE.DIF) Переход к 52 Переход к 51
52 А1 = 1.
IF (UJ1.LT.0.) А1 = — 1.
YEM=Y(J+1) + (Y(J)—Y(J +1)) * (U(NP3) +A1*DIF— U(J+ 1))/(U(J) — -U(J+1))
Переход к 46
45 YEM=Y(NP3)
46 YL = YEM—YIP Возврат Окончание
MASS
С Применима к b в случае непроницаемой стенки Подпрограмма MASS(XU.XD.AM)
АМ = 0.
Возврат
Окончание
OUTPUT
Подпрограмма OUTPUT
COMMON/GEN/PBI,AMI,AME,DPDX,PREF (2) ,PR (2) ,Р (2) ,DEN,
1MU,XU,XD,XP,XL,DX,IiNTG,CSALFA
1 /V/U (43) ,F (2,43) ,R (43) ,RHO (43) ,OM (43), Y (43)
1/C/SC (43) ,AU (43) ,BU (43) ,CU (43) ,A (2,43) ,B (2,43) ,C (2,43)
COMMON /l/ak.almg 1Ц1 /YL,UMAX,UMIN,FR,YIP,YEM 1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD 1/B/BETA,GAMA(2) ,TAUI,TAUE,AJI (2) ,AJE (2) ,IND1 (2) ,INDE (2)
IF(NTG.NE.l) Переход к 15 Запись (6,49) (ОМ(1),1 = 1,NP3
49 Формат (24Н1) величины со есть/(11F10.4))
15 Продолжение
IF (FLOAT INTG— 1) /5..NE.FLOAT ((INTG— 1) /5)) RETURN Запись (6,51)XU
Запись (6,51) (Y(1),1 = 1,NP3),(U(I),I = 1,NP3),(F(1,I),1= 1,NP3)
51 Формат (1P12E10.2)
Возврат
Окончание
PRE
Подпрограмма PRE(XU,XD,DPDX)
COMMON /PR/UGU,UGD 1/V/U (43) F (2,43) ,R (43) ,RHO (43) ,OM (43) ,Y (43) 1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD С Здесь UGU и UGD обозначены скорости внешнего течения в сечениях XU и XD
DPDX= (UGU + UGD) * (UGU—UGD) *.5*RHO(NP3)/(XD—XU)
Возврат
Окончание
RAD (Пример 1)
С Применима к плоским течениям
Подпрограмма RAD(X,R1,CSALFA)
CSALFA= 1.
Rl = l.
Возврат
Окончание
RAD (Пример 2)
С Применима к радиальной пристеночной струе Подпрограмма RAD(X,R1,CSALFA)
CSALFA=0.
R1 = Y
Возврат
Окончание
RAD (Пример 3)
С Применима к осесимметричным слою смешения и струе Подпрограмма RAD(X,R1,CSALFA)
COMMON/GEN/PEI,AMI,AME,DPDX,PREF (2) ,PR (2) ,Р (2) ,DEN,
1AMU,XU,XD,XP,XL,DX,INTG,CSALFA
1/V/U (43) ,F (2,43) ,R (43) ,RHO (43) ,OM(43) ,Y (43)
1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD
CSALFA= 1
IF(KIN.EQ.3) Переход к 17 IF(X.EQ.O.) Переход к 15
Rl = R(1) * (;R(1) —2.*AMI* (X—XP)/(RHO(l) *U(1)))
IF (Rl.LT.O.)'Rl =0.
R1 = SQRT (Rl))
Возврат
15 Rl = 50.
Возврат 17 Rl =0.
Возврат
Окончание
READY
Подпрограмма READY
COMMON/GEN/PEI ,AMI ,AME, DPDX, PREF (2), PR (2) ,P (2), DEN, AMU,
1XU,XD,XP,XL,DX,INTG,CSALFA
1/V/U (43) ,F (2,43) ,R (43) ,RHO (43) ,OM (43) ,Y (43)
1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,K:IN,KASE,K'RAD
1/B/BETA,GAMA (2) ,TAUI,TAUE,AJI (2) ,AJE (2) ,INDI (2) ,INDE (2)
CALL DENSTY
CALL RAD (XU,R (1) ,CSALFA)
Y Вблизи I-границы Переход к (71,72,73),KIN
71 Y(2) = (l. + BETA) *OM(3) *4./( (3 *RHO(2) + RHO(3)) * (U(2) +U (3))) Переход к 74
72 Y(2) = 12.*OM(3)/((3.*RHO(2) + RHO(3)) * (U(2) +U(3) +4.*U(1))) Переход к 74
73 Y(2) = ,5*OM(3) / (RHO(l) *U( 1))
74 Y(3) =Y(2) +.25*OM(3) * (l./(RHO(3) *U(3)) +2./(RHO(3) *U(3) +
+ RHO (2) *U(22)))
Y. Для промежуточных точек сетки DO 50 1 = 4,NP1
50 Y(I) =Y(I—1) +.5* (OM(I)—OM(I—.1)) * (l./(RHO(I) *U(I)) +
М./(RHO(I—.1) *U(I— 1)))
Y Вблизи Е-границы
Y(NP2) = Y(NP1) + .25*(OM(NP2) — OM(NPl)) * (l./(RHO(NPl) *
1U (N'Pil)) +'2.1 (RHO(NP1) *U(NPH) +RHO (NP2) *U(NP2)))
Переход к (81,82,83), КЕХ
81 Y(NP3) = Y(NP2) + (l.+BETA) * (OM(NP2)—OM(NPl)) *
,14./((iRHO (N'Pl) + 3.*RHO (NiP'2)) * (U (NP 1) + U(NP2)))
Переход к 84
82 Y (NP3) = Y (N'P'2) + 12.* (OM(NP2) —OM(NPl)) / ((RHO(NPl) +
43.*RHO(NP2)) * (U(NP2) +U(NPil) +4.*U(N'P3)))
Переход к 84
83 Y(NP3) = Y(NP2) + .5* (OM(NP2)—OM(NPl))/(RHO(NP3) *U(NP3))
84 IF (CSALFA.EQ.O..OR.KRAD.EQ.O) Переход к 51 DO 52 I=2,NP3
52 Y(I) =2.*Y(I) *PEI/(R(1) +SQRT(R(l) *R(I) +2.*Y(I) *
PEI*CSALFA))
Переход к 56
51 DO 54 1=2,NP3
54 Y(l) = PEI*Y (1)/R(1)
56 Y(2) =2.*Y(2)—Y(3)
Y (NP2) =2.*Y (NP2)—Y (NPI)
Вычисление радиусов DO 57 1 = 2,NP3 IF(KRAD.EQ.0)R(I)=R(1)
