УДК 621.56
Математическое моделирование процесса охлаждения хладоносителя системой замороженных шаров
Д-р техн. наук Е. В. СЕМЕНОВ [email protected]
Московский Государственный Университет Технологий и управления им. К. Г.Разумовского (Первый Казачий Университет) 109004, г. Москва, ул. Земляной вал, 73 Д-р техн. наук Б. С. БАБАКИН Московский Государственный Машиностроительный Университет (МАМИ) канд. техн. наук М. И. ВОРОНИН Московский Государственный Университет Пищевых Производств Канд. техн. наук А. Г. БЕЛОЗЁРОВ, канд. техн. наук С. Б. БАБАКИН Всероссийский научно-исследовательский институт холодильной промышленности (ВНИХИ)
Проблемы интенсификации производственных процессов и рационального использования сырья связаны с холодильной обработкой продуктов с целью сохранения их высоких биохимических показателей и снижения микробиологической порчи. Научное обоснование выбора способа и технических средств проведения процесса холодильной обработки определяет конечное товарное качество продуктов и эффективность работы предприятия в целом. Процесс охлаждения жидких сред относится к наименее изученной стадии обработки продукта холодом, а используемые для этой цели технические средства часто не отвечают современным требованием. В настоящее время перспективным направлением совершенствования техники и технологии охлаждения жидких сред является использование для охлаждения последних емкостей различной геометрии заполненных хладоносителем. Полученные результаты позволят оптимизировать процесс охлаждения жидких сред, снизить затраты на энергию и оборудование. Проблема охлаждения жидкостной среды путем интродуцирова-ния в данную среду коллектива замороженных шаров формулируется в рамках краевой задачи для уравнения теплопроводности. Решение данной задачи, полученное по методу осреднения, в области реальных значений параметров процесса по объекту исследования используется в качестве основы для численного эксперимента по моделированию протекания процесса охлаждения водной среды.
Ключевые слова: жидкостная система, охлаждение, теплопередача, шар, краевая задача, эффективность процесса.
Информация о статье
Поступила в редакцию 13.05.2016, принята к печати 24.10.2016 doi: 10.21047/1606-4313-2016-15-4-74-79 Ссылка для цитирования
Семенов Е. В., Бабакин Б. С., Воронин М. И., Белозёров А. Г., Бабакин С. Б. Математическое моделирование процесса охлаждения хладоносителя системой замороженных шаров // Вестник Международной академии холода. 2016. № 4. С. 74-79.
Mathematical modeling of thermostating liquid cooling process by the system of frozen ballons
D. Sc. E. V. SEMENOV
Moscow State University of Technologies and Management Named after K. G. Razumovskiy, Russia, Moscow
D. Sc. B. S. BABAKIN Moscow State University of Mechanical Engineering, Russia, Moscow Ph. D. M.I. VORONIN Moscow National University of Food Production, Russia, Moscow Ph. D. A.G. BELOZEROV, Ph.D S. B. BABAKIN All-Russian Research Institute of the refrigerating industry, Russia, Moscow
Problems of intensification of production processes and the efficient use of raw materials related to the cold treatment products in order to preserve their high biochemical parameters and reducing microbiological spoilage. The scientific rationale for the choice of the method and technical means of the cold treatment process determines the final quality of commodity products and the efficiency of the enterprise as a whole. The process of cooling liquid media refers to the least explored stage of cold treatment of the product, and used for this purpose technical devices often do not meet modern requirements. Currently, promising direction of improving technology and cooling technology is the use of liquid media for
cooling the latter containers filled with different geometry coolant. The results will allow to optimize the process of cooling liquid environments, reduce energy and equipment costs. The problem of cooling the liquid medium by introducing into the medium band of frozen balls is formulated within the framework of a boundary value problem for the heat equation. The solution to this problem, obtained by the averaging method in the field of real values on the object of study ofprocess parameters used as the basis for the numerical modeling of the process medium flow cooling the aqueous medium. Keywords: liquid, cooling, heat transfer, ball, boundary value problem, the process efficiency.
Одной из важных проблем интенсификации производственных процессов и рационального использования сырья, в частности, в молочной промышленности является охлаждение молока и молочных продуктов, что способствует сохранению его биологических свойств, предотвращает размножение микрофлоры в продукте. От того, насколько обоснован способ и технические средства проведения процесса охлаждения зависит товарное качество продукции и эффективность работы предприятия в целом. Вместе с тем до сих пор процесс охлаждения сырого молока относится к наименее изученной стадии обработки продукта холодом, а используемые для этой цели технические средства не в полной мере отвечают современным требованиям.
В настоящее время к одному из перспективных направлений совершенствования техники и технологии переработки молочной продукции относится ее охлаждение хладоносителем с интродуцированной в него системой замороженных до низких температур, наполненных жидкостью (хладоносителем) шаров. Данный технологический способ, в принципе, используют в смежных отраслях промышленности, в частности, в металлургии [1].
Результаты проведенных исследований выявили преимущества способа охлаждения водной среды с помощью замороженных шаров перед другими: отмечена высокая интенсивность процесса, снижение энергозатрат. Применительно к технологии молочного производства использование предложенного способа позволит в целом повысить количественные и качественные показатели эффективности процесса обработки продукта. В то же
Рис. 1. Схема расположения шаров в водяной рубашке
время имеется острая потребность в дальнейших углубленных исследованиях по этой важной научно-технической проблеме.
Постановка задачи
Пусть в некотором рабочем объеме имеют жидкостную систему, включающую жидкую и твердую фазы: жидкостную среду и коллектив замороженных шаров с хладоносителем. Предполагая, что группа шаров статистически равномерно распределена в жидкостной системе, в качестве модели процесса теплопередачи от жидкостной среды к шарам выбирают изолированный шар радиусом R. К данному шару из прилегающей к нему области R < г < l подается теплота из жидкостной среды (рис. 1, 2), где г — радиальная координата, 2l — среднее расстояние между центрами шаров (рис. 2).
При этом, если V и V — суммарный объем шаров и жидкости в рабочем объеме, соответственно, а с — объемная концентрация шаров в жидкостной системе, то, согласно принятой геометрической модели жидкостной системы «жидкостная среда + система шаров», половина расстояния l между двумя соседними шарами составит (рис. 3)
I = Я\[с 0<с<1. (1)
Полагают, что подвод теплоты к шару радиусом R реализуется симметричным образом, и, кроме того, для выбранной расчетной модели, с удовлетворительной точностью выполнены все, обычно принимаемые при анализе процесса теплопереноса в сплошном теле, допуще-
Рис. 2. Схема к расчету процесса размораживания шара
Замороженный шар
э т
д2Т 2 Э7\
Эт
дг
дг
а = U су,
э Г(/,т) Эг
= 0 (0 < т < да),
С этой целью предварительно вводят замену переменной
Т = и + Тк, (8)
в результате чего граничные условия (5), (6) становятся однородными по переменной и — приведенной температуре
и (Я, т) = 0, (9)
Рис. 3. Схема к расчету процесса теплопередачи в жидкостной системе
ния по физико-механическим свойствам изучаемого объекта [2-9].
Решение задачи
В качестве потока q теплоты в направлении оси г в принятой одномерной модели теплопереноса в сферической системе отсчета принимают
ди(1,т) дг
= 0, 0 < т < да.
(10)
Начальное условие принимает форму
и (г, 0) = ДТ ^ < г < I), (11)
где ДТ = Т - Т > 0 — перепад температур по объему жидкостной среды в начальный период времени.
При этом уравнение (3) по переменной и вследствие (8) сохраняет свой вид
Эи ,Э2и 2 ди, от дг г дг
(12)
(2)
В целях упрощения решения задачи, начальное условие (11) заменяют на усредненное по прилегающей к шару области г О (Я, I) значение
где Т—температура; X — коэффициент теплопроводности.
Тогда в качестве исходного соотношения, описывающего кинетику распределения температуры в жидкостной среде, выбирают, как обычно, отнесенное к сферическим координатам, с началом координат в центре шара, уравнение нестационарной теплопроводности [7] (рис. 1).
1 i
-ju(r,0)dr = AT .
I R n
(13)
В свою очередь, усредняя левую часть уравнения (12) по интервалу R < г < I, имеют
(3)
1-RidT
(14)
где Т — температура; т — время; г — радиальная координата; а — коэффициент температуропроводности,
(4)
здесь с — удельная теплоемкость; у — плотность шара.
За начальное условие, выбираемое как равномерное по объему жидкостной среды (воды), принимают
Т (r, 0) = Тс = const (R < r < l). (5)
Граничное условие на поверхности шара
Т (R, т) = Тк (0 < т < да), (6)
где Тк — криоскопическая температура воды (температура плавления льда).
Граничное условие симметричности поля температуры между соседними шарами [7]
где ф(т) — функция, подлежащая определению.
В результате уравнение (12) заменяют приближенным
д2и 2 дих „ „
(ТТ + -Т") = ф(Т). (15)
дг г дг
Так как правая часть уравнения (15) зависит лишь от т, то, интегрируя его дважды по г, находят
u (r, т) = D + D2 / r + г2ф (т) / 6,
(16)
(7)
где Б1, Б2 — произвольные величины, зависящие от т.
Согласуя (16) с граничными условиями (9), (10), получают
где I связана с R зависимостью (1).
Хотя решение задачи (3)-(7) известно, оно малопригодно для практических расчетов, так как имеет вид ряда. Поэтому, в некоторых случаях, применяют приближенный метод решения данной задачи, аналогичный методу осреднения [8].
и поэтому вместо (16) имеют
u (r, т) = у (r) ф (т),
где \|/(г) = —
6 R
- +—+—. Зг 6
(17)
(18)
Подставляя выражение (17) в уравнение (14), приходят к дифференциальному уравнению относительно искомой функции ф
1 \ ¿ср
1-Я{ йX
или в— = шр, йх
где 5 =
(19)
(20)
ах.
м(г,х) = £)3^(г)ехр(—).
а
а Л
чХг)=г-!-<0,
согласно (8) и (25) производная (22) по переменной г положительна
эт
дг
д и А Г г3-/3
ах.
дг ЪВ (В < 0, г < I),
•ехр(—)>0 В
Поток тепла от жидкостной среды к поверхности шара согласно (2), (26) вычисляют по зависимости
.. .Эи(Д,х) Ы ТВ?-1Ъ .ах. 2
) ъв я2 еМ~в)' Вт/м2. (27)
Количество тепла Q, подводимого к единице поверхности шара за время т расчитывают, интегрируя (27) по данному временному интервалу
^ Г /-\1 Ц13-Яъ)„ ,ах
[1-ехр(—)], Дж / м2. (28)
о
Интегрируя (19), получают
ф (т) = Б3ехр (ат / В), (21)
где Б3 — произвольная постоянная.
Поэтому, в соответствии с выражениями (17), (21) имеют
Следовательно, количество тепла, передаваемого шару от жидкостной среды, за счет теплопроводности в течение периода времени 1 составляет
Я* = бу*
(29)
(22)
где * — площадь поверхности шара, буд вычисляют согласно (28).
И поэтому, вследствие (28), (29) имеют
С целью определить значение Б , подставляют выражение (22) в (13)
/ л,
или Б3В = АТ, (23)
где В определяют по уравнению (20).
В результате чего, вследствие (22), (23) окончательно имеют
4тЛ(/3 -Я3)АТ
[1-ехрф]. 3 а В
(30)
Если Т0 и Тк — соответственно, температура замороженного шара и криоскопическая температура, то количество бш теплоты в шаре составляет
бш=т К Т - То) + бфД Дж
или
(24)
бш = 4 / 3*Я3у [сш (Тк - То) + бфудЪ Дж.
(31)
Как видно, что найденное по (24) выражение приведенной температуры и (г, т) имеет простую и удобную для количественного и качественного анализа мультипликативную форму. А именно, приведенная температура пропорциональна перепаду температуры АТ, и, поскольку, согласно расчету, В < 0, температура жидкостной среды убывает по экспоненте с течением времени т, и, тем быстрее, чем больше коэффициент коэффициент температуропроводности а.
Помимо этого, учитывая, что
где т1 — масса хладагента в шаре; у — плотность хладагента; сш — удельная теплоемкость; бфуд — удельная теплота плавления льда.
Пусть п — удельная теплота шара, расходуемая на охлаждение жидкостной среды такая, что
П = бш / б* (0 < П< 1),
(32)
(25)
где б* рассчитывается по уравнению (30), бш — по (31).
Поэтому, если необходимо рассчитать период времени охлаждения жидкостной среды, когда шар, отдает П величины своего холода, то тогда разрешают относительно т уравнение (32) или
бш - б* П =
(33)
(26)
то отсюда следует заключение, что температура жидкостной среды возрастает в интервале Я<г<1 и, следовательно, убывает по направлению к шару.
Отмеченные особенности процесса теплопередачи находятся в согласии с физической стороной описываемого явления.
Численный эксперимент
Рассматривается пример расчета процесса охлаждения жидкостной среды за счет размораживания коллектива шаров.
Пусть исходная жидкостная среда характеризуется параметрами: начальная температура Т = 290 К; криоскопическая температура Т = 273 К; плотность у = 1000 кг/м3; коэффициент теплопроводности = 0,56 Вт/(мК); удельная теплоемкость с = 4,2105 Дж/(кг-К) [10].
15
т
10
5
0
.2
Рис. 4. Зависимости периода обработки (р от коэффициента удельной теплоты шара (И) для различных значений радиуса шаров (Я) и коэффициента объемной концентрации шаров (с) в жидкостной системе: с = 0,01: 1 — Я = 1 см; 2 — Я = 2 см; с = 0,1: 3 — Я = 1 см, 4 — Я = 2 см
В качестве исходных геометрических и физико-механических параметров процесса для жидкостного хладагента (льда) в шаре выбирали (по порядку величин): радиус шара Я = 0,05 м; плотность льда у = 900 кг/м3; температура замороженного шара Т0 = 261 К; удельная теплоемкости сл = 2100 Дж/(кг-К); удельная теплота плавления Qфул = 3,35-105 Дж/кг.
На базе (33) для принятых значений параметров процесса и радиуса шаров Я = 1 см, Я = 2 см, а также значений коэффициента объемной концентрации шаров в жидкостной системе с = 0,01; с = 0,1, в виде графиков кривых рис. 4 показаны зависимости периода т обработки от коэффициента п — удельной теплоты шара, расходуемой на охлаждение воды.
Из вида графиков, показанных на рис. 4 выявляется согласие полученных результатов количественного моделирования исследуемого процесса с физическим смыслом задачи. А именно, возрастание периода обработки жидкостной среды когда, с одной стороны, увеличивается коэффициент п, т. е. отдача тепла шарами растет, а также, с другой стороны, когда увеличивается размер шара вследствие уменьшения его удельной поверхности.
Кроме того, естественно, период обработки жидкостной среды сокращается при увеличении объемной концентрации шаров (кривые 1,3 и 2,4 рис. 4).
Выводы
Таким образом:
— для обоснования базирующейся на явлении теплопередачи от жидкостной среды к коллективу замороженных шаров инновационной технологии охлаждения жидкостной среды предлагается аналитический аппарат по прогнозированию процесса теплопереноса в шаре;
— с целью количественного моделирования проблемы переноса теплоты в шаре, как сплошном теле заданным размером и физико-механическими параметрами процесса, используется закон сохранения энергии в фор-
ме уравнения теплопроводности, решение которого отыскивается приближенным способом;
— исходя из осесимметричной модели процесса переноса теплоты в прилегающей к шару жидкостной среде, на базе найденного решения получены удобные для инженерных расчетов и прогнозирования протекания процесса аналитические зависимости для определения распределения температуры по радиусу и по времени;
— на основе найденных, учитывающих все основные параметры задачи, явных зависимостей по распределению температуры в жидкостной среде в области реальных значений параметров процесса, проведено численное моделирование процесса охлаждения данной среды;
— полученные расчетные данные по моделированию процесса охлаждения воды адекватны физическому смыслу исследуемой проблемы.
Литература
1. Виноградов В. В., Шакуров А. Г., Тяжельникова И. Л., Виноградова Е. П., Есенбеков В. С. Математическое моделирование охлаждения шлакового расплава системой металлических шаров // Журнал технической физики. 2015. Том 8. Вып. 12. С. 21-25.
2. Семенов Е. В., Бабакин Б. С., Воронин М. И., ЛасароМ. Моделирование процесса генерации льда и инея на поверхности воздухоохладителя // Вестник Международной академии холода. 2009. № 4. С. 5-9.
3. Семенов Е. В., Бабакин Б. С., Воронин М. И., Морехон Л. Моделирование процесса инееобразования на поверхности прибора охлаждения // Вестник Международной академии холода. 2010. № 1. С. 36-39.
4. Белозеров Г. А., Бабакин Б. С., Макаров Б. А. Математическое моделирование продолжительности процесса замораживания и плавления эвтектического раствора в аккумуляторах холода // Известия Калининградского государственного технического университета. 2011. № 23. С. 141147.
5. Бредихин А. С., Семенов Е. В. и др. К вопросу о моделировании процесса размораживания мяса в вакууме. Ч. I // Хранение и переработка сельхозсырья. 2011. № 4. С. 25-29.
6. Бредихин А. С., Семенов Е. В. и др. К вопросу о моделировании процесса размораживания мяса в вакууме. Ч. II // Хранение и переработка сельхозсырья. 2011. № 5. С. 14-16.
7. Лыков А. В. Теория теплопроводности. — М.: Высшая школа, 1967. 600 с.
8. Битюков В. К., Тихомиров С. Г., Саввин С. С., Арапов Д. В. Математическая модель охлаждения оборотной воды в градирне с механической тягой. // Вестник Воронежского государственного университета инженерных технологий. 2014. № 1. С. 51-55. D0I:10.20914 / 2310-1202-2014-1-51-55
9. Лыков А. В. Тепломассобмен. — М.: Энергия, 1978. 480 с. 10. Напалков Г. Н. Тепло-массоперенос в условиях инееобразования. — М.: Машиностроение, 1983. 190 с.
References
1. Vinogradov V. V., Shakurov A. G., Tyazhel'nikova I. L., Vino-gradova E. P., Esenbekov V. S. Mathematical modeling of cooling of slag fusion with system of metal spheres. Zhurnal tekh-nicheskoi fiziki. 2015. Vol. 8. No. 12. p. 21-25. (in Russian)
0 0.05 0.1 0.15 0
h
физика
79
2. Semenov E. V., Babakin B. S., Voronin M. I., Lasaro M. Modeling of process of generation of ice and hoarfrost on an air cooler surface. VestnikMezhdunarodnoi akademii kholoda. 2009. No 4. p. 5-9. (in Russian)
3. Semenov E. V., Babakin B. S., Voronin M. I., Morekhon L. Modeling of process of an ineeobrazovaniye on the surface of the device of cooling. Vestnik Mezhdunarodnoi akademii kholoda. 2010. No 1. p. 36-39. (in Russian)
4. Belozerov G. A., Babakin B. S., Makarov B. A. Mathematical modeling of duration of process of freezing and melting of the eutectic solution in cold accumulators. Izvestiya Kaliningrad-skogo gosudarstvennogo tekhnicheskogo universiteta. 2011. No 23. p. 141-147. (in Russian)
5. Bredikhin A. S., Semenov E. V. i dr. To a question of modeling of process of defrosting of meat in vacuum. P. I. Khranenie i pererabotka sel'khozsyr'ya. 2011. No 4. p. 25-29. (in Russian)
6. Bredikhin A. S., Semenov E. V. i dr. To a question of modeling of process of defrosting of meat in vacuum. P. II. Khranenie i pererabotka sel'khozsyr'ya. 2011. No 5. p. 14-16. (in Russian)
7. Lykov A. V. Theory of heat conductivity. Moscow, Vysshaya shkola, 1967. 600 p. (in Russian)
8. Bitiukov V. K., Tikhomirov S. G., Savvin S. S., Arapov D. V. The mathematical model of cooling recycled water in a cooling tower with mechanical traction. Proceedings oof the Voronezh State University of Engineering Technologies. 2014; (1):51-55. (in Russian)
9. Lykov A. V. Heatmass exchange. Moscow. Energiya, 1978. 480 p. (in Russian)
10. Napalkov G. N. Warm mass transfer in the conditions of an ineeobrazovaniye. Moscow. Mashinostroenie, 1983. 190 p. (in Russian)
Поздравляем с юбилеем Георгия Автономовича Белозерова
29 октября 1946 г. исполнилось 70 лет Георгию Автономовичу Белозерову, доктору технических наук, главному научному сотруднику Всероссийского научно-исследовательского института холодильной промышленности (ВНИХИ), вице-президенту и действительному члену Международной академии холода.
После окончания Московского института химического машиностроения Белозеров Г. А. трудится на ряде крупных предприятий страны, таких как НПО «Криогенмаш», ВНИИторгмаш, НПО «Агро-холодпром». С 1988 г. по 2006 г. Георгий Автономович работает заместителем директора по научной работе ВНИИ холодильной промышленности (ВНИХИ), а затем занимает должность директора ВНИХИ, которым он успешно руководил до июня 2016 г.
Научные труды Белозерова Г. А. широко известны в России и за ее пределами, среди них девять монографий, свыше 50 авторских свидетельств, сотни публикаций по ресурсо- и энергосберегающим технологиям холодильной обработки и хранения пищевых продуктов, по созданию и внедрению в индустрию холода холодильно-технологического оборудования, отвечающего современному научно-техническому уровню.
Белозеров Г. А. — член редколлегий ведущих научно-технических журналов: «Холодильная техника», «Мир мороженого и быстрозамороженных продуктов», «Вестник Международной академии холода», «Все о мясе», «Холодильный бизнес». Является вице-президентом Союза мороженщиков России, членом подкомитета Торгово-промышленной палаты, членом комиссии Д2 Международного института холода (Париж), имеет государственные награды, медали ВДНХ, почетные грамоты Ми-нобрнауки Российской Федерации, Россельхозакадемии, правительства г. Москвы.
Коллеги по холодильному сообществу знают Георгия Автономовича как душевного, доброжелательного человека, обладающего выдающимся организаторским талантом, неиссякаемой энергией, творческой работоспособностью.
Президиум Международной академии холода, редколлегия журнала «Вестник МАХ» сердечно поздравляют Вас, Георгий Автономович, с юбилеем!
Искренне желаем Вам крепкого здоровья, энергии, творческого задора, профессиональных достижений, благополучия в семье!