Научная литература
booksshare.net -> Добавить материал -> Биология -> Патанкар С. -> "Тепло и массообмен в пограничных слоях " -> 62

Тепло и массообмен в пограничных слоях - Патанкар С.

Патанкар С., Сполдинг Д. Тепло и массообмен в пограничных слоях — М.: Энергия , 1971. — 128 c.
Скачать (прямая ссылка): teplomassoobmenvpogransloyah1971.djvu
Предыдущая << 1 .. 56 57 58 59 60 61 < 62 > 63 64 65 66 67 68 .. 76 >> Следующая

CS = SC(I)*(U(I+1)*U(I+1)—U(I)*U(I))/(0M(I+'1)__OM(I))
CS=CS—SO(I—11)*(U(I)*U(I) —U(I—.1) *U(I— 1))/(0M(I)— OM(I—1))
SC = (1.— 1 ,/P REF (J)) *cs/ (ОМ (I +1) —ОМ (I— 1))
DS = 0.
Возврат
Окончание
VEFF
Подпрограмма VEFF (I,IP 1,EMU)
COMMON/GEN/PEI, AMI, AME,DPDX,PR'DF (2) ,PR (2) ,P (2) ,DEN,
1AMU,XU,XD,XP,XL,DX,IiNTG,CSALFA
1 /V/U (43) ,F (2,43) ,R (43) ,RHO (43) ,QM(43) ,Y (43)
1/I/N,NP1.NP2.NP3,NEQ,NPH.KEX,KIN,KASE,KRAD COMMON /L/AK.ALMG 1/LI/YL,UMAX,UMIN,FR,YIP,YEM Эта подпрограмма попользует гипотезу пути смешения AL = ALMG*YL IF(KASE.EQ.2) Переход к 66 IF(KIN.EQ.l) YM= (Y(I) + (IP 1)) *.5 IF(KEX.EQ.1)YM = Y(NP3)— ,5*(Y(I)+Y(1P1))
IF (YM.LT.AL/AK) AL = A'K* YM ,66 EMU = .5* (RHO(1) + RHO(IP1)) *AL*AL*ABS((U(I) —U(IP 1))/(Y(I) —Y(IP 1)) Возврат Окончание
VISCO
Функция VISCO(I)
COM'MON/GE\/PEI,AMI,AME,DPDX,PREF (2) ,PR (2) ,P (2) ,DEN,AMU,XU,XD
1XP,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
VISCO = AMU*(F(l,I)/F(l,NP3))**.76
Возврат
Окончание
WALL
Подпрограмма WALL
COMMON /GEN/РЕ I,AMI,AME,DPDX,PREF (2) ,PR (2) ,P (2),DEN,AMU,XU, 1XD,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 1/B/BETA,GAMA(2) ,TAUI,TAUE,AJI (2) ,AJE (2) ,INDI (2) ,INDE (2)
COMMON L/AK.ALMG Вычисление (3 для Е-границы IF(KEX.NE.l) Переход к 15 YI = Y(NP3)—.5* (Y(NP1) +Y(NP2)
UI = .5* (U(NP2) + U(NP1))
PH = .25* (3,*RHO (NP2) + RHO (NP1))
RE = RH*UI*YI/VISCO(NP3)
FP = DPDX*YI(RH*UI*UI)
AM=—АМЕ/ (RH * U I)
CALL WF1 (RE,FR,AM,S)
BETA = SQRT(ABS(S + FP + AM))/AK TAUE = S * RH * UI * UI IF(NEQ.EQ.I) Переход к 36 Вычисление величин у для Е-границы DO 35 J= 1,NPH
CALL WF2(RE,FP,AM,PR(J),PREF(J),P(J),SF)
GAMA(J) = (SF + AM) *PREF(J)/(AK*AK*BETA)
IF(INDE (J) .EQ.l) AJE (J) =SF*RH*UI* (F (J,NP2) +F(J,NP1)—2*
F(J,NP3))*,5
35 Продолжение
36 IF (KIN.NE.l)
Возврат
Вычисление (3 для I-границы 15YI = .5*(Y(2)+Y(3))
UI = ,5* (U(2) +U(3))
RH=.25* (3.*RHO (2) + RHO (3))
RE = RH*UI*YI/VISCO(l)
FP=DPD*YI/(RH*UI*UI)
AM = AMI/(RH*UI)
CALL WF 1 (RE,FP,AM,S)
BETA = SQRT(ABS (S + FP + AM))/АК TAUI = S*RH*UI*UI IF(NEQ.EQ.l)
Возврат
С Вычисление величин (3 для 1-границы DO 38 J= 1,NP2
CALL WF2(RE,FP,AM,PR,J),PREF(J),P(J),SF)
GAMA(J) = (SF+AM) *PREF(J)/(AK*AK*BETA)
IF(INDI (J) .EQ.l) AJI (J) = SF*RH*UI* (2.*F(J,l)—F(J,2) —F(J,3)) *.5
38 Продолжение Возврат Окончание
WF1
Подпрограмма WF1 (R,F,AM,S)
COMMON /L/AK,ALMG l/WL/SYO,AKS,RT,FT,AMT
AKS=AK*AK RT = R*AKS
ST= 1 ,/RT—,1561*RT* * (—.45) + .08723 *RT* * (—.3) +
+ .03713*RT* * (—.18)
STO = ST
IF(F.EQ.O.) Переход к 15 FT = F/AKS
FM = 1.—4 * FT * RT/ (583. + RT * * 2.5) * * .4 IF (FM.LT.O) FM = 0.
ST = ST FM* * 1.6 Переход к 16
15 IF(AM.EQ.O.) Переход к 16 AMT=AM/AKS
AMM = 1.—AMT/ (7.74 * RT * * (—1.17) + .956* *RT* * (—.25))
ST = ST*AMM**4
16 S = ST*AKS Возврат Окончание
WF2
Подпрограмма WF2(R,F,AM,PR,PRT,P,S)
COMMON /L/AK.ALMG l/WL/STO,AKS,RT,FT,AMT
ST1 = STO/(l. + P*SQRT (SAO)) lF(F.EQ.O) Переход к 15 SSEP = 1.725*RT* * (—.3333) * (P +
+ 6.8) **(—1.165)
FD = .25*FT*RT/(1. + .0625*RT)
ST1 = ST1 *¦( 1.—FD) +FD*SSEP 15 ST = ST1/PRT S=ST*AKS Возврат Окончание
ПРИЛОЖЕНИЕ II
ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ
ПП-1. КРАТКОЕ ОПИСАНИЕ ОПЫТОВ И ЦЕЛЕЙ ЭКСПЕРИМЕНТАЛЬНОГО ИССЛЕДОВАНИЯ
На недостаток экспериментальных данных, пригодных для проверки справедливости соотношений для потоков стенки, указывалось в гл. 4. Основной целью выполненной нами экспериментальной работы является устранение в некоторой степени этого недостатка. В результате измерены диффузионные потоки на стенке, а также профили скорости и концентрации при конечной интенсивности массопереноса через стенку.
Исследовалась радиальная пристеночная струя.
Схема опыта, приведенная на рис. ПП-1, аналогична использованной Бейкером [Л. Б]. На эту работу Бейкера, где содержится подробное описание экспериментального стенда, мы будем в дальнейшем ссылаться неоднократно. Течение в радиальной пристенной струе действительно двумерно и характеризуется нулевым продольным градиентом давления.
В данном случае радиальная пристенная струя формируется на пористой плоской поверхности, через которую возможен массоперенос с конечной скоростью. В наших опытах к вторичному воздуха, проходившему через пористую пластину, подмешивался газ-индикатор гелий. Поверхность пористой пластины была гидродинамически гладкой, ,и свойства движущейся среды оставались почти однородными.
Предыдущая << 1 .. 56 57 58 59 60 61 < 62 > 63 64 65 66 67 68 .. 76 >> Следующая

Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

Есть, чем поделиться? Отправьте
материал
нам
Авторские права © 2009 BooksShare.
Все права защищены.
Rambler's Top100

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed