Моделирование химической кинетики. Моделирование молекул, химических процессов и реакций

Цель работы

1. Ознакомиться с методами построения кинетических моделей гетерогенных химических реакций.

2. Составить кинетическую модель гетерогенной химической реакции в соответствии с заданным механизмом.

3. Выбрать численный метод и разработать программу расчета.

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

Кинетика гетерогенных химических реакций

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

1) равноценность всех участков поверхности катализатора и независимость энергии хемосорбции от степени заполнения поверхности различными адсорбентами;

2) неизменность катализатора и независимость его свойств от состава реакционной смеси и ее воздействия на катализатор;

      равновесное распределение энергии.

Формальным аналогом кинетического закона действующих масс для элементарных процессов на твердых поверхностях является закон действующих поверхностей (ЗДП) .

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

Пусть протекает элементарная химическая реакция

При этом все вещества вступают во взаимодействие из адсорбированного состояния. Обозначим z i – долю поверхности, занятойi -м адсорбированным веществом. Тогда, в соответствии с законом действующих поверхностей, скорость необратимой реакции (1.38) можно записать как

,

где W – скорость химической реакции;

k – константа скорости;

–доля поверхности, занятой i -й адсорбированной частицей;

–доля свободной поверхности;

i – стехиометрические коэффициенты стадий;

– изменение числа молей при протекании химической реакции.

Если не все вещества вступают во взаимодействие из адсорбированного состояния, а реагируют непосредственно из газовой фазы, то в более общем виде выражение закона действующих поверхностей записывается следующим образом:

,

где – парциальные давления (концентрации)-веществ, реагирующих из газовой фазы;

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

n ,m – количество веществ, адсорбированных на поверхности катализатора и реагирующих из газовой фазы.

Пример: пусть протекает адсорбция водорода на активном центре катализатора Z с образованием адсорбированного поверхностного соединения ZH 2:

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

Последующие исследования показали существенную ограниченность этих предположений. Тем не менее Хиншельвудом, Швабом, Хоугеном, Ватсоном и другими на их основе получены уравнения, удовлетворительно описывающие кинетический эксперимент в определенном интервале изменения параметров.

Типовая формула кинетического уравнения, соответствующего этим предположениям, имела вид

,

где k – константа скорости;

С i – концентрацияi -го реагента газовой среды;

–константа равновесия стадии адсорбции i -го компонента;

i – стехиометрический коэффициент.

Наиболее общее описание кинетики сложных реакций дано в теории стационарных реакций Хориути – Темкина .

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

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

В теории стационарных реакций вводятся понятия «пробег стадии», «пробег по маршруту» и «скорость реакции по базисному маршруту». Под числом пробегов стадии понимается разность числа актов элементарной реакции в прямом и обратном направлениях. Тогда скорость простой реакции равна числу ее пробегов за единицу времени в единичном реакционном пространстве. Один пробег по маршруту означает, что произошло столько пробегов каждой из стадий, каково ее стехиометрическое число для данного маршрута. В том случае, когда образование молекулы промежуточного вещества в одной из стадий скомпенсировано расходованием этой молекулы в иной стадии, реализуется стационарный режим реакции. Если в ходе этой стадии образуется не конечный продукт, а новое промежуточное вещество, то и оно должно расходоваться в другой стадии. Полная компенсация образования и расходования промежуточных веществ и означает завершение пробега по какому-либо из маршрутов.

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

Условие стационарности элементарных стадий химических реакций можно записать следующим образом:

где
– скорости элементарных стадий (s -й, прямой и обратной);

–скорость по маршруту Р ;

s -й стадии, маршрутаР .

На основании уравнения (1.44) получим уравнение, которое называется уравнением стационарных реакций :

где
,
... – скорости по маршрутам;

, . .. – скорости элементарных стадий в прямом и обратном направлениях;

–стехиометрический коэффициент i- й стадии поj -му маршруту.

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

Исследование кинетики гетерогенной химической реакции

Рассмотрим пример сложной гетерогенной химической реакции – реакции гидрокрекинга толуола.

Исходными данными являются:

    детальный механизм гетерогенной химической реакции, где Z – активные центры на поверхности катализатора; Z H 2 и т. д. – адсорбированные промежуточные соединения:

1. H 2 + ZZH 2

2. ZH 2 + C 7 H 8 ZC 7 H 8 ∙H 2

3. ZC 7 H 8 ∙H 2  Z + C 6 H 6 +CH 4

C 7 H 8 + H 2  CH 4 + C 6 H 6

    начальные концентрации веществ и константы скоростей, которые равны

С Н2 (0) = 0,6 мольн. доли;

С С7Н8 (0) = 0,4 мольн. доли;

;

;
;

;

;
.

Запишем скорости элементарных стадий механизма по закону действующих поверхностей:

;
;

;

; .

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

При решении системы дифференциальных уравнений (1.48) можно использовать численные методы Эйлера и Рунге-Кутта. Примеры результатов расчетов кинетики гетерогенной химической реакции приведены на рис. 1.5, 1.6. Программа расчета кинетики химических реакций приведена в Приложении Б.

Цель работы

  • 1. Ознакомиться с методами построения кинетических моделей гетерогенных химических реакций.
  • 2. Составить кинетическую модель гетерогенной химической реакции в соответствии с заданным механизмом.
  • 3. Выбрать численный метод и разработать программу расчета.
  • 4. Исследовать динамику изменения концентраций реагирующих веществ реакции и промежуточных соединений.

Кинетика гетерогенных химических реакций

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

  • 1) равноценность всех участков поверхности катализатора и независимость энергии хемосорбции от степени заполнения поверхности различными адсорбентами;
  • 2) неизменность катализатора и независимость его свойств от состава реакционной смеси и ее воздействия на катализатор;
  • 3) равновесное распределение энергии.

Формальным аналогом кинетического закона действующих масс для элементарных процессов на твердых поверхностях является закон действующих поверхностей (ЗДП) .

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

Пусть протекает элементарная химическая реакция

При этом все вещества вступают во взаимодействие из адсорбированного состояния. Обозначим zi - долю поверхности, занятой i-м адсорбированным веществом. Тогда, в соответствии с законом действующих поверхностей, скорость необратимой реакции (1.38) можно записать как

где W - скорость химической реакции;

k - константа скорости; - доля поверхности, занятой i-й адсорбированной частицей; - доля свободной поверхности; i - стехиометрические коэффициенты стадий;

Изменение числа молей при протекании химической реакции.

Если не все вещества вступают во взаимодействие из адсорбированного состояния, а реагируют непосредственно из газовой фазы, то в более общем виде выражение закона действующих поверхностей записывается следующим образом:

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

n, m - количество веществ, адсорбированных на поверхности катализатора и реагирующих из газовой фазы.

Пример: пусть протекает адсорбция водорода на активном центре катализатора Z с образованием адсорбированного поверхностного соединения ZH2:

тогда на основании ЗДП скорость данной элементарной химической реакции можно записать как

В качестве основного фактора, определяющего кинетические зависимости, вначале рассматривался фактор вытеснения, «борьбы» компонентов реакционной смеси за места на поверхности катализатора. При этом принималось дополнительное предположение о высокой скорости адсорбционных и десорбционных стадий по сравнению с собственно химическими превращениями. Последующие исследования показали существенную ограниченность этих предположений. Тем не менее Хиншельвудом, Швабом, Хоугеном, Ватсоном и другими на их основе получены уравнения, удовлетворительно описывающие кинетический эксперимент в определенном интервале изменения параметров. Типовая формула кинетического уравнения

где k - константа скорости;

Сi - концентрация i-го реагента газовой среды;

Константа равновесия стадии адсорбции i-го компонента;

i - стехиометрический коэффициент.

Наиболее общее описание кинетики сложных реакций дано в теории стационарных реакций Хориути - Темкина .

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

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

В теории стационарных реакций вводятся понятия «пробег стадии», «пробег по маршруту» и «скорость реакции по базисному маршруту». Под числом пробегов стадии понимается разность числа актов элементарной реакции в прямом и обратном направлениях. Тогда скорость простой реакции равна числу ее пробегов за единицу времени в единичном реакционном пространстве. Один пробег по маршруту означает, что произошло столько пробегов каждой из стадий, каково ее стехиометрическое число для данного маршрута. В том случае, когда образование молекулы промежуточного вещества в одной из стадий скомпенсировано расходованием этой молекулы в иной стадии, реализуется стационарный режим реакции. Если в ходе этой стадии образуется не конечный продукт, а новое промежуточное вещество, то и оно должно расходоваться в другой стадии. Полная компенсация образования и расходования промежуточных веществ и означает завершение пробега по какому-либо из маршрутов.

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

Условие стационарности элементарных стадий химических реакций можно записать следующим образом:

где - скорости элементарных стадий (s-й, прямой и обратной); - скорость по маршруту Р;

Стехиометрический коэффициент s-й стадии, маршрута Р.

На основании уравнения (1.44) получим уравнение, которое называется уравнением стационарных реакций:

где, ... - скорости по маршрутам;

Скорости элементарных стадий в прямом и обратном направлениях;

Стехиометрический коэффициент i-й стадии по j-му маршруту.

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

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

Особенности кинетики гомогенных и гетерогенных реакций

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

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

Реакция окисления углерода в системе шлак-металл, как пример гетерогенной реакции

Реакция

Примером может служить реакция окисления углерода в системе шлак – металл применительно к ванне мартеновской или электросталеплавильной печи

Три стадии реакции

Здесь можно выделить, по меньшей мере, три стадии:

  1. диффузия кислорода из шлака в металл к месту реакции (граница раздела: металл – газовый пузырь, незаполненные поры на подине или поверхности кусков руды и извести);
  2. химическая реакция между кислородом и углеродом металла на границе раздела упомянутых фаз;
  3. выделение газообразного продукта реакции из металла.

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

Описание диффузии и массопереноса

Диффузия

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

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

а переходя к бесконечно малым приращениям и скорости диффузии (потоку массы через единичную площадь)

получаем уравнение

(3.57) описывающее стационарную диффузию и называемое первым законом Фика.

Диффузия системы с распределёнными параметрами позакон Фика

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

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

Условия применимости молекулярной диффузии

Необходимо подчеркнуть, что уравнения (3.57) и (3.58) относятся к молекулярному переносу в неподвижной среде и справедливы для изотермического процесса и случая, когда диффузия данного компонента не зависит от диффузии других компонентов.

Формула Стокса-Эйнштейна

В этих условиях зависимость коэффициента диффузии от температуры , вязкости среды и радиуса диффундирующих молекул определяется формулой Стокса – Эйнштейна:

(3.59) где

И – газовая постоянная и число Авогадро.

Турбулентная диффузия

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

Например, значение коэффициента атомарной диффузии в неподвижном расплавленном железе при 1500 – 1600°C составляет – . Величина же коэффициента турбулентной диффузии в мартеновской ванне, зависящая от скорости обезуглероживания, составляет 0,0025 -0,0082 , а в конвертерном процессе 2,0 -2,5, т. е. на три порядка выше.

Диффузия с учетом влияния конвекции

С учетом влияния конвекции уравнение диффузии принимает следующий вид:

(3.60) где – скорость переноса вещества, м/с.

Чаще, в случаях преобладающего влияния турбулентной диффузии, используется эмпирическое уравнение вида

– диффузионный поток;

– разность концентраций;

– коэффициент массоотдачи (турбулентной диффузии).

Эмпирическое уравнение для турбулентной диффузии

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

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

Число Рейнольдса

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

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

При значении <2300 поток является ламинарным. В этой области для описания диффузии могут использоваться уравнения (3.57) или (3.60). Область значений 2300<<10000 является переходной. Здесь, в зависимости от степени развития турбулентности и наличия ламинарного слоя, целесообразно использовать уравнения (3.60) или (3.61).

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

Уравнения кинетики гомогенных реакций

Скорость реакции

Скорость реакции представляет собой производную от концентрации по времени

Молекулярность реакции

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

В прямом направлении пропорциональна концентрациям реагирующих веществ и представляется уравнением

(3.63) где

– константа скорости, имеющая смысл при .

Порядок реакции

Определение

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

Остановимся на виде кинетических уравнений в зависимости от порядка реакции.

Реакция нулевого порядка

При реакциях нулевого порядка скорость постоянна во времени

(3.64) После интегрирования получаем

– постоянная интегрирования, имеющая смысл начальной концентрации при =0.

Таким образом, в рассматриваемом случае концентрация реагента линейно убывает во времени.

Реакция первого порядка

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

Кинетическое уравнение имеет вид:

(3.65) а его решение

показывает, что концентрация исходного компонента экспоненциально уменьшается во времени (рис. 3.2).

Рис. 3.2 Изменение концентрации и ее логарифма во времени при реакциях первого порядка

Решение этого уравнения можно представить и в другом виде, более удобном для определения константы скорости реакции. В результате разделения переменных и выбора пределов интегрирования

При температуре

получаем решение

из которого можно выделить, что линейно зависит от времени. Если опытные данные укладываются на прямую линию (см. рис. 3.2), то это указывает на первый порядок реакции. По углу наклона прямой определяется величина .

Реакция второго порядка

Схема реакции второго порядка имеет вид

Или, например,

А скорость реакции описывается уравнением

(3.66) которое при одинаковых концентрациях принимает вид

После разделения переменных и интегрирования соотношения

получаем соотношение

(3.67) которое может быть использовано для определения . Если начальные концентрации реагентов неодинаковы и равны соответственно и, а концентрация продукта в момент составляет , то получаем уравнение

Логарифмирование которого дает

(3.68)

Обратная реакция

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

(3.69) По мере расходования реагентов и и образования продукта скорость прямой реакции уменьшается, а увеличивается. При суммарная скорость равна нулю, наступает равновесие. Тогда

или

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

Влияние температуры на скорость химической реакции

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

– энергия активации.

После интегрирования, при условии, что , получаем

– постоянная, имеющая смысл логарифма константы скорости при бесконечной температуре ().

Это соотношение можно представить также в виде

(3.73)

Энергия активации

Величину можно определить по тангенсу угла наклона прямой (3.72), построенной в координатах , для чего необходимо измерить константы скорости при различных температурах.

Физический смысл энергии активации и механизм химических реакций можно объяснить на основе теории активных столкновений.

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

Рис. 3.3 К вопросу об энергии активации

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

Более подробное изложение молекулярной кинетики, нашедшее дальнейшее развитие в теории переходного состояния, выходит за рамки данного пособия.

Взаимосвязь массопереноса и кинетики в гетерогенных реакциях

Наглядное представление о взаимосвязях массопереноса и кинетики в гетерогенных процессах дает приведенная на рис.3.4 схема обобщенной модели.

Поверхности раздела фаз

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

Рис. 3.6 Диффузионный пограничный слой

  1. – твердое вещество
  2. – диффузионный пограничный слой
  3. – жидкость

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

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

Стадии гетерогенных реакций

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

  1. диффузия частиц исходных веществ к поверхности раздела фаз (реакционной зоне);
  2. адсорбция реагентов на поверхности;
  3. химическая реакция на поверхности;
  4. десорбция продуктов реакции на поверхности раздела фаз;
  5. диффузия этих продуктов от реакционной зоны вглубь одной из фаз.

Стадии 1 и 5 относятся к диффузионным, а стадии 2 – 4 к кинетическим.

Кинетическое сопротивление гетерогенной реакции

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

(3.74) где

– константа скорости суммарного (наблюдаемого) процесса;

– константа скорости кинетической стадии;

– константа скорости (коэффициент диффузии) диффузионной стадии.

Стадия, имеющая наибольшее сопротивление, является лимитирующей.

Особенности процессов в кинетической области

Рассмотрим основные особенности процессов в кинетической области:

Первые три особенности могут наблюдаться и в случае нахождения процесса в переходной области. Четвертый признак является основным экспериментальным подтверждением нахождения процесса в кинетической области.

Особенности процессов в диффузионной области

Основные особенности процессов в диффузионной области:

  1. процесс первого порядка;
  2. слабая зависимость скорости процесса от температуры и от величины поверхности раздела фаз;
  3. 3) резкое влияние на скорость процесса гидродинамических и аэродинамических условий процесса.

Наиболее важным признаком нахождения процесса в диффузионной области являются первая и третья особенности.

Растворение извести как пример гетерогенного процесса

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

Теоретическая часть

Цель работы

Ознакомиться с методами построения кинетических моделей гомогенных химических реакций.

Построить кинетическую модель заданной химической реакции.

Исследовать влияние температуры на выход продуктов и степень превращения.

кинетический химический реакция гомогенный

Кинетика гомогенных химических реакций

Скорость химической реакции есть изменение числа молей реагентов в результате химического взаимодействия в единицу времени в единице объема

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

Константа скорости химической реакции является функцией температуры, и зависимость её от температуры выражается законом Аррениуса:

Порядок выполнения работы

Составить в соответствии с заданной схемой реакций на основании закона действующих масс кинетическую модель.

Разработать алгоритм расчёта составленной кинетической модели с использованием численного метода Эйлера.

Разработать программу расчета кинетики с учетом температурной зависимости констант скорости гомогенной химической реакции.

Обсудить результаты. Сделать выводы по работе.

Практическая часть

Задание №11

А → В ↔ С

Сa0=0,8 моль/л; Cb0=Cc0=0моль/л;

h1=20; K10=0.3; K20=0.25; K30=0.18; E1=10; E2=8; E3=9; T1=393; T2=453; R=8.31;lb2;mas=array of real;tempk=10.0; h=1.0; Ca0=0.8; Cb0=0; Cc0=0; K10=0.3; K20=0.25; K30=0.18; E1=10; E2=8; E3=9; T1=393; T2=453; R=8.31; h1=20;C, f: mas; i:integer; T, k1, k2, k3, temp:real;:=T1;T<=T2 do begin:=k10*exp(-E1/(R*T));:=k20*exp(-E2/(R*T));:=k30*exp(-E3/(R*T));("T=", T:5:2, " ",

"k1=", k1:5:5, " ",

"k2=", k2:5:5, " ",

"k3=", k3:5:5, " ");:=0; C:=Ca0; C:=Cb0; C:=Cc0;temp<=tempk do begin:=-k1*C;:=K1*C-K2*C+K3*c;:=K2*C-K3*C;i:=1 to 3 do C[i]:=C[i]+h*f[i];("C=", C:5:2, " ",

"C=", C:5:2, " ",

"C=", C:5:2, " ",

"temp=", temp:5:2);:=temp+h; end;

T:=T+h1; end;.

Температура 393 К

k1=0,29908k2=0,24939k3=0,17950Концентрация С1ВремяКонцентрация С2ВремяКонцентрация С3Время0,5600,240000,3910,3510,0610,2820,3920,1420,1930,430,2130,1440,3940,2740,0950,3950,3250,0760,3760,3660,0570,3770,3970,0380,3680,4180,0290,3590,4390,02100,35100,4410

Температура 413 К

k1=0.29913 k2=0.24942 k3=0.17953Концентрация С1ВремяКонцентрация С2ВремяКонцентрация С3Время0.56 00.24 00.00 00.39 10.35 10.06 10.28 20.39 20.14 20.19 30.40 30.21 30.14 40.39 40.27 40.09 50.39 50.32 50.07 60.37 60.36 60.05 70.37 70.39 70.03 80.36 80.41 80.02 90.35 90.43 90.02 100.35 100.44 10

k1=0.29917 k2=0.24944 k3=0.17955 Концентрация С1ВремяКонцентрация С2ВремяКонцентрация С3Время0.56 00.24 00.00 00.39 10.35 10.06 10.28 20.39 20.14 20.19 30.40 30.21 30.14 40.39 40.27 40.09 50.39 50.32 50.07 60.37 60.36 60.05 70.37 70.39 70.03 80.36 80.44 80.02 90.35 90.43 90.02 100.35 100.41 10

Температура 453 К

k1=0.29920 k2=0.24947 k3=0.17957 Концентрация С1ВремяКонцентрация С2ВремяКонцентрация С3Время0.56 00.24 00.00 00.39 10.35 10.06 10.28 20.39 20.14 20.19 30.40 30.21 30.14 40.39 40.27 40.09 50.39 50.32 50.07 60.37 60.36 60.05 70.37 70.39 70.02 80.35 80.41 80.02 90.35 90.44 90.03 100.36 100.43 10

Вывод

В ходе лабораторной работы, составили кинетические модели гомогенной химической реакции; построили кинетическую модель химической гомогенной реакции; рассчитали концентрации реагирующего, промежуточных веществ и продукта реакции при различных температурах и времени протекания реакции. Построили зависимости изменения концентрации от времени, при увеличении времени протекания реакции концентрация исходного вещества уменьшается, промежуточных веществ увеличиваются, а концентрация продукта реакции - возрастает. Рекомендации по условиям проведения химической реакции с целью получения максимального выхода целевого продукта С: температура 393 К и время 10 с.

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

Моделирование кинетики гомогенных химических реакций

Этапы развития химической кинетики

Химическая кинетика как наука о скоростях химических реакций начала формироваться в 50-70 гг. XIX в.

В 1862-1867 гг. норвежские ученые Гульдберг и Вааге дали начальную формулировку закона действующих масс : при протекании химической реакции:

Концентрация i-го вещества, ;

Константа скорости;

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

Химическая кинетика в полной мере была сформулирована в работах Вант-Гоффа и Аррениуса в 80-х гг. XIX в.; был разъяснен смысл порядков реакций и введено понятие энергии активации. Вант-Гофф ввел понятия моно-, би- и полимолекулярных реакций:

где n - порядок реакции.

Вант-Гофф и Аррениус, развивший его идеи, утверждали, что температура не есть причина реакции, температура - причина изменения скорости реакции:

(Аррениус в 1889 г.), (1.3)

где А - предэкспоненциальный множитель;

Е - энергия активации;

R - газовая постоянная;

Т - температура.

С 1890 г. величина энергии активации стала универсальной мерой реакционной способности превращающихся веществ. Таким образом, в период 1860-1910 гг. была создана формальная кинетика. Ясность и немногочисленность основных постулатов отличают химическую кинетику периода Вант-Гоффа и Аррениуса.

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

Однако для химика и до настоящего времени остаются наиболее важными две концепции:

Закон действующих масс как закон простой реакции.

Сложность механизма химической реакции.

Основные понятия химической кинетики

Кинетика гомогенных химических реакций

Скорость химической реакции есть изменение числа молей реагентов в результате химического взаимодействия в единицу времени в единице объема (для гомогенных реакций) или на единице поверхности (для гетерогенных процессов) :

где W - скорость химической реакции, ;

V - объем, м3;

N - число молей;

t - время, с.Согласно уравнению (1.4), вводя концентрацию, получим

где С - концентрация, моль/м3,

Для реакций, идущих при постоянном объеме, второе слагаемое в уравнении (1.6) равно нулю и, следовательно,

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

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

где W - скорость химической реакции;

Константа скорости;

Концентрации исходных веществ, ;

Соответствующие стехиометрические коэффициенты в брутто-уравнении химической реакции.

Уравнение (1.8) справедливо для элементарных реакций. Для сложных реакций показатели степени в уравнении (1.8) называются порядками реакции и могут принимать не только целочисленные значения.

Константа скорости химической реакции является функцией температуры, и зависимость от температуры выражается законом Аррениуса:

где - предэкспоненциальный множитель;

E - энергия активации, ;

Т - температура, К;

R - газовая постоянная, .

Рассмотрим гомогенную реакцию

где a, b, c, d - стехиометрические коэффициенты.

Согласно закону действующих масс (1.8) скорость этой реакции запишется следующим образом:

Между скоростями реакции по отдельным компонентам (обозначим их WA, WB, WC, WD) и общей скоростью реакции W существует зависимость

Отсюда вытекают следующие выражения:

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

Кинетические уравнения

Кинетические уравнения связывают скорость реакции с параметрами, от которых она зависит. Наиболее важными из этих параметров являются концентрация, температура, давление, активность катализатора.

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

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

где Wi - скорость химической реакции, моль/м3с;

Ni - текущее количество i-го компонента реакционной смеси, моль;

V - объем реакционной смеси или слоя катализатора (объем реактора), м3;

S - поверхность катализатора, м2;

0 - удельная поверхность катализатора, м2/м3;

t - время, с.

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

где ni - мольный расход i-го компонента реакционной смеси, моль/с;

Скорость подачи реакционной смеси, м3/с;

Время контакта, с.

Для реактора непрерывного действия полного смешения, при установившемся режиме,

где ni0 - начальное количество i-го компонента реакционной смеси, моль/с.

На практике обычно измеряют скорость изменения мольной концентрации Сi (моль/м3; моль/л).

Для реактора периодического действия

Для реактора непрерывного действия

где - объемная скорость подачи реакционной смеси, м3/с.

Если реакция не сопровождается изменением объема, то для реактора идеального вытеснения

Для реактора непрерывного перемешивания

где xi - степень превращения, ;

Среднее время пребывания, = V/, с.

Методы решения кинетических уравнений

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

Для решения дифференциальных уравнений - интегрирования - применяются:

· табличный метод (с использованием таблиц интегралов) - применяется для простейших дифференциальных уравнений;

· аналитические методы применяются для решения дифференциальных уравнений первого порядка;

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

Численные методы

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

Например, уравнение скорости химической реакции описывается уравнением

где СА - концентрация вещества, моль / л;

Время, с.

При малых t можно приближённо принять, что

величину называют шагом интегрирования. Решая уравнение (1.23), получим общую формулу Эйлера

где - правая часть дифференциального уравнения (например,

Задав начальные условия: при t = 0 С = С0, величину шага интегрирования h, а также параметры уравнения, с помощью формулы (1.24) можно провести пошаговый расчёт и получить решение данного уравнения (рис. 1.1).

Рис. 1.1. Графическая иллюстрация метода Эйлера

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

Величина шага интегрирования выбирается, исходя из достижения минимального времени счёта и наименьшей ошибки вычислений.

Общие представления одношаговых методов решения обыкновенных дифференциальных уравнений

Пусть имеется дифференциальное уравнение

удовлетворяющее начальному условию

Требуется найти решение задачи (1.25), (1.26) на отрезке . Выполним разбиение отрезка точками

Этот набор точек называют сеткой, а точки xi (i = i, n) - узлами сетки.

Одношаговые численные методы дают приближения yn к значениям точного решения y(xn) в каждом узле сетки xn на основе известного приближения yn-1 к решению в предыдущем узле xn-1. В общем виде их можно представить так :

Для явных одношаговых методов функция F не зависит от yn+1.

Обозначая

явные одношаговые методы будем записывать также в виде

Явные методы типа Рунге-Кутта

Идея данного метода основана на вычислении приближённого решения y1 в узле x0 + h в виде линейной комбинации с постоянными коэффициентами :

Числа выбираются так, чтобы разложение выражения (1.29) по степеням h совпадало с разложением в ряд Тейлора:

Это эквивалентно следующему. Если ввести вспомогательную функцию

то ее разложение по степеням h должно начинаться с максимально возможной степени:

Если можно определить эти постоянные так, чтобы разложение имело вид (1.32), то говорят, что формула (1.29) с выбранными коэффициентами имеет порядок точности s.

Величина

называется погрешностью метода на шаге, или локальной погрешностью метода, а первое слагаемое в выражении (1.32)

называется главным членом локальной погрешности метода.

Доказано, что если q = 1, 2, 3, 4, то всегда можно выбрать коэффициенты так, чтобы получить метод типа Рунге-Кутта порядка точности q. При q = 5 невозможно построить метод типа Рунге-Кутта (1.29) пятого порядка точности, необходимо брать в комбинации (1.29) более пяти членов.

Исследование кинетики гомогенных химических реакций

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

Пусть протекают химические реакции

На основании закона действующих масс запишем уравнения скоростей химических реакций и составим кинетическую модель:

где СА, СВ, СС, СD - концентрации веществ, моль/л;

ki - константа скорости i-й химической реакции первого порядка, с-1; (для реакций второго порядка размерность константы; для реакций третьего порядка размерность константы);

Wi - скорость i-й химической реакции, моль/лс; t - время реакции, с.

Размещено на http://www.allbest.ru/

Рис. 1.2. Блок-схема расчета кинетики гомогенной химической реакции методом Эйлера

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

Блок-схема расчета кинетики гомогенной химической реакции методом Эйлера приведена на рис. 1.2.

Примеры программ расчета кинетики гомогенных химических реакций приведены в Приложении А. Результаты исследования на математической модели (1.35) влияния температуры на степень превращения исходного реагента и на концентрацию веществ представлены на рис. 1.3, 1.4.

Полученные результаты позволяют сделать вывод об оптимальном времени проведения процесса с целью получения целевого продукта. Математическая модель (1.35) также позволяет исследовать влияние состава сырья на выход продуктов реакции.

Необходимо учитывать, что скорость химической реакции зависит от температуры, поэтому, чтобы использовать кинетическую модель (1.35) для исследования процесса при различных температурах, необходимо ввести зависимость константы скорости химической реакции от температуры по уравнению Аррениуса (1.9).

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

где ai - коэффициенты Рунге-Кутта, которые рассчитываются по следующим формулам:

Литература

кинетика химический гомогенный

1. Панченков Г. М., Лебедев В. П. Химическая кинетика и катализ. - М.: Химия, 1985. - 589 с.

2. Яблонский Г. С., Быков В. И., Горбань А. И. Кинетические модели каталитических реакций. - Новосибирск: Наука, 1983. - 254 с.

3. Кафаров В. В. Методы кибернетики в химии и химической технологии. - М.: Химия, 1988. - 489 с.

4. Кравцов А. В., Новиков А. А., Коваль П. И. Методы анализа химико- технологических процессов. - Томск: изд-во ТПУ, 1994. - 76 с.

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

6. Мойзес О. Е., Коваль П. И., Баженов Д. А., Кузьменко Е. А. Информатика: учеб. пособие. В 2-х ч. - Томск, 1999. - 150 с.

7. Турчак Л. И. Основы численных методов. - М.: Наука, 1987. - 320 с.

8. Офицеров Д. В., Старых В. А. Программирование в интегрированной среде Турбо-Паскаль. - Минск: Беларусь, 1992. - 240 с.

9. Бесков В. С., Флор К. В. Моделирование каталитических процессов и реакторов. - М.: Химия, 1991. - 252 с.

10. Руд Р., Праустниц Дж., Шервуд Т. Свойства газов и жидкостей

/ под ред. Б. И. Соколова. - Л.: Химия, 1982. - 591 с.

11. Танатаров М. А. и др. Технологические расчеты установок переработки нефти. - М.: Химия, 1987. - 350 с.

12. Жоров Ю. М. Термодинамика химических процессов. - М.: Химия, 1985

13. Расчеты основных процессов и аппаратов нефтепереработки: справочник / под ред. Е. Н. Судакова. - М.: Химия, 1979. - 568 с.

14. Кафаров В. В. Разделение многокомпонентных систем в химической технологии. Методы расчета. - М.: Московский химико-технологический институт, 1987. - 84 с.

Размещено на Allbest.ru

Подобные документы

    Методы построения кинетических моделей гомогенных химических реакций. Исследование влияния температуры на выход продуктов и степень превращения. Рекомендации по условиям проведения реакций с целью получения максимального выхода целевых продуктов.

    лабораторная работа , добавлен 19.12.2016

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

    реферат , добавлен 16.11.2009

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

    презентация , добавлен 30.04.2012

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

    презентация , добавлен 11.08.2013

    Рассмотрение превращения энергии (выделение, поглощение), тепловых эффектов, скорости протекания химических гомогенных и гетерогенных реакций. Определение зависимости скорости взаимодействия веществ (молекул, ионов) от их концентрации и температуры.

    реферат , добавлен 27.02.2010

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

    презентация , добавлен 25.10.2014

    Термодинамика и кинетика сложных химических реакций. Фазовые превращения в двухкомпонентной системе "BaO-TiO2". Классификация химических реакций. Диаграммы состояния двухкомпонентных равновесных систем. Методы Вант Гоффа и подбора кинетического уравнения.

    курсовая работа , добавлен 19.05.2014

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

    контрольная работа , добавлен 04.04.2009

    Общее понятие о химической реакции, ее сущность, признаки и условия проведения. Структура химических уравнений, их особенности и отличия от математических уравнений. Классификация и виды химических реакций: соединения, разложения, обмена, замещения.

    реферат , добавлен 25.07.2010

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