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

BernnuxJBryWT/Proceedmgs of VSUET, Ж 80, № 2, 2018-

Оригинальная статья/Original article_

УДК 664

DOI: http://doi.org/10.20914/2310-1202-2018-2-18-22

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

механических колебаний

Юлиан И. Шишацкий 1 Алексей М. Барбашин 1

Сергей А. Никель 1 sergei.nickel@yandex.ru

1 Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия Реферат. В нашем случае твердым телом является сырье растительного происхождения - люпин, измельченный в крупку, а экстрагентом -подсырная сыворотка. Турбулентная обстановка в аппарате создавалась наложением низкочастотных механических колебаний, которые оказывают значительное влияние на характеристики гидромеханических, массообменных и тепловых процессов. Эту особенность необходимо учитывать при расчетах экстракционных аппаратов. Сформулированы основные допущения для решения поставленной задачи. Записано уравнение движения одиночной частицы, которое приводится в ряде работ (Соу, Хинце, Чен, Протодьяконов и др.). Оно справедливо при мгновенных значениях параметров. Выписано более простое уравнение, описывающее движение дисперсной частицы, а также тензоры временной корреляции с последующим их разложением в интеграл Фурье. Далее, учитывая определение тензоров, показаны зависимости для расчета интенсивности хаотического движения сплошной и диспергированной фаз, а также получено конечное выражение, показывающее соотношение интенсивностей движения фаз. Коэффициент турбулентной диффузии каждой из фаз пропорционален интенсивности хаотического движения соответствующей фазы. Поэтому записанное конечное уравнение для соотношения фаз позволяет оценить соотношение коэффициентов турбулентной диффузии жидкой и диспергированной фаз в экстракционном аппарате. В нашем случае отношение плотностей рт/рж составляет 1,1. Поскольку плотности люпина и подсырной сыворотки количественно разнятся, то следует ожидать некоторого увеличения относительной скорости движения фаз, что приведет к увеличению скорости массоотдачи. Интенсивности хаотического движения фаз не будут одинаковыми, равно как и коэффициенты турбулентной диффузии. Таким образом, рассмотренный случай движения одиночной частицы в турбулентном потоке сложен и может быть решен только при достаточно серьезных допущениях, сформулированных ниже.

Ключевые слова: механические колебания, турбулентное перемешивание, сферическая частица люпина, подсырная сыворотка, математическая модель

Mathematical model of movement of a single spherical lupine particle in _the extractor using low-frequency mechanical vibrations_

Yulian I. Shitshatskii 1 Aleksei M. Barbashin 1

Sergei A. Nikel_1 sergei.nickel@yandex.ru

1 Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia Summary.In our case, the solid body is the raw material of plant origin-lupine, crushed into grits, and the extractant is the cheese whey. The turbulent situation in the apparatus was created by the imposition of low-frequency mechanical vibrations, which have a significant impact on the characteristics of hydro-mechanical, mass transfer and thermal processes. This feature must be taken into account in the calculation of the extraction apparatus. The basic assumptions for the solution of the problem are formulated. The equation of motion of a single particle, which is contained in a number of works (Sow, an introduction, Chen, Protodyakonov, etc.). It is true in the instant values of the parameters. A simpler equation describing the motion of the dispersed particle and time correlation tensors with their subsequent decomposition into the Fourier integral are written. Further, taking into account the definition of tensors, the dependences for the calculation of the intensity of the chaotic motion of continuous and dispersed phases are shown, and the final expression is obtained, showing the ratio of the intensities of the phases. The coefficient of turbulent diffusion of each phase is proportional to the intensity of the chaotic motion of the corresponding phase. Therefore, the written finite equation for the phase ratio allows to estimate the ratio of the turbulent diffusion coefficients of the liquid and dispersed phases in the extraction apparatus. In our case, the ratio of the density of Hg / Hg is 1.1. Since the density of lupine and cheese whey differ quantitatively, we should expect some increase in the relative velocity of the phases, which will increase the rate of mass transfer. The intensities of the phases chaotic motion will not be the same, as well as the coefficients of turbulent diffusion. Thus, the case of motion of a single particle in a turbulent flow is complex and can be solved only under sufficiently serious assumptions formulated below. Keywords :mechanical vibrations, turbulent mixing, spherical lupine particle, subsurface serum, mathematical model

Турбулентное перемешивание двухфазной системы осуществляется в различных аппаратах, например, в экстракторах. В нашем случае твердым телом являлось сырьё растительного происхождения люпин, измельчённый в крупку, а экстрагентом - подсырная сыворотка. Молочно-растительный экстракт люпина (последний

Для цитирования Шишацкий Ю.И., Барбашин А.М., Никель С.А. Математическая модель движения одиночной сферической частицы люпина в экстракторе с помощью низкочастотных механических колебаний // Вестник ВГУИТ. 2018. Т. 80. № 2. С. 18-22. doi: 10.20914/2310-12022018-2-18-22

характеризуется высокой массовой долей белков - более 30 % на сухое вещество) может быть использован в различных отраслях промышленности, в частности, как сырье для функциональных продуктов питания. Турбулентная обстановка в экстракторе обеспечивалась наложением низкочастотных механических колебаний.

For citation

Shishatsii Yu.I., Barbashin A.M., Nickel S.A. Mathematical model of movement of a single spherical lupine particle in the extractor using low-frequency mechanical vibrations. Vestnik VGUIT [Proceedings of VSUET]. 2018. vol. 80. no. 2. pp. 18-22. (in Russian). doi:10.20914/2310-1202-2018-2-18-22

<Ветшк&ТУИТ/Фгоиг£^ о/ЧЯиЕТ, Т. 80, № 2, 2018-

Последние оказывают значительное влияние на характеристики гидромеханических, массооб-менных и иных процессов, что следует учитывать при расчётах экстракционных аппаратов.

Принятые допущения:

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

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

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

Уравнение движения одиночной частицы приводят С. Соу, И.О. Хинце, С.М. Чен. И др. И.О. Протодьяконов записал его в виде [3]

1 И3

— па рр-6

аю„

1 6

. 3 др 1 3

— па--1--па р

1

- 3пу<ар (ю -Юр) -а (ю - Юр)

дг 12

(1)

+—

а

2 '' {

а/

'(ю-юр)

где Юр , рр - 7-я составляющая скорости частицы и её плотность; о)7, р - 7-я составляющая скорости жидкости и давление в точке пространства, в которой находится частица в момент времени 1 V - вязкость жидкости.

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

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

для ускорения так называемой присоединённой массы частицы, обусловленной тем, что движущаяся частица вовлекает в движение окружающую её сыворотку. Четвёртый член учитывает влияние нестационарности течения на движение крупки (сила Бассэ). В целом правая часть характеризует полную силу, действующую на частицу в потоке.

Поскольку считается, что частица не влияет на движение, из (1) градиент давления можно исключить. Используя уравнение Навье - Стокса, описывающее движение жидкой фазы [2]:

(и V) и = -—Ур + мДи , Р

где и - вектор скорости; Ди - оператор Лапласа; V - кинетическая вязкость жидкости, V = ц / р, из уравнения (1) получается следующее равенство:

—лс?р„= —пvdр(ю -Юр) + 6 Ир Л 4 7 р'

1 ,3 ( дю дю д ю

+—пО р| —L + юк—7- -V-—

6 1 д/ дх дх, дх.

1 ,3 0 (ю-Ю?7), 3. /- 'гО^-^

(2)

+—пО р 12

а/

- +—а гр4яу^~

]/ - /'

О/',

В уравнении (2) правая часть в совокупности представляет собой силу, действующую на частицу. Её правомерно считать случайной силой. Это справедливо, учитывая сформированные выше допущения применительно к турбулентной обстановке в экстракторе.

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

у-

1 дюг,

2 д/

9

+—

д

3 дю

2 д/

&

О/

■—Г (ю-ЮР

'(ю-юр)

(3)

п / лЦ-Р '

í0

где у = рр / р.

Статистическая связь соответственно между значениями скорости в моменты времени 11 и 1 для дисперсной частицы, дисперсной частицы и сплошной фазы, а также сплошной фазы характеризуются соответственно тензорами временной корреляции:

(1/) = Юр, (/1)Юр, (/) р7 (/1,/) = Юр, (/1)ю (/)

(/l, /) = Юр1 (/1 )Ю(/1).

<Ветшк&ТУИТ/Фгоиг£^ о/Т. 80, № 2, 2018В этих формулах черта означает операцию статистического осреднения.

Умножая почленно уравнение (3) на а. (^ )

и на а.. (т) , а также осредняя, получим:

г+1 На (* )=3 ср (* )-

7 2 I 1 У ' 2 С* 'Л '

18у

[Р (*)- 0. (*)]■

(4)

+ 9 уф [ Р ( ^^ ( "Я ( С\ у/7-7

г+1 \-0& (* )=3(* )-7 2 I Су3 У ' 2 С* ,Л '

18^

[*з (*)- Р3 (*)] +

(5)

+ , Ы с[* ^("Я й*,.

с V п 8 77-7

Здесь 8 = 1 - 11, 8' = 1' - 11, 8 и 8' - переменные. Вместо тензоров Рщ, R1j, Q1j подставляем в уравнения (4) и (5) их разложения в интеграл Фурье. Тогда получим:

1= }¿а- (-"/у + )1 Оз (-) =

Ы}¿«е-1«*{3(-")Р& (-) + ^[&(-)-&")]} +

сСа

(6)

(7)

-75= -

+ЩтС-7Ш(•)-а (а)]; ^2=1(-<«)(,+ 21& —) =

чЫ«-" {I (-«) а(а)- 1Т [& (")-& -)]}+

+Щ ТГ-71 Ж & (-)-Р")]

Меняя порядок интегрирования и вводя новую переменную 1 = 8 - 8', последний член уравнения (6) запишем иначе:

9Ш1Ые- (-")[& ("-&(а)]} Т-

Вычислим последний интеграл в этом уравнении, используя известные значения определённых интегралов [1]:

С С

Г—со8—т= [—81П — Т = -у/п / 2а.

0 V т 0 V т

Тогда этот член принимает следующий вид:

9 (V г Са

с V 2 ¿>/2=

е--(1 + 1)(-1л-)[&. (а)-Оз (а)].

Используя изложенную процедуру вычислений, последний член уравнения (7) примет вид:

9 V г Сае

-(1-1)Н"[&(а}-&(О} (8)

Л 2 -г 42=

Учитывая уравнения (8) и (9), можно заключить, что каждое из уравнений (6), (7) представляет собой равенство двух интегралов Фурье, откуда следует равенство преобразований Фурье. Поэтому получаем:

Г + 2 ] & (-) = 2 а (а) + 1 ^ а (а)-&(а)] + ^

+ С^1 + ()[& (а) -&(а)];

Г + ))& (-) = )&(а)-((а)-&")] +

(10)

"'Л* (а)-&(а)].

Теперь из уравнения (9) и (10) исключаем тензор Р. (а) . После этого получим соотношение,

которое связывает преобразования Фурье тензоров временной корреляции сплошной и диспергированной фаз:

& (-) = Г(а)О. (а)

Г (-) =

1 9

■ — + —

2 С \ 2а

9 V. 18у С V 2а С2а

(11) (12)

3 9 V 2 9 V + 18у

2 С V 2- + С\ 2а С2а

И.О. Хинце [4] установил, что интенсивности хаотического движения сплошной и диспергированной фаз, соответственно, определяются соотношениями:

1.2 = о) (X), 12р1=ар1 (Г) при 1 = 1,2,3.

Эти соотношения можно переписать иначе, если учесть определения тензоров временной корреляции Я. (*) иQ¡j (5) [3]

а)Са;

1 г ~ = * (°)=-/Н * л

■42Л -г

1 г___

Ои(} 011 (-)Са. л/2л -Г

II =

(13)

(14)

Выражая (13) через величину &")

согласно найденному соотношению (12) для интенсивности хаотического движения сплошной фазы имеем:

12 = 31 у[)Л

1 г ____

= } Г (а) <2 й (а)

(15)

<Ветшк&ТУИТ/Фгоиг£^ о/Т. 80, № 2, 2018-

Определим явный вид этого соотношения. И.О. Хинце с достаточной точностью выразил тензор временной корреляции следующей формулой:

в„ (*)=, (16) где X - некоторые постоянные.

Находим преобразование Фурье этой функции:

£ = II .. (17)

л; + ю

Пренебрегая силой Бассэ (последний член правой части уравнения (1)), для функции Г(ю) получим:

'2 ] + ( (18)

Г И =

Y-

9 ( 18v 4 + 1 Id 2m

Выражая в соотношении (16) величины Q" (ю) и Г(ю), согласно (17) и (18), после некоторых вычислений имеем:

V ^ А-1

f

11, г 2__т2

1+|Г+-

1 X Id2

27v

1+

Xtd2 12v

(19)

где IЛ и I . - интенсивности хаотического дви-

J1 Рг

жения сплошной (жидкой) и диспергированной (твёрдой) фаз.

ЛИТЕРАТУРА

1 Ильин В.А., Кукина А.В. Высшая математика. М.: Изд-во Проспект, 2004. 600 с.

2 Протодьяконов И.О., Люблинская И.Е., Рыжков А.Е. Гидродинамика и массообмен в дисперсных системах жидкость-твердое тело. Л.: Химия, 1987. 336 с.

3 Протодьяконов И.О., Сыщиков Ю.В. Турбулентность в процессах химической технологии. JL: Наука, 1983. 318 с.

4 Хинце И.О. Турбулентность. Ее механизм и теория. М.: Гос. изд. физико-математической литературы, 1963. 680 с.

5 Шишацкий Ю.И., Буданов АВ., Никель С.А., Власов Ю.Н. Влияние наложения низкочастотных механических колебаний на эффективность экстрагирования // Вестник ВГУИТ. 2018. № 1. С. 25 - 29.

6 Пищиков Г.Б., Лазарев В. А., Шихалев С.В. Метод оценки интенсивности пространственного смешения микроорганизмов в биореакторах непрерывного действия // Вестник Воронежского государственного университета инженерных технологий. 2017. №79(3). С. 169-173.

7 Celis С., da Silva L. F. F. Lagrangian mixing models for turbulent combustion: review and prospects // Flow, Turbulence and Combustion. 2015. V. 94. №. 3. P. 643-689.

8 Watanabe Т., Nagata K. Mixing model with multi-particle interactions for Lagrangian simulations of turbulent mixing //Physics of Fluids. 2016. V. 28. №. 8. P. 085103.

9 Li L. J. et al. A modified turbulent mixing model with the consideration of heat transfer between hot buoyant plume and sidewalls in a closed stairwell // International Journal of Heat and Mass Transfer. 2015. V. 84. P. 521-528.

10 Barmparousis C., Drikakis D. Multidimensional quantification of uncertainty and application to a turbulent mixing model // International Journal for Numerical Methods in Fluids. 2017. V. 85. №. 7. P. 385-403.

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

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

Физический смысл заключается в том, что если плотности твердой и жидкой фаз количественно разнятся, то следует ожидать некоторого увеличения относительной скорости их движения, что приведет к увеличению коэффициента массопередачи. Интенсивности хаотического движения фаз не будут одинаковыми, равно как и коэффициенты турбулентной диффузии [5-10]. Из сказанного следует, что рассмотренный случай движения одиночной частицы в турбулентном потоке сложен и может быть решен только при достаточно серьёзных допущениях, сформулированных выше.

REFERENCES

1 Ilyin V.A., Kukina A.V. Vysshaya matematika [Higher mathematics] Moscow, Publishing House Prospekt, 2004. 600 p. (in Russian)

2 Protodyakonov Ю., Lublinskaya IE., Ryzhkov A.E. Gidrodinamika I massoobmen v dispersnykh sistemakh [Hydrodynamics and mass transfer in liquid-solid disperse systems] Leningrad, Chemistry, 1987. 336 p. (in Russian)

3 Protodyakonov I.O., Syshchikov Yu.V. Tusrbulentnost' v protsessakh khimicheskoi tekhnologii [Turbulence in the processes of chemical technology] Leningrad, Nauka, 1983. 318 p. (in Russian)

4 Khintse I.O. Turbulentnost' [Turbulence. Its mechanism and theory] Moscow, State ed. physico-mathematical literature, 1963. 680 p. (in Russian)

5 Shishatsky Yul, BudanovA.V.,Nikel SA, VlasovYuN. Effect of the imposition of low-frequency mechanical oscillations on the extraction efficiency. Vestnik VGUIT. [Proceedings of VSUET] 2018. no. 1. pp. 25 - 29. (in Russian)

6 Pishchikov G.B., Lazaiev V.A., Shikhalev S.V. A method of evaluating the intensity of spatial mixing of the microorganisms in the bioreactors, continuous. Vestnik VGUIT. [Proceedings of VSUET] 2017. no. 79(3). pp. 169-173. (in Russian)

7 Celis C., da SilvaL. F. F. Lagrangian mixing models for turbulent combustion: review and prospects. Flow, Turbulence and Combustion. 2015. vol. 94. no. 3. pp. 643-689.

8 Watanabe T., Nagata K. Mixing model with multi-particle interactions for Lagrangian simulations of turbulent mixing. Physics of Fluids. 2016. vol. 28. no. 8. pp. 085103.

9 Li L. J. et al. A modified turbulent mixing model with the consideration of heat transfer between hot buoyant plume and sidewalls in a closed stairwell. International Journal of Heat and Mass Transfer. 2015. vol. 84. pp. 521-528.

10 Barmparousis C., Drikakis D. Multidimensional quantification of uncertainty and application to a turbulent mixing model. International Journal for Numerical Methods in Fluids. 2017. vol. 85. no. 7. pp. 385-403.

BeemHunJBry^T/Proceedings of VSUET, Т. 80, № 2, 2018-

СВЕДЕНИЯ ОБ АВТОРАХ Юлиан И. Шишацкий д.т.н., профессор, кафедра физики, теплотехники и теплоэнергетики, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия Алексей М. Барбашин к.т.н., доцент, кафедра физики, теплотехники и теплоэнергетики, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия Сергей А. Никель к.т.н., доцент, кафедра физики, теплотехники и теплоэнергетики, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия, sergei.nickel@yandex.ru

КРИТЕРИЙ АВТОРСТВА

Все авторы в равной степени принимали участие в написании рукописи и несут ответственность за плагиат КОНФЛИКТ ИНТЕРЕСОВ Авторы заявляют об отсутствии конфликта интересов.

ПОСТУПИЛА 30.03.2018 ПРИНЯТА В ПЕЧАТЬ 26.04.2018

INFORMATION ABOUT AUTHORS Yulian I. Shitshatskii Dr. Sci. (Engin.), professor, physics, heat engineering and heat power engineering department, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia

Aleksei M. Barbashin Cand. Sci. (Engin.), associate professor, physics, heat engineering and heat power engineering department, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia Sergei A. Nikel Cand. Sci. (Engin.), associate professor, physics, heat engineering and heat power engineering department, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia, sergei.nickel@yandex.ru

CONTRIBUTION

All authors equally participated in writing the manuscript and responsible for the plagiarism

CONFLICT OF INTEREST The authors declare no conflict of interest.

RECEIVED 3.30.2018 ACCEPTED 4.26.2018