CFD-моделирование сгустителя красных шламов

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

Рубрика Физика и энергетика
Вид статья
Язык русский
Дата добавления 30.04.2018
Размер файла 1,6 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

CFD-моделирование сгустителя красных шламов

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

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

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

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

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

Рис. 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.

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

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

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

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

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

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

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

(1.1)

где е = (1 - ц) - пористость, c - индекс стесненного осаждения (индекс Ричардсона-Заки), - скорость осаждения Стокса, ц - концентрация твердых частиц.

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

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

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

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

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

(1.2)

(1.3)

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

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

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

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

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

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

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

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

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

(1.4)

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

(1.5)

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

(1.6)

Коэффициент сопротивления, CD, определяется как [3]:

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

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

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

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

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

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

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

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

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

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

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

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

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

В ANSYS Fluent:

· заданы граничные условия;

· включена сила гравитации при расчете движения масс частиц;

· для входа частиц использовалась поверхность inlet;

· включена инжекция твердой фракции заданного диапазона распределения частиц с помощью модели DPM (Discrete phase model) по закону Розина-Раммлера;

· заданы свойства рабочей среды и твердой фракции, описанные в работе ранее;

· для отображения результатов отображаем привязку по ID с поверхностью выхода;

· анализируем полученные результаты.

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

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

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

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

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

· поведения сгустителя;

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

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

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

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

пользователь пульпа сгуститель шлам

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.

Размещено на Allbest.ru

...

Подобные документы

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

    презентация [2,1 M], добавлен 24.10.2013

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

    презентация [797,7 K], добавлен 15.02.2015

  • Определение расчетной мощности электродвигателя. Выбор материалов червяка и червячного колеса. Определение допускаемых напряжений изгиба. Выбор коэффициента диаметра червяка. Уточнение передаточного числа. Расчет клиноременной передачи, ведущего шкива.

    курсовая работа [1,0 M], добавлен 18.07.2014

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

    курсовая работа [502,5 K], добавлен 07.12.2015

  • Деформация как изменение взаимного положения частиц тела, связанное с их перемещением относительно друг друга, ее причины и механизмы. Виды: растяжение, сжатие, кручение, изгиб и сдвиг. Основные факторы, влияющие на жесткость и прочность твердого тела.

    презентация [1,3 M], добавлен 26.01.2014

  • Дуализм в оптических явлениях. Недостатки теории Бора. Дифракция частиц, рассеяние микрочастиц (электронов, нейтронов, атомов) кристаллами или молекулами жидкостей и газов. Опыты по дифракции частиц. Корпускулярно-волновой дуализм микрочастиц вещества.

    презентация [4,8 M], добавлен 07.03.2016

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

    лабораторная работа [275,9 K], добавлен 29.08.2015

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

    статья [794,6 K], добавлен 07.02.2014

  • Броуновское движение как беспорядочное движение микроскопических видимых, взвешенных в жидкости или газе частиц твердого вещества. Формула Эйнштейна, ее справедливость. Причина броуновского движения, его особенности, хаотичность и интенсивность.

    презентация [932,4 K], добавлен 14.01.2015

  • Явление рассеяния света. Воздействие частиц вещества на световые волны. Понятие рэлеевского рассеяния и частицы пигмента. Относительный показатель преломления частиц и среды. Увеличение количества отраженного белого света. Исчезновение насыщения цвета.

    презентация [361,6 K], добавлен 26.10.2013

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

    дипломная работа [2,2 M], добавлен 21.01.2016

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

    презентация [127,2 K], добавлен 29.09.2013

  • Основы динамики вращения твёрдого тела относительно неподвижной и проходящей через него оси, кинетическая энергия его частиц. Сущность теоремы Гюгенса-Штейнера. Расчет и анализ результатов зависимости момента инерции шара и диска от массы и радиуса.

    курсовая работа [213,6 K], добавлен 02.05.2012

  • Математическая модель и решение задачи очистки технических жидкостей от твердых частиц в роторной круговой центрифуге. Система дифференциальных уравнений, описывающих моделирование процесса движения твердой частицы. Физические характеристики жидкости.

    презентация [139,6 K], добавлен 18.10.2015

  • Явление дифракции частиц. Структурные и магнитные характеристики вещества. Разложение волн по их частотному спектру. Свободное движение частицы. Волновой вектор монохроматической волны. Применение дифракции частиц для изучения физических объектов.

    реферат [109,6 K], добавлен 21.12.2016

  • Изучение корпускулярной концепции описания природы, сущность которой в том, что все вещества состоят из молекул - минимальных частиц вещества, сохраняющих его химические свойства. Анализ молекулярно-кинетической теории газа. Законы для идеальных газов.

    контрольная работа [112,2 K], добавлен 19.10.2010

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

    презентация [723,6 K], добавлен 02.10.2013

  • Коэффициенты диффузии, ступенчатые поверхности. Алгоритм Метраполиса, метод Монте-Карло, парциальное и среднее покрытие, термодинамический фактор. Диффузия системы взаимодействующих частиц. Зависимость среднего покрытия от химического потенциала.

    дипломная работа [1,2 M], добавлен 10.12.2013

  • Основной закон динамики вращательного движения твердого тела относительно неподвижной оси. Изучение методических рекомендаций по решению задач. Определение момента инерции системы, относительно оси, перпендикулярной стержню, проходящей через центр масс.

    реферат [577,9 K], добавлен 24.12.2010

  • Классификация, основные характеристики и методы разделения неоднородных систем. Их роль в химической технологии. Основные параметры процесса разделения жидких неоднородных систем. Осаждение в поле действия сил тяжести и под действием центробежных сил.

    контрольная работа [404,8 K], добавлен 23.06.2011

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.