Телефон: 8-800-350-22-65
WhatsApp: 8-800-350-22-65
Telegram: sibac
Прием заявок круглосуточно
График работы офиса: с 9.00 до 18.00 Нск (5.00 - 14.00 Мск)

Статья опубликована в рамках: XL Международной научно-практической конференции «Естественные и математические науки в современном мире» (Россия, г. Новосибирск, 07 марта 2016 г.)

Наука: Информационные технологии

Секция: Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей

Скачать книгу(-и): Сборник статей конференции

Библиографическое описание:
Исаев С.М., Тлебаев М.Б. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОТОЧНОГО РЕАКТОРА С МЕШАЛКОЙ (CSTR) // Естественные и математические науки в современном мире: сб. ст. по матер. XL междунар. науч.-практ. конф. № 3(38). – Новосибирск: СибАК, 2016. – С. 32-40.
Проголосовать за статью
Дипломы участников
У данной статьи нет
дипломов

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОТОЧНОГО РЕАКТОРА С МЕШАЛКОЙ (CSTR)

Исаев Саги Махсутулы

д-р техн. наук, проф., Таразский государственный университет имени М.Х. Дулати,

Республика Казахстан, г. Тараз

Тлебаев Манат Бейшенович

д-р техн. наук, проф., Таразский государственный университет имени М.Х. Дулати,

Республика Казахстан, г. Тараз

 

MATHEMATICAL MODEL OF CONTINUOUS STIRRED TANK REACTOR (CSTR)

Sagi Issayev

master’s degree Taraz State University named after M.Kh. Dulati,

Kazakhstan, Taraz

Manat Tlebayev

doctor of engineering sciences, Professor, Taraz State University named after M.Kh. Dulati,

Kazakhstan, Taraz

 

АННОТАЦИЯ

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

ABSTRACT

This paper deals with basic simulation studies on of the common used devices in chemical industry Continuous Stirred Tank Reactor (CSTR). Simulations on mathematical models has several advantages over the experiment on a real model or system. The mathematical model is developed from material balances. Simulation is very important and popular tool now a days, when computation speed of computers increases exponentially every day numerical mathematics is used for steady-state analysis and dynamic analysis.

 

Ключевые слова: проточный реактор с мешалкой; температура; концентрация; интенсивность подачи; жидкость; реактор; энергия.

Keywords: сontinuous stirred tank reactors; temperature; concentration; flow rate; liquid; reactor; energy.

 

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

Рассмотрим и выведем математическую модель проточного реактора с мешалкой (CSTR) и рубашкой. Параметры математической модели реактора с рубашкой следующие:

Fi – интенсивность подачи, расход;

CAi – концентрация подачи;

Ti – температура подачи.

 

Рисунок 1. Проточный реактор с мешалкой

 

Предположим, что подача имеет Fi, CAi, Ti. Индекс i обозначает input (введение) начальное значение, А – начальная концентрация.

Пар, который выпускается в конце, – это конечный пар, который имеет следующие параметры: F, CA, T. До того, как выведем математическую модель, мы должны знать единицы измерения разных потоков, таких как Fi, F, FС , все эти единицы – расходы подачи, или объемная интенсивность подачи. CA, CAi – молярная концентрация, моль/объем [2, с. 45].

Предложим, что:

  1. Идеальное смешивание (Т и СА исходного одинакова) с Т и СА субстрат, который внутри реактора. Идеальная смесь означает, что в каждой точке реактора Т и СА одинаковы, а также Т и СА при выпуске одинаковы с реактором.
  2. r – плотность жидкости; СА – теплоемкость. r, CA – константы.
  3. А®В экзотермическая реакция первой степени. Для того, чтобы убрать экзотермическое тепло в реакторе, была использована смазочно- охлаждающая жидкость.
  4. Реактор идеально изолирован, что означает, что нет потери тепла (не выделяется (теряется) тепло в окружающую среду).
  5. Смазочно-охлаждающая жидкость идеально смешана в реакторе.
  6. Мы не будем рассматривать баланс энергии рубашки.

Это предположения, рассматриваемые в системе. На их основе мы выведем формулы. Для начала рассмотрим закон сохранения массы (massbalance) [3, с. 56]:

V – объем жидкости; r – плотность; Vr – масса.

Дифференцируем массу и получаем следующее дифференциальное уравнение:

®   (1) – закон сохранения массы.

Закон сохранения массы компонента (составной части, компонент А)

Чтобы представить скорость изменения, мы дифференцируем

 – скорость исчезновения компонента А;

заменяем формулой (1), после сокращаем одинаковые члены, затем делим на :

Рассмотрим уравнением Аррениуса:

– скорость реакции;  – фактор частоты (константа перед экспонентой);  – энергия активации;  – универсальная газовая константа. Заменяем  в уравнении уравнением Аррениуса и получаем следующее [4, с. 100]:

                      (2)

– уравнение закона сохранения массы компонента.

Далее рассмотрим закон сохранения энергии [5, с. 65]:

              (3)

Q= – энергетический баланс. Входные переменные:– переменная загрузки; Q,(F) – управляющие переменные, исходные переменные – V,,T. Все три переменные рассмотрим, как накопление, их можно также назвать параметры состояния [6, с. 100].

Среди входных переменных есть манипулируемые переменные. Для примера CSTR системы мы рассмотрим объем жидкости (таблица 1).

Таблица 1.

Объем жидкости

            MV

CV

 F

V

Q

T

 

 

CV – контрольные переменные;

MV – управляющая переменная:

LV – переменная загрузки;

Это разработка простой CSTR системы, здесь мы увидели разные переменные, которые были вовлечены в этот CSTR пример.

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

 

Рисунок 2. Блок-схема

 

Программа на MatLab:

function xdot=cstr1(t,x)

global u

% Входной (1):

% Температура охлаждающей рубашки (К)

Tc = u;

% Состояние (2):

% Концентрация в ПРМ (моль / м ^ 3)

C = x(1,1);

%Температура в ПРМ (К)

T = x(2,1);

% Параметры:

% Объемный расход (м ^ 3 / сек)

Fi = 100;

% Объем ПРМ (м ^ 3)

V = 100;

% Плотность смеси A-B (кг / м ^ 3)

rho = 1000;

% Теплоемкость смеси A-B (Дж / кг-К)

Ca = .239;

% Тепло реакции на А-> В (Дж / моль)

mdelH = 5e4;

% Е – энергия активации в Аррениус Уравнение (Дж / моль)

% R – универсальная газовая постоянная = 8,31451 Дж / моль-К

EoverR = 8750;

% Предварительно экспоненциальный фактор (1 / сек)

k0 = 7.2e10;

% U – общий коэффициент теплопередачи (Вт / м^2-К)

% А – площадь, это значение является специфичным для расчета U (м ^ 2)

UA = 5e4;

% Подача Концентрация (моль / м ^ 3)

Caf = 1;

% Подача температура (K)

Tf = 350;

% Вычислить х точка:

xdot(1,1) = (Fi/V*(Caf - Ca) - k0*exp(-Eover/R*T)*Ca);

xdot(2,1) = (Fi/V*(Tf - T) + mdelH/(rho*Cp)*k0*exp(-Eover/R*T)*Ca + UA/V/rho/Cp*(Tc-T));

global u

% Устойчивое состояние Начальные условия для Состоянии

Ca_ss = 0.87725294608097;

T_ss = 324.475443431599;

x_ss = [Ca_ss;T_ss];

% Устойчивое состояние Начальное условие для контроля

u_ss = 300;                                                        

% Открыть Петля шаг Изменить

u = 290;

% Финал Время (сек)

tf = 5;

[t,x] = ode15s('cstr1',[0 tf],x_ss);

% Разобрать состояние ценности

Ca = x(:,1);

T = x(:,2);

% Участок результаты

figure (1);

plot (t,Ca);

figure (2)

plot (t,T);

 

Результат программы:

 

Рисунок 3. График концентрации

 

Рисунок 4. График температуры

 

В соответствии с проведенными расчетами были получены графики зависимости концентрации подачи субстрата и времени относительно времени подачи субстрата и подогрева.

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

 

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

1.  Брайнес Я.М. Введение в теорию и расчеты химических и нефтехимических реакторов. – 2-е изд., пер. и доп. – М.: Химия, 1976. – 232 с.

2.  Касаткин А.Г. Основные процессы и аппараты химической технологии: Учебник для вузов. – 10-е изд., стереотипное, доработанное. Перепечатано с изд.1973 г. – М.: ООО ТИД «Альянс», 2004. – 753 с., ил.

3.  Крамерс Х., Вестертерп К. Химические реакторы. – М.: Химия, 1967. – 264 с.

4.  Попов Ю.В., Корчагина Т.К., Панчехин В.А. Химические реакторы (теория химических процессов и расчет реакторов): учебное пособие. – Волгоград: ВолгГТУ, 2013. – 240 с.

5.  Самойлов Н.А. Моделирование в химической технологии и расчет реакторов: монография. – Уфа: ООО «Монография», 2005. – 224 с. 

Проголосовать за статью
Дипломы участников
У данной статьи нет
дипломов

Оставить комментарий

Форма обратной связи о взаимодействии с сайтом
CAPTCHA
Этот вопрос задается для того, чтобы выяснить, являетесь ли Вы человеком или представляете из себя автоматическую спам-рассылку.