NGL debutanizer with a bottoms recycle — a PENG-ROBINSON process flowsheet
The flowsheet a real CAPE-OPEN interchange deck describes, built on the three specifications that deck states and this product can now solve. An NGL cut (16% propane / 13% isobutane / 26% n-butane / 45% n-pentane, 40 mol/s at 15.5 bar) mixes with a bottoms recycle, is preheated to exactly 10 mol% vapour, and is fractionated in a 22-stage column with a 1.40 → 1.45 MPa pressure profile; 15% of the bottoms returns to the feed and the loop is torn and converged.
The preheater is specified by the VAPOUR FRACTION it must leave behind (outlet_vapor_fraction = 0.10), not by a duty or an outlet temperature — which is how a column preheater is always specified, because the feed-stage flash is what the number is for. The solver roots the duty for it: 0.404 MW, landing the feed at 373.6 K. The column carries a real PRESSURE PROFILE rather than one pressure: trays resist vapour flow, so the condenser sits at 14.00 bar and the reboiler at 14.50 bar and each stage boils at its own value. The cut is set the way the deck sets it, by the absolute distillate rate: 20.833 mol/s (75 kmol/h).
Result: an LPG distillate at 345.3 K carrying 30.7% propane, 23.7% isobutane, 43.6% n-butane and 2.0% n-pentane, against a 411.9 K bottoms of 91.8% n-pentane with 6.9% n-butane. Reboiler 1.021 MW, condenser 1.003 MW. The recycle settles at 3.382 mol/s — a fixed point of the loop's own balance, specified nowhere.
The finding worth carrying, because it is a trap this deck sets. The only composition the deck states anywhere sits on its recycle stream, and that stream is 96% n-pentane — only 11.6 mol% propane + butanes. Read it as a guide to the FEED and the deck's own 75 kmol/h distillate rate stops being achievable: a clean overhead would need a feed of at least 179 mol/s, so at the 43.4 mol/s this flowsheet runs, 76% of the distillate would have to be n-pentane by MATERIAL BALANCE alone — whatever the reflux, with no solver involved. Underwood agrees from the other side, putting minimum reflux at 8.8 for a 95%/98% key recovery on that composition against 1.465 on the NGL feed shipped here. That stream is the BOTTOMS, not the feed, and on the real feed the deck's reflux of 2.0 is 1.37x minimum: a sensible design point, so the deck is self-consistent after all. The same lesson as a distillate ratio set equal to a feed's light fraction — the cut specification has to fit inside what the feed carries.
Bounds, stated rather than implied. The fresh feed's rate and composition are ours — the deck states neither, and the recycle estimate it does state is for a scale it never gives. The feed is mixed ON-STAGE (the whole feed joins its stage's balance); the deck declares Above-Stage, RadFrac's default, which differs by about 0.02 K on a feed this lightly vaporised. The reboiler is a kettle. The recycle returns at the bottoms pressure, so the mixed feed sits at 14.50 bar and the preheater's 0.5 bar drop lands it at the column's top pressure, 0.24 bar below its own feed stage — a real plant closes that with a feed pump the deck does not show. The tear is seeded with the converged circulation, which cuts the solve from 109 s to 32 s; a rough seed converges to the same fixed point, slower.
45 views 0 forks
- Rigorous PENG-ROBINSON thermodynamics, solved by the same engine every simulation runs on.
- 4 unit operations modeled: M 101, E 101, C 101, S 101.
- Focus areas: Debutanizer, Pressure profile, Vapour-fraction preheater, Recycle, CAPE-OPEN import.
- Thermodynamics
- PENG-ROBINSON
- Components
- propane, isobutane, n_butane, n_pentane
- Unit operations
- M 101E 101C 101S 101
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("ngl-debutanizer-bottoms-recycle")
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
Helium liquefaction (Linde-Hampson cycle)
A Joule-Thomson liquefaction loop: helium is compressed, aftercooled, precooled well below its ~40-51 K JT-inversion temperature, then throttled through a valve into a flash drum that draws off liquid helium while the cold vapor recycles back to the compressor suction — the same JT-cycle physics as the platform's LNG/refrigeration loops, run on a fluid whose inversion curve makes precooling mandatory before the throttle can cool it at all.
LOHC hydrogen release (methylcyclohexane dehydrogenation)
A liquid organic hydrogen carrier (LOHC) power train: methylcyclohexane (MCH) is preheated and catalytically dehydrogenated (MCH → toluene + 3H2, endothermic, equilibrium-limited to ~90% single-pass conversion) to release hydrogen for downstream power generation. A flash drum separates the H₂-rich gas from the liquid; a distillation column then recovers the unreacted MCH overhead and recycles it to the reactor feed, closing a real tight liquid recycle loop, while spent toluene leaves the bottoms (to be re-hydrogenated back to MCH off-site, closing the supply-chain loop — out of scope for this dehydrogenation-side flowsheet). Screening-fidelity note: MCH and toluene are genuinely close-boiling (relative volatility ~1.4 under this thermo package at column conditions), so the recycle carries real toluene along with the recovered MCH rather than a sharp cut — a production column would use more stages to purify it further; this shows the recycle topology and its convergence, not an optimized column design.
Recycle loop
A mixer/heater/splitter loop with a tear stream — exercises Wegstein recycle convergence.
Propane refrigeration cycle
A single-stage vapor-compression refrigeration loop: propane vapor is compressed, condensed to subcooled liquid, throttled through a valve, and superheated back to vapor in the evaporator before recycling — a small makeup feed and purge close the mass balance (inspired by ChemSep's Refrigeration_* reference cases).
Propane (R-290) heat pump heating a hydronic loop
A single-stage vapor-compression heat pump delivering heat to a space-heating water loop — the same cycle as a refrigerator, read from the hot end. Propane evaporates at 4.7 atm (273 K, an air- or ground-source evaporator at 0 °C), is compressed to 21 atm (condensing at 333 K, 60 °C) and condenses against the heating water in a real two-sided exchanger, so the delivered heat is a water stream you can read: 40 °C return warmed to 58 °C. A small makeup feed and purge close the refrigerant mass balance around the recycle, as in the refrigeration examples. Two numbers define the machine: 74 kW delivered for 21 kW of compressor work, a heating COP of 3.5 — 63% of the Carnot limit for this 60 K lift at the compressor's 0.75 isentropic efficiency. Read COP on the flowsheet as the user variable cop_heating, or sweep the condensing pressure: from 18 to 27 atm the COP falls 3.8 → 3.1 as the lift grows, the trade every heat pump makes between water temperature and electricity. Honest scope: the compressor is a fixed isentropic efficiency, the condenser outlet is a specified subcooling, and there is no defrost, part-load or refrigerant-charge model. Below ~18 atm the specified 326 K condenser outlet sits above saturation and the cycle stops condensing — the sweep starts where the physics does.
Ammonia refrigeration cycle
A single-stage vapor-compression refrigeration loop using ammonia — the classic industrial (rather than HFC) refrigerant — condensed to subcooled liquid, throttled through a valve, and superheated back to vapor in the evaporator before recycling (inspired by ChemSep's Refrigeration_Ammonia-30C reference case).