Pressure-swing azeotropic distillation (methanol/acetone) — a NRTL process flowsheet
Methanol and acetone form a minimum-boiling azeotrope whose composition shifts noticeably with pressure — enough that two ordinary columns run at different pressures can cross the azeotrope without an entrainer. Column 1 (low pressure) takes pure acetone bottoms and a near-azeotropic overhead; that overhead is compressed into Column 2 (high pressure), whose new, shifted azeotrope composition lets pure methanol come off the bottoms (the recycle of Column 2's now off-azeotrope overhead back to Column 1 is reported as a product here rather than closed, the same open-loop simplification the acetone-water example already uses). Luyben, Ind. Eng. Chem. Res. 2008, 47, 2696.
346 views 0 forks
- Rigorous NRTL thermodynamics, solved by the same engine every simulation runs on.
- 2 unit operations modeled: 2× Hpcol, COMP.
- Focus areas: Pressure-swing distillation, Azeotrope, Two-column train.
- Thermodynamics
- NRTL
- Components
- methanol, acetone
- Unit operations
- 2× HpcolCOMP
Opens in a new tab, loaded straight into the app — no setup.
Read the step-by-step guideReproduce this exact result from Python — the real client.get_example() → run_and_wait() path, not a mockup.
from flowsim.sdk import FlowSimClient
client = FlowSimClient()
example = client.get_example("pressure-swing-methanol-acetone")
sim = client.create_simulation(example["title"], example["flowsheet"])
result = client.run_and_wait(sim["id"])
print(result["status"]) # "converged"
streams = client.streams(sim["id"])Related models
Ethanol–water distillation
An 8-stage column concentrating ethanol overhead toward the azeotrope (the headline demo).
Pressure-swing ethanol dehydration (Gᴱ mixing rule)
Ethanol–water is the classic azeotrope, and pressure-swing distillation breaks it without an entrainer: the azeotrope moves with pressure, so a low-pressure column and a high-pressure column pass each other's azeotropic distillate and each recovers a pure product. The whole process only works if the property package tracks that shift — which is exactly where a conventional package choice falls between two chairs. This flowsheet runs the high-pressure column at 15 bar on pr-mhv1: Peng-Robinson with an MHV1 excess-Gibbs mixing rule, so the cubic equation of state gets its attraction parameter from NRTL's excess Gibbs energy instead of from a single binary interaction constant. Switch the thermo package (Solver menu) and compare the predicted azeotrope: | package | 1 atm | 15 bar | valid at 15 bar? | |---|---|---|---| | NRTL | 0.891 | 0.802 | no — γ-φ is a low-pressure formulation (~10 bar) | | Peng-Robinson (kij) | 0.586 | 0.613 | yes, but a kij cannot represent this azeotrope | | pr-mhv1 | 0.949 | 0.798 | yes | (mole fraction ethanol; the repo's DECHEMA-validated 1 atm anchor is 0.894.) At 15 bar pr-mhv1 lands within 0.005 of NRTL while remaining a genuine equation of state, whereas plain Peng-Robinson is off by ~0.19 and puts the azeotrope in the wrong place entirely. Selecting nrtl here also trips the applicability guard, which warns that the activity model is past its pressure ceiling and names the fix. The flowsheet demonstrates the mechanism on itself. Drop the column pressure to 1 atm and re-run, changing nothing else: the solve fails with SPEC_THERMODYNAMICALLY_IMPOSSIBLE, because at atmospheric pressure the requested bottoms purity sits beyond the azeotrope and no column can reach it. At 15 bar the same specification converges and the bottoms leaves at x_EtOH ≈ 0.924 — past the atmospheric azeotrope of 0.894, which is precisely the composition an atmospheric column cannot cross.
Ethanol dehydration (pervaporation)
A hydrophilic pervaporation membrane pulls water out of near-azeotropic ethanol — crossing the 89 mol% ethanol-water azeotrope that ordinary distillation cannot. Permeate is ~95% water; the retentate is dried past the azeotrope.
Extractive distillation (ethanol/water with ethylene glycol)
Near-azeotropic ethanol/water (85 mol% ethanol) can't be pushed past the 89 mol% azeotrope by ordinary distillation. A heavy, high-boiling entrainer (ethylene glycol) fed near the top of the column raises water's relative volatility enough to pull overhead ethanol past the azeotrope in one pass; a second column then strips water overhead from the glycol, regenerating it as the bottoms product. The classic extractive-distillation textbook case (Luyben, Ind. Eng. Chem. Res. 2006, 45, 4625).
Ethanol-water column at 6 atm (SRK + MHV1)
The ethanol-water column everyone knows, run at 6 atm on the Gᴱ-mixing-rule package: Soave-Redlich-Kwong whose mixture energy parameter comes from NRTL's excess Gibbs energy through the MHV1 rule (Michelsen 1990) instead of a kij — the package built for polar, hydrogen-bonding mixtures at pressure, where a γ-φ activity model is past its ~10 bar ceiling and a plain cubic with van der Waals mixing has no hydrogen bonding at all. A 10 mol% ethanol feed through 16 stages at a reflux ratio of 2.5 gives a 66.7 mol% distillate and an ethanol-free bottoms, 2.6 MW on the reboiler. Read it against NRTL, which at 6 atm is still inside its own range: the two packages put the pressure-shifted azeotrope in different places — 0.88 mol fraction ethanol here, 0.84 on NRTL at 6 atm (0.97 against 0.89 at 1 atm) — so this card shows where the MHV1 rule sits, not a number to prefer over the fitted binary at low pressure. The cut is kept well short of the azeotrope on purpose: at a distillate ratio of 0.12 the profile pinched against it and the column reported exactly that. Each solve is about a minute, because every stage K-value is a Gᴱ mixing-rule evaluation.
n-Hexane / 1-butanol column with no fitted binary (UNIFAC)
A binary the databank has no regressed NRTL parameters for, so the column runs on predictive original UNIFAC from the two molecules' groups alone: 50/50 n-hexane and 1-butanol, 20 stages, reflux 2, 45% distillate. UNIFAC predicts the minimum-boiling azeotrope at 97 mol% hexane and 341.7 K, and the column lands its distillate exactly there (96.7%) with an 88% butanol bottoms — the azeotrope, not the stage count, is what caps the overhead purity, and that is a prediction from group contributions with no data on this pair behind it. Read it as a screening result: original UNIFAC's mean error against fitted binaries in this tree is 0.08 in ln gamma, and an alkane/alcohol pair is in its well-behaved range.