Numerical investigation of the vortex formation in a liquid metal under the action of disk agitator.pdf Известно, что измельчение структуры литого металла ведет к повышению его физико-механических свойств. Для измельчения структуры сплавов их подвергают модифицированию путем введения перед заливкой в литейную форму специальных добавок - микрочастиц, которые становятся дополнительными центрами кристаллизации, или изменяют поверхностное натяжение расплава на границе с зародышем кристаллизации. Для обеспечения равномерного распределения добавок по сплаву применяют специальные устройства-смесители, которые создают турбулентное поле течения внутри жидкого металла и тем самым интенсифицируют диффузионные процессы, способствующие равномерному распределению микрочастиц. Для организации процесса перемешивания в химической технологии нашли широкое применение смесительные устройства, состоящие из емкости, в которую помещают перемешиваемые компоненты (жидкости, сыпучие материалы и т.д.), и мешалки-завихрителя, являющиеся рабочим органом смесительного устройства [1-4]. По своему устройству и числу оборотов мешалки-завихрители разделяют на тихоходные: до 100 об/мин (лопастные, якорные, рамные, ленточные, шнеко-вые) и быстроходные: свыше 100 об/мин (пропеллерные, турбинные, дисковые, специальные лопастные) [2, 4]. Тихоходные завихрители создают, как правило, окружной поток жидкости, а быстроходные, в зависимости от количества и места расположения лопастей, дисков, пропеллеров, могут создавать радиальный поток, осевой поток или радиально - осевой поток жидкости. Выбор того или иного типа смесителя зависит от свойств смешиваемых жидкостей, их объема и времени перемешивания. Решение задачи гидродинамики процесса перемешивания позволяет судить о картине течения в смесительном устройстве и предсказать степень однородности получаемой смеси, а также определить параметры турбулентности, которые оказывают влияние на процесс перемешивания жидкости и микрочастиц. Авторами [5] проведено численное исследование влияния известных k-е-моделей турбулентности на прогнозируемые значения параметров потока в смешивающем устройстве и на основе сравнений с экспериментальными данными показано, что эти модели способны адекватно предсказывать параметры турбулентного потока. Исследование влияния несимметрично расположенной мешалки в смешивающем устройстве на гидродинамику процесса перемешивания было выполнено в работе [6] с помощью численного моделирования с привлечением реализованной Л-е-модели турбулентности. Авторам [7, 8] удалось получить хорошее соответствие между результатами численного моделирования с привлечением Л-е-модели турбулентности и экспериментальными данными при определении формы поверхности свободного вихря и при определении содержания газа в смесительном устройстве [8]. В данной работе рассматривается гидродинамика процесса перемешивания жидкого (расплавленного) металла в тигле под действием завихрителя, выполненного в виде трех перфорированных дисков с размещенными на них штифтами [9], рис. 1, а также исследуется влияние скорости вращения завихрителя на картину течения в тигле и основные параметры смесительного устройства. Рис. 1. Форма завихрителя: 1 - ось; 2-4 - перфорированные диски; 5, 6 - штифты; 7 - отверстия; 8 - штифты центрального диска Fig. 1. Agitator construction: 1, axis; 2-4, perforated disks; 5,6, pins; 7, openings; 8, pins of the central disk Математическая модель и метод решения Основные допущения, используемые при описании процесса перемешивания, следующие: жидкость несжимаемая; процесс перемешивания нестационарный, изотермический, турбулентный; течение осесимметричное и закрученное. Для описания положения границы жидкий металл - воздух используется модель VOF, для описания турбулентного течения привлекается классическая Л-е-модель. Определяющая система уравнений, записанная во вращающейся со скоростью за-вихрителя системе координат, описывает: - изменение объемной доли воздуха d02 + V(a2^r ) = 0; (1) - изменение объемной доли жидкого металла a1 = 1 -a2 ; (2) - изменение количества движения смеси ^ + V(pV.V) + proxV = -Vp + V-feff + FnH +Fg ; (3) dt s - изменение турбулентной кинетической энергии М+V(pkvr) = v dt ц + ^л | Vk ct, + Gk-ре; (4) - изменение скорости диссипации турбулентной кинетической энергии 2 ц+^^Ь. |V8 d^ + V(p8Vr) = V + ci8 jGk -C2.pk> (5) где р = р1а1 +p2a2 - плотность смеси; р1 - плотность жидкого металла; р2 -плотность воздуха; ю - угловая скорость вращения завихрителя; Vr = V - ю х r -скорость смеси в системе координат, связанной с завихрителем; r - расстояние от оси вращения до рассматриваемой точки; F^ = 2стк-Р-Va1 - сила поверхно- Pl +Р2 стного натяжения; к = V(n/|n|) - кривизна поверхности раздела фаз; n = Va1 -нормаль к поверхности раздела фаз; ст - коэффициент поверхностного натяжения; Fg =(-pg, 0) - сила тяжести; Teff = т + Tturb - эффективный тензор напряжений; Tturb = f-lturb (VV + VVT)- pkU - тензор турбулентных напряжений; т = ^(VV + VVT) - тензор вязких напряжений; ц = ц1а1 + ц2a2 - динамическая вязкость смеси; цturb = рСцk2 / 8 - турбулентной вязкость; Gk = 2цturb S : S - производство турбулентной кинетической энергии вследствие сдвигового течения (действие градиентов скорости); S = 2(vV + VVT) - тензор деформаций. Параметры в уравнениях (3), (4) следующие: Cl8 = 1.44; С28 = 1.92; Сц = 0.09; CTk = 1.0; ст8 = 1.3. На твердых поверхностях тигля и завихрителя задаются условия прилипания, на открытой границе тигля - «мягкие» граничные условия. Для решения системы уравнений (1), (3) - (5) используется метод конечных объемов, причем уравнение (1) решается по явной схеме с использованием метода геометрической реконструкции [10] для определения потоков на гранях ячеек, находящихся вблизи поверхности раздела фаз. Уравнения (3) - (5) решаются по неявной схеме, в которой для аппроксимации конвективных членов уравнений используется противопоточная разностная схема второго порядка аппроксимации. Для нахождения поля давления используется алгоритм PRESTO. Согласование поля давления и поля скоростей осуществляется по методу SIMPLE [11]. Значения параметров турбулентности на твердых поверхностях определялись в соответствии с моделью EWT (Enhanced Wall Treatment). Решение поставленной задачи осуществлялось с помощью пакета программ Ansys-Fluent [12]. Расчеты проводились при следующих значениях параметров: плотности жидкого металла и воздуха р1 = 2700 кг/м3, р2 = 1.2 кг/м3, динамическая вязкость жидкого металла и воздуха ц1 = 10-3 Пас, ц2 = 1.810-5 Пас, поверхностное натяжение жидкого металла ст = 0.9 Н/м, скорость вращения завихрителя n = 1003000 об/мин. В начальный момент времени высота уровня жидкого металла составляла половину высоты тигля. Высота тигля - 0.3 м, диаметр - 0.2 м. Обсуждение результатов Для определения основных безразмерных критериев, от которых будет зависеть решение системы уравнений (1) - (5), введем масштаб длины d - внешний диаметр завихрителя, масштаб времени ю-1, где ю - угловая скорость вращения завихрителя, масштаб плотности р1 - плотность жидкого металла. Тогда, согласно теории размерностей, масштаб давления выразится как p1(rod)2, масштаб скорости - rod. Запись уравнения (2) в безразмерном виде показывает, что безразмерными критериями рассматриваемой задачи являются: Fr = ro2d/g - число Фруда, We = р1ю2d3/ст - число Вебера, Re = p1rod2 /число Рейнольдса. Безразмерной энергетической характеристикой завихрителя служит модифицированное число Эйлера, Eu = P/(p^d5) представляющей отношение мощности завихрителя P к мощности вращающейся жидкости [1, 2]. К группе вышеперечисленных безразмерных критериев относится геометрический параметр - отношения характерных размеров завихрителя и тигля, d/D. При изменении скорости вращения завихрителя от 100 об/мин до 3000 об/мин число Re изменяется в диапазоне 1.8-105 - 5.4 106, число We - в диапазоне 1.7 102 - 1.5 105, число Fr - в диапазоне 0.9-805. Указанные значения параметров говорят о том, что инерционные силы намного превосходят силы вязкой диссипации и силы поверхностного натяжения, а также больше силы тяжести. Поэтому течение жидкого металла в тигле является существенно турбулентным. Влияние силы поверхностного натяжения и силы тяжести может играть существенную роль в областях, где течение заторможено - вблизи стенок тигля. Для исследования сходимости численного решения была проведена серия расчетов по определению момента сил, действующего на завихритель со стороны жидкого металла, с использованием четырех разностных сеток, размер ячеек которых уменьшался в два раза от 2 до 0.25 мм. Скорость вращения завихрителя задавалась равной 1500 об/мин. Система уравнений (1) - (5) решалась до тех пор, пока численное решение не сходилось к стационарному. Значения момента сил М в зависимости от размеров ячеек h и количества ячеек N разностной сетки приведены в таблице. За «точное» значение момента сил принималось значение, полученное на разностной сетке с размером ячеек 0.25 мм. Из таблицы следует, что имеет место сходимость численного решения, поскольку относительная погрешность в определении момента сил уменьшается при увеличении числа ячеек. Дальнейшие расчеты проводились на разностной сетке, число ячеек которой равно 119216, а размер ячеек 0.5 мм. Параметры Размер ячеек h, мм 2 1 0.5 0.25 N ячеек 8001 29804 119216 476864 М, Н-м 0.2016 0.2151 0.2211 0.2226 (AM/M>100 % 10.44 3.52 0.69 С ростом скорости вращения завихрителя глубина воронки за счет центробежных сил, отбрасывающих жидкий металл на периферию, увеличивается, рис. 2. При скорости вращения 500 об/мин поверхность жидкого металла достигает верхнего перфорированного диска, при скорости 1000 об/мин - среднего диска. При скорости вращения 1500 об/мин и выше поверхность жидкого металла достигает нижнего перфорированного диска, и площадь поверхности соприкосновения жидкого металла с элементами завихрителя уменьшается. Те штифты, которые соприкасаются как с жидким металлом, так и с воздухом, способствуют проникновению воздуха в жидкий металл, создавая предпосылки для кавитации. Рис. 2. Положение уровня жидкого металла (1 - n = 500; 2 - 1000; 3 - 1500; 4 - n =3000 об/мин) Fig. 2. Level position of a liquid metal (n = (1) 500, (2) 1000, (3) 1500, and (4) 3000 rpm) Для скоростей вращения завихрителя от 100 до 1500 об/мин течение жидкого метала в тигле выходит на установившийся режим. При этом минимальное значение момента сил и его максимальное значение для заданной скорости вращения завихрителя слабо отличаются друг от друга, рис. 3. При скоростях вращения свыше 1500 об/мин течение становится существенно нестационарным и момент сил изменяется в пределах от минимального значения (обозначено крестиками) до своего максимального значения (обозначено треугольниками). Следует отметить, что учет силы поверхностного натяжения практически не оказал существенного влияния на значение момента сил завихрителя (кружочки). На рис. 4 показаны характерные картины течения жидкого металла при различных скоростях вращения завихрителя. При скорости вращения 500 об/мин жидкий металл полностью покрывает завихритель и в окрестности штифтов завихрителя образуются шесть мелких вихрей, а вблизи боковой стенки - один крупный (рис. 4, а). При увеличении скорости вращения до 1000 об/мин жидкость наполовину закрывает завихритель и количество вихрей в окрестности штифтов сокращается до четырех, а в окрестности боковой стенки вихрь распадается на два (рис. 4, b). Дальнейшее увеличение скорости вращения ведет к тому, что поверхность соприкосновения жидкого металла с завихрителем уменьшается, и внутри области формируются два относительно крупных вихря (рис. 4, c-f). 01 X2 Д3 04 -A Д t i V 1 / ft J X ч > / /4 i 1 ***** 1000 2000 n, об/мин 3000 0 Рис. 3. Зависимость момента сил от числа оборотов завихрителя (1 - без учета сил поверхностного натяжения; 2 - нижняя граница; 3 - верхняя граница; 4 - с учетом сил поверхностного натяжения) Fig. 3. The force moment as a function of the agitator revolution number (1, excluding the surface tension forces; 2, bottom boundary; 3, top boundary; 4, taking into account the surface tension forces) 0.5 0.4 S 0.3 к ^ 0.2 0.1 0 Рис. 4. Линии тока в жидком металле (а - n = 500; b - 1000; с - 1500; d - 2000; e - 2500; f - 3000 об/мин) Fig. 4. Stream lines of a liquid metal, n = (a) 500; (b) 1000; (с) 1500; (d) 2000; (e) 2500; and f) 3000 rpm Оценка размеров вихрей, вызванных турбулентными пульсациями, показывает, 3 I 0 25 1 5 что их значение, изменяющееся от колмогоровского масштаба (v /е) до k / е [13], не превосходит 1 см и убывает с ростом скорости вращения завихрителя до 2000 об/мин, рис. 5. Размер крупных турбулентных вихрей остается неизменным при увеличении скорости вращения завихрителя от 2000 до 3000 об/мин и лежит в пределах от 5 до 6 мм. 10 9 8 I 7 i 6 5 4 1000 2000 и, об/мин 0 3000 Рис. 5. Зависимость размера турбулентных вихрей от числа оборотов завихрителя Fig. 5. Size of turbulent eddies as a function of revolution number of the agitator Проведенные расчеты показывают, что увеличение скорости вращения завихрителя приводит к росту момента сил M, действующего на него со стороны жидкого металла, и росту турбулентной диффузии жидкого металла D. После достижения скорости вращения 1500 об/мин темпы роста турбулентной диффузии снижаются, что связано с уменьшением поверхности соприкосновения завихрителя и жидкого металла, рис. 6. Увеличение скорости вращения завихрителя ведет к падению работы внешних сил по отношению к кинетической энергии вращающегося жидкого металла (рис. 7), описываемое следующей зависимостью критерия мощности от числа Рейнольдса Eu = 1.183exp(-4.298-10-7Re). и, об/мин Рис. 6. Зависимость турбулентной диффузии от числа оборотов завихрителя Fig. 6. Turbulent diffusion as a function of revolution number of agitator Рис. 7. Зависимость критерия мощности от числа Рейнольдса Fig. 7. Power number vs Reynolds number Re-10-6 Заключение Показано, что увеличение скорости вращения завихрителя ведет к интенсивному росту турбулентной диффузии в жидком металле до тех пор, пока поверхность раздела металл - воздух не начнет взаимодействовать с нижним перфорированным диском завихрителя. Критерий мощности экспоненциально убывает с ростом числа Рейнольдса (скорости вращения завихрителя). Условия, при которых элементы завихрителя погружены полностью в жидкий металл, более благоприятны для формирования мелких вихрей, чем условия, при которых поверхность раздела жидкий металл - воздух взаимодействует с элементами завихрителя. Следует ожидать, что самые крупные вихри, сопоставимые по величине с размерами завихрителя, способны неравномерно перемешивать жидкий металл и внедренные в него микрочастицы, тогда как мелкие вихри, хотя и энергетически более слабые, могут обеспечить равномерное перемешивание, создавая участки равномерности в несколько десятков микрон.

Скачать электронную версию публикации