Моделирование теплообмена в воздушно-испарительном охладителе при воздействии электростатического поля

УДК 621.573

Моделирование теплообмена в воздушноиспарительном охладителе при воздействии электростатического поля

Канд. техн. наук В.В.КИРЕЕВ ОАО «Мясокомбинат «Тулунский»

The use of electrostatic field (ESP) in the cooling of heat carriers in an air-evaporative chiller allows to increase the heat transfer coefficient and reduce consumption of water.

The offered calculation algorithm as obtained on the basis of the analysis of the constraint equations of the surface energy, the intensity of the ESP, radius of drops and their electric and thermophysical characteristics can be used as the basis of engineering principles of refrigerating systems construction.

Применение электростатического поля (ЭСП) при охлаждении теплоносителей в воздушно-испарительном охладителе позволяет увеличить коэффициент теплоотдачи и сократить расход воды. Струя диспергированной жидкости (СДЖ) в ЭСП может двигаться при определенной разности потенциалов в присутствии электродов в зоне охлаждения. Диспергирование жидкости в ЭСП происходит более эффективно [1,2]. Большие капли, попадая в ЭСП, «взрываются» с образованием мелких, которые приобретают электрический заряд. При этом увеличивается напряженность поля, уменьшаются внутреннее давление капель и соответственно поверхностное натяжение.

С уменьшением внутреннего давления капель увеличивается эффективность воздействия аэродинамических сил, что резко повышает силу «взрыва» капель.

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

Наличие ЭСП, с одной стороны, будет приводить к росту скорости капель, а с другой — к деформации поля скоростей полидисперсной системы капель, вследствие чего усиливается возмущающее воздействие капель на поверхность пленки жидкости, стекающей по трубе, и интенсифицируется процесс теплоотдачи [ 1,2,4]. Благодаря воздействию ЭСП на режим охлаждения изменяется поверхностная энергия капель, что может привести к их охлаждению [8].

Для аппроксимации экспериментальных данных по псевдоламинарному режиму течения пленки предложена степенная зависимость [7]:

N11 = С Яе^ \¥е"' Ргл ва^ Оа‘<х Ргг

где N11 - число Нуссельта;

Яе - число Рейнольдса;

\¥е - число Вебера;

Гг — число Фруда; ва - число Галилея;

Рг - число Прандтля.

Яе = со ШV .

пл ж ' '

Скорость стекающей пленки определяется как

со = С./ко а

ПЛ ф Г ж пл4 тр '

где 8 - толщина стекаемой пленки.

пл

^ =/Г2 р / Кг -

где]х - плотность орошения, кг/(см2 с); х - определяющий размер;

с — коэффициент поверхностного натяжения заряженной капли;

/?о3 — радиус заряженной капли.

Для псевдотурбулентного режима [7]

Nu = CRe°'5We0’22 Fr°-75“^ Ga°’17 Ga0-17 Pr05 x(Pr /Pr )0-25,

где С = 5,1 при ф = 90°; С = 1,5 при ф = 65°. Зависимость (5) справедлива для условий:

10 < Ле< 545; 1,5107< \Уе<2,6105;

4,13 -104 < Рг < 11,8-104; 21 < Оад < 53;

4,6-108<ва <9,2210".

7 х 7

Значение коэффициентов ф, С, Яе, т, «для псевдолами-нарного режима течения пленки в уравнении (5) приведены в таблице; 1= —0,25; 7*= 0,33; г — 0,5; <? = 0,17.

ф Re min Re max С т п

90° 10 220 - 230 3,36-10-* 1,36 0

65° 10 165-360 5,МО-5 0,83 0,327

15° 44 425 4,3-Ю 0,15 0,725

Радиус распиливаемых капель в зависимости от потенциала, который подводится к штуцеру распылителя, при известной напряженности определяется по уравнению

'«.И-

Анализ моделей и принятых в них допущений дает основание предположить, что характер теплообмена в ЭСП будет подобным, если в качестве допущения принять, что процессы электролиза в стекающей пленке жидкости практически не влияют на теплообмен и гидродинамику течения. В этом случае ЭСП будет воздействовать только на гидродинамику течения СДЖ и стекающей пленки жидкости.

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

1. Наличие ЭСП увеличивает деформацию поля скоростей капель СДЖ, что следует из уравнения

Ж = е0 {Я 3 Е2[ 1+2

где - скорость капель без струи в ЭСП;

Е — напряженность поля; е — диэлектрическая постоянная жидкости;

Цг — коэффициент динамической вязкости газа, из которого видно, что с ростом 7?о3 капли ее скорость в ЭСП увеличивается, т.е. крупные капли будут еще быстрее «обгонять» мелкие.

2. Одноименный заряд капель снижает эффективность коагуляции и способствует диспергированию крупных капель [5].

3. Вектор напряженности ЭСП направлен по нормали к поверхности охлаждения. В этом случае скорость капель определяется достаточно просто по уравнению

IV =

к х с эк с эк ' ' 3 у/

где IVк— скорость капель;

V/ - скорость капель, результирующая в струе без ЭСП;

ф-угол, образованный осью струи и вектором напряженности ЭСП.

4. Рассматривается одномерная модель течения СДЖ с допущением асимметричности струи.

Наличие ЭСП будет способствовать уменьшению уноса капель жидкости за поверхность охлаждения, поэтому вводим коэффициент

^о)- Лоо I Л**)-

Плотность орошения при отсутствии ЭСП определяется как

где (7 — массовый расход воды через форсунку;

/ — орошаемая поверхность.

Зависимость плотности орошения в ЭСП _/ будет

иметь более сложный вид и зависеть от напряженности электростатического поля:

Е =£//Л,1п

где £7 — приложенное напряжение;

/?,— радиус электрода;

Я 2 — радиус трубы.

Задаваясь в первом приближении значением г] = 1,2...1,5, в дальнейшем уточним его величины, используя экспериментальные данные. Тогда плотность орошения в ЭСП будет определяться как

7э

Для определения режима стекания пленки (псевдо-ламинарного или псевдотурбулентного) вычисляется критическое число Рейнольдса [6]:

Яе = А^е°'73-10-7,

Крит 7 V '

где К- 0,31 при ф = 90°; К— 0,22 при ф = 65° и значение координаты X , соответствующее переходной границе,

И *ЧР /Л«-

Для проведения комплекса исследований по разработанным моделям процесса охлаждения нами создана программа для ЭВМ. Алгоритм расчетов построен таким образом, что позволяет проводить многовариантные расчеты (численное моделирование) путем задания определенного набора исходных данных и ограничений.

Влияние ЭСП на интенсивность теплоотдачи неоднозначно. Как показывает анализ полученной нами системы уравнений, критерии Фруда, Вебера и Галилея, входящие в уравнение (1) и (5), имеют различный ха-

Рис. 1. Зависимость критерия Нуссельта от угла распыла охлаждающей жидкости

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

рактер зависимости от напряженности ЭСП и первоначального диаметра распыливаемых капель. Как и следовало ожидать, наиболее «чувствительны» к действию ЭСП мелкие капли радиусом менее 10-5м. С ростом напряженности ЭСП в большей степени начинают проявляться эффекты электростатического диспергирования [8], связанные с уменьшением сил поверхностного натяжения жидкости.

На рис. 1 представлены зависимости критерия Нус-сельта от радиуса капель /?о и угла ср. Наибольшее влияние размера капель сказывается на больших углах атаки. Так, при ф =60. ..90° и уменьшении радиуса от Ло = 10~4 до Я о = 10~5 м критерий Нуссельта возрастает более чем в 2 раза. При меньших значениях ф,

а, Вт /

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

1 — в ЭСП; 2 — испарительное охлаждение;

3 — водой; 4 — воздухом

т.е. когда ось струи становится параллельной оси трубы, изменение радиуса капель практически слабо влияет на эффективность теплоотдачи. Это можно объяснить тем, что при «атаке» струи по нормали к поверхности капли оказывают существенное возмущающее воздействие на поверхность стекающей пленки, увеличивая эффективность теплоотдачи.

На рис. 2 представлены зависимости коэффициента теплоотдачи от скорости охлаждающей среды и напряженности электрического поля. Наибольшего значения а достигает при скорости охлаждающей среды v = 1,6...1,8 м/с и напряженности поля Е — 30 кВ/м.

На рис. 3 представлены результаты теоретического исследования зависимостей коэффициентов теплоотдачи от скорости воздуха для различных методов охлаждения.

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

Предложенный алгоритм расчета может быть применен: при разработке и конструировании конденсаторов испарительного охлаждения; в оросительных камерах систем кондиционирования воздуха; в градирнях охлаждения оборотной воды.

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

1. Абросимов Ю.Г., Исаченко В.П., Солодов А.П. Об устойчивости поверхности жидкости в электрическом поле //Труды МЭИ. 1971.Вып.81.

2. Абросимов Ю.Г., Исаченко В.П., Солодов А.П. Неустойчивость поверхности жидкости в электрическом поле//Докл. научно-техн. конф. по итогам НИР за 1968 - 1969 гг., МЭИ, 1969.

3. Амелин А.Г. Теоретические основы образования тумана при конденсации пара.-М.: Химия, 1972.

4. Бабакин Б.С. Электротехнология в холодильной промышленности.-М.: Агропромиздат, 1990.

5 Бейтуганов М.Н., Белгороков Л.Г. Изучение взаимодействия капель в электрическом поле //Атмосферное электричество: Тр. III Всесоюзн.симпозиума.-Л.: Гидрометеоиз-дат, 1988.

6. Гогонин И.И., Дорохов А. Р., Пирогова О.Ю. Теплообмен на струях в пространстве между трубами горизонтального пучка конденсатора//Тепломассообмен-ММФ-96. Т. 4,

ч. 2.-Минск, 1996.

7. Исаченко В. П. и др. Струйное охлаждение. -М.: Энергоатомиз-дат, 1987.

8. Пажи Д.Г., Галустов B.C. Основы техники распыления жидкостей.-М.: Химия, 1984.