0% found this document useful (0 votes)
28 views12 pages

Non-Isothermal Reactor Design Basics

This document discusses energy balances for non-isothermal chemical reactors. It presents the general energy balance equation and defines terms like heat of reaction. It then derives specific energy balance equations for closed and open systems, as well as for batch, continuous stirred-tank, and plug flow reactors. The equations relate the heat flow, temperatures, heat capacities, heat of reaction, and other terms. Constant heat capacities are also assumed in some cases to simplify the equations.
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
28 views12 pages

Non-Isothermal Reactor Design Basics

This document discusses energy balances for non-isothermal chemical reactors. It presents the general energy balance equation and defines terms like heat of reaction. It then derives specific energy balance equations for closed and open systems, as well as for batch, continuous stirred-tank, and plug flow reactors. The equations relate the heat flow, temperatures, heat capacities, heat of reaction, and other terms. Constant heat capacities are also assumed in some cases to simplify the equations.
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

IQ2011B 

Chemical reactor design

Module 3: Non-isothermal
reactors
IQ2011B – Chemical reactor design

Mass balance general equation


The molar balance for species j is:

𝑑𝑁
𝐹 𝐹 𝑟 𝑑𝑉
𝑑𝑡

1
IQ2011B Chemical reactor design

General energy balance equation


𝑑𝐸
𝑚 𝐸 𝑚 𝐸 𝑄 𝑊
𝑑𝑡
1
𝐸 𝑈 𝐸 𝐸 𝑚𝑢 𝑈 𝑚𝑔ℎ
2
𝑊 𝑊 𝑊
𝐻 𝑈 𝑊 𝑈 𝑝𝑉
Ignoring potential and kinetic energy, and shaft work:
𝑑𝑈
𝑚 𝐻 𝑚 𝐻 𝑄
𝑑𝑡

Energy Balance (Closed Systems)


The first Law of Thermodynamics for Closed Systems states:
Heat Flow
Work of the
Total to the
System on the
Energy System
Surroundings

dE = dQ  dW
gc is the factor to convert
from mass to force
1𝑣 𝑔𝑧 kg m lb ft
𝑑 𝑈 𝑔 1 32.2
2𝑔 𝑔 Ns lb s

At constant volume dU = dQ and also


H = U + pV and for liquids and solids pV is negligible so Q ≈ ∆H

2
IQ2011B Chemical reactor design

Energy Balance (Open Systems)


For open systems, the first law is expressed:
Specific energy Specific energy
associated with associated with
System of the mE mS
System

dEsist   1 v 2 gz   1 v 2 gz 
QW  E  Ĥ 
S m   m S  Ĥ   
dt Entradas  2 g c g c Salidas  2 g c gc 
Heat Shaft
flow work

Steady state and disregarding kinetic and potential energy:


𝑄 𝑚 ∑ 𝐻 𝑚 ∑ 𝐻

Energy Balance (Open Systems)


Due to steady state mE = mS and using molar flows
𝑄 𝐹𝐻 𝐹𝐻

And according to the stoichiometric table


𝑖
𝐹 𝐹 𝐹 𝑥
𝑎
 i   i 
 
Q   Fi 0  FA 0x A Ĥi    Fi 0  FA 0x A Ĥi
Salidas  a  Entradas  a 

3
IQ2011B Chemical reactor design

Energy Balance
Heat of reaction @ T
Rearranging:

c d b 
 
Q F i0 Ĥi  F i0 Ĥi  FA 0x A  ĤC  ĤD  ĤB  Ĥ A 
Salidas Entradas a a a 

And with:
c d b
ĤRXN ( T)  ĤC  ĤD  ĤB  ĤA
a a a

Energy Balance

Simplifying:
 
Q  F (Ĥ  Ĥ
i
i0 i i0 )  FA 0x A ĤRXN ( T)

If there is no phase change: a + bT + cT2 + dT3

T
Ĥi ( T)  Ĥ0fi ( TR )  
TR
Cpi dT

4
IQ2011B Chemical reactor design

Energy Balance
Replacing:
Ĥi ( T)  Ĥi 0 ( T)  Ĥ0fi ( TR )  CpidT  Ĥ0fi ( TR )  CpidT
T Ti 0

  TR   
TR 

T Ti 0
 
TR
CpidT  
TR
Cpi dT

T
 
Ti 0
CpidT

Energy Balance
Replacing

F 
T
 
Q Cpi dT  FA 0x A ĤRXN ( T)
i0
Ti 0
i

And the heat of reaction is expressed:

   
TR T
ĤRXN  Cpi dT  Ĥ0RXN  CpidT
T TR
Re activos Pr oductos

5
IQ2011B Chemical reactor design

Energy Balance
Taking into account the coefficients:
i
 
T
ĤRXN  Ĥ0RXN  Cpi dT
a TR

Convective Heat of
heat flow reaction

 i 
   
T T
 
Q Fi 0 Cpi dT  FA 0x A Ĥ0RXN  Cpi dT
i
Ti 0  a TR 

Heat flow to the


surrounding area

Energy Balance
Using Constant Cp's:
  F C p ( T  T )  F x Ĥ0   i C p ( T  T )
Q i0 i i0 A0 A  RXN i R 
 a 
Simplifying:

  F C p ( T  T )  F x Ĥ0   C p( T  T )
Q i0 i i0 A0 A RXN R  
𝑖
∆𝐶𝑝 𝐶𝑝
𝑎

6
IQ2011B Chemical reactor design

Energy Balance in BRs


For a batch reactor, instead of F molar flows, we have moles N:

Sensitive and
latent heat of
the mixture

i i0 A0 A RXN 
 t  N C p ( T  T )  N x Ĥ0   C p ( T  T )
Q i i R 
Heat flow with Heat of reaction
the
surroundings

This would be a BdeE accumulated from the beginning to the end


of the reaction, where Q = UA(Tm – T)

Energy Balance in BRs


For an instant of time t:
 dt  N Cp dT  N dx Ĥ
Q i i A0 A RXN

But from the BdeM, NA0dxA = (– rAV)dt:


  N Cp dT  ( r V )Ĥ
Q i i A RXN
dt

dT UA( Tm  T)  ( rA V )ĤRXN



dt N iCpi

7
IQ2011B Chemical reactor design

Energy Balance in CSTRs


Taking the equation:

i i0 A0 A RXN 
  F C p ( T  T )  F x Ĥ0   C p ( T  T )
Q i0 i R 
where Q = UA(Tm – T) and A is the heat transfer area.

Energy Balance in PFRs


Taking the equation:
UDL( Tm  T)  Fi 0 C pi ( T  Ti 0 )  FA 0x A ĤRXN
where Q = UA(Tm – T). Again, this is a cumulative BdeE from the
input to the output of the PFR. On a dV:
UdA ( Tm  T)  FiCpi dT  FA 0dx A ĤRXN

8
IQ2011B Chemical reactor design

Energy Balance in PFRs


But FA0dxA = – rAdV. Rearranging:
UDdz( Tm  T)  FiCpi dT  ( rA A Tdz )ĤRXN

where AT = D2/4. Dividing between dz and simplifying:

dT UD( Tm  T)  ( rA A T )ĤRXN



dz FiCpi

Batch
It can be analyzed with the following simplifications:
a) You have a constant temperature or a temperature program.
b) You have a specified heat flow (zero, constant, or a function of
time.
c) The system operates neither isothermal nor adiabatic.

9
IQ2011B Chemical reactor design

a) Specified temperature
To calculate the heat flow, Qtotal is cleared and derived with
respect to time:

𝑑𝑄 𝑑 ∑ 𝑁 𝐶𝑝 𝑇 𝑇 ∆𝐻 𝑁 𝑋
𝑄 𝑡
𝑑𝑡 𝑑𝑡

 H0RXN  CpT  TR 

This gives us the heat flow program required to maintain the


specified temperature conditions.

b) Specified heat flow


The BdeE is a nonlinear equation relating T and XA that must be solved simultaneously with
the EDO of the BdeM. In the case of an adiabatic batch reactor:
𝛥𝐻rxn 𝑁 𝑋 ∑𝑁 Cp 𝑇 𝑇

Clearing T as a function of XA:


Δ𝐻 𝑁 𝑋
𝑇 𝑇
Σ𝑁 𝐶𝑝

When XA tends to 1, it is known as the adiabatic elevation of temperature.

Δ𝐻 𝑁 Complete
𝑇 𝑇 Conversion
Σ𝑁 𝐶𝑝

10
IQ2011B Chemical reactor design

c) Non-isothermal, non-adiabatic
In this case the BdeM equation cannot be solved independently
of the BdeE and only applies a simultaneous solution of the
EDOs.

𝑑𝑋 𝑉
𝑟
𝑑𝑡 𝑁

𝑑𝑇 ∆𝐻 𝑟 𝑉 𝑈𝐴 𝑇 𝑇
𝑑𝑡 ∑ 𝑁 𝐶𝑝

CSTR

i i0 A0 A 
  F C p ( T  T )  F x Ĥ0   C p ( T  T )
Q i0 RXN i R 
where T is the operating temperature of the reactor, Cp
represents the effect of the T on the HRXN, Q = UA(Tm – T)
and A is the heat transfer area.

11
IQ2011B Chemical reactor design

PFR
Mass Balance:
xA
dx A dx A D2
V  FA 0
0  rA

dz 4FA 0
( rA )

Energy Balance: Momentum Balance:

dT UD( Tm  T)  ( rA A T )ĤRXN 𝑑𝑝 𝜌𝑣


 2𝑓
dz FiCpi 𝑑𝑧 𝐷𝑔

Transfer of Energy to the Surroundings


Additionally, if there is heat transfer with the surroundings
(cooling/heating fluid), the BdeE can be considered in the heat
transfer medium:

dTm UD( T  Tm )

dz m
 ECpE

Where mE is the mass flow of the medium and CpE its specific
heat.

12

Common questions

Powered by AI

The heat transfer area (A) plays a crucial role in non-isothermal reactor designs as it determines the rate at which heat can be transferred between the reactor contents and the surroundings. This heat transfer affects the reactor’s energy balance by influencing the temperature gradient and the overall heat flow into or out of the system. In steady-state conditions for both CSTR and PFR, the energy balance is impacted by the equation Q = UA(Tm - T), where U is the overall heat transfer coefficient, A is the heat transfer area, Tm is the temperature of the surroundings, and T is the reactor temperature. A larger area enhances heat transfer, potentially stabilizing temperature fluctuations .

In open reactor systems, energy balances account for material flow across system boundaries, which involves considering the enthalpy of incoming and outgoing streams. The energy balance formula incorporates terms for flow enthalpies, heat transfer, and work done on/by the system as 𝑄̇ = ∑(𝐹𝐻̂)𝑖𝑛 − ∑(𝐹𝐻̂)𝑜𝑢𝑡, indicating exchanges with surroundings. Conversely, in closed systems, energy balances focus on changes in internal energy (dU) since there is no mass transfer across the system boundary. The first law of thermodynamics is expressed as dU = dQ - dW, highlighting internal changes due to heat transfer and work done, with negligible impact from material flow enthalpy .

In non-isothermal reactors, the heat of reaction at temperature T is integrated into energy balance calculations by considering the change in enthalpy due to the reaction progress at a given temperature. This is expressed as Δ𝐻̂𝑅𝑋𝑁(T) in energy equations, which accounts for heat released or absorbed during the reaction. It is included with other terms like the heat capacity integrated over the temperature range to evaluate the total heat effects. These calculations use the formulation 𝐻̂(T) - ∫𝐶𝑝dT to estimate how changes in temperature affect the overall energy structure of the reaction, influencing both the product and the reaction’s thermal behavior .

The non-isothermal, non-adiabatic system requires a more complex solution approach because it involves simultaneous changes in both temperature and reaction kinetics, which influence each other dynamically. Unlike isothermal or adiabatic systems where either temperature is constant or heat exchange is considered negligible, non-isothermal, non-adiabatic systems must solve coupled differential equations for both energy and material balances (BdeM and BdeE). This integration requires considering heat losses/gains and their interaction with the reaction kinetics, making it necessary to perform simultaneous solutions of ordinary differential equations reflecting both concentration and temperature profiles .

The adiabatic temperature rise for complete conversion in a batch reactor is determined by assuming no heat loss to the surroundings and using the relation (T - T₀) = (-ΔĤrxn)N₀ΣCp⁻¹. This represents the maximum temperature increase achievable due to the entire reaction heat being used to raise the mixture's temperatu... as it transforms from reactants to products. It is significant as it defines the thermal limit of the system under adiabatic conditions, crucial for designing safety measures and understanding potential thermal effects that could impact reaction rates, stability, and mechanisms in industrial applications .

The assumption of negligible pV work (pressure-volume work) is valid in energy balance calculations for liquids and solids because these phases typically exhibit minimal volume changes under pressure variations at constant temperature conditions. In these states, the pV product is small and doesn’t significantly impact the overall energy calculations as it does in gaseous phases where volume changes are substantial under pressure differences. Therefore, for liquids and solids, the internal energy change can be approximated primarily by changes in enthalpy, simplifying the energy balance equation to Q ≈ ΔH .

The stoichiometric table approach aids in simplifying energy balance equations by directly linking reactant conversion to heat flow calculations via stoichiometric relationships. It allows for the systematic representation of moles and the related enthalpy terms, thereby transforming complex reaction networks into straightforward algebraic expressions. By defining changes in molar flow rates through stoichiometric coefficients (e.g., F_i = F_i₀ - νA F_A₀ xʁA) and integrating these into enthalpy balances, it enables quick adjustment of balance equations to reflect changes in conversion and reaction conditions, facilitating easier energy accounting .

Performing energy balance calculations for batch reactors requires using time-dependent integration over the process duration, as there is no continuous flow. Unlike continuous flow reactors (CSTRs and PFRs), where steady-state assumptions often simplify balance equations leading to distinct input and output terms, batch reactors involve temporally cumulative energy exchanges. Equations must reflect changes over time, such as integrating specific heat over temperature changes and considering reaction progression (e.g., Q = ∫CpNdT + ∫r֦ dV) to capture temporal variations in both composition and thermal conditions .

Assuming constant heat capacities (Cp) simplifies energy balance calculations by treating the heat absorption or release per unit temperature change as invariant across the temperature range. This assumption allows for straightforward integration of temperature-dependent terms in energy equations, expressed compactly as ∆Cp(∆T), and is justified in small ranges or cases with minimal temperature dependence. However, it can lead to inaccuracies in calculating precise thermal states in systems with significant temperature variation or where Cp varies substantially with temperature, necessitating considering empirical temperature-dependent Cp expressions for a more accurate model .

The factor gc enters the energy balance equations in reactor design to convert units of mass flow into force, particularly when dealing with potential and kinetic energy changes. It ensures consistency between mass and force unit systems, especially critical in calculations involving British units where gc equals 32.2 lb·ft/lbf·s². This conversion factor prevents dimensional inconsistencies in equations that involve gravitational forces and energy expressions, ensuring that kinetic and potential energies are correctly accounted for within the broader thermodynamic framework .

You might also like