Внастоящее время перспективным способом безмашинного охлаждения является применение аккумуляторов холода в авторефрижераторах при внутригородских перевозках
| Вид материала | Документы |
- Система охлаждения двигателя автомобиля, 231.89kb.
- Заочный Государственный Университет Внастоящее время все большее применение в разработке, 64.47kb.
- Математическое моделирование процесса каталитического получения четырехвалентного урана, 28.23kb.
- Путь повышения эффективности лекций, 69.69kb.
- Доклад внастоящее время индустрия туризма является одной из наиболее динамично развивающихся, 88.22kb.
- Шкаф автоматического управления системой охлаждения кочегаров И. А, ведущий конструктор, 52.15kb.
- Правила составления актов при перевозках грузов железнодорожным транспортом общие положения, 308.75kb.
- Внастоящее время школа является важнейшим фактором ускорения социально-экономического, 1473.01kb.
- Т н. Прохоренко Н. Н. Обратные термодинамические циклы и холодильные установки. Термодинамика, 190.29kb.
- Внастоящее время в промышленности, в сельском хозяйстве, в армии и на флоте находят, 318.64kb.
УДК 664.7
Математическое моделирование продолжительности
процесса замораживания и плавления эвтектического
раствора в аккумуляторах холода
Г.А. Белозеров*, Б.С. Бабакин**, Б.А. Макаров***
*Государственное научное учреждение «Всероссийский научно-исследовательский холодильный институт», г. Москва, ул. Костякова, 12
**Московский государственный университет прикладной биотехнологии, г. Москва, ул.Талалихина, 33, e-mail: babakin@bk.ru
***Открытое акционерное общество Центральный научно-исследовательский институт «Курс», г. Москва, ул. Авиамоторная, 10
Длительность процесса кристаллизации эвтектического раствора неограниченной пластины могла бы быть представлена как сумма трех компонентов: длительность процесса охлаждения эвтектического раствора жидкой фазы 1 первоначальной температуры t0 до начала кристаллизации; длительность процесса кристаллизации 2 и длительность процесса предварительного охлаждения твердой фазы 3 до некоторой конечной температуры tк в центре пластины.
математическое моделирование, замораживание, плавление, эвтектический раствор, аккумулятор холода
В аккумуляторных системах охлаждения используют теплоту плавления эвтектических растворов. В настоящее время перспективным способом безмашинного охлаждения является применение аккумуляторов холода в авторефрижераторах при внутригородских перевозках. Применение аккумуляционного охлаждения в авторефрижераторах создает экологическую чистоту окружающей среды.
За рубежом данный вид охлаждения широко используется, в частности, фирма «Carrier Transicold» (Франция) предлагает серию установок с машинно-аккумуляционным охлаждением “Vatna” для изотермических автофургонов объемом от 4 до 23 м3 с температурой замерзания эвтектического раствора минус 32°С. Установки предназначены для перевозки мороженого и быстрозамороженных продуктов и рассчитаны на 14 ч непрерывной работы при температуре окружающего воздуха 30°С.
Мировой лидер в производстве эвтектических плит – фирма «FIT s.p.a.» из Италии предлагает на российском рынке серии эвтектических плит моделей EBS и EFR.
Ниже рассматривается математическая модель продолжительности процесса замораживания (плавления) аккумуляционных охлаждающих приборов (см. рис. 1, 2).
Приближённое определение продолжительности процесса кристаллизации эвтектического раствора неограниченной пластины принято при граничных условиях третьего рода.
Рис.1. Температурное поле в плоской неограниченной пластине
при ее охлаждении
Fig.1. The temperature field in a plane infinite plate at cooling
Продолжительность процесса кристаллизации эвтектического раствора неограниченной пластины можно представить как сумму трёх составляющих: продолжительности процесса охлаждения жидкой фазы эвтектического раствора τ1 от начальной температуры t0 до начала процесса кристаллизации; продолжительности процесса кристаллизации τ2 и продолжительности процесса переохлаждения твёрдой фазы τ3 до некоторой конечной температуры tk в центре неограниченной пластины.
Если пренебречь конвективным теплообменом, то продолжительность процесса охлаждения жидкой фазы τ1 от начальной температуры t0 до температуры t3 кристаллизации эвтектического раствора на внешней поверхности может быть определена из решения одномерного нестационарного уравнения теплопроводности:
, 
. (1)При граничных начальных условиях
;
;
, (2)где аж, λж – соответственно коэффициенты температуропроводности и теплопроводности жидкой фазы; α – коэффициент теплоотдачи на внешних поверхностях неограниченной пластины; tо.с – температура окружающей среды, оС.
Решение уравнения (1) при граничных условиях (2) имеет следующий вид:
, (3)где
; μn – тень корня характеристического уравнения, 
; 
– число Фурье; 
– критерий Bio для жидкой фазы эвтектического раствора.Для больших чисел Фурье можно ограничиться первым членом ряда (3), и тогда приближённо продолжительность процесса охлаждения жидкой фазы можно будет определить по следующей зависимости:
,где
.Точное значение τ1 можно определить путём численного решения уравнения (3) относительно числа Фурье Foж при
и 
.Наибольшую сложность представляет определение продолжительности процесса τ2 кристаллизации эвтектического раствора.
Для приближённого определения продолжительности процесса кристаллизации эвтектического раствора воспользуемся вторым методом Л.С. Лейбензона.
Предположим, что в каждый момент времени распределение температур в твёрдой фазе определяется по линейному закону (квазистационарная модель)
,где ξ – координата границы раздела фаз.
Рис.2. Температурное поле в плоской неограниченной пластине
при замораживании эвтектического раствора
Fig.2. The temperature field in a plane infinite plate at freezing eutectic solution
Текущее значение энтальпии твёрдой фазы будем определять по следующей зависимости:
, (4)где wT3 – энтальпия твёрдой фазы эвтектического раствора при температуре кристаллизации; СPT – теплоёмкость твёрдой фазы.
Аналогичным образом определим текущее значение энтальпии жидкой фазы:
, (5)где wжз – энтальпия жидкой фазы эвтектического раствора при температуре кристаллизации; Срж – удельная теплоёмкость жидкой фазы.
Уравнение теплового баланса для каждого момента времени можно записать следующим образом:
. (6)В уравнении (6) первый член характеризует изменение во времени теплосодержания жидкой фазы; второй – изменение во времени теплосодержания твёрдой фазы; третий – изменение во времени теплового потока на внешних поверхностях пластины.
Для простоты в этом уравнении начальный момент времени принят равным нулю.
Определяем текущее значение первого и второго интегралов в уравнении (6):
. (7)Величину второго слагаемого приближённо определим следующим образом:
,где Qж – общее количество теплоты, которое необходимо отвести от жидкой фазы, чтобы охладить последнюю до температуры кристаллизации от состояния, в котором находилась жидкая фаза в момент начала кристаллизации:
,где Tж – пространственное распределение температуры в момент времени τ1 начала процесса кристаллизации, которое можно определить приближённо по формуле:
.Следовательно:
.Таким образом, первое слагаемое в уравнении (6) приближённо определяется следующим выражением:
. (8)Второе слагаемое в уравнении (6) определяется при линейном законе изменения распределения температуры по толщине твёрдой фазы:
После подстановки выражений (8) и (9) для первого и второго слагаемых в уравнение (10) и произведя дифференцирование, получим следующее уравнение, определяющее скорость перемещения границы раздела фаз:
(10)В уравнении (10) Н – скрытая теплота фазового перехода (кристаллизации), определяемая следующим образом:
.Интегрируя уравнение (10), получим формулу для определения продолжительности процесса кристаллизации эвтектического раствора:
или в безразмерной форме
. (11)В формуле (11) первое слагаемое определяет продолжительность процесса, необходимого для собственно кристаллизации и охлаждения жидкой фазы до температуры кристаллизации; второе – определяет продолжительность процесса, необходимого для переохлаждения твёрдой фазы.
В конце продолжительности процесса кристаллизации, согласно принятому допущению, распределение температуры по толщине твёрдой фазы будет определяться линейной зависимостью:
.Для определения продолжительности процесса переохлаждения твёрдой фазы τ3 необходимо решить уравнение теплопроводности (5) при следующих граничных и начальных условиях:
,
, (12)
.При этом в уравнении теплопроводности (1) коэффициент для жидкой фазы aж необходимо заменить на коэффициент температуропроводности твёрдой фазы эвтектического раствора аТ.
Это решение имеет следующий вид:
. (13)Коэффициенты Dn разложения начального распределения температуры определяются по формуле
Для больших значений числа Фурье продолжительность процесса переохлаждения твёрдой фазы τ3 можно определить, учитывая только первый член ряда Фурье (13):
или 
.Точное значение продолжительности процесса переохлаждения твёрдой фазы τ3 можно определить путём численного решения уравнения (13) относительно числа Фурье F0T3 при t(0,τ) = t1.
Таким образом, продолжительность всего процесса кристаллизации эвтектического раствора будет равна:
или в безразмерной форме
.Аналогичным образом может быть рассмотрена обратная задача - плавление эвтектического раствора неограниченной пластины.
Approximate computation of the time period of the eutectic solution freezing or melting in the cold accumulator
G.A. Belozerov, B.S. Babakin, B.A. Makarov
The duration of crystallization process of the eutectic solution of unrestricted plate could be presented as sum of three components: duration of the cooling process of the eutectic solution liquid phase τ1 from initial temperature t0 to the crystallization beginning; duration of the crystallization process τ2 and duration of the subcooling process of the solid phase τ3 to some final temperature tk at the center of the plate.
approximate computation, freezing, melting, eutectic solution, cold accumulator