Математическая модель неизотермического течения высоковязких сред в каналах матрицы экструдера Текст научной статьи по специальности «Физика»

УДК 532

Старший преподаватель А.С. Сидоренко,

(Военный учебно-научный центр военно-воздушных сил «Военная воздушная академия им. проф. Н.Е. Жуковского и Ю.А. Гагарина») кафедра общепрофессиональных дисциплин. тел. 8(904)210-17-90 E-mail: sas1.vrn@mail.ru

доцент А.И. Потапов

(Воронеж. гос. ун-т. инж. технол.) кафедра машин и аппаратов пищевых производств

тел. 8(906)586-75-97

E-mail: a.i.potapov@rambler.ru

Senior lecturer A.S. Sidorenko,

(Russian Air Force Military Educational and Scientific Center "Air Force Academy named after Professor N.E. Zhukovsky and Y.A. Gagarin") Departmentof All-Professional Disciplines. phone 8(904)210-17-90 E-mail: sas1.vrn@mail.ru

associate professor A.I. Potapov

(Voronezh state university of engineering technologies) Department of food production machines. phone 8(906)586-75-97 E-mail: a.i.potapov@rambler.ru

Математическая модель неизотермического течения высоковязких сред в каналах матрицы экструдера

Mathematical model non-isothermal flow highly viscous media channels matrix extruder

Реферат. Рассматривается одномерное установившееся течение высоковязкой среды в цилиндрическом канале с учетом диссипации и зависимости коэффициента вязкости от температуры. Предположено, что на сравнительно малых интервалах температуры изменение коэффициента динамической вязкости с достаточной степенью точности можно принять линейным. В основу модели были положены уравнения гидродинамики жидкости и теплопереноса. В поставленной задаче температура стенки канала принимается постоянной. Получено приближенное решения рассматриваемой задачи, в соответствии с которым распределение скорости, давления и температуры ищется в виде разложения по степеням безразмерной поперечной координаты. Рассмотрен частный случай, когда в соотношениях распределение скорости, давления и температуры допустимо ограничиться следующим числом членов разложения: для скорости - первые 3, для давления - первые 2, для температуры - первые 5. Получены выражения для определения профиля температуры среды в канале и определения характеристик диссипативного разогрева. Для моделирования процесса теплопере-носа высоковязких сред разработана программа для персональных электронно-вычислительных машин. Расчет проведен по экспериментальным данным исследования течения расплава зерновой смеси гречихи и сои для скорости нагрузки 0,08 мм/с. Методом машинного эксперимента осуществлена проверка полученных решений на адекватность реальному процессу теплопереноса. Анализ результатов указывает, что при малых значениях длины канала влияние диссипативной функции проявляется главным образом у стенки. При увеличении приведенной длины это явление распространяется на все сечение канала. При большой длине канала профиль температур целиком определяется диссипацией. В случае теплообмена, обусловленного только теплотой трения, форма кривых распределения температур является следствием взаимодействия эффектов нагрева за счет вязкого сдвига с эффектами охлаждения за счет теплопроводности. Отклонение расчетных данных от экспериментальных по абсолютному значению не превышало 12 %.

Summary. We consider a one-dimensional steady flow of highly viscous medium in a cylindrical channel with Dissipation and dependence of the viscosity on the temperature. It is assumed that a relatively small intervals of temperature variation of the dynamic viscosity with a sufficient degree of accuracy can be assumed to be linear. The model was based on the equations of hydrodynamics and the heat transfer fluid. In the task channel wall temperature is assumed constant. An approximate solution of the problem, according to which the distribution of velocity, pressure and temperature is sought in the form of an expansion in powers of the dimensionless transverse coordinates. A special case, when the ratio of the velocity distribution, pressure and temperature is allowed to restrict the number of terms in the expansion as follows: for speed - the first 3 to the pressure - the first two for the temperature - the first 5. The expressions to determine the temperature profile of the medium in the channel and characterization dissipative heating. To simulate the process of heat transfer highly viscous media developed a program for personal electronic computers. The calculation was performed using experimental research data melt flow grain mixture of buckwheat and soybeans for the load speed of 0.08 mm / s. The method of computer simulation carried out checks on the adequacy of the solutions to the real process of heat transfer. Analysis of the results indicates that for small values of the length of the channel influence dissipation function appears mainly at the wall. By increasing the reduced length of this phenomenon applies to all section of the channel. At high temperature profile along the channel length is determined entirely by dissipation. In the case of heat transfer due to fric-tional heat only, the form of curves of temperature distribution is a consequence of the interaction effects of heating due to viscous shear effects cooling by conduction. The deviation from the experimental data calculated by the absolute value does not exceed 12%.

Ключевые слова: математическое моделирование, неизотермическое течение, конвективный теплоперенос

Keywords: mathematical modeling, isothermal flow, convective heat transfer

© Сидоренко А.С., Потапов А.И., 2015

Несмотря на то, что теоретическому исследованию неизотермического течения высоковязких сред с учетом диссипации энергии посвящен ряд работ, однако до настоящего времени эти вопросы изучены явно недостаточно. Осложняющими факторами моделирования подобных процессов является зависимость физико-механических свойств расплавов экструдата от температуры. В обычных условиях наиболее сильно изменяется с температурой коэффициент динамической вязкости. Поэтому ограничиваемся учетом изменения одной лишь вязкости, полагая остальные физические свойства постоянными [1-4].

Рассмотрим одномерное плоское установившееся течение высоковязкой жидкости в цилиндрическом канале. Канал имеет длину L и радиус R .

Введем цилиндрическую систему координат. Начало отсчета продольной оси z системы совместим с входным сечением канала.

Рассмотрим движение жидкости за счет разности давлений на входе Р0 и на выходе из канала РL.

В поставленной задаче температура стенки канала принимается постоянной и равной Т№ .

Одномерное течение и теплоперенос высоковязких сред с учетом общепринятых допущений описывается системой уравнений неразрывности, гидродинамики и конвективного теплопереноса с учетом диссипации механической энергии, которые могут быть записаны в цилиндрической системе координат следующим образом:

ди _ п д/(Т) ди дР _

-— 0;----— 0 ;

дг дг дг дг

1 Аи у и 1-дР — о;

г дг ^ дг ^ дг

ттдТ 1 д ( дТ^ (ди рсИ — — Х--[г —1 + /(Т)[ ——

дz г дг ^ дг) ^ дг Запишем граничные условия:

(1)

— 0;

аи

"дГ

И о — 0;

1г—R

г—0

дт

дг

т|

— 0;

г—0

— Т •

г—^ " ^ ;

т(г)|2—0 — X твх^; 2—0 j—0

(2)

где р - плотность среды; г, 2 - цилиндрические координаты; Р-давление среды в канале; И-скорости частиц среды; с - удельная теплоемкость среды в канале; X - коэффициент теплопроводности среды; /(Т) -коэффициент динамической вязкости, зависящий от температуры; Т • -заданные коэффициенты в пред-вхЛ

ставлении профиля температуры во входном сечении канала, определяемые из дополнительных условий(принятый здесь вид граничного условия для температуры при 2 равном нулю, как это будет следовать из дальнейшего рассмотрения обуславливается видом предполагаемого приближенного решения для температуры в канале); J - степень полинома, используемого для аппроксимации начального профиля температуры на входе в канал.

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

/(т) —/0 +/1 (тхар - т^ (3)

где /0, /1 - эмпирические коэффициенты, определяемые экспериментально для каждой жидкости в отдельности; Тхар- некоторая характерная температура процесса.

Для перехода к безразмерному виду введем следующие переменные, параметры, функции и критерии:

г' — В — ^ ; 2' — Ей — ^ 15 £ Ь PV12,

' — Твх^ — Твх,0 , — Т — Твх,0 , — ; Т —

, j т Т ^ т т

Ткрит — твх,0 Ткрит — твх,0

_ тхар твх,0 _ V» ОЛ.и / л\

Тхар— т —^г;Ту— т ; (4) , Р — PL , и - 0

Р' —-— ; и' — —; 0 —-^-;

Р0 — PL ^ пяЧ

/1 (Ткрит — Твх,0 ) / С

Vls ———-—; Рг ——

/ 0 Х

Tw Твх,0

—Ь — РЬ,

1

2

п 2pVmR ^

Re = m ; Ec =

V,

*Т

m

крит

- Т

вх,0

)"

где г' - текущий радиус канала; D - вспомогательный геометрический параметр, характеризующий отношение гидравлического диаметра канала к его длине; z' - безразмерная

продольная координата; Т,Тхар,Ту -значения

I

температуры; Ткрит - безразмерное значение критической температуры; P' -давление среды в канале; U' -скорость среды; Q - расход среды через канал; Vis - параметр характеризующий отношение характерной вязкости к величине масштабной вязкости; Pr, Re, Ec, Eu -критерии Прандтля, Рейнольдса, Эккерта, Эйлера, соответственно.

Здесь и далее верхние штрихи указывают на безразмерную форму записи соответствующих величин.

I

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

I

в качестве Ткрит может быть принята, например, температура фазового перехода.

С учетом (4) основные уравнения (1) могут быть записаны в следующей форме

dU dz'

= 0;

1 + V,s(4p -Т'))= 0,

dP D dU d h ... „ .. _.

+ EURe77+ МТх--Т " =0;

дг EURe дг dz

Eu Re D dP

4

1 д r dr

dz

(l + Vis(fxap - T )) r

dü' dr

(5)

Pr Re и dT

D dz'

4 i _d_

D2 r dr

dT_ dr'

+

+

^[l + -V

D2

dr

С учетом (4) граничные условия (2) принимают вид:

dU

dr U

= 0;

r'= 0 = 0;

dr

dr Т

= 0;

г = 0

г = 0

= Т

г = 0

w '

Т

P

'(г')

J ' 'j

z=0 j=0 ,j

(6)

= 0.

z =1

В общем случае точное решение системы уравнений (5) с граничными условиями (6) представляется затруднительным. В этой связи рассмотрим подход к построению приближенного решения рассматриваемой задачи, в соответствии с которым распределение скорости, давления и температуры будем искать в виде разложения по степеням безразмерной поперечной координаты г':

'(г')=|. Uj(r j

(7)

P

Т

(8)

(9)

и (г 1= Е и

j=0

>'(гУ )=Е р (г')(г'} ;

j=0

'("У)= 2 ^(г')(г' j=0

Подставляя (7) - (9) в (5), приходим после соответствующих преобразований к системе обыкновенных дифференциальных уравнений

относительно коэффициентов разложения (г ), Т (г ), как, функций продольной координаты г ,

и постоянных коэффициентов j = 0,1,....

Рассмотрим наиболее простой частный случай, когда в соотношениях (7) - (9) допустимо ограничиться следующим числом членов разложения:

- для скорости - первыми тремя;

- для давления - первыми двумя;

- для температуры - первыми пятью.

В рамках рассматриваемого приближения с учетом допущения о параболическом профиле скорости коэффициенты разложения в (7) определяются из выражений:

' О ' ' '

=у; и = 0; и2 =-Ц). (10)

2

r

В итоге после соответствующих преобразований приходим к следующей системе обыкновенных дифференциальных уравнений относительно Т0, Т', Т2, Т3, Т_4, Р0 и Р[ как функций продольной координаты:

1 + ^(4р — Т0 ))и

— Т0))и2 — -Еи5е° др0-

16

дг

РгЯеБИ0 дТ0 — 16Т2;

(11)

(

Р^еБ

И,

дг

дТ2

с2

^ + И2-

дТ0 !

дг

— 64

Т4 + ЕсРг1 + VlsT

(1 + У^р )(и2 )2

Р1— 0; Т1— 0; Т0 + Т2 + Т4 — Т^ .

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

После некоторых преобразований система уравнений (11) сводится к одному дифференциальному уравнению относительно распределения давления в канале:

■5 ! Л ! !

+ + В2^ — в

4')3 1 ^

3;

(12)

здесь для краткости записи приняты обозначения:

в Рг2Яе3Б3ЕиИ' ;В

в0 —-о-и0 ; В1 —

162Vls

(

В — 64Ем 5е В

3РгЕиЯе2Б2 16Vls

Л

V У18и0

+ Ее Рг ип

)

Вз — 64[ ТЛ„ — Т

1хар

1

Принимая во внимание функциональные связи между коэффициентами разложения для давления и температуры, вытекающие из (11), граничные условия для уравнения (12) с учетом (6) записываются в виде уравнения (13).

После нахождения распределения давления вдоль оси продольного канала, температурное поле определяли из (9) с учетом уравнений (11).

дР0

дг

г — 0

16И2 ЕиЯеБ

(1 + У^(гхяр — Т,

хар Твх,0

а2Р0

4?

—0;Р0

(13)

— 0

г—1

г —0

Получено аналитическое решение уравнения (12):

Р0

С1 cos

( А^

2

V )

+ С2 sln

(Аг^

2

V )

(14)

2В3 + С3 +

В0 В1

где для краткости записи приняты обозначения:

В2 В3

С1

( , Л

дР0 В0

дг г'— 0 В1

V )

А2 +(В2/В3 )

,2

С 2 — 1

2 А

(

дР0

дг

г — 0

—В0+^С!

В1 2В3 1

Л

В2

С3 — —е

2В3

С1с°Н у) + С2^п(АА

В0

В1

А — Л4В!/В3 — (В2/В3)2 .

Средняя температура в канале описана уравнением:

Т'ат^'):

1

| г'Т'(г',г')и'(г')аг'

_0_

1

| г'и'(г')аг'

(15)

Для моделирования процесса теплопере-носа высоковязких сред разработана программа на ПЭВМ. Расчет проведен по экспериментальным данным исследования течения расплава зерновой смеси гречихи и сои для скорости нагрузки 0,08 мм/с [5].

Исходные данные и результаты расчета приведены в таблице. На рисунке 1 приведено распределение температурного поля расплава смеси в зависимости от длины канала.

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

х

В

2

г

х е

4

32,0 z,mm 16,0

г 458 453 448 443 1-438 433

T,K

0,00 0,33 0,66 0,99 г,мм

Рисунок 1. Распределение температурного поля расплава смеси зерна гречихи и сои

Методом машинного эксперимента осуществлена проверка полученных решений на адекватность реальному процессу теплопере-носа. Отклонение расчетных данных от экспериментальных по абсолютному значению не превышало 12 % (таблица 1).

Т а б л и ц а 1

Исходные данные и результаты

расчета

Параметр

Радиус канала, м Длина канала, м

Температура на входе в канал, К Температура стенки канала, К Характерная температура, К Критическая температура, К Скорость потока в канале, м/с Избыточное давление на выходе, кПа

Эмпирический коэффициент Ц0, кПас

Эмпирический коэффициент Ц1, кПас

Теплопроводность расплава, Вт/(мК)

Удельная теплоемкость расплава, Дж/(кгК)

Плотность расплава, кг/м3 Расход среды через канал, кг/м3 Число Vis Число Прандтля Число Рейнольдса Число Эккерта

Число Эйлера_

Значение

0,0011 0,0320 433 453 443 459 0,8 10-4

31.00

11.1

0,240

0,220

1600 1200 1,15710-8

0,562 8,073107 1,18910-5 6,0110-8 1,167104

Выводы:

1. Выявлен характер изменения поля температур высоковязких сред в различных сечениях цилиндрического канала.

2. Полученная модель позволяет с достаточной точностью (±12 %) рассчитать не только значения температур высоковязких сред в цилиндрическом канале, но и провести анализ влияния основных параметров системы на характеристики диссипативного разогрева материала.

ЛИТЕРАТУРА

1 Василенко В.Н., Остриков А.Н., Ряжских В.И. Математическая модель течения двух вязкопла-стичных сред в формующем канале экструдера при коэкструзии // Вестник Воронежского государственного университета инженерных технологий. 2012. № 2. С. 64-67.

2 Остриков А.Н., Татаренков Е.А., Попов А.С., Василенко В.Н. Моделирование течения расплава биополимера в динамической матрице экструдера // Хранение и переработка сельхозсырья. 2011. № 8. С. 66-69.

3 Остриков А.Н. Математическое моделирование течения аномально-вязких сред в каналах экс-трудеров. Воронеж : Изд-во ВГУ, 2010. 237 с.

4 Комов А.А., Потапов А.И., Тарарыкова И.В, Шахов С.В. Математическое описание процесса микрофильтрации суспензии в трубчатом канале // Современные наукоемкие технологии. 2014. № 5-1. С. 164-165

5 Кретов И.Т., Попов Е.С., Потапов А.И., Попов Д.С Математическое моделирование процесса микрофильтрации // В книге: Материалы LI отчетной научной конференции преподавателей и научных сотрудников ВГУИТ за 2012 г. 2012. С. 42.

REFERENCES

1 Vasilenko V.N., Ostrikov A.N., Rya-zhskikh V.I. A mathematical model of the flow of two viscoplastic media in the form of a channel-le extruder at coextrusion. Vestnik VGUIT. [Proceedings of VSUET], 2012, no. 2, pp 64-67. (In Russ.).

2 Ostrikov A.N., Tatarnikov E.A., Popov A.S., Vasilenko V.N. Modeling of biopolymer melt extruder in the dynamic matrix. Khranenie i pererabotka sel 'khozsyr 'ya. [Storage and processing of agricultural-hozsyrya], 2011, no. 8, pp. 66-69 (In Russ.).

3 Ostrikov A.N. Matematicheskoe modeliro-vanie techeniya anomal'no-vyazkikh sred [Mathematical modeling of abnormally viscous media channels extruders]. Voronezh, VGU, 2010. 237 p. (In Russ.).

4 Komov A.A., Potapov A.I., Tararykova I.V., Shakhov S.V. The mathematical description of the process of suspension in a tubular micro channel. Sovremennye naukoemkie tekhnologii. [Modern high technologies], 2014, no. 5-1, pp 164-165(In Russ.).

5 Kretov I.T., Popov E.S, Potapov A.I., Popov D.S. Mathematical modeling of microfiltration. Materialy LI otchetnoi nauchnoi konferentsii prepo-davatelei VGUIT. [In the book: Materials LI reporting conference of teachers and researchers VSUET for 2012]. 2012, pp. 42 (In Russ.).