УДК 664.03
Аналитическое решение уравнений теплопроводности для абразивной обработки пищевого сырья
Д-р техн. наук Г. В. АЛЕКСЕЕВ, Б. А. ВОРОНЕНКО, В. А. ГОЛОВАЦКИЙ
Санкт-Петербургский национальный исследовательский университет ИТМО Институт холода и биотехнологий 191002, Санкт-Петербург, ул. Ломоносова, 9
В работе приводится математическое описание теплового процесса при первичной обработке пищевого сырья на стадии очистки. Полученное аналитическое решение может быть предложено для инженерных расчетов процесса теплового воздействия на пищевое сырье при его переработке.
Ключевые слова: аналитическое решение, уравнение теплопроводности, пищевое сырье, стадия очистки.
Analytical decision of the equations heat-conducting for abrasive
processing food cheese
G. V. ALEKSEEV, B. A. VORONENKO, V. A. GOLOVATSKY
National Research University of Information Technologies, Mechanics and Optics Institute of Refrigeration and Biotechnologies 191002, St. Petersburg, Lomonosov str., 9
The paper presents a mathematical description of the thermal process in primary treatment of food raw material in the purification step. The resulting analytical solution can be offered for engineering calculations, process heat influence on food raw material for reprocessing.
Keywords: analytical solution, heat equation, food raw materials, purification step.
В ряде работ [1, 2] авторами показано, что обрабатываемые на стадии очистки пищевые материалы (картофель, крупа, зерно) подвергаются значительным температурным воздействиям, влияющим на качество конечного продукта и производительность аппаратов. Однако многие тепловые процессы при первичной обработке сырья не изучены.
В значительной части оборудования, применяемого в пищевой промышленности, используется абразивный метод воздействия на объект. При этом предполагается непрерывность контакта рабочих органов с материалом. Актуальной является проблема получения и аналитического исследования
математических моделей, описывающих влияние термических процессов на пищевое сырье на указанном этапе его обработки.
Тепловой эффект при снятии кожуры овощей или оболочки зерновых может быть описан следующим образом при рассмотрении непрерывного контакта обрабатываемого продукта с абразивом при вращении продукта или абразива.
Клубень картофеля (или зерно) можно представить в виде тела основной геометрической формы - шара, состоящего из двух сферических тел: кольца оболочки (кожуры) и внутреннего шара (собственно ядра объекта) (см. рис.)
Рис. OR2 = R2 -радиус внутреннего шара; R2R] -радиус внешнего шара
Теплофизические свойства этих двух составных тел различны.
Клубень вращается со скоростью v = 2;r(^ ^2) ^ где р^=()1у - толщина
оболочки, Т- период вращения. Вследствие трения абразива о картофель (зерно) температура наружной поверхности тела - составного шара -изменяется по экспериментально найденному закону h (~R\ >т) = h - (?2 - h) ®ф( ■-кт). (1)
Очевидно, что для достаточно продолжительного времени обработки пищевого продукта температура поверхности, а, следовательно, и всего тела не может превысить максимального значения t2.
Так как толщина оболочки значительно меньше радиуса шара (картофеля, зерна), т.е. ORl = Rl « OR2 = R2 (длина полоски кольца /» /( ), то для описания температурного поля оболочки можно использовать уравнение теплопроводности для неограниченной пластины. Кроме того, принимаем в первом приближении следующее допущение: толщину оболочки считаем
неизменной. Таким образом, математически задачу нагрева вращающегося сферического тела можно сформулировать следующим образом: требуется решить систему дифференциальных уравнений теплопроводности
Э^(г,т) d2t1(r,r) dtl
дт
- а
дг2
-V-
дг
(2)
(г > 0, -i?j < г < 0); dt2(r,T)
<Эг
= а
^ d2t2 (г, г) 2 dt2 (г, г) Л
Sr
г дг
(3)
(г > 0, 0<r<R2);
при начальном условии
?i(r,0) = r2(>,0) = ?i, (4)
граничном условии первого рода (1) на внешней границе составного шара, граничном условии четвертого рода (равенстве температур и тепловых потоков) на границе оболочки и ядра *i(0 ,T) = t2(0,T);
-к
5^(0,г) =_^ dt2(0,r)
дг
дг
(5)
(6)
условии симметрии
dt2{R2,r)
дг
= 0
и физической ограниченности температуры в центре шара
t2(R2,0 < 00
Аналитическое решение краевой задачи (1) - (8) методом обратного интегрального преобразования Лапласа позволило получить распределение полей температур в составном шаре в следующем безразмерном виде:
гг Г V 17 Л tl(r’T)~tl ! <KPd)
TY (X, Fo) = -X- = 1 - exp
<p(Pd)
(7)
(8)
t2 tx
Pe
— (1 + X)- PdFo
V 2 .
+
00 A
+ X—exP „=1 <Pn
Pe i
— (\ + X)-(^ +
f т>„\2 \
Pe
v2 у
Fo
(9)
TJXJ’o)
exp
t2(rx)-ti ^
Pd
Pe
Pd sin( y]ka Pd X)
12 ^1
yjkaPdX(p{Pd)
Pe \
— -PdFo + /
v 2 J JkaX^
кг Am sin(V^(//„ + 0т) )X)
2^---------*-----------------exP
<Pn
2 Л
где
jun - последовательные положительные корни характеристического уравнения
tq/u = -
Р .
п „
Pe
inn <!ЛП < (2п +1)—);
(11)
^(Рй?)=
^Pe^
V 2 у
-Pd ch
f Pe^
V z У
-Pd X
\ ( Pe
-----s/z
2
V
^Pe^
V 2 у
-PdX
1
2
2
2
2
<p{Pd) -
гРел
L 2 J
■ Pd ch
\ f
Pe ,
н----sh
2
V
гРел
L 2 J
Pd
A, = 2^пМ (jun cos(//„X) - Ц- sin( дйХ);
■*721 \f~n
COS jLln
4,2 =
<Pn =
2 цп2 Pd
cos//„
f T>0\2\
Pe
V2y
(Pt/ - {/лп
( Pe^
V2y
' 2 9
)(1+!T<A
Pe
fPe^
v 2 у
));
Проведенный численный эксперимент показал адекватность решений (9) и (10) реальным процессам.
Разработанная математическая модель в виде аналитического решения соответствующей краевой задачи теплопроводности дает возможность прогнозировать и управлять температурным полем тела сферической формы (клубня картофеля, зерна, крупы), предотвращая перегрев продукта, и тем самым влиять на его качество.
Полученное решение позволяет решить обратную задачу по нахождению времени, необходимого для достижения нужной температуры в любой точке обрабатываемого пищевого материала.
Обозначения
tx(r,T)~ температура оболочки тела, °(\К\
(t2(г,т)-температура ядра тела °С,К;
t\= tmn = const - минимальная температура, равная начальной температуре тела; t2 = /|гах = const - наибольшая допустимая температура тела;
г - текущая координата (по диаметру) обрабатываемого тела, м; т- время, с;
к- положительная постоянная, зависящая от скорости вращения
исследуемого тела, 1/с;
а1 и а2 - коэффициенты температуропроводности оболочки и ядра
соответственно, м2/с;
\ и /^-коэффициенты теплопроводности, вт/(мК);
vRx
Ре = —- - число Пекле; a
kR
Pd = —— - число Предводителева;
a,
2
2
2
2
2
2
^ а, т .
Fo = -hr - число Фурье;
X = -- безразмерная координата;
Ri
к
a
ах
а2
Список литературы
1. Алексеев Г. В.,Головацкий В. А. и др. Патент на изобретение № 2228795 от 20.05.2004. Устройство для измельчения пищевого сырья.
2. Головацкий В. А. Совершенствование процессов и аппаратов для переработки пищевого сырья. -СПб.: НИЭУиД, 2008.-123с.