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

МЕХАНИКА

УДК 533:532

Н. Л. Великанов, В. А. Наумов, С. И. Корягин

ТЕЧЕНИЕ ГАЗА С ДИСПЕРСНОЙ ПРИМЕСЬЮ В ЦИЛИНДРИЧЕСКОМ КАНАЛЕ ПРИ ДОЗВУКОВОЙ СКОРОСТИ

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

Задача определения динамических характеристик процессов подачи инертных газов, содержащих дисперсные примеси, в ковши и конвертеры используется для повышения точности дозирования инертных газов. Соответствующая краевая задача не имеет аналитического решения. Рассмотрена динамика стационарного течения совершенного газа, содержащего дисперсные примеси, в вертикальном канале. Представлены результаты решения краевой задачи численным методом в среде Mathcad. Показаны итоги исследования изменений параметров газа вдоль канала при различных расходах. Графики представлены в размерных переменных. В частности, исследованы распределенная плотность дисперсной фазы, температура несущей и дисперсной фазы, числа Маха M, Рейнольдса Re, Нуссельта Nu, Стокса St.

In this article, we present a method and an algorithm for calculating the dynamic characteristics of gas flow in a cylindrical channel at a subsonic speed, in the presence of dispersed particles in the gas.

The problem of identifying the dynamic characteristics of dispersed particle-containing inert gas feed into ladles and the converter is used to improve the accuracy of inert gas dosing. The relevant boundary value problem has no analytical solution. We consider the dynamics of a steady flow of a perfect gas containing dispersed particle in a vertical channel. Numerical methods are used to obtain a solution to the boundary value problem in the Mathcad environment. We present the results of a study of changes in gas parameters along the channel at different gas flow rates. The graphs show dimensional variables. We consider the distributed density of the dispersed phase, the temperatures of the carrier and dispersed phases, the Mach number M, the Reynolds number Re, the Nusselt number Nu, and the Stokes number Stk.

Ключевые слова: дисперсная примесь, число Рейнольдса, течение газа.

Keywords: dispersed particles, Reynolds number, gas flow.

Ведение

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

81

© Великанов Н. Л., Наумов В. А., Корягин С. И., 2018

Вестник Балтийского федерального университета им. И. Канта. Сер.: Физико-математические и технические науки. 2018. № 2. С. 81 — 91.

82

гии, подаваемой для перемешивания металла в ковше падающим потоком и нагнетаемым газом, должно поддерживаться на непрерывно высоком уровне. Одним из способов достижения этого является увеличение количества закачиваемого газа в течение второго периода работы отвода [3].

Новые исследования в области производства чугуна и стали должны быть ориентированы на компактные, гибкие и экономичные технологии, обеспечивая при этом более высокое качество конечного продукта [4].

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

Затронуты пять направлений [4]:

1) производство стали из жидкого чугуна;

2) производство стали из лома или предварительно восстановленной руды;

3) новые технологии плавки лома;

4) прямое восстановление;

5) новые технологии.

Ведется работа по усовершенствованию существующих доменных печей с использованием вдувания пылеугольного топлива с целью сокращения потребления кокса. Фактические ограничения не способствуют созданию новых коксовых мощностей [4].

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

Как показали авторы статьи статьи [8], для оптимизации химических реакторов применима вычислительная гидродинамика.

В работе [9] на основе уравнения Пуассона — Больцмана в цилиндрических координатах рассчитана проводимость одной заряженной нанотрубки, заполненной электролитом.

Экспериментально исследованы динамика и параметры бифуркации ударной волны, обусловленные взаимодействием отраженной ударной волны с пограничным слоем, в аргоне, воздухе и водородно-азотной смеси для чисел Маха 1,3—3,5 в ударной трубке диаметром 76 мм. Измерения проводились при постоянной плотности газа за отраженной ударной волной [10].

1. Исходные уравнения математической модели

Ранее [1] нами было рассмотрено установившееся течение идеального газа в вертикальном канале длины L, диаметра D, L ^ D, а также установлена температура стенки Tw канала. Поток предполагался высокоскоростным, дозвуковым, давление на выходе газа равно давлению среды (обратному давлению). Силы гравитации не учитывались.

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

Уравнение сохранения массового расхода газа G1 с учетом постоянства плошади поперечного сечения трубы S:

Gi = pfWfS = const,

G

Gf = —S" = pfWf = const, где pf, Wf — средние по сечению плотность и скорость газа соответственно; S = %D ; D — диаметр трубы.

Уравнение состояния совершенного газа с учетом Tf < 1000 К:

p = PfRTf,

где R — газовая постоянная, равная универсальной газовой постоянной, деленной на молекулярную массу; Tf — термодинамическая температура газа (средняя по сечению трубы); Р — давление. Координата X направлена по оси трубы в направлении течения.

Система уравнений неизотермического течения газовзвеси в канале (одномерное приближение) при малой объемной доле дисперсной фазы имеет следующий вид [1; 7]:

dWf dP Wf

p Wit=-dxf p W -F,

dTf dWf Wf pCVWt—L = -P-f- + Cfpf — Wf + QW -Q,

f fV f dX dX f f 2D f W

dW Wf

p W-f = -C p —- + F,

f dX 4fPf 2D

83

йТ„

р С Щ —- = о,

Р- - - ах

где Щ-, Т- — скорость и термодинамическая температура дисперсной фазы; ¥ — сила межфазного взаимодействия (на единицу объема); О — интенсивность межфазного теплообмена в единице объема; С/ — коэффициент гидравлических потерь на трение.

2. Расчетные формулы

__Так как длина трубы Ь достаточно велика: Ь/О ^ 100, — можно не

84 учитывать эффект начального участка трубы и рассчитывать коэффициент гидравлических потерь на трение по формуле Блазиуса:

С / = 0,3164 Ие"0,25. Тепловой поток через стенки трубы рассчитывается по формуле [8]

Qw = DРfCPWf (Tw - Tf),

где С/р — теплоемкость газа при постоянном давлении.

Число Стэнтона БЬ в трубах при числах Маха М < 0,8 можно найти по формуле А. А. Гухмана [1]:

St = 0,0167(Ие Рг)"0,18 ,

ТЩ

где ТЩ — термодинамическая температура стенки трубы; Тм — средняя

по сечению трубы температура торможения [2]:

Щ/

Тм = Т/ + 0,5 С-.

С/р

При одномерном движении газовзвеси можно рассчитывать силу межфазного взаимодействия, умножая силу каждой из частиц на их счетную концентрацию, по формуле [1; 7]:

F = 3*Svpff(Re5)(Wf -Wr)np, ß = ^Г, . (1)

pp

В англоязычной научной литературе, как правило, стандартную кривую сопротивления аппроксимируют зависимостью [1; 11]

С0 = <!

R

24

-(1 + 0,15Re0'687) приReg< 1000,

Res (2)

0,44 при 1000 < Reg < 2 • 105.

Из (2) следует, что функция, учитывающая отклонение закона гидродинамического сопротивления в (1) от линейного для Reg < 1000, есть

f (Reg) = 1 + 0,15 Re

0,687

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

дисперсной фазы и объем частицы ©1

n _р-/p-

- ©1 л53/б'

Сила межфазного взаимодействия для Res < 1000 может быть вычислена по формуле

F _ ß/(Re5)(Wf -W;)Р;, ß_1|^Г, X_Pf.

Динамический коэффициент формы частиц будем рассчитывать, как в [12; 13]:

Г1(§), при Re < 0,2,

Г3(§, Re), при 0,2 < Re < 1162,

Г2(§), при Re > 1162,

r&Re) _■

Г Д) _ 1 + 0,348(2, -1), ВД _ 10,0 - 9,0/

85

Г3(1, Re) _ГД) + (ГД)-ГД))

Re-0,2 1161,8

0,44

где % — геометрический коэффициент формы. Интенсивность межфазного теплообмена

Q _ßT (T/ -T; )Р; , ßT _

6Xv Nu

V^-' К - 52 pr

Без учета сжимаемости число Нуссельта Ми частицы можно рассчитать по формуле Маршалла — Ранца [1]:

Ми = 2 + 0,6Ке0'5Рг0'33.

3. Окончательные уравнения математической модели

При числах Маха 0,1 < М < 0,7 можно воспользоваться следующим приближением [7]:

ХТ Nu А„ \Wf - W;\

NuM _-, М5 -

M 1 + 3,42M5Nu 5 a

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

W-dX = ßT(Tf -T;),

dW; W- (G/T/R )

W-- _-C —— + ß/(Re5)

- dX - 2D 5

-W P -

(3)

86

Форма уравнений неизотермического течения газа из-за обменных членов изменится следующим образом:

С

fp

( GlZR'] (GtR} 1 —

-7 /

v P

Tf

f G)TfR 2D P

+ F

f GfR 2

CfP + Tf

V p V

dP _ dX ~

4GfR ~(Qw - Q )f *

с

fp

G 2TfR 1--f—f—

( GR \

v P v

dT t

f

(GfTfR)2 (G2fTfR С

+ f

P 2D

Qw - Q

dX

/

G 2fTfR л

1 —Ц—

Граничные условия системы (3) — (6):

T(0) _ To, Wp(0) _ Wpo, Tp(0) _ Tpo, P(L) _ PK.

(5)

(6)

(7)

4. Результатах расчетов

Краевая задача (3) —(7) решалась численным методом в среде Mathcad. На рисунках 1—12 представлены результаты расчетов при различных расходах дисперсной примеси и следующем базовом наборе значений параметров:

G1 _ 0,06 кг/с, То _ 300 К, Tw _ 800 К, L _ 5 м, D _ 0,02 м, б _ 0,2 мм, | _ 1,3; Рк _ 3 105 Па, р° _ 2600 кг /м3, Wp0 _ 1 м/с, Tp0 _ 300 К.

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

•

—-

"" ~ <- - 1 * ~ , • *

кПа

420

380

340

300

0.5

1.5

2.5

3.5

4.5 X.

Рис. 1. Изменение давления вдоль канала при различных расходах частиц: 1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,12 кг/с; 4 - Сг = 0,18 кг/с

Распределенная плотность (концентрация) дисперсной фазы на рисунке 2 уменьшается вдоль канала из-за роста скорости частиц (см. рис. 4). При этом как скорость, так и температура у дисперсной фазы заметно меньше, чем у несущей фазы (см. рис. 11 и 12).

Рис. 2. Изменение распределенной плотности дисперсной фазы

вдоль канала при различных расходах С2: 1 - С2 = 0,0001 кг/с; 2 - С2 = 0,06 кг/с; 3 - С2 = 0,12 кг/с; 4 - С2 = 0,18 кг/с

Щ

м/с 110

80

50

20

Л-'*

Г1*

0

0.5

1

1.5

2.5

3.5

4.5 X.

Рис. 3. Изменение скорости несущей фазы вдоль канала при различных расходах С2: 1 - С2 = 0,0001 кг/с; 2 - С2 = 0,06 кг/с; 3 - С2 = 0,12 кг/с; 4 - С2 = 0,18 кг/с

87

Щ

м/с 75

50

25

0

У--'-'

'—* 1

0.5

1

1.5

2.5

3.5

4.5 Х.1

Рис. 4. Изменение скорости дисперсной фазы вдоль канала при различных расходах С2: 1 - С2 = 0,0001 кг/с; 2 - С2 = 0,06 кг/с; 3 - С2 = 0,12 кг/с; 4 - С2 = 0,18 кг/с

О 0.5 1 1.5 2 2.5 3 3.5 4 4.5 X, м Рис. 5. Изменение чисел Маха вдоль канала при различных расходах С2:

1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,12 кг/с; 4 - Сг = 0,18 кг/с

Рис. 6. Изменение чисел Стокса вдоль канала при различных расходах Сг: 1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,1г кг/с; 4 - Сг = 0,18 кг/с

Рис. 7. Изменение чисел Нуссельта вдоль канала при различных расходах С2: 1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,^ кг/с; 4 - Сг = 0,18 кг/с

89

Рис. 8. Изменение чисел Рейнольдса частиц вдоль канала при различных расходах С2: 1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,12 кг/с; 4 - Сг = 0,18 кг/с

Рис. 9. Изменение чисел Маха (относительных) вдоль канала при различных расходах С2: 1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,12 кг/с; 4 - Сг = 0,18 кг/с

Рис. 10. Изменение динамического коэффициента формы частиц

вдоль канала при различных расходах С2: 1 - Сг = 0,0001 кг/с; 2 - Сг = 0,06 кг/с; 3 - Сг = 0,12 кг/с; 4 - Сг = 0,18 кг/с

90

620

540

460

380

300

*

_ *■

о

0.5

1

1.5

2.5

3.5

4.5 X.

Рис. 11. Изменение температур несущей (Т) и дисперсной (Тр) фаз вдоль канала

при расходе С2 = 0,06 кг/ с

щ.

м/с 120

90 60 30

Щ

, # * * » — ^" в ^ • ^

Г

о

0.5

1

1.5

2.5

3.5

4

4.5 X.

Рис. 12. Изменение скоростей несущей (И/) и дисперсной (Ир) фаз вдоль канала

при расходе С2 = 0,06 кг/ с

Заключение

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

Список литературы

1. Великанов Н. Л., Наумов В. А., Корягин С. И. Течение газа в цилиндрическом канале при дозвуковой скорости // Вестник Балтийского федерального университета им. И. Канта. Сер.: Физико-математические и технические науки. 2018. № 1. С. 96-103.

2. Корнеев С. В. Применение инжекционных технологий в металлургическом производстве // Литье и металлургия. 2011. № 2. С. 152—159.

3. Svyazhin A. A., Krushke E. Mixing of the metal and mass transfer during the tapping of a semifinished product from an oxygen converter // Metallurgist. 2005. Vol. 49, iss. 11-12. P. 432-438.

4. Faure H. New technologies for the production of iron and steel // Revue de métallurgie. Cahiers d'informations techniques. 1994. Vol. 91, iss. 1. P. 49 — 60.

5. Koria S. C. Principles and applications of gas injection in steelmaking practice // Scandinavian journal of metallurgy. 1993. Vol. 22, iss. 5. P. 271 — 279.

6. Харлашин П. С., Мохаммед А. К., Харин А. К., Куземко Р. Д. Влияние концентрации порошка на течение газовой взвеси в торкрет-фурме 160-т конвертера // Сталь. 2015. № 4. С. 21 — 25.

7. Валуев Д. В. Внепечные и ковшовые процессы обработки стали в металлургии. Томск, 2010.

8. Svyazhin A. A., Krushke E. Mixing of the metal and mass transfer during the tapping of a semifinished product from an oxygen converter // Metallurgist. 2005. Vol. 49, iss. 11 — 12. P. 432—438.

9. Faure H. New technologies for the production of iron and steel // Revue de métallurgie. Cahiers d'informations techniques. 1994. Vol. 91, iss. 1. P. 49 — 60.

10. Koria S. C. Principles and applications of gas injection in steelmaking practice / / Scandinavian Journal of Metallurgy. 1993. Vol. 22, iss. 5. P. 271 — 279.

11. Crowe C. T., Sommerfeld M., Tsuji Y. Multiphase Flows with Droplets and Particles. Florida, 1998.

12. Sommerfeld M., Huber N. Experimental analysis and modelling of particle-wall collisions // International Journal of Multiphase Flow. 1999. Vol. 25. P. 1457—1489.

13. Великанов Н. Л., Наумов В. А., Космодамианский А. С. Моделирование осаждения твeрдых частиц в пульпопроводе // Наука и техника транспорта. 2011. № 2. С. 69 — 78.

14. Великанов Н. Л., Наумов В. А., Примак Л. В. Осаждение частиц взвесей в воде // Механизация строительства. 2013. № 7. С. 44 — 48.

Об авторах

Николай Леонидович Великанов — д-р техн. наук, проф., Балтийский федеральный университет им. И. Канта, Россия.

Е-mail: monolit8@yandex.ru

Владимир Аркадьевич Наумов — д-р техн. наук, проф., Калининградский государственный технический университет, Россия.

Е-mail: van-old@rambler.ru

91

Сергей Иванович Корягин — д-р техн. наук, проф., Балтийский федеральный университет им. И. Канта, Россия. Е-mail: SKoryagin@kantiana.ru

The authors

Prof. Nikolay L. Velikanov, I. Kant Baltic Federal University, Russia. E-mail: monolit8@yandex.ru

Prof. Vladimir A. Naumov, Kaliningrad state technical university, Russia. E-mail: van-old@rambler.ru

Prof. Sergey I. Koryagin, I. Kant Baltic Federal University, Russia. E-mail: SKoryagin@kantiana.ru