Методика расчета пластинчато-ребристых конденсаторов-испарителей, страница 10

QQ=qek*Fk*pdk   !кол-во выд.теплоты в кан.конд.

dTst=QQ*dst/(2*(b1-2*b3)*hr*alst)   !температурный напор в проставочном листе

!9. ––––––––––––––––––––––––––––––––––––––– Расчет канала кипения ––––––––––––––––––––––––––––––––––––––––––––

eee=1.

qei=QQ/Fi   !эфф.плотность Q в канале кип.

do while (eee>.001)   !––––––––––––-начало цикла по qei–––––––––-

!массовые доли компонентов жидкости на входе

call dolm(xi1,xi2,xxi1,xxi2)

!10.   Температура и свойства кипящей смеси на входе

call temp(pi,Ti,xi1,xi2)

call ssmes(pi,xi1,xi2)

Re=w0*dei*ros/amus   !число Рейнольдса в испарителе

xxi=4.*hr*qei/(ros*rs*dei*w0)   !массовая доля пара на выходе из испарителя

ckr=1./xxi   !кратность циркуляции

Go=Rsl*rs*1.e3/AMs/Q   ! отн. расх. ж. на вых. из кан.кип.

!   Объемные доли компонентов на выходе и входе в каналы кипения

call olet(pi,xi1,xi2,a1,a2)

cc=a1*xi1+a2*xi2+1.-xi1-xi2

c1=a1/cc; c2=a2/cc

cc=ckr*(1.+Go)

xi1v=xi1*cc/(ckr*c1+Go*(ckr+c1-1.))   !объемная доля на выходе из испарителя

xi2v=xi2*cc/(ckr*c2+Go*(ckr+c2-1.))   !- " -

xi1n=xi1v*(ckr+c1-1.)/ckr   !объемная доля на входе в испаритель

xi2n=xi2v*(ckr+c2-1.)/ckr   !- " -

!11.   Усредненные параметры смеси в каналах испарителя

h=6100.   !

pis=pi+.5e-4*h*ros*hr*9.81   !среднее давление в каналах кипения

xi1s=.5*(xi1n+xi1v); xi2s=.5*(xi2n+xi2v)   !средние объемные доли комп. в жидкости

call temp(pis,Tis,xi1s,xi2s)

call ssmes(pis,xi1s,xi2s)

RM=w0*Fsi*ros   !массовый расход жидкости на входе в парогенератор

Rem=qei*dei*ros/(rs*amus*rops)   !модифиц. Re потока массы от теплообм.пов-ти

!12.   Температурный напор на стороне кипения

dTi=.15*QQ*Re**1.1*(dei/hr)**1.1/(RM*cps*Rem**.75)   !

ali=qei/dTi   !коэффициент теплоотдачи при кипении

ami=sqrt(2*ali/(alst*di))   !параметр ребра канала кипения

pdri=.5*ami*hli; pdri=tanh(pdri)/pdri   !КПД ребра канала кипения

pdi=1-(1-pdri)*Fri/Fi   !КПД поверхности кипения

!13.   Новое значение плотности теплового потока

qein=QQ/Fi/pdi

eee=(qein-qei)/qein

qei=qein

enddo  !––––––––––––––––––окончание цикла по qei–––––––––––-

!   свойства тройной смеси в испарителе

TSi=TS; AMsi=AMs; AMpsi=AMps; rosi=ros; ropsi=rops

rsi=rs; ssi=ss; amusi=amus; amupsi=amups

!14. ––––––––––––––––––––––––––––––––––––––––– Баланс температур ––––––––––––––––––––––––––––––––

call temp(pi,Tir,xi1v,xi2v)

if(xi1<.965) Tir=Tis

dTr=Tk-Tir   !располагаемый температурный напор

dTs=dTk+dTst+dTi   !полученный суммарный температурный напор

if(dTs>dTr)then; qeka=qek; else; qeki=qek; endif   !метод половинного деления

enddo   !–––––––––––––––––-окончание цикла по qek–––––––

!15. ––––––––––––––––––––––––––––––––––-Параметры конструкции ––––––––––––––––––––––––––––––––––––––––––––––

Fke=Q/qek   !площадь эффективной поверхности конденсатора

nk=Fke/(Fk*pdk)+1   !число каналов конденсации

nk1=(b2-(hli+2.*db))/(hli+hlk+2.*dst)   !число каналов конденсации в пакете

np=nk/nk1+1   !число пакетов в конденсаторе–испарителе

!16. ––––––––––––––––––––––––––– Гидравлика, опускная система –––––––––––––––––––––––––––––––––––––––––––

w0vp=4*hr*qei/(dei*rsi*ropsi)   !приведеная скорость пара на вых. из испарителя

w0v=w0*(1.-(w0vp/w0)*(ropsi/rosi))   !приведеная скорость смеси на вых. испарит.

beta=w0vp/(w0v+w0vp)   !объемное расходное паросодержание на вых. испарит.

Vop=Fsi*w0v*(nk1+1.)*3600.   !объемный расход жидкости в опускной системе

cmb=10000.; if((cmb*hst**1.5)<10.)then

cmb=6400.; if((cmb*hst**1.5)>5.)then

cmb=2800./(1.-720*hst**1.5); endif; endif   !коэффициент расхода перелива

b=Vop/(cmb*hst**1.5)   !параметр перелива

u0=.35*sqrt(9.81*(rosi-ropsi)*dei/rosi)   !скорость всплытия пузырьков в исп.(fi<=.7)

c0=1.061*w0**(-.065)   !параметр распределения (fi<=.7)

cc=1./(c0*(1.-ropsi/rosi))

w0p=.7*(c0*w0-u0)/(1.-.7/cc)   !приведенная скорость пара при fi=.7

u01=.35   !скорость всплытия пузырьков в испарителе (fi>.7)

c01=1.03   !параметр распределения (fi>.7)

cc1=u0+c0*w0; cc2=u01+c01*w0

fis=cc*(1.-alog10(1.+w0p/(cc*cc1))*cc*cc1/w0vp)   !средн.истинное паросодерж. исп.

if(w0p<w0vp)then; fis=fis+cc*(w0p/w0vp-1.)+cc*(1.-w0p/w0vp- &  ! alog10((1.+w0v/(cc*cc2))/(1.+w0p*c01/(c0*cc*cc2)))*cc*cc1/w0vp)   !паросодерж.

endif