MaximaLabs
All unit operations

Membrane distillation simulation

Governing equations

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

J=Cm[psat(Tfm)awpsat(Tpm)]J = C_m\left[p^{sat}(T_{fm})\,a_w - p^{sat}(T_{pm})\right]
q=hf(TfTfm)=JΔHvap+kmδ(TfmTpm)=hp(TpmTp)q = h_f\left(T_f - T_{fm}\right) = J\,\Delta H_{vap} + \frac{k_m}{\delta}\left(T_{fm} - T_{pm}\right) = h_p\left(T_{pm} - T_p\right)
TPC=TfmTpmTfTp\mathrm{TPC} = \frac{T_{fm} - T_{pm}}{T_f - T_p}
n˙dist=JA(solvent only — vapour alone crosses)\dot{n}_{dist} = J\,A \quad\text{(solvent only — vapour alone crosses)}
JJ
vapour flux across the membrane [mol/m^2/s]
CmC_m
membrane distillation coefficient — measured per module [mol/(m^2 s Pa)]
psatp^{sat}
saturation vapour pressure of the solvent [Pa]
awa_w
water activity of the feed — salt lowers it, and with it the driving force [-]
TfT_f
bulk feed temperature [K]
TpT_p
bulk permeate temperature [K]
TfmT_{fm}
temperature at the feed side of the membrane, below the bulk [K]
TpmT_{pm}
temperature at the permeate side of the membrane, above the bulk [K]
hfh_f
feed-side film heat-transfer coefficient [W/m^2/K]
hph_p
permeate-side film heat-transfer coefficient [W/m^2/K]
km/δk_m/\delta
conductive loss through the membrane per unit area [W/m^2/K]
ΔHvap\Delta H_{vap}
latent heat carried by every mole that crosses [J/mol]
TPC\mathrm{TPC}
temperature polarisation coefficient — the fraction of the measured driving force that reaches the membrane; real modules run 0.4-0.7 [-]
AA
membrane area [m^2]
n˙dist\dot{n}_{dist}
distillate rate [mol/s]

Parameters

membrane_area [m^2, REQUIRED — the module size, which cannot be inferred]; permeate_temperature [K, the cold side]; optional membrane_coefficient [mol/(m^2 s Pa), default 2.78e-5 = the mid-range literature 5e-7 kg/(m^2 s Pa)], feed_heat_transfer_coefficient / permeate_heat_transfer_coefficient [W/m^2/K, default 3000 — these set temperature polarisation and therefore the real flux], membrane_conductance [W/m^2/K, default 15, conductive loss through the membrane], solvent [default 'water']. Outlets in order: distillate, concentrate.

Example flowsheets that use it

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