MaximaLabs
All equations

Adiabatic mixer energy balance

Material & energy balance

Two streams mix adiabatically; the outlet temperature is the enthalpy-weighted average of the inlets — the linear building block of every energy balance.

Governing equation

n˙1cp(T1Tout)+n˙2cp(T2Tout)=0\dot n_1 c_p (T_1 - T_{\mathrm{out}}) + \dot n_2 c_p (T_2 - T_{\mathrm{out}}) = 0
ToutT_{\mathrm{out}}
Mixed-outlet temperature (K) — the unknown.
n˙1,n˙2\dot n_1,\,\dot n_2
Inlet molar flows (mol s⁻¹).
T1,T2T_1,\,T_2
Inlet temperatures (K).
cpc_p
Molar heat capacity (J mol⁻¹ K⁻¹), taken constant.

Derivation

  1. 1
    inn˙jhj(Tj)=n˙outh(Tout)\sum_{\mathrm{in}} \dot n_j h_j(T_j) = \dot n_{\mathrm{out}} h(T_{\mathrm{out}})

    Steady-state energy balance on an adiabatic mixer: enthalpy in equals enthalpy out.

  2. 2
    h(T)h(Tref)=cp(TTref)h(T) - h(T_{\mathrm{ref}}) = c_p (T - T_{\mathrm{ref}})

    With a constant heat capacity the enthalpy is linear in temperature about a reference.

  3. 3
    n˙1cp(T1Tout)+n˙2cp(T2Tout)=0\dot n_1 c_p (T_1 - T_{\mathrm{out}}) + \dot n_2 c_p (T_2 - T_{\mathrm{out}}) = 0

    Substituting and cancelling the reference gives the balance; T_out is the flow-weighted mean inlet temperature.

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.

∂/∂Tout
Adiabatic energy balance

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

Symbolic derivative
F1Tout=\dfrac{\partial F_{1}}{\partial\, \mathrm{Tout}} =
(n˙1+n˙2)cp-(\dot n_1 + \dot n_2)\,c_p
Evaluated at the operating point: -300

For the residual Adiabatic energy balance with respect to Tout. 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.

Tout = 337.5n1 = 3n2 = 1cp = 75T1 = 350T2 = 300

References

  • Felder & Rousseau, Elementary Principles of Chemical Processes, 3rd ed., Ch. 8.

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