MaximaLabs
All equations

CSTR mole balance (2nd-order)

Reaction

The steady-state mole balance on a continuous stirred-tank reactor with a second-order rate law — outlet concentration from a single nonlinear balance.

Governing equation

V˙(C0C)kC2V=0\dot V (C_0 - C) - k\,C^2\,V = 0
CC
Outlet (= tank) concentration (mol m⁻³) — the unknown.
V˙\dot V
Volumetric flow rate (m³ s⁻¹).
C0C_0
Inlet concentration (mol m⁻³).
kk
Second-order rate constant (m³ mol⁻¹ s⁻¹).
VV
Reactor volume (m³).

Derivation

  1. 1
    dNdt=V˙C0V˙C+rAV=0\frac{dN}{dt} = \dot V C_0 - \dot V C + r_A V = 0

    Mole balance on the well-mixed tank at steady state: in − out + generation = 0.

  2. 2
    rA=kC2r_A = -k\,C^2

    Second-order rate law; the tank is at the outlet composition because it is perfectly mixed.

  3. 3
    F(C)=V˙(C0C)kC2V=0F(C) = \dot V (C_0 - C) - k\,C^2 V = 0

    A single nonlinear (quadratic) equation for the outlet concentration.

Jacobian

The matrix Jij=Fi/xjJ_{ij} = \partial F_i / \partial x_j of the residual vector with respect to the unknowns, at the representative operating point. Every value below is computed by the exact forward-mode autodiff engine the solver uses — click a non-zero cell to see the symbolic derivative, each one verified against that engine in our test suite. This is the block a Newton step factorizes.

∂/∂C
Steady-state mole balance

Shaded cells are non-zero (the sparsity pattern); blanks are structural zeros.

Symbolic derivative
F1C=\dfrac{\partial F_{1}}{\partial\, \mathrm{C}} =
V˙2kCV-\dot V - 2 k C V
Evaluated at the operating point: -4

For the residual Steady-state mole balance with respect to C. This exact expression is asserted equal to the autodiff value in the test suite.

Operating point

Representative operating point — a real solution of this instance.

C = 0.66667Fv = 2C0 = 1k = 0.5V = 3

References

  • Fogler, Elements of Chemical Reaction Engineering, 5th ed., Ch. 2 & 5.

Stop fighting legacy software. Build your first flowsheet in 60 seconds.