Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

THE BASIC OF CALCULATION OF HYDRODYNAMIC PARAMETERS IN THE MOTION OF FLUID WITH MONODISPERSE OF COARSE SOLID PARTICLES IN VERTICAL PIPES

Kondratev A.S. 1 Nha T.L. 1
1 Moscow State University of Mechanical Engineering (MAMI)
Model calculation of two-phase flow is based on the assumption A. Fortier, that if the solid particle size approximately equal or exceeds the linear scale turbulence, the particles move along the trajectory, which does not affect the surrounding small turbulent fluid structure. Under this model the particle velocity is determined by the Lagrangian method, considering the process of movement established. Is accepted that the force of gravity, the Archimedes force, Magnus force and Safman force acts on the particles. The values of the Magnus forces and and the Safman forces calculated based on known literature correction functions. It is also assumed that the process of the radial distribution of solids is determined by the one-dimensional diffusion equation. The coefficient of turbulent diffusion of solid particles is defined as a product of the turbulent diffusion of fluid multiplied by a correction factor depending on the weight of the solid particles. Turbulent friction factor in the two-phase medium is determined by on the basis of the formula with the introduction Prandtl correction for the dispersed solids. The last two correction depending previously been successfully used in the calculation of the horizontal two-phase mixtures. As we considered monodisperse particles, as it was assumed that collisions between particles does not occur. The calculated ratio, for determining the local hydrodynamic flow parameters (concentration distribution of solids particles and liquid velocity), on the basis which is calculated of the integral characteristics of the two-phase flows (average concentration, the average velocity of the two-phase medium) during the motion of fluid mixtures and monodisperse large solid particles in vertical pipes. A comparison with published experimental data and the method of calculation, based representations of turbulent kinetic theory of granular flows using ANSYS CFX-4-4 software. The results of the comparison showed that the particle density of 2650 kg/m3 with a size of 470 microns with the solid phase volume fraction of 0,087 and 0,278 and an average speed of 2,6 m/s somewhat better agreement with experiment is observed when using the published method of calculation. For particles with a size of 1700 mm at a volume fraction of solids 0,085 and 0,177, and velocity 2,8 m/s the best agreement with experiment is observed when using the proposed method of calculation and the known method of calculation gives the dependence is qualitatively different from the experimental data.
monodisperse of coarse solid particles
the concentration of solids
friction loss
the volume flow
1. Zdanov V.G., Stapov V.M. Opredelenye effektivnoy vjazkosty konsentrirovannykh supenzjy // Kolloid. Zurnal. 1998. T. 60. no 6. pp. 771–778.
2. Kondratiev A.S. Dvizhenye monodispersnykh tviordykh chatics v potoke newtonovskoy zidkosti v vertykalnykh trubak. Sbornik. Statey. Problemy aksiomatiki v gidrogarodinamike. Vol. 21. M.: Kompanya Sputnik +. 2010, pp 15–24.
3. Kondpatiev A.S., Svydko. P.P. Opredelenye skorocty dvizenye gidrosmecsey v pripovierhnocthoy zone tetrenya. Sbornik. statey. Problemy aksiomatiki v gidrogarodinamike. 2014. no. 28. pp 344–354.
4. Kondpatiev A.S. Raschot raspredelenya konsentrasya polimodalnoy smesi tviordykh chactis proizvolnoy formy v newtohnovskoy zidkosti v gorizontalnykh trubak // Vesti. MGOU. Ser. Technol. 2013. no 2. pp 10–20.
5. Fortje A. Mekhanika suspenzjy. M.: Mir. 1971. 264 p.
6. Kharlamov A., Chara Z., Vlasak P. Experimental investigation of Magnus force acting on smooth sphere at high Reynolds number. Colloquium Fluid Dynamics. 2008. Inst. Thermomechanics AS CR. Prague. October 22–24. 2008. 11 p.
7. Krampa-Morlu F.N., Bergstrom D.J., Bugg J.D., Sanders R.S., Schaan J. Numerical Simulation of Dense Coarse Particle Slurry Flows in a Vertical Pipe. 5th ICMF. Yokohama. Japan. May 31 – June 3. 2004. pp. 1.
8. Krampa F.N. Two-Fluid Modelling of Heterogeneous Coarse Particle Slurry Flows. Univer. Saskatchewan. Doct. Disser. 2009. 242 p.
9. McLaughlin J.B. Inertial migration of a small sphere in linear shear flows // J. Fluid Mechan. 1991. V. 224. pp. 261–274.
10. Summer R.J., McKibben M.J., Shook C.A. Concentration and velocity distributions in turbulent vertical slurry flows // Ecoulements Solide-Liquide. 1990. V. 2(2). pp. 33–42.

В работе [3] построена физико-математическая модель движения мелких монодисперсных твердых частиц в потоке ньютоновской жидкости в вертикальных трубах, основанная на представлении, что линейный масштаб турбулентности λt существенно больше диаметра частиц d, то есть λt >> d [1]. В данной работе рассматривается другой предельный случай, а именно, принимается, что размер частиц порядка или более линейного масштаба турбулентности, то есть λt ≤ d [1]. Величину максимального линейного масштаба λt определим выражением, приведенным в [1]:

λt = 0,03R, (1)

где R – характерный размер потока в поперечном направлении – радиус трубы.

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

Силы Fp, FA, Fg и Fr определяются следующими выражениями:

kond01.wmf

kond02.wmf

kond03.wmf

kond04.wmf, (2)

где d – диаметр частицы; kond05.wmf – составляющая градиента давления, действующая на твердую частицу в направлении движения, в рассматриваемом случае вверх; r и rp – плотность жидкости и материала твердых частиц, соответственно; g – ускорение силы тяжести; Uf и Up – скорости движения жидкости и частиц в направлении движения потока вверх; f – коэффициент стесненности, учитывающий стесненное движение жидкости в межчастичном пространстве; cfu – коэффициент гидравлического сопротивления сферической частицы; коэффициента стесненности f определяется выражением [3]:

kond06.wmf

где j – объемная доля частиц твердой фазы; cfu – коэффициент гидравлического сопротивления сферической частицы.

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

kond07.wmf. (3)

Коэффициент гидравлического сопротивления cfu определяется по формуле Озеена:

kond08.wmf (4)

где kond09.wmf – число Рейнольдса, определенное по фактической скорости стесненного движения жидкости в межчастичном пространстве kond10.wmf, диаметру частицы d и кинематической вязкости двухфазного потока, определяемой по формуле

kond11.wmf

где kond12.wmf – кинематическая вязкость несущей жидкости; числа 0,65 и 1,675 – эмпирические коэффициенты [4].

Подставляя приведенные выражения в (3), получим

kond13a.wmf

kond13b.wmf (5)

где kond14.wmf – динамическая скорость, D – диаметр трубопровода.

Если число Рейнольдса значительно отличается от 4,5, то в выражении (5), вместо постоянного коэффициента 4,5 следует использовать величину kond15.wmf, и расчет по формуле (5) ведется путем итераций.

В работе [3] показано, что если использовать гипотезу Прандтля о распределении турбулентного трения с введением поправки на присутствие в потоке твердых частиц, то выражение для профиля скорости записывается в виде

kond16.wmf. (6)

где NΔ – относительная скорость двухфазной среды на расстоянии толщины приповерхностной зоны течения; Δ – толщина приповерхностной зоны движения двухфазной среды; y – расстояние от поверхности нижней образующей трубы;

kond17.wmf

длина пути смешения жидкости; y – расстояние от поверхности трубы; kpj – коэффициент, определяемый отношением усредненного квадрата пульсационной скорости частицы к усредненному квадрату пульсационной скорости жидкости [3] при характерной частоте, определяющей теоретическое значение нестационарной скорости осаждения частиц под действием силы тяжести [1].

Величина NΔ определяется зависимостью [5]

kond18.wmf, (7)

где kond19.wmf – среднее значение для поверхности трубы действительная степень шероховатости, которой неизвестна; φΔ – объемная доля частиц на расстоянии Δ от стенки трубы; φmax – максимальная объемная доля твердых частиц, соответствующая предельно шероховатой поверхности. Например, если сферические частицы и расположены в виде простой кубической решетки, то kond20.wmf. Если же частицы имеют форму откалиброванных песчинок, например, кубической формы, то φmax = 1. В формуле (7) второй член в правой части учитывает проскальзывание (скачок скорости) на «поверхности шероховатого квазитвердого тела», образованного движущимся придонным слоем твердых частиц ближайшего к поверхности, со стороны, обращенной к поверхности трубы [5].

Величина Δ определяется как сумма толщин подслоев жидкости на поверхности стенки трубы и диаметра частиц d: kond21.wmf где kond22.wmf.

В работе [3] показано, что при установившемся осесимметричном движении взвесенесущего потока, распределение объемной концентрации твердой фазы в поперечном направлении определяется выражением

kond23.wmf, (8)

где j0 – объемная концентрация твердой фазы на оси трубы; Vp – скорость перемещения частиц в горизонтальном направлении; Dp – коэффициент диффузии твердых частиц в жидкости.

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

kond24.wmf; (9)

где jm – средняя объемная концентрация частиц твердой фазы.

Выражение коэффициента турбулентной диффузии для частиц Dp можно представить в виде зависимости [2]

kond25.wmf (10)

где коэффициент диффузии жидкости Df определяется эмпирическим выражением

kond26.wmf. (11)

В поперечном, горизонтальном направлении на движущиеся частицы действует сила Магнуса:

kond27.wmf, (12)

где kond28.wmf – градиент скорости жидкости при стесненном обтекании твердых частиц; у – расстояние от поверхности; kM – поправочная функция для силы Магнуса.

Сила Сафмана также действует в горизонтально направлении:

kond29.wmf; (13)

где kS – поправочная функция для силы Сафмана.

На основе проведенного анализа рекомендуется использовать следующие поправочные зависимости.

Для величины kM для диапазона 0,5 < Rep < 140×103; 0,1 < Γ < 7 [6];

kond30.wmf, (14)

где a = 1,333 – 0,061×Г + 0,029×Г2, b = (5,9 + Г)/(39 + 4,6×Г), с = 1/(25 + 21×Г).

kond34.wmf поступательное число Рейнольдса частицы; kond35.wmf – вращательное число Рейнольдса; ω – угловая скорость вращения частицы относительно окружающей среды; в случае если жидкость в поперечном сечении потока не вращается, то принимается, что kond36.wmfа kond37.wmf – число вращения.

Для величины FS, при использовании в качестве параметра, определяющего диапазон использования рекомендуемых зависимостей, величины kond38.wmf, где kond39.wmf – число Рейнольдса частицы, определенное по скорости сдвига жидкости, используются следующие зависимости [9]:

kond40.wmf (15)

Из последней зависимости следует, что величина kS знакопеременна, а значит, направление действия силы Сафмана может изменяться на противоположное направление.

Сила гидродинамического сопротивления в горизонтальном направлении, подобно выражению, приведенному в (2), определяется выражением

kond41.wmf, (16)

где Vp – скорость движения частицы в горизонтальном направлении; cfv – коэффициент гидравлического сопротивления, определяемый по скорости Vpf; рассчитывается по формуле (4) с соответствующей заменой величины Рейнольдса Repu на kond42.wmf.

Скорость движения твердых частиц в поперечном направлении Vp определим из условия, что в поперечном направлении к оси потока на частицы действуют силы Магнуса FM и Сафмана FS противоположно направленной, силой гидродинамического трения Frv:

kond43.wmf. (17)

Отметим, что при записи правой части выражения (17) учитывалось, что в установившемся режиме при движении частиц от стенки трубы в направлении к оси, навстречу ему перемещается такой поток частиц за счет диффузии, вследствие чего результирующая скорость жидкости в поперечном направлении равна нулю, то есть принимается, что частицы перемещаются в неподвижной жидкости. Подставляя выражения (12)÷(16), с учетом (4), в (17) и разрешая его относительно Vp, получим

kond44.wmf

kond45.wmf (18)

Если второй член в аналоге формулы (4) существенно отличается от 4,5, то также используется рассмотренный выше итерационный процесс вычислений.

При известном распределении скорости жидкости, рассчитанном по формуле (6), рассчитываются распределения скоростей твердых частиц по формуле

kond46.wmf

kond47.wmf. (19)

Объемный расход жидкой и твердой фаз рассчитывается по формулам:

kond48.wmf, (20)

kond49.wmf (21)

Суммарный объемный расход взвесенесущего потока равен kond50.wmf.

Средняя объемная скорость равна

kond51.wmf. (22)

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

kond52a.wmf

kond52b.wmf (23)

В выражении (23) знак + относится к случаю движения двухфазного потока снизу вверх, а знак – к движению в обратном направлении.

При заданных градиенте давления и средней объемной доли твердых частиц и их размере приведенные соотношения позволяют рассчитать все локальные характеристики взвесенесущего потока монодисперсных твердых частиц, включая профили скорости жидкости и твердых частиц и поперечное распределение объемной доли частиц, и интегральные величины двухфазного потока: среднюю скорость и суммарную объемную долю твердых частиц. Порядок вычислений следующий. Задается произвольное значение концентрации твердой фазы, например, на оси потока, j0. Далее итеративным путем определяются все характеристики двухфазного потока на некотором малом расстоянии dy. После нахождения локальных значений скоростей и объемных долей частиц, определяются интегральные показатели двухфазного потока Um и jm. При необходимости расчет повторяется при новом значении j0. Расчет завершается при малом (допустимом) отличии значений jm от заданного значения.

В случае если движение взвесенесущего потока направлено вниз, то силы, определяемые градиентом давления в жидкости Fp и тяжести Fg, направлены вниз, а сила Архимеда FA – вверх, и этом случае частицы опережают жидкость, в соответствии с чем преобразуются приведенные выше выражения.

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

В рамках сформулированной проблемы проанализируем экспериментальные данные работы [10], в которой приведены опытные данные по движению смеси воды с частицами песка размером 470 мкм и 1700 мкм в трубопроводе с внутренним диаметром 0,040 м. В опытах измерялось распределение объемной доли частиц j и их скорости Up в радиальном направлении при различных средних объемных долях твердой фазы и при различных скоростях движения смеси. Величина максимального линейного масштаба турбулентности. λt, определенная по формуле (1), равна λt = 0,03×0,020 = 600 мкм. Поскольку в обоих случаях частицы имеют линейный размер, близкий или превосходящий максимальный линейный масштаб турбулентности λt, может быть использован рассмотренный метод расчета. В данных расчетах также принималось, что отсутствует механическое трение твердых частиц о поверхность стенки трубы.

В таблице приведены основные интегральные характеристики экспериментов и соответствующие расчетные величины, которые определялись по изложенной методике. Плотность материала частиц rp = 2650 кг/м3, жидкости ρ = 1000 кг/м3, кинематическая вязкость жидкости n = 10-6 м2/с.

Сопоставление приведенных опытных данных с расчетными профилями распределения объемной доли твердой фазы и скоростей частиц ранее выполнялось в работах [7, 8]. Используемый метод расчета основывался на использовании различных физических моделях турбулентности основанных на кинетической теории гранулированных потоков с использованием коммерческого программного продукта ANSYS CFX-4-4.

Основные интегральные характеристики экспериментов и соответствующие расчетные величины

№ п/п

d, мкм

Опыт

Расчет

φm

Um, м/с

φm

Um,м/с

Umах, м/с

1

470

0,087

2,6

0,0867

2,65

3,09

2

470

0,278

2,6

0,277

2,65

3,03

3

1700

0,085

2,8

0,0849

2,85

3,45

4

1700

0,177

2,8

0,176

2,77

3,28

kond1a.tif kond1b.tif

Рис. 1. Сравнение измеренных и результатов расчетов для частиц с d = 470 мкм и средней объемной долей частиц jm = 0,087: (a) скорости твердых частиц и (б) объемной доли частиц. Точки – опыт [10]; сплошная кривая – расчет [8]; пунктирная кривая – наш расчет

kond2a.tif kond2b.tif

Рис. 2. Сравнение измеренных и результатов расчетов для частиц с d = 470 мкм и средней объемной долей частиц jm = 0,278: (a) скорости твердых частиц и (б) объемной доли частиц. Точки – опыт [10]; сплошная кривая – расчет [8]; пунктирная кривая – наш расчет

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

Сравнение опытных данных по профилю скорости частиц Up показывает, что расчет по рассмотренной методике и методу работ [7, 8] достаточно хорошо совпадают. В обоих случаях расчетные распределения объемной доли твердой фазы j отличаются от экспериментальных данных в значительно большей мере. Опытные данные по всему поперечному сечению трубы близки к постоянному значению примерно равному среднему значению, φm = 0,087. Ни один из методов не отражает локальное уменьшение величины φ при значении y/R = 0,23. Метод расчета работ [7, 8] лучше согласуется с опытными данными, хотя и различие между изложенным методом и опытными значениями φ в области 0,2 ≤ y/R ≤ 1 не превосходит 15 %.

На рис. 2 приведено распределение объемной доли твердой фазы и скорости частиц для условий второго опыта.

Так же, как в предыдущем случае, сравнение опытных данных по профилю скорости частиц Up показывает, что расчет по рассмотренной методике и методу работ [7, 8] достаточно хорошо совпадает. Также имеется достаточно удовлетворительное соответствие между обоими расчетными методами и опытными значениями j. Ни один из методов не отражает локальное увеличение величины φ при значении y/R = 0,1.

kond3a.tif kond3b.tif

Рис. 3. Сравнение измеренных и результатов расчетов для частиц с d = 1700 мкм и jm = 0,085: (a) скорости твердых частиц и (б) объемной доли частиц. Точки – опыт [10]; сплошная кривая – расчет [8]; пунктирная кривая – наш расчет

kond4a.tif kond4b.tif

Рис. 4. Сравнение измеренных и результатов расчетов для частиц с d = 1700 мкм и средней объемной долей частиц jm = 0,177: (a) скорости твердых частиц и (б) объемной доли частиц. точки – опыт [10]; сплошная кривая – расчет [8]; пунктирная кривая – наш расчет

На рис. 3 приведено распределение объемной доли твердой фазы и скорости частиц для условий третьего опыта.

Сравнение опытных данных по профилю скорости частиц Up показывает, что расчет по рассмотренной методике и метод работ [7, 8] несколько лучше совпадает с опытными данными, хотя и расчет по рассмотренной методике достаточно хорошо совпадает с опытными данными. Сравнение опытных данных по распределению объемной доли твердой фазы j показало, что расчетные зависимости по методике работ [9, 10] в области 0,3 ≤ y/R ≤ 1 с погрешностью около 50 % отличаются от опытных данных. Качественный характер расчетной зависимости практически полностью повторяет аналогичные зависимости на рис. 1, а и 2, а. В этой области потока различия между опытными и расчетными значениями j по изложенной методике существенно меньше, хотя в области 0,9 ≤ y/R ≤ 1 также весьма значительно. Ни один из методов не отражает локальное увеличение величины j при значении y/R = 0,75.

На рис. 4 приведено распределение объемной доли твердой фазы и скорости частиц для условий второго опыта.

Также как в предыдущем случае сравнение опытных данных по профилю скорости частиц Up показывает, что расчет по рассмотренной методике и методы работ [7, 8] несколько лучше совпадают с опытными данными, хотя и расчет по рассмотренной методике достаточно хорошо совпадает с опытными данными.

Качественный характер расчетной зависимости j по методике работ [7, 8] практически полностью повторяет аналогичные зависимости на рис. 1, а, 2, а и 3, а и достаточно значительно отличается от опытных данных. Расчет по рассмотренной методике значительно лучше согласуется с опытными данными во всем диапазоне 0,05 ≤ y/R ≤ 0,1. Следует отметить, что и в этом случае ни один из методов не отражает локальное увеличение величины j при значении 0,6 ≤ y/R ≤ 0,75.

Проведенный сравнительный анализ показал, что рассмотренная методика дает лучшее соответствие с опытными данными с более крупными частицами, когда отношение величин (d/λt) = (1700/600) = 2,83 существенно больше единицы, хотя и для частиц меньшего диаметра (d/λt) = (470/600) = 0,783 результаты сравнения опытных и расчетных данных достаточно удовлетворительны. В работе [8] в связи с большим отличием опытных и расчетных данных по j, для частиц большого размера, отмечается, что развитая авторами модель расчета не в состоянии учесть особенности движения частиц большого размера. Эта констатация, как нам представляется, свидетельствует о том, что использование достаточно сложных моделей турбулентности двухфазного потока не обеспечивает необходимой общности, а, как и любая другая более простая модель турбулентности двухфазного потока, требует «настройки» на процессы, которые являются определяющими в том или ином конкретном случае.