Аналитическое решение задачи совместного тепломассопереноса при инфракрасном нагреве колбасных изделий

УДК 637.523.27:536.25

Аналитическое решение задачи совместного тепломассопереноса при инфракрасном нагреве

колбасных изделий

Д.т.н. Вороненко Б.А., д.т.н. Пеленко В.В., аспирант Стариков В.В.

Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий

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

Ключевые слова: ИК-нагрев, тепломассоперенос, математическая модель, система дифференциальных уравнений, начальные и граничные условия.

Технологический процесс тепловой обработки колбасных изделий как сложных капиллярно-пористых коллоидных тел, состоит из стадий подсушки, обжарки и варки. Интересующая нас первая стадия (подсушка) заключается в прогреве, главным образом, поверхности продукта (Н.С. Михеевой показано, что углубление поверхности испарения имеет место с начала процесса сушки [1]) в среде с низкой относительной влажностью (до 10%). При подсушке с поверхности колбасного изделия удаляется влага смачивания (механически удерживающаяся вода), что способствует равномерной прокраске поверхности и диффузии в продукт коптильных веществ при последующей обжарке.

Как известно [2-4], использование энергии электромагнитного поля инфракрасного (ИК) диапазона является одним из эффективных путей интенсификации тепловой обработки пищевых материалов, позволяющим значительно сократить длительность тепловой обработки и повысить качество готовых изделий.

При составлении математической модели процесса термической обработки колбасных изделий в коптильных установках с ИК-нагревом конвективным теплообменом от окружающего воздуха пренебрегаем, причем, чем выше температура источника радиационного излучения, тем обоснованнее эта предпосылка [2, 8].

Лучистый поток, проникая в продукт, затухает по экспоненциальному закону [5].

При толщине продукта I >10 мм и высоком значении коэффициента поглощения лучистый поток быстро затухает по мере проникновения в продукт

[1], и можно считать, что вся энергия отдается поверхности, а в нагреве внутренних слоев не участвует. В этом случае, особенно для периода постоянной скорости сушки — стадии подсушки (от начала процесса до

времениТх), затухание лучистого потока, проникающего в продукт, может быть описано параболическим законом [6]:

г2 -I2

= -¥-1г

Принимаем колбасный батон за тело основной геометрической формы — бесконечный цилиндр, однородный и изотропный, теплофизические характеристики которого в процессе копчения остаются неизменными.

При ИК-облучении нагрев происходит с поверхности, а внутренним слоям колбасного батона теплота передается теплопроводностью и влагопровод-ностью. На основе многих экспериментальных данных [3] не будем учитывать процесс термодиффузии влаги, т.е. изменение влагосодержания за счет термовлагопроводности, и учтем также сравнительно медленное изменение влагосодержания со временем при ИК-нагреве [7].

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

Э ^ Г Э2 £ 1 Э £ ^ ер Э и г)

Jr = 3q

з 2 + + ;

V Эг г Эг) сч Эт CJ0 Э u

Э2 u 1 Эй Э г2 г Э г )

Эт

(0< т < т;0< г< R) при следующих условиях:

t

u

3t (0,т) = Эи (0,т) = 0; (6) Эг Эг

t

«

^ + 9.

Здесь: (2) — уравнение теплопереноса с внутренним источником тепла, определяемым зависимостью (1); (3) — уравнение влагопереноса; (4) и (5) — начальные условия, задающие равномерное распределение температуры и

влагосодержания в теле в момент начала процесса нагрева; (6) и (7) — условия симметрии и физической ограниченности температуры и влагосодержания; (8) — граничное условие второго рода, задающее поток теплоты через поверхность тела; (9) — граничное условие второго рода, описывающее влагообмен тела с окружающей средой.

Аналитическое решение задачи (2) — (9) получено методом интегрального преобразования Лапласа в следующем безразмерном виде:

Т

м с 0 ( V пр'

13 768

/

V

X2 -1 | Ра - — X4 + — X2 пр 8 ) 16 32

- I 2 0 {Уп:^ ) е-м\Ра

тП Jo

еКаЬи

2

/

л/Р^Л {4Рё2 X) VI

Jo

и

Ьи

X

Лл/Р^

Л

Ьи

- Рё, Ра

Э 2

[4Рё1 X) 2 Jl^jPd1

е~ рё1ра + I

п=1

1

еКаЬи

1 - М

( .2 \

1

Мп

Р±

е~М Ра +

\ " "2 У У

+ -

еКаЬи

~МП ЬиРа

М Л

V Рё2 У

)

J0

для Ьи=1

т(X, Ра) =

<(г,т)- <0

<мс - <0

= 2 Л + -

Ра

Ц 2 J0 ^ ) е-м2Ра

П=1 Мп J0

2

УРё!J0 ^урё! X) J1 )

768

пр

13 ._ Г X2 - !! Ра - — X4 + — X2 -

8 У 16 32

е

- Рё1Ра

+I

п=1

J0

í 2 Л

м1

е Ра

V Рё1 У

J0

в( X ,Ра) = 2

- 1 [Рё7

2 4 Ьи

J 0

Рё 2

X

V Ьи

J1

е - Рё2 Ра

+ I

Jo

\ /

' 1 М2Ьи > Рё

_еРа

J0

2 У

;

Мп —

последовательные положительные корни уравнения

Jl(М) = 0

л

г

е

\

\

/

1

п = 1

Анализ полученных решений (1) и (12) показывает, что благодаря быстрому увеличению абсолютной величины последовательного ряда значений характеристических корней т (уравнение (13)),

т <т <  <т< •••,

а, следовательно, быстрого уменьшения экспоненциальных сомножителей

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

Так, формулой удобной для инженерных расчетов полей потенциала влагопереноса, будет следующее выражение, полученное из (12) при п = 1:

и

= 1 - 2

и

Л

\Рй1 Ьи

X

Ьи

Л

Ьи

е~ Рё^Р° +

2Л

сТМ?Р°

т Ьи

Л

2 У

Среднее значение безразмерной величины потенциала переноса влаги для упрощенного выражения 0(Х, Бо) при п = 1

в( Р°) =

и г

и

и

= 21 Xв

'о о

Ограничиваясь одним первым членом ряда, получаем приближенно формулу для расчета времени, необходимого для достижения нагреваемым телом определенного влагосодержания:

в

т =

2 к„

Среднее значение безразмерного потенциала упрощенного выражения Т(Х, Бо) при п = 1

' t (хт)-Л =

температуры для

Т(Р°) = 21 ХТ (X, Р°) ёх =

Р°

О - X)

2 + 2 Рё 1 Р°

Для нахождения времени доведения продукта до необходимой средней

объемной температуры, запишем в следующем виде:

( \

т =

1

2 к.

1 \ м с 0 /

^о Х-

я.

я

2 Л

п р

у У

1

о

о

Г

1

Выводы

Поставлена и аналитически решена задача совместного тепло- и массопереноса при инфракрасном нагреве колбасного батона цилиндрической формы.

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

Разработанные аналитические решения дают возможность определить время, необходимое для достижения нагреваемым продуктом определенной температуры и влагосодержания.

Обозначение

1;(г, т) — температура, К, °С; ^ начальная температура; м максимальная

лУ т _ ОМ-О/ А_ г 1(уРо) _ /О -г) температура среды в камере; м /Км ™ - безразмерная

температура; и(г, т) - влагосодержание, кг влаги/кг абс. сух. вещества; и начальное влагосодержание;

6(Х О _ (и -и(гт))/

/и - безразмерное влагосодержание; г - текущая

координата, м; Х = г/Я - безразмерная координата; Я - радиус поперечного

сечения цилиндрического колбасного батона, м;

т - время, с;т время окончания первого периода термической обработки

- подсушки; ач коэффициент температуропроводности, м2/е; е - коэффициент фазового перехода (0 < е < 1); р - удельная теплота фазового перехода, Дж/кг;

Сд удельная теплоемкость, Дж/

проводности, м2/с; т предельная глубина проникновения инфракрасного

1 -

излучения в материал, м; д коэффициент теплопроводности, Вт/(мК);

д(т) _ к^чУоЯАíe-+ Якру,"ое"

экспериментально найденная

2- Чт (Т) _ к270 " 0 Яе

- к2т

плотность теплового потока на поверхности тела, Вт/м2; т 2 0 0

- экспериментально установленная плотность потока массы вещества через

поверхность тела, кг/(м2-с);7о - плотность абсолютно сухого вещества, кг/м3;

к 1 - коэффициент, характеризующий убывание температуры поверхности тела по экспоненциальному закону, 1/с; к2- коэффициент сушки, 1/с;

Ро _ аТ/2- 1и _ ат/а -

/Я критерий гомохронности (число Фурье); / ад критерий

Ко _ гАу(с-

(число Лыкова) взаимосвязи массо- и теплопереноса; / ( дА) число

Ж=, ж=

Коссовича; / 4 / 4 - числа Предводителева (теплообменное и

массообменное); и ^ - функции Бесселя от вещественного аргумента 7

первого рода и нулевого порядка соответственно; и°- мощность ИК-источника,

Вт/м3.

Литература

1. Гинзбург А.С. Основы теории и техники сушки пищевых продуктов. - М.: Пищевая промышленность, 1973. - 528 с.

2. Гинзбург А.С. Инфракрасная техника в пищевой промышленности. - М.: Пищевая промышленность, 1966. - 376 с.

3. Рогов И. А., Некрутман С.В. СВЧ и ИК-нагрев пищевых продуктов. - М.: Пищевая промышленность, 1976. - 21° с.

4. Ильясов С.Г., Красников В.В. Физические основы инфракрасного излучения пищевых продуктов. - М.: Пищевая промышленность, 1978. - 359 с.

5. Лыков А.В. Тепло и массообмен в процессах сушки. - М.: Госэнергоиздат, 1956. - 464 с.

6. Белобородов В.В., Вороненко Б.А. Решение задачи нагрева тел в электромагнитном поле сверхвысоких частот //Журнал прикладной химии АН СССР, «Наука» - Ленинградское отделение, № 1°, 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.