Journal of Siberian Federal University. Engineering & Technologies 2 (2010 3) 159-170
УДК 669.295
Методы эффективного управления технологическим процессом электролитического получения алюминия
Т.В. Пискажова*
Сибирский федеральный университет, 660041 Россия, Красноярск, пр. Свободный, 79 1
Received 7.05.2010, received in revised form 28.05.2010, accepted 8.06.2010
На основе математической модели процесса разработаны методы управления технологией электролиза алюминия, принесшие экономический эффект. В статье представлены алгоритмы оптимальной стабилизации химического состава электролита и программа-имитатор процесса для поддержки принятия технологических решений.
Ключевые слова: математическая модель процесса, алгоритмы управления, виртуальный электролизер
Введение
Задача управления технологией электролиза алюминия состоит из двух основных частей:
1. Отыскание набора целевых значений для технологических параметров процесса, способных обеспечить хорошие технико-экономические показатели.
2. Осуществление мероприятий по достижению этих целей.
Первая часть задачи управления осложняется отсутствием единого функционала качества процесса, выраженного от параметров технологии. Себестоимость алюминия, расход энергии и сырья зависят не менее чем от 100 основных переменных электролиза, свойств сырья, конструкции, и не менее половины этих величин взаимозависимы. Поэтому необходим инструмент, позволяющий рассчитывать и отображать изменения технологических параметров при воздействиях на процесс.
Важной составляющей второй части задачи управления является разработка алгоритмов Автоматической системы управления технологическим процессом (АСУТП) для стабилизации технологических параметров вблизи назначенных целей.
Сущность процесса электролиза состоит в следующем. В прямоугольных ваннах (электролизерах), футерованных углеродистыми материалами, находится на подине слой расплавленного алюминия, выше - слой электролита (рис. 1). На стенках ванны образуется настыль из застывшего электролита. Температура процесса около 960 °С. Сверху в ванну опущен уголь-
* Corresponding author E-mail address: Piskazhova@yandex.ru
1 © Siberian Federal University. All rights reserved
Рис. 1. Поперечный разрез алюминиевого электролизера с самообжигающимся анодом
ный анод, частично погруженный в электролит. Катодом служит расплавленный алюминий. Электролитическое разложение глинозема, растворенного в электролите, происходит в междуполюсном зазоре анод - металл.
При прохождении тока через токоподводящие элементы и сам электролизер возникают значительные электромагнитные силы, создающие стационарные магнитогидродинамические (МГД) явления в расплавах внутри ванны, а также нестационарные возмущения поверхности раздела.
При прохождении тока через электролит выделяется большое количество тепла, которое, с одной стороны, необходимо для поддержания реакции электролиза, с другой - температуры различных слоев и областей электролизера должны быть постоянны во времени и оптимальны, обязательно наличие защитного слоя настыли на бортах электролизера.
Из-за различного химического состава расплавленного и застывшего электролита в электролизере присутствует взаимосвязь теплового и массового балансов. Поэтому управление составом электролита эффективно осуществляется только вместе с решением тепловой задачи.
Математическая модель процесса
Таким образом, математическое описание процесса должно содержать следующие разделы.
• Модель МГД-явлений в расплавах: расчет деформации поверхности раздела, поля скоростей расплавов.
• Тепловую динамическую модель электролизера: расчет температур выделенных областей и поверхностей, расчет изменения толщины настыли, уровня спекания анода, тепловых затрат на реакции и нагрев сырья.
• Модели массового и электрического балансов: поступление и расход сырья, производство металла, расчет падения напряжения в электролите, прогноз анодных эффектов.
Модель для управления построена на обыкновенных дифференциальных уравнениях и зависимостях, полученных из предварительных расчетов трехмерными моделями; модель служит для принятия решений в реальном времени.
В рамках одной статьи невозможно рассмотреть все перечисленные выше подмодели, отметим здесь найденное нами решение проблемы оценки МГД-стабильности электролизера для модели управления и приведем некоторые уравнения центральной модели процесса -тепловой.
Нарушение МГД-стабильности, развитие длинноволновых возмущений на поверхности раздела расплавов является серьезной технологической проблемой и возникает при снижении междуполюсного расстояния (МПР) меньше критического значения. При этом слишком большое МПР приводит к перегреву электролизера и повышенному расходу энергии.
Запас МГД-стабильности зависит как от магнитных полей, конструкции электролизера, так и от меняющихся в ходе технологической практики уровня металла, плотности электролита, силы тока, длины настыли на подине. Моделирование волнообразования в электролизерах является пространственной задачей с пересчетом магнитных полей и не может быть использовано в модели для управления.
Обозначим через: Ьх, Ьу - длину и ширину электролизера, 1С - ток серии, рм - плотность металла, рэ - плотность электролита, Нм - высоту металла, Бг - вертикальную компоненту магнитного поля. Волновые вектора
к = (кх,к ) и к' = (к'х,кУ) ,кх =
Lx
, я , п , к y =-n , кх = — m ,
y L y Lx
,, л ,
y = Tn
Lyy
т, п, т, п - целые числа.
На основе критерия Бояревича - Роме рио [1] на ми в работе [2] была получена зависи мое ть критического МПР от те хнологических параметров электролизера в следующем виде:
где
4/)=
мпрК
= 40
2 • Wo
n m - nm
LiL y
( Pm-PS )HM
(/ (m '+m,ri+n,/ )- l(m'-m,ri-n,/ ))-
(1)
n' m + nm'
' LLy
(i (m'+m, n'-n, /) - l(tn'-m, n'+n,/))
gn
ЛГ
Lí LÍ
y
Lí L2„
y
Lí
Lí
Представим коэффициенты Фурье Д»г,п,/) вертикалпной составляющей магнитного поля с учетом длины настыли в виде
, А 4 , А . rn-x . ту
1 ím,n,l i=-J B í x, y, l i sin-sin—— dxdy,
L L j-T L L
x y l x -y
где l - длина настыли. С помощью пространствпнных моделей магнитное поле Bz(x,y,l) было
рассчитано при различных значениях длины настыли. Зависимость коэффициентов Фурье Bz
от длины настыли была аппроксимирована многочленом следоющего вида:
Y13
Коэффициенты многочлена найдены методом наименьших квадратов и для каждого типа электролизера хранятся в виде матрицы в программе.
Для каждой конструкции электролизера трехмерным моделированием были определены неустойчивые пары волновых чисел т, пит', п'.
Таким образом, по формуле(1) возможно определение критического междуполюсного расстояния, или, что то же самое, критического напряжения для электролизера в зависимости от технологических параметров. 1В раббт е [2] при водятся расчеты критических МПР для разных конструкций электролизеров.
Остановимся на методикл тепловлк расчетов. Для описания динамики средних температур металла, электролита, слоев футеровки и анода используются сосредоточенные (нульмер-ные) модели. По токи тепла через границы исследуемых областей опр еделяются по одномерным приближениям с учетом фазовых переходов.
На рис. 2 выделены области для расчета средних температур электролизера с самообжигающимся анодом.
Теплообмен в областях глинозема, корей, электролита, металле, подовых блоков, теплоизоляции описывается уравнениями относительно средних температур:
Здесь Т - средние температуры соответствующих слоев, О - выделение тепла при прохождении тока в электролите, <2ан - теплообмен с анодом, с/ - удельные теплоемкости соответствующих материалов, М/[1 - массы соответствующих областей, QГ - потери тепла на нагрев и растворение глинозема, 2Р - расход тепла на основную реакцию, - тепловой поток от внутреннего источника тепла в катоде. Потоки тепла через границы рассматриваемых слоев находятся по выписанным ниже формулам:
бэ-и = ТЭ-Я • аэ-н(Т.э ~ ) - тепловой по ток из электролита в настыль;
Тшк - температура ликвидуса электролита;
бэ-м =^э-м • кэ-м(Тэ ~тм) - тепловой пото 1С из электролита в металл;
О-кор-г = ^кор-г ■ Ккср-г (Ткср ~тг) -от коржи 1С глинозему на корке;
Ом-к =5м-к ■ км-к(ТМ -Тк) - из металла в катод;
Ом-н = ?>м-н ■ ам-н(Тм ~тшк) - из металла в настыль;
Ог-ат = Яг-ат ■ аг-а»^!^ "Тат) - с поверхности глинозема в атмосферу, 7>|г=0 - температура поверхности глинозема;
Ок-дн = $к-дн' кк-дн(Тк ~Тди) -из подовых блоков в теплоизолирующее дно;
Одн-ат = $дн-ат-Ядн-ат(Гдн\л, -Т™) - поток с наружной поверхности дна в атмосферу.
смМ,
а +а» -б.э-н -е э-м -о. г-аР -еэ.
еэ -м - <2м - н - а
(2)
Анодная часса
Спеченный анод
По!' руженный Глинозем, Корка
анод Электролит гарнис- £ 1 ^ £ о о
Эгектралит сэж г1
Настыль З1
Металл к к Б £ £
Подовые блоки
Теплоизоляция, стальной кожух
Атмосфера
Рис. 2. Схематичное представление электролизера для тепловых расчетов
- площади теплообмена соответствующих областей. Коэффициенты теплоотдачи аГ-ат, адн-ат, аЭ-н, аМ-н считаются заданными, а коэффициенты теплопередачи ккор-Г, кЭ-М, кМ-к, кк-дн вы-
числяются по формулам к1 _] = 2/
Ч+ 4 ^ Л ^ ,
, где Д,, Д,- - толщины соседних областей, к,, к,- - их
теплопроводности.
Подобная система теплообмена выписывается для самообжигающегося анода в работе [3] и рассчитывается высота спеченного анода.
В работе [4] для расчета теплопотерь через бортовую футеровку и настыль приведена система уравнений, подобная системе (2), и с использованием условия фазового перехода получена формула (3) для динамического расчета толщины настыли:
зч:1= 0.5
РнЬн
(
+, 0.25
А-а
РнЦ
Тэ ~ТлЭик)
+ка т
т Vлик н/ • РнЬн
(3)
Здесь рн, Ьн, кн - плотност1>, удельная теплота плавления, теплопроводность настыли, 5н -толщина настыли, Тя - средняя температура слоя настыли, Д/ - расчетный шаг по времени. Формулы для определения толщ ин настыл и в металле и электролите отражают связь массового и теплового баланса электролизера, что необходимо для управления химическим составом электролита.
Начальное стационарное распределение температур находится совместным решением систем теплообмена с нулевой правой частью во всех областях. В работе [4] представлены результаты расчетов различных управляющих воздействий на тепловые параметры электролизера.
Модель растворения и потребления глинозема с учетом свойств сырья описывалась нами в [5]. Используемые в модели процесса уравнения электрического и массового баланса были приведены в [6].
2
Рис. 3. Схема работы динамической модели
Совокупность систем дифференциальных уравнений для расчета средних температур и теплообмена в выделенных областях, а также уравнения массового баланса решаются численно методом Эйлера.
Порядок работы динамической модели в целом представлен на рис. 3.
Программа, содержащая динамическую модель, для расчета получает начальные значения технологических переменных (состав и температура электролита, уровни расплавов, значения силы тока и напряжение на электролизере и т. п.) и по заданному типу электролизера использует параметры конструкции электролизера (размеры шахты и анода, толщины и теплофизиче-ские свойства футеровки, характеристики магнитного поля). Осуществляется расчет начального состояния, затем при воздействиях рассчитываются динамические изменения параметров. Выходными данными могут быть все рассчитываемые моделью переменные - температуры областей и поверхностей электролизера, химический состав электролита, толщина настыли, высота спеченного анода, запас МГД стабильности (разность между рабочим и критическим напряжением) и т.п.
Так как математическая модель разработана для использования реальными управляющими системами, то было уделено повышенное внимание ее идентификации. Для этого использовались измерения тепловых балансов разных конструкций электролизеров, стационарные расчеты трехмерными пакетами, проводились комплексные активные эксперименты. Сравнение динамических откликов модели и объекта приведено в работе [6].
Управление химическим составом электролита
Программа «Стабилизация состава электролита на базе математической модели» [7] является предметом ежедневного промышленного применения и рассчитывает для каждого элек-
Запрос значений технологических параметров по электролизеру из базы данных
Расчет значений КО, технологических параметров к текущему моменту
1
Формирование таблицы расчетных значений технологических переменных
г
Расчет добавок АШ3. Оптимизация по общей массе и распределению сырья по дням
Вывод результата расчета добавок АШз. Передача команд исполнительным механизмам АПФ
Рис. 4. Схема работы программы «Стабилизация состава электролита на базе математической модели»
тролизера добавки сырья А№3 для корректировки криолитового отношения (КО) электролита -важнейшего параметра химического состава.
Криолитовое отношение равно отношению числа молей NaF/AlF3 в расплаве. Анализы КО выполняются 1-2 раза в неделю для каждого электролизера. Впервые в практике АСУТП электролиза мы использовали в программе встроенный динамический имитатор для расчета прогноза КО между анализами и оптимального выбора суточной дозы сырья А№3.
Необходимость расчета прогноза КО обусловлена:
• зависимостью изменений КО от теплового баланса ванны;
• инерционностью КО;
• зависимостью потерь А№3 от параметров электролиза и свойств сырья.
Например, если повышать напряжение электролизера, то повысится температура электролита, как следствие, расплавится настыль, КО повысится, так как криолитовый модуль настыли выше, чем электролита, и мы должны отреагировать на это увеличением подачи А№3. Последствия увеличенной отдачи А№3, изменений напряжения сказываются полностью через 3-5 суток.
На рис. 4 приведена схема работы алгоритмов программы. После запроса необходимых данных программа рассчитывает КО и другие параметры электролизера на текущий момент, затем по этим рассчитанным параметрам составляет задание по добавке сырья на три последующих дня.
£Ц;таСшиыция («МВ ЗЛСКТрО№"13
Вижзд кибсрхдлтуы (Икнет гга щуппв | Риме г по элет ролилфу |
I_I РЯП КД |рлотви»111- с гшипплпп. % [их [ЭЬ.«
Рв£ч*т
Корпус П. I рцигш "лп ии. ЬВ ?0.
11РЧЛТЬ ЙППНКЛ ЗДОДАШЙ «я дпгр^лк.у —>
J
J
сдр|ч1нсш10 номиий а]г) л ол —г пр^ть | 1упнять | г^ринять |
м ЛГрр икшвд т<ж юо Прог.ки Темп ал та <СЗ 0(4 к слу*£ы (нОС) № ггер дРЭ (И »41 (*Г) 4-1 -Гп Нодидм (ЧРЮП
Итчяпо 7Л6 гм »я.; 1*1.4 474
350 3 4Я1 356 14И н 1« 17 £ ИЛ
ЭЙ не 2А 3.411 9Н 35 И 1Ь.) к.г 1Л ЫА
36* т 2» 2366 2121 18.1 1П 123
347 255 гя1 1.ве 956 г.я 411 17.3 11Л 113
34 ив а 3.4Ш № 05.21 ял ИЗ зи 3<С
35& ? 45 3457 9Я1 Я!ЛП 1В.4 ■л ЕМ
356 ? 45 7» 3 7В1 жп ига А с II Л
357 г-15 3 417 ВД 37 77 ш ня зи 35.1
358 2.48 3.41 3.4« 45$ 18.1 (и 1П 123
ли ик 1 м гя 953 i4.ua »Л М.4 ял 223
1рнпечонгч1 НйЩнТ е и - л |ку г ллн р<иСч&1 о ни: А. Г- * Д||>| | онццки (ки.'чш □ кшуку С_иг|. Г и ш -^л иь н'н райхе* □ ^¡ртИОНбЖнЧН ОР1ШПЩ ШМ4Н 1м4и4мЛВДЫ0№ЩЛ1*МТС«»*№ОА1Т Поел» рпгифтл. Пипсииптл иплпштьАМу ЛРЗ 1 169ЛМ|Й 11пн ийЭн^ага к ний йОаи & ПриЙия ганна* ЛОЗ^а аьшечини шиич ЦНТйгъ шивпенноя ируииуй <рППШН ! к ! 1ИТГ РУ'ППКу Рр1|| 1ЧТ1. Л^^ ГЛлрЯ! 1ГччИ Гц |"ЧГ1ч.:к"| 0 Г'.П ОчС^ Я™ ОЫОМИ ГМГр^ПКу "1ИТГ НЫ^Пку ГГлЧПТ*.
Рис. 5. Примеры расчетов программой стабилизации состава электролита
Программа должна ликвидировать разницу между прогнозом КО на сегодняшний день и целевым значением КО. Для этого используется оптимизационная процедура выбора общей массы добавки сырья и распределения этого сырья по дням. От разных способов добавки сырья меняется прогнозируемый тренд КОфакт в целевой функции
2 (КОЧелъ - КОфакт )2
СКО = 1 —
N
Ш1П ,
(4)
где N - общее число расчетных шагов.
Для оптимизации применяется метод двумерного градиента с одномерным спуском. В пространстве двух параметров - массы сырья и распределения этого сырья на суточные порции - определяется градиент, и в направлении антиградиента осуществляется «спуск» до тех пор, пока целевая функция уменьшается, или до заданной точности. Результаты расчетов выводятся либо в виде специальной таблицы для технологов корпусов, не оборудованных системой автоматической подачи фторсолей (АПФ) - рис. 5 , либо передаются в виде команд на исполнительные механизмы АПФ.
На рис. 5 приведен интерфейс программы с выполненным расчетом сырья. Например, на электролизер №350 рассчитано добавить по 18 кг А№3 на каждый день, что близко к нормативному расходу (последний столбец), а на №353 добавки в первый и второй день 43 и 27 кг значительно превышают нормативный расход. Это связано с отличием рассчитанного значения КО от цели. В столбце «Прог. КО» - указывается расчет прогноза КО на текущий момент, в столбце
«Тек. КО» результат последнего измерения КО, в столбце «Цел. КО» - целевое значение параметра по каждому электролизеру.
О результатах внедрения программы сообщается в [8], [9]. Благодаря использованию встроенной математической модели программа смогла заменить технолога в составлении ежесуточного задания на каждый электролизер по добавке фторида алюминия - это было основной задачей разработки алгоритмов. Но программа принесла еще экономический эффект, выразившийся:
1. В улучшении качества управления - снижении среднеквадратичного отклонения КО от целевого значения не менее чем на 0,03 единицы.
2. В снижении удельного расхода А№3 в среднем на 1 кг на тонну алюминия. На некоторых опытных корпусах экономия доходила до 3 кг на тонну алюминия. Это составляет от 3 до 8 % удельного расхода сырья.
Программа внедрена на Красноярском, Новокузнецком и Хакасском алюминиевых заводах и функционирует в составе АСУТП. При этом обеспечивается выполнение технологических требований к стабильности параметра КО. Хакасский завод оборудован системой автоматической подачи фторсолей, что позволило полностью автоматизировать управление КО на этом заводе. Стандартное отклонение КО на ХАЗе не превышает 0,06 единиц, что является очень хорошим показателем качества управления. Расход сырья А№3 на всех производствах удовлетворяет техническим нормам.
Поддержка принятия технологических решений
На основе математической модели процесса нами разработана Учебно-Консультационная программа « Виртуальный электролизер» [10]. Программа позволяет посредством специального интерфейса рассчитывать и просматривать влияние управляющих воздействий или случайных факторов на параметры процесса и используется:
• для обучения персонала;
• для расчетов технологических параметров и технико-экономических показателей при планируемых изменениях процесса;
• для выбора регламентов ведения процесса.
На рис. 6 представлено моделирование программой изменений температуры электролита, концентрации глинозема и МПР при работе системы автоматической подачи глинозема (АПГ).
Внизу на графике приведен период АПГ, который говорит о частоте срабатывания дозатора глинозема: меньше период - больше глинозема попадает в электролизер, и наоборот. Видно, что при недостатке глинозема температура растет, при полном отключении питания в 06.00 произошел резкий рост температуры из-за произошедшего анодного эффекта. Это электрохимическое явление характеризуется мгновенным увеличением напряжения на электролизере приблизительно в 10 раз.
Справа от графиков в окне программы отображаются текущие расчетные значения перечисленных параметров.
На рис. 7 приведен графический результат решения технической задачи, состоявшей в подборе величины снижения заданного напряжения и добавочного питания фтористыми солями
Рис. 6. Расчет изменений температуры электролита и концентрации глинозема при работе АПГ в течение суток
Сначала:
Рис. 7. Среднесуточные изменения параметров при увеличении силы тока в течение месяца
при увеличении силы тока. При этом нужно было сохранить тепловой баланс электролизера, в частности температуру электролита и криолитовое отношение. Видно, что вначале при увеличении силы тока температура и КО выросли, но затем при снижении напряжения и увеличении добавки AlF3 были стабилизированы.
Такой подбор на реальных электролизерах приводит к расстройству технологии. В работах [11], [12] даются примеры расчетов технологических регламентов, выбора целевых значений и показывается эффективность такого способа управления.
Использование программного комплекса «Виртуальный электролизер» в итоге позволило:
1. По рассчитанным регламентам провести повышение токовой нагрузки на нескольких сериях электролиза без нарушения теплового режима электролизеров.
2. Обосновать выгодность изменения химического состава электролита, определить регламенты переходов к новым составам с наименьшими затратами.
3. Выбрать режимы системы автоматического питания глиноземом, наиболее устойчивые к изменению его свойств. Оценить возможные улучшения управления от внедрения новых алгоритмов АПГ.
Заключение
Разработаны новые методы и программные инструменты для управления электролизом алюминия, позволяющие эффективно стабилизировать химический состав электролита, выбирать целевые значения параметров процесса и регламенты проведения изменений технологии.
Список литературы
1. Bojarevics V. and Romerio M.V. Long wave instability in liquid metal-electrilyte interface in aluminium electrolysis cells: a generalization of Sele's criterion // European Journal of Mechanics, B/ Fluids, Vol. 13, N 1, 1994, 33-56.
2. Коростелев И.Н., Пискажова Т.В., Проворова О.Г., Синельников В.В. Разработка методики использования критерия устойчивости Бояревича-Ромерио в алгоритмах АСУТП электролиза алюминия // Вестник КрасГУ Серия Физ.-матем. науки. 2005, Вып 3. - С. 118-124.
3. Белолипецкий В.М., Пискажова Т.В., Емельяшин М.В. Тепловая динамическая модель для управления технологией электролиза алюминия // Вестник Сибирского государственного аэрокосмического университета им. М.Ф.Решетнева. 2006. Вып. 5 (12). С. 19-24
4. Белолипецкий В.М., Пискажова Т.В. Тепловая динамическая модель для автоматического управления технологией производства алюминия // Вычислительные технологии. 2008. Т. 13. Вестник КазНУ, №3(58). С. 267-274.
5. Berezin A.I., Isaeva L.A., Belolipetsky V.M., Piskazhova T.V., Sinelnikov V.V. A model of Dissolution and Heating of Alumina Charged by Point-Feeding System in «Virtual Cell» program // Light Metals. - 2005. - Р. 151-156.
6. Yurkov V.V., Mann V.C., Piskazhova T.V., Nikandrov K.F. Virtual Aluminum Reduction Cell // Light Metals. - 2001. - Р.1259.
7. Манн В.Х., Юрков В.В., Пискажова Т.В., Требух О.А., Никандров К.Ф. СВИДЕТЕЛЬСТВО об официальной регистрации программы для ЭВМ №2003612185 «Стабилизация состава электролита на базе математической модели», РОСПАТЕНТ 23.09.2003.
8. Пискажова Т.В., Манн В.Х., Юрков В.В., Никандров К.Ф., Требух О.А. Виртуальный электролизер и его практическое применение на ОАО КрАЗ // Алюминий Сибири - 2002: сборник статей. Красноярск. - С. 60-65
9. Piskazhova T.V., Mann V.C. The Use of a Dynamic Aluminum Cell Model // JOM, Vol. 58, № 2, February 2006, pp. 48-52.
10. Никандров К.Ф., Пискажова Т.В., Синельников В.В., Манн В.Х. СВИДЕТЕЛЬСТВО об официальной регистрации программы для ЭВМ №2003610551 «Учебно-консультационная программа «Виртуальный электролизер» «, РОСПАТЕНТ 28.02.2003.
11. Пискажова Т.В., Своевский А.В., Попов Ю.А., Бузунов В.Ю. Использование программы «Виртуальный электролизер» для анализа и выбора технологических регламентов при поднятии тока // Алюминий Сибири - 2003: сборник статей. Красноярск, 2003. - С. 132-136
12. Колесов М.С., Пискажова Т.В. Разработка технологических регламентов с помощью программного обеспечения «Виртуальный электролизер» на ОАО НкАЗ // Технико-Экономический Вестник РУСАЛа. 2004. №9. С. 11-13.
Efficient Methods for Control of Technological Aluminum Reduction Process
Tatyana V. Piskazhova
Siberian Federal University, 79 Svobodny, Krasnoyarsk, 660041 Russia
On the basis of mathematical model of process management methods by aluminum reduction technology, brought economic benefit are developed. In article algorithms of optimum stabilization of electrolyte composition and the program - the simulator of process for support of acceptance of technological decisions are presented.
Keywords: mathematical model of process, control algorithms, Virtual cell.