Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

ОСНОВЫ РАСЧЕТА ГИДРОДИНАМИЧЕСКИХ ПАРАМЕТРОВ ПРИ ДВИЖЕНИИ ЖИДКОСТИ С МОНОДИСПЕРСНЫМИ КРУПНЫМИ ТВЕРДЫМИ ЧАСТИЦАМИ В ВЕРТИКАЛЬНЫХ ТРУБАХ

Кондратьев А.С. 1 Ньа Т.Л. 1
1 Московский госудаственный машиностроительный университет (МАМИ)
Модель расчета двухфазного потока основывается на предположении А. Фортье, что если размер твердых частиц примерно равен или превосходит линейный масштаб турбулентности потока, то частицы движутся по траекториям, на которые не влияет мелкая турбулентная структура окружающей жидкости. В рамках этой модели скорость частицы определяется по методу Лагранжа, считая процесс движения установившимся. Принимается, что на частицы действует сила тяжести, сила Архимеда, сила Магнуса и сила Сафмана. Величины сил Магнуса и Сафмана рассчитываются с учетом известных по литературе поправочных функций. Принимается также, что процесс радиального распределения твердых частиц определяется одномерным уравнением диффузии. Коэффициент турбулентной диффузии твердых частиц определяется как произведение турбулентной диффузии жидкости, умноженной на поправочный коэффициент, зависящий от массы твердой частицы. Турбулентное трение в двухфазной среде определяется на основе формулы Прандтля с введением поправки на влияние дисперсной твердой фазы. Последние две поправочные зависимости ранее успешно использовались при расчете горизонтальных двухфазных смесей. Поскольку рассматривались монодисперсные частицы, предполагалось также, что столкновения между частицами не происходит. Получены расчетные соотношения, позволяющие определить локальные гидродинамические параметры потока (распределение концентрации твердой фазы, скорости частиц и жидкости), на основании которых рассчитываются интегральные характеристики двухфазного потока (средняя концентрация, средняя скорость движения двухфазной среды) при движении смесей жидкости и монодисперсных крупных твердых частиц в вертикальных трубах. Проведено сравнение с опубликованными экспериментальными данными методом расчета, основаным на представлениях кинетической теории турбулентных гранулированных потоков с использованием коммерческого программного продукта ANSYS CFX-4-4. Результаты сравнения показали, что для частиц плотностью 2650 кг/м3 с размером 470 мкм при объемной доле твердой фазы 0,087 и 0,278 и средней скорости 2,6 м/с несколько лучшее соответствие с опытом наблюдается при использовании опубликованного метода расчета. Для частиц с размером 1700 мкм при объемной доле твердой фазы 0,085 и 0,177 и скорости 2,8 м/с лучшее соответствие с опытом наблюдается при использовании предложенного метода расчета, а известный метод расчета дает зависимости, качественно отличающиеся от опытных данных.
монодисперсные крупные твердые частицы
концентрация твердой фазы
потери на трение
объемный расход
1. Жданов В.Г., Старов В.М. Определение эффективной вязкости концентрированных суспензий // Коллоид. журн. – 1998. – Т. 60. № 6. – С. 771–778.
2. Кондратьев А.С. Движение монодисперсных твердых частиц в потоке ньютоновской жидкости в вертикальных трубах. Сб. ст. Проблемы аксиоматики в гидрогазодинамике. Вып. 21. – М.: Компания Спутник +, 2010. – С. 15–24.
3. Кондратьев А.С., Швыдько П.П. Определение скорости движения гидросмесей в приповерхностной зоне течения. Сб. ст. Проб. аксиомат. в гидрогазодин. – 2014. – № 28. – С. 344–353.
4. Кондратьев А.С. Расчет распределения концентрации полимодальной смеси твердых частиц произвольной формы в ньютоновской жидкости в горизонтальных трубах // Вестн. МГОУ. Сер. Техн. и технол. – 2013. – № 2. – С. 10–20.
5. Фортье А. Механика суспензий. – М.: Мир, 1971. – 264 с.
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. – Р. 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. – P. 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). – P. 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, для частиц большого размера, отмечается, что развитая авторами модель расчета не в состоянии учесть особенности движения частиц большого размера. Эта констатация, как нам представляется, свидетельствует о том, что использование достаточно сложных моделей турбулентности двухфазного потока не обеспечивает необходимой общности, а, как и любая другая более простая модель турбулентности двухфазного потока, требует «настройки» на процессы, которые являются определяющими в том или ином конкретном случае.


Библиографическая ссылка

Кондратьев А.С., Ньа Т.Л. ОСНОВЫ РАСЧЕТА ГИДРОДИНАМИЧЕСКИХ ПАРАМЕТРОВ ПРИ ДВИЖЕНИИ ЖИДКОСТИ С МОНОДИСПЕРСНЫМИ КРУПНЫМИ ТВЕРДЫМИ ЧАСТИЦАМИ В ВЕРТИКАЛЬНЫХ ТРУБАХ // Фундаментальные исследования. – 2016. – № 9-1. – С. 35-42;
URL: http://fundamental-research.ru/ru/article/view?id=40691 (дата обращения: 25.01.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074