Green ammonia synthesis — a PENG-ROBINSON process flowsheet
An equilibrium reactor runs N₂ + 3H₂ ⇌ 2NH₃ on green-hydrogen syngas (Keq illustrative) — the new-energy leapfrog workflow.
167 views 0 forks
- Rigorous PENG-ROBINSON thermodynamics, solved by the same engine every simulation runs on.
- 1 unit operations modeled: RX.
- Focus areas: Green ammonia, Haber–Bosch, Equilibrium reactor.
- Thermodynamics
- PENG-ROBINSON
- Components
- n2, h2, ammonia
- Unit operations
- RX
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("green-ammonia")
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
Ammonia synthesis loop (Haber-Bosch, ChemSep casebook)
The classic industrial ammonia loop: fresh syngas joins recycled unconverted gas, reacts over an equilibrium reactor (N₂ + 3H₂ ⇌ 2NH₃) at synthesis conditions, is chilled to condense high-purity liquid ammonia, and the remaining vapor splits into a recycle (back to the loop) and a purge — the purge exists specifically to bleed off the inert argon a single-pass reactor could never consume, which would otherwise concentrate in the recycle forever.
Ethane steam cracker
A simplified thermal-cracking train: an equilibrium reactor converts ethane to ethylene + H₂ at furnace conditions (1100 K), then a quench, compression, and a cold flash strip the H₂ tail gas before a C₂ splitter separates ethylene product from unconverted ethane. An equilibrium reactor is used here rather than the kinetic reactor because the latter's concentration term always uses a liquid-phase density (a real gap for a vapor-phase reaction at 1100 K) — noted, not fixed, out of scope for this example. The C₂ splitter needs many stages and high reflux since ethylene/ethane relative volatility is modest. From the ChemSep casebook (ethane cracker).
Toluene hydrodealkylation (HDA)
Toluene + H₂ react to benzene + methane over an equilibrium reactor (900 K), then a quench, a high-pressure separator, and a letdown-valve-plus-flash stabilizer strip dissolved H₂/methane before the aromatics split — the same SEP-LETDOWN-SEP2 degassing pattern used for methanol synthesis and CO₂-from-natural-gas elsewhere in this gallery. The final benzene/toluene split uses a shortcut separator rather than a rigorous distillation column: even the small H₂/methane trace left after stabilization has a pure-component bubble point far below any real column temperature (H₂ at 1.5 bar: ~22 K), which corrupts the column's initial T-profile guess and the inside-out/Newton MESH solvers land on a physically nonsensical low-temperature root for this feed. A fix narrowed to just that seed was tried and reverted: it also changes the wide-boiling classification and the seed fed into the (otherwise-fine) inside-out solver for other columns, and broke the already-shipped methanol-synthesis example — the same cross-example fragility already hit twice this session trying to patch the shared thermo core, so it's noted honestly here rather than chased further. From the ChemSep casebook (toluene HDA).
Quench-converter ammonia synthesis loop
The multi-bed quench converter BASF first installed in 1942, which displaced the tube-cooled converters before it, running inside a complete 200 bar synthesis loop. Ammonia synthesis is equilibrium-limited and exothermic, so the catalyst wants to run hot for rate and cool for equilibrium, and no single adiabatic bed can do both. The quench converter's answer is mechanical: split the loop gas, send 55% through the interchanger into bed 1 and inject the other 45% as cold shots between the beds, so each bed heats up adiabatically and each quench pulls it back down. What it computes: the sawtooth, 673 → 766 K in bed 1, quenched to 689 K and out at 760 K, quenched to 712 K and out at 766 K; 12.5 mol% ammonia at the converter exit, in the 12-18% band real quench converters deliver, at 20.9% nitrogen conversion per pass; 44.6 mol/s of 97.8% liquid ammonia off the 250 K separator; and argon and methane held at 3.3% and 4.2% in the loop by a 3% purge, with the recycle closed by a Wegstein tear. The acid test is the counterfactual: one adiabatic bed producing that same 12.51 mol% ammonia exit ends at 834.7 K — 562 C, past the sintering limit of promoted-magnetite catalyst — while the three-bed quench converter delivers the identical duty with a 766 K (493 C) peak.
High-recycle ammonia loop (equation-oriented)
A tight, high-recycle ammonia synthesis loop built to show why a modern simulator solves recycles the way legacy sequential-modular tools can't. Fresh syngas (N₂ + 3H2, with argon inert) mixes with a large recycle, reacts to only ~10% per pass, chills so ammonia condenses out as product, and the unreacted gas recycles — a recycle-to-fresh ratio of ~5:1, with argon building up until a small purge balances it. Run this in Equation-Oriented mode (Solver menu > Mode > Equation-oriented). In the default sequential-modular mode the solver tears the recycle and iterates Wegstein ~78 times to close the loop; the equation-oriented solver instead makes every inter-unit stream a global unknown and closes all ~35 of them in one simultaneous Newton solve — the same simultaneous approach AVEVA SimCentral / gPROMS / IDAES use, and the reason tightly coupled recycles that crawl (or stall) in sequential-modular converge cleanly here.
CO₂ methanation (e-fuels / power-to-gas)
The Sabatier reaction (CO₂ + 4H₂ ⇌ CH₄ + 2H₂O) converts captured CO₂ and green H₂ into synthetic natural gas at a 97% single-pass-equivalent conversion (real plants stage several adiabatic beds with intercooling to reach this; lumped into one reactor here), then a cooler and knockout drum condense the reaction water from the SNG — the e-fuels loop that closes the loop on captured carbon instead of just storing it.