Тепло и массообмен в пограничных слоях - Патанкар С.
Скачать (прямая ссылка):
OM(,NP3) =1.
IF(NEC.EQ.l) RETURN DO 69 J=,l, NPH IF(KEX.EQ.1)INDE(J) = 1 IF(KIN.EQ.l) INDI(J) = 1 69 Продолжение Возврат Окончание
COEFF Подпрограмма COEFF
COMMON/GEN/PEI,AMI,AME,DPOX,PREF (2) ,PR (2) ,P (2),DEN,AMU, lXU,XD,XP,XL,DX,rNTG,CSALFA 1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD l/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)
1/C/SC (43) ,AU(43) ,BU(43) ,CU (43) ,A(2,43) ,B (2,43) ,C (2,43)
COMMON /L/АК, ALMG
Размерность G1 (43) ,G2 (43),G3(43) ,D(2,43) ,S 1 (43) ,S2 (43) ,S3(43)
С Вычисление малых С-величин DO 98 I = 2,NP 1 RA = .5*(R(I+1)+R(I))
PH = .5* (RHO (I +il) + RHO (I))
U'M = .5*(U(I+1)+U(I))
CALL VEFF(I,I+,l,EiMU)
98 SC(I)=RA*RA*RH*UM*EMU/(PEI*PEI)
С Член конвекции
SA = R(1) * AMI/PEI
SB=(R(N'P3)*AME—R(1)*AMI)/PEI
DX = XD—XU
DO 71 1=3,NP1
OMD = ОМ (I + 1) —ОМ (I— 1)
P2 = .25/DX P3—P2/OMD
Pl = (OM(I+'l)— OM(I))*P3 P3=.(OM(I)—OM(I—1))*P3 P2 = 3.*P2 Q = SA/OMD R2 =—SBA.25 R3 = Ri2/OMD
R1 =—(OM(I + 1) + 3.*OM(I))*R3 R3= (OM(I—<1) +3.*OM(I)) *,R3 GI (I) =P1 +Q + R1 G2(I)=P2 + R2 G3(I)=P3—Q + R3
CU(I)=—P1*U(I+I1)—P2*U(I)—P3*U(I—1)
С Член диффузии
AU(I) =2./OMD
BU(I)=SC(I— l)*AU(I)/(OM(I)— OM(I—1))
AU(I) =SC(I) *AU(I)/(OM(I+1)—OM(I)
IF(NiEQ.EQ.I) Переход к 33 DO 34 J= l.NPH
C(J,I)=— P1*F(J,I+1)—P2*F(J,I)— P3*F(J,I—1)
CALL SOURCE(J.I,CS,D(J,I))
С (J,I) =—С (J,I) +CS—F (J, I) *D(J,I)
A(J,I) =AU(I)/PREF(J)
В (J,I) =BU (I)/PREF (J)
34 Продолжение Г Член источника в уравнении движения 33 SI(I)=DPDX*DX
S2 (I) =P2*S1 (I)/(RHO(I) *U(I))
S3(l) =P3*S 1 (I)/(RHO (I—1) *U(I—1))
S1(1)=P1*S1(I)/(RH0(I + 1)*U(I+'1)
CU(1)=—CU(I)— 2 *(S1(I)+S2(I)+S3(I)
Sl(I)=Sl(I)/(I+'l)
S2(I) = S2 (I) /U (I)
S3(I)=S3(I)/U(I—1)
71 Продолжение С Коэффициенты в окончательном виде DO 91 1 = 3,NP1
RL= 1./ (G2(I) +AU (I) +BU(I)—S2(I))
AU(I) = (AU(I) +SI (I)—G1 (I))*RL BU(I) = (BU(I)+S3(I)—G3(I))*RL
91 CU(I) =CU(I)*RL IF(NEQ.EQ.l) Переход к 76 DO 92 J=1,NPH
DO 92 1=3,NP1 RL = 1 -/(02 (I) + A(J,I) +B(J,I)—D(J,I))
A(J,I) = (A(J,I)— G1(I))*RL B(J,I) = (B(J,I)— G3(I))*RL
92 С (J,I) =C (J,I) *RL
76 CALL SLIP
Возврат
Окончание
CONST Подпрограмма CONST
COMMON/GEN/PEI,AMI,AME,DPDX,PREF (2) ,PR (2) ,P (2),DEN,AMU,XU,
11XD,XP,XL,DX,INTG,CSALFA
COMMON/L/AK.ALMG
1/LI/YL,UMAX,UMIN.FR,YIP,YEM
AK = -435
ALMG=.09
FR = .01
PREF(l) =.9
PR(1)=.7
PRRAT = PR (1)/PREF(1)
P (1) =3.68* (PRRAT—1.) *PRRAT** (—.25)
DEN = .075 AMU = .000012 Возврат Окончание
DENSTY
Подпрограмма DENSTY
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) l/I/N,NPa,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD DO 45 1 = 1,NP3
45 RHO (I) =DEN*F(1,NP3)/F (4,1)
Возврат
Окончание
ENTRN
Подпрограмма ENTRN
COMMON /GEN/PEI, AMI, AME, DPDX, PREF(2), PR (2), P (2), DEN, AMU, 1XU,XD,XP,XL,DX,IiNTG CSALFA COMMON /L/AK.ALMG l /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 l/Ll/YL,UMAX,UMIN,FR,VIP,YEN - Эта подпрограмма использует гипотезу пути смешения Переход к (71,72,73),KIN
71 Переход к 74
72 AMI = 8.*RHO(l) * ((ALMG*YL)/(Y(2) +Y(3))) * *2*ABS (U(2) +
+ U(3)-2.*U(1))
Переход к 74
73 AMI = 0.
74 Переход к (81,82,83) КЕХ
81 Возврат
82 AME = 8.*RH0(NP3)*((ALMG*YL)/(Y(NP1)+Y(NP2) —
12.*Y (NP3))) **2*ABS(U(NP1) +U(NP2)—2.*U(NP3))
Возврат
83 AME = 0.
Возврат
Окончание
FBC (Пример 1)
С Предписывает ступенчатое изменение температуры стенки Подпрограмма FBC (X,J,IND,AJFS)
IND=il
AJFS = 100. + 5.*X
Возврат
Окончание
FBC (Пример 2)
С Предписывает однородный поток
Подпрограмма FBC (X,J,IND,AJFS)
IND = 2 AJFS = 50.
Возврат
Окончание
LENGTH
Подпрограмма LENGTH
COMMON/GEN/PEI, AMI, AME,DPDX,PREF(2),PR (2), DEN, AMU,
UXU,XD,XP,XL,DX,IiNTG,CSALFA
1 /V/U (43) ,F (2,43) ,R (43) ,RHO (43) ,OM (43) ,Y (43)
1/I/N,NP1,NP'2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD 1/Ll,YL, UMAX, UMIN,FR,YIP,YEM С Отыскание максимальной и минимальной скоростей
40 UMAX = U (1)
UMIN = U(1)
DO 41 J=3,NP3 IF(J.EQ.NP2) Переход к 41
IF (U (J) .GT.UMAX) UMAX=U (J)
IF(U (J) .LT.UMIN) UMIN = U (J)
41 Продолжение
DIF=ABS (UMAX—UMIN) *FR С Поиск вблизи I-границы IF (KIN.NE.2)
