CFD МОДЕЛИРОВАНИЕ СГУСТИТЕЛЯ КРАСНЫХ ШЛАМОВ

Научная статья
DOI:
https://doi.org/10.18454/IRJ.2016.54.011
Выпуск: № 12 (54), 2016
Опубликована:
2016/12/19
PDF

Федорова Э.Р.

Аспирант, Санкт-Петербургский Горный Университет, Санкт-Петербург

CFD МОДЕЛИРОВАНИЕ СГУСТИТЕЛЯ КРАСНЫХ ШЛАМОВ

Аннотация

Разработана CFD модель сгустителя, описывающая поведение флокулированной суспензии по всему объему сгустителя. В модель дополнительно были заложены функции, описывающие стесненное осаждение и эффективное сжатие твердого вещества при помощи встроенных и добавленных пользователем моделей через UDF. Трехмерная модель учитывает распределение частиц по размерам в питающей пульпе сгустителя и действие перемешивающих граблин.  Результатом данной работы является CFD модель сгустителя, разработанная в программном пакете ANSYS Fluent.

Ключевые слова: сгуститель, красный шлам, CFD моделирование, стесненное осаждение, эффективное сжатие, распределение частиц по размерам.

Fedorova E.R.

Postgraduate student, Saint-Petersburg Mining University, Saint-Petersburg

CFD MODELING OF THE RED MUD THICKENER

Abstract

The present paper is devoted to development of CFD thickener model which describes behavior of flocculated suspension in all volume of the unit. Functions for description of hindered sedimentation and effective compressive yield stress of the solid have been additional included to the model by means of built-in and user-added models using UDF. The three-dimensional model takes into account the particle size distribution in the thickener feed pulp and mixing action of rakes. The result of this work is CFD model of the thickener which develops in ANSYS Fluent software package.

Keywords: thickener, red mud, CFD modeling, hindered sedimentation, effective compressive yield stress, particle size distribution (PSD).

Введение

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

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

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

Объект управления

Радиальный одноярусный сгуститель красных шламов диаметром 30 м и общей высотой 3,5 м (цилиндрическая  часть – 2,15 м, коническая – 1,35 м). В центре чана расположена разгрузочная воронка с патрубками для разгрузки сгущенного продукта, по борту чана – кольцевой сливной желоб. В центре чана установлено загрузочное устройство  для приема пульпы из питающей коробки – питающий колодец с поперечной решеткой внутри для улавливания случайных крупных предметов и для гашения скорости потока поступающей пульпы. Гребковая рама сгустителя состоит из двух радиально расположенных граблин с прикрепленными к ним гребками в количестве 20 шт. Рама при помощи крестовин и тяг укреплена на конце центрального вала. Гребки установлены под определенным углом, позволяющим перемещать осадок к центру. Для предотвращения поломок приводного механизма и гребковой рамы при перегрузках вал вместе с граблинами может перемещаться на некоторую высоту в вертикальном направлении.

Геометрия рассматриваемого сгустителя представлена на рис.1 ниже:

24-11-2016-09-42-43

Рис.1 – Геометрия сгустителя и питающего колодца

 

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

 

Таблица 1 –  Параметры технологического процесса

№ п/п Параметр Значение параметра Единица измерения параметра
1 Рабочая температура процесса 102
2 Объемная скорость подачи пульпы в питании сгустителя 300 - 430 м3/ч
3 Объемная скорость откачки сгущенного продукта 80 - 100 м3/ч
4 Концентрация твердого в пульпе питания 1,59 – 2,00 % об.долей
5 Концентрация твердого в сгущенном слое 2,83 – 6,00 % об.долей
6 Плотность твердой фазы 3200 кг/м3
7 Плотность жидкой фазы 1240 кг/м3
8 Вязкость суспензии при рабочей температуре 2,10Е-03 Па*с
 

Гранулометрический состав исследуемого образца красного шлама

Преобладающий размер агрегатов представлен на рис. 2.

24-11-2016-09-46-49

Рис. 2 –  Гранулометрический состав исследуемого красного шлама

 

Состав исследуемого шлама полидисперсный, в основном твердое представлено мелкими классами, из которых класс – < 0,005 мм – около 70 %, содержание в гидросмеси частиц класса 0-0,004 мм способствует структурообразованию, а значит, для данных шламов методы осветления должны сопровождаться процессами флокуляции и агрегации частиц.

Распределение размера частиц на входе в сгуститель

Процесс агрегации частиц в процессе флокуляции можно описать с помощью модели баланса популяций (Population balance model), которая описывает динамику числа частиц заданного размера от процессов соударения частиц, слипания частиц, распада частиц и распределения распада частиц. Уравнение баланса популяций было реализовано в программном пакете MATLAB. Полученный в процессе моделирования баланса популяций диапазон флокулируемых частиц использовался как начальный диапазон распределения размера частиц на входе в сгуститель.

Для описания дисперсного состава твердой фракции в питающем потоке использовалась формула  Розина-Раммлера, заложенная в программном пакете ANSYS Fluent.

Скорость стесненного осаждения

Основное уравнение, описывающее осаждение сферических частиц в жидкой среде – уравнение Стокса. Для описания зависимости скорости осаждения Us от содержания твердой фракции в концентрированной суспензии со взвешенными частицы используется уравнение Ричардсона и Заки 1954 года:

 24-11-2016-09-49-01    (1.1)

где ε = (1 - φ) — пористость, c – индекс стесненного осаждения (индекс Ричардсона-Заки), 24-11-2016-09-50-16 – скорость осаждения Стокса, φ - концентрация твердых частиц.

Фактор стесненного осаждения используется для определения относительной скорости между твердой и жидкой фазой, когда концентрация твердой фазы относительно высока (превышает критическую концентрацию) [2].

Скорость Стокса и индекс стесненного осаждения были получены экспериментальным путем и равны 0,0175 м/с и 89,39 соответственно.

Эффективное сжатие твердого вещества

Для концентраций твердых веществ выше критической концентрация или гелевой точки, 24-11-2016-09-51-37, образованная структура частиц передает силы непосредственно между частицами в осадке.

Предел текучести при сжатии описывается следующими выражениями [1]:

24-11-2016-09-52-11  (1.2)

24-11-2016-09-53-17  (1.3)

Значения параметров  и σn были взяты из опубликованной ранее литературы и равняются 2 Па и 6,5 соответственно. Значение критической концентрации или гелевой точки получено экспериментальным путем и равно 0,0297 об. долей.

 Метод  Эйлера-Лагранжа

На сегодня существует два подхода расчета многофазных потоков: метод Эйлера-Лагранжа и метод Эйлера-Эйлера.

Модель дискретной фазы (DPM) по методу  Эйлера-Лагранжа строится путем решения усредненных по времени уравнений Навье-Стокса для жидкой фазы, которую рассматривают, как континиум, в то время как дисперсная фаза описывается путем отслеживания траектории большого числа частиц через расчетные поля течения. Дисперсная фаза может обмениваться с жидкой фазой импульсами, массой и энергией [3].

В данной модели сделано допущение: дисперсная фаза имеет малую объемную долю (менее 10-12 %), несмотря на приемлемую большую массу частиц  (m частицы ≥ m жидкость). Модель дискретной фазы является единственной многофазной моделью, которая позволяет определить распределение частиц в питании. Расчет траектории частиц происходит в определенном интервале времени.

Поведение частиц в турбулентных потоках

Дисперсию частиц за счет турбулентности в жидкой фазе можно предсказать с помощью стохастической модели отслеживания или облачной модели частиц. Модель стохастического отслеживания (случайное блуждание) включает влияние мгновенных турбулентных пульсаций скорости на траектории частиц за счет использования стохастических методов. Облачная модель частиц отслеживает статистическую эволюцию облака частиц вокруг средней траектории. Концентрация частиц в облаке представлена гауссовой вероятностной плотностью распределения вокруг средней траектории. В обеих моделях частицы не имеют никакого прямого влияния на процесс генерации или диссипации турбулентности в непрерывной фазе.

Уравнения движения частиц

Баланс сил, действующих на частицу, прогнозирует траекторию частиц дискретной фазы путем интегрирования баланса сил на частицу, который записывается в лагранжевой системе отсчета. Этот баланс силы уравнивает инерцию с силами, действующими на частицу, и может быть записан (для направления в декартовой системе координат) как [3]:

24-11-2016-09-55-33  (1.4)

где 24-11-2016-09-59-04 - это сила сопротивления на единицу массы частицы [3].

24-11-2016-09-59-36  (1.5)

Где U – скорость жидкой фазы, Up – скорость частицы,  – вязкость жидкости, ρ – плотность жидкости, ρp – плотность частицы, dp – диаметр частицы, Re – число Рейнольдса, которое определяется как [3]:

24-11-2016-10-01-40  (1.6)

Коэффициент сопротивления, CD, определяется как [3]: 24-11-2016-10-02-28

где a1, a2, a3 –  константы, применяемые для гладких сферических частиц в рамках нескольких диапазонов Re [3].

Также коэффициент сопротивления по Хайдеру и Левеншпилю можно выразить как [3]:

24-11-2016-10-03-55

где α – коэффициент формы, s – это площадь поверхности сферы, имеющей такой же объем, как и частица, и является фактической площадью поверхности частицы [3].

В качестве дополнительных сил, действующих на частицу, 24-11-2016-10-05-31, могут выступать силы, действующие на частицы, которые возникают из-за вращения системы отсчета. Эти силы возникают при моделировании потоков во вращающихся системах отсчет: градиент давления, термофоретическая, броуновское движение, сила Саффмана, например.

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

Результаты

Для построения геометрии использовался стандартный модуль Design modeler. Модель построена из двух основных частей – стационарный объем (корпус сгустителя) и динамический объем (призматический элемент – граблины) (рис. 3).

24-11-2016-10-06-35

Рис. 3 –  Геометрия сгустителя

 

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

Для получения качественной расчетной сетки импортировалась 1/4 часть модели. Для построения использовался модуль Mesh (рис. 5). Области перехода плоскостей и пересечений выполнены с меньшим шагом сетки. На рис. 4 указаны симметричные области, вход и выход потока.

24-11-2016-10-07-28

Рис. 4 – Вход, выход потока и симметричные области

24-11-2016-10-08-25

Рис. 5 – Расчетная сетка объекта управления

 

Также были выделены области контакта двух объемов (рис. 6). Данная область контакта позволит объединить два объема во время вращения одного из объемов.

24-11-2016-10-09-17

Рис. 6 – Область контакта двух объемов – основного корпуса сгустителя и вращающихся граблин

 

В ANSYS Fluent:

  • заданы граничные условия;
  • включена сила гравитации при расчете движения масс частиц;
  • для входа частиц использовалась поверхность inlet;
  • включена инжекция твердой фракции заданного диапазона распределения частиц с помощью модели DPM (Discrete phase model) по закону Розина-Раммлера;
  • заданы свойства рабочей среды и твердой фракции, описанные в работе ранее;
  • для отображения результатов отображаем привязку по ID с поверхностью выхода;
  • анализируем полученные результаты.

Выводы

Результатом работы является трехмерная модель сгустителя (рис. 7), по результатам которой можно оценить концентрацию твердой фракции в каждой точке исследуемого объекта управления, что невозможно сделать на реальном объекте.

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

24-11-2016-10-10-18

Рис. 7 –  Распределение твердой фракции в сгустителе

 

Для построения такой модели в Fluent требуется дополнительно описать параметры стесненного осаждения и уплотнения осадка с помощью UDF (User Defined Functions).

Модель может использоваться для оценки:

  • поведения сгустителя;
  • распределения твердой фракции в объеме аппарата и в верхнем осветленном сливе, что является одним из целевых параметров при управлении процессом сгущения красных шламов;
  • влияния диапазона распределения твердой фракции в питании сгустителя или диапазона флоккул/кластеров, образовавшихся в процессе флокуляции в питающем стакане.

Модель требует доработки. Для дальнейшей отладки работы модели необходимо описать процесс флокуляции в питающем стакане непосредственно  в среде Fluent.

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

  1. Fernando Concha A. Solid–Liquid Separation in the Mining Industry// Fluid Mechanics and Its Applications. Vol. 105. P. 209.
  2. Usher S., Scales P. Steady state thickener modelling from the compressive yield stress and hindered settling function // Chemical Engineering Journal. 2005. No. 111. P. 253-261.
  3. ANSYS, Inc. ANSYS FLUENT. User's Guide. 2003.