CC BY
73
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И
ЧИСЛЕННЫЙ АНАЛИЗ ШАХТНОЙ ПЕЧИ ИЗВЕСТКОВОГО ПРОИЗВОДСТВА
Е. П. Протасова, доцент Нижегородского государственного архитектурно-строительного
университета;
Н. Т. Пузиков, к.т.н., доцент кафедры «Механика» Нижегородский государственный инженерно-
экономический институт
На российском рынке, как показывает анализ, основными конкурентными преимуществами печных систем являются низкий удельный расход топлива, компактность установок и приемлемая стоимость оборудования [4]. Высокие расценки на импортные печи не позволяют рассчитывать на широкое их внедрение на территории России, а качество российских альтернатив недостаточно.
Большое количество работ, проводимых в настоящее время, посвящено совершенствованию старых и разработке новых конструкций отечественных шахтных печей [2,3,5]. Существует около 150 действующих российских патентов в области оборудования и способов производства извести, однако целый ряд мероприятий, эффективность использования которых установлена опытом работы некоторых разновидностей шахтных печей зарубежной конструкции [7,8], в отечественной практике шахтного обжига известняка, применения пока не нашел. Повышенное внимание уделяется и вопросам автоматизации установок обжига известняка [6]. Реализация создаваемых этими мероприятиями возможностей, разработка новых эффективных систем отопления, устройств загрузки и выгрузки, точная режимная настройка шахтных
печей требуют глубокого изучения теплотехнических особенностей шахтного обжига и их пристального анализа методами физического и математического моделирования, компьютерного расчета. Недостаточно изученными остаются и некоторые вопросы, касающиеся кинетики диссоциации известняка, аэродинамики шахты, условий движения материала, эффективности сжигания топлива, стойкости футеровок, режимной оптимизации шахтных печей. Все это сдерживает инновационную динамику и увеличивает остроту проектных рисков при реализации новых конструкторских идей в области технологии и оборудования шахтного обжига известняков.
С учетом вышеизложенного задачей исследования является синтез прогнозирующей математической модели шахтной печи с газовым отоплением на базе современных представлений о кинетике диссоциации карбонатных пород, механизме тепломассообменных процессов, аэродинамики газового потока и движения гранул, модели сгорания газообразного топлива (метана).
При анализе процессов, протекающих в шахтной печи, их условно можно разделить на два уровня: микро- и макроуровень, поэтому и общую модель шахтной печи целесообразно разбить на две части: микро- и макроки-нетическую [1]. Совокупность факторов, определяющих скорость физико-химических процессов в малых областях шахтной печи, можно отнести к микрокинетике. Микроки-нетическая группа включает уравнения состояния единичной гранулы материала, реакции горения в газовой фазе. На эффекты микроуровня оказывают влияние крупномасштабные тепловые, диффузионные,
гидродинамические явления, которые относятся к макрокинетической группе. Их специфика зависит от конструкции печи, способа и условий ввода топлива, воздуха, структуры плотного слоя и других факторов.
С точки зрения системного анализа можно выделить пять базовых уровней, на которых моделируются процессы, протекающие в шахтной печи обжига известняка:
1) физико-химическая кинетика (реакции горения топлива, испарение влаги, диссоциация карбонатов);
2) протекание процессов в отдельно взятой грануле материала (теплоперенос, диффузия, напорное течение в узких поровых каналах);
3) обстановка в пограничном слое вокруг гранулы (гидродинамика обтекания, теплообмен газами, диффузия);
4) характер взаимодействия в пределах рассматриваемой зоны аппарата (преимущественно гидродинамика);
5) обстановка в масштабе аппарата (гидродинамика печи, движение материала, теплоперенос).
Ввиду тесной связи перечисленных явлений решение задачи моделирования шахтной печи должно осуществляться совокупно на всех уровнях иерархии. Такой подход позволяет создать завершенную модель аппарата и объединить общие закономерности теплофизических процессов с конкретной технологией, режимными параметрами и конструкцией шахтной печи. Для этого были изучены и применены следующие методы:
1) оценки порозности слоя;
2) теплообмена в плотном фильтруемом слое;
3) анализа аэродинамики газового потока и гидродинамического сопротивления шихты;
4) численного моделирования движения материала;
5) расчета горения газообразного топлива;
6) диссоциации кускового известняка;
7) многопараметрической оптимизации с ограничениями.
В работе синтезирована динамическая математическая модель шахтной печи обжига известняка, структурированная по иерархическому принципу. Модель высшего ранга создана для осесимметричной геометрии. Ее
обобщенное уравнение можно записать в простой форме
- ^ \fiPfPi ]+тг ]+~ 1гтт ] - -г
ot ох г ог ох
7г¥г
д(р1
1д_ г дг
Нг¥г
д(р{
дг
= 5,
_ дх.
где t, x и г - соответственно время, вертикальная и радиальная координаты печи, с, м; i - индекс,
указывающий на принадлежность параметров твердой или газовой фазе; - объемная доля фазы; - зависимая
переменная; - коэффициент переноса (диффузионный
коэффициент); - источниковый член.
В таблице 1 дается расшифровка величин, превращающих уравнение (1) в систему балансовых уравнений неразрывности, массы, энтальпии и инерции, отражающих состояние гранулированной твердой и газообразной фаз. В таблице приняты следующие обозначения: g и 5 - индексы, указывающие на принадлежность параметра газовой или твердой фазе; gi - массовая доля j-го компонента в твердой фазе, причем i меняется от 1 до 5 (Н20, С02, 02, К2, СН4); к - энтальпия, Дж/кг; и и V -соответственно вертикальная и горизонтальная составляющие вектора скорости потока, м/с; Sj и -соответственно суммарный и покомпонентный перенос массы из твердой фазы в газовую с испаряющейся водой в зоне сушки 8=1 и выделяющимся углекислым газом в зоне обжига 8=2, кг/(м3с); vik - стехиометрический коэффициент г-го компонента газа в к-ой реакции; Як - скорость к-ой реакции, кг/(м3с); Qgm - перенос тепла между газом и материалом в результате конвекции, Вт/м3; Нк и И8д -соответственно изменение энтальпии газа в результате кой реакции в газовой фазе и выделения вещества из материала, Дж/кг; р и р5 - давление в газовой и твердой
фазе соответственно, Па; Тё и Т^0 - соответственно фактическая и опорная температура, К; ^ - плотность газа при опорной температуре, кг/м3; ё - коэффициент температурного расширения газа, К-1; g - ускорение свободного падения, м/с2; и - вертикальная и радиальная составляющая объемной силы ^ взаимодействия газа и материала, Н/м3; , -
стехиометрический коэффициент пропорциональности, причем , = 1 для зоны сушки, у = 100/44 для зоны диссоциации; g и ,, - соответственно динамический коэффициент вязкости газа и приведенный коэффициент вязкости гранулированной среды, Пас. Указанная система решается при граничных условиях I и II рода.
Таблица 1
Подстановочные выражения для уравнения (1)
Балансовое уравнение 1 5
неразрывности g 1 0 5,
массы вещества g & gDgi к
энтальпии g в,т +11{~АНк К + Е • я,
импульса X g Ug g - е(др / дх)+ер^ (Г, - Г,0) - ^
г g Vg g -б{др1дг + /и^ё1 г1)-Ргё5
неразрывности £ (1-) 1 0 - 5,
массы вещества £ (1-) & 0
энтальпии £ (1-) К - 08» + Е1 Л// ^КГ8гЪ Л//< ' 51
импульса X £ (1-) Щ £ (1 -£){р^-др5 /дх)
г £ (1-) £ -(1 -£){дрх/дг + /л^3/г2)
Состояние газовой фазы в областях печи, не заполненных материалом, может определяться по уравнениям, составленным для плотного слоя, если при использовании выражений из таблицы 1 принять: Sj = St=j = 0, Qgm = 0, FgS = 0.
Коэффициент межфазного переноса тепла при расчете Qgm определяется по известному критериальному выражению Ранца-Маршала. Движение материала описывалось моделью механики жидкости, адаптированной к механике гранулированных сред Жангом и Ю (Zhang и Yu). Силовое взаимодействие между газом и материалом учтено с помощью закона Эргуна.
Для описания процесса диссоциации кускового известняка разработана усовершенствованная
математическая модель, возникшая на базе работ В. Г. Пивоварова, В. Моффата и Р. В. Валмсли. В модели совместно рассмотрены процессы теплопроводности, конвекции и диффузии, кинетика и физические условия реакции диссоциации карбоната (рис. 1). Разработанная модель является динамической, что позволяет при переменных значениях внешних факторов получать данные о состоянии гранулы в каждый момент времени.
В ее основе лежит дифференциальное уравнение 2-го порядка для одномерной симметричной задачи
2
д (р.
— = Y
dt
д д>, | 2 д<р,
у дг2 г дг
-г—
г
dw, ди
и —- + (р, — дг дг
где ¡иг- соответственно временная и радиальная координаты, м; - зависимая переменная; и -
коэффициенты уравнения; и - скорость газа в пористой оболочке оксида кальция, м/с; / - индекс, указывающий на принадлежность параметров оксидной оболочке I = 1 и ядру I = 2 гранулы.
°° - A&JL* - ” • '.А:=JüliL i Jl •. НУ" ____________________1_. .Р ■:: »t.w lJ
Рис. 1. Тепломассообменные процессы в грануле известняка при диссоциации:
R - радиус гранулы, м; rf - радиус фронта диссоциации, м;
T- температура, К; с - концентрация CO2, моль/м3; p - давление, Па; Qgm - тепловой поток от газов к материалу, Вт; нижние индексы: g - газ; s - поверхность гранулы; f -фронт реакции; с - центр гранулы.
Система балансовых уравнений, отражающих состояние гранулы известняка, составляется на базе таблицы 2, аналогично предыдущему случаю. В таблице приняты следующие обозначения: g - плотность газа, находящегося в порах оболочки, кг/м3; T - температура, К; си p -соответственно концентрация углекислого газа в объеме пористой оболочки и давление, моль/м3 и Па; aef, ef и Cef -соответственно эффективные значения температуропроводности, плотности и удельной теплоемкости оболочки, м2/с, кг/м3 и Дж/(кгК); g и Cg и g - плотность, удельная теплоемкость динамический коэффициент вязкости газа, находящегося в порах, кг/м3, Дж/(кгК) и Пас; D -коэффициент диффузии, м2/с; k - проницаемость пористой структуры, м2.
Таблица 2
Подстановочные выражения для уравнения (2)
Уравнение i
Тепловой баланс 1 T aef gCglefCef
2 T a2 0
Матер иальный баланс 1 с D 1
Силовой баланс 1 Р k/g 0
Скорость газа в порах определяется по уравнению неразрывности с допущением, что газ несжимаем. Давление во фронте гранулы формулируется с помощью граничного условия II рода. Кинетика реакции разложения известняка рассчитывается по известной методике с учетом величины равновесного давления CO2 и температурной зависимости для константы скорости. Условие сохранения тепловой энергии во фронте реакции учитывает сток тепла при диссоциации. Также рассчитывается скорость смещения фронта реакции, объемная и массовая степени диссоциации гранулы. Описанные модели реализованы в пакете Comsol Multiphysics 3.5 a.
Для эффективного решения задачи оптимизации математическая модель шахтной печи приведена к стационарному одномерному виду и реализована совместно с моделью диссоциации гранулы и блоком многофакторной оптимизации в программной среде Matlab R2006a. В разработанной программе «LSF Wizard» при поиске решения используется метод прогонки (рис. 2). Модель идентифицирована на основе промышленных данных.
Проверка достоверности модели проводилась на основе данных о работе шахтных печей, предоставленных ОАО «Силикатный завод №1» (г. Нижний Новгород) и ЗАО «Борский силикатный завод» (г. Нижний Новгород),
на базе НПФ ООО «Промэнергогаз-2». Совокупность процессов, протекающих в шахтной печи, удовлетворительно описывается предложенной системой уравнений, отклонение расчетных значений содержания СаО в извести от экспериментальных и производственных данных не превышает 2 %. Математическая модель, учитывающая тепломассообмен в оксидной оболочке, химическую кинетику процесса диссоциации карбоната кальция и нестационарный характер влияния внешних факторов удовлетворительно описывает разложение гранул известняка различной крупности в широком диапазоне температур и концентраций С02 в газовой фазе.
I 3) блок обработки Выда ча результатов решяшя
результатов в аиде таблиц и графиков
Рис. 2. Алгоритм решения математической модели
171
Список литературы
1. Кафаров, В. В. Системный анализ процессов химической технологии: основы стратегии / В. В. Кафаров, И. Н. Дорохов. - М.: Наука, 1976. - 498 с.
2. Косников, А. А. Разработка математических моделей кинетики окислительно-восстановительных процессов шахтной свинцовой плавки: автореф. дис. канд. техн. наук: 05.16.03 / А. А.Косников; Ин-т металлургии и обогащения. - Алма-Ата, 1991. - 24 с.
3. Котельников, Ю. В. Комовая известь. Оборудование для обжига извести [Электронный ресурс] / Ю. В. Котельников. - Красково: Спецпромтех, 2006. - Режим доступа : http://www.spsss.ru.
4. Материалы отраслевого семинара производителей извести (пос. Угловка, октябрь 2005 г.) / ОАО «УИК». -Угловка [пос.]: УИК, 2005.
5. Шубин, А. Ф. Прямоточно-противоточная шахтная печь для производства активной извести / А.Ф. Шубин [и др.] // Сталь. - 2001. - № 5 - С. 80-81.
6. Constantinescu, R. MRAC strategy for the temperature
profile control of a lime kiln [Electronic resource] / R. Constantinescu. - Bellingham : IDEAS Simulation Inc. - Режим доступа: http://www.ideas-simulation.com/ down-
loads/LimeКiln.pdf.
7. Piringer, H. Process optimization on Maerz lime kilns / H. Piringer // ZKG Int. - 2005. - Vol. 58, № l. - Р. 41-48.
8. Predescu, L. Improving the Vertical Lime Kiln / L. Predescu // Steel Times Int. - 1997. - Vol.
