Например, Бобцов

Решение задачи механической очистки пищевых сред

УДК 665+602.1:519
Решение задачи механической очистки пищевых сред
Д.т.н. Вороненко Б.А., д.т.н. Пеленко В.В., аспирант Поляков С.В.
Санкт-Петербургский государственный университет низкотемпературных и пищевых технологий
К.т.н. Марков В.Н.
Всероссийский научно-исследовательский институт жиров РАСХН
Представляет аналитическое решение интегрально-дифферинциального уравнения, являющегося математической моделью одного из процессов разделения двухфазных систем. Это решение может быть использовано, в частности, для инженерных расчётов механической очистки сырья при производстве кондитерского крема и первичной очистки нерафинированных растительных масел.
Ключевые слова: аналитическое решение, двухфазная система, механическая очистка, крем кондитерский, нерафинированное масло.
В технике широко распространен и используется для разделения двухфазных или многофазных систем (жидкость – твердое тело, газ – твердое тело, жидкость – жидкость, твердое тело – твердое тело, жидкость – газ – твердое тело и т.д.) гидромеханический процесс осаждения твердых (или жидких) частиц в жидкой (газовой) среде. Физические характеристики движения, например, твердых частиц в среде жидкости зависят от реологических свойств системы, от факторов, связанных с условиями обтекания [1]. В случае осаждения мелкодисперсных твердых частиц в жидкости или газе основной характеристикой процесса является скорость осаждения.
В молочной промышленности для механической очистки, являющейся составной частью водоподготовки при производстве крема кондитерского, применяется процесс отстоя, при котором частицы жира, как наиболее легкие, всплывают, выделяясь из молока. Физические закономерности, которым подчиняются процессы осаждения и отстоя, сходны и отличаются только направлением движения частиц, отделяемых от среды [2].
В производстве растительных масел важнейшим этапом очистки нерафинированных масел является первичная очистка – удаление примесей нежирового характера (мелких частиц жмыха, фосфолипидов, воскоподобных соединений) [3].
Движение твердых частиц в жидкости может быть описано с помощью уравнений Навье-Стокса.
При малых скоростях движения потока жидкости характерен ламинарный (слоистый) режим течения. При ламинарном движении слои жидкости скользят

один относительно другого не перемешиваясь. В условиях установившегося

(стационарного) движения скорость потока W при ламинарном режиме постоянна
в каждой его точке, т.е. W = f (x, y, z), где x, y, z - текущие координаты.

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

осаждении в жидкости частиц малых размеров. При значениях числа Рейнольдса

Re = W0dρ = W0d µν

< 1 (W0

- скорость потока, обтекающего частицу (вдали от нее);

d - диаметр частицы; ρ - плотность жидкости; µ - динамический коэффициент

вязкости; ν - кинематический коэффициент вязкости), уравнения движения

Навье-Стокса упрощаются и из них выводится в [1] формула полной

кинетической силы сопротивления, проявляющаяся только в движущемся потоке:

Fk = 6πµRW0 ,

(1)

где R — радиус шарообразной частицы.

Уравнение (1) выражает закон Стокса.

Полная сила сопротивления при обтекании частицы (шара) пропорциональна

импульсу (количеству движения) и площади лобового сечения f:

Fk

= C ρW02 2

f

(формула Ньютона).

(2)

Отсюда коэффициент пропорциональности:

C

=

2Fk ρW02

f

.

(3)

Учтя f = πR2 и Fk из (1), получим [1]:

C

=

2 ⋅ 6πRµW0 ρW02 πR 2

=

12µ ρW0 R

=

24 Re

=

ξ.

(4)

Уравнение (4) выражает закон сопротивления среды при ламинарном режиме

(Re < 1).

При этих условиях скорость осаждения частицы описывается формулой

Стокса:

Woc

=

d (2 ρтв − ρ)g
18µ

,

(5)

которая справедлива для области 10−4 < Reoc < 2 [1].

В (5) ρтв — плотность частицы.

Для нестационарных условий (неустановившемся движении) и при движении

твердых тел в жидкости близкой плотности формула Стокса (1) заменяется для

силы сопротивления движущемуся в жидкости телу формулой Буссинеска [4-8]:

Fсопр

(t)

=

−6πµRW

(t)



2 3

πρR3

dW (t) dt



 6πµR


R

t


dW

(τ)

πν 0 dτ





2R2

 W (o)

t − τ πνt 

(6)

где t — время.

(В [4] и в предыдущих изданиях монографии [5] под знаком интеграла
ошибочно указана скорость W , а не ускорение dW (τ) (или W '); решение

уравнения (6) дано только в [4] для ошибочного случая). Предполагается, что вдали от частицы жидкость неподвижна. В формуле Буссинеска первый член представляет собой стационарную
формулу Стокса (силу сопротивления Стокса) (1), второй – инерционную составляющую силы сопротивления, учитывающую присоединенную массу жидкости (шара) [Движение шара в жидкости можно рассматривать как происходящее в пустоте, если только к массе шара присоединить дополнительную массу, равную половине массы жидкости в объеме шара], третье слагаемое – так называемая сила Бассе. Сила Бассе обычно включается в выражение для силы сопротивления движущейся частице, если сила межфазного взаимодействия зависит от предыстории движения частиц при их нестационарном движении. Иногда слагаемыми типа силы Бассе пренебрегают [6], что значительно упрощает описание движения дисперсной частицы, однако при этом величина относительной погрешности при вычислении скорости осаждения твердой частицы может достигать 20% [8].

Из формулы Буссинеска (6) следует, что при t = 0 и W (o) ≠ 0 сила сопротивления становится бесконечной: невозможно шару, погруженному в вязкую жидкость, мгновенно сообщить конечную скорость. Поэтому следует

принять W (o) = 0 . В случае шара, приведенного импульсивно из состояния покоя в состояние

поступательного, прямолинейного и равномерного движения со скоростью W (o) ,

из (6) следует величина силы сопротивления, равная [5]:

Fсопр

=

6πµRW0

1 

+

R  ,
πνt 

(7)

из которой при t → ∞ вновь получается формула Стокса [1].

При вертикальном падении частицы в форме шара уравнение его движения

будет (основной закон динамики – второй закон Ньютона):

F cопр

+

mg

+

FA

=

m

dW dt

или

(8)

4 3

πρ1R

3

dW dt

= F сопр

+

4 3

π(ρ1



ρ)R 3

g

,

(9)

где FA — выталкивающая сила (сила Архимеда); ρ1 — плотность материала шара;

g — ускорение свободного падения; mg — сила тяжести.

Проектируя (6) на направление движения с учетом (8), получим интегро-

дифференциальное уравнение [4]:

dW (t) + α2W (t) +



t


dW

(t

)



dt π 0 dτ

dτ t−τ

=

f0 ,

(10)

где введены следующие обозначения:

α2

=



2R

2

 

ρ1

+

ρ  2

;



=

9 µρ

2R 

ρ1

+

ρ 2

 

;

f0

=

g(ρ1 − ρ) .

ρ1

+

ρ 2

Для решения уравнения (10) применим один из методов операционного

исчисления – метод интегрального преобразования Лапласа. Он заключается в

том, что изучается не сама функция (оригинал), а её видоизменение

(изображение). Это видоизменение – преобразование – производится при помощи

умножения на некоторую экспоненциальную функцию и интегрирования в

определённых пределах.
Преобразование Лапласа функции y = f (τ) (оригинала) состоит в умножении

ее на e−sτ и интегрировании в пределах от 0 до ∞ [4, 9]:


∫ f (τ)e−sτdτ = F (s),

0
( f (τ) L→ F(s)

где s = ξ + iη - некоторая комплексная величина; i = − 1 ; F (s) - изображение функции f (τ).

Преобразованное уравнение (например, дифференциальное, интегро-

дифференциальное) в области изображений значительно проще исходного, часто

это обычные алгебраические уравнения. Затем находят искомую функцию в

изображениях и затем совершают обратный переход (по обратному

преобразованию Лапласа) к оригинальной искомой функции.

Если

W

(t

)

L
⋅ →VL

(s

)

,

то

dW (τ)
dt

L
⋅ →

sWL

(s

)



W

(0)

;

поскольку W (0) = 0 ,

то

dW (t
dt

)



sWL

(s

)

=

sWL

;

t


dW

(τ)



0 dτ

dτ t−

τ

L
⋅ →

sWL



1= s

[по теореме умножения изображений (теореме Бореля):

sW L

t
∫ ϕ(t) f

(t



τ)dτ

=

t
∫ ϕ(t



τ) f

(τ)dτ

=

f



ϕ

L
⋅ →

F

(s

)



Φ(s

)

,

00
где f ∗ τ - свёртка непрерывных функций f (t) и ϕ(t); F (s) и Φ(s) - изображения функций f (t) и ϕ(t)].

Таким образом, после применения к уравнению (10) преобразования Лапласа,

получим изображающее уравнение:

sWL + α2WL + 2β

sWL =

f0 , s

(11)

откуда задача сводится к разысканию искомой функции для изображения

скорости:

( )WL

=

s

s

+

f0 2β s

+

α2

.

(12)

Обозначим s = x и разложим дробь (12) на простейшие

неопределенных коэффициентов. В результате получим:

( )WL

=

x2

x2

f0 + 2βx + α2

=

x2 (x



f0
x1 )(x



x1 )

=

=

f0 α4

 

α

2

 x2



2β x

+



2β2 2

+ 2β β2 −

β2 − α2 +
α2 (x − x1)

α2

+

2β2 + 2β β2 − α2 − α2
2 β2 − α2 (x − x2 )

 , 

где x1 и x2 — корни квадратного трехчлена x2 + 2βx + α2 :

методом (13)

x1 = −β − β2 − α2 ; x2 = −β + β2 − α2 .

Здесь возможны два случая.

При β > α , т.е. ρ > 1,6ρ1, оба корня х1 и х2 действительны и отрицательны,

плотность жидкости больше плотности тела (газовые пузырьки).

В этом случае

( ) ( )WL(s) =

f0 α4

 

α

2

 s



2β s

+



2β2 + 2β β2 2 β2 − α2

− α2 + s − x1

α2

+

2β2 2

+ 2β β2 −

β2 α2

− α2 s−

− α2 x2

 . 

(14)

По обратному преобразованию Лапласа (используя соответствующие

таблицы изображений функций [9, 10]) получаем окончательное решение для

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

( )W

(t )

=

f0 α4

 α2 



2β + − 2β2 + 2β β2 − α2 + α2 

πt

2 β2 − α2



1 πt

+

x1e x1t erfc



x1

t  + 

( )+ 2β2 − 2β β2 − α2 − α2 

2 β2 − α2



1 πt

+

x2e x2terfc



x2

t

 

=

( )=

f0 α4

 

α

2



2β − α2 + 2β πt 2

β2 β2

− α2 − α2



2β2

∑2 (− 1)i xie xiterfc
i =1



xi

t

 

,

(15)

где

erfz = 1 − erfz =

2


∫e



z

2

dz

;

πz

erfz

=

2

z


e



z

2

dz

=

Φ(z)

π0

— функция ошибок Гаусса [11].

(16)

Разложение в ряд функции Φ(z) для малых значений z = −xi t (i = 1, 2)
будет:

Φ(z) =

2

z


∑∞ (− 1)n

z 2n dz

=

π 0n=0

n!

2 π

∑∞ (−
n=0

1)n

z 2n +1
n!(2n +

1)

=

2 π



z 1



z3 1!3

+

z5 2!5



... ,

(17)

для больших значений z:

Φ(z) =1−

1 π

e−z2

 1 z



1 2z3

+

1⋅3 22 ⋅ z5



1⋅3⋅5 23 ⋅ z7

+ ... .

(18)

При ρ < 1,6ρ1, т.е. α > β (плотность жидкости меньше плотности тела) выражение для скорости осаждения частиц имеет вид (соответствует

поставленной задаче удаления нежировых примесей):

( )W (t) =

f0 α2



  1

+

−α
e

2

1− 2

β α

 2 

t



 

   ⋅

1



2

β α

2 

  α 

1



 

β α

 

2

A1

+

2β α2





A2

 

cos







t

+

( )

+

  

2β α2





A1



1



2

β α

 

2

α

1



 

β α

 

2





A2

 

sin







t  ,

(19)

где

( )A1 =

2 π

e

−β

2

t

  

α

1



 

β α

2 

ϕ1

+

βϕ2

  



α

1



 

β α

 

2

erfc

β

t

;

( )A2 =

2 π

e−β2t

  

βϕ1



α

1



 

β α

2 

ϕ2

  



βerfc

β

t

;

α 1− β 2 t

( )ϕ1 =

α


e y2 sin 2β t y dy ;

0

α 1− β 2 t

( )ϕ2 =

α


e y2 cos 2β t y dy .

0

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

Преобразованное интегро-дифференциальное уравнение решено методами математической физики. Полученное аналитическое решение является математической моделью одного из процессов разделения двухфазных систем, что может быть, в частности, использовано для инженерных расчетов механической очистки сырья при производстве кондитерского крема и первичной очистки нерафинированных растительных масел.
Список литературы
1. Романков П.Г., Курочкина М.И. Гидромеханические процессы химической технологии. –Л.: Химия, 1982. –288 с.
2. Кук Г.А. Процессы и аппараты молочной промышленности. –М.: Пищевая промышленность, 1973. –768 с.
3. Вороненко Б.А., Марков В.Н., Кунилова Т.М. Аппаратурное оформление процесса первичной очистки растительных масел// Межвузовский сб.науч.трудов «Теория и практика разработки ресурсосберегающего пищевого оборудования», СПб: ГОУ ВПО «СПбГУНиПТ», 2006. –С.51-56.
4. Лурье А.И. Операционное исчисление и его приложение к задачам механики. – М. –Л.: Гостехиздат, 1950. –432 с.
5. Лойцянский Л.Г. Механика жидкости и газа. Учеб. для вузов. –7-е изд. – М.: Дрофа, 2003. – 840 с.
6. Протодьяконов И.О., Чесноков Ю.Г. Гидромеханические основы процессов химической технологии. –Л.: Химия, 1987. – 360 с.
7. Ландау Л.Д., Лифшиц Е.М. Механика сплошных сред. –М.: Гостехиздат, 1953. – 788 с.
8. Наумов В.А. Влияние силы Бассе на динамику твердых частиц в придонном слое. // Сб.науч.трудов «Теоретические и практические аспекты применения инженерной физико-химической механики с целью совершенствования и интенсификации технологических процессов пищевых производств», Моск.госуд.университет прикладной биотехнологии, -М.: 2002. – С.257-260.
9. Мартыненко В.С. Операционное исчисление.–К.:Высш.шк.,1990.–359 с. 10. Лыков А.В. Теория теплопроводности. –М.: Высшая шк., 1967. –600 с. 11. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и
произведений. –М.: Физматгиз, 1963. –1100 с.