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

write(1,*)'Параметр ребра канала конденсации:       ',amk

write(1,*)'КПД ребра:                               ',pdrk

write(1,*)'КПД пов. конд.:                          ',pdk

write(1,*)'Q конд (=Q исп):                         ',QQ

write(1,*)'Темп. напор в проставочном листе:   dTst=',dTst

write(1,*)

write(1,*)'              ИСПАРИТЕЛЬ'

write(1,*)'              =========='

write(1,*)'Массовые доли жидкости на входе:     Xi1=',xxi1

write(1,*)'                                     Xi2=',xxi2

write(1,*)'Темп. смеси на входе в прогенератор:  Ti=',Ti

write(1,*)'Re:                                      ',Re

write(1,*)'Массовая доля пара на выходе:            ',xxi

write(1,*)'Кратность циркуляции:                    ',ckr

write(1,*)'Относительный расход сливаемой жидкости: ',Go

write(1,*)'Объемн. доли на выходе испарителя:  xi1v=',xi1v

write(1,*)'                                    xi2v=',xi2v

write(1,*)'Объемн. доли на входе в испаритель: xi1n=',xi1n

write(1,*)'                                    xi2n=',xi2n     

write(1,*)'Среднее давление в испарителе:           ',pis

write(1,*)'Средняя температура в испарителе:        ',Tis

write(1,*)'       Свойства тройной смеси'

write(1,*)'       ----------------------'

write(1,*)'Температура на линии насыщения:          ',TSi

write(1,*)'Молярная масса в жидкой фазе:            ',AMsi

write(1,*)'Молярная масса в паровoй фазе:           ',AMpsi

write(1,*)'Плотность в жидкой фазе:                 ',rosi

write(1,*)'Плотность в паровoй фазе:                ',ropsi

write(1,*)'Теплота испарения:                       ',rsi

write(1,*)'Поверхостное натяжение:                  ',ssi

write(1,*)'Динамическая вязкость жидкости:          ',amusi

write(1,*)'       ---------------------'

write(1,*)'Массовый расход на входе:                ',RM

write(1,*)'Модифицированное Re потока от пов-ти:    ',Rem

write(1,*)'Температурный напор кипения:         dTi=',dTi

write(1,*)'Коэффициент теплоотдачи в канале кипения:',ali

write(1,*)'Параметр ребра канала кипения:           ',ami

write(1,*)'КПД ребра:                               ',pdri

write(1,*)'КПД поверхности кипения:                 ',pdi

write(1,*)'Эффективная плотность q в канале кипения:',qei

write(1,*)

write(1,*)'           БАЛАНС ТЕМПЕРАТУР'

write(1,*)'           ================='

write(1,*)'Суммарный перепад температур:        dTs=',dTs

write(1,*)'Действит. эффект. плотность q конд.: qek=',qek

write(1,*)

write(1,*)'           ГЕОМЕТРИЧЕСКИЕ ПАРАМЕТРЫ'

write(1,*)'           ========================'

write(1,*)'Эффективная поверхность конденсации:     ',Fke

write(1,*)'Число каналов конденсатора:              ',nk

write(1,*)'Число каналов в 1 пакете:                ',nk1

write(1,*)'Число пакетов:                           ',np

write(1,*)

write(1,*)'           ГИДРАВЛИКА, ОПУСКНАЯ СИСТЕМА'

write(1,*)'           ============================='

write(1,*)'Привед. скорость пара на вых. испарителя:',w0vp

write(1,*)'Привед. скорость смеси на там же:        ',w0v

write(1,*)'Объемн. расходное паросодержание там же: ',beta

write(1,*)'Объемн. расход жидкости в опускной сист.:',Vop

write(1,*)'Средн. истиное объемн. паросодерж в исп.:',fis

write(1,*)'Скорость жидкости в опускной системе:    ',wop

write(1,*)'Сечение опускной системы:                ',fop

write(1,*)'Относит. уровень жидкости в опускн.сист.:',hop

write(1,*)'Движущее давление циркуляции:            ',pdv

write(1,*)'Сопротивление на экономайзерном участке: ',dpe

write(1,*)'Сопротивление жидкости на входе в исп.:  ',dp

write(1,*)'Интеграл psi:                            ',cks

write(1,*)'Сопротивление трения двухфазному потоку: ',dpd

write(1,*)'Истинное объемн. паросод. на вых. исп.:  ',cfi

write(1,*)'Потери давления на ускорение:            ',dpu

write(1,*)'Сопротивление на выходе из испарителя:   ',dpv

write(1,*)'Сопротивление на столб жидкости:         ',dps

write(1,*)'Гидравлическое сопротивл. подъемн. части:',dppod

return

end

================================================================

!Основная подпрограмма для вычисления температуры

subroutine temp(p,T,x1,x2)

x3=1-x1-x2

B=-364.65*x3-345.48*x2-302.82*x1+x3*x2*(13.44-7.4*(x2-x3)+2.*(x2-x3)**2)+x3*x1* &