Моделирование - Лекция
.pdf
Расчет однократного испарения
Задан состав питания zi , степень испарения f и давление P. Рассчитать температуру и состав паровой yi и жидкой xi фаз.
N |
N |
z |
|
|
||
xi |
|
i |
1 |
|||
|
|
|||||
|
|
1 K |
|
1 f |
||
i 1 |
i 1 |
i |
|
|||
|
|
|
||||
и |
K |
i
y x
i i
Задан состав питания zi , температура T и давление P. Рассчитать степень испарения f и состав паровой yi и жидкой xi фаз.
N |
K |
|
1 z |
|
0 |
||
|
|
i |
|
|
|
i |
|
i 1 |
1 K |
i |
1 f |
|
|||
|
|
|
|
|
|
||
Задача сводится к решению нелинейного уравнения с одной неизвестной.
2015 |
Мат. моделирование ХТ |
51 |
Химическое равновесие
Одной из важнейших и традиционных вычислительных задач химии и химической технологии является задача вычисления равновесного состава реагирующей смеси. Актуально для быстрых реакций, в которых смесь, как правило, находится в состоянии термодинамического равновесия.
Равновесная концентрация задает предел, к которому стремятся концентрации компонентов с течением времени.
• Получение водорода для нужд ХТ |
А |
|
|
(паровая конверсия метана) |
B |
• Получение аммиака из водорода и азота |
R |
|
A, B |
|
S |
Очень важно для энергетики и пищевой промышленности.
2015 |
Мат. моделирование ХТ |
52 |
Характерный вид кинетических кривых
|
|
|
2015 |
Мат. моделирование ХТ |
53 |
Равновесный состав смеси
А
B
|
Реактор |
A, B |
R |
|
S
A B R S R S A B
A B R S
Равновесные концентрации компонентов реакционной смеси зависят от внешних условий (например, от температуры или давления).
Химическая реакция протекает не всегда быстро, даже с катализатором. Далеко не всегда конечный состав реагирующей смеси совпадает с равновесным составом. Равновесная концентрация задает предел, к которому стремятся концентрации компонентов с течением времени.
Задача: предсказать равновесный состав смеси при заданных условиях.
2015 |
Мат. моделирование ХТ |
54 |
Расчет равновесного состава идеальной смеси
Константа равновесия является термодинамической величиной и связана с энергией Гиббса:
RT ln K p GT0 HT0 T ST0
Расчет константы равновесия по справочным данным. случай существенно сложнее, хотя данное уравнение справедливо.
|
|
R |
|
S |
|
|
|
|
|
|
|
|
|
R |
|
S |
|
|
A |
|
|
|
|
|
|
|
|||
Для идеальной смеси |
K p |
A |
B |
|
|||||
|
|
|
|
|
|
||||
|
|
|
A |
|
B |
|
|
рав н |
i |
|
|
|
|
|
|
|
|
|
|
Неидеальный по-прежнему
|
|
|
|
|
i |
i |
рав н |
Для идеального газа |
K p ( газ ) Pi i |
|
|
yRvR ySvS |
PvR vS vA vB |
|
равн |
yAvA yBvB |
|||||
|
i |
|
|
2015 |
Мат. моделирование ХТ |
55 |
Расчет равновесного состава неидеальной смеси
RT ln K p GT0 HT0 T ST0
Для неидеальных смесей необходимо заменить концентрацию и давление на активность и фугитивность, соответственно.
|
a |
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
R |
|
S |
|
|
aA |
|
|||
K p |
|
|
|
|
|
|
|
|||||
|
|
|
R |
|
S |
|
|
|
i |
Ai |
|
|
|
aA |
aB |
|
|
|
|||||||
|
|
A |
|
B |
|
|
|
i |
равн |
|||
|
|
|
|
|
|
|
|
|
равн |
|
|
|
Поскольку коэффициент фугитивности показывает отклонение отклонения от идеального газа, то фугитивность (компонента газовой смеси):
f |
i |
P |
|
Py |
i |
|
i i |
i |
|
Коэффициент фугитивности может быть рассчитан из уравнения состояния:
|
|
|
P |
|
v |
|
1 |
|
|
ln |
|
|
|
|
i |
|
|
dP |
|
i |
|
|
|
||||||
|
|
RT |
|
P |
|
||||
|
|
|
|
|
|||||
|
|
|
0 |
|
|
|
|
|
|
2015 |
|
Мат. моделирование ХТ |
56 |
||||||
Пример – паровая конверсия СО
CO H O CO H |
|||
|
ZnO/CuO |
|
|
2 |
|
2 |
2 |
Имеется по 1 молю угарного газа и воды (всего 2 моля), T=500 K.
Константа скорости реакции при данной температуре по справочным данным равняется 148,4.
Константа равновесия через состав смеси и общее давление:
|
|
|
yvR yvS |
PvR vS vA vB |
yCO |
yH |
2 |
148,4 |
||||
K |
|
R |
S |
|
|
2 |
|
|
||||
p |
yvA yvB |
y |
|
y |
|
|
|
|||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
CO |
H |
O |
|
|||||
|
|
|
A |
B |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
Ключевым компонентом |
является CO, |
поэтому |
вводится неизвестная |
|||||||||
величина xCO – прореагировавшее количество СО. Далее мы выражаем количество всех веществ через эту величину.
2015 |
Мат. моделирование ХТ |
57 |
Пример – паровая конверсия СО
|
|
|
|
|
ZnO/CuO |
|
|
CO H2O CO2 H2 |
|||
Частицы |
Начало |
Конец |
Доля в смеси, yi |
|
|
|
|
СО |
1 |
1-х |
(1-x)/2 |
|
|
|
|
H2O |
1 |
1-х |
(1-x)/2 |
|
|
|
|
CO2 |
0 |
х |
x/2 |
H2 |
0 |
х |
x/2 |
|
|
|
|
Всего |
2 |
2 |
1 |
|
|
|
|
yCO yH |
|
|
x / 2 x / 2 |
|
x2 |
148,4 |
||||
y |
y |
|
|
2 |
148,4 1 x / 2 1 x / 2 |
1 x2 |
||||
|
2 |
|
|
|
|
|
|
|
|
|
CO |
|
H |
O |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2015 |
Мат. моделирование ХТ |
58 |
Пример – паровая конверсия СО
CO H O CO H |
||||
|
|
ZnO/CuO |
|
|
2 |
2 |
2 |
||
|
x |
2 |
|
|
|
|
148,4 |
|
|
|
|
2 |
|
|
|
1 x |
|
|
|
0.99326 x2 |
2 x 1 0 |
|
||
Корни квадратного уравнения: x = 0.924 и x = 1.089
Физический смысл имеет только первый корень, т.к. количество превратившегося вещества не может быть больше его исходного количества.
Таким образом, в данных условиях протекания реакции 92,4% угарного газа превращаются в водород.
2015 |
Мат. моделирование ХТ |
59 |
Несколько стехиометрических уравнений
CH |
H O CO 3H |
|
4 |
2 |
2 |
CH |
2H O CO |
4H |
4 |
2 2 |
2 |
Алгоритм расчета равновесного состава в случае нескольких стехиометрических уравнений остается практически неизменным. Единственное отличие заключается в том, что мы вводим для каждого стехиометрического уравнения свою неизвестную (степень превращения).
Задача сводится к решению системы нелинейных алгебраических уравнений. Численное решение системы, иногда достаточно высокого порядка, может быть получено при помощи любого стандартного математического пакета (Excel, MathCad, MathLab и т.д.).
2015 |
Мат. моделирование ХТ |
60 |
