📑 Table of Contents (Tap to view sections)
1. Geometric Foundations of Partially Filled Horizontal Vessels
In process engineering design and facility operations, horizontal cylindrical vessels—including three-phase production separators, flare knock-out drums, reflux accumulators, and liquid storage bullets—rarely operate completely full. Instead, they operate with a gas-liquid interface located at some liquid level height $h$ from the bottom invert.
Calculating the exact liquid volume $V(h)$, the remaining vapor disengagement volume $V_{vapor}(h)$, the liquid cross-sectional area $A_L(h)$, and the vapor disengagement interfacial area $A_{int}(h)$ is required for:
- Calibrating differential pressure (DP) cell transmitters and guided-wave radar (GWR) level instruments to read true engineering volumes (barrels, cubic meters, or gallons) rather than simple linear percentage height.
- Evaluating true liquid residence and retention times for oil-water gravity settling (API Spec 12J).
- Determining effective vapor space velocity to verify the Souders-Brown droplet settling criteria across the length of the vessel.
Geometric volume formulations comply with ASME Boiler & Pressure Vessel Code (BPVC) Section VIII, Division 1, API Specification 12J (Oil and Gas Separators), and GPSA Engineering Data Book (Section 7).
2. Exact Analytical Derivations for Cylindrical Shells
Consider a horizontal cylinder of inside radius $R = D/2$ and cylindrical shell length $L$, filled with liquid to height $h$ ($0 \le h \le D$).
2.1 Cross-Sectional Area of a Circular Segment
The cross-sectional area of liquid $A_L(h)$ is obtained by integrating circular slices or via trigonometry of the circular sector:
The subtended half-angle $\theta$ (in radians) from the vertical centerline to the liquid surface contact point is:
$$\cos\theta = \frac{R - h}{R} = 1 - \frac{h}{R} \implies \theta = \arccos\left(1 - \frac{2h}{D}\right)$$The area of the circular segment occupied by liquid is:
$$A_L(h) = R^2 (\theta - \sin\theta \cos\theta) = R^2 \left[\arccos\left(1 - \frac{h}{R}\right) - \left(1 - \frac{h}{R}\right)\sqrt{1 - \left(1 - \frac{h}{R}\right)^2}\right]$$Expanding in terms of vessel diameter $D$ and liquid level $h$:
$$A_L(h) = \frac{D^2}{4} \arccos\left(1 - \frac{2h}{D}\right) - \left(\frac{D}{2} - h\right)\sqrt{h(D - h)}$$The total cylinder cross-sectional area is $A_{total} = \frac{\pi D^2}{4}$, and the remaining vapor cross-sectional area is:
$$A_v(h) = A_{total} - A_L(h)$$2.2 Shell Liquid Volume
For the straight cylindrical portion of length $L$, the partial volume is simply the cross-sectional area multiplied by length:
$$V_{shell}(h) = A_L(h) \cdot L$$2.3 Interfacial Surface Area (Vapor-Liquid Contact Width)
The horizontal chord width $W(h)$ of the liquid surface at height $h$ is given by:
$$W(h) = 2 \sqrt{R^2 - (R - h)^2} = 2\sqrt{h(D - h)}$$The interfacial gas-liquid surface area available for vapor disengagement is:
$$A_{int}(h) = W(h) \cdot L = 2 L \sqrt{h(D - h)}$$Maximum interfacial surface area occurs at the centerline ($h = D/2 = R$), where $W_{max} = D$ and $A_{int, max} = D \cdot L$.
3. Formed Head Volumes (2:1 Ellipsoidal, Hemispherical & ASME F&D)
Pressure vessels terminate in formed dished heads that contribute significant volume, especially in vessels with low $L/D$ ratios ($L/D < 3$). The three most common ASME head geometries are:
| Head Geometry | Depth Ratio ($b/D$) | Total Volume (Both Heads) | Application & Characteristics |
|---|---|---|---|
| 2:1 Ellipsoidal (ASME) | $b = D/4 = 0.25$ | $V_{2heads} = \frac{\pi D^3}{12} \approx 0.2618 D^3$ | Industry standard for medium to high pressure ($P > 15\,\text{bar}$). |
| Hemispherical | $b = D/2 = 0.50$ | $V_{2heads} = \frac{\pi D^3}{6} \approx 0.5236 D^3$ | Highest pressure capability; minimum wall thickness; high fabrication cost. |
| ASME Torispherical (F&D) | $b \approx 0.169$ | $V_{2heads} \approx 0.177 D^3$ | Economical for low-pressure vessels and atmospheric tanks ($P < 15\,\text{bar}$). |
3.1 Partial Volume of 2:1 Ellipsoidal Heads
For a pair of 2:1 ellipsoidal heads in a horizontal vessel, the partial volume as a function of dimensionless liquid depth $y = h/D$ ($0 \le y \le 1$) is derived from the integration of horizontal elliptical slices:
$$V_{heads}(h) = V_{2heads, total} \cdot \left[3 y^2 - 2 y^3\right] = \frac{\pi D^3}{12} \cdot \left[3 \left(\frac{h}{D}\right)^2 - 2 \left(\frac{h}{D}\right)^3\right]$$This cubic formulation is exact for horizontal 2:1 ellipsoidal heads. At $h = D/2$ ($y = 0.5$), $3(0.5)^2 - 2(0.5)^3 = 0.75 - 0.25 = 0.50$, correctly confirming that the heads are exactly half-full at the centerline.
3.2 Total Vessel Partial Volume Equation
The complete partial liquid volume of a horizontal separator with formed 2:1 ellipsoidal heads is the sum of the cylindrical shell and both heads:
$$V_{total}(h) = V_{shell}(h) + V_{heads}(h) = A_L(h) \cdot L + \frac{\pi D^3}{12} \left[3 \left(\frac{h}{D}\right)^2 - 2 \left(\frac{h}{D}\right)^3\right]$$4. Dimensionless Fill Fraction Reference Table
Process engineers frequently utilize dimensionless fill ratios ($h/D$ vs. $V/V_{total}$) for rapid field estimation and DCS level configuration:
| Level Fraction ($h/D$) | Shell Area Fraction ($A_L/A_T$) | Head Volume Fraction ($V_{head}/V_{head,T}$) | Total Volume ($L/D = 3.0$) | Operating Significance |
|---|---|---|---|---|
| $0.05$ (5%) | $0.0187$ | $0.0073$ | $0.0175$ (1.75%) | Low-Low Level Trip (LAHH/LALL cutoff); minimal volume |
| $0.15$ (15%) | $0.0941$ | $0.0608$ | $0.0908$ (9.08%) | Typical Low Liquid Level (LLL); bottom sump buffer |
| $0.30$ (30%) | $0.2523$ | $0.2160$ | $0.2487$ (24.87%) | Low Normal Operating Level (NLL) in gas-dominated drums |
| $0.50$ (50%) | $0.5000$ | $0.5000$ | $0.5000$ (50.00%) | Vessel Centerline: Maximum surface area $W = D$ |
| $0.70$ (70%) | $0.7477$ | $0.7840$ | $0.7513$ (75.13%) | High Liquid Level (HLL); vapor disengagement begins choking |
| $0.85$ (85%) | $0.9059$ | $0.9393$ | $0.9092$ (90.92%) | High-High Liquid Level Trip (HHLL); shutdown trigger |
Notice that at 15% liquid level height ($h/D = 0.15$), the vessel holds only 9.08% of total volume! Conversely, going from 40% to 60% height ($20\%$ change in level) accounts for 27.3% of total volume. Assuming linear level-to-volume response causes massive errors in residence time and inventory control.
5. Horizontal Vapor Space Capacity & Gas Residence Time
In a horizontal vessel, gas enters one end and flows horizontally toward the outlet nozzle while liquid drops settle downward:
5.1 Vapor Passage Velocity
The average horizontal superficial gas velocity across the vapor disengagement space is:
$$v_g = \frac{Q_g}{A_v(h)} = \frac{Q_g}{A_{total} - A_L(h)}$$To avoid bulk interfacial liquid stripping and re-entrainment, $v_g$ must not exceed the critical re-entrainment velocity:
$$v_{crit} = K_{horiz} \sqrt{\frac{\rho_l - \rho_g}{\rho_g}}$$where $K_{horiz} \approx 0.07\text{--}0.09\,m/s$ ($0.23\text{--}0.30\,ft/s$) for vessels without demister pads, and $0.12\text{--}0.15\,m/s$ when an internal horizontal vane or mesh pad is installed.
5.2 Gas Residence Time
The time available for a droplet to fall from the top of the vessel to the liquid interface is the gas residence time:
$$t_{residence} = \frac{L_{eff}}{v_g} = \frac{L_{eff} \cdot A_v(h)}{Q_g}$$where $L_{eff}$ is the effective horizontal separation distance between the inlet distributor and the outlet mist extractor (typically $L_{eff} \approx 0.80\text{--}0.85 \times L_{T-T}$).
6. Step-by-Step Worked Engineering Example
A production separator has the following design dimensions and operating parameters:
- Inside Diameter: $D = 2.40\,\text{m}$ ($2400\,\text{mm}$ / 94.5 in)
- Shell Tangent-to-Tangent Length: $L = 7.20\,\text{m}$ ($L/D = 3.0$)
- Head Type: 2:1 Ellipsoidal heads on both ends
- Current Operating Level: $h = 0.96\,\text{m}$ from bottom invert
- Gas Volumetric Flow: $Q_g = 8.50\,\text{m}^3\text{/s}$ at operating conditions
- Liquid Flow: $Q_l = 85.0\,\text{m}^3\text{/h}$
Step 1: Calculate Dimensionless Fill Ratio
$$y = \frac{h}{D} = \frac{0.96}{2.40} = 0.400\quad (40.0\%)$$Step 2: Calculate Circular Cross-Sectional Areas
Total shell cross-sectional area:
$$A_{total} = \frac{\pi D^2}{4} = \frac{\pi \times (2.40)^2}{4} = 4.5239\,\text{m}^2$$Liquid cross-sectional area $A_L(h)$:
$$\cos\theta = 1 - \frac{2 \times 0.96}{2.40} = 1 - 0.80 = 0.20 \implies \theta = \arccos(0.20) = 1.3694\,\text{rad}\quad (78.46^\circ)$$ $$A_L = \frac{(2.40)^2}{4} \times 1.3694 - \left(\frac{2.40}{2} - 0.96\right) \sqrt{0.96 \times (2.40 - 0.96)}$$ $$A_L = 1.44 \times 1.3694 - (0.24) \times \sqrt{0.96 \times 1.44} = 1.9719 - 0.24 \times 1.1758 = 1.9719 - 0.2822 = 1.6897\,\text{m}^2$$Liquid area fraction: $A_L / A_{total} = 1.6897 / 4.5239 = 0.3735$ (37.35%).
Vapor cross-sectional area: $A_v = 4.5239 - 1.6897 = 2.8342\,\text{m}^2$ (62.65%).
Step 3: Calculate Component & Total Liquid Volumes
Shell liquid volume:
$$V_{shell} = A_L \times L = 1.6897\,\text{m}^2 \times 7.20\,\text{m} = 12.166\,\text{m}^3$$Formed heads liquid volume (2:1 ellipsoidal, both ends):
$$V_{heads} = \frac{\pi D^3}{12} \left[3 y^2 - 2 y^3\right] = \frac{\pi \times (2.40)^3}{12} \left[3(0.40)^2 - 2(0.40)^3\right]$$ $$V_{heads} = 3.6191 \times [3(0.16) - 2(0.064)] = 3.6191 \times [0.48 - 0.128] = 3.6191 \times 0.352 = 1.274\,\text{m}^3$$Total liquid inventory currently held in separator:
$$V_{total, liq} = 12.166 + 1.274 = 13.440\,\text{m}^3\quad (\approx 84.5\,\text{barrels})$$Step 4: Calculate Liquid Residence Time
$$t_{retention} = \frac{V_{total, liq}}{Q_l} = \frac{13.440\,\text{m}^3}{85.0\,\text{m}^3\text{/h}} = 0.1581\,\text{h} = 9.49\,\text{minutes}$$Per API Spec 12J, standard retention time for light to medium oil-water separation is 5 to 10 minutes. At 40% level ($h=0.96\,\text{m}$), the vessel comfortably satisfies retention criteria.
Step 5: Verify Vapor Space Velocity & Residence Time
Horizontal gas passage velocity:
$$v_g = \frac{Q_g}{A_v} = \frac{8.50\,\text{m}^3\text{/s}}{2.8342\,\text{m}^2} = 3.00\,\text{m/s}$$Gas residence time across effective separation distance ($L_{eff} = 0.85 \times 7.20 = 6.12\,\text{m}$):
$$t_{gas} = \frac{6.12\,\text{m}}{3.00\,\text{m/s}} = 2.04\,\text{seconds}$$7. Level Transmitter Calibration (Strapping Tables)
To program automated DCS/SCADA systems, process engineers translate the continuous $h/D$ relationship into a 21-point strapping table (0% to 100% in 5% increments). Modern smart radar and magnetostrictive transmitters can store the non-linear Strapping polynomial directly:
$$V_{liquid} = C_1 h^3 + C_2 h^2 + C_3 h$$Applying strapping tables prevents false high-level alarms during sudden production rate surges, smoothing out the response of downstream feed pumps.
8. ChemProCal Integration
Perform complete vessel geometry and partial volume calculations using:
- Separator Utilities: Instant partial cylinder and head volume computation for any level height.
- Horizontal Separator Sizing Tool: Complete multi-phase vessel sizing per API 12J.
- Vertical Separator Sizing Tool: Comprehensive vertical vessel geometry and nozzle calculations.
Separator Utilities
Apply this methodology directly in the ChemProCal calculator.
Open Calculator →Live Separator Souders-Brown Capacity Estimator
Adjust parameters below to test the methodology equations in real time before running full simulations: