CC BY
22
УДК 664.7
Математическая модель процесса намораживания льда на сферической поверхности применительно для аккумуляторов холода
Д-р техн. наук И. Е. ЛОБАНОВ Московский авиационный институт 125993, г. Москва, Волоколамское шоссе, 4 Д-р техн. наук Б. С. БАБАКИН, Р. Б. АЙТИКЕЕВ, канд. техн. наук М. И. ВОРОНИН, С. Б. БАБАКИН [email protected]
Московский государственный университет пищевых производств 109316, г. Москва, ул. Талалихина, 33
Предложенное аналитическое решение квазистационарной задачи о намораживании на внешней сферической поверхности с граничными условиями третьего рода на внутренней поверхности и граничными условиями третьего рода на внешней поверхности, более общее, чем существующие, которое необходимо использовать для решения задачи о намораживании льда на рабочих поверхностях аккумуляторов холода.
Ключевые слова: воздухоохладитель, моделирование, аккумуляторы холода, охлаждаемая поверхность, тепло-и массообмен, фазовый переход, сферическая поверхность.
Mathematical model of the ice build-up process on a spherical surface for cold accumulators
D. S^ I. E. LOBANOV
Moscow Aviation Institute 125993, Russia, Moscow, Volokolamskoe Shosse, 4 D. Sc. B. S. BABAKIN, R. B. AITIKEEV, Ph. D. M. I. VORONIN, S. B. BABAKIN
Moscow National University of Food Production 125080, Moscow, Volokolamskoye shosse, 11
An analytical solution is suggested for the quasi-stationary problem of ice build-up on an outer spherical surface with third-order boundary conditions on the inner surface and third-order boundary conditions on the outer surface. More general than the existing ones, this method should be applied to solve the task of ice formation on the working elements of cold accumulators.
Keywords: air cooler, modelling, cold accumulators, cooled surface, heat and mass transfer, phase transition, spherical surface.
В рамках данной работы рассматривается задача затвердевания, т. е. переход вещества из жидкого состояния в твердое условно рассматривается без учета особенностей формирования кристаллов, что следует использовать для решения задачи о намораживании на аккумуляторах холода. Изучается проблема только с учетом тепловых взаимодействий между твердой и жидкой фазами и телами, которые находятся в тепло-
вом контакте между ними. Поскольку решения частных случаев данной задачи, представленные ранее, были получены независимым образом, авторами ставится задача получения решения общего вида для намораживании на сферической поверхности, как на внутренней, так и на внешней, для граничных условий первого и третьего родов. Процесс намораживания на внешней поверхности сферы (см. рис.) описывается следующими дифференциальными уравнениями одномерной нестационарной теплопроводности в сферических координатах
ciPi = 1 1
'эч
Ът1
2 dt
Ът
"т0 < т < т0 + R ;
ср ât =À
dr2
"r,. + R < r < r
2 dt_ r dr R + X
(1)
(2)
где т — время; с — теплоемкость намороженного слоя; р — плотность намороженного слоя; X — коэффициент теплопроводности намороженного слоя; t — температура намороженного слоя; с'1 — теплоемкость материала стенки; р1 — плотность материала стенки; Х1 — коэффициент теплопроводности материала стенки; t1 — температура стенки; г — текущая радиальная координата, отсчитываемая от центра сферы; г0 — внутренний радиус полой сферы; R — радиальная толщина сферы; £ — радиальная толщина намороженного слоя.
Граничные условия:
= t
л dt, I . .
1 gT I r = * = «0( tn " t0);
(3)
1 Г кр
к
Схематическое изображение модели для решения задачи о намораживании внешней сферической поверхности с граничными условиями IIIрода на внутренней поверхности
дг'
Ч |т=г0 + Я+Х 1 кр'
а X
=г0 + К+%= «ж «ж - + Рг — >
Ч г=г„ + Я Цг=Г0 + Я
1 I _ 1 ^ I
1 э7 I - +* _ 1 ^ I - _*+* •
Стационарное распределение температуры в полой сфере при одномерной теплопроводности выглядит следующим образом:
С
1 = СХ--2-,\/г0+Я<г<г0+11 + Х;
tl = С3 -^ "г0 < г < г0 + R г
(9) (10)
где С, i = 1, 2, 3, 4 — константы.
Вышеуказанные константы находятся из граничных условий и условий сопряжения. Сначала необходимо составить систему уравнений:
1 I Г= 1 = «о( *п - {о) Ю11 С=а0 ( *п - Ч); (11)
дг
Ч г = г0 + Я ¿гр Ю ¿гр С 3
С
Г, + Я
г=г0 + Я ¿гр Ю ¿гр С1
С
го + Я
= г ю г = С,-
кр 7 кр 1
С
г=г0 + К+Х кр кр
Г + К + Х
(12)
(13)
(14)
(5)
(6)
Вышеуказанные константы получаются из решения следующей системы уравнений:
, = С -- С ■
'гр ^ 1
где / — криоскопическая температура; аж — коэффициент теплоотдачи со стороны жидкости; ¿ж — температура жидкости; г — удельная теплота фазового перехода; /0 — температура холодильного агента; а0 — коэффициент теплоотдачи со стороны холодильного агента; / — температура стенки со стороны холодильного агента.
Условия сопряжения:
= ^(Х);
г = с, -
кр 1
а
го + Я
с
Г, + * + Х
(15)
1 -2-=а о (ч - О;
с
г = С3 - с
гр 3
Ч + *
(7)
(8)
В развернутой форме константы представлены в виде
(16)
Ч _
С2 =
X
г0 X + г1 + 2 г0 Я + ЯХ +Я2
г -
кр
Задача намораживания решается при граничных условиях третьего рода на внутренней поверхности (г = г0) и граничных условиях третьего рода на внешней поверхности сферы (г = г0 + R + £), поскольку эти условия являются более общими, чем рассмотренные ранее.
В данном исследовании основной целью является отыскание именно аналитических решений. Вышеуказанный подход, базирующийся на методе замещения, ранее был применён к решению различного типа задач и является вполне обоснованным.
Следовательно, основным допущением, принятым при решении задачи замораживания на сферической стенке является то, что температура во всех промежутках — [г; г0 + R] и [г; г0 + R + £] — замещается температурами, получающимися при стационарной одномерной теплопроводности в полой сфере.
В дальнейшем, решение задач намораживания получается без дополнительных допущений и, в рамках уже принятых допущений, может считаться точным аналитическим решением.
ДХ+г02 + 2 г0 Я + г0Х+Я2
¿гр'
с = Угр Г0 +1 ¿гр Х + «0 Г0 ¿п - «0 Г0 ¿0
3 1( го + X)
г2
(17)
(18) (19)
Таким образом, стационарное температурное поле в полой сфере будет выжаться следующим образом:
, Уо+'кр * + К, Х - ^гр _
X
Чр ^ОХ + Тр К Х+ ^кр^О2 +2 УкрД + Ч*кр X _ гХ
- (гр ч2-2^гр я+^кр II2+ЩрЛ - ЧрЯ2
гХ '
"г0 + Я < г < г0 + Я + Х; (20)
о
о
Х
4 =
^ гр^О 1 гр п 0^0 0 ,, ^ Г02
У
М'-о+Л)
Далее следует найти температуру tп из условия
= ^:
— ... ЦЛ'п Ч) л '
Решив последнее уравнение, получим
К
1 Лгр Г0 + 1 Л/ + 1 0а 0Г0Я
г0 + 11 Я + а 0 г0 Я
В дальнейшем следует найти температуру t из условия сопряжения (7). Производные температур равны:
Э?!
дг соп
дг
дг
дг
(
-а0 Ч
К
(го + я + Х).
%г2
Эти же производные в точке сопряжения равны:
+я = -( го+*)( К - 0('0
У./о+У^+^ад*
Кго + \К + аогоК
(г+я+Х)
,2 '
Эг
Э?! I
г=г„ + я /
2 '
= -«о
'о-
ч
Уц/р +
\г0+Х1Я + аХ)г0Я
Условие сопряжения (7) для нахождения температуры t :
£,(г0+К)2
= -«о
ч-
Угр^о +А.1^,Я + Цх0г0Д Ус+М+ <¥•(>* .
Решение (28) дает температуру t :
С =
¿кр Ч,2 Я + «Л ¿крЧЧ Д + г0 ДО ! + 2а 0 г0 + 2 ^ Яг0 \Я + ^1X0, И- + ДО ^ + ^ 2Я1 1г021 1 +а 0 1г02а 0Я + Г01Х11 + 1Я 2 о0 г0 + 1 г0Я,Я + г01£оо Я +1X1 Л + 1Я211
(21) (22)
(23)
(24)
(25)
(26)
(27)
(28) (29)
Продолжительность процесса намораживания тн слоя толщиной £ детерминируется из решения дифференциального уравнения (6) после подстановки в него значений соответствующей производной и значения температуры ^ из выражения (29):
Оо + Я)(7гр-гкр)
^ | г=г0+*+5 - -Чг0 + «)<Лр ~ 4р ) ^
г=г0+Д+5
(30)
После преобразований, решение вышеуказанного дифференциального уравнения относительно продолжительности процесса намораживания тн слоя толщиной £ будет выражаться следующим интегралом:
(■х а X _
(31)
т„ = -
РГ /о
+Ф (с - ОАа
[(г0П 1 + а011 г02 + г^ао Л + ХХ1Я )Х+^02111 + 211 г0 Я + 1а 0 г02 Я + г^ а0 Л 2 + 1Л1Тг 2) ]( г0 + Я + X)
-а (г - г )
ж V ж кр /
Запишем окончательный интеграл (31) в каноническом виде:
гх ¿X
"РГ То
К1
к2-х тКз -х+К4
где
К1 = (го + Я)(4р - и) 11аго2; К2 = го 111 + ао 11 го2 + го 1ао Я + 111 Л;
Кз = (го + Л)(2го 1ао Я + 2111Л + 2го 111 + ао 11 го2); К4 = 1 (го + Я)2 (11 го + 11Я + ао го Я); К5 = - аД4-4р).
После проведения интегрирования и соответствующих преобразований, получим окончательное выражение относительно времени намораживания т слоя толщиной 4 в каноническом виде:
t„ =
рг
К
-X+2KJ0
где
0
arctg(K5(2K2X+K3^) -arctg( К,К30 1
(33)
+4 KJLJt, - K^IQ'
Данное решение является канонической и компактной записью аналогичного решения, полученного ранее, но в отличие от последних, как показывает анализ, имеет обобщенный вид, к которому могут быть приведены аналогичные решения для различных граничных условий и поверхностей намораживания (как для внешней, так и для внутренней сферы). Как видно из окончательных выражений, решения задачи внешнего намораживания при граничных условиях третьего рода на внутренней поверхности и граничных условиях третьего рода на внешней поверхности сферы и при граничных условиях первого рода на внутренней поверхности и граничных условиях третьего рода на внешней поверхности различаются только комплексами К.
Анализ полученных решений при граничных условиях третьего рода на внутренней поверхности сферы и граничных условиях третьего рода на внешней поверхности сферы показал, что они переходят в решения при граничных условиях первого рода на внутренней поверхности и граничных условиях третьего рода на внешней поверхности при условии гп —> г0Ла 0 ^ что является дополнительным доказательством правильности обоих решений. Полученное в работе аналитическое решение квазистационарной задачи о намораживании (задачи Стефана) на внешней сферической поверхности устанавливает непосредственную зависимость между продолжительностью намораживания и толщиной намороженного слоя в зависимости от геометрических и режимных параметров. В существующих работах показано, что при прочих равных условиях толщина намороженного слоя больше на цилиндрической поверхности, чем на плоской.
Сравнительные расчеты, полученных в исследовании решений, показывают, что при прочих равных условиях толщина намороженного слоя на внешней сфери-
ческой поверхности больше, чем на плоской и цилиндрической поверхностях.
Предложенное аналитическое решение квазистационарной задачи о намораживании на внешней сферической поверхности с граничными условиями третьего рода на внутренней поверхности и граничными условиями третьего рода на внешней поверхности, более общее, чем существующее, которое необходимо использовать для решения задачи о намораживании льда на рабочих поверхностях аккумуляторов холода.
Список литературы
1. Семенов Е. В., Бабакин Б. С., Воронин М. И., Ла-саро М. Моделирование процесса генерации льда и инея на поверхности воздухоохладителя // Вестник Международной академии холода. 2009. № 4. С. 5-9.
2. Белозеров Г. А., Бабакин Б. С., Макаров Б. А. Математическое моделирование продолжительности процесса замораживания и плавления эвтектического раствора в аккумуляторах холода // Известия Калининградского государственного технического университета. 2011. № 23. С. 141-147.
3. Семенов Е. В., Бабакин Б. С., Воронин М. И., Мо-рехон Л. Моделирование процесса инееобразования на поверхности воздушного прибора охлаждения // Вестник Международной академии холода. 2010. № 1. С. 36-39.
4. Семенов Е. В., Бабакин Б. С., Воронин М. И., Выгодин В. А., Бабакин С. Б. Моделирование процесса пневмокриоэлектросепарации сырья биологического происхождения // Вестник Международной академии холода. 2013. № 2. С. 62-66.
References
1. Semenov E. V. Babakin B. S., Voronin M. I., Lasaro M. Vestnik Mezdunarodnoj akademii holoda. 2009. No. 4. pp. 5-9.
2. Belozerov G. A. Babakin B. S., Makarov B. A. Izvestija Kaliningradskogo gosudarstvennogo tehnicheskogo universiteta. 2011. No. 23. pp. 141-147.
3. Semenov E. V., Babakin B. S., Voronin M. I., Morehon L. Vestnik Mezhdunarodnoj akademii holoda. 2010. No. 1. pp. 36-39.
4. Semenov E. V, Babakin B. S., Voronin M. I., Vygodin V A., Babakin S. B. Vestnik Mezhdunarodnoj akademii holoda. 2013. No. 2. pp. 62-66.
