Математическое моделирование продолжительности процесса замораживания и плавления эвтектического раствора в аккумуляторах холода

Размещено на http://www.allbest.ru

Размещено на http://www.allbest.ru

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

За рубежом данный вид охлаждения широко используется, в частности, фирма «Carrier Transicold» (Франция) предлагает серию установок с машинно-аккумуляционным охлаждением “Vatna” для изотермических автофургонов объемом от 4 до 23 м3 с температурой замерзания эвтектического раствора минус 32°С. Установки предназначены для перевозки мороженого и быстрозамороженных продуктов и рассчитаны на 14 ч непрерывной работы при температуре окружающего воздуха 30°С.

Мировой лидер в производстве эвтектических плит - фирма «FIT s.p.a.» из Италии предлагает на российском рынке серии эвтектических плит моделей EBS и EFR.

Ниже рассматривается математическая модель продолжительности процесса замораживания (плавления) аккумуляционных охлаждающих приборов (см. рис. 1, 2).

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

Рис. 1. Температурное поле в плоской неограниченной пластине при ее охлаждении

Продолжительность процесса кристаллизации эвтектического раствора неограниченной пластины можно представить как сумму трёх составляющих: продолжительности процесса охлаждения жидкой фазы эвтектического раствора ф1 от начальной температуры t0 до начала процесса кристаллизации; продолжительности процесса кристаллизации ф2 и продолжительности процесса переохлаждения твёрдой фазы ф3 до некоторой конечной температуры tk в центре неограниченной пластины.

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

,. (1)

При граничных начальных условиях:

;

;

, (2)

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

Решение уравнения (1) при граничных условиях (2) имеет следующий вид:

, (3)

где ; мn - тень корня характеристического уравнения, ; - число Фурье; - критерий Bio для жидкой фазы эвтектического раствора.

Для больших чисел Фурье можно ограничиться первым членом ряда (3), и тогда приближённо продолжительность процесса охлаждения жидкой фазы можно будет определить по следующей зависимости:

,

где:

.

Точное значение ф1 можно определить путём численного решения уравнения (3) относительно числа Фурье Foж при:

и .

Наибольшую сложность представляет определение продолжительности процесса ф2 кристаллизации эвтектического раствора.

Для приближённого определения продолжительности процесса кристаллизации эвтектического раствора воспользуемся вторым методом Л.С. Лейбензона.

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

,

где о - координата границы раздела фаз.

Рис. 2. Температурное поле в плоской неограниченной пластине при замораживании эвтектического раствора

Текущее значение энтальпии твёрдой фазы будем определять по следующей зависимости:

, (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.

Таким образом, продолжительность всего процесса кристаллизации эвтектического раствора будет равна:

или в безразмерной форме:

.

Аналогичным образом может быть рассмотрена обратная задача - плавление эвтектического раствора неограниченной пластины.

Размещено на Allbest.ru