§ 4.3. Уравнение материального баланса для элементарного объема проточного химического реактора Прежде чем рассмотреть отдельные типы химических реакторов в соответствии с приведенной классификацией, составим уравнение материального баланса по произвольному участнику реакции – веществу J – для элементарного объема произвольного проточного химического реактора и элементарного промежутка времени.
Рассмотрим поток жидкости, протекающей через реактор. О ходе химического процесса в реакторе будем судить по изменению молярной концентрации вещества J в жидкости сJ. Так как, в общем случае, в реакторе имеет место то, или иное распределение концентрации сJ по объему, а в каждой произвольно выбранной точке еще и распределение концентрации во времени, то считается, что сJ является функцией четырех переменных: трех пространственных координат х, у, z и времени τ: сJ = сJ(х, у, z, τ).
В соответствии с рассмотренными в § 4.2 требованиями к размеру элементарного объема и значению элементарного промежутка времени выберем в качестве элементарного промежутка времени бесконечно малый интервал dτ (dτ 0), а в качестве элементарного промежутка пространства – параллелепипед с бесконечно малыми сторонами dx, dy и dz и объемом dV = dxdydz (рис. 4.1).
В уравнении материального баланса по веществу J должны быть отражены, как указано выше, все изменения, которые произойдут за время dτ с веществом J при прохождении потоком элементарного объема. Эти изменения могут быть связаны с тремя причинами: конвективным переносом, диффузионным переносом и химической реакцией.
Рис. 4.1. Элементарный объем химического реактора:
1, 1/– конвективный и 2, 2/– диффузионный потоки соответственно на входе в элементарный объем и на выходе из него Конвективный перенос, или перенос импульса, вызван движением потока со скоростью u в результате какого-либо внешнего воздействия (например, из-за перепада давления, созданного насосом или компрессором). При макроскопическом движении жидкости каждый данный ее участок передвигается как целое с неизменным составом, и в результате происходит чисто механическое перемешивание: хотя состав каждого передвигающегося участка жидкости может оставаться неизменным (если нет химической реакции) в каждой неподвижной точке пространства (неподвижном элементарном объеме), концентрация жидкости будет со временем меняться. Охарактеризовать конвективный перенос можно изменением импульса единицы объема жидкости с, и.
Диффузионный перенос вызван наличием неравномерного распределения вещества J в пространстве. Вследствие выравнивания концентрации молекулярным переносом веществ реакционной смеси из одного участка жидкости в другой также происходит изменение состава внутри элементарного объема. Охарактеризовать диффузионный перенос можно в соответствии с законами Фика изменением диффузионного потока вещества J, равного D grad cJ (D – коэффициент диффузии).
Протекание химической реакции в элементарном объеме – неотъемлемая часть любого химического процесса. Расход или образование вещества J в ходе химической реакции пропорционален скорости реакции wrJ.
Алгебраическая сумма всех этих трех изменений должна быть равна накоплению (положительному или отрицательному) вещества J в элементарном объеме, т. е. изменению количества вещества J, находящегося внутри элементарного объема, за тот промежуток времени, для которого составляется материальный баланс.
Запишем теперь отдельные составляющие уравнения материального баланса.
Количество вещества, попадающее за время dτ в элементарный объем с конвективным потоком, можно рассматривать как сумму составляющих потока, которые войдут через отдельные грани параллелепипеда. В направлении оси z через грань dx dy за время dτ войдет cJ иz dx dy dt моль вещества J.
Аналогично через грань dy dz войдет cJ ux dy dz dτмоль вещества J, а через грань dx dz cJ uy dx dz dτмоль J.
Суммарно с конвективным потоком в элементарный объем будет внесено
cJ(uzdxdy + uxdydz + uydxdz)dτ. (4.1)
При прохождении элементарного объема произойдет изменение импульса единицы объема (так как в общем случае и сJ и скорость и имеют неравномерное распределение в пространстве). В результате количество вещества J, которое будет вынесено за тот же промежуток времени dτ через противоположные грани параллелепипеда, составит:
Суммарно по всем осям:
cs(uzdxdy + uxdydz + uydxdz)dτ + [u grad cj + cj div u] dxdydzdτ, (4.2)
где .
После вычитания выражения (4.2) из (4.1) получим (с учетом того, что для несжимаемой жидкости div u = 0) изменение количества вещества в элементарном объеме в результате конвективного переноса за время dτ:
∆nJ, конв= –и grad сJdVdτ. (4.3)
Аналогично получим член уравнения материального баланса, описывающий изменение количества вещества J в результате диффузионного переноса. Диффузионный поток на входе в параллелепипед через грань dxdy (в направлении оси z)в соответствии с первым законом Фика равен
При прохождении потока через элементарный объем произойдет изменение градиента концентрации dcJ/dz на величину (d2cJ /dz2)dz,следовательно, диффузионный поток на выходе из параллелепипеда через противоположную грань составит
Изменение количества вещества J в результате диффузионного переноса через все грани параллелепипеда за время dτ
(4.4)
Расход вещества на химическую реакцию (или его образование в ходе химической реакции) внутри элементарного объема dV заэлементарный промежуток времени dτ пропорционален скорости реакции wrJ(она определяется концентрацией вещества JcJ,установившейся внутри элементарного объема), объему dV ивремени dτ:
. (4.5)
Следует отметить, что в соответствии с формальным правилом о знаках при составлении кинетических уравнений (см. §3.2) ∆nJ,хрположительно, если вещество J – реагент, и отрицательно, если J – продукт. Поэтому для сохранения физического смысла в уравнение материального баланса член ∆nJ,хрдолжен всегда входить со знаком «минус».
Накопление вещества J за время dτ внутри элементарного объема может произойти в результате приращения концентрации сJ при изменении времени на величину dτ. Это изменение концентрации равно (дcj/дτ)dτ. Соответственно накопление вещества в элементарном объеме dV
(4.6)
Таким образом, уравнение материального баланса по веществу J в соответствии с выражениями (4.3)–(4.6) можно записать как:
или, сократив все его члены на dFdx,
(4.7)
Уравнение (4.7) достаточно полно описывает химический процесс, протекающий в любом химическом реакторе (при его выводе не было принято никаких допущений об его применимости только к какому-то одному определенному типу химических реакторов). В нем отражен перенос импульса (первый член уравнения), диффузионный перенос (второй член) и протекание химической реакции (третий член).
Уравнение (4.7) вместе с уравнением теплового баланса, учитывающим явления теплопереноса в элементарном объеме реактора, составят полную математическую модель реактора. Таким образом, будет решен вопрос и о небольшом числе уравнений, составляющих математическую модель, и об ее полноте.
Однако уравнение (4.7) слишком сложно для решения (дифференциальное уравнение второго порядка в частных производных). Следовательно, реальный путь создания математических моделей, пригодных для решения практических инженерных задач по расчету и проектированию химических реакторов, заключается в упрощении математической модели, которое можно провести для различных частных случаев.
В соответствии с такой концепцией рассмотрим математические модели различных типов реакторов:
реакторов для гомогенных процессов, работающих в изотермическом режиме; в неизотермическом режиме; реакторы для гетерогенных процессов.
|