Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Моделирование - Лекция

.pdf
Скачиваний:
52
Добавлен:
13.03.2016
Размер:
2 Мб
Скачать
☆

Расчет однократного испарения

Задан состав питания 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