Sulfolane extractive distillation of aromatics — a PSRK process flowsheet
Extractive distillation recovering benzene and toluene from a stabilized reformate's C₆-C₇ non-aromatics (represented by methylcyclohexane) using sulfolane, after Figure 10.2/10.6 of T. Brouwer (PhD thesis, TU Twente, 2021). A vacuum solvent-recovery column then splits the aromatics from the sulfolane (reported as its own product rather than recycled). Both columns converge (ED 44 iterations, SR 98, residuals ~1e-5): methylcyclohexane leaves overhead at 91.3 mol%, the aromatics product carries no sulfolane at all (0.00000 mole fraction) at 95.2% recovery of the feed's benzene and toluene, and the recovery column returns sulfolane 99.96% pure at 490 K.
Both columns were originally mis-specified in the same way, and it is worth seeing once because the symptom looks like a thermodynamic failure rather than an arithmetic one. distillate_to_feed is a fraction of a column's TOTAL feed, and both were set above the amount of light key actually present, so each column was forced to drag its heavy key overhead to make up the flow. The ED asked for 0.15 of 140 mol/s = 21 mol/s of overhead when only 17.5 mol/s of non-aromatics exists, so 3.5 mol/s of aromatics had to come over — that was the 19.4% benzene lost to the raffinate. The recovery column asked for 0.42 of its bottoms against roughly a quarter of it being aromatics, so ~20 mol/s of solvent had to come over — that was the 42 mol% sulfolane in the product. Setting each to the light key's real share (0.125 and 0.265) removes both.
The property method is the whole story here, and it is why this example ran unconverged for a long time. Extractive distillation exists because of solvent selectivity, so a cubic with van der Waals mixing rules — Peng-Robinson as this was originally written — cannot represent the one effect the column depends on; it was being asked to converge on a model that did not contain the physics. A plain activity model cannot be used either: sulfolane boils at 558 K, so a reboiler hot enough to strip it drives benzene past its 562 K critical point, where a gamma-phi formulation has no answer at all. What the system needs is both at once, which is exactly what a Ge mixing rule provides — PSRK is a cubic (no supercritical ceiling) whose mixing is driven by UNIFAC (real selectivity). pr-mhv1 converges here too, in 47 iterations.
Modeling assumptions & limitations
- 1And these are trades rather than defects. Sending exactly the non-aromatic flow overhead means the ED's imperfect split leaves ~9% of the methylcyclohexane in the bottoms, where it reports to the aromatics product at 4.7 mol% — off-spec for a benzene product, and the lever that would clear it is a sharper ED (more reflux, or the higher solvent-to-feed ratio a real sulfolane unit runs). The solvent is still reported as its own product rather than recycled, so the loop is open and the makeup/purge balance a real unit turns on is not modelled. Reflux, not stages, is what the recovery column responds to here: 24 stages at reflux 2.0 still hits the iteration cap where 16 stages at 2.5 converges.
166 views 0 forks
- Rigorous PSRK thermodynamics, solved by the same engine every simulation runs on.
- 1 unit operations modeled: 2× SR.
- Focus areas: Extractive distillation, Sulfolane, Aromatics, Reformate.
- Thermodynamics
- PSRK
- Components
- benzene, toluene, methylcyclohexane, sulfolane
- Unit operations
- 2× SR
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("sulfolane-extractive-distillation")
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
MCHT extractive distillation with phenol
Extractive distillation of methylcyclohexane (MCH) and toluene using phenol as the selective solvent, adapted from Tiverios and Van Brunt (Ind. Eng. Chem. Res. 2000, 39, 1614). Phenol raises toluene's relative volatility away from MCH enough to split an otherwise close-boiling pair; a solvent-recovery column then splits toluene from the phenol (reported as its own product rather than recycled).
Sulfolane liquid-liquid extraction of aromatics
Liquid-liquid extraction of aromatics from a refinery catalytic-reformer stream with sulfolane, after Figure 10.1 of T. Brouwer (PhD thesis, TU Twente, 2021). A countercurrent extraction cascade pulls benzene and toluene preferentially into the sulfolane-rich extract, leaving a methylcyclohexane-rich raffinate.
Reformate splitter
A single column cutting catalytic reformate into a light benzene/hexane overhead and a toluene/xylene bottoms — the IECR 50, 5680 configuration from the ChemSep casebook.
Reformate splitter via dividing-wall column
A refinery catalytic-reformer stream (benzene/toluene/o-xylene plus a non-aromatic n-heptane proxy) is split into three purified aromatics cuts by a single thermally-coupled dividing-wall column instead of two ordinary columns in series — the same DWC economics as the BTX example, applied to a genuine reformate splitter duty. Dejanovic, Matijasevic, Jansen, Olujic, "Designing a Packed Dividing Wall Column for an Aromatics Processing Plant," Ind. Eng. Chem. Res. 2011, 50, 5680.
Aromatics column (Strigle)
A simple aromatics splitter recovering benzene overhead from a toluene/xylene-heavy feed, as described by R. Strigle (Gulf Publishing, 1987) — a classic packed-column textbook design case.
Dilute acetic acid recovery by extractive distillation
Recovering acetic acid from a dilute aqueous stream (e.g. cellulose acetate or terephthalic acid plant wastewater) using diisobutyl ketone (DIBK) as an extractive entrainer. Because acetic acid and water form no azeotrope but have a relative volatility that flattens toward 1 as the mixture concentrates, straight distillation needs an impractically tall column and heavy reboil duty; the entrainer raises water's relative volatility so the extractive column takes water overhead in far fewer stages, and a second column then strips the entrainer from the acid bottoms for recycle.