Статья: Модель системы охлаждения индукционного устройства

ФГОУ ВО Сибирский федеральный университет

Модель системы охлаждения индукционного устройства

Кинев Е.С., Тяпин А.А., Гришко Г.С.

Аннотация

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

Ключевые слова: индукционное устройство, система охлаждения, испаритель, конденсатор, теплопередача, тепловой поток.

Induction device cooling system model

Abstract

The paper deals with the calculation model of the cooling system of an induction device with power windings, characterized by increased efficiency. Evaporative-condensation cooling system is used as a base for building the model. During the simulation, the authors investigated the properties of three types of coolant; the characteristic features that determine their properties with respect to the ability of increased heat transfer were established. The study of the thermal modes of the model made it possible to determine the nature of the influence of the depth and width of the condenser, as well as the air flow rate on the amount of the extracted heat flux for the used working fluids.

Keywords: induction device, cooling system, evaporator, condenser, heat transfer, heat flow.

Введение

Специфика индукционного оборудования заключается в относительно невысокой тепловой эффективности, поскольку тепловой КПД индукторов, предназначенных для перемешивания расплавленных металлов не превышает значений 20-25%. Для повышения нагрузки нередко применяют обмотки из полой медной трубки, с проточным водяным охлаждением. Другим способом охлаждения обмоток является принудительный обдув, который нередко является предпочтительным по причине безопасности производства, особенно при использовании индукторов в МГД-устройствах на печах с жидким алюминием [1]. Однако интенсивность воздушного охлаждения весьма мала и не позволяет существенно увеличивать плотность тока в обмотках. Аналогичными свойствами обладают индукторы, предназначенные для экструзии алюминия. В связи с изложенным, поиск и разработка новых эффективных систем охлаждения технологических машин, является актуальной задачей.

С увеличением нагруженности оборудования в технологическом процессе увеличивается выделение тепла, которое нужно рассеивать. Концептуально любая система охлаждения построена на одном или нескольких способах переноса теплоты [2], [3]. Отвод тепла во внешнюю среду происходит путем теплопроводности, излучения, естественной или принудительной конвекции, а так же кипения-конденсации [4]. Развитие индукционного оборудования идет в направлении повышения удельных энергетических параметров при сохранении размеров и массы [5]. Это требует от системы охлаждения максимальной производительности при минимальных габаритах. Причем под производительностью, нередко, понимают отводимую от объекта мощность при возможно меньшем температурном напоре.

Среди существующих средств охлаждения, требованию повышенной эффективности теплоотвода в большей степени отвечает испарительно-конденсационная система [6], в основе которой лежит испарение (кипение) жидкости на охлаждаемой поверхности с поглощением тепла и последующая конденсация пара в теплообменнике с отдачей тепла промежуточному теплоносителю. Однако, у проектировщика, как правило, отсутствует информация о характере влияния основных параметров системы на ее эффективность в зависимости от типа рабочей жидкости [7]. При разработке подобных систем охлаждения это является серьезной проблемой, поэтому оценка рабочих свойств жидкостей имеет в себе и практическое значение. Другие аспекты технического проектирования индукторов с испарительно-конденсационной системой охлаждения здесь не рассмотрены.

Постановка задачи

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

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

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

Допущения

В расчетной модели принята следующая система допущений и ограничений: нестационарные процессы нагрева и охлаждения отсутствуют; теплообмен с окружающей средой происходит в теплообменнике; температура кипения сравниваемых теплоносителей задана при нормальном атмосферном давлении; не конденсирующиеся газы в системе отсутствуют; конденсатор установлен вертикально (j_с = 0 рад); влажность окружающего воздуха равна нулю; отвод тепла от нагреваемых поверхностей происходит при кипении рабочей жидкости в большом объеме.

Решение

охлаждение индукционный испарительный конденсационный

Многофазные среды и течения, теории их расчета и моделирования в основном рассматриваются применительно к паропроизводящим установкам (паровые котлы, ядерные реакторы и т. д.). Достаточно большой объем исследований в этой области накоплен в нашей стране и за рубежом [8], [9]. Отдельным направлением является разработка теоретических основ проектирования и расчета различных типов теплообменных аппаратов [10]. На их основе существующих теорий авторами предложен подход, позволяющий оценить эффективность многофазной системы в качестве системы охлаждения технологического оборудования.

Цикл работы рассматриваемой системы охлаждения выглядит следующим образом:

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

2. Транспорт пара к конденсатору (теплообменнику) предусмотрен с помощью трубопроводов.

3. Конденсация пара на поверхности конденсатора-теплообменника и передача тепла во внешнюю среду или промежуточному теплоносителю.

4. Возврат конденсата из теплообменника обратно в бак (свободная циркуляция теплоносителя).

5. Обеспечение чистоты и требуемого удельного электрического сопротивления теплоносителя (например, деионизация, если в качестве теплоносителя используют воду).

При моделировании теплового режима индукционного устройства, в качестве исходных параметров приняты следующие: выделяемая тепловая мощность в катушках - Q_kt = 100 кВт; площадь тепловыделяющих элементов S_kt = 1,472 м²; давление на линии насыщения - p_r = 101325 Па; температура воздуха на входе в конденсатор - t_sr0 = 293,15 К; тип поверхности теплообмена в конденсаторе - пучки плоских оребрённых труб со сплошными гофрированными ребрами по типу ПлСР-3 [6]; высота конденсатора l_c = 0,57м; глубина конденсатора Y_c = 0,045 м; ширина конденсатора X_c = 0,395 м; расход воздуха вентилятором w_v = 24000 м³/ч.

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

Рис. 1 - Эскиз расчетной области для построения модели

На рисунке приняты следующие обозначения. Герметичная немагнитная емкость с диэлектрическим теплоносителем, заполненная до уровня H, обозначена 1. В теплоноситель полностью погружены силовые обмотки индуктора 2, имеющие площадь охлаждаемой поверхности S_kt. Для фазового превращения теплоносителя предназначен бокс конденсатора 3 в виде параллелепипеда с известными размерами, сообщающийся с основной емкостью парожидкостными технологическими патрубками 4.

Удельный тепловой поток на поверхности катушек, Вт/м²:

(1)

где Q_kt - величина тепловой мощности, выделяемой на катушках, Вт; S_kt - площадь поверхности тепловыделяющих элементов, контактирующая с жидкостью, м².

Расчетное значение коэффициента теплоотдачи для воды (Вт/м²ЧК) проводят по известной формуле [6]:

(2)

где p_r - давление на линии насыщения, МПа.

Для модели системы охлаждения с органическими рабочими жидкостями (спирт, фреон) коэффициент теплоотдачи может быть определен с использованием эмпирического выражения [7+5]:

(3)

где r - плотность жидкости, кг/м³; r_p - плотность пара, кг/м³; л - коэффициент теплопроводности рабочей жидкости, Вт/мЧК; н - коэффициент кинематической вязкости жидкости, м²/с; у - поверхностное натяжение жидкости, Н/м; T_s - температура на линии насыщения (кипения), К.

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

(4)

где r - теплота парообразования, Дж/кг; ц_с - угол наклона бокса конденсатора к горизонту, рад; l_c - высота конденсатора, м; Дt - температурный напор между паром и стенкой конденсатора, К; м - коэффициент динамической вязкости жидкости, ПаЧс.

В качестве промежуточного теплообменника 3 (рис. 1) для сброса тепла, выделяющегося на катушках, в окружающую среду выбран компактный водо-воздушный теплообменник с развитой поверхностью теплообмена. Для определения параметров теплообменника и процесса передачи тепла промежуточному теплоносителю в расчетной модели была использована методика, описанная в [12], [13] с адаптацией к условиям работы рассматриваемой системы охлаждения. Теплоотдача и гидравлическое сопротивление теплообменника описано экспериментальными зависимостями вида [13]:

(5)

(6)

где f_voz - коэффициент сопротивления потоку; - безразмерный параметр, характеризующий теплоотдачу и представляющий собой произведение числа Стэнтона на число Прандтля в степени 2/3.

В рассматриваемой модели зависимости (5) и (6) с помощью известных методов регрессионного анализа представлены в виде функции логарифмического вида:

(7)

(8)

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

(9)

где ДP_pod - потери давления на подъемном участке, Па; ДP_c - потери давления в конденсаторе, Па; ДP_op - потери давления в опускном участке, Па, g - ускорение свободного падения, м/с².

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

(10)

где t_sr1 - температура воздуха на выходе, К; t_sr0 - температура воздуха на входе, К.

В модели основные теплофизические свойства воздуха - коэффициент динамической вязкости, теплоемкость, теплопроводность и плотность определены в зависимости от температуры. Теплоотдача к воздуху задана критерием Стэнтона:

(11)

Коэффициент теплоотдачи со стороны воздуха, Вт/м² К:

(12)

где G_v - массовая скорость потока воздуха, кг/сЧм²; c_pv - теплоемкость воздуха, Дж/кгЧК.

На интегральный коэффициент теплопередачи существенное влияние оказывает оребрение поверхности со стороны воздуха. Параметры эффективности оребрения включены в модель с помощью общепринятых зависимостей. Интегральный коэффициент теплопередачи (отнесенный к теплопередающей поверхности со стороны воздуха), Вт/м²ЧК:

(13)

где м₀ - общая эффективность развитой поверхности на стороне воздуха; Ш_cw - отношение поверхности теплоотдачи на стороне конденсата к полному объему, м²/м³; Ш_c - отношение поверхности теплообмена на стороне воздуха к полному объему, м²/м³.

Тепловая мощность при отводе тепла от конденсатора в окружающий воздух, Вт:

(14)

где V_c - объем конденсатора, м³.