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

Размещено на http://www.allbest.ru/

Институт Энергетики АН Молдовы

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

Шит Б.М.

Аннотация

эжектор углерод насос диоксид

Приведен расчет двухфазного эжектора для теплового насоса на диоксиде углерода. Метод расчета базируется на методе, разработанном S.M. Kandil, W.E. Lear, S.A. Sherif, модифицированным с учетом требования задания коэффициента эжекции.

Ключевые слова: двухфазный эжектор, диоксид углерода, тепловой насос.

Rezumat

Este prezentat calculul ejectorului bifazic pentru pompa de cгldurг pe dioxid de carbon. Metoda se bazeazг pe metoda, elaboratг de S. Kandil, W. Lear, S. Sherif, modificatг luвnd оn consideraюie valoarea prescrisг a coefiюientului de ejecюie.

Cuvinte-cheie: ejector bifazic, dioxid de carbon, pompa de cгldurг.

Abstract

It is presented the calculus of the two-phase ejector for carbon dioxide heat pump. The method of calculus is based on the method elaborated by S.M. Kandil, W.E. Lear, S.A. Sherif, and is modified taking into account entrainment ratio as the input for the calculus.

Key words: two-phase ejector, carbon dioxide, heat pump.

Задачи расчета двухфазных эжекторов «жидкость-газ» описаны во многих работах. Это работы Kandil S.M. [1], Lear W.E. [2], Дейч М.Е. [3], Соколов Е.Я. єi Зингер Н.М. [4], Цегельский В.Г. [5], Васильев Ю.Н.[6], Спиридонов Е.К. [7] и многих других авторов.

В работах [1,2] приведена методология расчета сверхзвукового двухфазного эжектора, когда заданы давление на выходе эжектора и конструктивные параметры эжектора, но не рассмотрен случай, когда задан коэффициент эжекции и требуется обеспечить максимальное давление на выходе эжектора. В работах [3-7] отсутствуют методы расчеты эжекторов, когда среда на выходе из сопла двухфазная и теплофизические параметры рабочей и эжектируемой среды переменные в зависимости от давления и температуры. В работах [6,7] расчет эжектора проводится при условии, что процесс в эжекторе изотермический.

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

Для расчета двухфазного эжектора нами использована в качестве базовой методика, описанная в [1,2]. При расчетах была использована программа REFPROP 7 разработанная “American National Institute of Standarts” (NIST).

Расчет эжектора

На рис. 1 представлена схема эжектора.

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

Рис. 1. Схема эжектора

При расчете используем методику, разработанную в [1,2]. Модификация методики [1,2] заключается в том, что при определении давления на выходе смесительной камеры и энтальпии потока на выходе смесительной камеры выбираем координаты параметров состояния потока на выходе камеры на кривой постоянной сухости смеси, находящейся под пограничной кривой. При этом задаемся давлением на выходе смесительной камеры с учетом исключения возникновения условий запирания потока и обеспечения максимального СОР теплового насоса. Чтобы получить параметры потока в критическом сечении сопла, задаются давлением в нем, и условием, что поток является изоэнтропическим. На основании закона о сохранении массы получаем:

,

Задаются давлением в критическом сечении и разностью энтальпий до тех пор, пока число Маха в сечении „nt” не станет меньшим 1. Процесс считается изоэнтропическим. . Параметры среды на выходе из сопла вычисляются, используя уравнения сохранения энергии для адиабатического процесса.

Скорость жидкости на выходе из сопла рассчитываем по формуле:

.

Вычислим отношение площадей из уравнения:

Уравнение (1.3) следует из уравнения (1.5).

,

Затем вычислим энтальпию в конце камеры смешения и по заданному значению сухости смеси в конце камеры смешения выбираем исходя из условия

Условия запирания потока по входу вычисляем по формуле:

Первое место, где может произойти запирание потока, - это сечение (см. рис. 1). В этом месте может произойти запирание потока, если противодавление падает, и достигает скорости звука. Такой режим называется «сверхзвуковой насыщенный режим». В этом случае изменяют , чтобы оно достигло значения давления торможения . В этом случае первичный поток расширяется в смесительной камере, так что вторичный поток ускоряется. При этом скорость вторичного потока перестает зависеть от условий по ходу потока.

Следующий анализ используется для вычисления . Уравнение моментов для контрольного объема, показанного на рис. 1, выглядит следующим образом:

.

Итерационный процесс начинается с задания значения , зная, что , что определяет состояние . Из уравнения энергии следует:

Тогда может быть вычислено как:

Если задаться условием и давлением , то состояние будет полностью определено. Скорость вычисляется по следующему уравнению:

.

Давление должно быть подобрано так, чтобы скорость достигла скорости звука.

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

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

.

Затем задаются давлением , чтобы вычислить .

.

Если координаты точки «» не находятся на линии постоянной влажности близкой к требуемой изменяем до получения необходимого результата. При этом варьируем величиной . Определим число Маха потока в сечении . Для этого воспользуемся формулами:

,

,

 скорость звука в газе на кривой насыщения при давлении ; , - показатели адиабаты диоксида углерода на линии насыщения при давлении и удельной энтальпии (определяем с помощью программы REFPROP.

Рассмотрим характеристики двухфазного потока в трубах и каналах. Общий массовый расход смеси жидкости и газа (kg/s) равен и является постоянной величиной, одинаковой в любом сечении канала. Значения массовых расходных газосодержаний лежат в пределах . Массовые расходы газа соответствуют объемным расходам (m³/s) єi . Их сумма является объемным расходом смеси. В отличие от массового расхода объемный расход смеси в общем случае переменен по длине эжектора и может изменяться от значения если в сечении движется только жидкость, до значения величины - если в сечении движется только газ. Массовое содержание газа (или степень сухости) равна отношению массового расхода газа к массовому расходу смеси. Объемное расходное газосодержание равно отношению объемного расхода газа к объемному расходу смеси:

.[8]

Следует напомнить, что в адиабатических (без подвода тепла) скачках сжатия происходит увеличение энтропии газа, а также то, что в прямом скачке уплотнения всегда сверхзвуковая скорость газа переходит в звуковую. Во многих случаях имеется два возможных решения для потока, выходящего из смесительной камеры. Первое решение - дозвуковая скорость. В этом случае условия на входе в диффузор устанавливаются равными условиям на выходе из смесительной камеры. Второе решение - это сверхзвуковое решение. Для этого случая анализ выполняется двумя способами. Сначала условия на входе диффузора устанавливаются равными условиям на выходе смесительной камеры, и производится изоэнтропический анализ для диффузора. Это решение, которое называют вторым решением, является теоретически верным и подтверждено экспериментально. Вторая методика анализа при сверхзвуковом режиме течения в смесительной камере требует использования условий Рэнкина-Гьюгонио (Rankin-Hugoniot relations) для моделирования ударной волны [1,2].

.

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

Как только дозвуковые или сверхзвуковые условия на входе в диффузор определены, проводится изоэнтропический анализ условий на выходе диффузора.

Скорость смеси на выходе диффузора вычисляется с использованием закона сохранения энергии:

.

Значения величин и выбирают таким образом, чтобы удовлетворяло уравнениям (1.23) и (1.24). Применительно к эжектору на диоксиде углерода становится очевидным, что оптимизация давления осуществляется путем выбора скоростей потоков на входах эжектора и его конструктивных параметров таким образом, чтобы точка с координатами и находилась на линии заданного массового газосодержания, находящейся под пограничной кривой, что может обеспечить заданное значение коэффициента эжекции.

Выводы

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

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

Литература

1. Kandil S.M., Lear W.E., Sherif S.A. Analysis of two-phase supersonic flow in jet pumps. Proceedings of FEDSM2005 (2005 ASME Fluids Engineering Division Summer Meeting and Exhibition June 19-32, 2005 Houston, TX, USA, FEDSM2005-77476.

2. Lear, W.E., Sherif, S.A., Steadham, J.M., 2000, “Design Considerations of Jet Pumps with Supersonic Two-Phase Flow and Shocks for Refrigeration and Thermal Management Applications,” International Journal of Energy Research, Vol. 24, pp. 1373-1389.

3.Дейч М.Е. Газодинамика двухфазных сред. - М.: Энергоиздат, 1981. - 471 с.

4. Соколов Е.Я., Зингер Н.М. Струйные аппараты. - М.: Энергия, 1970.-287 с.

5.Цегельский В.Г. Двухфазные струйные аппараты. - М.: Изд-во МГТУ им Н.Э. Баумана, 2003. -408 с.

6.Васильев Ю.Н.. Теория двухфазного газожидкостного эжектора с цилиндрической камерой смешения. //Лопаточные машины и струйные аппараты. - М.: Машиностроение, 1971. - Вып.5-С.175-261.

7. Спиридонов Е.К. Энергетический анализ жидкостно-газовых течений. Вестник ЮУрГУ, №1, 2003. Серия «Машиностроение», вып.1 с. 143-150.

8. Благов Э.Е. Критическое отношение давлений и критическая скорость при течении однокомпонентной вскипающей жидкости через сужающие устройства. Теплоэнергетика, 2005, №6, с.56-66.

9. Шит М.Л., Балануцэ А.П., Шит Б.М. Промышленная теплонасосная установка на диоксиде углерода с испарителями, работающими на различных температурных уровнях и при переменной нагрузке. Problemele Energeticii Regionale, N2, 2010, p.47-58.

Шит Борис Михайлович - программист Лаборатории энергоэффективности и систем управления ИЭ АНМ. Область научных интересов: проектирование систем управления в экономических и технических системах, сетевое программирование; автор или соавтор более 20 статей и одного патента. E-mail: boris_fld@rambler.ru

Приложение

Пример расчета

Рассмотрим расчет эжектора для теплового насоса, описанного в [9].

Исходные данные: коэффициент эжекции ;

Свойства рабочего потока:

Свойства эжектируемого потока:

Подбираем и с учетом того, что процесс между состояниями „pi” и „nt” является изоэнтропическим и чтобы удовлетворяло уравнениям (А1.1) и (А1.2). Координаты состояния «nt» будут следующими:

Вычислим:

Принимаем:

.

Как показали расчеты, значение величины очень чувсвительно к изменению величин и () и поэтому должно быть стабилизировано с помощью регулируемого сопла.

Выберем давление и температуры среды на выходе из сопла и определим теплофизические свойства среды в сечении «ne».

Скорость газа в сечении «ne»

Определим отношение сечений по формуле:

Зададимся значением коэффициента эжекции По формуле (1.15) в которой индексы „me” заменены индексами „ne”(см. [8]) вычислим скорость звука в сечении „ne”: Так как режим работы эжектора дозвуковой.

Выберем давление в сечении , при условии, что Тогда скорость потока в сечении “se” будет равна:

Отношение сечений alegem egal cu . Зададимся давлением

Тогда:

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

Вычислим значение величины .

Вычислим критическую скорость звука по формуле:

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

Проверяем выполнение условия . Если скорость смеси меньше скорости звука, возвращаемся в начало расчета и изменяем или геометрические параметры эжектора или скорости сред на входе в эжектор. Так как критическая скорость звука, вычисленная по формуле (А1.9) больше скорости , вычисленной по формуле (А1.9), то режим течения является дозвуковым и вычисление чисел Фабри не требуется. Тем не менее, вычисление чисел Фабри показало, что при выбранных параметрах эжектора запирание потока не происходит. Размещено на Allbest.ru