
- •Моделирование физических процессов Учебное пособие
- •Содержание
- •3. 8. Математическая модель массообмена ……………………………. 54
- •5.1.1. Формулы конечного дифференцирования …………………… 74
- •5.6.2. Метод Ньютона ………………………………………………. 95
- •Введение
- •VII этап
- •2. Структурные математические модели
- •Основы структурного моделирования
- •2.2. Типовые элементы
- •2.2.1. Способы математического описания звеньев и систем.
- •2.2.2. Типовые динамические звенья.
- •2.3. Схемы соединения типовых элементов
- •2.4. Структурные модели теплообменного оборудования
- •2.4.1. Составление математических моделей тепловых объектов.
- •3. Физические математические модели
- •3.1 Использование математических моделей
- •Под углом 57 к горизонту со скоростью 20 м/с
- •Изменение температуры при остывании тела
- •3.2 Математическая модель теплопроводности
- •. Математическая модель гидродинамики
- •3.4 Математическая модель теплообменника
- •3.5 Математическая модель парогенератора
- •С естественной циркуляцией
- •Математическая модель горения
- •Математическая модель образования токсичных
- •3. 8. Математическая модель массообмена
- •Дифференциальные уравнения диффузионного пограничного слоя
- •Моделирование теплоотдачи диффузией.
- •3.9. Математическая модель парового котла
- •3.9.1. Математическая модель динамических процессов.
- •3.9.2. Регулирование давления пара и тепловой нагрузки
- •3.10. Пример построения математической модели объекта
- •4. Планирование эксперимента
- •Полный трехфакторный эксперимент
- •Численные методы
- •5.1. Интерполяция
- •5.1.1. Формулы конечного дифференцирования
- •5.1.2. Сплайны
- •Способы получения краевых условий:
- •5.2. Решение систем алгебраических линейных уравнений
- •Метод прогонки. Метод прогонки используется для решения систем специального вида
- •5.3. Решение систем алгебраических нелинейных уравнений
- •5.3.1. Метод итераций
- •5.3.2. Метод Ньютона
- •5.4. Решение дифференциальных уравнений
- •5.4.1. Постановка задачи. Разностные схемы.
- •5.4.2. Разностные схемы. Метод Рунге – Кутта.
- •Пример решения дифференциального уравнения методом Рунге-Кутта.
- •Результаты решения уравнения методом Рунге-Кутта
- •5.5. Решение краевых задач методом прогонки
- •5.5.2. Распространение тепла в пластине радиатора.
- •Сравнение метода прогонки с точным решением
- •5.6. Решение краевых задач методом пристрелки
- •5.6.1. Метод пристрелки.
- •5.6.2. Метод Ньютона. Рассмотрим граничную задачу, определяемую дифференциальным уравнением второго порядка
- •5.7. Решение краевых задач методом конечных разностей
- •5.7.1. Метод конечных разностей.
- •5.7.2. Линейные дифференциальные уравнения.
- •Моделирующие программы
- •Смеси газов в реакторе
- •Перемещение свободной поверхности жидкости
- •Скоростей течения испаряющегося метана
- •Рекомендованный библиографический список
- •Пасько Петр Иванович «моделирование физических процессов»
- •346360, Г. Волгодонск, ул. Ленина 73/94
-
Математическая модель образования токсичных
ВЕЩЕСТВ
Образование и превращения соединений серы. Сера входит в состав мазута и угля в виде сложных сероорганических соединений. При сжигании топлива сера из этих соединений через множество промежуточных реакций может переходить в газовую фазу в виде диоксида и триоксида серы SO2 и SO3, сероводорода Н2S, а также других соединений и радикалов: SO, SH, S2O, S, S2. Доля каждого соединения зависит в основном от избытка воздуха и от температуры. На начальном участке горения топлива, когда воздух еще не успевает перемешаться с топливом, а температура не высока (t < 1100 °С), практически вся сера переходит в Н2S. Далее сероводород окисляется по реакции
2Н2S + 3О2 = 2Н2О + 2SO2
и при нормальном количестве воздуха в продуктах сгорания его не остается. Однако до сгорания он может вызывать коррозию экранов топки.
Из двух оксидов серы, образующихся при ее газификации, наиболее опасен для котла триоксид серы SО3, т.к., соединяясь с парами воды, он дает пары серной кислоты: Н2О + SО3 = Н2SО4. Конденсируясь на поверхностях нагрева в зоне низких температур, она вызывает их сильную коррозию.
Как показывают исследования, количество SО3 в газах после первичной газификации серы может существенно увеличиваться. Есть два пути такого увеличения: окисление SО2 в SО3 в гомогенных реакциях и каталитическое окисление SО2 в SО3 на горячих поверхностях нагрева.
Образование оксидов азота. При сжигании органического топлива в паровых котлах непосредственно в топочной камере образуется практически только монооксид азота (NО = 95 – 100 %, остальное – диоксид азота NО2). Далее по газовому тракту котла происходит доокисление NО в NО2. Этот процесс продолжается далее в атмосфере за счет озона, в итоге в NО2 превращается порядка 40 – 70 % первоначально образовавшегося в топке NО. Механизм образования NОХ в процессе горения включает более ста реакций и до конца еще не изучен.
3. 8. Математическая модель массообмена
Движущей силой массопереноса могут быть различные эффекты, обусловленные разностью концентраций вещества, разной восприимчивостью к действию электрических или магнитных сил, разностью плотностей компонентов смеси и т.д.
Если масса переносится механизмом диффузии, то математическая формулировка такой задачи вполне аналогична задаче теплопроводности.
.
В стационарной постановке такая задача описывается уравнением Лапласа.
При конвективном переносе массы следует решать систему уравнений, которая должна быть дополнена соответствующими граничными условиями. Среди них должно присутствовать уравнение, определяющее источники и стоки массы. Источниками массы могут, в частности, быть химические реакции. Причем здесь необходимо строго отличать гомогенные и гетерогенные (с участием материала стенки) химические реакции. Гомогенные реакции следует учитывать, вводя соответствующие "источниковые" члены непосредственно в дифференциальные уравнения, а гетерогенные реакции следует учитывать с помощью граничных условий.
Задачи массопереноса намного сложнее задач конвективного теплообмена, так как увеличивается число уравнений, а коэффициенты переноса следует принимать зависящими от температуры и концентрации химических компонентов, что делает задачу нелинейной.
Для
массообмена использует следующую
систему обозначений. Масса обозначается
символом
.
Производная по времени обозначается
точкой над символом
.
Производные по координатам обозначаются
одним штрихом
,
производные по поверхности – двумя
штрихами
,
а производные по объему – тремя штрихами
.
Для рассматриваемой величины символ
обозначает массовый расход (точка
сверху) через единицу поверхности (два
штриха) и представляет собой, таким
образом, скорость массопереноса,
считается положительной, если полный
поток массы направлен в рассматриваемую
фазу. Так же вводится понятия
"рассматриваемой" и "соседней"
фаз. Вектор плотности полного потока
массы обозначается
и включает все компоненты смеси.
Вследствие диффузии скорости отдельных
компонентов могут различаться. Вектор
скорости смеси будем обозначать
.
Плотность смеси будем обозначать .
Проекции скорости будут обозначаться
.
– конвективная составляющая плотности
потока j-го
компонента на поверхности раздела фаз
и под ней.
– мощность источника j-го
компонента, образующегося в некоторой
точке пограничного слоя вследствие
химической реакции, кг/(м3
с). Для этой величины выполняется
равенство
.
– массовая концентрация химического
элемента
,
кг элемента
/кг
смеси.
– стехиометрическая постоянная элемента
в соединении j,
кг элемента
/кг
j-го
компонента.
С учетом принятого правила знаков получаем выражение для формального определения энтальпии переносимого вещества
.