Технические науки /5. Энергетика

М.н.с Абайдуллин Б.Р.

Исследовательский центр проблем энергетики Федерального бюджетного учреждения науки Казанского научного центра Российской академии наук

Исследование влияния реологических параметров на возникновение критического режима теплообмена в коаксильном реакторе полимеризации

Аннотация

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

Ключевые слова:

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

Введение

В работе [3] рассматривались установившееся течение и теплообмен неньютоновской жидкости в коаксиальном канале бесконечной длины, со сформировавшимся профилем температуры. На практике чаще всего наблюдается изменение профиля температуры по длине, и в связи с этим наибольший интерес представляет исследование влияния основных безразмерных параметров на условия возникновения критических режимов теплообмена в канале с учетом изменения температуры по длине и по радиусу. В данном случае численная реализация поставленных задач усложняется в силу того, что необходимо учитывать конвективный теплоперенос, особенности изменяющегося неньютоновского характера поведения среды, убывание концентраций при условии химической реакции, поэтому в литературе можно встретить небольшое количество работ по данной тематике [2]. Исследование распределений температуры и компонент вектора скорости на начальном участке коаксиального канала в зависимости от интенсивности химического и диссипативного источников при условии преобладания одного из них уже рассмотрены в [4]. Исследования этих характеристик проведены и от некоторых других основных определяющих параметров в [4] и [5].

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

При построении математической модели стационарного тепломассобмена при течении неньютоновской жидкости в коаксиальном реакторе, как и в [5], приняты следующие допущения: течение ламинарное, осесимметричное, стационарное; теплофизические характеристики жидкости, такие как плотность, теплопроводность и теплоемкость меняются в рассматриваемом диапазоне температур незначительно и приняты постоянными; силы тяжести входят в уравнения неявно через избыточное давление; перенос тепла вдоль основного направления движения за счет теплопроводности много меньше вынужденного в этом же направлении; профиль скорости почти мгновенно подстраивается под профиль температуры (); присутствует химический источник теплоты, распределенный по всему объему канала, заданный реакцией дробного порядка. На границе задаются гидродинамические условия прилипания, тепловые и концентрационные граничные условия первого рода.

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

(1)

где , и .

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

(2)

Для системы уравнений (1) граничные условия задаются в виде

1.

2. (3)

3.

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

На входе в канал профиль компоненты вектора скорости (то есть функция ) является решением задачи

, (4)

при этом , а составляющая градиента давления в направлении оси определяется из выражения

, (5)

где - расход жидкости через поперечное сечение канала.

В качестве реологической модели используется модель Кутателадзе-Хабахпашевой [1] для структурно вязкой жидкости, подробно описанная в [4].

Для приведения системы уравнений (1) к безразмерному виду иcпользуются те же обозначения что и в [3],[4].

Приведя вязкость, как и в [4], к безразмерному виду, получим

, (6)

, (7)

, ,(8)

Подставляя (6) в (1) и переходя к безразмерным параметрам, как и в [4], имеем

(9)

После перехода к безразмерным параметрам граничные условия первого рода (3) примут вид

1. ,

2. , (10)

3. , .

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

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

Результаты

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

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

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

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

а б

Рис.1. Критические поверхности зависимостей и .

Параметру в реологической модели в [3] соответствует параметр из реологической модели в [3], параметру параметр , а параметру параметр . Поэтому на рис. 1 приведены критические поверхности, содержащие именно эти параметры.

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

Как следует из графика на рис 1 б, критические значения параметра не изменяются при изменении параметра .

Литература

1. Кутателадзе С.С., Хабахпашева Е.М., Попов В.И., К гидродинамике жидкостей с переменной вязкостью // ПМТФ. 1966. №1, С. 45-49.

2. Франк-Каменецкий Д.А. Теплопередача и диффузия в химической кинетике // 3-е изд. М.: Наука. 1987. С. 320-358.

3. Абайдуллин Б.Р. Влияние геометрии бесконечного коаксиального канала и реологических характеристик жидкости на критические значения температуры, скорости и вязкости. // Труды Академэнерго №3, 2006. С. 16-27.

4. Абайдуллин Б. Р. Исследование влияния основных определяющих параметров на критический теплообмен при течении неньютоновской жидкости на начальном участке коаксиального канала. // Труды Академэнерго №1, 2007. С. 13-24.

5. Абайдуллин Б.Р. Зависимости характеристик тепломассобмена от параметра Франка Каменецкого при критическом режиме теплообмена неньютоновской жидкости в реакторе. // Труды Академэнерго № 4, 2012. С. 20-28.

Абайдуллин Булат Равилевич – млад. науч. сотр. лаб. ТФИ Исследовательского центра проблем энергетики Федерального бюджетного учреждения науки Казанского научного центра Российской академии наук (г. Казань, 420056, ул. Компрессорщиков, д.31), 8(843) 239-82-83