ross=rosi-(rosi-ropsi)*fis !средняя плотность двухфазной смеси в испарителе
wp=1.5*(ssi*(rosi-ropsi)*(9.81/ross)**2)**.25 !скор.вспл.пузырей в большом объеме
cw=.9032 !
wop=cw*wp !скорость жидкости в опускной системе
fop=Vop/(3600.*wop) !площадь сечения опускной системы
!17. ––––––––––––––––––––––––––- Цикл – поиск hop методом половинного деления ––––––-–––––––––
hopi=.2; hopa=1.
dpp=10.
do while(abs(dpp)>.1)
hop=.5*(hopi+hopa) !относит.уровень жидколсти в ОП задаем =0,2...1,0
!18. ––––––––––––––––––––––––––––––––––––––– Баланс давлений –––––––––––––––––––––––––––––––––––––––––––––––––––––
!ode=f(hop) !f(hop) – аппроксимация графиков рис. 5.
ode=.3087 !относительная длина экономайзерного участка
pdv=hr*9.81*(rosi*(hop-ode)-ross*(1.-ode)) !движущее давление циркуляции
cg=.5 !коэффициент гидравл.сопротивления на входе в испаритель
cgm=1. !коэффициент истинного гидравлического сопр.на выходе из испарителя
dp=cg*rosi*w0**2/2. !сопротивление жидкости на входе в испаритель
if(Re>2300.)then;cla=.3164/Re**.25
else;cla=64./Re; endif !коэффициент гидравлического сопр.
dpt=.5*cla*rosi*w0**2/dei !сопротивление трения жидкости на единицу длины
dpe=dpt*hr*ode !сопротивление на экономайзерном участке
im=1000*xxi
if(Re>2300.)then;ck=.25;else;ck=1.; endif !коэфф. для вычисления xx и cks
cc=sqrt(ropsi/rosi)*(amusi/amupsi)**ck/2.; cks=0.
do i=1,im
z=i/1000.
xx=cc*(1./z-1.)**((2.-ck)/2.) !параметр Мартинелли
do j=2,20;if (xxx(j-1)<=xx.and.xxx(j)>=xx)then !аппроксимация (fi ж) рис. 6
fig=fff(j-1)+(fff(j)-fff(j-1))*(xx-xxx(j-1))/(xxx(j)-xxx(j-1))
endif; enddo
cks=cks+(1.-z)**(2.-ck)*fig**2. !интеграл в ф-ле (55)
enddo
cks=1.e-3*cks/xxi !Кси
dpd=dpe*cks*hr*(1.-ode) !сопротивление трения двухфазному потоку
cc=w0+w0vp*(1.-ropsi/rosi)
if(w0p>w0vp)then; cfi=w0vp/(c01*cc+u01) !истинное объемн.паросодерж.вых.исп.
else; cfi=w0vp/(c0*cc+u0); endif паросодержание на вых. исп.
dpu=(rosi*w0v)**2*(xxi**2/(cfi*ropsi)+(1.-xxi)**2/((1.-cfi)*ropsi)-1./rosi) !потери на уск.
dpv=cgm*rosi*w0**2*(1.+w0vp/w0v*(1.-ropsi/rosi))/2. !сопр. двухф.потоку на вых.
dps=9.81*rosi*hst !потери давления на преодоление столба жидкости
dppod=dp+dpe+dpd+dpu+dpv+dps
dpp=pdv-dppod !полезное давление циркуляции – сведение баланса
if(dpp>0.)then;hopa=hop; else;hopi=hop; endif
enddo !–––––––––––––––-конец цикла по hop (пока не достиенем dpp=0.)––––––––––-
19. ––––––––––––––––––––––––––––––––––––––––––––- Результаты расчета ––––––––––––––––––––––––––––––––––––––––––––
write(1,*);write(1,*)
write(1,*)' РЕЗУЛЬТАТЫ РАСЧЕТА'
write(1,*)
write(1,*)' КОНДЕНСАТОР'
write(1,*)' ==========='
write(1,*)'Объемн. доли компонент на входе: xk1v=',xk1v
write(1,*)' xk2v=',xk2v
write(1,*)'Относительные летучести на входе: a1v=',a1v
write(1,*)' a2v=',a2v
write(1,*)'Темп. конденсации в верхнем сеч.: Tkv=',Tkv
write(1,*)'Объемн. доли компонент на выходе: xk1n=',xk1n
write(1,*)' xk2n=',xk2n
write(1,*)'Относительные летучести на выходе: a1n=',a1n
write(1,*)' a2n=',a2n
write(1,*)'Средн.состав пленки - массовые доли: X1=',xxk1
write(1,*)' X2=',xxk2
write(1,*)'Средняя температура пленки: Tk=',Tk
write(1,*)' Свойства тройной смеси'
write(1,*)' ----------------------'
write(1,*)'Tемпература на линии насыщения: ',TSk
write(1,*)'Молярная масса в жидкой фазе: ',AMsk
write(1,*)'Молярная масса в паровoй фазе: ',AMpsk
write(1,*)'Плотность в жидкой фазе: ',rosk
write(1,*)'Плотность в паровoй фазе: ',ropsk
write(1,*)'Теплота фаз. перехода пар-жидк.: ',rsk
write(1,*)'Поверхостное натяжение: ',ssk
write(1,*)'Динамическая вязкость жидкости: ',amusk
write(1,*)' ---------------------'
write(1,*)'Температурный напор конденсации: dTk=',dTk
write(1,*)'Коэфф. теплоотдачи в канале конд.: ',alk
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.