Аналитическое решение задачи совместного тепломассопереноса при инфракрасном нагреве колбасных изделий
УДК 637.523.27:536.25
Аналитическое решение задачи совместного тепломассопереноса при инфракрасном нагреве
колбасных изделий
Д.т.н. Вороненко Б.А., д.т.н. Пеленко В.В., аспирант Стариков В.В.
Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий
Использование энергии инфракрасного излучения является одним из эффективных путей интенсификации тепловой обработки пищевых материалов, позволяющим значительно сократить длительность тепловой обработки и повысить качество готовых изделий. На основе предпосылок предложена математическая модель совместного тепломассопереноса тепловой обработки колбасного батона цилиндрической формы в коптильной установке с ИК-нагревом.
Ключевые слова: ИК-нагрев, тепломассоперенос, математическая модель, система дифференциальных уравнений, начальные и граничные условия.
Технологический процесс тепловой обработки колбасных изделий как сложных капиллярно-пористых коллоидных тел, состоит из стадий подсушки, обжарки и варки. Интересующая нас первая стадия (подсушка) заключается в прогреве, главным образом, поверхности продукта (Н.С. Михеевой показано, что углубление поверхности испарения имеет место с начала процесса сушки [1]) в среде с низкой относительной влажностью (до 10%). При подсушке с поверхности колбасного изделия удаляется влага смачивания (механически удерживающаяся вода), что способствует равномерной прокраске поверхности и диффузии в продукт коптильных веществ при последующей обжарке.
Как известно [2-4], использование энергии электромагнитного поля инфракрасного (ИК) диапазона является одним из эффективных путей интенсификации тепловой обработки пищевых материалов, позволяющим значительно сократить длительность тепловой обработки и повысить качество готовых изделий.
При составлении математической модели процесса термической обработки колбасных изделий в коптильных установках с ИК-нагревом конвективным теплообменом от окружающего воздуха пренебрегаем, причем, чем выше температура источника радиационного излучения, тем обоснованнее эта предпосылка [2, 8].
Лучистый поток, проникая в продукт, затухает по экспоненциальному закону [5].
При толщине продукта ℓ > 10 мм и высоком значении коэффициента поглощения лучистый поток быстро затухает по мере проникновения в продукт
[1], и можно считать, что вся энергия отдается поверхности, а в нагреве
внутренних слоев не участвует. В этом случае, особенно для периода
постоянной скорости сушки — стадии подсушки (от начала процесса до
времениτ 1 ), затухание лучистого потока, проникающего в продукт, может быть
описано параболическим законом [6]:
w(r )
=
w0
r2
−
λ2 пр
R2
−
λ2 пр
( )0 < λпр ≤ r ≤ R
(1)
Принимаем колбасный батон за тело основной геометрической формы —
бесконечный цилиндр, однородный и изотропный, теплофизические
характеристики которого в процессе копчения остаются неизменными.
При ИК-облучении нагрев происходит с поверхности, а внутренним слоям
колбасного батона теплота передается теплопроводностью и влагопровод-
ностью. На основе многих экспериментальных данных [3] не будем учитывать
процесс термодиффузии влаги, т.е. изменение влагосодержания за счет
термовлагопроводности, и учтем также сравнительно медленное изменение
влагосодержания со временем при ИК-нагреве [7].
В соответствии с указанными допущениями нагрев колбасного батона в
коптильной установке с ИК-нагревом может быть описан следующей краевой
задачей связанного тепло- и массопереноса: требуется решить систему
дифференциальных уравнений в частных производных в цилиндрических
координатах
∂t ∂τ
=
aq
∂ 2t ∂r 2
+
1
∂t
+
ερ
r ∂r cq
∂u + w(r) ; ∂τ cqγ 0
(2)
∂u ∂τ
=
∂ 2u ∂r 2
+
1 r
∂u ∂r
(3)
(0 < τ < τ1; 0 < r < R)
при следующих условиях:
t (r, 0) = t0 = const;
(4)
u(r,0) = u0 = const;
∂t (0,τ ) = ∂u (0,τ ) = 0;
∂r ∂r
t (0,τ ) < ∞;u (0,τ ) < ∞;
(5) (6) (7)
−λq
∂t
( R,τ
∂r
)
+
q
(τ
)
−
(1−
ε
)
qm
(τ
)
=
0;
(8)
amγ
0
∂u
( R,τ
∂r
)
+
qm
(τ
)
=
0
(9)
Здесь: (2) — уравнение теплопереноса с внутренним источником тепла,
определяемым зависимостью (1); (3) — уравнение влагопереноса; (4) и (5) —
начальные условия, задающие равномерное распределение температуры и
влагосодержания в теле в момент начала процесса нагрева; (6) и (7) — условия
симметрии и физической ограниченности температуры и влагосодержания;
(8) — граничное условие второго рода, задающее поток теплоты через
поверхность тела; (9) — граничное условие второго рода, описывающее
влагообмен тела с окружающей средой.
Аналитическое решение задачи (2) — (9) получено методом интегрального
преобразования Лапласа в следующем безразмерном виде:
T ( X , Fo) =
t(r,τ ) − t0 tмc − t0
= 2 1 +
2(1
Po −X
2 пр
)
13 768
−
X
2 пр
−
1 8
Fo
−
1 16
X4
+
1 32
X2
−
∑−
∞ n =1
2J0 (µnX )
µ
4 n
J
0
(
µ
n
)
e−
µ
2 n
F
o
+
ε KoLu 2(1 − Lu )
( )Pd2 J 0 Pd2 X − J1 Pd2
Pd2 Lu
J
0
J
1
Pd2
Lu
Pd2
Lu
X
e− Pd2Fo
−
( ) ∑−
Pd1 J0 2J1
Pd1 X Pd1
e− Pd1Fo
+
∞
1
n =1
1 −
µ
2 n
Pd1
−
(1 −
ε KoLu
Lu) 1 −
µ
2 n
Pd2
e−
µ
2 n
Fo
+
+
(1 −
ε KoLu
Lu) 1−
µ
2 n
Pd2
e−
µ
2 n
LuFo
J0(µn X )
J0 (µn )
для Lu=1 ( am = aq )
T (X , Fo) =
t(r,τ ) − t0 tмc − t0
= 2 1+
2(1
Po −X
2 пр
)
13 768
−
X
2 пр
−
1 8
Fo
−
1 16
X4
+1 32
X2
−
∑ ( )−
∞ n =1
2J0 µn X
µ
4 n
J
0
(
µ
n
)
e−
µ
2 n
Fo
−
( )Pd1 J 0 Pd1 X ( ) ∑J1 Pd1
( )e − Pd1Fo
+
∞ n =1
J0 µnX
1−
µ
2 n
Pd1
J0 (µn )
e
−
µ
2 n
Fo
(10) (11)
θ
(X
, Fo)
=
2
1 −
1 2
Pd2 Lu
J
0
J
1
∑P d 2
X
( )L u
Pd2
( )L u
e − Pd2Fo
∞
+
n =1
1 −
J0 µnX
µ
2 n
L
u
Pd2
J0
e
−
µ
2 n
F
o
;
µn
(12)
µn — последовательные положительные корни уравнения
J1(µ) = 0 .
(13)
Анализ полученных решений (1) и (12) показывает, что благодаря
быстрому увеличению абсолютной величины последовательного ряда значений
характеристических корней µn (уравнение (13)),
µ1 < µ2 < K < µn < K ,
а, следовательно, быстрого уменьшения экспоненциальных сомножителей
e ,−
µ
2 n
Fo
бесконечные
суммы,
входящие
в
решения,
сходятся
достаточно
быстро.
Поэтому, начиная с определенного значения числа Фурье, с заранее заданной
степенью точности из всего ряда разложения можно удовлетвориться одним-
двумя первыми членами. Такое упрощение имеет большое практическое
значение из-за существенного сокращения объема расчетной работы и
возможности представления общих решений в удобной для практического
применения форме.
Так, формулой удобной для инженерных расчетов полей потенциала
влагопереноса, будет следующее выражение, полученное из (12) при n = 1:
u(r,τ u0
)
=1−
2
Pd2 Lu
J
0
J
0
Pd2
X
Lu
Pd2
Lu
e− Pd2 Fo
+
1−
2J0 (µ1X )
µ12 Lu Pd2
J
0
(
µ1
)
e − µ12 Fo
.
(14)
Среднее значение безразмерной величины потенциала переноса влаги для
упрощенного выражения Ө(Х, Fo) при n = 1
∫ ( )θ ( F o ) =
u 0 − u (τ )
1
=2
X θ ( X , Fo)dx = 2
1 − e− Pd2Fo
u0 0
.
(15)
Ограничиваясь одним первым членом ряда, получаем приближенно
формулу для расчета времени, необходимого для достижения нагреваемым
телом определенного влагосодержания:
τ= θ
2k2 .
(16)
Среднее значение безразмерного потенциала температуры для
упрощенного выражения Т(Х, Fo) при n = 1
∫T ( F o ) =
1
2
XT (X , Fo)dx =
t(x,τ ) − t0
=
( )0 t мс − t0
Po 1 − X пр
2 + 2Pd1Fo
.
(17)
Для нахождения времени доведения продукта до необходимой средней
объемной температуры, запишем в следующем виде:
τ
=
2 k1
1
(t мс
−
t0
)
t
−
t0
−
w0R 2
λq
1 −
λ2 пр
R2
.
(18)
Выводы
Поставлена и аналитически решена задача совместного тепло- и массопереноса при инфракрасном нагреве колбасного батона цилиндрической формы.
Полученные решения позволяют рассчитать поля температуры и влагосодержания, усредненные значения соответствующих потенциалов переноса, темп нагрева, расход тепла в процессе теплообмена, а также получить формулы, удобные для инженерных расчетов.
Разработанные аналитические решения дают возможность определить время, необходимое для достижения нагреваемым продуктом определенной температуры и влагосодержания.
Обозначение
t(r, τ) — температура, К, °C; t0 − начальная температура; tмс − максимальная
температура
среды
в
камере;
∆t =tмс −t0;
T(X,
Fo)
=
(t(x,τ)
−t0) (tмс
−t0)
–
безразмерная
температура; u(r, τ) – влагосодержание, кг влаги/кг абс. сух. вещества; и0 −
начальное влагосодержание;
θ(X,
Fo)
=
(u0
−и(r,τ)) u0
–
безразмерное
влагосодержание;
r
–
текущая
координата, м; Х = r/R – безразмерная координата; R – радиус поперечного
сечения цилиндрического колбасного батона, м;
τ – время, с;τ1 − время окончания первого периода термической обработки
– подсушки; aq − коэффициент температуропроводности, м²/c; ε – коэффициент
фазового перехода (0 < ε < 1); ρ – удельная теплота фазового перехода, Дж/кг;
cq − удельная теплоемкость, Дж/(кг·K); am− коэффициент потенциало- (влаго-)
проводности, м²/c; λm − предельная глубина проникновения инфракрасного
излучения в материал, м; λq − коэффициент теплопроводности, Вт/(м·К);
q (τ ) = k1cqγ 0 R∆te −k1τ + Rk2 ργ 0u0e− k2τ – экспериментально найденная
плотность теплового потока на поверхности тела, Вт/м²; q m (τ ) = k 2γ 0u 0 R e − k2τ
– экспериментально установленная плотность потока массы вещества через
поверхность тела, кг/(м²·c); γ 0 – плотность абсолютно сухого вещества, кг/м³;
k 1 - коэффициент, характеризующий убывание температуры поверх-
ности тела по экспоненциальному закону, 1/с; k 2 - коэффициент сушки, 1/с; Fo = aqτ R2 − критерий гомохронности (число Фурье); Lu = am aq − критерий (число Лыкова) взаимосвязи массо- и теплопереноса; Ko = r∆u (cq∆t) − число
Коссовича; Pd1 = k1R2 aq , Pd1 = k2R2 aq - числа Предводителева (теплообменное и массообменное); J0(z) и J1(z) - функции Бесселя от вещественного аргумента z первого рода и нулевого порядка соответственно; w0 − мощность ИК-источника, Вт/м³.
Литература
1. Гинзбург А.С. Основы теории и техники сушки пищевых продуктов. – М.: Пищевая промышленность, 1973. – 528 с.
2. Гинзбург А.С. Инфракрасная техника в пищевой промышленности. – М.: Пищевая промышленность, 1966. – 376 с.
3. Рогов И.А., Некрутман С.В. СВЧ и ИК-нагрев пищевых продуктов. – М.: Пищевая промышленность, 1976. – 210 с.
4. Ильясов С.Г., Красников В.В. Физические основы инфракрасного излучения пищевых продуктов. – М.: Пищевая промышленность, 1978. – 359 с.
5. Лыков А.В. Тепло и массообмен в процессах сушки. – М.: Госэнергоиздат, 1956. – 464 с.
6. Белобородов В.В., Вороненко Б.А. Решение задачи нагрева тел в электромагнитном поле сверхвысоких частот //Журнал прикладной химии АН СССР, «Наука» - Ленинградское отделение, № 10, 1984. – С. 2276-2282.
7. Вороненко Б.А., Пеленко В.В., Стариков В.В. Постановка задачи тепломассопереноса процесса горячей сушки рыбы // Межвузовский сборник научных трудов «Ресурсосберегающие технологии и оборудование пищевой промышленности», СПб, СПбГУНиПТ, 2006. – С. 71- 75.
8. Стариков В.В., Вороненко Б.А. Применение ИК-нагрева при копчении // Межвузовский сборник научных трудов «Теория и практика разработки и эксплуатации пищевого оборудования», СПб, СПбГУНиПТ, 2007. – С. 12-17.
Analytical solution of a combined heat and mass transfer problem for infrared heating of sausage goods
B.A. Voronenko DSc, V.V.Pelenko DSc, V.V.Starikov graduate student
Saint-Petersburg State University of Refrigeration & Food Engineering
Use of infrared radiant energy is an effective way to enhance heat treatment of food materials, allowing both to curtail heat treatment period and to upgrade finished product quality. On the basis of prerequisites a mathematical model of combined heat and mass transfer is proposed for heat treatment of a long cylinder sausage loaf in an IR-heated smoker.
Key words: IR-heating, heat and mass transfer, mathematical model, differential equation set, initial and boundary conditions.
Аналитическое решение задачи совместного тепломассопереноса при инфракрасном нагреве
колбасных изделий
Д.т.н. Вороненко Б.А., д.т.н. Пеленко В.В., аспирант Стариков В.В.
Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий
Использование энергии инфракрасного излучения является одним из эффективных путей интенсификации тепловой обработки пищевых материалов, позволяющим значительно сократить длительность тепловой обработки и повысить качество готовых изделий. На основе предпосылок предложена математическая модель совместного тепломассопереноса тепловой обработки колбасного батона цилиндрической формы в коптильной установке с ИК-нагревом.
Ключевые слова: ИК-нагрев, тепломассоперенос, математическая модель, система дифференциальных уравнений, начальные и граничные условия.
Технологический процесс тепловой обработки колбасных изделий как сложных капиллярно-пористых коллоидных тел, состоит из стадий подсушки, обжарки и варки. Интересующая нас первая стадия (подсушка) заключается в прогреве, главным образом, поверхности продукта (Н.С. Михеевой показано, что углубление поверхности испарения имеет место с начала процесса сушки [1]) в среде с низкой относительной влажностью (до 10%). При подсушке с поверхности колбасного изделия удаляется влага смачивания (механически удерживающаяся вода), что способствует равномерной прокраске поверхности и диффузии в продукт коптильных веществ при последующей обжарке.
Как известно [2-4], использование энергии электромагнитного поля инфракрасного (ИК) диапазона является одним из эффективных путей интенсификации тепловой обработки пищевых материалов, позволяющим значительно сократить длительность тепловой обработки и повысить качество готовых изделий.
При составлении математической модели процесса термической обработки колбасных изделий в коптильных установках с ИК-нагревом конвективным теплообменом от окружающего воздуха пренебрегаем, причем, чем выше температура источника радиационного излучения, тем обоснованнее эта предпосылка [2, 8].
Лучистый поток, проникая в продукт, затухает по экспоненциальному закону [5].
При толщине продукта ℓ > 10 мм и высоком значении коэффициента поглощения лучистый поток быстро затухает по мере проникновения в продукт
[1], и можно считать, что вся энергия отдается поверхности, а в нагреве
внутренних слоев не участвует. В этом случае, особенно для периода
постоянной скорости сушки — стадии подсушки (от начала процесса до
времениτ 1 ), затухание лучистого потока, проникающего в продукт, может быть
описано параболическим законом [6]:
w(r )
=
w0
r2
−
λ2 пр
R2
−
λ2 пр
( )0 < λпр ≤ r ≤ R
(1)
Принимаем колбасный батон за тело основной геометрической формы —
бесконечный цилиндр, однородный и изотропный, теплофизические
характеристики которого в процессе копчения остаются неизменными.
При ИК-облучении нагрев происходит с поверхности, а внутренним слоям
колбасного батона теплота передается теплопроводностью и влагопровод-
ностью. На основе многих экспериментальных данных [3] не будем учитывать
процесс термодиффузии влаги, т.е. изменение влагосодержания за счет
термовлагопроводности, и учтем также сравнительно медленное изменение
влагосодержания со временем при ИК-нагреве [7].
В соответствии с указанными допущениями нагрев колбасного батона в
коптильной установке с ИК-нагревом может быть описан следующей краевой
задачей связанного тепло- и массопереноса: требуется решить систему
дифференциальных уравнений в частных производных в цилиндрических
координатах
∂t ∂τ
=
aq
∂ 2t ∂r 2
+
1
∂t
+
ερ
r ∂r cq
∂u + w(r) ; ∂τ cqγ 0
(2)
∂u ∂τ
=
∂ 2u ∂r 2
+
1 r
∂u ∂r
(3)
(0 < τ < τ1; 0 < r < R)
при следующих условиях:
t (r, 0) = t0 = const;
(4)
u(r,0) = u0 = const;
∂t (0,τ ) = ∂u (0,τ ) = 0;
∂r ∂r
t (0,τ ) < ∞;u (0,τ ) < ∞;
(5) (6) (7)
−λq
∂t
( R,τ
∂r
)
+
q
(τ
)
−
(1−
ε
)
qm
(τ
)
=
0;
(8)
amγ
0
∂u
( R,τ
∂r
)
+
qm
(τ
)
=
0
(9)
Здесь: (2) — уравнение теплопереноса с внутренним источником тепла,
определяемым зависимостью (1); (3) — уравнение влагопереноса; (4) и (5) —
начальные условия, задающие равномерное распределение температуры и
влагосодержания в теле в момент начала процесса нагрева; (6) и (7) — условия
симметрии и физической ограниченности температуры и влагосодержания;
(8) — граничное условие второго рода, задающее поток теплоты через
поверхность тела; (9) — граничное условие второго рода, описывающее
влагообмен тела с окружающей средой.
Аналитическое решение задачи (2) — (9) получено методом интегрального
преобразования Лапласа в следующем безразмерном виде:
T ( X , Fo) =
t(r,τ ) − t0 tмc − t0
= 2 1 +
2(1
Po −X
2 пр
)
13 768
−
X
2 пр
−
1 8
Fo
−
1 16
X4
+
1 32
X2
−
∑−
∞ n =1
2J0 (µnX )
µ
4 n
J
0
(
µ
n
)
e−
µ
2 n
F
o
+
ε KoLu 2(1 − Lu )
( )Pd2 J 0 Pd2 X − J1 Pd2
Pd2 Lu
J
0
J
1
Pd2
Lu
Pd2
Lu
X
e− Pd2Fo
−
( ) ∑−
Pd1 J0 2J1
Pd1 X Pd1
e− Pd1Fo
+
∞
1
n =1
1 −
µ
2 n
Pd1
−
(1 −
ε KoLu
Lu) 1 −
µ
2 n
Pd2
e−
µ
2 n
Fo
+
+
(1 −
ε KoLu
Lu) 1−
µ
2 n
Pd2
e−
µ
2 n
LuFo
J0(µn X )
J0 (µn )
для Lu=1 ( am = aq )
T (X , Fo) =
t(r,τ ) − t0 tмc − t0
= 2 1+
2(1
Po −X
2 пр
)
13 768
−
X
2 пр
−
1 8
Fo
−
1 16
X4
+1 32
X2
−
∑ ( )−
∞ n =1
2J0 µn X
µ
4 n
J
0
(
µ
n
)
e−
µ
2 n
Fo
−
( )Pd1 J 0 Pd1 X ( ) ∑J1 Pd1
( )e − Pd1Fo
+
∞ n =1
J0 µnX
1−
µ
2 n
Pd1
J0 (µn )
e
−
µ
2 n
Fo
(10) (11)
θ
(X
, Fo)
=
2
1 −
1 2
Pd2 Lu
J
0
J
1
∑P d 2
X
( )L u
Pd2
( )L u
e − Pd2Fo
∞
+
n =1
1 −
J0 µnX
µ
2 n
L
u
Pd2
J0
e
−
µ
2 n
F
o
;
µn
(12)
µn — последовательные положительные корни уравнения
J1(µ) = 0 .
(13)
Анализ полученных решений (1) и (12) показывает, что благодаря
быстрому увеличению абсолютной величины последовательного ряда значений
характеристических корней µn (уравнение (13)),
µ1 < µ2 < K < µn < K ,
а, следовательно, быстрого уменьшения экспоненциальных сомножителей
e ,−
µ
2 n
Fo
бесконечные
суммы,
входящие
в
решения,
сходятся
достаточно
быстро.
Поэтому, начиная с определенного значения числа Фурье, с заранее заданной
степенью точности из всего ряда разложения можно удовлетвориться одним-
двумя первыми членами. Такое упрощение имеет большое практическое
значение из-за существенного сокращения объема расчетной работы и
возможности представления общих решений в удобной для практического
применения форме.
Так, формулой удобной для инженерных расчетов полей потенциала
влагопереноса, будет следующее выражение, полученное из (12) при n = 1:
u(r,τ u0
)
=1−
2
Pd2 Lu
J
0
J
0
Pd2
X
Lu
Pd2
Lu
e− Pd2 Fo
+
1−
2J0 (µ1X )
µ12 Lu Pd2
J
0
(
µ1
)
e − µ12 Fo
.
(14)
Среднее значение безразмерной величины потенциала переноса влаги для
упрощенного выражения Ө(Х, Fo) при n = 1
∫ ( )θ ( F o ) =
u 0 − u (τ )
1
=2
X θ ( X , Fo)dx = 2
1 − e− Pd2Fo
u0 0
.
(15)
Ограничиваясь одним первым членом ряда, получаем приближенно
формулу для расчета времени, необходимого для достижения нагреваемым
телом определенного влагосодержания:
τ= θ
2k2 .
(16)
Среднее значение безразмерного потенциала температуры для
упрощенного выражения Т(Х, Fo) при n = 1
∫T ( F o ) =
1
2
XT (X , Fo)dx =
t(x,τ ) − t0
=
( )0 t мс − t0
Po 1 − X пр
2 + 2Pd1Fo
.
(17)
Для нахождения времени доведения продукта до необходимой средней
объемной температуры, запишем в следующем виде:
τ
=
2 k1
1
(t мс
−
t0
)
t
−
t0
−
w0R 2
λq
1 −
λ2 пр
R2
.
(18)
Выводы
Поставлена и аналитически решена задача совместного тепло- и массопереноса при инфракрасном нагреве колбасного батона цилиндрической формы.
Полученные решения позволяют рассчитать поля температуры и влагосодержания, усредненные значения соответствующих потенциалов переноса, темп нагрева, расход тепла в процессе теплообмена, а также получить формулы, удобные для инженерных расчетов.
Разработанные аналитические решения дают возможность определить время, необходимое для достижения нагреваемым продуктом определенной температуры и влагосодержания.
Обозначение
t(r, τ) — температура, К, °C; t0 − начальная температура; tмс − максимальная
температура
среды
в
камере;
∆t =tмс −t0;
T(X,
Fo)
=
(t(x,τ)
−t0) (tмс
−t0)
–
безразмерная
температура; u(r, τ) – влагосодержание, кг влаги/кг абс. сух. вещества; и0 −
начальное влагосодержание;
θ(X,
Fo)
=
(u0
−и(r,τ)) u0
–
безразмерное
влагосодержание;
r
–
текущая
координата, м; Х = r/R – безразмерная координата; R – радиус поперечного
сечения цилиндрического колбасного батона, м;
τ – время, с;τ1 − время окончания первого периода термической обработки
– подсушки; aq − коэффициент температуропроводности, м²/c; ε – коэффициент
фазового перехода (0 < ε < 1); ρ – удельная теплота фазового перехода, Дж/кг;
cq − удельная теплоемкость, Дж/(кг·K); am− коэффициент потенциало- (влаго-)
проводности, м²/c; λm − предельная глубина проникновения инфракрасного
излучения в материал, м; λq − коэффициент теплопроводности, Вт/(м·К);
q (τ ) = k1cqγ 0 R∆te −k1τ + Rk2 ργ 0u0e− k2τ – экспериментально найденная
плотность теплового потока на поверхности тела, Вт/м²; q m (τ ) = k 2γ 0u 0 R e − k2τ
– экспериментально установленная плотность потока массы вещества через
поверхность тела, кг/(м²·c); γ 0 – плотность абсолютно сухого вещества, кг/м³;
k 1 - коэффициент, характеризующий убывание температуры поверх-
ности тела по экспоненциальному закону, 1/с; k 2 - коэффициент сушки, 1/с; Fo = aqτ R2 − критерий гомохронности (число Фурье); Lu = am aq − критерий (число Лыкова) взаимосвязи массо- и теплопереноса; Ko = r∆u (cq∆t) − число
Коссовича; Pd1 = k1R2 aq , Pd1 = k2R2 aq - числа Предводителева (теплообменное и массообменное); J0(z) и J1(z) - функции Бесселя от вещественного аргумента z первого рода и нулевого порядка соответственно; w0 − мощность ИК-источника, Вт/м³.
Литература
1. Гинзбург А.С. Основы теории и техники сушки пищевых продуктов. – М.: Пищевая промышленность, 1973. – 528 с.
2. Гинзбург А.С. Инфракрасная техника в пищевой промышленности. – М.: Пищевая промышленность, 1966. – 376 с.
3. Рогов И.А., Некрутман С.В. СВЧ и ИК-нагрев пищевых продуктов. – М.: Пищевая промышленность, 1976. – 210 с.
4. Ильясов С.Г., Красников В.В. Физические основы инфракрасного излучения пищевых продуктов. – М.: Пищевая промышленность, 1978. – 359 с.
5. Лыков А.В. Тепло и массообмен в процессах сушки. – М.: Госэнергоиздат, 1956. – 464 с.
6. Белобородов В.В., Вороненко Б.А. Решение задачи нагрева тел в электромагнитном поле сверхвысоких частот //Журнал прикладной химии АН СССР, «Наука» - Ленинградское отделение, № 10, 1984. – С. 2276-2282.
7. Вороненко Б.А., Пеленко В.В., Стариков В.В. Постановка задачи тепломассопереноса процесса горячей сушки рыбы // Межвузовский сборник научных трудов «Ресурсосберегающие технологии и оборудование пищевой промышленности», СПб, СПбГУНиПТ, 2006. – С. 71- 75.
8. Стариков В.В., Вороненко Б.А. Применение ИК-нагрева при копчении // Межвузовский сборник научных трудов «Теория и практика разработки и эксплуатации пищевого оборудования», СПб, СПбГУНиПТ, 2007. – С. 12-17.
Analytical solution of a combined heat and mass transfer problem for infrared heating of sausage goods
B.A. Voronenko DSc, V.V.Pelenko DSc, V.V.Starikov graduate student
Saint-Petersburg State University of Refrigeration & Food Engineering
Use of infrared radiant energy is an effective way to enhance heat treatment of food materials, allowing both to curtail heat treatment period and to upgrade finished product quality. On the basis of prerequisites a mathematical model of combined heat and mass transfer is proposed for heat treatment of a long cylinder sausage loaf in an IR-heated smoker.
Key words: IR-heating, heat and mass transfer, mathematical model, differential equation set, initial and boundary conditions.