How to simulate thermal (non-catalytic) hda — broad aromatics feedstock

A real thermal-HDA licensor's process (non-catalytic, tolerant of a mixed aromatics feedstock and inert impurities) modeled as a genuine **competing-and-sequential reaction network** in a real plug-flow reactor, not a fixed-Keq equilibrium reactor like the simpler hda-toluene-dealkylation example. The main toluene dealkylation (toluene + H₂ → benzene + methane) uses **real published kinetics**: (-r_toluene) = k0*exp(-Ea/RT)*C_toluene*C_H2^0.5 (order 1 in toluene, order 0.5 in H₂, the classic 1.5-order thermal-HDA rate law), k0 = 5.67e9 (m3/mol)^0.5/s, Ea = 228,200 J/mol, heat of reaction -49,000 J/mol — transcribed from Shull & Hixon (Ind. Eng. Chem. Process Des. Develop. 5, 147, 1966) as tabulated in Rase, H.F., 'Chemical Reactor Design for Process Plants, Vol. 2: Case Studies and Design Data' (Wiley-Interscience, 1977), Case Study 104 'Toluene Dealkylation', pp. 36-42 (unit-converted from the source's imperial (ft, lb-mol, degR) basis to SI). That same source explicitly states 'rate equations for the side reactions have not been reported' — so the coking/heavies side reaction (2 toluene → biphenyl + H₂, the real reaction that caps aromatic yield below 100% and drives decoking intervals) and the xylene demethylation feeding the shared toluene pool (o-xylene + H₂ → toluene + methane — the real single-methyl-loss pathway, not a lumped double-dealkylation straight to benzene, so this is a genuine reactions-in-series-and-parallel network, not three independent reactions) both stay illustrative order-of-magnitude Arrhenius parameters, tuned only to land in the source's cited ~98% main-reaction selectivity at 2:1 H₂:toluene — honestly labeled as such, not fabricated citations. H₂S/CO₂/ammonia ride through the flowsheet as inert impurities — exactly what 'impurity tolerance' means in a mass balance, not a special chemistry path. A real fired heater (rigorous combustion stoichiometry + stack energy balance, the same unit op used in the naphtha-steam-cracker example) brings the mixed fresh-plus-recycle feed up to reaction temperature — not a feed specified as already at 950 K. The H₂-rich high-pressure separator vapor is a real recycle loop: fresh feed joins recycled gas at a mixer ahead of the furnace, and past the separator the vapor splits into a recompressed recycle (85%) and a genuine purge (15%) that bleeds off the methane a single pass could never consume — the same role a purge plays in the ammonia-synthesis-loop example. Quench, the high-pressure separator, and a letdown-valve-plus-flash stabilizer strip dissolved H₂/methane before the aromatics split. The reactor is genuinely **adiabatic**, not isothermal: temperature is a real energy-balance-coupled ODE state (the adiabatic PFR mode) that rises along the reactor length as the exothermic reactions release heat — starting from the furnace's real ~900 K outlet (in the source's cited 894-922 K inlet range) and climbing as the reaction proceeds, self-consistent to within ~0.01% of the reactor's own required external duty (which should be ~0 for a truly adiabatic reactor — verified, not assumed). One honest simplification remains: the final benzene/toluene/xylene-derived-benzene split still uses a shortcut separator rather than a rigorous column, for the same documented reason as hda-toluene-dealkylation: trace H₂/methane corrupts the column's wide-boiling initial-guess classification, and a fix narrowed to just that seed broke the shipped methanol-synthesis example when tried previously.

Frequently asked questions

What does the Thermal (non-catalytic) HDA — broad aromatics feedstock model simulate?

A real thermal-HDA licensor's process (non-catalytic, tolerant of a mixed aromatics feedstock and inert impurities) modeled as a genuine **competing-and-sequential reaction network** in a real plug-flow reactor, not a fixed-Keq equilibrium reactor like the simpler hda-toluene-dealkylation example. The main toluene dealkylation (toluene + H₂ → benzene + methane) uses **real published kinetics**: (-r_toluene) = k0*exp(-Ea/RT)*C_toluene*C_H2^0.5 (order 1 in toluene, order 0.5 in H₂, the classic 1.5-order thermal-HDA rate law), k0 = 5.67e9 (m3/mol)^0.5/s, Ea = 228,200 J/mol, heat of reaction -49,000 J/mol — transcribed from Shull & Hixon (Ind. Eng. Chem. Process Des. Develop. 5, 147, 1966) as tabulated in Rase, H.F., 'Chemical Reactor Design for Process Plants, Vol. 2: Case Studies and Design Data' (Wiley-Interscience, 1977), Case Study 104 'Toluene Dealkylation', pp. 36-42 (unit-converted from the source's imperial (ft, lb-mol, degR) basis to SI). That same source explicitly states 'rate equations for the side reactions have not been reported' — so the coking/heavies side reaction (2 toluene → biphenyl + H₂, the real reaction that caps aromatic yield below 100% and drives decoking intervals) and the xylene demethylation feeding the shared toluene pool (o-xylene + H₂ → toluene + methane — the real single-methyl-loss pathway, not a lumped double-dealkylation straight to benzene, so this is a genuine reactions-in-series-and-parallel network, not three independent reactions) both stay illustrative order-of-magnitude Arrhenius parameters, tuned only to land in the source's cited ~98% main-reaction selectivity at 2:1 H₂:toluene — honestly labeled as such, not fabricated citations. H₂S/CO₂/ammonia ride through the flowsheet as inert impurities — exactly what 'impurity tolerance' means in a mass balance, not a special chemistry path. A real fired heater (rigorous combustion stoichiometry + stack energy balance, the same unit op used in the naphtha-steam-cracker example) brings the mixed fresh-plus-recycle feed up to reaction temperature — not a feed specified as already at 950 K. The H₂-rich high-pressure separator vapor is a real recycle loop: fresh feed joins recycled gas at a mixer ahead of the furnace, and past the separator the vapor splits into a recompressed recycle (85%) and a genuine purge (15%) that bleeds off the methane a single pass could never consume — the same role a purge plays in the ammonia-synthesis-loop example. Quench, the high-pressure separator, and a letdown-valve-plus-flash stabilizer strip dissolved H₂/methane before the aromatics split. The reactor is genuinely **adiabatic**, not isothermal: temperature is a real energy-balance-coupled ODE state (the adiabatic PFR mode) that rises along the reactor length as the exothermic reactions release heat — starting from the furnace's real ~900 K outlet (in the source's cited 894-922 K inlet range) and climbing as the reaction proceeds, self-consistent to within ~0.01% of the reactor's own required external duty (which should be ~0 for a truly adiabatic reactor — verified, not assumed). One honest simplification remains: the final benzene/toluene/xylene-derived-benzene split still uses a shortcut separator rather than a rigorous column, for the same documented reason as hda-toluene-dealkylation: trace H₂/methane corrupts the column's wide-boiling initial-guess classification, and a fix narrowed to just that seed broke the shipped methanol-synthesis example when tried previously.

Which thermodynamic method does it use?

The PENG-ROBINSON property package, over toluene, o_xylene, benzene, biphenyl, h2, methane, h2s, co2, ammonia — already selected. You can switch the method on the canvas before running.

Which unit operations are in the flowsheet?

It chains MIX, Furnace, RX, Quench, 2× SEP2, Recycle Split, Recycle COMP, Letdown, COL. 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. Thermal (non-catalytic) HDA — broad aromatics feedstock runs entirely in your browser on MaximaLabs — free, no install, no license. Open the model to load it live and run the deterministic solver.