Тепло и массообмен в пограничных слоях - Патанкар С.
Скачать (прямая ссылка):
IF(KRAD.NE.O)R(I) = R(1) +Y(I) *CSALFA
57 Продолжение Возврат Окончание
SLIP
Подпрограмма SLIP
COMMON/GEN,PE I,AMI,AME,DPDX,PREF(2) ,PR (2) ,P (2) ,DEN, l'AiMU,XU,XD,XP,XL,DX,INTG,CSALFA 1/I/N,NP1,NP2,NP3,NEQ,NPH,KEX,KIN,KASE,KRAD 1/V/U (43) ,F (2,43) ,R (43) ,RHO (43) ,OM(43) ,Y (43)
1/B/BETA,GAMA (2) ,TAUI,TAUE,AJI (2) ,AJE (2),INDI(2),INDE (2) COMMON /l/ak.almg
1 /C/SC (43) ,AU (43) yBU (43) ,CU (43) ,A (2,43) ,B (2,43) ,C (2,43)
Коэффициенты скольжения вблизи I-границы для уравнения движения CU (2) =0.
CU (NP2) =0.
Переход к (71,72,73),KIN
71 BU (2) =0.
AU(2)=1./(1. + 2.*BETA)
Переход к 74
72 SQ = 84.*U(1)*U(1) — 12.*U(1)*U(3)+9.*U(3)*U(3) BU(2)=8.*(2.*U(1)+U(3))/(2.*U(1)+7.*U(3) + SQRT(SQ))
AU(2> = 1.—ВЬТ(2)
Переход к 74
73 BU(2) =0.
CALL VEFF (2,3,EMU)
AK1 = 1./DX—DPDX/(RHO(l)*U(l)*U(l))
AK2=—U(l) *AK1 +DPDX/(RHO(l) *U(1))
AJ = RHO (1) *U (1) *.25* (Y (2) + Y (3)) * *2/EMU IF (KRAD.EQ.O) Переход к 75 AU(2)=2./(2.+AJ*AKl)
CU (2) =—,5*AJ*AK2*AU (2)
Переход к 74
75 CU(2) = 1./12.+3.*AJ*AK1)
AU (2) =CU(2) * (2,—AJ*AK1)
CU (2) =— CU(2)*4.*AJ*AK2
Коэффициенты скольжения вблизи Е-границы для уравнения движения
74 Переход к (81,82,83) КЕХ
81 AU(NPL) =0.
BU(NP2) = l./(l. + 2.*BETA)
Перехот к 84
82 SQ = 84.*U(NP3)*U(NP3)—12.*U(NP3)*U(NP1) +9.*
U(NP1)*U(NP1)
AU(NP2) ='8 *(2.*U(NP3)+U(NP1))/(2.*U(NP3)+7.*U(NP1) +
+ SQRT (SQ))
BU (NP2) = 1.—AU (NP2)
Переход к 84
83 AU(NP2) =0.
CALL VEFF (NP1.NP2,EMU)
BK1 = l./DX—DPDX/(RHO(NP3) *U(NP3) *U(NP3))
BK2 =—U (NP3) *BK1 + DPDX/(RHO(NP3) *U(NP3))
BJ=RHO(NP3) *U(NP3) *.25* (2.*Y(NP3)—Y(NP1)—Y(NP2)) * *2/EMU CU(NP2) = l./(2.+3.*BJ*BKl)
BU(NP2) =CU(NP2) * (2,—BJ*'BK1)
CU(NP2) =—CU(NP2) *4.*BJ*BK2
84 IF(NEQ.EQ.l)
Возврат
Коэффициенты скольжения вблизи I-границы для других уравнений DO 54 J=l, NPH C(J,2) =0.
C(J,NP2) =0.
Переход к (41,42,43),KIN
41 CALL FBC (XD,J,INDI (J) ,Q I)
IF(INDI(J).EQ.l) Переход к 61 АЛ (J) =QI
A(J,2) = 1.
В (J,2) =0.
C(J,2) =6.* (l. + 2.*BETA) *PREF(J) *AJI (J)/(AK*AK*BETA* (l. + BETA) * 1 fl.+BETA) * (3 *iRHO(2) +RHO(3)) *U(3))
Переход к 44 61 F(J,1)=QI
A (J,2) = (1. + BETA—GAMA (J))/(l. + BETA+GAMA (J))
В (J,2) = 1.—A(J,2)
Переход к 44
42 A(J,2) = (U(2)+U(3)— 8.*U(1))/(5.*(U(2)+U(3))+8.*U(1))
GF= (1.—PREF (J)) / (l. + PREF (J))
A(J,2) = (A(J,2) +GF)/(1.+A(J,2)*GF B(J,2) = 1.—A(J,2)
Переход к 44
43 B(J,2) =0.
CALL SOURCE(J,l,CD,DS)
АК1 = l./DX—DS AK2=—AK1*F(J,1)—CS AJF=AJ*PREF(J)
IF(KRAD.EQ.O) Переход к 45 A(J,2)=2./(2.+AJF*AKl)
C(J,2) =—,5*AJF*AK2*A(J,2)
Переход к 44
45 C(J,2) = l./(2. + 3.*AJF*AKl)
A(J,2) = C(J,2) * (2,—2JF*AK1)
C(J,2)=—C(J,2)*4 *AJF*AK2 С Коэффициенты скольжения вблизи Е-границы для других уравнений
44 Переход к (51,52,53),КЕХ
51 CALL FBC (XD,J,INDE(J) ,QE)
IF(INDE(J).EG.l) Переход к 31 AJE(J) =QE
В (J,NP2) =1.
A(J,NP2) =0.
С (J,NP2) =—8.* (l. + 2.*BETA) *PREF(J) *AJE (J) /(AK*AK*BETA* (1.+ 1BETA) * (l. + BETA) * (RHO(NPl) + 3.*RHO(NP2)) *U(NP1))
Переход к 54
31 F(J,NP3)=QE
В (J.NP2) = (1. + BETA—GAMA (J))/(l. + BETA + GAMA (J))
A(J,NP2) =1.—B(J,NP2)
Переход к 54
52 B(J,NP2) = (U(NP2)+U(NPl)-8.*U(NP3))/(5.*(U(NP2)+U(NPl))+8*
U (NP3))
GF= (1.—PREF(J))/(1. + PREF(J))
B(J,NP2) = (B(J,NP2)+GF)/(1.+B(J,NP2)*GF)
A (J,NP2) = 1.—В (J.NP2)
Переход к 54
53 A(J,NP2) =0.
CALL SOURCE (J,NP3,CS,DS)
BK1 = l./DX—DS
BK2=— BK1*F(J,NP3)— CS
BJF=BJ*PREF(J)
C(J,NP2) =l./(2. + 3.*iBJF*BKl)
В (J,NP2) =C(J,NP2) * (2.—BJF*BK1)
С (J,NP2) =— С (J.NP2) *4.*;BJF*BK2
54 Продолжение Возврат Окончание
SOLVE
Подпрограмма SOLVE (A.B,C,F,NP3)
С Решаются уравнения вида
С F(I)=A(I)*F(I + 1)+B (I) *F(I—1) +C(I)
С Для 1 = 2, NP2
Размерность A(NP3) ,В (NP3) ,С (NP3) ,F(NP3)
NP2=NP3—1 B(2)=B(2)*F(1)+C(2)
DO 48 I = 3,NP'2 T= l./(l.—В (I) *A(I—1))
A (I) = A (I) *T
48 B(I) = (B(I)*B(I—1)+C(I))*T DO 50 1 = 2,NP2 J=NP2—1 + 2
50 F(J)=A(J)*F(J + 1)+B(J)
Возврат
Окончание
SOURCE (Пример 1)
С Для сохранения массы химически инертного вещества Подпрограмма SOURCE (J,I,CS,DS)
CS = 0.
DS =0.
Возврат
Окончание
SOURCE (Пример 2)
С Для сохранения энтальпии торможения
С Предостережение — следует использовать согласованные между собой единиц» Подпрограмма SOURCE (J) ,I,CS,DS)
C0MM0N7GEN/PEI,AMI,AME,DPDX,PREF (2) ,PR (2) ,P (2),DEN,
1AMU,XU,XD.XP,XL,DX,INTG,CSALFA
1/V/J (43) .F (2,43) ,R (43) ,RHO (43) ,0M (43) ,Y (43)
I/C/SC (43) ,AU (43) ,BU (43) ,CU (43) ,A (2.43) ,B (2,43) ,C (2,43)
