Моделирование процесса распыления расплава в газовой струе

УДК 669-492

МОДЕЛИРОВАНИЕ ПРОЦЕССА РАСПЫЛЕНИЯ РАСПЛАВА В ГАЗОВОЙ СТРУЕ

П.А. Лыков, Р.М. Байтимеров, Е.В. Сафонов, А.О. Шульц

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

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

Увеличение числа областей применения порошковых материалов, обусловленное существенным экономическим эффектом от замены деталей, произведенных традиционными методами, на спеченные, привело к тому, что в начале XXI века ежегодный прирост продукции порошковой металлургии составлял 5 % в США и 13,5 % в странах Западной Европы. Одним из самых распространенных методов производства металлических порошков является распыление металлических расплавов в газовой струе.

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

1) изменением агрегатного состояния, изменением химического состава и свойств распыляемого материала в результате взаимодействия с распыляющим газом;

2) высокими скоростями процессов охлаждения и дробления металла;

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

Расчет

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

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

Для расчета потока воздуха использовалась стандартная модель турбулентности. Воздух

при этом рассматривался как идеальный газ.

Результаты моделирования скорости газового потока и давления распылительного узла

установки распыления металлов УРМ-001 (рис. 1) показали, что предложенная конструкция форсунки позволяет получать газовую струю с максимальной скоростью порядка 320 м/с (рис. 3).

Рис. 2. Схема распылительного узла установки распыления металлов УРМ-001:

1 - распылительная форсунка, 2 - канал подачи расплава, 3 - линия оси распылительной форсунки, 4 - периферийная линия оси канала подачи расплава

Рис. 3. График изменения скорости газовой струи на выходе из сопла: 1 - вдоль оси распылительной форсунки, 2 - вдоль периферийной линии оси канала подачи расплава

Минимальную скорость струи, при которой величина напора газа позволяет дробить капли размером ёк, определяли [2] из формулы

=

V

^ •Рг • ак

где О/ - коэффициент аэродинамического сопротивления капли; рг - плотность газа, кг/м3; Ж - скорость газовой струи, м/с; сж - коэффициент поверхностного натяжения металла, Н/м; йк - диаметр капли, м.

При определении критической скорости дробления учитывалось влияние изменения давления (рис. 4) и температуры (рис. 5) газа на его плотность.

Рис. 4. График изменения давления газовой струи на выходе из сопла: 1 - вдоль оси распылительной форсунки, 2 - вдоль периферийной линии оси канала подачи расплава

Title

" г-=

У1

/

I—I—,—,—,—| ,—,—,—,—| ,—,—,—,—|—,—,—.—, |—,—,—,—,—,—,—,—,—,—|

О 0,05 0,1 0,15 ОД 0,25 0,3

Рис. 5. График изменения температуры газа на выходе из сопла: 1 - вдоль оси распылительной форсунки, 2 - вдоль периферийной линии оси канала подачи расплава

Из рис. 6 видно, что дробление расплава до капель размером менее 50 мкм будет происходить при скорости струи газа выше 264 м/с.

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

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

потоком, достаточно велик, а время существования мало. Поэтому охлаждения металла в этот период не происходит [2].

Размер первоначально образующейся капли в условиях обдува струи расплава газовым потоком вычислялся по формуле

= 3Д

ж

где Лст - радиус струи металла, м; We - критерий Вебера,

=

Ж2 • Рг • ¿к

Время дробления капли радиусом г0, согласно [1], составит

I

V

т = ж

др

Рж • Г0

2 • ож

где Го - радиус капли, м; рж - плотность жидкого металла, кг/м3.

а

н

о

а

о

м

Л

и

н

и

о

а

о

и

Диаметр капли, мкм Рис. 6. Зависимость критической скорости струи газа от диаметра капли

Сферическая капля образуется из пряди, по форме близкой к цилиндру, при этом время сфе-роидизации капли составит [2]

1,785ж • Гк •Цж

Тсф =

где гк - радиус частицы после затвердевания, м; ц* - динамическая вязкость расплава, Нс/м2. Время охлаждения конечной капли до температуры кристаллизации [2]

См • Рж • ¿к_ !п 'м 'г 'пл - 'г

6а

где см - удельная теплоемкость металла, Дж/°С; 'м - температура металла, °С; 'г - температура распыляющего газа; 'пл - температура кристаллизации металла, °С; а - коэффициент теплоотдачи для конечной капли

а =

^г • Nu

¿к

где Хг - коэффициент теплопроводности газа, Вт/(м°С); Nu - критерий Нуссельта. Критерий Нуссельта [3] вычислялся по формуле № = 2 + 0,03 • Рг0,33Яе0,54 + 0,35 • Рг0,356Яе0,58, где Рг - критерий Прандтля; Яе - критерий Рейнольдса.

о

ж

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

Яе = Ж • ¿к •рг ;

ц

г

Рг = ,

^г

где цг - динамическая вязкость газа, Н с/м ; сг - удельная теплоемкость газа, Дж/°С.

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

^ = Ки-ХгТ , ^•

-(кр - 'г )Р

к

где ^кр - удельная теплота кристаллизации, Дж/кг; т - масса капли, кг; Е - площадь поверхности

2

капли, м .

Расчеты показали, что время дробления струи составило 7,95-Ю-4 с, время сфероидизации капли - 4,04 10-7 с, время охлаждения металла до температуры плавления - 4,48-Ю-4 с, а время выделения скрытой теплоты кристаллизации - 1,3110-3 с. Таким образом, отношение времени остывания металла к суммарному времени деформации составляет 2,21, что позволяет не учитывать фазовые переходы при моделировании процесса дробления.

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

Эксперимент

На основе расчетов и результатов моделирования была разработана схема установки распыления металлов УРМ-001 (патент №110312). Распыление производилось при следующих параметрах: температура металла 1250 °С, расход газа 2500 л/мин, давление наддува рабочей камеры плавильного агрегата 0,1 атм, удельный расход газа 1,5 л/г.

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

/

Рис. 7. Процесс распыления жидкого расплава на установке УРМ-001

35

20 60 100 140 180 220 260 300

Диаметр, мм

Рис. 8. Гранулометрический состав материала: 1 - расчет, 2 - эксперимент

Вывод

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

Представленная работа выполнена при поддержке Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы» (государственный контракт № 14.B37.21.0759).

Литература

1. ANSYS CFX-Solver Theory Guide: Release 13.0 //ANSYS, Inc. - 2010. - 390p.

2. Нечипоренко, О. С. Распыленные металлические порошки / О. С. Нечипоренко, Ю.И. Найда, А.Б. Медведовский. - Киев: Наук. Думка, 1980. - 240 с.

3. Кацнельсон, Б.Д. Исследование конвективного теплообмена между частицами и потоком в нестационарных условиях / Б.Д. Кацнельсон, Ф.А. Тимофеева // Тр. ЦКИТИ, кн. 12, 1949.

Лыков Павел Александрович. Инженер УНИД, Южно-Уральский государственный университет (Челябинск). Тел.: (351) 267-92-65; pavel2504@ya.ru.

Байтимеров Рустам Миндиахметович. Аспирант кафедры электромеханики и электромеханических систем, Южно-Уральский государственный университет (Челябинск). Тел.: 89507259421; rustam_shurmach@mail.ru.

Сафонов Евгений Владимирович. Кандидат технических наук, доцент, и.о. заведующего кафедрой «Двигатели летательных аппаратов», декан Аэрокосмического факультета, Южно-Уральский государственный университет (Челябинск). Тел.: (351) 267-92-65; e-safonov@yandex.ru.

Шульц Алексей Оттович. Инженер кафедры «Двигатели летательных аппаратов», ЮжноУральский государственный университет (Челябинск). Тел.: (351) 267-92-65; shao75@mail.ru.

Bulletin of the South Ural State University Series “Mechanical Engineering Industry” __________2013, vol. 13, no. 2, pp. 148-154

NUMERICAL SIMULATION OF ATOMIZATION OF MOLTEN METAL IN GAS SPRAY

P.A. Likov, South Ural State University, Chelyabinsk, Russian Federation, pavel2504@ya.ru,

R.M. Baitimerov, South Ural State University, Chelyabinsk, Russian Federation, rustam_shurmach@mail.ru,

E.V. Safonov, South Ural State University, Chelyabinsk, Russian Federation, e-safonov@yandex.ru, A.O. Schulz, South Ural State University, Chelyabinsk, Russian Federation, shao75@mail.ru

A mathematical calculation and numerical simulation of atomization of molten metal in gas spray have been carried out using ANSYS. Mathematical model is based on fact, that molten metal is considered as dispersed phase is gas flow. The results of simulation are validated with experimental data, a good agreement was found.

Keywords: atomization of molten metal, metal micropowder, numerical simulation.

Поступила в редакцию 17 мая 2013 г.