Расчетные исследования пневматической системы стабилизации температуры

УДК 536.74; 62-522.7

А. А. Крутиков, К. Е. Д е м и х о в, А. В. Чернышев, О. В. Белова, И. А. Смирнов

РАСЧЕТНЫЕ ИССЛЕДОВАНИЯ ПНЕВМАТИЧЕСКОЙ СИСТЕМЫ СТАБИЛИЗАЦИИ ТЕМПЕРАТУРЫ

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

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

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

Ьенжвятвв

\ жидкость

Рис. 1. Схема пневматической системы термостабилизации

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

Рассматриваемое устройство представляет собой открытую термодинамическую систему, границы которой проницаемы как для потока вещества, так и для теплового взаимодействия, и может рассматриваться как объект с сосредоточенными параметрами или как объект с распределенными параметрами. Математическое описание тепловых процессов в приближении объекта с сосредоточенными параметрами не позволяет получить распределение физических величин по объему системы. Следовательно, невозможно отследить неравномерность теплового поля объекта термостабилизации, что является существенным недостатком этого метода. Для получения температурных полей необходимо рассматривать устройство в приближении объекта с распределенными параметрами. Такое приближение более точно описывает реальные физические процессы, в том числе переходные процессы, протекающие в устройствах.

Постановка задачи. Принципиальная схема пневматической системы стабилизации температуры включает в себя следующие основные элементы: источник теплоты, объект теплового воздействия и рабочее тело. В исследуемой системе стабилизации температуры источником теплоты служит электрический нагреватель, а объектом теплового воздействия — жидкость, которая расположена в сосуде объемом 0,2 мл. В качестве рабочего тела в пневматических системах используется воздух окружающей среды.

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

Рис.2. Расчетная область П:

$1 и 52 — вход и выход из расчетной области; 53 — внешние границы расчетной области; $4 и $ — границы раздела между жидкой средой и твердым телом, а также воздухом в сосуде

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

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

В итоге расчетная область будет состоять из газообразного рабочего тела Пь жидкости П2, воздуха П3, сосуда П4, стенок П5 и электрического нагревателя П6.

Математическая модель. При построении математической модели приняты следующие допущения: жидкая среда считается ньютоновской и несжимаемой; движение жидкой среды ламинарное; все рассматриваемые материалы однородны по составу и изотропны по теплофизическим свойствам; физические параметры жидкой среды

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

Система уравнений для описания процессов гидродинамики и теплообмена жидкой среды включает в себя [1-4]:

уравнения движения для жидкой среды

д dp

—

d dp

—

d dp

—

уравнение неразрывности для жидкой среды

f + div

уравнение энергии для жидкой среды

d

—

где p — плотность, кг/м3; p — давление, Па; ц — динамический коэффициент вязкости, Па-с; c — теплоемкость при постоянном давлении, Дж/(кг-К); A — коэффициент теплопроводности, Вт/(м-К); и — вектор скорости; u,v,w — скорость в проекциях на оси координат x, y и z, м/с; T — температура, К; t — время, с; SZi — член, учитывающий массовые силы, которые действуют в направлении оси z, Н/м3.

В области Q (рабочее тело) не учитывается свободная конвекция теплоты и, следовательно, Sz1 = 0. При расчете областей Q2 (жидкость) и Q3 (воздух) учитывается свободная конвекция теплоты в приближении Буссинеска [1-4]:

Szi = -poiÇPTi

где p0i — плотность жидкой среды при температуре T0i; fîTi — температурный коэффициент объемного расширения жидкой среды.

При расчете областей Q2 и Q3 под давлением, обозначаемым символом "p", понимается давление жидкой среды за вычетом гидростатического давления.

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

теплопроводности для твердого тела получается из уравнения энергии для жидкой среды при условии, что и = 0 [1-4]:

д

—

где STi — член, учитывающий количество теплоты, выделяемое внутренними источниками в единице объема среды в единицу времени, Вт/м3.

В области П 6 (электрический нагреватель) происходит выделение теплоты вследствие протекания электрического тока и, следовательно, ST 6 = Q.B областях П 4 (сосуд) и П 5 (стенки) отсутствуют внутренние источники теплоты и, следовательно, ST4 = ST5 = 0.

Областью определения системы уравнений (1)-(7) является область П (см. рис. 2), состоящая из подобластей П ^ (i = 1, 2,..., 6) и определяемая положением радиус-вектора в декартовой системе координат r (x,y,z,t) £ П или в цилиндрической системе координат R (r,ö,z,t) £ П.

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

u

где (R) £ Si.

На выходе S2 из расчетной области П задано нулевое изменение скорости и нулевое изменение температуры рабочего тела в направлении выхода:

du

где (R) £ S2.

На внешних границах S3 расчетной области П задано нулевое значение скорости рабочего тела для уравнений движения и адиабатическое условие теплообмена для уравнения энергии:

dT

u

где (r) £ S3.

На границах раздела S4 между жидкой средой и твердым телом задано нулевое значение скорости жидкой среды для уравнений движения:

u

где (r) £ S4.

На границе раздела S5 между жидкостью П2 и воздухом П3 в сосуде задано условие непроницаемости для уравнений движения и условие

идеального теплового контакта для уравнения энергии:

w

где (г) е S5.

Начальные условия. В начальный момент времени t = 0 жидкая среда считается неподвижной, т.е. и = 0. Также в начальный момент времени задано равномерное распределение температуры по объему жидкой среды и твердого тела Т = 293 К.

Условием окончания счета является достижение времени расчета заданного значения t = 60 с.

Решение. Представленная система дифференциальных уравнений (1)-(7) с граничными условиями (8)-(12) решается численными методами. Для получения дискретного аналога используется метод контрольного объема в конечно-элементной постановке. Методы составления дискретных аналогов и методы решения полученной системы алгебраических уравнений приведены в работах [5-7].

Исходные данные. Физические свойства материалов при средней температуре Т = 333 К [8,9] сведены в таблицу.

Таблица

Физические свойства используемых материалов

Область Материал Свойства материалов

р, кг/м3 ßT,1/K Л, Вт/(м-К) с, Дж/(кг-К) р, Па-с

П1 Воздух 1,060 - 0,029 1005 20,1 • 10-6

П2 Вода 998,2 4,98 • 10-4 0,659 4179 469,9 • 10-6

П3 Воздух 1,205 2,69 • 10-3 0,029 1005 20,1 • 10-6

П4 Полипропилен 900 — 0,2 1927 —

П5 Алюминиевый сплав Д16Т 2700 227 920 —

Пб Нихром 8800 — 12,6 448 —

Объемный расход воздуха в системе составляет G = 0,1 • 10-3 м3/с. Площадь входного сечения S = 0,39 • 10-3 м2. Скорость рабочего тела на входе в расчетную область и = G/S = 0,26 м/c. Температура рабочего тела на входе в расчетную область T = 293 K.

Электрическая мощность, подводимая к нагревателю, P = 10 Вт. Объем материала нагревателя V = 3,76 • 10-8 м3 (общая длина проволоки нагревателя 0,5 м, диаметр проволоки нагревателя 0,3 • 10-3 м). Количество теплоты, выделяемое в нагревателе в единице объема за единицу времени, Q = P/V = 2,66 • 108 Вт/м3.

Рис. 3. Сетка расчетной области

Толщина стенки сосуда 0,4 • 10-3 м, толщина металлических стенок 0,3 • 10-3 м. Объем жидкости в сосуде 22 • 10-9 м3.

В результате разбиения расчетной области на конечные элементы получена расчетная сетка (рис. 3). Число элементарных ячеек составляет 1043 472, число узлов 1081923. Расчет проводился с шагом 0,01 с, что составило 6000 шагов по времени.

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

Результат численного исследования теплового поля жидкости и рабочего тела в процессе нагрева показан на рис. 5, из которого видно, что во время прогрева температура рабочего тела значительно отличается от температуры жидкости в сосуде. Следовательно, во время переходных процессов будет затруднительно регулировать температуру жидкости по температуре рабочего тела.

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

На базе созданной математической модели проведены численные исследования нестационарных тепловых процессов в элементах системы.

Рис.4. Распределение температур (а) и скоростей (б) в расчетной области на 30-й секунде прогрева

4 \

/ С

/ \ \ \

\ 2

/ /V \ 3

и 1

10 15 20

25 30 35 Время, с

40 45 50 55 60

Рис. 5. Результаты численного исследования поля температур в процессе прогрева:

1 и 2 — температура зоны максимального и минимального прогрева жидкости; 3 — среднемассовая температура жидкости; 4 — температура рабочего тела в непосредственной близости от сосуда

Результаты численного исследования демонстрируют, что во время переходных процессов скорость нагрева рабочего тела достигает 10 K/с, в то время как скорость нагрева реакционной смеси не превышает 2 K/с. В то же время разброс температур в микрообъеме реакционной смеси (22мкл) достигает 7K.

Работа выполнена при частичной поддержке Госконтракта № 02.467.11.3005 с Федеральным агентством по науке и инновациям.

СПИСОК ЛИТЕРАТУРЫ

1. Исаченко В. П., Осипова В. А., Сукомел А. С. Теплопередача. - М.: Энергоиздат, 1981. - 416 с.

2. Самарский А. А., Вабищевич П. Н. Вычислительная теплопередача. - М.: Едиториал УРСС, 2003. - 784 с.

3. Себеси Т., Брэдшоу П. Конвективный теплообмен / Пер. с англ. - М.: Мир, 1987. - 592 с.

4. Лойцянский Л. Г. Механика жидкости и газа: Учеб. для вузов. - М.: Дрофа, 2003. - 840 с.

5. Чернышев А. В., Белова О. В. Метод решения сопряженной задачи конвективного теплообмена на примере термостатирующего устройства // Вестник МГТУ им. Н.Э. Баумана. Сер. "Машиностроение". - 1998. - № 4. -C. 77-87.

6. Кюрджиев Ю. В. Моделирование рабочих процессов, разработка и модернизация пневматических систем и агрегатов с учетом образования конденсата рабочего тела // Автореф. дис.... канд. техн. наук. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. - 163 с.

7. П а т а н к а р С. В. Численное решение задач теплопроводности и конвективного теплообмена при течении в каналах / Пер. с англ. Е.В. Калабина; Под ред. Г.Г. Янькова. - М.: Изд-во МЭИ, 2003. - 312 с.

8. Теплотехнический справочник. В 2 т. / Под общ. ред. В.Н. Юренева, П.Д.Лебедева. - М.: Энергия, 1975. - 1976.

9. Физические величины: Справочник / А.П.Бичев, Н.А.Бабушкина, А.М. Братковский и др. / Под ред. И.С. Григорьева, Е.З. Михайлова. - М.: Энер-гоатомиздат, 1991. - 1232 с.

Статья поступила в редакцию 26.06.2007

Алексей Александрович Крутиков родился в 1981 г., окончил МГТУ им. Н.Э. Баумана в 2004 г. Аспирант кафедры "Вакуумная и компрессорная техника" МГТУ им. Н.Э. Баумана. Специализируется в области вычислительной гидродинамики и теплообмена.

A.A. Krutikov (b. 1981) graduated from the Bauman Moscow State Technical University in 2004. Post-graduate of "Vacuum and Compressor Technology" department of the Bauman Moscow State Technical University. Specializes in the field of computing hydrodynamics and heat exchange.

Константин Евгеньевич Демихов — д-р. техн. наук, профессор, первый проректор по научной работе МГТУ им. Н.Э. Баумана, зав. кафедрой "Вакуумная и компрессорная техника". Автор более 150 научных трудов в области вакуумной и компрессорной техники.

K.Ye. Demikhov — D. Sc. (Eng.), professor, first pro-rector for scientific work, head of "Vacuum and Compressor Technology" department of the Bauman Moscow State Technical University. Author of more than 150 publications in the field of vacuum and compressor technology.

Андрей Владимирович Чернышев родился в 1954 г., окончил МВТУ им. Н.Э. Баумана в 1977 г. Канд. техн. наук, доцент кафедры "Вакуумная и компрессорная техника" МГТУ им. Н.Э. Баумана. Автор ряда научных работ в области теплофизики, теории и расчета пневматических систем.

A.V. Chernyshov (b. 1954) graduated from the Bauman Moscow Higher Technical School in 1977. Ph. D. (Eng.), assoc. professor of "Vacuum and Compressor Technology" department of the Bauman Moscow State Technical University. Author of some publications in the field of thermal physics, theory and design of pneumatic systems.

Ольга Владимировна Белова родилась в 1971 г., окончила МГТУ им. Н.Э. Баумана в 1994 г. Канд. техн. наук, доцент кафедры "Вакуумная и компрессорная техника" МГТУ им. Н.Э. Баумана. Автор ряда научных работ в области вычислительной гидродинамики и теплообмена.

0.V. Belova (b. 1971) graduated from the Bauman Moscow State Technical University in 1994. Ph. D. (Eng.), assoc. professor of "Vacuum and Compressor Technology" department of the Bauman Moscow State Technical University. Author of some publications in the field of computing hydrodynamics and heat exchange.

Игорь Алексеевич Смирнов родился в 1947 г. Зав. лабораторией ГНЦ РФ — ИМБП РАН, канд. техн. наук. Автор ряда научных работ в области теоретических и экспериментальных исследований искусственной газовой среды.

1.A. Smirnov (b. 1947) - head of laboratory of "IMBP" of RAS, Ph. D. (Eng.). Author of some publications in the field of theoretical and experimental study of artificial gas medium.

ЖУРНАЛ "ВЕСТНИК МОСКОВСКОГО ГОСУДАРСТВЕННОГО ТЕХНИЧЕСКОГО УНИВЕРСИТЕТА имени Н.Э. БАУМАНА"

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

Журнал "Вестник МГТУ им. Н.Э. Баумана" в соответствии с постановлением Высшей аттестационной комиссии Федерального агентства по образованию Российской Федерации включен в перечень периодических и научно-технических изданий, в которых рекомендуется публикация основных результатов диссертаций на соискание ученой степени доктора наук.

Главный редактор журнала "Вестник МГТУ им. Н.Э. Баумана" — ректор МГТУ им. Н.Э. Баумана, чл.-кор. РАН, д-р техн. наук, профессор И.Б. Федоров.

Журнал издается в трех сериях: "Приборостроение", "Машиностроение", "Естественные науки" с периодичностью 12 номеров в год.

Подписку на журнал "Вестник МГТУ им. Н.Э. Баумана" можно оформить через агентство "Роспечать".

Подписка по каталогу "Газеты, журналы" агентства "Роспечать"

Индекс Наименование серии Объем выпуска Подписная цена

Полугодие 3 мес. 6 мес.

72781 "Машиностроение" 2 250 500

72783 "Приборостроение" 2 250 500

79982 "Естественные науки" 2 250 500

Адрес редакции журнала "Вестник МГТУ им. Н.Э. Баумана": 105005, Москва, 2-я Бауманская ул., д. 5.

Тел.: (495) 263-62-60; 263-60-45. Факс: (495) 261-45-97. E-mail: press@bmstu.ru