How to simulate btx aromatics train
Benzene / toluene / p-xylene split in a two-column train — 99+% purity on all three products. From the ChemSep casebook (CScasebook_BTX).
Also known as: BTX unit, benzene toluene xylene separation, aromatics train.
- 1Open the ready-made model
Open the "BTX aromatics train" model in the MaximaLabs workspace — no install, no license. It loads live on the canvas, ready to edit and run.
- 2Confirm the thermodynamics
This process is modeled with the PENG-ROBINSON property package over benzene, toluene, p_xylene — already selected, so the phase equilibrium and enthalpy are physically consistent from the first run.
- 3Review the flowsheet
The flowsheet chains 2× TOLU. Every block is a real, solvable unit op you can reconfigure on the canvas.
- 4Run the simulation
Click Run. The deterministic solver converges the material and energy balances (recycles included) and fills the live stream table — the AI never invents a number.
- 5Read the results and iterate
Inspect the converged streams, tweak a spec, and re-run — or ask the AI copilot to explain a result or diagnose a failed solve in plain English.
- Thermodynamics
- PENG-ROBINSON
- Components
- benzene, toluene, p_xylene
- Unit operations
- 2× TOLU
Opens live on the canvas — free, no install.
Explore the model & flowsheetFrequently asked questions
- What does the BTX aromatics train model simulate?
- Benzene / toluene / p-xylene split in a two-column train — 99+% purity on all three products. From the ChemSep casebook (CScasebook_BTX).
- Is "BTX aromatics train" the same as a BTX unit?
- Yes — this model covers what is also called BTX unit, benzene toluene xylene separation, aromatics train. It runs the real process on the rigorous solver, so you can size and study it directly.
- Which thermodynamic method does it use?
- The PENG-ROBINSON property package, over benzene, toluene, p_xylene — already selected. You can switch the method on the canvas before running.
- Which unit operations are in the flowsheet?
- It chains 2× TOLU. Every block is a real, solvable unit operation you can reconfigure, add to, or remove.
- Do I need to install software or buy a license?
- No. BTX aromatics train runs entirely in your browser on MaximaLabs — free, no install, no license. Open the model to load it live and run the deterministic solver.
More guides like this
BTX separation via dividing-wall column
Benzene/toluene/ethylbenzene split into three purified products by a single rigorous dividing-wall column (prefractionator + two coupled main columns) instead of two ordinary columns in series — the thermally-coupled route real DWC retrofits use.
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.
Natural gas fractionation train
A four-column NGL fractionation train — demethanizer, deethanizer, depropanizer, debutanizer — recovering pipeline-spec sales gas plus ethane, propane, and butane products from wellhead-pressure raw natural gas, after Luyben (Ind. Eng. Chem. Res. 2013, 52, 10741).
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 extractive distillation of aromatics
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.
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.