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

ТЕПЛОЭНЕРГЕТИКА МЕТАЛЛУРГИЧЕСКОЙ ОТРАСЛИ

УДК 542.49 (669. 784: 536.2.072) DOI: 10.18503/1995-2732-2016-14-2-116-125

СОВЕРШЕНСТВОВАНИЕ И ИДЕНТИФИКАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПРОКАЛКИ УГЛЕРОДНЫХ МАТЕРИАЛОВ В БАРАБАННОЙ ВРАЩАЮЩЕЙСЯ ПЕЧИ

Герасименко Т.Е., Мешков Е.И., Томилина A.A.

Северо-Кавказский горно-металлургический институт (государственный технологический университет), Владикавказ, Россия

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

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

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

Приведены полученные с применением разработанных алгоритмов и компьютерной программы данные вычислительного эксперимента прокалки нефтяного кокса. Параметрической идентификацией с применением данных промышленной печи подтверждена адекватность усовершенствованной математической модели. Отклонения расчетных параметров от параметров реального процесса составляют менее 5%. Исследования выполнены на кафедре «Теория и автоматизация металлургических процессов и печей» СевероКавказского горно-металлургического института в 2014-2015 годах. Математическая модель и программа эффективно применяются для исследования процессов вычислительными экспериментами. Усовершенствованная комплексная математическая модель процессов термообработки в барабанной вращающейся печи и разработанная компьютерная программа для ее реализации могут быть использованы для исследования широкого круга технологических процессов термообработки как инертных, так и неинертных материалов в барабанных вращающихся печах.

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

Введение

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

Эффективным методом решения этой задачи является математическое моделирование тепловой

© Герасименко Т.Е., Мешков Е.И., Томилина A.A., 2016

работы печи и протекающих физико-химических процессов. Математическое моделирование процессов термообработки в барабанных вращающихся печах все чаще используется как в России [1-4], так и за рубежом [5, 6].

Методика моделирования и исследования численными экспериментами обжига в режиме противотока инертного материала - шамота на основе зонального метода описана в [1]. Этот метод основан на широко применяемых нелинейных уравнениях зонального теплового баланса

К ^.т4+ёк,?)+, - б* = о,

к

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

Вт/К4; Т - зональная температура, К; gJd и g0, - коэффициенты конвективного обмена, Вт/К, Вт; ^ , суммарный зональный (результирующий) тепловой поток, Вт.

Методы исследования

С целью применения зонального метода печь по длине разбивается на определенное рациональное число п участков (рис. 1). На каждом участке выделяется по три зоны: на поверхности обжигаемого материала (материал), на поверхности футеровки печи (футеровка) и объёмная газовая зона (газ). Выделяются также две воображаемые поверхностные зоны по торцам барабана печи (торец 1 и торец 2). Зоны участков имеют следующие номера: , = 1,..., п - зоны материал; , = п +1, ..., 2п - футеровки; , = 2п +1, 2п + 2 - торцов и , = т,..., I - газа, где т = 2п + 2 - количество поверхностных зон, а I = т + п - общее количество зон.

7=1 17=2 ' } 11 >=л-1| 7=п

Рис. 1. Схема выделения участков и зон в печи: 1, 2,...,], п - номераучастков;

1, 2, п, т, I - номеразон; I - футеровка, II - материал, III - газ;

R - внутренний радиус печи, р - центральный угол сегмента материала

Для каждой зоны составляется уравнение теплового баланса (1). Система этих нелинейных уравнений (например, для двадцати участков составляется 62 уравнения), где неизвестными являются зональные температуры, решается итерационным методом. После задания первого приближения зональных температур на каждом шаге итерационной процедуры уточняется значение этих температур, одновременно рассчитываются все величины, входящие в уравнения математической модели, и составляется новая система уравнений. Для получения конечных значений зональных температур обычно требуется 7-9 шагов итерационной процедуры.

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

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

Учитывая, что при прокалке углеродных материалов в барабанной вращающейся печи, кроме горения топлива и изменения свойств прокаливаемого материала, горят летучие вещества, образующиеся в результате пиролиза органики [7], и прокаливаемый материал, прототип - математическая модель [1] дополнена нами (рис. 2) подсистемой уравнений физико-химических превращений [2].

Рис. 2. Структура комплексной математической модели

Система уравнений комплексной математической модели прокалки углеродных материалов приведена в работе [8].

Коэффициенты радиационного обмена в уравнении (1) определяются по известным формулам [1]. Наибольшую трудность при этом составляет расчёт угловых коэффициентов излучения в рабочем пространстве барабанной вращающейся печи. Вопросы их автоматизированного расчёта рассмотрены в публикациях [9,10].

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

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

Результаты исследования и их обсуждение

Настоящая работа направлена на совершенствование математической модели теплообмена с целью учета возникающих конвективных тепловых потоков при термообработке углеродных материалов. Рассмотрим вывод машинно-ориентированных алгоритмов расчёта коэффициентов конвективного обмена,/ и g01 на примере объёмных зон (рис.3).

Рис. 3. Схема конвективных тепловых потоков, источников и стоков для объемных зон: I - номер объемной газовой зоны;^ - номер смежной зоны на поверхности материала; ]2 - номер смежной зоны на поверхности футеровки; ]з - номер предыдущей по ходу газа смежной объемной зоны; м, ф, Т, л - индексы, обозначающие материал, футеровку, топливо и летучие соответственно; а - коэффициенттеплоотдачи, Вт/(м2-К); в - массовый расход, кг/с; сР- удельная теплоемкость при постоянном давлении, Дж/(кг-0С); Г - температура, °С; Р - площадь поверхности зоны, м2; Ав - изменение массового расхода, кг/с; 0 - мощность теплового потока, Вт

В соответствии с рис. 3 и с учетом замены температуры ^ на (Т-273) суммарный конвективный тепловой поток объемных зон, кроме первой по ходу газа, выражается следующим образом:

б; =а-Е; ? - ?г) + (?* - ?г) + (Тз - 273) + (<л + с^ \ьо*\)(?? - 273)--сгр А (?г - 273) .

В результате разделения переменных получим

б = (а; Е; + < л \ач\+^ | до; |) ? + а * ^ + 3 ед - (+ «; Е + а % Л) Тг --273(с;,^ -с;,,с; + с;,л|дс; |+с;,л\ьо-|).

В соответствии с этим выражением алгоритмы расчёта коэффициентов конвективного обмена уравнения (1) для объёмных зон (,=т+2, ..., I) имеют вид

= -273 («^Л - +с; л 1д^ |+),

gki ~

ет + с; ^ ^+^

а% Р?2

-

при к Ф j 1,к Ф } 2,к Ф У3,к Ф,;

при к =

при к = у2;

при к = Л3;

при к =, .

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

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

- Первой по ходу газа объемной зоны

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

0 при к Ф ]\,к Ф ]2,к Ф 2п +1,к Ф

а^ + ¡ДО^ + с;,я|дО;| при к = ]{;

% ] к = ]2;

gki ~

а

<2 + с; ,2 „+1 + с;,2„+о - * к = 2п+1

при

(сгр о +а л р; + а* ) при к =1,

где п - индекс, обозначающий подсос воздуха в горячую головку печи.

- Поверхностных зон материала (1=2, ..., п-1), кроме первой и последней

= -273Г

V P, ]2 ]2 Р^ 1 ) 1 | Р,Ч 1 \)

Здесь и далее верхний знак относится к режиму прямотока, а нижний - к режиму противотока, и ' -индексы, обозначающие начало и конец зоны соответственно; ]1 и ]2 - номера зон смежных с 1-й зоной объемной зоны и предыдущей зоны материала соответственно.

0, при к Ф ]1,к Ф ]2,к Ф1,к Ф1 +1;

±с". О" /2,

Р, ]2 ]2 / ,

/„ . ч/ ПРИ к = ]г;

±(с". С" -)/2-ам-

V Р, ]2 ]2 Р,1 1 Н 11

-с" 1доЛ1 - с1доЧ ПРИ к=1;

Р ,М М Р ,М 1 ' _ ''м^м /'-)

+СР/2, при к = 1 +1;

Мт^М а1 ^ ,

при к = ]\.

Вывод алгоритмов прокалки в режиме противотока газа и материала выполнен при подаче материала в зону с номером п, а топлива в зону 2п+1. - Первой и п-й зон материала (1=1, 1=п): для режима прокалки в прямотоке

g0, = ^РТт - 273

где Р - производительность печи по материалу в загрузке, кг/с; Тт - температура этого материала, К. g0 = -273

о0п \ р,п п-1 р,п п р,п п р,п п / 5 4 '

Совершенствование и идентификация математической модели

Герасименко Т.Е., Мешков Е.И., Томилина А.А.

§кк

" М /~1М /л Ср,п-1°п-1/2,

(с'м СМ1 - 3с"м Ом )/2-

у р,п п-1 р,п п }

-а*Е " - сл АОл - см АО * ,

п п р ,п п р ,п п '

м рм ХпГп

для режима прокалки в противотоке

?0| = -273(с"м.О" - с'".в," - с". (ДО/ I - с". (ДО." I) . >01 \ р,1 1 р,1 1 р,1 1 р,1 1 / -

при к Ф ]1, к Ф п -1, к Ф п; при к = п — 1;

при к = п ; при к = ^

<к1 _

(" о " м/~~1м \ /Ъ _ М ТТЛ

СРЛ°2 - 3сРЛО1 V2 ~а1р1

(^ + СрО )/2,

МТ7М а1 р1 ,

при к Ф ]1, к Ф1, к Ф 2; при к = 1;

при к = 2; при к = ]

ОмТт - 273 (с'"Р - с'м Ом - с" |ДОл| - с" |ДОм I).

р,п п \ р,п р,п п р ,п п р,п п / '

§кп

" М /~1М /л ~Ср,п-1Оп-1/2,

-( см 1Ом1/2 + а " Е м

^ р, п-1 п-1 I п п

при к Ф ]1, к Ф п -1, к Ф п; при к = п — 1;

+ср ,п 1АО,

м рм ХпГп

при

к = п ;

при к = ]\.

<к1 ~

- Поверхностных зон футеровки (1=п +1,., 1=2п) при к Ф1,к Ф1 +п + 2; при к = 1; при к = 1 + п + 2.

0,

- а? , а? рф,

- Поверхностных зон на торцах печи для зоны на первом торце

=стр О ттт + с"роете+спро птп - 273 (сррр+с;ои+с;ов + ср Оп),

(19) 121

-сг Ог

I p,i i

при к Ф m +1; при к = m +1;

- для зоны на втором торце

= -273вг

n у p,n p,i У' ^ '

Ski ~

при

к Ф n, к ФI;

г<г

cP,nGn, при к = n;

при к = i.

Алгоритмы (4)-(18) в отличие от прототипа учитывают конвективные тепловые потоки, возникающие при массопереносе в газовые зоны выделяющихся летучих веществ и сгорающего прокаливаемого материала, а алгоритмы (19)-(22) - не учтенные в прототипе потоки газов, участвующие в теплообмене торцовых зон.

Алгоритмы расчета коэффициентов конвективного обмена (4)-(22) использованы с целью совершенствования математической модели для синтеза системы зональных уравнений теплового баланса, которая в случае прокалки углеродных материалов имеет следующий вид:

yiaT + ST) + s0 + м |ag-| - = 0,

к

z( aT + ST)- & = 0,

к

i aT - й? = 0,

к

К aj: + ST) + s0 + q,fT| ag^ +

к

(+1 -км)QPH " |ag"i + QPH,л |AG?|,

i=1, ..., n; i=n+1, ..., 2n;

i= 2n+1, m; i=m, ..., l,

где км - доля теплового потока от горения материала, потребляемая на его нагрев; Qn - тепловой поток потерь.

С целью реализации комплексной математической модели [8] и проведения численных экспериментов прокалки углеродных материалов разработана компьютерная программа в операционной среде визуального программирования Delphi XE2. Данными таких экспериментов являются поле зональных температур и значения двадцати параметров в виде распределения их по длине печи. При этом обеспечивается возможность варьирования пятидесяти входных технологических и конструктивных параметров.

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

Значения технологических контролируемых параметров промышленного процесса и вычислительного эксперимента были одинаковыми: время прокалки 60 мин, производительность по загрузке 3,1 кг/с и расход природного газа 0,12 кг/с. Все величины технологических параметров получены усреднением за сутки 96 значений и являются их математическими ожиданиями. Некоторые результаты вычислительного эксперимента при этих условиях приведены на рис. 4, а данные идентификации - в таблице.

Качество прокаленного кокса оценивается по величине его удельного электросопротивления, которое в процессе прокалки снижается с 10000 мкОм • м до соответствующего значения техническим требованиям в 1000 мкОм • м.

ей

Рн &

О

е з

о

2к 10"

1.5x10"

1x10

500

„— ' 1 ' * ! 1

/ < ( Лч ! \ -

ч ♦ л ♦ ч ♦ ч ч ч ч ч ч ч ч { * * * $ * -

\ /

\ / * / >*

0.1

0.08

0.06

0.04

0.02

10

20

30

40

-Температура материала

.....Температура газа

Расход угара - дифференциальная функция

-----Расход летучих - дифференциальная функция

Удельное электросопротивление материала, 0,00001 *мкОм*м

и Я "С: Е

к £

В

ей О. Н

к

о

о

5

и о X о а

Си

ч. X :

\ \ • *• \*

Ч \ 1 ч - ' 1 V \ "

* * » ^ .....____.. \ \ \\ \ч \

0.1

0.05

(Я

кГ ш

К

§

с

о

е-

«

о

х

а ^

Си

10

20

0

30

-Концентрация кислорода

.....Угар - интегральная функция

Сгоревшая часть летучих - интегральная функция

-----Несгоревшая часть летучих - интегральная функция

Летучие - интегральная функция ..... Топливо

Рис. 4. Распределение параметров по длине печи, работающей в режиме противотока

Данные параметрической идентификации математической модели

Наименование параметра Размерность Значение параметра Отклонение расчетного параметра, %

промышленной печи численного эксперимента

математическое ожидание дисперсия

Температура: - кокса - отходящего газа °С 1336 428 1339 0,22

°С 710 1666,26 712 0,28

Угар кокса кг/с т/ч 0,23 0,828 - 0,225 0,81 -2,17

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

Выводы

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

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

1. Арутюнов В.А., Бухмиров В.В., Крупенников С.А. Математическое моделирование тепловой работы промышленных печей / под науч. ред. В.А. Арутюнова. М.: Металлургия, 1990. 240 с.

2. Мешков Е.И. Развитие теории и совершенствование технологии производства графитированной электродной продукции на основе математических моделей массо- и теплопереноса: автореф. дис. ... д-ра. техн. наук. Владикавказ: СКГМИ (ГТУ), 2009.

3. Арунянц Г.Г., Рутковский А.Л. Математическое моделирование в задачах проектирования систем управления сложными объектами. Калининград: Изд-во ФГОУ ВПО КГТУ, 2011. 304 с.

4. Салихов З.Г., Арунянц Г.Г., Рутковский А.Л. Системы оптимального управления сложными технологическими объектами. М.: Теплоэнергетик, 2004. 240 с.

5. Khan J.A., Pal D. and Morse J.S. Numerical modeling of a rotary kiln incinerator// Hazardous Waste & hazardous Materials. 1993, 10 (1), pp. 81-95.

6. Leger C.B., Cundy V.A. and Sterling A.M. A three dimensional detailed numerical model of a field-scale rotary kiln incinerator // Environmental Science & Technology. 1993, 27, pp. 677-690.

7. Химическая технология твёрдых горючих ископаемых: учебник для вузов / под ред. Г.Н. Макарова и Г.Д. Харламповича. М.: Химия, 1986. 496 с.

8. Мешков Е.И., Герасименко Т.Е., Ковалёва М.А. // Цветная металлургия. 2012. №4. С. 53-56.

9. Мешков Е.И // Изв. вузов. Цветная металлургия. 2006. № 4. С. 71-79.

10. Мешков Е.И. // Цветные металлы. 2007. № 5. С. 36-39.

Материал поступил в редакцию 10.12.15. INFORMATION ABOUT THE PAPER IN ENGLISH

DOI:10.18503/1995-2732-2016-14-2-116-125

IMPROVEMENT AND IDENTIFICATION OF THE MATHEMATICAL MODEL OF CARBON MATERIALS CALCINATION IN A ROTARY KILN

Gerasimenko Tatiana Evgenyevna - Ph.D. (Eng.), Assosiate Professor, North Caucasian Institute of Mining and Metallurgy (The State Technological University), Vladikavkaz, Russia. E-mail: gerasimenko_74@mail.ru/

Meshkov Evgeniy Ivanovich - D.Sc. (Eng.), Professor, North Caucasian Institute of Mining and Metallurgy (The State Technological University), Vladikavkaz, Russia. E-mail: eimeshkov@gmail.com.

Tomilina Aleksandra Aleksandrovna - Postgraduate Student, North Caucasian Institute of Mining and Metallurgy (The State Technological University), Vladikavkaz, Russia.

Abstract. The work is aimed at expanding the application and improving numerical experiment based studies of heat treatment processes extensively applied to non-inert materials in rotary kilns.

The purpose of the research is to develop new computer oriented algorithms for calculation of convective heat transfer coefficients used to build heat balance equations applied in zone method simulation of heat transfer processes that take place in kilns. The article describes a case of carbon materials calcinated in parallel-flow and counterflow modes and demonstrates an algorithm building procedure for the convective heat transfer coefficients identified first for volumetric gas zones and after that for surface zones. Unlike the algorithms used before in heat treatment simulation of inert materials, these new algorithms account for the convective flows resulting from stripping of volatiles from the treated material and the interaction of the material and the volatiles with the gas flow. The algorithms designed for volumetric zones are identical for the parallel-flow and counterflow modes. As for the surface zones, some of the equations included in the algorithms are different.

The article presents numerical experiment data for a petroleum coke calcination process which were obtained with the help of the developed algorithms and a computer software. Correctness of the improved mathematical model was verified by means of parametric identification based on actual kiln data. Deviation of the calculated parameters from the actual process parameters is less than 5%. The research was carried out by the Department of Theory and Automation of Metallurgical Processes and Furnaces of the North Caucasian Institute of Mining and Metallurgy in 2014 and 2015. The mathematical model and the software are effectively applied in numerical experiment based research. The advanced comprehensive mathematical model designed for heat treatment processes of a rotary kiln and the software developed for the above model can be used to study a broad range of

heat treatment processes applied to both inert and noninert materials in rotary kilns.

Keywords: Carbon materials, calcination, rotary kiln, mathematical model, algorithm, convective heat transfer, identification, numerical experiment.

References

1. Arutyunov V.A., Bukhmirov V.V., Krupennikov S.A. Matematich-eskoe modelirovanie teplovoy raboty promyshlennykh pechey [Mathematical modeling of the thermal performance of industrial furnaces]. Ed. Arutyunov V.A. M.: Metallurgiya. 1990, 240 p.

2. Meshkov E.I. Developing the theory and improving the technology for production of graphite electrodes based on mathematical models of mass and heat transfer. Aftoreferat dissertatcii doktora tekhnicheskikh nauk [Abstract of a doctoral dissertation]. Vladikavkaz: North Caucasian Institute of Mining and Metallurgy, 2009.

3. Arunyants G.G., Rutkovskiy A.L. Matematicheskoe modelirovanie v zadachakh proektirovaniya sistem upravleniya slozhnymi obyektami [Mathematical modelling in design of management systems for complex objects]. Kaliningrad, Publ. Kalinigrad State Technical University. 2011, 304 p.

4. Salikhov Z.G., Arunyants G.G., Rutkovski A.L. Sistemy optimal-nogo upravleniya slozhnymi tekhnologicheskimi obyektami [Adequate management system for complex technological objects]. Moscow: Teploenergetik. 2004, 240 p.

5. Khan J.A., Pal D. and Morse J.S. Numerical modeling of a rotary kiln incinerator. Hazardous Waste & Hazardous Materials. 1993, 10 (1), pp. 81-95.

6. Leger C.B., Cundy V.A. and Sterling A.M. A three dimensional detailed numerical model of a field-scale rotary kiln incinerator. Environmental Science & Technology. 1993, 27, pp. 677-690.

7. Chemical technology of solid combustible minerals. Uchebnik dlya vuzov [Textbook for Universities]. Edited by Makarov G.N., Kharlampovich G.D. Moscow: Khimiya. 1986, 496 p.

8. Meshkov E.I., Gerasimenko T.E., Kovalyova M.A. Tcvetnaya metallurgiya [Non-ferrous metallurgy]. 2012, no.4, pp. 53-56.

9. Meshkov E.I. izv. Vuzov. Tcvetnaya metallurgiya [Proceedings of the Russian universities. Non-Ferrous Metallurgy]. 2006, no.4, pp. 71-79.

10. Meshkov E.I. Tcvetnye metally[Non-Ferrous Metals]. 2007, no. 5, pp. 36-39.

Герасименко Т.Е., Мешков Е.И., Томилина A.A. Совершенствование и идентификация математической модели прокалки углеродных материалов в барабанной вращающейся печи // Вестник Магнитогорского государственного технического университета им. Г.И. Носова. 2016. Т.14. №2. С. 116-125. doi:10.18503/1995-2732-2016-14-2-116-125

Gerasimenko T.E., Meshkov E.I., Tomilina A.A. Improvement and identification of the mathematical model of carbon materials calcination in a rotary kiln. Vestnik Magnitogorskogo Gosudarstvennogo Tekhnicheskogo Universiteta im. G.I. Nosova [Vestnik of Nosov Magnitogorsk State Technical University]. 2016, vol. 14, no. 2, pp. 116-125. doi:10.18503/1995-2732-2016-14-2-116-125