Однофазные и многофазные математические модели электролиза алюминия тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Анпилов, Сергей Валерьевич

• Анпилов, Сергей Валерьевич
• кандидат физико-математических наук
• 2011

Оглавление диссертации кандидат физико-математических наук Анпилов, Сергей Валерьевич

Введение.

Постановка задачи.

Обзор литературы.

Обзор работы.

Основные результаты.

Глава 1.Математическая модель МГД-процессов в алюминиевом электролизере.

§ 1.1 Рабочая гипотеза модели модели.

§ 1.2 Дифференциальные уравнения масс и импульсов для каждой фазы.

§ 1.3 Математическое моделирование давления.

§ 1.4 Математическое моделирование токов.

§ 1.5 Математическое моделирование электромагнитного поля.

§ 1.6 Полная математическая модель.

Глава 2. Численный метод решения.

§ 2.1 Разделение на этапы по физическим процессам.

§ 2.2 Численный метод решения задачи по этапам.

§ 2.3 Численное исследование разностного метода.

Глава 3. Анализ численных экспериментов.

§ 3.1 Тестовые и верификационные расчёты.

§ 3.2 Моделирование некоторых процессов, возникающих при промышленном производстве алюминия.

Моделирование замены выгоревших анодных блоков.

Оценка вклада индуцированных токов.

Численное исследование влияния формы настыли на мгд-стабильность ванны.

§ 3.3 Сравнение численных экспериментов, проведенных по однофазной и двухфазной математическим моделям.

Введение диссертации (часть автореферата) на тему «Однофазные и многофазные математические модели электролиза алюминия»

Постановка задачи

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

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

- настыль. Она служит защитой от коррозии внутренней поверхности ванны электролитом. Нижняя поверхность ванны — неровная, неровности обусловлены образованием отложений по бокам и на дне ванны (т.н. гарнисаж). Верхняя поверхность жидкости ограничена погруженными в электролит анодами прямоугольной формы. Геометрические размеры рассматриваемого электролизёра составляют 8.9 м в длину, 3.7 м в ширину и

0.65 м в высоту. Рабочая температура колеблется в пределах 940-980°С, при таких температурах плотность алюминия (2.546-103 кг/м3) превосходит плотность криолита (2.112-103 кг/мг), что обуславливает их разделение на два слоя. На рис. 1 изображен схематичный срез алюминиевого электролизёра с погруженными анодами. гарннсаж аноды криолит свободная поверхность катод

Рис. 1. Схематичное изображение электролизёра -поперечный разрез.

Основной процесс, происходящий на катоде - восстановление трехвалентного алюминия А1 + Зе = А1, на угольном аноде происходит разряд кислородосодержащих ионов с образованием окислов углерода СО и СО2. Суммарная реакция, происходящая в электролизере, может быть описана уравнением:

А1203 + ХС = 2А1 + (2Х-3)СО + (3-Х)С02

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

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

Основным показателем работы электролизёра является выход по току -отношение количества практически полученного алюминия к количеству, которое должно выделиться на катоде согласно закону Фарадея. Считается, что первопричина снижения выхода по току - окисление растворенного алюминия анодными газами:

Основные потери энергии происходят в электролите, что* побуждает уменьшать междуполюсное расстояние (МПР). При этом на поверхности раздела металл-электролит легко возникают нежелательные возмущения или волны. При некоторых условиях наблюдается рост амплитуд этих, волн, который называют неустойчивостью, или магнитогидродинамической (МГД) нестабильностью. На рис.2 приведено схематичное изображение участка рабочей зоны электролизёра. Даже простое незатухающее колебание поверхности раздела должно быть устранено, так как при нем возрастает массоперенос растворенного алюминия от катода в межполюсное пространство, где он снова окисляется. МГД-нестабильность является основным препятствием увеличения выхода по току, и её устранение является одним из самых тонких мест в технологическом процессе производства первичного алюминия.

Рис.2. Схематичное изображение рабочей зоны.

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

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

Технология электролиза с обожженными анодами подразумевает периодическую замену анодов, в* процессе чего могут происходить нежелательные возмущения поверхности раздела, приводящие к снижению эффективности производства. При, замене выгоревших анодов, которые выводят парами по мере их выгорания, происходит перераспределение токов (суммарная величина которых порядка нескольких, сотен кА) по анодам, это вызывает перераспределение электрического тока и, соответственно, электромагнитных сил в расплаве. Это, в свою очередь, вызывает изменение распределения скоростей, в электролите и металле, образование вихрейг в средах, что вызывает возмущение' границы раздела* сред металла* и электролита. Часто* при этом возникает МГД-нестабильность, зазора между анодом и жидким^ металлом! становится) меньше допустимого, «3 см (значение1 минимального допустимого МПР р для каждой электролизной ¡ванны индивидуально).

За последние два десятилетия появилось много, работ, посвященных исследованию механизмов, вызывающих волнение на поверхности металла. Само по себе существование различных подходов, дающих объяснения,этому явлению, свидетельствует о том, что исчерпывающего ответа* нет. Математические модели имеют ограниченные области применения. Создание и использование их зависит от поставленной задачи.

Важным параметром, влияющим на эффективность, работы электролизера, является межполюсное расстояние (МПР). МНР — это расстояние между подошвой анодов, и поверхностью металла. Чем больше межполюсное расстояние, тем больше выход по току. Если межполюсное расстояние больше 6 см, то выход по току растет не значительно. Если межполюсное расстояние меньше 3 см, то выход по току резко падает, так как алюминий попадает в реакционную зону, где происходит его обратное окисление. Поэтому для большинства ванн оптимальное МНР лежит в диапазоне 4—6 см.

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

Таким образом, при электролизе криолито-глиноземных расплавов на выход по току влияет:

• Температура рабочей зоны

• Межполюсное расстояние

• Плотность тока

• Состав электролита

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

В [45] подробно описаны неполадки, которые могут возникать в технологическом процессе производства первичного алюминия. К основным неполадкам в работе ванн авторы относят следующие:

• Холодный ход ванны

• Выпадение глинозема на подину ванны

• Затяжной, «негаснущий» анодный эффект

• Горячий ход ванны

Холодный ход ванны. наблюдается если ванна получает недостаточное количество тепла, его причинами может быть пониженная сила тока серии, малое МПР, а также большое количество,1 металла в ванне, опасность холодного хода ванны заключается в увеличении криолито-глинозёмной смеси при понижении; температуры и как результат. — выпадения; её: на дно ванны. Выпадение глиноземаг на подину ванны, или «отравление ванны глиноземом» также. может происходить, если в ванну загружено глинозема больше того количества, которое может раствориться в электролите. :

Анодный эффект возникает при снижении концентрации глинозема в электролите: до 1-2%, поэтому его гасят внесением, глинозема; вI ванну. Елавный признак анодного эффекта — резкое; почти мгновенное, повышение напряжения в ванне. Основной: причиной: анодного, эффекта является ухудшение смачиваемости поверхности: анода: электролитом из-за увеличения? поверхностного * натяжения на границе: раздела; электролит—анод, т.е. на поверхности анодашоявляется! сплошная прослойка из пузырьков газа продукта электролиза. Растет сопротивление на границе электролит-анод, что влечет за- собой потери тока; резкое повышение напряжения и резкое перераспределение токов; по анодам. Для каждой ванны существует критическая плотность тока. Чем она выше, тем реже*, возникает анодный эффект. Величина критической плотности тока зависит главным образом? от состава электролита и температуры. Затяжной анодный эффект характеризуется тем что он длится несколько часов; и его не- удается ликвидировать обычным приемом — загрузкой в электролит очередной порции глинозема. Если затяжной анодный эффект не прекратить, то неизбежен прорыв электролита и металла и даже полный выход ванны из строя.

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

• Чрезмерно малое МПР (избыточное тепло выделяется в рабочей зоне при обратном окислении алюминия);

• Сильно возросшее напряжение из-за увеличения МПР;

• Недостаточное количество алюминия в ванне;

• Анодный эффект;

• Перегрев электролита на каком-то участке (местный перегрев) при неравномерном распределении тока. Главная причина — сильный перепад поверхности раздела металл—электролит. При этом плотность тока растёт и электролит перегревается, происходит усиленное окисление алюминия и повышается температура ванны.

Сильный перепад МНР возникает, когда:

• Нарушена параллельность между нижней поверхностью анодов и дном ванны из-за перекоса анодов после выливки алюминия и перетяжки рамы анодов.

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

• Анод замыкается с алюминием через обломившийся кусок анода или кусок подового блока, который не вынули из ванны.

• Анод касается настыли.

Поскольку физические эксперименты по оптимизации работы электролизера затруднительны, а замеры полей при температуре близкой к 1000 °С на практике дают погрешность порядка 15%, возрастает роль проведения численного эксперимента на основе математического моделирования. Таким образом, математическое моделирование фактически является единственным способом исследования, и визуализации динамических процессов на границе раздела сред электролит-алюминий. В силу высокой температуры и химической агрессивности среды проведение экспериментальных замеров основных характеристик процесса электролиза алюминия в промышленной электролизной ванне крайне затруднено. Поэтому достаточно адекватное математическое моделирование позволяет получить достоверную информацию о форме границы раздела сред, конфигурации распределения скоростей металла и электролит, электромагнитных полей и электрических токов в средах, а также зависимости этих величин от геометрии ванны, конфигурации анодов и условий проведения процесса электролиза. Это даёт возможность выработать рекомендации по оптимальному режиму проведения процесса электролиза и выбору формы рабочего пространства ванны с целью увеличения выхода алюминия по току.

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

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

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

Обзор литературы

В математическом моделировании алюминиевых электролизеров можно выделить два крупных направления:

1. Получение основных полей (магнитного, электрического), формы стационарной поверхности металла и скоростей в установившемся стационарном течении [18, 21, 31, 33].

2. Изучение колебаний поверхности металла и получение критериев устойчивости [8-13, 30].

Как правило, исследование устойчивости проводится на основании приближения «мелкой воды» и линеаризации в окрестности стационарного решения с последующей оценкой вкладов членов в уравнении движения. Полученные уравнения отличаются правой частью, авторы работ по-разному учитывают возмущающие электромагнитные силы, также уравнения отличаются граничными условиями [8, 9].

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

Первый из них [21], основан на двумерных уравнениях Навье-Стокса< и к-е модели турбулентности; учитывается только движение в горизонтальной плоскости, не учитывается вертикальный перенос импульса, т. е. трение относительно узких слоев жидкости о дно ванны, нижнюю поверхность анода и между собой. Для нахождения скоростей, давления и формы поверхности металла решается усредненная по времени система уравнений Навье-Стокса. В [21] приведены результаты расчета скоростей в металле и электролите и формы поверхности металла для ванн с различными схемами ошиновки. В работе [21 ] также описана модель для расчета электрического и магнитного поля в алюминиевом электролизере. Для нахождения распределения тока в электролизере авторы работы используют закон Ома. Электрический потенциал находится из уравнения Лапласа для потенциала. Авторы работы [21] делают предположение, что величина компоненты электрического тока, направленной вдоль длинной стороны электролизера мала по сравнению с вертикальной компонентой и с компонентой, направленной вдоль короткой стороны электролизера. Уравнение на электрический потенциал решается в плоскости, в поперечном сечении электролизера. Магнитное поле находится с помощью закона Био-Савара-Лапласа, интегрирование ведется по объему, занятому токонесущими элементами системы. Магнитное поле рассчитывается без учета ферромагнетиков.

Второй подход — это рассмотрение движения металла, электролита, распределения температуры в вертикальном разрезе, представляющий'собой систему уравнений теплоэлектропереноса, Максвелла и Навье-Стокса [31]. Этот подход не позволяет получить распределение температур и* линии циркуляции в планарной плоскости, тем не менее, конвекция учитывается с помощью эффективного коэффициента теплопроводности.

Третья, наиболее распространенная в настоящее время модель Моро-Эванса [33]. Одна из главных идей в формулировке этой модели состоит в том, что основной интерес представляет описание циркуляции средних слоев жидкостей. Как следствие этого в уравнении движения в горизонтальной плоскости учитывается трение слоев друг о друга. Следующий шаг - это оценка значений вклада каждого члена в уравнение для моментов движения. В результате разность между электромагнитным силами и градиентом давления уравновешивается силами трения, которые предполагаются пропорциональными скорости. Первая двумерная модель в приближении «мелкой воды» была предложена в работе [3], в данной модели используется нелинейное волновое уравнение для описания поверхности невязкой жидкости. Нелинейность обусловлена электромагнитной силой. К сожалению, в работе не приводится вывод этого уравнения.

В работе [33] рассматривается электролизер бесконечных размеров. В основе модели лежат уравнения Навье-Стокса. Используется классическая техника.линейной теории гидродинамической устойчивости, за исключением того, что для моделирования течений на середине каждого слоя предполагается линейный закон трения в приэлектродной области. Авторы отмечают, что электромагнитная неустойчивость может развиваться, если плотность, тока в жидком алюминии достаточно велика. Она может порождать крупномасштабные волны (1,5-2 м); распространяющиеся в направлении горизонтального тока. Возрастание уровня»турбулентности (или возрастание коэффициента трения) может стабилизировать возмущения. Это предполагает существование некоторого механизма насыщения, и требует развитие нелинейных моделей, способных различать условия насыщения от условий развития хаоса. В работе [33] получен результат, который противоречит [4]. А именно: рассматривается электролизер со сплошным анодом. Система^ запитана вертикальным электрическим проводником, расположенным выше анодного блока. Какова бы ни была длина проводника, система (в рамках используемой модели) неустойчива. Если добавить однородное поле в г-направлении, то система становится устойчивой- в случае, когда В: -компонента магнитного поля превышает некоторое критическое значение 10 “3 7л-. Авторы объясняют это различие тем, что-они учли индуцированный электрический ток, в то- время- как в статьях, опубликованных ранее, индуцированным1 током пренебрегают. Уравнения выводятся при следующих предположениях:

• Движения несжимаемых жидкостей подчиняется уравнениям Навье-Стокса.

• Электромагнитное поле подчиняется квазистационарным уравнениям Максвелла.

• Электрические токи вне электролизера, в анодных и катодных блоках известны и не зависят от времени.

• Флуктуации вокруг стационарного движения настолько маленькие, что уравнения можно линеаризовать в окрестности стационарного решения.

• Изменения электромагнитного поля настолько маленькие, что ротор электрического поля равен нулю.

• В линеаризованной системе для флуктуаций можно пренебречь вязкими членами.

• При расчетах дополнительно предполагалось:

• Стационарное поле скоростей - нулевое.

• Стационарная поверхность — плоская.

• Аноды и катод представляют из себя два слоя.

В [6] выводится волновое уравнение методом разложения по малому параметру. В качестве малого параметра взяты: отношение глубины слоев к горизонтальным размерам ванны, отношение амплитуды волны к глубине, отношение проводимости электролита к проводимости металла. Получены уравнения, аналогичные уравнениям в работе [3], усиленные наличием анодного слоя. В работе [12] демонстрируются области справедливости подходов на основе теории «мелкой воды», хорошо описываются неустойчивости, порожденные вертикальным магнитным полем. Авторы строят трехмерное интегро-дифференциальное уравнение для описания колебаний в металле, осредняя его по глубине слоя, получают новое двумерное волновое уравнение. Тщательно моделируются аноды, приведена трактовка граничных условий. Из-за этого коэффициент в волновом уравнении отличается от [26]. Новое волновое уравнение было выведено из теории «мелкой воды» в работе [27] и в других работах этих же авторов. В этом уравнении сила Лоренца выражена непосредственно в терминах движения жидкости. Приводится общий энергетический критерий для того, чтобы установить, какие движения могут быть неустойчивыми. Это могут быть как бегущие, так и стоячие волны. Приводится достаточное условие стабильности, которое не требует решения основных уравнений.

Позднее эта модель получила свое развитие в работе [28], где задавалось идеализированное распределение магнитного поля и электрического тока в жидкой среде, что позволило получить решение гидродинамических уравнений в виде аналитических выражений. В" работе [28] сделано предположение, что распределение электромагнитных величин при продольном расположении электролизеров в серии двухмерно и не зависит от координаты х вдоль, длинной* стороны ванны, аналогичное предположение для токов сделано в работе [21]. Для4 расчета магнитного поля« от соседнего’ ряда электролизеров,* соседний* ряд представляется- как линия с током, равным по модулю току серии. Магнитное поле от соседнего ряда электролизеров и от ошиновки считается* по ^закону Био-Савара-Лапласа для прямого отрезка с током*. Для нахождения магнитного, поля* от токов; текущих внутри электролизера; из всех возможных решений уравнения

Ч2В = 0 выбираются те, которые соответствую простейшему распределению плотности тока в жидкости - как линейной функции координат. Индуцированные токи в модели авторами [28] не учитываются. В работе [29] авторы описывают методику расчета распределения тока и магнитных полей для конкретной конструкции электролизера. При этом авторы не приводят данных об идентификации модели. Далее в- ряде работ, например [20, 21] также используется модель Моро-Эванса, для расчета магнитного поля используются замеры магнитного поля по периметру ванны. В этих работах, как и в модели Бояревича [9], также делается переход от трехмерной задачи к двухмерной, он делается путем усреднения трехмерных уравнений (Др = о ) по высоте каждого слоя. Получая в результате усредненный по высоте каждого слоя потенциал, как функцию от горизонтальных координат, можно найти горизонтальные токи. Вертикальный же ток находится с помощью предположения, что трехмерный потенциал равномерно изменяется по вертикали в пределах каждого слоя. Влияние индуцированных токов не учитывается на этапе расчета электрического потенциала, а учитывается после расчета скоростей и магнитного поля при пересчете электрического тока. Магнитное поле от токов, текущих внутри ванны, считается путем решения дифференциальных уравнений с граничными условиями -экспериментальными замерами. Наблюдается совпадение рассчитанных траекторий движения металла с экспериментальными, полученными в [27].

В работе [34] авторы исследуют циркуляцию электролита и металла в алюминиевом электролизере. Для описания совместного движения металла и электролита в работе применяется система, уравнений Навье-Стокса. Значения электромагнитных сил задаются аналитически по формулам, не зависящим от вертикальной координаты. При таком задании сил вертикальная компонента скорости в слое металла близка к нулю, что свидетельствует о горизонтальной конвекции алюминия. С этим связан тот факт, что отклонение положения поверхности металла от начального 5 задавалось прямой) в работе составляет величину ~10" м.

В работе [35] как и в работе [34] авторы исследуют топологию движения расплава в алюминиевом электролизере. Однако в работе [35] используется другой подход. В [35] предлагается концепция оценки гидродинамики электролизера, позволяющая получить информацию о характере движения расплава. В этой работе говорится, что поле скорости напрямую связано с ротором электромагнитной силы. Если расчетным или экспериментальным путем найти распределение плотности тока и индукцию магнитного поля в расплаве, то можно качественно определить топологию течения. При этом имеется в виду, что траектории движущихся объемов расплава должны быть замкнуты вследствие сИу V = 0. Предлагаемая концепция позволяет без расчета гидродинамики расплава составить представление о характере его движения при изучении схемы ошиновки, что дает технологам возможность оценить последствия предлагаемых решений. формулы векторного анализа, можно получить выражение для ротора электромагнитной силы из которого следует, что завихренность силы возникает в том случае, если плотность тока меняется вдоль направления магнитного поля или магнитное поле меняется вдоль направления плотности тока. Так как толщина слоя электролита и жидкого алюминия много меньше горизонтальных размеров ванны, то основная циркуляция будет происходить в горизонтальной плоскости. Ответственной за эту циркуляцию является вертикальная составляющая ротора силы, из чего следует, что завихренность электромагнитной силы, а значит и горизонтальная циркуляция расплава, вызываются неравномерностями в распределении вертикальной компоненты плотности тока в объеме расплава. Таким образом, в [35] говорится, что, если ликвидировать неравномерность вертикальной компоненты плотности^ тока, то в любом внешнем магнитном поле, т.е. при любой схеме ошиновки, расплав будет оставаться практически неподвижным.

Согласно оценкам [35] значение плотности индуцированных токов в электролите пренебрежимо мало по сравнению с общим током. В металле для электролизера с током серии 175 кА величина индуцированных токов сравнима с общим током.

Используя уравнения электродинамики сИу у = 0, сИу В = 0 и известные

В работе [35] описаны различные модификации ошиновки электролизера для оптимизации распределения магнитного поля и увеличения выхода по току. Отмечается, что важнейшей мерой повышения эффективности работы, электролизера является уменьшение скоростей движения расплава в межэлектродном пространстве. Основным способом уменьшить скорости циркуляции расплава является снижение электродинамических сил уменьшением индукции магнитного поля и распределение ее так, чтобы вместо больших контуров циркуляции расплава образовалось несколько с малыми скоростями и перекосами уровня. Для компенсации' магнитного поля тока соседнего ряда электролизеров, расположенных в. два' ряда применяется ошиновка, создающая встречный компенсирующий поток. Состояние настылей на подине определяет величину и направление горизонтальных токов. Если край настыли у борта за проекцией* края анода, то горизонтальные токи максимальны по величине и протекают в сторону борта. Если настыли заходят под анод, то горизонтальные токи протекают к центру. Когда граница круто падающей настыли находится под. краем анода [1][35], то горизонтальные токи минимальны. При' этом металл как бы зажимается- в пространстве, ограниченном проекцией, анода, увеличивается гидравлическое сопротивление, обеспечивается более спокойное состояние расплава.

Тем не менее, никакие модели, не могут претендовать на точное описание процессов в электролизере. Даже использование пакетов программ известных фирм [11][42] показывает значительное расхождение в вычисленных и измеренных значениях г-компоненты магнитного поля. Задача состоит в том, чтобы на основании модели можно было судить о том, как изменения конструкции и технологических параметров отражаются на количественных и качественных характеристиках процессов, происходящих в электролизере.

Количественной характеристикой неустойчивости в алюминиевом электролизере является технологический параметр «уровень шума», который представляет собой амплитуду флуктуации напряжения на ванне и измеряется, в вольтах. Существуют анодные (высокочастотные) шумы и катодные (низкочастотные) шумы. Появление анодных шумов происходит в результате изменений в газовом слое под рабочей поверхностью анода. Поверхность анода подвергается воздействию выделяющихся из-под него газов, которые содержат СО и СО2, а также 02, содержащемся в воздухе. Эти газы, вступая во взаимодействие с углеродом вертикальных граней анода, приводят к их повышенному износу. Поэтому периферия анода обгорает быстрее, и расстояния между концами штыря и подошвой анода на периферии становится меньше, чем в центральной зоне, что может привести к сильному перепаду поверхности раздела металл-электролит. Выделяющиеся пузырьки газов влияют на распределение токов в аноде, так как налипают на нижнюю поверхность и понижают коэффициент электропроводности в этой области. Так же эти пузырьки вызывают мелкую рябь на поверхности металла. Появление катодных шумов происходит в результате колебаний поверхности, металла, обусловленных гравитационными и электромагнитными силами.

Обзор работы

• В главе 1 предложена трёхмерная двухфазная модель для описания магнитогидродинамических процессов, протекающих в алюминиевом электролизёре.

• В главе 2 описан численный метод решения предложенной системы уравнений.

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

Основные результаты

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

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

• Численно исследована МГД-стабильность ванны при замене различных пар выгоревших анодов и зависимость МГД-стабильности работы электролизёра от формы рабочего пространства ванны для заданной конфигурации анодов.

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

Особую благодарность автор выражает академику РАЕН, доктору физико-математических наук, профессору физического факультета МГУ им. М.В. Ломоносова Р.Н. Кузьмину за постоянное внимание, консультации по физической постановке задачи и полезные советы в ходе выполнения диссертационной работы.

Автор благодарен преподавателям кафедры вычислительных методов МГУ им. М.В. Ломоносова профессору А.П. Фаворскому, профессору Г.Г. Еленину, профессору A.B. Гулину, профессору В.Б. Андрееву, профессору С.И. Мухину, профессору Н.В. Соснину, доценту Терновскому В.В. за обсуждение работы и ценные замечания.

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

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Анпилов, Сергей Валерьевич

Основные результаты работы:

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

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

• Численно исследована МГД-стабильность ванны при замене различных пар выгоревших анодов и зависимость МГД-стабильности работы электролизёра от формы рабочего пространства ванны для заданной конфигурации анодов.

Заключение

Список литературы диссертационного исследования кандидат физико-математических наук Анпилов, Сергей Валерьевич, 2011 год

1. Овчинников В.В., Проворова О.Г., Пингин В.В., Пискажева Т.В. Математические модели и МГД-явления в электролизере Содерберга // Цветные металлы. 1997. № 1. С. 61-63.

2. Н. П. Савенкова, А. В. Шобухов, С. В. Анпилов, Р. Н. Кузьмин «Математическое моделирование физико-технологического процесса электролиза» // Прикладная физика. 2009. №6. С.43-51

3. Urata N., Mori К., Ikeuchi Н. Behavior of Bath and Molten Metal in Aluminium Electrolytic Cell // Light Metals Japan. 1976. P. 573-583.

4. Sele T. Instabilities of the Metal Surface in Electrolytic Cells // Light Metals. 1977. P. 7-24.

5. Moreau R., Ewans J.W. An analysis of hydrodynamics of aluminum reduction cells // Jour. Electrochem. Soc. 1984. Vol. 131, No.10. P. 22512259.

6. Sneyd A.D., Wang A. Interfacial Instabilities due to MHD Mode Couplingin Aluminium Reduction Cells // Jour. Fluid Mech. Д994. Vol.263. P.343-359. '

7. Bojarevics V.V., Romerio M.V. Long waves instability of liquid metal-electrolyte interface in aluminium electrolysis cells: a generalization of Sele’s criterion //Eur. Jour. Mech. В /Fluids. 1994. Vol.13, No.l. P. 33-56.

8. L. Leboucher, V. Bojarevics & K. Pericleous The shallow water approximation applied to the aluminium electrolysis process // Centre for Numerical Modelling and Process Analysis University of Greenwich, 1998.

9. M. Dupuis, V. Bojarevics and J. Freibergs, “Demonstration Thermo-Electric and MHD Mathematical Models of a 500 kA A1 Electrolysis cell: Part 2”, Light Metals, TMS, (2004), 453-459.

10. М. Dupuis, V. Bojarevics and J. Freibergs, “Demonstration Thermo-Electric and MHD Mathematical Models of a 500 kA A1 Electrolysis cell”, Proceedings of the 42nd Conference on Light Metals, CIM, (2003), 3-20.

11. M. Dupuis and I. Tabsh, “Thermo-electric Coupled Field Analysis of Aluminum Reduction Cells using the ANSYS® Parametric Design Language”, Proceeding of the ANSYS® Fifth International Conference, volume 3, 17.80-17.92, (1991).

12. M. Dupuis and I. Tabsh, “Thermo-electromagnetic Modelling of a Hall

13. Heroult Cell”, Proceeding of the ANSYS® Sixth'International Conference, volume 4, 9.3-9.13, (1994). .

14. Проворова О.Г., Пингин B.B., Овчинников B.B., Пискажева Т.В., Горин Д.А. Математические модели физических полей в электролизере Содерберга//Магнитная гидродинамика. 1998. Т.34. С. 375-385.

15. Самарский А.А., Попов Ю.П. Разностные методы решения задач газовой динамики. — М. : Наука, 1992.

16. Самарский А.А., Гулин А.В. Численные методы математической физики. — М. : Научный мир, 2000.

17. Н. П. Савенкова, А. В. Шобухов, С. В. Анпилов, Р. Н. Кузьмин «Математическое моделирование физико-технологического процесса электролиза» // Прикладная физика. 2009. №6. С.43-51.

18. Головизнин В.М., Самарский А.А. Некоторые свойства разностной схемы «Кабаре» // Математическое моделирование. 1998. Т. 10, № 1. С. 101-116.

19. Головизнин В.М., Самарский А.А. Разностная аппроксимация конвективного переноса с пространственным расщеплением временной переменной // Математическое моделирование, 1998. Т. 10, № 1. С. 86100. '

20. Gerbeau J.-F., Lelievre Т., Le Bris С. Simulations of MHDI flows with moving interfaces // Jour. Comput. Phys. 2003. Vol. 184. P. 163-19k

21. Нигматулин P.И. Основы механики гетерогенных сред. -М.: Наука,1978. :

22. Дорохов ИЛ Г., Кафаров В.В., Нигматулин Р:И. Методы механики сплошной» среды для: описания* многофазных многокомпонентных смесей с химическими реакциями и процессами тепло и массопереноса. Прикладная математика и механика. 1975. Т.39, No.3. С. 485-496:

23. Нигматулин Р.И. Мелкомасштабные течения; • и поверхностные эффекты в: гидродинамике многофазных сред: Прикладная математика и механика. 1971. Т.35, №3, С. 450-463

24. Самарский A.A., Попов Ю.П. Разностные методы' решения задач газовою динамики; 1-М::;Наука, 1980:

25. Gerbeau, J. Metal pad roll instabilities Text.,/ J; Gerbeau,.T. Lelievre, G. Le Bris, N. Ligonesche, C. Vanvoren // Light Metals. — 2002.

26. Grjotheim, K. Introduction to Aluminum Electrolysis |Text. / K. Grjotheim,

27. H. Kvande // Dusseldorf: Aluminium-Verlag. — 1993. Р.260Л

28. Moraru, A. Current field in an aluminum electrolysis cell Text. / A. Moráru, A. Panaitesku, A. Grisu // Light Metals. — 2005. — P.469-474.

29. Lindsay, R.I. Applications of new stability criteria to industrial cell design Text. / R.I: Lindsay, P.A. Davidson // Light Metals. 1997. - P.423-428.

30. Lympany, S.D. The hall-heroult cell: some designs alternatives examined bya mathematical model Text. / S.D. Lympany, J.W. Evans // Metall. Trans.

31. B. 1983. -March. - Vohl4B;-P.63-70. .

32. Dupuis M. Weakly coupled thermo-electric and MHD mathematical models of an aluminium electrolysis cell Text. / M. Dupuis, V. Bojarevics // Light Metals. — 2005. — P.449-454.

33. Deslöux J:, Flueck M., Romerio M.V. Modelling for Instabilities in Hall-Heroult Cells: Mathematical; . and' Newmerical Aspects //

34. Magnetohydrodynamics in Process Metallurgy. 1991. P. 107-110.

35. Василевский, О.И. Численное исследование течений- расплавов в . алюминиевом электролизере Текст. / О.И.' Василевский^ ВЖ Иванов,

36. В:А. Крюковский* С.А; Щербинин // Цветные металлы. 1989: - №9. -. . С.50-54. . ' ' : ■ . , . ■

37. Горбачев; Е.В. Об описаниш МГД явлений Bi алюминиевых электролизерах различной мощности Текст.! / Е.В. Горбачев, Э.В. Щербинин // Цветные металлы. 1990. - №3. — С.47-52.

38. Самарский A.A., Гулин A.B. Численные методы* математической физики. -—М.: Научный мир, 2000.

39. Коростелев. И.Н.,. Патрин Г.С., Пингин В.В., Проворова О.Г., Шайдулин E.Pi’ “Магнитогидродинамические явления в электролизере”

40. Ch. Droste, М. Segatz and'D.Vogelsang; Magnetohydrodynamics Instability Analysys in Reduction Cells, Light Metals^ 1999

41. Исследование магнитогидродинамических явлений в- алюминиевых электролизерах Текст.: отчет о НИР / Красноярск, 1993. — 50с.

42. Алаторцев A.B., Кузьмин Р.Н., Проворова О.Г., Савенкова Н.П. Исследование МГД-нестабильности в алюминиевом электролизере // Сб.тезисов XI международной конференции «Математика. Компьютер. Образование». Дубна, 2004. С.89.

43. Минцис М.Я., Поляков П.В., Сиразутдинов Г.А. «Электрометаллургия алюминия» Новосибирск: Наука. - 2001.

44. Р.И. Нигматулин «Мелкомасштабные течения и поверхностные эффекты в гидромеханике многофазных сред» // Прикладная математика и механика, 1975, том 35. С.451-463

45. С Н Droste PHOENICS' Applications in the Aluminium Smelting Industry // PHOENICS User Conference 2000 Proceedings, C.05/00-22

46. Сушков А.И., Троицкий И. А. Металлургия Алюминия. —1. М. : Металлургия, 1965. ч

47. Содержание диссертации изложено в следующих работах:

48. М.-Ижевск: "Регулярная и хаотическая динамика": сб. тезисов, 2008, том1, С.113.

49. Н. П. Савенкова, А. В. Шобухов, С. В. Анпилов, Р. Н. Кузьмин «Математическое моделирование физико-технологического процесса электролиза» // Прикладная физика; 2009; №6; с.43-51.

50. Н.П. Савенкова, Р.Н. Кузьмин, О.Г. Проворова, A.B. Шобухов, С.В. Анпилов, С.А. Складчиков «Двумерная и трехмерная математические подели электролиза алюминия» // Динамика сложных систем, 2009, №2,1. С.53-59.

51. Савенкова Н.П., Анпилов С.В., Кузьмин Р.Н., Проворова О.Г., Пискажова Т.В. «Двухфазная трёхмерная модель алюминиевого электролизёра» // Прикладная физика, 2011, №6. С.25-33.

52. Савенкова Н.П., Анпилов С.В. «Двухфазная трёхмерная модель МГД-стабильности алюминиевого электролизёра» // Тихоновские чтения, сб. тезисов, М.: МАКС Пресс, 2011, С 70-71.

53. Анпилов С.В. Двухфазная 3-0 модель алюминиевого электролизёра // Препринт, М.: МАКС Пресс, 2011, 18с.

54. Савенкова Н.П., Анпилов С.В., Кузьмин Р.Н., Проворова О.Г., Пискажова Т.В. «Двухфазная 30 модель мгд-явлений алюминиевого электролизёра» // Сборник докладов третьего международного конгресса «Цветные металлы 2011». Красноярск,. - С. 282-286.