close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
УДК 664.696.9
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕЧЕНИЯ АНОМАЛЬНО ВЯЗКИХ ЖИДКОСТЕЙ
ПРИ ПОЛУЧЕНИИ КОЭКСТРУДИРОВАННЫХ КОМПОЗИЦИЙ
А.А. Шевцов, Л.И. Лыткина, Е.А. Острикова
Разработана математическая модель движения аномально вязкой жидкости (расплава зерновой смеси) и изменения давления в цилиндрическом кольцевом зазоре матрицы экструдера. Приведены уравнения для двух случаев: с и
без проскальзывания расплава на стенках канала
Ключевые слова: математическая модель, экструзия, аномально вязкие жидкости
Экструзия является прогрессивным способом получения качественных продуктов питания, в том числе функционального назначения.
Особое место среди экструзионных продуктов
занимают коэкструдированные, состоящие из
зерновой оболочки и находящейся внутри
начинки [1, 2]. Используемые в качестве начинок жиросодержащие смеси относятся к аномально вязким (неньютоновским жидкостям).
Для их дозировки в матрицу для формирования
коэкструдированного продукта используется
насос высокого давления, для правильного
подбора которого необходимо знать характер
изменения давления и скорости расплава зерновой смеси (экструдата), формирующего зерновую оболочку при выходе из экструдера. Течение расплава экструдата в каналах с кольцевым зазором в матрице экструдера с зоной ввода начинки определяет не только качество готового продукта, но и производительность экструдера.
Целью работы явилась разработка математической модели течения расплава зерновой смеси
в кольцевом канале формующей матрицы экструдера с введением в центральную зону жиросодержащей начинки и получение информации
об изменении величины давления экструдата в
любой точке в кольцевого канала формующей
матрицы экструдера.
В основе математической модели процесса
течения псевдопластичной жидкости лежат законы сохранения массы, количества движения
и энергии [3]. Приняты следующие допущения:
течение является полностью установившимся
(ни в одной точке канала не наблюдается переходных режимов); силы инерции пренебрежи
Шевцов Александр Анатольевич – ВГУИТ, д-р техн.
наук, профессор, тел. (473) 255-65-11
Лыткина Лариса Игоревна – ВГУИТ, д-р техн. наук,
профессор, тел. (473) 255-65-11
Острикова Елена Александровна – ВГУИТ, аспирант,
e-mail: power-of-inspiration@yandex.ru
46
мо малы по сравнению с силами трения; течение ламинарное и изотермическое; является
полностью развитым; расплав несжимаем, его
плотность постоянна; течение происходит при
отсутствии внешних сил; влиянием силы тяжести можно пренебречь; расплав прилипает к
твердым стенкам канала (на стенках отсутствует проскальзывание).
Рассмотрим процесс течения псевдопластичной жидкости (расплава зерновой смеси) в
канале с кольцеобразным поперечным сечением (с внутренним радиусом R1 , внешним радиусом R2 и длиной L ) (рис. 1).
Рис. 1. Расчетная схема течения расплава и начинки
в канале матрицы
Предположим, что в нем имеет место равновесие сил, действующих на кольцевой массовый элемент с толщиной слоя dr , движущийся
со скоростью vz . В связи с тем, что равновесие
количества движения сводится к равновесию
действующих сил (это является следствием несжимаемости расплава и допущения, что расплав течет по прямолинейным параллельным
траекториям с постоянной скоростью), можно
записать:
p
dz;
z

(1)
(r  dr )  (r )  dr ,
r
где r, z – координаты, r – текущий радиус, м,
p – давление, Па.
Так как течение является полностью развившимся, то градиент давления можно считать
постоянным р / z  p /L . Отбросив все члены высшего порядка, получаем следующее
дифференциальное уравнение:
p  d  1 
(2)
 

( r ).
L r dr r r
В результате интегрирования уравнения (2)
получаем
С
р
(3)
 (r ) 
r 1 .
2L
r
Чтобы решить уравнение (3) для данного
случая, необходимо сделать предположение,
что напряжение сдвига  принимает нулевое
значение при r   r (рис. 1), где скорость vz
принимает максимальное значение (vz )max . Тар
ким образом, С1   ( R2 )2 . Подставляя C1 в
2L
уравнение (3), получаем
p R2  r
 R2 
 r  
(4)
 
.
2 L  R2
r 
Для течения ньютоновской жидкости скорость сдвига вычисляется по формуле [1]:
dv
p R2  r
 R2 
 r    z  
(5)
 
,
dr
2 L  R2
r 
где  – вязкость, Па∙с.
Проинтегрировав уравнение (5), получаем
2

p R2  r 
R 
   2  ln  2   C2  . (6)
vz ( r )  
2 L  R2 

 r 

Значения двух неизвестных величин  и
С2 можно определить, используя следующие
граничные условия: при r  k R1vz  0 ; при
r  R1vz  0 , здесь k  R1 /R 2 – соотношение
внутреннего и внешнего радиусов кольцевого
зазора.
Подставив граничные условия в уравнение
(6), получаем
1  ( R1 /R2 )2
2 2 
;
(7)
С2  1 .
ln( R1 /R2 )1
Распределение скоростей можно выразить
формулой:
p( z  dz )  p( z ) 
2
R22 p   r  1  ( R1 /R 2 ) 2  r  

vz ( r ) 
1   
ln    . (8)
4 L   R2  ln  R1 /R 2 1  R2  


Когда r  R , получаем выражение для
максимальной скорости потока (vz )max :
(vz )max 


 1  ( R /R )2   
R22 p 
1  ( R1 /R 2 )2 
1
2

 
1

ln
1 
.
 2ln  R1 /R 2 1   
4 L  2ln  R1 /R 2 1 

 


(9)
Средняя скорость получается из уравнения
(8):
R22 p 1  ( R1 /R 2 )4 1  ( R1 /R 2 ) 2 

 . (10)

4 L 1  ( R1 /R 2 ) 2 ln  R1 /R 2 1 


Объемный расход V получается умножением выражения для средней скорости на площадь поперечного сечения кольцевого зазора:
vz 
2
R 
V  R (1   1  ) vz 
 R2 
2 2

 R1  
. (11)

(1    )
4

R
 R1 
R24 
1
 2

(1    ) 
 p
1
8L 
R
 R1 
 2

ln




 R2 


Среднее время пребывания расплава в канале обратно пропорционально средней скорости
2
2
1
2
  R 4
 R1  
1
1    1    
R 
8 L2   R2 
t  L / vz  2 
  2   . (12)
2
R2 p
R 
1   R1 
ln  2  
  R2 
 R1  

Анализ формул для расчета объемного
расхода V показывает, что все они могут быть
записаны в обобщенной форме. Обобщенная
форма уравнения объемного расхода для ньютоновских жидкостей:
(13)
V  K p /  ,
где K – гидравлическая пропускная способность матрицы (зависит от геометрической
формы канала матрицы [3]).
Применив степенной закон Оствальда-де
Виля, получим уравнение для объемного расхода
(14)
V  K ' p m ,
где K ' 
 ( R2  R1 )( R2  R1 )m 2
2m  2 ( m  2 )Lm
– гидравличе-
ская проводимость матрицы (величина K ' является функцией геометрии канала матрицы и
индекса течения m);  – коэффициент текучести; m – индекс течения.
Однако формулы (13) и (14) применимы
только для расчета течения ньютоновских жидкостей. Расчетные формулы для псевдопластичных материалов при использовании степенного закона становятся громоздкими. Поэтому удобнее использовать метод характерной
47
вязкости [3], позволяющий применять уравнения для ньютоновских жидкостей для псевдопластичных материалов.
Согласно данному методу, если известна
некоторая характерная точка в канале, то для
этой точки можно вычислить эффективную
скорость сдвига  на основе объемного расхода V и характерной вязкости  для известной
кривой вязкости. Затем полученное значение
следует подставить в уравнения для ньютоновских жидкостей и определить перепад давления. Таким образом, можно получить уравнение течения в кольцевом канале матрицы, аналогичное уравнению (13), но использующее характерные параметры.
При выводе уравнений, описывающих течение псевдопластичной жидкости в кольцевом
канале матрицы экструдера, использовалось
предположение об отсутствии проскальзывания
на стенках канала, т. е. предполагалось, что
жидкость, текущая в канале матрицы, прилипает
к внутренней поверхности стенок, что равносильно граничному условию v = 0 на стенках
канала (рис. 2, а). Однако при превышении критического значения расплав может проскальзывать вдоль стенок с конечной скоростью v .
Критические значения напряжений сдвига для
расплава зерновой смеси составляют  крит 
0,12…0,14 Н/мм2, выше них наблюдается проскальзывание продукта вдоль стенок канала
матрицы (рис. 2, б).
Рис. 2. Профили скорости расплава зерновой смеси при
различных граничных условиях: а – прилипание продукта
к стенке; б – проскальзывание с конечной скоростью
Э. Уландом [4] предложена математическая модель течения, учитывающая явление
проскальзывания на стенках и представляющая
собой уравнение равновесия между силами
вязкости и силами трения для элементарного
объема жидкости при течении через трубу:
R dp
(15)
 ст  2
  p fc ,
ражение изменения давления по длине цилиндрического канала матрицы:
 2

(16)
p  pL exp  G ( L  z )  .
 R2

где G – вязкость расплава зерновой смеси, Па∙с.
Решая совместно (15) и (16), получаем
2 f

(17)
ст   p f c exp  c  L  z   .
 R2

Анализ уравнения (17) показывает, что при
скольжении расплава вдоль стенки напряжение
сдвига  ст на стенке изменяется по всей длине
канала (рис. 3) и для кольцевого канала матрицы экструдера напряжения сдвига равно
 cт 
R2  p  r
1  ( R1 / R2 )2 R2  ,



2 L  R2
2 ln( R2 / R1 ) r 
(18)
где R1 – внутренний радиус кольцевого канала
матрицы, м; L – длина цилиндрического кольцевого канала, м; r – текущее значение радиуса кольцевого канала матрицы, м;  p – перепад давления в канале матрицы, Па.
Из уравнения (17) следует, что сила трения
возрастает с удалением от выхода из канала
матрицы. Вероятно, что сила трения при этом
возрастет настолько, что проскальзывание станет невозможным ( z< z1 ). В этом случае
напряжение сдвига на стенке  ст будет меньше,
чем напряжение сдвига, необходимое для преодоления сил трения. Поэтому на начальном
участке течения наблюдается прилипание расплава к стенкам канала ( 0<z<z1 ). Дальше, на
участке z1 <z<L , начинает развиваться проскальзывание на стенках, сопровождаемое
сдвиговым течением жидкости (рис. 3).
Из совместного решения уравнений (17) и
(18) получаем
R2  1  1  ( R1 / R2 )2 R22  
z1  L 
ln 
r 
  .(19)
2 f c  2 L f c  2ln ( R2 / R1 ) r  
2 dz
где R2 – внешний радиус кольцевого канала
матрицы, м; f c – коэффициент трения при проскальзывании; р – давление в рассматриваемом
элементарном объеме, Па.
Интегрирование уравнения (15) с учетом
граничных условий p  pL при z  L (где pL –
давление в конце канала матрицы, Па) дает вы-
48
Рис. 3. Изменение давления и напряжения сдвига на стенке канала по его длине при проскальзывании расплава на стенках
Скольжение расплава по стенкам на протяжении всей длины канала будет всегда
наблюдаться при z1 = 0. Необходимый для этого объемный расход Q может быть получен из
формулы
m
 2 fc L 
 R23 
Q
 pL f c exp 
 .
m  3 
 R2  
(20)
Давление в точке отрыва расплава от стенки z1 вычисляется из уравнений (19) и (16):
m
  R23 
(21)

 .
  m  3Q  
Давление на входе в фильеру po определяется из уравнения
p0  p1
dp
(22)
 .
z1
dz
dp
p
С учетом того, что
, а также на

dz
L
основании уравнений (19) и (20), получаем
1
p1 
fс
m
m 


  R23   2 L 1 
1   R23    
p0  


   1  ln 
.

 PL f c   m  3 Q    
  m  3 Q   R2 f c 

  



(23)
При полном прилипании расплава к стенкам по всей длине канала давление p0H вычисляется по следующей формуле:
m
p0 H
2 L   R23 


 .
R2   m  3 Q 
(24)
Это значение больше давления, наблюдающегося в случае, когда на некотором участке
канала наблюдается проскальзывание расплава
на стенках (см. рис. 3). В. Микаэли [5] предложил следующую формулу для определения поp
терь давления
при течении псевдоплаL
стичной жидкости в кольцевом канале матрицы
Дальнейшее преобразование (27) позволяет
получить выражение для vcт
vcт 
2 f m

R
Q
   PL f c m 2 exp  c  L  z  . (28)
2
R2
m3
 R2

Из уравнения (28) видно, что максимальное значение vcт наблюдается в конце канала
матрицы, а в точке отрыва потока z1 (переход
от прилипания к проскальзыванию) vcт = 0 (см.
уравнение (19).
На рис. 4 показано развитие профиля скорости на участке канала матрицы экструдера
z1  z  L. Для рассматриваемого здесь примера на участке ( r / R  0,425) наблюдается
чисто сдвиговое течение, для которого справедливо граничное условие vz  vст  0 на
стенке канала, тогда как осевое течение устанавливается на выходе из канала (где vz  vст ).
Очевидно, что существует компонента скорости в радиальном направлении r , которой
можно пренебречь, так как vz vст [4, 5]. Таким образом, установлено, что явление проскальзывания расплава по стенке наиболее выражено к выходу из канала, а в противоположном направлении его влияние уменьшается.
Это означает, что проскальзывание на стенках
начинается в зоне выхода из канала. Однако
существенным препятствием применения представленного метода является невозможность
точного определения коэффициента трения f c
для расплавов зерновых смесей и проблема
корректного выбора скорости проскальзывания
vcт на внутренней стенке канала матрицы.
m
p
  DH m  2 
  m 1
(25)
 .
L  2 ( m  2 )Q 
На основании уравнения (17) с учетом степенного закона Оствальда-де Виля    m , а
также условия   cт r / R получаем следующее уравнение:
m

 f

dv
r 
   z    pL f c  exp  c  L  z  . (26)
dr
R2 

 R2

Интегрируя уравнение (26) с учетом граничного условия vz  vст при r  R2 , получим
уравнение для профиля скоростей:
m
 p f  R m 1  r m 1
2 f m

vz  vcт    L c 
exp  c  L  z  . (27)
m 1
 R2 
 R2

Рис. 4. Развитие профиля скорости в канале матрицы
для z1  z  L : R = 0,040 м; Q = 12,010-3 м3/с; p = 6,5
МПа; f c = 0,2; m = 2,3
Следовательно, при разработке конструкции матрицы экструдера для получения коэкструдированных продуктов, имеющих тенденцию к проскальзыванию на стенках канала мат-
49
рицы, полученные решения (24) и (28) являются достоверными, т. к. получены на основе учета явления проскальзывания расплава на стенках канала. Они дают более точные значения
перепада давления и скорости проскальзывания
vcт по сравнению с наблюдающимися в используемых экструдерах.
Был проведен машинный эксперимент с использованием системы символьной математики
Maple для среды плотностью  = 1230 кг/м3 и
динамической вязкостью  = 13680 Пас и значением индекса течения m  1, 3 . Объемный
расход жидкости составлял V = 0,421∙10-5 м3/с.
Геометрия области течения – цилиндрический
кольцевой канал длиной L = 25∙10-3 м.
Расчет осевой скорости проводился для
N = 5 участков. Потери давления являются
важной характеристикой течения жидкости в
цилиндрическом кольцевом канале. На рис. 5
приведен график изменения потерь давления
вдоль оси канала.
в цилиндрическом кольцевом канале матрицы по оси z
Проведен сравнительный анализ характера
течения расплава зерновой смеси при различных значениях индекса течения m  1, 3 . Проверка адекватности представленной модели показала, что результаты вычислительного эксперимента согласуются с экспериментальными
данными, погрешность изменялась в пределах
18 %, что говорит о возможности использования полученных результатов в проектировании
экструдеров.
Литература
1. Остриков А. Н., Василенко В. Н., Соколов И. Ю.
Коэкструзионные продукты: новые подходы и перспективы. М.: ДеЛи принт, 2009. 232 с.
2. Лыткина Л. И., Шевцов А. А., Дранников А. В.,
Клейменов А. И. Техника и технология тепловых и механических процессов в задачах энергосбережения на комбикормовых заводах. Воронеж: ВГТА, 2011. 304 с.
3. Остриков А. Н., Абрамов О. В., Василенко В. Н.,
Попов А. С. Математическое моделирование течения
аномально вязких сред в каналах экструдеров. Воронеж:
Изд-во ВГУ, 2010. 237 с.
4. Uhland E. Modell zur Beschreibung des Fliepens
Wandgleitender Substanzen durch Disen, Itlieol. Acta, 15.
1976. Pp. 30-29.
5. Микаэли В. Экструзионные головки для пластмасс и резины. Конструкция и технические расчеты. Пер.
с англ.: Под ред. В. П. Володина. СПб: Профессия, 2007.
472 с.
Рис. 5. Изменение давления  Р при течении экструдата
Воронежский государственный университет инженерных технологий
MATHEMATICAL MODEL OF QUASI-VISCOUS LIQUIDS FLOW IN THE PROCESS
OF COEXTRUDED COMPOSITIONS OBTAINING
A.A. Shevtsov, L.I. Lytkina, E.A. Ostrikova
Mathematical model of quasi-viscous liquid (molten grain mixture) flow and pressure variation in cylindrical annular
space of extruder die is developed. Equations for cases either with or without melt sliding on channel walls are discussed
Key words: mathematical model, extrusion, quasi-viscous liquid
50
Документ
Категория
Без категории
Просмотров
13
Размер файла
464 Кб
Теги
аномально, жидкостей, вязких, математические, pdf, коэкструдированных, получения, модель, течение, композиций
1/--страниц
Пожаловаться на содержимое документа