Конденсаторы для динамических систем
В настоящее время расчет конденсаторов проводится обычно по среднему для всей теплообменной поверхности значению коэффициента теплопередачи, которое определяется по сумме тепловых сопротивлений с обеих сторон охлаждающей поверхности. Если тепловые сопротивления со стороны хладоагента и стенки могут быть рассчитаны относительно легко и достаточно точно, то определение коэффициентов тепло — и массоотдачи со стороны парогазовой смеси вызывает большие затруднения. Критический обзор опубликованной по этому вопросу литературы приводит к заключению, что имеющиеся опытные данные пока недостаточны для получения надежных количественных зависимостей, позволяющих определять интенсивность процессов при конденсации пара из парогазовых смесей различного состава и для реальных условий работы теплообменных аппаратов.
При конденсации пара, содержащего даже сравнительно небольшую примесь некондеисирующегося газа, интенсивность процесса в отличие от случая чистого пара определяется не столько скоростью отвода выделяющейся при конденсации теплоты, а главным образом скоростью переноса вещества (пара) к охлаждаемой поверхности, на которой происходит конденсация. Поэтому в случае конденсации пара из парогазовой смеси правильнее говорить о массоотдаче, чем о теплоотдаче, или же, когда разность температур смеси и теплообменной поверхности такова, что отдачей теплоты нельзя пренебречь, — о совместно протекающих процессах тепло — и массоотдачи. Главным источником трудностей, возникающих при иссследовании задачи о конденсации пара из парогазовой смеси является взаимное влияние процессов тепло — и массообмена.
Процесс конденсации пара из парогазовой смеси, как и всякий процесс взаимодействия между твердым телом и омывающей его средой, по физической сущности представляет собой перенос носителей энтальпии, количества движения и вещества в направлении, нормальном к поверхности тела (такое перемещение перпендикулярно поверхности тела принято называть поперечным). Указанное взаимодействие выражается в том, что поверхность получает некоторое количество теплоты и соответственно массы. При таком комплексном процессе, соединяющем в себе явления переноса теплоты и вещества в движущейся среде, возникают дополнительные эффекты, отсутствующие при раздельном их протекании. Совместные процессы оказываются вследствие этого взаимосвязанными.
Существенным для рассматриваемой задачи является еще тот факт, что твердая стенка, непроницаемая для инертного компонента парогазовой смеси, служит поверхностью стоков для конденсирующейся фазы. Вследствие наличия градиента концентраций возникает перенос массы путем диффузии в направлении, перпендикулярном поверхности раздела фаз. Этот диффузионный поток активного вещества должен сопровождаться встречной диффузией инертного компонента. Но, с другой стороны, очевидно, что в силу отмеченного свойства полупроницаемости теплообменной поверхности реальное существование при стационарных условиях постоянного видимого потока инертного компонента в направ — 1G* 243
лении нормали к поверхности является невозможным. Эти два взаимно противоположных требования удовлетворяются благодаря тому, что возникает дополнительный конвективный поток смеси, направленный навстречу диффузионному потоку инертного компонента и компенсирующий его (эффект Стефана).
Нейтрализуя диффузионный перенос инертного газа, поток смеси увеличивает перенос пара по сравнению со случаем, когда имеют место только молекулярная и турбулентная диффузия. Система диффузионных уравнений и граничные условия для рассматриваемого случая должны, следовательно, учитывать наличие поперечного потока активного компонента смеси, а также наличие потока Стефана.
Попытка найти строгое аналитическое решение наталкивается на непреодолимые трудности. Вместе с тем для ряда практически важных задач получено хорошее приближение. Так, С. С. Кутате — ладзе и А. И. Леонтьев [5.12] провели теоретический анализ процесса, в котором применили приближенные методы расчета трения, тепло — и массообмена, основанные на использовании интегральных соотношений импульсов, энергии и диффузии.
Применительно к рассматриваемому явлению тепломассоперено — са эффективными являются методы обобщенного анализа, применение которых позволяет даже в наиболее сложных случаях обнаружить и представить в явной форме связи, скрытые в уравнениях. Конкретная форма этой связи, представленная уравнениями подобия, устанавливается опытным путем.
Л. Д. Берман, анализируя систему дифференциальных уравнений процесса на базе физических представлений, предложил систему чисел подобия, позволяющую обобщить экспериментальные данные для процессов тепломассоотдачи от парогазовой смеси к поверхности конденсации [5.11]:
В общем случае количественная мера результирующего эффекта одновременно протекающих процессов тепло — и ■ массообмена при конденсации пара из парогазовой смеси может быть выражена в виде условного суммарного коэффициента теплоотдачи, определяемого из равенства
Aw А
асум Q 1 ‘■ем ‘ гр
где а — коэффициент теплоотдачи; |3Р — коэффициент массоотдачи, отнесенный к градиенту парциальных давлений; г — теплота фазового перехода; р — парциальное давление; R•—термическое сопротивление; индексы: см — смесь; ст — стенка; п — пар; гр — на границе раздела фаз; к — конвективное; пл—-пленка конденсата.
Так как при конденсации водяного пара теплота фазового перехода значительно превышает количество теплоты, переданной конвекцией от смеси к поверхности раздела фаз, то ак можно пренебречь. Кроме того, можно пренебречь также термическим сопротивлением Ягр, так как добавление в пар инертного газа существенно понижает роль сопротивления на границе раздела фаз.
Особого внимания заслуживает оценка термического сопротивления пленки конденсата Япл. Вопрос, очевидно, не возникает, когда разделение фаз происходит с отсосом образующегося конденсата через пористую стенку. Обычно же определение перепада температуры в пленке конденсата требует проведения большого числа трудоемких опытов при конденсации движущегося чистого пара. Однако, как показали визуальные наблюдения авто* ров, из парогазовой смеси с параметрами, какие обычно имеют место на выходе изТЭ, осуществляется капельная конденсация, в связи с чем парциальное давление пара у поверхности раздела фаз можно определить по температуре стенки. Таким образом, анализ зависимости (5.30) показывает, что с точки зрения инженерной практики для обобщения опытных данных по тепло — и массо — обмену при конденсации пара в присутствии неконденси — рующегося газа в теплообменных аппаратах ЭХГ достаточно знать закономерности изменения коэффициента массоотдачи и соответственно диффузионного числ^ Нуссельта.
Получение искомой зависимости для коэффициента массоотдачи облегчается, если проводить обобщение экспериментальных данных по рекомендации Л. Д. Бермана в относительной форме, используя в качестве масштабных множителей диффузионное и тепловое числа Нуссельта.
Тогда практически независимо от режима течения можно представить в относительной форме и при условии Срп/Ср=Яп/Ят=: с учетом тройной аналогии
Nu _ Niip Nu0 Nudo ®г»
где ng=Ap/p — безразмерная разность парциальных давлений; Ар — разность парциальных давлений активного компонента; р — давление смеси.
В выражении (5.31) Nu и NuD — тепловое и диффузионное числа Нуссельта с учетом поперечного потока массы; Nu0 и NuD0 — тепловое и диффузионное числа Нуссельта для данных условий при исчезающе малом влиянии поперечного потока — масштабные множители. Масштабные множцтели отражают влияние особенностей гидродинамической обстановки (формы поверхности и условий ее обтекания), Nu0 и Nu^o определяются в условиях полной аналогии, что позволяет определить
их по имеющимся в литературе данным для теплообмена при отсутствии массообмена или же при малой его интенсивности, его — объемное содержание газа у стенки; Єгоо — объемное содержание газа в основной массе парогазовой смеси.
Соотношения Кутателадзе — Леонтьева п Бермаиа позволяют значительно сократить число необходимых опытов, решать широкий круг задач с учетом поперечного потока вещества, установить границы аналогии теплоті массообмена и при принятых допущениях хорошо согласуются между собой.
/ — данные работы 15.23]; 2 — данные работы [5.24]; 3— по уравнениям Л. Д. Бермана. |
Опытные данные работы [5.14] были получены в основном па паровоздушной смеси, имеющей различные параметры по давлению, температуре, концентрации воздуха при поперечном обтекании одиночного цилиндра (Re=230-:-2500), турбулентном течении внутри трубы и гравитационно-вязкостном режиме течения в гидродинамическом начальном участке круглой горизонтальной трубы.
На рис. 5.16 в качестве примера представлена обобщенная зависимость, полученная при конденсации водя — ‘. 246
його пара из движущейся паровоздушной смеси при пб — перечном обтекании одиночной трубки и на вертикальной плоской пластине при турбулентном движении смеси внутри прямоугольного канала. Опытные данные удовлетворительно обобщаются уравнениями, полученными на базе уравнения подобия
Nuo/NuTO = 0,7Ь7° ,4 при 0,1 < sr/rcg < 1; (5.32) NuD/Nufl0= 0,7Ь70’^-0Л при а, Ле> 1. (5.33)
Кроме аналитических решений задач тепломассообмена и зависимостей, полученных на основании физических представлений о процессе, можно отметить эффективный метод численного решения системы уравнений пограничного слоя, предложенный С. В. Патанкером и Д. Б. Сполдингом [5.13].
Метод Патанкера—-Сполдинга позволяет решать пространственные задачи совместного тепломассообмена конечно-разностным методом с использованием ломаных профилей; независимой переменной поперек слоя при этом служит безразмерная функция тока.
Перейдем теперь к рассмотрению задачи массообмена в’конденсаторах ЭХГ.
Почти все опытные данные по конденсации пара в присутствии неконденсирующегося газа получены на паровоздушной смеси или смеси с другими тяжелыми газами, и применение рекомендованной обобщенной зависимости для случая конденсации из пароводородной смеси без соответствующей экспериментальной проверки является некорректным, так как сведения о влияний рода газа на интенсивность рассматриваемого процесса весьма разноречивы. Реальные скорости пароводородной смеси в конденсаторах ЗХГ (Re<400) значительно отличаются от имевших место в опытах большинства исследователей (Re>1000). Кроме того, в некоторых ЭХГ применяются малогабаритные теплообменники, процессы в которых проходят при гидродинамически и термически нестабилизированном течении смеси. В литературе практически отсутствуют рекомендации по методике их расчета.
Такого рода теплообменные аппараты для ЭХГ исследованы авторами. Конденсация осуществлялась на продольно обтекаемых трубках с наружным диаметром 3 и длиной 125 мм, расположенных в два ряда по 17 штук,
разделенных перегородкой и заключенных в узкий прямоугольный кожух.
При конденсации пара из парогазовых смесей с большим объемным содержанием газа (ег>0,85), которое, как правило, имеет место в теплообменниках ЭХГ, можно считать справедливой приближенную аналогию между тепло — н массообменом, и для ег/яг> >2,5 при условии jRu/Rv= = 1 справедливо соотношение Nud/Nudo^I — Для выявления величины Nux»o=f(Re, Ргд), определяемой по аналогии между раздельно протекающими процессами тепло — и массообмена, были проведены опыты по конвективному теплообмену на воздухе. Опытные данные, представленные на рис. 5.17, в диапазоне чисел Рейнольдса Re=25-i—350 обобщены уравнением подобия
Nu0=0,048 Re°-8Pr°i33. (5.34)
По аналогии с уравнением (5.34) для диффузионного числа Нуссельта запишем
Nudo=0,048R e°’sP ід0’33. (5.35)
С целью определения возможности применения зависимости (5.35) в расчетах конденсации паров воды из смесей с легкими и тяжелыми газами были проведены четыре серии опытов (воздух, аргон, водород и гелий). Опытные данные с разбросом ±15% удовлетворительно аппроксимируются зависимостью (рис. 5.18)
NuD= 0,048Re°’8Pr^33 (5.36)
Здесь же для сравнения представлены данные работ [5.16, 5.17], полученные при конденсации пара из влажного воздуха в прямоугольных каналах в переходной области чисел Рейнольдса Re> 1000. (Можно заметить, что применение рекомендованных в этих работах расчетных зависимостей приводит к значительному увеличению необходимой тенлообменной поверхности.)
Авторами был испытан малогабаритный конденсатор с поперечным обтеканием шахматного пучка трубок. Трубный пучок был изготовлен из 73 трубок наружным диаметром 4 и длиной 104 мм. Опытные данные удовлетворительно обобщаются зависимостью
(5.37)
Полученная зависимость может быть рекомендована для расчета теплообменных аппаратов аналогичной конструкции. Порядок расчета поверхностного аппарата для конденсации пара из парогазовой смеси может быть следующим.
Обычно в постановке технической задачи известны исходные данные: давление и температура влажного водорода на входе в аппарат, равновесная концентрация электролита, рабочий ток ТЭ, расход охлаждающей жидкости и ее начальная температура. Априори задается конфигурация поверхности теплообмена, выбираемая в зависимости от назначения аппарата. При этом для трубчатых аппаратов задаются диаметром трубок, размером сечений каналов н т. д. Определяется удельная тепловая нагрузка.
Рассчитывается коэффициент теплоотдачи со стороны хладоагента, и определяется температура поверхности конденсации.
По известным формулам для теплообмена определяются Nu0 и Num для данного значения числа Рейнольдса,
По формулам (5.34), (5.36), (5.37) вычисляются значения Nile.
По вычисленным значениям NuD находятся удельная тепловая нагрузка и требуемая поверхность тепломассообмена.
Полученные значения сравниваются с заданными.
Расчет ведется методом последовательных приближений.