GASCO
VC Fired Heater — API 560 Design Suite
Engineering Calculation & Estimation Tool  ·  Internal Use Only
v16.0
Verify all results before submission

VC Fired Heater

Vertical Cylindrical Fired Heater — API 560 Design Modules

Version 1.0  ·  M14 Combustion  ·  M15 Radiant Section (Lobo-Evans)  ·  M6 Tube Wall (API 530)
Governing codes: API STD 560, 5th Ed. + Add.1 (2021) + Add.2 (2023)  ·  API RP 530, 7th Ed. (2015)  ·  GPSA Engineering Data Book, 14th Ed.  ·  Internal Evaluation Build — verify all results against primary standards before use in proposals

▶ STANDARD: | ★=HTRI/FRNC-5
▲ REAL:

M14 — Combustion Stoichiometry & Heat Balance

Calculates flue gas composition, MW, LHV, thermal efficiency and heat balance from fuel composition and combustion conditions. Outputs feed M15.

API 560 §6.2 · GPSA 14th Ed. Ch.27
JANAF Tables (Chase, NIST 1998)
Borger, Hyd.Proc. Sept 1994
⚙ Design Case Studies — Click any case to load all parameters
All module inputs (M1, M14, M15, M10, M9) cascade automatically. Use M25 → Run All Cases for side-by-side comparison table.
A — Standard Hot-Oil (Therminol 66, 155→220°C) · API 560 §6 sizing basis
B — HTRI Xfh / FRNC-5 PC / HYSIS Calibrated · Therminol 66 · 28–40 MW range
★ = HTRI Xfh v7 + FRNC-5 PC v3.1 calibrated  ·  ⚠ = h/D below API 560 §6 minimum 1.5 — review geometry
C — Real Project Cases · Loaded from Client Datasheets (Gasco supplied) · As-built geometry
Sources: Q13946=50191C-DS-03-MA-5044 · P30339B=14-01-34-1607 · Q9446=2971-5-13-0061 · Q10539N=MAEE012-S0207-0235 · B-0105=Methanex · Waitsia=45954-DS-1400-M-62320
Project identification
Process requirements
Site conditions
Project summary
Heat & mass balance
Design basis — applicable codes & standards
Fuel composition (mole %)
Total100.00 % ✓
Quick presets — click to load, then adjust mole %:
Enter mole percentages. Click Normalise if total ≠ 100%.
Combustion conditions
Reference temperature T_ref = 15°C (LHV/net calorific value basis per API 560 §6.2.3)
Thermal Efficiency (LHV basis) — API 560 §6.2.3
Flue gas composition (mole %)
Heat balance (per kg fuel)
Combustion parameters — key outputs for M15
Fuel consumption & flue gas quantities
▦ Adiabatic Flame Temperature & Sulphur / Acid Dew Point
AFT computed by iterative energy balance on flue gas (JANAF Cp polynomials). SO₂ from H₂S combustion; H₂SO₄ ADP via Verhoff-Bournstein (1971). Stack temperature check against ADP margin.
References & equations ▾
LHV = Σ(yᵢ × ΔHc_i) / MW_fuel [kJ/kg fuel]

Standard molar enthalpies of combustion (NIST/JANAF): CH₄ = 802,625; C₂H₆ = 1,427,820; C₃H₈ = 2,043,400; n-C₄H₁₀ = 2,658,000; H₂ = 241,818; CO = 282,989; H₂S = 518,000 kJ/kmol. Source: Chase, M.W., NIST-JANAF Thermochemical Tables, 4th Ed., J.Phys.Chem.Ref.Data, Monograph No.9, ACS/AIP, 1998.

Cp_fg = Σ(yᵢ × Cp_i) / MW_fg [kJ/kgK] where Cp_i are JANAF polynomial fits in T[K]

JANAF polynomial fits (kJ/kmol·K): CO₂: 37.14+0.0294T−8.74×10⁻⁶T²; H₂O: 32.22+0.00193T+1.06×10⁻⁵T²; N₂: 28.90+0.00167T; O₂: 29.52+0.00406T; Ar: 20.80 (monatomic). Source: Perry's Chemical Engineers' Handbook, 8th Ed., Table 2-172 (Green & Perry, McGraw-Hill, 2008). Also: GPSA Engineering Data Book, 14th Ed., Ch.27, Fig.27-5.

η_thermal = (LHV + Q_air_sensible − Q_stack − Q_rad_loss) / LHV × 100 [%]

Efficiency on LHV (net calorific value) basis per API STD 560, 5th Ed. §6.2.3. Radiation + casing heat loss = 1.5% of heat release for well-insulated heaters (API 560 §6.2). Source: API STD 560, 5th Edition (Feb 2016) + Addendum 1 (2021) + Addendum 2 (Dec 2023). Also: Borger, G.G., 'Approximate Formulae for Fired Heater Design', Hydrocarbon Processing, September 1994, pp.100-106.

Dry air composition: O₂=20.95%, N₂=78.09%, Ar=0.93% (mole%, dry basis)

Standard dry air composition per GPSA Engineering Data Book, 14th Ed., Fig.27-1. Molecular weight of dry air = 28.97 kg/kgmol. Actual air flow = stoichiometric × (1 + EA/100). Flue gas MW and composition calculated by atom balance on fuel + air.

Thermal inputs — from M14
Geometry — quick presets & design variables
Quick size presets (click to load)
Pitch = 2.00×OD (minimum per API 560)
▮ Flow Direction — Radiant Section
Counter-current: maximum tube metal T at bottom zone (near burners) where ΔT is highest. This is the conservative design case — API 560 preferred for stress/creep calculations.
Iterate pitch, N & L until all checks pass and Q_R ≥ Q_required. Use M10 Burner module to check flame clearances.
Calculating…
Plan view — tube & burner layout
● Tubes  •  TCD  •  Burner circle  •  D_int
Elevation — radiant section geometry
Geometry computed
Lobo-Evans radiation model
Heat transfer results
Flux & volumetric checks
API 560 / IS 13705 Compliance Checks
Radiant Section Design Summary
Cross-check: Area method vs Lobo-Evans
▦ HTRI Xref-Style Radiant Zone Gradient Analysis
Zone-by-zone temperature, heat flux and tube metal temperature profile using fired heater methodology. The radiant section is divided into N vertical zones from burner floor (bottom, Zone 1) to bridgewall arch (top, Zone N). Flue gas temperature declines from T_flame → T_bridgewall. Local Lobo-Evans radiation is applied at each zone. Both co-current and counter-current process fluid routings are solved and compared.
References & equations ▾
Q_R = σ · (α · A_cp) · F · (T_g⁴ − T_w⁴) [W] — Lobo-Evans fundamental equation

Lobo, W.E. & Evans, J.E., 'Heat Transfer in the Radiant Section of Petroleum Heaters', Transactions AIChE, Vol.35, 1939, pp.743-778. The foundational reference for all fired heater radiant section calculation. σ = 5.67×10⁻⁸ W/m²K⁴ (Stefan-Boltzmann). Tube absorption factor α = 0.97 for single tube row against refractory wall, pitch/diameter = 2.0, from Wimpress, R.N., 'How to Rate Fired Process Heaters', Hydrocarbon Processing, October 1963, pp.115-126.

A_cp = N_tubes × pitch × L_eff [m²] Cold plane area (tubes in circular arrangement)

Cold plane area = projected tube bank area (plane through tube centrelines). For circular VC arrangement the cold plane is the developed cylindrical surface. Effective (absorbed) area = α × A_cp where α from Wimpress (1963) effectiveness chart. Source: GPSA Engineering Data Book, 14th Ed., Chapter 13, Fig.13-2.

L_m = 0.60 × D_int [m] Mean beam length for vertical cylindrical heater

Mean beam length formula for cylindrical firebox per Wimpress, R.N. (1963). The mean beam length represents the characteristic path length of radiation through the gas volume. For other configurations: box heater = 0.60×(V_rad)^(1/3); general = 3.5×V/A_surface. Reference also in Hottel, H.C. & Sarofim, A.F., 'Radiative Heat Transfer', McGraw-Hill, 1967, Chapter 6.

ε_g ≈ 0.53·(1−exp(−1.4·PL^0.7))·(1100/T_g)^0.45 — curve-fit of Lobo-Evans emissivity

Analytical approximation of the Lobo-Evans (1939) flue gas emissivity chart. P = sum of CO₂+H₂O partial pressures [atm]. Temperature correction exponent per Borger (1994). For rigorous design, use the Hottel-Sarofim weighted-sum-of-grey-gases (WSGG) model as implemented in rigorous zone model software (rigorous zone model software technical literature). The WSGG model accounts for spectral band overlap between CO₂ and H₂O absorption bands across the full temperature range 500–2000 K.

F = ε_g / [ε_g + (1−ε_g) / (1 + A_w/αA_cp)] — Mekler-Fairall exchange factor

Overall radiant exchange factor accounting for gas emissivity and refractory re-radiation. Mekler, L. & Fairall, R.S., 'Better Heater Performance by Improved Design', Petroleum Refiner, June 1952, Vol.31, pp.101-110. Also tabulated in GPSA Engineering Data Book, 14th Ed., Ch.13, Fig.13-3. This is equivalent to the Hottel crossed-string formula for a single enclosure with partially absorbing gas.

h/D = H_rad / D_int ≤ 2.75 (API 560 §6.3.5 — mandatory)

API STD 560, 5th Edition, February 2016, §6.3.5: 'Vertical cylindrical heaters shall be designed with a maximum h/D of 2.75 [minimum 1.50]'. This ensures adequate combustion volume and prevents excessive height that would create temperature gradients requiring multi-zone analysis. h = H_rad; D = D_int (internal firebox diameter).

q_max ≤ 1.5 × q_avg (API 560 §6.1.4 — circumferential flux non-uniformity)

API STD 560, 5th Ed. §6.1.4: maximum heat flux density must not exceed 1.5× the average. F_max ≈ 1.30 for single tube row, s/d=2.0 (Wimpress 1963 circumferential non-uniformity chart). Higher F_max (1.4–1.5) occurs at tight pitch — confirm with burner supplier for flame impingement assessment. Also see Jegla, Z. et al., 'Numerical Analysis of Heat Transfer in Radiant Section of Fired Heaters', Chemical Engineering Transactions, Vol.39, 2014.

Tube specification
Design metal temperature (DMT) inputs
Thickness allowances (API 530 §4.4)
⚠ API 530 design life = 100,000 hours (§4.2.3). Schedule table shown is for NPS 6" (OD 168.3mm). Scale wall thickness for other NPS sizes using API 530 Eq.3.1/3.2 with the correct OD.
Calculating…
DMT build-up (API 530 §4.2)
T_DMT = T_fluid_max + ΔT_film + ΔT_wall + ΔT_safety
Thickness calculation steps
NPS 6" Schedule Selection — ASME B36.10M (OD = 168.3 mm, all schedules)
Schedule Wall (mm) ID (mm) Wt (kg/m) Status Notes

Source: ASME B36.10M-2018 (CS/alloy) and B36.19M-2018 (stainless). Selected schedule = first standard schedule with t_nominal ≥ t_specified.

Material allowable stresses — API 530 Annex B
Temp (°C) σ_a elastic (MPa) σ_r 100kh creep (MPa) Notes
API 530 / API 560 Compliance Checks
References & equations ▾
T_DMT = T_fluid_max + q_avg/h_i + q_avg·t_wall/k_metal + ΔT_safety [°C]

API RP 530, 7th Ed. (2015) §4.2.2. Minimum safety adder = 15°C. Inside film h_i from Dittus-Boelter: Nu = 0.023·Re⁰·⁸·Pr⁰·⁴ (heating), valid Re>10,000 (Bergman, T.L. et al., 'Fundamentals of Heat and Mass Transfer', 7th Ed., Wiley, 2011, Eq.8.60). Tube wall k_metal: CS A106-B = 45 W/mK; Cr-Mo alloys (T11/T22/T9) = 30 W/mK; austenitic SS (304H) = 16 W/mK (ASME BPVC Sec.II Part D).

t_min_elastic = P·d_o / (2·σ_a·f_t + 0.8·P) [mm] — API 530 Eq.3.1 (elastic design)

API RP 530, 7th Ed. (2015) §4.3.1. Identical formulation in ISO 13704:2021. Weld joint factor f_t = 1.0 for seamless tubes (mandatory per API 560 §7.1.4). σ_a = allowable stress = min(Yield_min/1.5, UTS_min/3) from API 530 Annex B material stress curves. Units: P in MPa, d_o in mm → t in mm. For P in kPa: divide by 1000 before substituting.

t_min_creep = P·d_o / (2·σ_r) [mm] — API 530 Eq.3.2 (creep-rupture design)

API RP 530, 7th Ed. (2015) §4.3.2. Applies when T_DMT exceeds the elastic limit temperature of the alloy. σ_r = minimum 100,000-hour rupture strength at T_DMT from API 530 Annex B rupture stress curves. Design life = 100,000 hours per §4.2.3. Governing thickness = max(t_el, t_cr). For remaining life / retirement thickness assessment see API 530 Annex A (Larson-Miller parameter).

t_specified = t_min + CA + t_min × (TOL/100) [mm] — API 530 §4.4

Corrosion allowance per API 560 §7.1.2 and process corrosion assessment. Mill undertolerance 12.5% for ASTM A106/A213/A335 seamless per ASTM A530 standard tolerances. Selected schedule must have t_nominal ≥ t_specified. Tube scheduling per ASME B36.10M-2018 (CS/alloy steel) and B36.19M-2018 (stainless). All material specifications are seamless construction only (API 560 §7.1.4 — no ERW or SAW tubes permitted in radiant section).