Текст подпрограммы теплового и гидравлического расчета
пластинчато-ребристого конденсатора-испарителя PlRKICalc
!Подпрограмма расчета пластинчато-ребристого конденсатора-испарителя
subroutine PlrkiCalc(trez)
!1. ––––––––––––––––––––––––––––––––––––––– Исходные данные ––––––––––––––––––––––––––––––––––––––––––––––––––––––
!свойства смеси компонентов
common /sv4/TS,AMs,AMps,rs,rsk,ros,rops,ss,amus,amups,als,cps
!исходные данные
common/dan/Q,pi,pk,xi1,xi2,yi1,yi2,yk1,yk2,Rsl,A,w0, &
hk,hst,co,hp,b1,b2,b3,hli,hlk,si,sk,di,dk,dst,db
!результаты
common/rez/Ti,Re,pis,TSi,dTi,ali,qei,Tk,TSk,dTk,alk,qek,dTst,QQ,dTs,Fke,nk,nk1,np, &
w0vp,w0v,wop,fop,hop,pdv,dpe,dp,dpd,dpu,dpv,dps,dppod
character(80)trez
DIMENSION xxx(20),fff(20) кусочно-линейная аппроксимация (рис.6)
DATA xxx/.2,.3,.4,.5,.6,.7,.8,1.,1.5,2.1,3.,4.,5.,6.,7.,8.,9.,10.,20.,30./
DATA fff/14.2,9.1,7.,6.08,5.27,4.62,4.06,3.66,3.2,2.85, &
2.6,2.24,2.04,2.,1.85,1.75,1.67,1.58,1.5,1.48/
call smes !константы для расчета свойств смеси
alst=145. !теплопроводность материала стенки
eps=1.e-4 !погрешность расчетов
!2. ––––––––––––––––––––––––––––––-–– Геометрические параметры ––––––––––––––––––––––––––––––––––––––––––––
dei=2*(si-di)*(hli-di)/(si+hli-2*di) !эквивалентный диаметр каналов испарителя
hr=hp-hk-b3 !рабочая высота пакета (в т.ч. парогенерирующего канала)
Fi=2*hr*(b1-2*b3)*(si+hli-2*di)/si !площадь раб. поверхности канала испарителя
Fk=2*hr*(b1-2*b3)*(sk+hlk-2*dk)/sk ! - " - канала конденсатора
Fri=2*hr*(b1-2*b3)*(hli-di)/si ! - " - ребер канала испарителя
Frk=2*hr*(b1-2*b3)*(hlk-dk)/sk ! - " - ребер канала конденсатора
Fsi=(b1-2*b3)*(si-di)*(hli-di)/si !площадь поперечного сечения канала испарителя
!3. ––––––––––––––––––––––––––––––––––––- Расчет канала конденсации ––––––––––––––––––––––––––––––––––––––––
! Объемные доли компонентов в конденсате на входе
xk1v=yk1; xk2v=yk2 !начальное приближение xк1в и xк2в
x1=0.; x2=0. !начальное xк1в и xк2в с предыдущей итерации
do while (abs(x1-xk1v)>=eps.or.abs(x2-xk2v)>=eps) !заголовок цикла, погрешн.
call olet(pk,xk1v,xk2v,a1v,a2v) !относительные летучестеи a1в, a2в
x1=xk1v; x2=xk2v !старые значения xк1в и xк2в
cc=(yk1/a1v+yk2/a2v+1-yk1-yk2)
xk1v=(yk1/a1v)/cc; xk2v=(yk2/a2v)/cc !новые значения xк1в и xк2в
enddo !конец цикла
!4. Промежуточный анализ результатов
call dolm(xk1v,xk2v,xxk1,xxk2) !массовые доли комп. в конденсате на входе
call temp(pk,Tkv,xk1v,xk2v) !температура конденсации в верхнем сечении Tкв
call ssmes(pk,xk1v,xk2v) !свойства смеси в конденсате на входе
!5. Объемные доли компонентов в конденсате на выходе
xk1n=yk1; xk2n=yk2 !начальное приближение xк1н и xк2н
x1=0.; x2=0. !начальное xк1н и xк2н с предыдущей итерации
do while (abs(x1-xk1n)>=eps.or.abs(x2-xk2n)>=eps) !заголовок цикла, погрешн.
do while (abs(x1-xk1v)>=eps.or.abs(x2-xk2v)>=eps) !заголовок цикла, погрешн.
call olet(pk,xk1n,xk2n,a1n,a2n) !относительные летучести a1н, a2н
x1=xk1n; x2=xk2n !старые значения xк1н и xк2н
cc=a1n*x1+a2n*x2+1-x1-x2
c1=a1n*cc; c2=a2n*cc
cc=Q+rsk*A
xk1n=yk1*cc/(Q+c1*rsk*A) !новое значение xк1н
xk2n=yk2*cc/(Q+c2*rsk*A) !новое значение xк2н
enddo
!6. Средний состав пленки конденсата
xk1=(xk1v+xk1n)/2.; xk2=(xk2v+xk2n)/2. !объемные доли
call dolm(xk1,xk2,xxk1,xxk2) !массовые доли
!7. Средняя температура и свойства пленки конденсата
call temp(pk,Tk,xk1,xk2)
call ssmes(pk,xk1,xk2)
TSk=TS; AMsk=AMs; AMpsk=AMps; rosk=ros; ropsk=rops; ssk=ss;amusk=amus
cdTk=(1./als)*(amus/ros/9.81)**.333*(hr/(rsk*amus))**.282 !комплекс для расчета dTк
!8. Цикл по эффективному q при конденсации
qeki=800.; qeka=4500. !мин и макс q
do while(abs(qeka-qeki)>.001) ) !––––––––––––-начало цикла по qek–––––––––-
qek=.5*(qeki+qeka)
dTk=(qek/.95)**1.282*cdTk !темп. напор на стороне конденсации
alk=qek/dTk !коэф. теплоотдачи конд.
amk=sqrt(2*alk/(alst*dk)) !параметр ребра кан.конд.
pdrk=.5*amk*hlk; pdrk=tanh(pdrk)/pdrk !КПД ребра канала конденсации
pdk=1-(1-pdrk)*Frk/Fk !КПД поверхности конденсации
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.