CC BY
16
УДК 653.093
Определение градиента скорости сдвига продукта в емкостном оборудовании на основании полуэмпирической теории турбулентного переноса
Канд. техн. наук Б.Л. НИКОЛАЕВ СПбГУНиПТ
On the basis of the half-empirical theory of turbulent carry the author offers the mathematical model, allowing to count the velocity gradient of shear in tanks with mixing devices.
При обработке вязких структурированных продуктов в емкостном оборудовании главным фактором, влияющим на разрушение их структуры, является градиент скорости сдвига продукта. В ряде работ предлагаются различные зависимости для определения градиента скорости [1-3,7]. Однако теоретические обоснования приводимых формул либо отсутствуют, либо содержат большое количество весьма спорных допущений.
Если мощность, расходуемую на перемешивание, и коэффициент теплоотдачи в имеющемся оборудовании или его модели можно определить опытным путем [7], то непосредственно измерить градиент скорости не представляется возможным. Поэтому он может быть только рассчитан.
В данной работе автор на основании полуэмпирической теории турбулентного переноса предлагает научно обоснованную методику расчета градиента скорости сдвига продукта у в емкостных аппаратах со скребковыми перемешивающими устройствами.
Подавляющее большинство пищевых продуктов представляют собой неньютоновские среды. Для таких сред существует общеизвестная формула
Шф = *т"4. (1)
где Цэф - эффективная вязкость, Па-с; к - коэффициент Оствальда, Па с; т - показатель неньютоновского поведения.
И.В. Доманский и В.Н. Соколов в работе [4] приводят уравнение
Ео = Й^У2, (2)
где Е0 - диссипация энергии у твердой стенки, Вт/м3.
В этой же работе приводится формула
= С)
где х ~ безразмерный коэффициент;
Е - средняя по объему жидкости диссипация энергии, Вт/м3.
Физический смысл коэффициента % авторами не раскрывается.
Из работы [4] известно, что
4И
£ = —= - 2 V nD Н
(4)
где И-мощность, расходуемая на перемешивание, Вт; V - объем резервуара, м3;
О - диаметр резервуара, м;
Н - высота резервуара, м.
Из уравнений (2) - (4) получаем
ЛN
ЦэфУ =х
kD2H
(5)
Последовательно преобразуя уравнение (5), имеем
4NX4
У2 =
^эф =
4jVx4
(6)
(7)
у2лй2Н '
Подставив в уравнение (7) формулу (1) и произведя серию последовательных преобразований, получим:
*г' = , ,
у nD Н
х4 =
х =
4ІУХ4
y2nD2h
4N
кІЇНку 4JV
т+1
\0,25
(8)
(9)
(10)
.т+1
_4Л^_ кй*Нк'
0,5 0.5/я—I
у =
4Л'Х< '
тФгНк
т+1
(11)
(12)
И.В. Доманский и В.Н. Соколов в работе [4] приводят следующее уравнение:
а V Рг
V ’
(13)
где а - коэффициент теплоотдачи на стороне продукта, Вт/(м2К);
X - коэффициент теплопроводности продукта, Вт/(мК);
V - коэффициент кинематической вязкости продукта, м2/с;
р - плотность продукта, кг/м3;
Рг = с Иэф/Х. - коэффициент Прандтля; с - теплоемкость продукта, Дж/(кг-К); у- средняя безразмерная разность температур. Известно [7], что
т-1
(14) Р Р
Подставив в уравнение (13) уравнения (10) и (14) и внеся х под знак корня, получим
аку
т—1
Рг
Хр
(15)
кй1Н!сут+1 Ь( Е
ч 4 N р‘ )
Подставив в уравнение (15) формулу (4), будем иметь
оку"-' Рг
Хр
,0,25 — •
V” »к
4Ы р2п02Н
(16)
Проведем серию последовательных преобразований уравнения (16):
аАу”"1 Рг
Хр
0,25 —
¥
акуг
Хр
Рг
к0'5 у
(17)
(18)
.0,5
ак у
Рг
»0,5
(19)
Хр - V
Подставив в уравнение (19) значение коэффициента Прандтля и произведя преобразования, получим:
^ СЦэф >,р0,5 уХ
Ок0-У'т~' сЬ/"-'
(20)
.0,5
Р V
а ск°'5у°'5т
Л5
Г =
ац/
(21)
(22)
(23)
уср0’5*0’5; '
И.В. Доманский и В.Н. Соколов [4] на основании по-луэмпирической теории турбулентного переноса предлагают формулу для расчета средней безразмерной разности температур у:
Пм пм
Л
¿Т]
1 8
о О---------1---
Рг V
(24)
где т\м= и*у^ - максимальное безразмерное расстояние от стенки;
м* - динамическая скорость, м/с; ум - расстояние от стенки до точки с наименьшей или наибольшей (в зависимости от направления теплового потока) температурой жидкости, м;
5 - коэффициент турбулентного обмена, м2/с; г| = и* у/\ - безразмерное расстояние от стенки; у - расстояние от стенки в направлении радиуса кривизны поверхности, м.
В работе [8] на основе гипотезы о вязком затухании турбулентных пульсаций у твердой стенки приводится уравнение, единое для всех значений т|:
0,5
1 + ^ 1 + 0,64ц3
1-ехр
26
-1
(25)
V 2
Для емкостного оборудования точка с наименьшей или наибольшей температурой продукта находится в центре резервуара. Поэтому
У. = 0/2. (26)
Динамическую скорость определим по формуле иф = со(0/2) = пйп, (27)
где со - угловая скорость вращения мешалки, с-1; п - частота вращения вала мешалки, с-1.
С учетом формул (26), (27) и (14) распишем в явном виде формулу для определения г)м:
Лм = {[пОп />/2\/Ы} р = [п&п р]/(2/су™-'). (28)
Подставив формулы (25) и (28) в уравнение (24) и расписав в явном виде критерий Прандтля, получим
я02пр яР'пр
2ктГ' .
г г________________________________________________
_ I [ Я/(с<:у"-,)-ЦИ-[И-0,64т?г(1-ехр(-п/26))гГ,}/2-1 по, ¥ л02пр/(2ку~') '
Подставив (29) в (23), окончательно получим
.(30)
В уравнении (30) все величины, кроме градиента скорости у, известны. Таким образом, мы имеем одно уравнение с одним неизвестным. Входящий в уравнение (30) двойной интеграл в явном виде взять невозможно. Соответственно получить из уравнения (30) явную зависимость для расчета у не представляется возможным. Поэтому уравнение (30) следует решать методом последовательных приближений. При современном уровне развития компьютерной техники это не представляет трудностей. Например, в среде MathCAD для решения уравнения (30) можно использовать стандартную процедуру root, дающую на выходе готовое значение у.
Таким образом, нами была разработана методика, позволяющая рассчитывать градиент скорости сдвига продукта у в емкостном оборудовании со скребковыми перемешивающими устройствами при следующих известных величинах: а, п, D, р, с, X, к, т. Предложенный способ был получен на основе полуэмпирической теории турбулентного переноса посредством строгих математических преобразований, выполненных без ка-ких-либо допущений и аппроксимаций. Поэтому погрешность при определении у по формуле (30) будет определяться только тем, с какой точностью замерены входящие в нее величины а, п, D, р, с, X, к, т.
В качестве примера рассмотрим, как использовать описанную выше методику определения градиента скорости сдвига в резервуаре с D = 3 м и количеством скребковых лопастей z = 2 при изменении п от 0 до
*DV
"Г ____________________________________dr}
1 J
г_ . - ^ 1(скуя~') + {1 + [1 + 0,64г}2 (1 - ехр(-Т) / 26))г ]°5} / 2 — 1
2 с-1 для пломбира шоколадного, реологические характеристики которого были исследованы нами ранее [5].
При температуре продукта 15 °С пломбир имеет следующие реологические характеристики: к=0,253 Па с; т = 0,787.
Расчеты будем производить для запатентованной нами емкости с многоконтактным очищающим устройством [6]. Как показали проведенные нами эксперименты на модельных средах и реальных продуктах, для данной емкости коэффициент теплоотдачи а на стороне продукта составил 80 % от теоретического. Расхождение между теоретическим коэффициентом, умноженным на 0,8, и реальным не превышало 6 %, что находится в пределах погрешности измерений. Такое близкое совпадение значений теоретического и реального коэффициентов теплоотдачи можно объяснить тем, что в предложенной нами модели пристенные слои счищаются полностью, и благодаря конструкции мешалки происходит полное перемешивание счищенных слоев с основной массой продукта. То есть два основных фактора, уменьшающих реальный коэффициент теплоотдачи по отношению к теоретическому, в предложенном агрегате сведены к нулю.
Теоретический коэффициент теплоотдачи в емкостном оборудовании со скребковыми перемешивающими устройствами рассчитывается по известной формуле
а = (2/-Утг )(срХяг)0,5. (31)
Для получения реальных значений коэффициента теплоотдачи в емкости правую часть уравнения (31) нуж-
Зависимость градиента скорости у от частоты вращения п
но умножить на 0,8. Таким образом, значения коэффициента теплоотдачи в [6] будут рассчитываться по формуле (32): а = 0,8(2/у[к jicpXnz)0'5. (32)
Значения, а, рассчитанные по формуле (32), и будут подставляться в уравнение (30).
На рисунке представлена зависимость (30) градиента скорости сдвига продукта от частоты вращения мешалки.
Как видно из рисунка, предлагаемая методика расчета градиента скорости сдвига достаточно точно отражает реальное положение вещей.
Список литературы
1. Барановский Н.В., Коваленко Л.М., Ястребенец-кий А.Р. Пластинчатые и спиральные теплообменники.-М.: Машиностроение, 1973.
2. Бегачев В.И., Брагинский JI.H., Глухов В.П. и др. О моделировании полимеризаторов со скребками // Тепло- и массообмен в неньютоновских жидкостях. — М.: Энергия, 1969.
3. Глуз М.Д., Павлушенко КС. О мощности, затрачи-
ваемой на перемешивание скребковыми мешалками //Теоретические основы химической технологии. 1969. Т. 3. № 5.
4. Доманский И.В., Соколов В.Н. Обобщение различных случаев конвективного теплообмена с помощью полуэмпирической теории турбулентного переноса // Теоретические основы химической технологии. 1968. Т. 2. № 5.
5. Николаев Б.Л. Математическая обработка реологических характеристик смесей мороженого // Материалы II Междунар. науч.-техн. конф. «Прогрессивные технологии и оборудование для пищевой промышленности». - Воронеж, 2004. 4.2.
6. Патент № 2002129724/20. Емкость с многоконтактным очищающим устройством / Б.Л.Николаев. Опубл. 20.06.2003, Бюл. ФИПС № 17.
7. Стренк Ф. Перемешивание и аппараты с мешалками. -Л.: Химия, 1975.
8. Vermeulen Т., Williams G., Lang Lois G. Interfacial area in liquid-liquid and gas-liquid agitation // Chem. Eng. Progress. 1955. V. 51.
