Постановка задачи исследования диффузионного перехода через границу шлак-металл в колонном реакторе и алгоритм ее решения Текст научной статьи по специальности «Компьютерные и информационные науки»

УДК 004.021 Дата подачи статьи: 02.03.16

Б01: 10.15827/0236-235Х.114.156-164 2016. Т. 29. № 4. С. 156-164

ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ ДИФФУЗИОННОГО ПЕРЕХОДА ЧЕРЕЗ ГРАНИЦУ ШЛАК-МЕТАЛЛ В КОЛОННОМ РЕАКТОРЕ И АЛГОРИТМ ЕЕ РЕШЕНИЯ

П.А. Сеченов, аспирант, pavesaS9@mail.ru; В.П. Цымбал, д.т.н.., профессор, tsumbal33@mail.ru; А.А. Оленников, к.т.н., доцент, tgsv-sibsin@mait.ru (Сибирский государственный индустриальный университет, ул. Кирова, 42, г. Новокузнецк, 654007, Росссия)

Предметом данного исследования являются задачи диффузионного перехода через границу шлак-металл и обезуглероживание частиц металла при прохождении через слой шлака. Авторами дана математическая постановка задач перехода через границу шлак-металл в рамках задачи имитационной модели гравитационного сепаратора струйно-эмульсионного реактора. Проведено сравнение объектно-ориентированных языков программирования (Dephi, Visual C#, ActionScript 3.0), выделены целевые задачи имитационной модели. Сделан сравнительный анализ имитационных моделей, созданных за последние 10 лет, в качестве критерия сравнения взята мерность пространства. Приведена также схема взаимодействия между классами имитационной модели и основным модулем. Особое внимание уделено классам шлака внизу, металла внизу, расчетов и основному модулю, так как именно в них реализуется поставленная задача. Для класса шлака описаны свойства (проценты массовых составляющих компонентов шлака, переменная для взаимодействия с классом расчетов) и методы (инициализация переменной класса, функции приращения и изменения состава шлака). В классе расчетов показаны функции определения балансного содержания FeO и скорости восстановления на границе шлак-металл, а затем алгоритм обезуглероживания железа шлака и его взаимодействие с другими классами программы. Для реализации задачи выбран метод имитационного моделирования (Монте-Карло), позволяющий при знании механизмов изучаемого объекта на основе разыгрывания случайности построить алгоритм. В результате реализации второй схемы - обезуглероживания частиц металла при прохождении через слой шлака - на имитационной модели были проведены опыты с одинаковыми начальными условиями и достаточным временем моделирования, равным 10 минутам.

Проведенные модельные исследования показали, что при реализации только первого механизма содержание углерода в конечном продукте составляет 4 %, а при реализации двух отмеченных выше схем - 3,5 %, что свидетельствует о возможности более гибкого управления процессом обезуглероживания.

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

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

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

Для решения задачи был выбран метод имитационного моделирования (Монте-Карло). Первоначально он использовался в основном для решения задач нейтронной физики, где традиционные численные методы оказались малопригодными, однако не утратил актуальности и на сегодняшний

день [2]. Большинство методов реализуют генерацию случайности по одной координате [3], реже встречаются задачи двухмерного [4] и трехмерного моделирования [5]. Анализ готовых средств имитационного моделирования [6, 7] показывает, что для их использования требуется знание инструментальной среды, а для гибкого взаимодействия - знание языка среды. Использование объектно-ориентированных языков программирования позволит создать более сложную систему имитационного моделирования. Авторы рассматривали объектно-ориентированные языки программирования Delphi [8], Visual C# [9] и ActionScript 3.0 [10, 11]. Язык ActionScript [12, 13] известен благодаря Flash-роли-кам и Flash-играм, в которых используется большое количество взаимодействующих объектов на сцене, поэтому он и был выбран для анимации движения частиц в колонном реакторе.

Одной из важнейших диссипативных структур, более подробно рассмотренных в [14, 15], является гравитационный сепаратор металла, шлака и газа, образующийся в верхней половине колонного реактора (зоны 6 и 7 на рисунке 1). Из-за влияния гравитационной составляющей и близкого к параболическому распределению вектора скорости дви-

Задачи подготовки исходной рабочей смеси

î à Ä 5---g .Y й-

и

__ о -

— в-

Неполное сжигание топлива и подготовка восстановительного газа

ИГ

Образование газовзвеси

ИГ

Газодинамический затвор

и создание потенциала _давления_

Целевые задачи имитационной модели колонного СЭР

6 Газ

Восстановительные процессы

на витающих взаимодействующих частицах шихтовых материалов в потоке газа

(1) (FeO)+{CO}={CO2}+[Fe] ^

(2) (FeO)+CiB={CO}+[Fe] 2

(3) {CO2}+Cie={CO} Конденсированная фаза

ИГ

X

\

Динамическая подушка

FeOr

ш

Газошлаковая эмульсия и формирование состава шлака

1 Г

в!-

Граница шлак-металл

(FeO)+[C]=[Fe]+{CO}

иг

Процессы в жидкой ванне

-Шлак—

m

Движение металла в пристенном слое

^зап = f (g )

I Fe, C]

[Fe, C]

[ Fe, C]

[Fe, C]-»- A

0

.1.

Задачи

М 1 } Восстановление оксидов железа, Ч

^—законы сохранения ..ч

1 2'' Разделение продуктов ^ К-—на металл, шлак и газ

! 2 Задача разделения процессов «V У

; в верхней и нижней частях

реактора в зависимости от \ дисперсионного состава материалов

i з )

Задачи

...................................................

з 1 ^ Расплавление частиц рудных материалов

Материальный и тепловой баланс шлака

«.. 4 J Задача перехода оксидов железа I

из шлака в металл

.......................................

^ 5 } Задача формирования состава металла, выделение и всплывание пузырьков оксида углерода

/■—V

Li;

Задача оценки макрокинетики и круговорота вещества в агрегате

Рис. 1. Целевые задачи имитационной модели колонного СЭР Fig. 1. Targets of a column jet emulsion reactor simulation model

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

Таким образом, с учетом противодействия гравитационных и аэродинамических сил в верхней половине колонного реактора (зона 6) происходит сепарация металла и шлака, а также образуется неравномерное (по плотности, газосодержанию и химсоставу) распределение параметров по высоте агрегата. Распределенность содержания оксидов железа по высоте, что играет важную роль в управлении соотношением восстановительных и окислительных процессов, обусловливается также термодинамически неравновесным характером процессов, в том числе интенсивными потоками вещества и энергии, поступающими снизу через соединительный канал 3 из реактора-осциллятора 2.

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

В процессе плавления частицы железной руды образуются частицы металла, шлака и газа. Частицы металла из-за гравитационных сил попадают в копильник (зона 8). Частицы шлака, не повстречавшие частицы газа, попадают на верхнюю границу металла и образуют слой шлака (зона 5). Данный слой шлака будет накапливаться до тех пор, пока высота шлака не достигнет соединительного канала (зона 3), который будет «сдувать» частицы шлака.

В данной статье основное внимание уделено алгоритмам двух следующих моделей:

- диффузионный переход на границе шлак-металл (зоны 5-8);

- процессы, протекающие при прохождении частицы металла через шлак.

Стоит также отметить, что ранее в имитационный модели был задействован только первый механизм [14].

Математическая постановка задачи диффузионного перехода на границе шлак-металл

На рисунке 2 показаны процессы, происходящие в зонах 5-8 колонного СЭР.

Сначала рассмотрим задачу диффузионного перехода на границе шлак-металл.

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

^ = Ршм [(ТеО)ш - ^еО)ршм ], (1)

где ^еО)ш - содержание FeO в шлаке; ^еО)ршм -содержание FeO, равновесное с металлом; Ршм -коэффициент массоотдачи из шлака в металл.

(FeO^

№

№

(2)

f (FeO)L 0,0023 где L - коэффициент распределения между шлаком и металлом; [0]рм - равновесное содержание кислорода в металле.

0,0025 • Pm

[0]рм = , ГГ1 CO, (3)

где Pco - давление в пузырьках CO; [С]м - содержание углерода в металле.

2. В результате диффузионной силы, определенной по формуле (1), осуществляется переход кислорода из шлака в металл, диссоциация FeO:

(FeO)^[Fe] + [O], а затем [O] + [C] = {CO}. (4)

Объемы компонентов углерода, железа и оксида углерода, образовавшиеся в результате диссоциации FeO, рассчитываются по формулам

12

V0 = — V0 Vc 72VFe

V0 = — V0

' Fû __' То

Vo =

V co

22,4

Vo V

М„

72 гс° 72 100

где Мм - масса металла.

3. Осуществляется подъем газовых пузырей СО, образовавшихся в результате диссоциации FeO по формуле (4).

Рис. 3. Взаимодействие между классами имитационной модели гравитационного сепаратора Fig. 3. The interaction between the classes of a gravity separator simulation model

Задача обезуглероживания частиц железа в слое шлака.

4. Из частицы железа, проходящей через слой шлака, выгорает углерод: Ашо = ткАС, где шк -масса капли; АС - процент обезуглероживания; АС= W•Тпр, где W - скорость обезуглероживания; тпр - время пребывания капли в слое шлака.

5. При выгорании углерода в шлак попадает часть оксида железа пропорционально массе сго-

72

ревшего углерода: -A^FeO = - ^c .

6. Осуществляется переход железа из шлака в металл:

AmFe = J6 Amc . (5)

7. Образуется газовый пузырь в результате сго-

22,4

рания углерода: AVco = Amc.

Структура имитационной модели гравитационного сепаратора

Прежде чем переходить к алгоритмам на границе шлак-металл, рассмотрим общую структуру имитационной модели (рис. 3).

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

телем, а отображение графиков - в соответствующем классе.

Стоит также отметить, что в основном модуле и классах конденсированных частиц и газов перерасчеты происходят 30 раз в секунду (например, определение нового положения частицы в колонном реакторе). На рисунке 4 показан алгоритм работы основного модуля, в котором выделены два этапа: ввод и инициализация данных, а также вызов основной функции МБ.

Ввод данных пользователем

г

function MF(e:Event) - основная функция (выполняется 30 раз в секунду)

Основные расчеты функции

Графическое отображение полученных данных (один раз в секунду)

Рис. 4. Алгоритм работы основного модуля Fig. 4. The basic module algorithm

В основных расчетах функции МБ() выполняется ряд задач:

- инициализация переменных и массивов;

- добавление новых частиц с заданной дискретностью;

- расчет плотности по газу;

- расчет изменений состава шлака и металла при диффузионном переходе на границе шлак-металл;

- образование частиц шлака, если высота шлака достигла входного отверстия потока газовзвеси;

- расплавление частицы и образование металла, шлака и газа;

- учет изменений металла и шлака при достижении частицей металла нижней границы шлака;

- удаление частиц, достигших выходного отверстия;

- расчеты упругих и неупругих ударов между частицами;

- определение плотности слоев по высоте канала;

- накопление статистических данных об исследуемом объекте и их отображение.

Алгоритмы процессов на границе шлак-металл

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

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

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

Как видно из рисунка 5, в функции определения нового состава шлака идет вызов функции класса расчетов, в котором по входным данным (массам и массовым долям заданного компонента) рассчитывается новая массовая доля в результате прихода или ухода какого-то вещества. Для шлака вызов данной функции происходит 12 раз, так как в шлаке учитываются 12 компонентов.

Рис. 5. Свойства и методы класса шлака внизу Fig. 5. Properties and methods of the "Slag bottom" class

Рис. 6. Алгоритм диффузионного перехода на границе шлак-металл

Fig. 6. The diffusion transition algorithm on the slag-metal border

Устройство класса металла внизу аналогично классу шлака внизу, за исключением компонентов, входящих в металл Fe, С, Мп, Р).

Алгоритм диффузионного перехода показан на рисунке 6.

Алгоритм диффузионного перехода на границе шлак-металл (рис. 6) состоит из определения массы шлака, удаления FeO из шлака, удаления С из копильника, добавления Fe в копильник, образования газового пузыря СО, визуального отображения изменения высоты шлакового слоя и отображения образовавшейся газовой частицы.

Рассмотрим более подробно алгоритм расчета изменения массы шлака (рис. 7). Для расчета массы шлака, прореагировавшего в результате диффузионного перехода, необходимо знать балансное содержание FeO (по формуле (2)) и скорость восстановления на границе шлак-металл (по формуле (1)). Эти две функции расположены в классе расчетов. Для знания балансового содержания FeO требуется передать в функцию содержание углерода. А для расчета скорости обезуглероживания передаются содержание FeO в шлаке и балансное FeO (формула (3)), рассчитываемое в соответствующей функции.

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

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

1) определение массы частицы железа в килограммах;

Рис. 7. Алгоритм определения массы прореагировавшего шлака на границе шлак-металл Fig. 7. The mass determining algorithm for unreacted slag in the slag-metal border

_mFeBefore = calc.vSphere(v[i].ball_r) * v[i].pt * kV;

_mFeBefore_rr = calc.vSphere(v[i].rr) * v[i].pt;

3

_dFe = calc.vSphere(v[i].ball_r) * newptMetal * kV; _dme.dmeComposition(v[i].mS, v[i].mFe, v[i].mC, v[i].mMn, v[i].mP, _dFe);

10

T

3 timeFe = calc.timeFeinSlak(v[i].pt, 2 * v[i].rr);

t

4 v[i].mC -= timeFe * 0.02;

i

5 firstMe.initMe(v[i].mS, v[i].mFe + 56 / 12 * _timeFe * 0.02, v[i].mC, v[i].mMn, v[i].mP, 0); newptMetal = firstMe.newPtMetalOne(firstMe);

V

6 v[i].ball_r = v[i].ball_r / v[i].rr * calc.rMP(calc.vSphere(v[i].rr) * v[i].pt * (1 + (56 - 12) / 12 * _timeFe * 0.02), newptMetal);

*

7 dDownSlak = 72 / 12 * mFeBefore * timeFe * 0.02; _dsl.dSlComposition(72 / 12 * _timeFe * 0.02, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -_dDownSlak); _sldown.newComposition(_dsl);

t

8 g.push(newGas(2, calc.rMP(_mFeBefore_rr * _timeFe * 0.02 * (28 / 12), 1.25 * 400) * _kmc, true)); add_Particle(_gx, downruda + 25, _gz, g);

_medown.newComposition(_dme);

Рис. 8. Алгоритм обезуглероживания частицы железа в слое шлака Fig. 8. Decarburization algorithm for iron particles in a slag layer

1

2

9

2) определение массы частицы в модельных единицах измерения;

3) определение времени пребывания частицы в шлаке по плотности и радиусу;

4) удаление углерода из частицы железа с учетом скорости выгорания углерода и времени пребывания частицы в железе;

5) инициализация выделившегося из шлака металла (формула (5)), вычисление новой плотности металла;

6) определение нового состава шлака с учетом того, что из него выделился БеО;

7) создание газовой частицы СО;

8) учет прироста массы железа, пришедшей из шлака;

9) определение нового состава металла в ко-пильнике.

Как видно из рисунка 7, в алгоритме используются данные о частицах железа, которые берутся из класса конденсированных частиц (пп. 1-6, 9). Принадлежность к частицам железа определяется по идентификатору.

В пп. 1-3 происходит обращение к классу расчетов, в котором рассчитываются масса частицы по радиусу и плотности (пп. 1 и 2) и время нахождения частицы металла в шлаке (п. 3).

Время растворения влияет на выгорание углерода из частицы металла (п. 4), определение выделившегося из шлака металла (п. 5), радиуса частицы металла (п. 6), состава шлака (п. 7) и объема частицы газа, образовавшейся в результате реакции (п. 8).

С учетом пришедшего из шлака металла (п. 9) и частицы металла определяется состав металла в ко-пильнике (п. 10).

Алгоритм обезуглероживания частицы металла в шлаке взаимодействует с другими классами программы: классом расчетов (пп. 1-3, 6, 8, 9), классом конденсированных частиц (пп. 1-6, 9), классом газовых частиц (п. 8), классом металла (пп. 5, 9), классом шлака внизу (п. 7), классом металла в ко-пильнике (п. 10).

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

мощностью компьютера. Такой подход обладает гибкостью, так как разработанные классы (или некоторые функции классов) могут быть использованы при разработке других задач. В результате проведенных модельных исследований показано, что при реализации только первого механизма содержание углерода в конечном продукте составляет 4 %, а при реализации двух отмеченных выше схем - 3,5 %, что свидетельствует о возможности более гибкого управления процессом обезуглероживания.

Литература

1. Цымбал В.П., Мочалов С.П., Рыбенко И.А. и др. Процесс СЭР - металлургический струйно-эмульсионный реактор. М.: Металлургиздат, 2014. 488 с.

2. Kroese D.P., Brereton T., Taimre T., Botev Z.I. Why the Monte Carlo method is so important today. WIREs Comput Stat, 2014, no. 6, pp. 386-392.

3. Милованов М.М. Разработка инструментальной системы имитационного моделирования бизнес-процессов предприятия // Моделирование, программное обеспечение и наукоемкие технологии в металлургии: сб. докл. III Всерос. науч.-практич. конф. Новокузнецк: Изд-во СибГИУ, 2011. С. 262-266.

4. Меликов А.В. Имитационная модель оценки эффективности стрельбы подразделений РСЗО // NB: Кибернетика и программирование. 2014. № 3. С. 91-102.

5. Гришко М.С., Безносюк С.А., Жуковский М.С. Компьютерное моделирование методом Монте-Карло процессов фор-

мирования кластеров никеля в нанопорах материалов // Изв. Алтайского гос. ун-та. 2007. N° 3. С. 71-75.

6. Карпов Ю.Г. Имитационное моделирование систем. Введение в моделирование с AnyLogic 5. СПб: БХВ-Петербург, 2009. 400 с.

7. Лычкина Н.Н. Современные технологии имитационного моделирования и их применение в информационных бизнес-системах и системах поддержки принятия решений // Имитационное моделирование. Теория и практика (ИММОД-2005): сб. докл. II Всерос. науч.-практич. конф. СПб: Изд-во ЦНИИТС, 2005. Т. 1. С. 25-31.

8. Кожемяченко В.И. Программирование в среде Delphi. В 2 ч. Новокузнецк: Изд-во СибГИУ, 2015. Ч. 1. 32 с.

9. Кожемяченко В.И. Программирование в среде Visual С#. В 2 ч. Новокузнецк: Изд-во СибГИУ, 2012. Ч. 2. 54 с.

10. Изучение ActionScript 3.0. URL: http://help.adobe.com/ ru_RU/as3/learn/as3_learning.pdf (дата обращения: 21.02.2016).

11. Справочник по AcionScript 3.0 для платформы Adobe Flash. URL: http://help.adobe.com/ru_RU/FlashPlatform/reference/ actionscript/3/ (дата обращения: 21.02.2016).

12. Лотт Дж., Шалл Д., Питерс К. ActionScript 3.0. Сборник рецептов. СПб: Символ-Плюс, 2007. 608 с.

13. Мук К. ActionScript 3.0 для Flash. Подробное руководство. СПб: Питер, 2009. 992 с.

14. Сеченов П.А. Алгоритм и программная реализация имитационной модели гравитационного сепаратора колонного струйно-эмульсионного реактора // Программные продукты и системы. 2015. № 3 (111). С. 214-219.

15. Цымбал В.П., Сеченов П.А., Оленников А.А., Падал-ко А.Г. Имитационное моделирование гетерогенного потока и проблема уравнения Навье-Стокса // Соврем. науч. достижения металлург. теплотехники и их реализация в пром-ти: сб. докл. Междунар. научн.-практ. конф. Екатеринбург: Изд-во УрФУ, 2015.С. 357-365.

Software & Systems Received 02.03.16

DOI: 10.15827/0236-235X.114.156-164 2016, vol. 29, no. 4, pp. 156-164

STATEMENT OF THE RESEARCH PROBLEM OF THE DIFFUSION TRANSITION ACROSS

THE SLAG-METAL BORDER IN A COLUMN REACTOR AND ITS SOLUTION ALGORITHM

P.A. Sechenov1, Postgraduate Student, pavesa89@mail.ru

V.P. Tsymbal1, Ph.D. (Engineering), Professor, tsumbal33@mail.ru

A.A. Olennikov1, Ph.D. (Engineering), Associate Professor, tgsv-sibsin@mail.ru

1 Siberian State Industrial University, Kirov St. 42, Novokuznechk, 654007, Russian Federation

Abstract. The subject of the research is the problems of the diffusion transition across a slag-metal border and metal particles decarbonization when passing through a slag layer. The authors give a mathematical problem statement through the slag-metal border problems within a simulation model of a jet emulsion reactor gravity separator. The article compares the object-oriented programming languages (Delphi, Visual C #, ActionScript 3.0). It allocates simulation model targets. There is a comparative analysis of simulation models created over the past 10 years. Dimension space is taken as a comparison test (one-, two- and three-dimensional). The authors give the scheme of interaction between simulation model classes and the main unit. Particular focus is on the classes: "slag at the botto", "metal at the bottom", calculations and the main unit since they implement the assigned task. The paper describes the properties of the slag class (percentages of slag components' weight and the variable for interaction with the settlement class) and methods (initialization of the class variable, increment function and slag composition change). The calculating class contains the functions determining FeO balanced content and the recovery rate on the slag-metal border. This article shows an algorithm of slag iron decarburization and its interaction with other program classes. To implement the task the authors chose the simulation method (Monte-Carlo). It allows constructing an algorithm with the knowledge of object mechanisms, based on random playing. As a result of the second scheme (metal particles decarbonization when passing through a slag layer) the authors conducted some experiments on the simulation model with the same initial conditions and sufficient modeling time (10 minutes). The simulation studies showed that when implementing only

the first mechanism carbon content in the final product is 4 %, and when implementing two above-mentioned circuits it is 3,5 %. It shows the possibility of more flexible management of a decarburization process.

Keywords: diffusion transfer, simulation model, ActionScript 3.0, gravity separator, the interaction of particles, algorithm, software implementation, jet-emulsion reactor, Monte-Carlo method.

References

1. Tsymbal V.P., Mochalov S.P., Rybenko I.A. Protsess SER - metallurgichesky struyno-emulsionny reactor [JER Process - Metallurgical Jet Emulsion Reactor]. Moscow, Metallurgizdat Publ., 2014, 488 p.

2. Kroese D.P., Brereton T., Taimre T., Botev Z.I. Why the Monte-Carlo method is so important today. WIREs Comput Stat. 2014, no. 6, pp. 386-392 (in Russ.).

3. Milovanov M.M. Development tool system simulation of business processes. Sb. dokladov III Vseross. nauch. -praktich. konf. "Modelirovanie, programmnoe obespechenie i naukoemkie tekhnologii v metallurgii" [Simulation Software and High Technologies in Metallurgy: Proc. 3rd All-Russian Science and Practice Conf.]. Novokuznetsk, SibSIU Publ., 2011, pp. 262-266 (in Russ.).

4. Melikov A.V. Simulation model of evaluating efficiency of Multiple Launch Rocket System units firing. Kibernetika iprogrammirovanie [Cybernetics and Programming]. 2014, no. 3, pp. 91-102 (in Russ.).

5. Grishko M.S., Beznosyuk S.A., Zhukovsky M.S. Computer simulation using Monte-Carlo method for nickel clusters formation processes in material nano-pores. Izvestiya Altayskogo gosudarstvennogo universiteta [The News of Altai State University]. 2007, no. 3, pp. 71-75 (in Russ.).

6. Karpov Yu.G. Imitatsionnoe modelirovanie sistem. Vvedenie v modelirovanie s AnyLogic 5 [Simulation systems. Introduction to Modeling with AnyLogic 5]. St. Petersburg, BHV-Peterburg Publ., 2009, 400 p.

7. Lychkina N.N. Modern technologies of simulation modeling and their application in information business systems and decision-making systems. Sb. dokl. II Vseross. nauch.-praktich. konf. "Imitatsionnoe modelirovanie. Teoriya i praktika" (IMMOD-2005) [Simulation. Theory and Practice: Proc. 2nd All-Russian Science and Practice Conf. IMMOD 2005]. St. Petersburg, TsNIITS Publ., 2005, vol. 1, pp. 25-31 (in Russ.).

8. Kozhemyachenko V.I. Programmirovanie v sredeDelphi [Programming in Delphi]. 2 parts. Part 1. SibSIU Publ., 2015, 32 p.

9. Kozhemyachenko V.I., Krasnoperov S.Yu. Programmirovanie v srede Visual C# [Programming in Visual C#]. In 2 parts. Part 2. Novokuznetsk Univ. SibSIU Publ., 2012, 54 p.

10. Izuchenie ActionScript 3.0 [Learning ActionScript 3.0]. Available at: http://help.adobe.com/ru_RU/as3/learn/ as3_learning.pdf (accessed February 21, 2016).

11. Spravochnikpo ActionScript 3.0 dlyaplatformy AdobeFlash [ActionScript 3.0 Reference for AdobeFlash Platform]. Available at: http://help.adobe.com/ru_RU/FlashPlatform/reference/actionscript/3/ (accessed February 21, 2016).

12. Lott J., Schall D., Piters K. ActionScript 3.0. Cookbook. O'Reilly, 2006 (Russ.ed.: St. Petersburg, Symbol plus Publ., 2007, 608 p.).

13. Moock C. ActionScript for Flash MX: The Definitive Guide. O'Reilly, 2002 (Russ.ed.: St. Petersburg, 2009, 992 p.).

14. Sechenov P.A. The algorithm and software implementation of a simulation model of the gravitational separator column jet emulsion reactor. Programmnyeprodukty i sistemy [Software & Systems]. 2015, no. 3, pp. 214-219 (in Russ.).

15. Tsymbal V.P., Sechenov P.A., Olennikov A.A., Padalko A.G. Simulation modeling of heterogeneous flow and the problem of the Navier-Stokes equations. Sb. dokl. Mezhd. nauch.-praktich. konf. "Sovremennye nauchnye dostizheniya metallurgicheskoy teplotekhniki i ikh realizatsiya v promyshlennosti" [Modern scientific achievements of metallurgical heat engineering and their industrial implementation: Proc. Int. Science and Practice Conf.]. Ekaterinburg, UrFU Publ., 2015, pp. 357-365 (in Russ.).

Примеры библиографического описания статьи

1. Сеченов П.А., Цымбал В.П., Оленников А.А. Постановка задачи исследования диффузионного перехода через границу шлак-металл в колонном реакторе и алгоритм ее решения // Программные продукты и системы. 2016. Т. 29. № 4. С. 156-164; DOI: 10.15827/0236-235X. 116.156-164.

2. Sechenov P.A., Tsymbal V.P., Olennikov A.A. Statement of the research problem of the diffusion transition across the slag-metal border in a column reactor and its solution algorithm. Programmnye produkty i sistemy [Software & Systems]. 2016, vol. 29, no. 4, pp. 156-164 (in Russ.); DOI: 10.15827/0236-235X.116.156-164.