Гидравлическое сопротивление при кипении хладагентов в трубах горизонтальных и вертикальных испарителей судовых холодильных машин
Аннотация и ключевые слова
Аннотация (русский):
Проведено сравнение потерь напора при движении кипящего хладагента в горизонтальных и вертикальных испарителях. Рассчитаны, сопоставлены и проанализированы значения всех составляющих гидросопротивления: потери на трение, местное сопротивление, на ускорение потока и нивелирный напор. Получены графические зависимости потери напора двухфазного потока от плотности смеси и высоты аппарата, а также влияния гидросопротивлений на изменение температуры кипения. Отмечено, что исследование эффективности вертикальных испарителей имеет две главные составляющие: определение теплоотдачи и гидросопротивлений при кипении рабочих тел в вертикальных трубах. Рассматривается вторая часть исследования. Вертикальные испарители имеют большую высоту, чем горизонтальные, что приводит к большим потерям при преодолении нивелирного напора большего столба двухфазного потока. Другие составляющие гидросопротивления в горизонтальном и вертикальном аппарате сопоставимы, хотя потери на трение и местные сопротивления в вертикальном могут быть и меньше. В вертикальном аппарате главные потери приходятся на преодоление нивелирного напора, остальные потери меньше, чем в горизонтальном, поэтому при конструировании вертикальных аппаратов нельзя пользоваться соотношением длины аппарата к его диаметру, которое рекомендуется для горизонтальных (L/d для горизонтального и вертикального равно 6). Соотношение для вертикальных аппаратов должно быть существенно меньше. Эксперименты, выполненные на водовоздушной модели, подтвердили, что L/d вертикальных испарителей должно быть в диапазоне 1,0÷2. Более точные значения этого параметра могут быть установлены при проведении экспериментов с хладагентом R410A.

Ключевые слова:
гидросопротивление, двухфазный поток, паросодержание, испаритель, хладагент
Текст
Текст (PDF): Читать Скачать

Введение

Наиболее эффективный и часто применяемый в холодильной технике затопленный испаритель мало применяется в судовых холодильных установках, т. к. требует для заполнения большую массу дорогостоящего фреона; кроме того, в условиях качки при свободном уровне жидкого хладагента в испарителе возможно попадание жидкости во всасывающий тракт компрессора, что приводит к серьезным авариям. Еще одним недостатком этих испарителей является возможность замерзания хладоносителя в трубах при нарушении технологического режима, что приводит к разрыву труб и выходу установки из строя на долгое время.

Применяемые в настоящее время испарители с внутритрубным кипением хладагента выпускаются только в горизонтальном исполнении. Они лишены вышеперечисленных недостатков, но у них ниже коэффициент теплопередачи, поэтому они имеют большую теплообменную поверхность, что увеличивает их габариты и стоимость. В этих испарителях невозможно достичь равномерного распределения хладагента во все трубы пучка, поэтому в некоторых нижних трубах будет избыток хладагента, что может привести к гидроудару в компрессоре, а в верхних – недостаток, что уменьшит эффективное использование теплообменной поверхности.

Вертикальные испарители с внутритрубным кипением фреона лишены указанных недостатков и предполагают высокую эффективность. У них минимальное заполнение хладагентом, замерзание хладоносителя в трубах не возможно, нет свободного уровня жидкости и лучшее распределение хладагента по трубам. Все это характеризует вертикальные испарители с внутритрубным кипением как высокоэффективные и безопасные аппараты с меньшей стоимостью, удобные в эксплуатации и занимающие меньшую площадь, что особенно важно в судовых условиях.

Исследование эффективности вертикальных испарителей имеет две главные составляющие: определение теплоотдачи и гидросопротивлений при кипении рабочих тел в вертикальных трубах. В настоящей работе рассматривается вторая часть исследования.

Вертикальные испарители, конечно, имеют большую высоту, чем горизонтальные, что приводит к большим потерям при преодолении нивелирного напора большего столба двухфазного потока. Другие составляющие гидросопротивления в горизонтальном и вертикальном аппаратах сопоставимы, хотя потери на трение и местные сопротивления в вертикальном могут быть и меньше.

Литературный обзор [1–3] свидетельствует, что с ростом паросодержания гидросопротивление существенно изменяется, при этом в двухфазном потоке оно значительно превышает сопротивление при движении жидкости с той же массовой скоростью. Определение потерь напора в испарителе является важным этапом для расчета аппарата, т. к. это влияет на температуру кипения фреона.

 

Расчет гидросопротивления двухфазного потока

Общий перепад давления между двумя сечениями канала складывается из перепадов, необходимых для преодоления нивелирного напора ∆pнив, потерь на трение ∆pтр, местных сопротивлений ∆pмс, а также потерь на ускорение жидкой и паровой фазы ∆pус:

p = ∆pнив + ∆pтр + ∆pмс + ∆pус.    

Схемы для расчета гидравлических сопротивлений в трубе горизонтального и вертикального испарителя и характер движения двухфазного потока в них приведены на рис. 1.

 

Описание: Безымянный

 

Рис. 1. Схемы для расчета вертикального и горизонтального испарителей:
а – вертикальный испаритель; б – горизонтальный испаритель; 1 – передняя крышка;
2 – задняя крышка; 3 – кожух; 4 – труба; 5 – уровнемер

 

В испарители хладагент подается после регулирующего вентиля, где дросселируется и частично превращается в пар. В условиях расчета массовое паросодержание на входе в испаритель составляет х = 0,15, что соответствует снарядному режиму движения. В трубах горизонтального фреонового испарителя выявлены следующие режимы: снарядный, волновой и расслоенный. Кольцевой режим при тепловых нагрузках, соответствующих работе фреоновых испарителей, не отмечен [1, 4, 5].

В трубах вертикального фреонового испарителя при паросодержании на входе x = 0,15 отмечается снарядный режим, который затем плавно переходит в кольцевой [4, 6, 7].

За базу расчета взят корпус горизонтального фреонового испарителя И-24 диаметром 350 мм, имеющего трубы длиной 2 м, диаметр трубы – 14 мм, среднее расстояние между трубами прямого и обратного хода хладагента s′ = 160 мм (символами ′ и ″ обозначаются параметры, относящиеся к жидкой и паровой фазе соответственно). Этот же корпус используется и в вертикальном аппарате, конструктивная разница аппаратов только в крышках.

Расчет выполнен для хладагента R410А при температуре кипения t0 = –10 оС и удельной тепловой нагрузке q = 5 кВт/м2.

Определяющим фактором для расчета потери напора является массовая скорость ωρ, кг/(м2·c), которая при движении в трубе равна массе хладагента, подаваемой в единицу потока площади сечения трубы в секунду:

                       

где хвх – массовое паросодержание на входе в испаритель после регулирующего вентиля, по диаграмме PI (давление – энтальпия) х = 0,15; f – площадь поверхности трубы, м2; S – площадь сечения трубы, м2; ω – скорость, м/с; ρ – плотность, кг/м3; r – теплота, парообразование, Дж/кг (кДж/кг).

Важным параметром двухфазного потока является истинное паросодержание , иллюстрирующее, какую часть сечения трубы занимает пар. В этих условиях по [2] среднее значение  = 0,87 при значении среднего массового паросодержания  = 0,57.

Тогда нивелирный напор:

                                                 

где lh – высота подъема двухфазного потока, м; g – ускорение свободного падения, м/с2;

pнив. гор = [0,87 · 20 + (1 – 0,87) · 1 209] · 9,81 · 0,16 = 275 Па; pнив. верт  = 3 441 Па.

Согласно расчету pнив. верт велико и значительно больше, чем ∆pнив. гор, и это подтверждает ранее высказанное предположение об определяющем влиянии высоты вертикального аппарата на потери в нем.

Потери давления на трения двухфазного потока [4]:

                                              

где ξ – коэффициент сопротивления потока (ξ = 0,03); ψ – коэффициент, учитывающий влияние потока по [1] (ψ = 1,5); dдиаметр трубы, м; l – длина трубы, м; ω0 – скорость циркуляции, м/с:

                                     

Тогда

В работе [5] представлен другой вид формулы для расчета ∆pтр:  

                                                     

где ψтр – коэффициент, учитывающий влияние потока (ψтр = 0,7).

Тогда

           

Отношение результатов расчетов по [1, 4] составило примерно 20 %.

Потери напора на трение в горизонтальном аппарате больше, чем в вертикальном, вследствие большего значения массовой скорости в нем.

Потери напора на преодоление местных сопротивления для двухфазного потока по [1, 5]:

                                                       

где x – массовое паросодержание в местном сопротивлении; ξм – коэффициент местного сопротивления (принимаем по [6]); ξм для внезапного сужения из крышки в трубу равен 0,5; ξм для внезапного расширения из трубы в крышку равен 1; ξм для поворота на 180°:

                                        

где R – средний радиус поворота потока.

Последовательность местных сопротивлений при движении хладагента в горизонтальном испарителе: сужение из передней крышки в трубу, расширение из трубы в заднюю крышку, поворот на 180°, сужение из задней крышки в трубу, расширение из трубы в переднюю крышку. В вертикальном испарителе: сужение из нижней крышки в трубу, расширение из трубы в верхнюю крышку.

Тогда потеря напора в местных сопротивлениях

                          

Потеря напора в местных сопротивлениях в вертикальных аппаратах меньше вследствие меньшего их числа в нем.

Потеря напора на ускорение рассчитывается как разность количества движения между рассматриваемыми сечениями канала. Потери напора на ускорение по [1]:

                                               

где  у1 – сечение в начале трубы; у2 – сечение в конце трубы.

Тогда потери напора на ускорение

                  

Общая потеря напора в горизонтальном испарителе длиной L = 2 м и диаметром d = 0,35 м:

pгор = ∆pнив + ∆pтр + ∆pмс + ∆pус = 275 + 121 + 63 + 31 = 490 Па.

Общая потеря напора в вертикальном испарителе высотой L = 2 м и диаметром d = 0,35 м:

p = ∆pнив + ∆pтр + ∆pмс + ∆pус = 3 441 + 15 + 6 + 8 = 3 470 Па.

Сравнение подтвердило, что главные потери в вертикальном испарителе – это потери на преодоление нивелирного напора, остальные потери несущественны, меньше, чем в горизонтальном аппарате. Естественно, что для вертикальных аппаратов нужно рекомендовать меньшие значения соотношения длины аппарата и его диаметра, чем те, что приняты в горизонтальных. Кроме того, ∆pнив. верт зависит от плотности двухфазного потока в трубе.

На рис. 2 приведены результаты расчета изменения потерь нивелирного напора в вертикальном аппарате от его высоты и плотности рабочего вещества .

 

 

Рис. 2. Зависимость ∆pнив от средней по длине трубы плотности хладагента при значениях L/d:

 1L/d = 6; 2L/d = 2; 3 L/d = 1

 

С уменьшением плотности смеси и высоты аппарата потери нивелирного напора в вертикальном аппарате значительно сокращаются, если учесть, что значения ρ смеси в первую очередь определяются плотностью теплового потока, то можно сделать вывод, что с увеличением тепловой нагрузки в вертикальном аппарате потери на преодоление нивелирного напора сокращаются, а поскольку это определяющая часть ∆p аппарата, то и общие потери напора будут меньше (рис. 2).

С уменьшением высоты аппарата ∆pнив заметно понижается. Линия 1 соответствует рекомендациям, принятым для горизонтальных аппаратов L/d = 6, что крайне не рационально. Если эту же теплообменную поверхность выполнить с длиной труб L = 0,6 м, при диаметре аппарата 0,60 м, то потери нивелирного напора сократятся в 3–4 раза.

На рис. 3 проиллюстрировано изменение температуры кипения хладагента R410A
от потери напора в испарителях.

 

 

0

р, Па

Описание: Безымянный

 

 

Рис. 3. Зависимость изменения температуры кипения от потери напора в аппарате
при различных температурах

 

Согласно графикам потери напора в высокотемпературных испарителях мало сказываются на температуре кипения, но при низких температурах t0 ниже –30 oC это влияние оказывается заметным, т. к. рекомендуемый температурный напор в этих аппаратах: 7–8 oC, повышение температуры кипения на 0,6–0,8 oC заметно снижает его, соответственно уменьшая производительность аппарата.

 

 

Заключение

Потери давления в горизонтальном испарителе с кипением хладагента внутри труб невелики, все составные части ∆p соизмеримы, и трудно выделить параметр, который бы дал существенное уменьшение ∆p при изменении конструкции аппарата.

В вертикальном аппарате главные потери приходятся на преодоление нивелирного напора, остальные потери меньше, чем в горизонтальном, поэтому при конструировании вертикальных аппаратов нельзя пользоваться соотношением длины аппарата к его диаметру, которые рекомендуются для горизонтальных. Так, в данном примере L/d для горизонтального и вертикального аппаратов равны 6. Это соотношения для вертикальных аппаратов должно быть существенно меньше.

Эксперименты, выполненные на водовоздушной модели, подтвердили, что L/d вертикальных испарителей должно быть в диапазоне 1÷3, более точно значения этого параметра будут уточнены при проведении экспериментов с хладагентом R410A.

В вертикальнотрубных аппаратах с ростом тепловой нагрузки уменьшаются общие потери напора, что не приводит к значительному повышению температуры кипения и уменьшению температурного напора в аппарате.

Список литературы

1. Букин В. Г., Кузьмин А. Ю., Минеев Ю. В. Обобщение экспериментальных данных по гидродинамике и теплообмену при течении двухфазных потоков внутри труб // Вестн. Астрахан. гос. техн. ун-та. 2006. № 6 (35). С. 108-115.

2. Малышев А. А., Киссер К. В., Зайцев А. В. Истинные параметры кипящих хладагентов в трубках и каналах // Вестн. Междунар. академии холода. 2017. № 2. С. 53-56.

3. Ейдеюс А. И., Никишин М. Ю., Кошелев С. В. Потери давления на трение при кипении хладагента в трубах // Вестн. Междунар. академии холода. 2014. № 1. С. 64-67.

4. Кутепов А. М., Стерман Л. С., Стюшин Н. Г. Гидродинамика и теплообмен при парообразовании. М.: Высш. шк., 1986. 447 с.

5. Данилова Г. Н., Богданов С. Н., Иванов О. П. и др. Теплообменные аппараты холодильных установок. Ленинград: Машиностроение, 1986. 303 с.

6. Идельчик И. Е. Справочник по гидравлическим сопротивлениям. М.: Книга по Требованию, 2012. 466 с.

7. Bukin V. G., Kuzmin A. Yu., Mineev Yu. V. Generalization of experimental data on a heat transfer and a flow friction at boiling of mixture ozonefriendly refrigerating agent R407C in a horizontal pipe with twisted tabulators // Science and Technology: International Journal of Scientific Articles “Association of Universities of Pre-Caspian States”. 2006. N. 4. P. 69-73.


Войти или Создать
* Забыли пароль?