Математическое описание распределения концентраций реагентов по длине реактора идеального вытеснения в случае последовательной реакции Текст научной статьи по специальности «Математика»

2015 Химическая технология и биотехнология № 4

УДК 66.974.434

Н.Р. Валиуллин, В.П. Плаксина

Пермский национальный исследовательский политехнический университет, Пермь, Россия

МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ РАСПРЕДЕЛЕНИЯ КОНЦЕНТРАЦИЙ РЕАГЕНТОВ ПО ДЛИНЕ РЕАКТОРА ИДЕАЛЬНОГО ВЫТЕСНЕНИЯ В СЛУЧАЕ ПОСЛЕДОВАТЕЛЬНОЙ РЕАКЦИИ

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

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

В связи с этим актуальными являются исследования, связанные с моделированием химических процессов, происходящих в реакторах различных типов. Наибольший интерес представляет собой распределение концентраций по объему реактора при заданных параметрах (констант скоростей, соотношений веществ, вступающих в химическое взаимодействие, температуры и др.), для получения продукта в заданном количестве. Этот метод имеет ряд преимуществ перед остальными: результативность, небольшие экономические затраты, осуществимость.

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

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

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

N.R. Valiullin, V.P. Plaksina

Perm National Research Polytechnic University, Perm, Russian Federation

MATHEMATICAL DESCRIPTION OF REAGENT CONCENTRATION DISTRIBUTION LENGTHWAYS THE IDEAL DISPLACEMENT REACTOR IN A CONSECUTIVE REACTION

One of the problems of the commercial chemico-techological process developmeut is the selection of the type and design of the chemical reactor.

The modeling (simulation) method which allows predicting the course of processes in devices of the specified design of any size and approaching calculation of their optimum hardware and technological parameters from a position of the general engineering principles of chemical technology is applied to studying, calculation and design of chemical reactors.

In this regard the research related to connected with modelling of chemical processes of various types occurring in reactors is timely. Distribution of concentration lengthways parameters (constants of speeds, ratios of the substances entering chemical interaction, temperatures and others), is important for the production of the product in the amount reguired. This method has a number of advantages as compared to other methods: productivity, low economic expenses, feasibility.

In this work the complete mathematical model of the chemical reactor is presented. By means of simplifications the authors have written a differential eguation describing the chemical process proceeding in the ideal displacement reactor working in isothermal and stationary mode sufficiently completely.The solution of this differential eguation in terms of functional and differential eguation theory is shown on a consecutive reaction (in general terms) taken as an example. The research resulted in finding the functional dependence of the length of the reactor on the degree of the initial substance transformation, drawing 2 sched-

ules with the specified values of speed constants and stoichiometric coefficients showing dependence of the intermediate substance transformation degree on the initial substance transformation degree

Keywords: mathematical model, concentration, reactor of ideal replacement, extent of transformation.

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

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

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

Математическая модель химического реактора имеет следующий вид [8]:

/7c

-div (cj • u) + DV2cj - Wj = -J ;

dt

rJ

(1) (2)

div(cJ • u) = u • grad cJ + cJ • div u;

„2 д2сг д2сг д2сг Ч2сг =—Г + —Г + - г

дх2 д12 &2

где сг - молярная концентрация вещества Г в потоке жидкости, протекающей через реактор; и - линейная скорость потока; Б - коэффициент молекулярной диффузии; - скорость химической реакции; I - время; I - длина реактора; Qвх - теплота входящего потока; Qвых - теплота выходящего потока; Qx.fi - теплота химической реакции; Qт.о - теплота, расходуемая на теплообмен с окружающей средой; Qф.пр - теплота фазовых превращений; Qнак - теплота накопления.

Упростим модель (1)-(2). Пусть гидродинамическая обстановка в реакторе соответствует реактору идеального вытеснения, работающему в изотермическом и стационарном режиме. Тогда из системы можно исключить уравнение теплового баланса (2), так как температура в реакторе постоянна. Кроме того, можно исключить из уравнения

(1) диффузионный перенос БУ2сГ исходного вещества Г, ввиду отсутствия в реакторе идеального вытеснения перемешивания в направлении оси потока и в каждом отдельно взятом сечении, перпендикулярном оси потока. Параметры процесса (концентрация, температура и т.д.) полностью выравнены. Частная производная концентрации по времени равна нулю при стационарном режиме работы реактора. Таким образом, система (1)-(2) запишется в виде

Нс

-и^ - ЖгГ = 0. (3)

В данном случае знак частной производной заменяется на выражение полного дифференциала, так как при стационарном режиме концентрация сг и линейная скорость вещества иу будут величинами, зависящими только от расстояния I от начала реакционной зоны реактора. Определим концентрацию вещества Г:

= ПГ , (4)

V

где пг - число молей вещества Г, проходящее через заданное сечение

потока в единицу времени; V - объем вещества, проходящего в единицу времени через заданное сечение.

Линейную скорость можно выразить формулой

V

S

где S - площадь сечения.

Подставляя (4) и (5) в уравнение (3), получим

иу =-, (5)

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

У —— у^В —— у2С,

где к1 и к2 - константы скоростей соответствующих реакций; у1 и у2 -

стехиометрические коэффициенты.

Пусть в начальный момент времени имеется только вещество У. На некотором расстоянии I от начала реакционной зоны реактора число молей вещества У составит п0У (1 - х). Согласно химическому уравнению должно возникнуть у1п0Ух молей вещества В, из которых п0Уу

молей превращается в вещество С. Таким образом, число молей вещества В, проходящего в единицу времени через сечение, находящегося на расстоянии I от начала реакционной зоны реактора, будет п0У (у1 х - у). Кроме того, согласно химическому уравнению образуется

у 2 п,

2"0У

у/ у1 молей вещества С. Суммируя эти величины, получаем

I п -

п0У

1 + (у1 -1) х +

-1

чу1 У

у

(7)

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

р¥ -I пЯТ , (8)

где р - общее давление смеси всех газов; V - общий объем смеси газов; I п - общее число молей всех участников реакции (как исходных, так

и конечных продуктов); Я - универсальная газовая постоянная; Т -температура смеси газов.

Выразим из формулы (8) общий объем смеси газов:

V =

£ пят

(9)

Концентрация любого компонента в смеси газов определяется как отношение числа молей этого компонента к общему объему смеси газов. Подставляя уравнение (9) в уравнение (4), определим концентрацию вещества у:

пгР = пог (1 - х) Р

сг =

£ пят £ пят

(10)

где пг = п0Г (1 - х) - количество вещества пвыраженное в долях от

исходного количества, еще не вступившего в реакцию.

Подставив уравнение (7) в уравнение (10), окончательно получим

сг =

(1 - х)Р

1 + (V! - 1)X + - 1 I у

V.

(11)

ЯТ

Аналогично определяем концентрацию вещества В:

Кх - У) Р

св =

1 + (V -1) х + -1 I у

V.

VV1

(12)

ЯТ

Уравнения (11) и (12) показывают зависимость концентраций веществ Г и В от расстояния до начала реакционной зоны. Давление р также является функцией расстояния I от начала реакционной зоны. Однако, учитывая пренебрежимо малые значения перепада давления по длине реакционной зоны, величину Р будем рассматривать как константу.

Используя выражение (6) в дифференциальном виде, запишем зависимости количества прореагировавшего вещества Г и образующегося вещества С, на расстоянии I от начала реакционной зоны, выраженные в долях от исходного вещества. Будем считать, что вторая реакция имеет первый порядок:

1 ёп/ 8 ё/

= К

(1 - х) Р

V.,

1 ёпС 8

— К2

1 + (V -1) х + -1 I у

У .

(V1х - У) Р

ЯТ

1 + (v1 -1) х +

^ _ 1 VV1 У

у

ЯТ

Так как пГ = пГ (1 - х), пС = v2п0Гу/V!, то

(1 - х) Р

и ёх = К 8 ё/ 1

1 + (V! - 1)х +

^ -1 VV1 У

у

ЯТ

п0Г ¿У = , V1 — 2 8 ё/ V 2

(vlх - У) Р

1 + (v1 -1) х +

^ -1 vVl У

у

ЯТ

Разделим уравнение (16) на уравнение (15), получим

йу К К v1x

—+—У = —— , ёх 1 - х 1 - х

К = ■

V 2

(13)

(14)

(15)

(16)

(17)

(18)

Для решения дифференциального уравнения зададим начальные условия:

йу К К v1x

—+-У = ——.

ёх 1 - х 1 - х

У(0) = 0

Запишем систему в операторном виде:

Г Ь(у (х)) = / (х), 1 У(0) = 0. Ее решение имеет вид [9]

(19)

(20)

у(х) - Ь-1(/(х)) - | С(х, 5)/(5^5. (21)

0

При фиксированном 5 С(х,5) = gs(х) выполняются условия [9]:

1. Ь ^(5)) = 0.

2. С (5, 5) =1.

3. С(х, 5) = 0, если х < 5.

Найдем решение однородного уравнения при заданном краевом условии:

4У К

—+-у - 0, „ ч

ах 1 - х (22)

Iу(5) - 1.

Решив эту систему, получим

У(х) - (1 - х)К (1 - 5)-К . (23)

Так как gS(х) = у(х), тогда функция Коши имеет вид

С(х,5) - (1-х)К (1 - 5)-К . (24)

Найдем решение системы (19):

у(х) -\С(х, 5)/(5)ё5 - 1 (1 - х)К (1 - 5)-К (] ё5 -

0 0 ^ 15 У

К у1(1 - х) у1(1 - х)К

= у1 +-

1 - К 1 - К

Задавая определенные значения констант скоростей к1 и к2 и сте-хиометрических коэффициентов у1 и у2, можно построить графики

(рис. 1, 2), показывающие зависимость степени превращения вещества В от степени превращения вещества У:

к2 - 1 . у -у - 1. К - к2у1 - 1. - ; у1 -У 2 - 1; К - - -— ;

к1 2 к^у 2 2

^ = 2; у1 -у2 - 1; К - - 2. к1 к1у2

й РО и о

в

и м

2

<и

в

а

Он

&

л к

н О

у(рс)

-1

1-

0,8 ■

0,6 ■

0,4 ■

0,2 ■

-1 -0,8 -0,6 -0,4 -0,2 -0,2 " ) 0,2 0,4 0,6 0,8 1

-0,4 ■

-0,6 ■

-0,8 ■

-1-

-1

Степень превращения вещества 3 Рис. 1. Зависимость степени превращения вещества В от степени превращения вещества 3 при К = -2 и — = V 2 = 1

03

& л

аз

н О

у(х)

\

\ 0,8 ■

\ 0,6 ■

\ 0,4 ■

X. 0,2 ■

-1 -0,8 -0,6 -0,4 -0,2 -0,2 " ) 0,2 0,4 0,6 0,8 1

-0,4 ■

-0,6 ■

-0,8 "

-1-

Степень превращения вещества 3

Рис. 2. Зависимость степени превращения вещества В от степени превращения вещества 3 при К = 2 и — = V 2 = 1

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

l = n°jRT j-v2(1 -K)ln(1 - x) - (v, -1 + K -Kv2)x + Sk1p(1 - K ) 1 2 1 2

, ч (1 - x)K , X 1 1 + (V2 --(V2 -V1) K J.

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

1. Мухленова И.П. Общая химическая технология. - Ч. 1. - М.: Высшая школа, 1984. - 256 с.

2. Кафаров В.В., Глебов М.Б. Математическое моделирование основных процессов химических производств. - М.: Высшая школа, 1991. - 367 с.

3. Натараев С.В. Системный анализ и математическое моделирование процессов химической технологии: учеб. пособие / под ред. В.Н. Блинчева; Иван. гос. хим.-технол. ун-т. - Иваново, 2007. - 80 с.

4. Панченков Г.М., Лебедев В.П. Химическая кинетика и катализ: учеб. пособие для вузов. - 3-е изд., испр. и доп. - М.: Химия, 1985. - 592 с.

5. Кондауров Б.П., Александров В.И., Артемов А.В. Общая химическая технология. - М.: Академия, 2005. - 336 с.

6. Бесков В. С. Общая химическая технология. - М.: Академкнига, 2005. - 452 с.

7. Стромберг А.Г., Семченко Д.П. Физическая химия. - М.: Высшая школа, 2006. - 527 с.

8. Беренгартен М.Г., Бондарева Т.И., Кутепов А.М. Общая химическая технология. - М.: Высшая школа, 1990. - 520 с.

9. Азбелев Н.В., Максимов В.П., Рахматуллина Л.Ф. Введение в теорию функционально-дифференциальных уравнений. - М.: Наука, 1991. - 206 с.

References

1. Mukhlenova I.P. Obshchaya khimicheskaya tekhnologiya [General chemical technology]. Moscow: Vysshaya shkola, 1984. 256 p.

2. Kafarov V.V., Glebov M.B. Matematicheskoe modelirovanie osnov-nykh protsessov khimicheskikh proizvodstv [Mathematical modeling of the basic processes of chemical production]. Moscow: Vysshaya shkola, 1991. 367 p.

3. Nataraev S.V. Sistemnyj analiz i matematicheskoe modelirovanie protsessov khimicheskoj tekhnologii [System analysis and mathematical modeling of processes of chemical technology]. Ivanovo, 2007. 80 p.

4. Panchenkov G.M., Lebedev V.P. Khimicheskaya kinetika i kataliz [Chemical kinetics and catalysis]. Moscow: Khimiya, 1985. 592 p.

5. Kondaurov B.P., Aleksandrov V.I., Artemov A.V. Obshchaya khimicheskaya tekhnologiya [General chemical technology]. Moscow: Akademiya, 2005. 336 p.

6. Beskov V.S. Obshchaya khimicheskaya tekhnologiya [General chemical technology]. Moscow: Akademkniga, 2005. 452 p.

7. Stromberg A.G., Semchenko D.P. Fizicheskaya khimiya [Physical chemistry]. Moscow: Vysshaya shkola, 2006. 527 p.

8. Berengarten M.G., Bondareva T.I., Kutepov A.M. Obshchaya khimicheskaya tekhnologiya [General chemical technology]. Moscow: Vysshaja shkola, 1990. 520 p.

9. Azbelev N.V., Maksimov V.P., Rakhmatullina L.F. Vvedenie v teoriyu funktsionalno-differentsialnykh uravnenij [Introduction to the theory of functional differential equations]. Moscow: Nauka, 1991. 206 p.

Получено 01.07.2015

Об авторах

Валиуллин Наиль Рафикович (Пермь, Россия) - студент 2-го курса, кафедра машин и аппаратов производственных процессов Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29; e-mail: [email protected]).

Плаксина Вера Павловна (Пермь, Россия) - кандидат физико-математических наук, доцент кафедры высшей математики Пермского национального исследовательского политехнического университета (614990, г. Пермь, Комсомольский пр., 29; e-mail: [email protected]).

About the authors

Nail R. Valiullin (Perm, Russian Federation) - Student, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation; e-mail: [email protected]).

Vera P. Plaksina (Perm, Russian Federation) - Ph.D. in Physical and Mathematical Sciences, Associate professor, Department of Higher Mathematics of the Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation; e-mail: [email protected]).