MaximaLabs
All unit operations

Polymerization reactor simulation

monomer -> polymer (Mn/Mw/PDI)

Governing equations

The exact equations the solver works for a polymerization reactor — the same math shown in the app's "Theory" panel, not a black box.

[R]=2fkd[I]/kt,Rp=kp[M][R][R^\bullet]=\sqrt{2f k_d [I]/k_t},\qquad R_p=k_p[M][R^\bullet]
ν=kp[M]2fkdkt[I],DPn=ν,    PDI2\nu=\dfrac{k_p[M]}{2\sqrt{f k_d k_t [I]}},\quad \overline{DP}_n=\nu,\;\; \text{PDI}\approx 2
Mn=DPnM0,Mw=PDIMnM_n=\overline{DP}_n M_0,\qquad M_w=\text{PDI}\cdot M_n
[R][R^\bullet]
radical concentration [mol/m³]
ff
initiator efficiency [0..1]
kdk_d
initiator decomposition rate constant
[I][I]
initiator concentration [mol/m³]
[M][M]
monomer concentration [mol/m³]
kpk_p
propagation rate constant
ktk_t
termination rate constant
RpR_p
polymerization rate [mol/m³/s]
ν\nu
kinetic chain length — how many monomers one radical adds before terminating
DPn\overline{DP}_n
number-average degree of polymerization
M0M_0
monomer molar mass [kg/mol]
Mn,MwM_n, M_w
number- and weight-average molecular weight [kg/mol]
PDI\text{PDI}
polydispersity M_w/M_n — ≈2 for this termination mechanism

Parameters

monomer [feed component]; free-radical CSTR: volume [m^3], initiator_conc [mol/m^3], temperature [K], monomer_mw [kg/mol], optional kp0/kp_ea, kt0/kt_ea, kd0/kd_ea Arrhenius constants, initiator_efficiency f — reports conversion, Mn, Mw, PDI

Example flowsheets that use it

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