Liquid ethylene (C₂H₄) is one of the highest-volume petrochemical feedstocks on the planet, yet its thermodynamic behavior at cryogenic temperatures makes accurate property estimation anything but trivial. Density alone can swing from 655 kg/m³ near the freezing point down to roughly 214 kg/m³ at the critical boundary — a shift of over 67% that directly governs vessel sizing, pressure relief design, and transfer-line hydraulics.

This calculator resolves that complexity into a single automated workflow. Provide a system temperature and a storage volume, and it returns the saturated liquid density, corresponding vapor pressure, total stored mass, specific gravity, specific volume, and an estimated latent heat of vaporization — all cross-referenced against high-accuracy NIST tabulated data.

Required Calculation Parameters

Before running the estimation, gather the following design inputs:

  • System Temperature — the operating temperature of the liquid ethylene, accepted in Kelvin (K), Celsius (°C), or Fahrenheit (°F). The valid thermodynamic range spans from the melting point at 104.0 K (−169.15 °C) to the critical temperature at 282.34 K (9.19 °C). Values outside this range are automatically clamped to the nearest boundary.
  • Temperature Unit Selection — choose between Kelvin, Celsius, or Fahrenheit. All internal computations are performed in Kelvin; conversion is handled automatically.
  • Liquid Volume — the volume of ethylene liquid stored in the vessel or pipeline segment, specified in cubic meters (m³), liters (L), US gallons (gal), or cubic feet (ft³).
  • Volume Unit — the measurement system for the liquid volume entry.
  • Density Unit — the preferred output unit for the saturated liquid density result: kg/m³, g/cm³, lb/ft³, or lb/gal.

Theoretical Foundation and Governing Equations

Saturated Liquid Density by Piecewise Linear Interpolation

The core density model is built on a 21-point lookup table derived from NIST Standard Reference Data for ethylene along the saturation curve. The data spans from 104.0 K to 282.34 K and captures the nonlinear relationship between temperature and liquid-phase density.

For any temperature $T$ falling between two adjacent data points $(T_i, \rho_i)$ and $(T_{i+1}, \rho_{i+1})$, the density is computed by linear interpolation:

$$\rho(T) = \rho_i + \frac{T - T_i}{T_{i+1} - T_i} \cdot (\rho_{i+1} - \rho_i)$$

This method avoids the error accumulation inherent in polynomial curve-fitting across the full temperature range, particularly near the critical point where $\frac{\partial\rho}{\partial T}$ increases sharply. By anchoring every segment to experimentally validated NIST values, the interpolation maintains sub-percent accuracy in all regions except the immediate vicinity of $T_c$.

Vapor Pressure via the Antoine Equation

The equilibrium vapor pressure of ethylene at a given temperature is calculated using the Antoine equation, one of the most widely used semi-empirical correlations in chemical thermodynamics:

$$\log_{10}(P) = A - \frac{B}{T + C}$$

For ethylene in this implementation, the Antoine constants are:

  • $A = 3.91382$
  • $B = 596.526$
  • $C = -16.47$

where $T$ is in Kelvin and $P$ is returned in bar. The Antoine equation is generally valid from approximately 119.8 K to 282.3 K for ethylene. This correlation is derived from regression of experimental vapor-pressure measurements and is standard practice as documented in Perry's Chemical Engineers' Handbook and the NIST Chemistry WebBook.

Latent Heat of Vaporization — Watson's Correlation

Estimating the energy required to vaporize liquid ethylene at temperatures away from the normal boiling point uses Watson's relation, a corresponding-states correlation first published by K. M. Watson in 1943:

$$\Delta H_v(T) = \Delta H_{v,\text{bp}} \cdot \left( \frac{T_c - T}{T_c - T_{\text{bp}}} \right)^{n}$$

The key parameters for ethylene are:

  • $\Delta H_{v,\text{bp}} = 482.0 \text{ kJ/kg}$ (latent heat at the normal boiling point of 169.4 K)
  • $T_c = 282.34 \text{ K}$ (critical temperature)
  • $T_{\text{bp}} = 169.4 \text{ K}$ (normal boiling point at 1 atm)
  • $n = 0.38$ (Watson's universal exponent)

The exponent 0.38, while sometimes substance-specific as shown by Viswanath and Kuloor (1967), provides engineering-grade accuracy for light hydrocarbons including ethylene. A critical property of this equation: as $T \to T_c$, the numerator approaches zero and $\Delta H_v \to 0$, correctly reflecting the physical reality that liquid and vapor phases become indistinguishable at the critical point.

Derived Properties

From the primary density result, several secondary properties are computed:

Total Mass — the product of liquid density and volume, converted to consistent SI units:

$$m = \rho \cdot V$$

Specific Volume — the reciprocal of density, expressed in liters per kilogram:

$$v = \frac{1}{\rho} \times 1000$$

Specific Gravity — the ratio of liquid ethylene density to the density of water at 4 °C (999.97 kg/m³):

$$SG = \frac{\rho_{\text{ethylene}}}{\rho_{\text{water}}}$$

Technical Specifications and Reference Data

The table below provides the complete NIST-derived dataset used for interpolation. These values represent saturated liquid ethylene at each temperature — i.e., the liquid in equilibrium with its own vapor at the corresponding saturation pressure.

Temperature (K)Temperature (°C)Density (kg/m³)Approx. Vapor Pressure (bar)Region
104.0−169.15655.45≈ 0.0001Near melting point
110.0−163.15647.85≈ 0.0004Deep cryogenic
120.0−153.15635.01≈ 0.003Deep cryogenic
130.0−143.15621.92≈ 0.015Deep cryogenic
140.0−133.15608.51≈ 0.056Cryogenic
150.0−123.15594.70≈ 0.16Cryogenic
160.0−113.15580.40≈ 0.40Cryogenic
170.0−103.15565.51≈ 1.08Near boiling point
180.0−93.15549.91≈ 2.03Low-pressure liquid
190.0−83.15533.43≈ 3.53Low-pressure liquid
200.0−73.15515.88≈ 5.72Moderate pressure
210.0−63.15496.99≈ 8.77Moderate pressure
220.0−53.15476.38≈ 12.8Moderate pressure
230.0−43.15453.53≈ 18.1Elevated pressure
240.0−33.15427.61≈ 24.7Elevated pressure
250.0−23.15397.16≈ 32.5High pressure
260.0−13.15359.10≈ 37.9Near-critical
270.0−3.15305.15≈ 43.8Near-critical
275.0+1.85268.50≈ 47.0Near-critical
280.0+6.85214.00≈ 49.7Near-critical
282.34+9.19214.2050.4Critical point

Key physical constants for ethylene (C₂H₄):

  • Molecular weight: 28.054 g/mol
  • Normal boiling point: 169.4 K (−103.75 °C) at 1 atm
  • Critical temperature: 282.34 K (9.19 °C)
  • Critical pressure: 50.4 bar (5.04 MPa)
  • Critical density: ~214.2 kg/m³
  • Melting point: 104.0 K (−169.15 °C)

Engineering Analysis and Real-World Application

How Temperature Governs Density — and Why It Matters

The dominant engineering insight from this calculator is the extreme sensitivity of liquid ethylene density to temperature as the system approaches the critical point. Between 104 K and 200 K — the typical deep cryogenic storage range — density drops at roughly 1.4 kg/m³ per degree Kelvin. This is manageable and predictable.

However, above 250 K, the rate of density change accelerates dramatically. Between 260 K and 280 K, density plummets from 359 to 214 kg/m³ — a collapse of 145 kg/m³ in just 20 K. For any process operating in this near-critical zone, even a modest temperature excursion can cause:

  • Rapid liquid expansion that overwhelms ullage space in storage tanks, potentially triggering pressure relief.
  • Hydraulic instability in transfer lines, as the fluid's compressibility increases by orders of magnitude.
  • Loss of net positive suction head (NPSH) at pump inlets, leading to cavitation.

Tank Sizing and Thermal Expansion

A standard engineering practice for cryogenic ethylene storage is to limit the liquid fill level to 85–90% of total vessel volume, preserving adequate vapor space for thermal expansion and boil-off management. The calculator's mass estimation directly supports this workflow.

For example, a 100 m³ vessel filled with saturated liquid ethylene at 170 K contains approximately 56,550 kg of product. If a process upset warms the liquid to 200 K, the same mass now requires roughly 109.7 m³ at the lower density of 515.88 kg/m³ — a volume increase of nearly 10%. This overflow scenario is a primary driver behind pressure relief valve sizing on ethylene storage tanks.

Interpreting the T / T_crit Ratio

The critical-point proximity indicator ($T / T_c$) provides an immediate sense of operational margin. Industry guidance generally considers:

  • Below 80% — the liquid is well within the cryogenic regime. Properties are stable and predictable. Standard correlations perform at their highest accuracy.
  • 80% to 95% — the liquid is entering the near-critical transition zone. Property gradients steepen, and small temperature changes produce large density shifts. Additional engineering caution is warranted.
  • Above 95% — the system is approaching or at the critical point. The distinction between liquid and vapor phases begins to vanish. Conventional liquid-handling assumptions may no longer apply.

Latent Heat and Energy Budgeting

The latent heat result ($\Delta H_v$) is essential for sizing boil-off gas (BOG) compressors, designing condenser duty in ethylene recovery trains, and calculating the refrigeration load needed to maintain storage temperature. At the normal boiling point, ethylene requires roughly 482 kJ/kg to fully vaporize — less than half the value for water, but still a significant energy penalty in large-scale operations.

As the liquid warms toward the critical temperature, $\Delta H_v$ falls steadily to zero. This means that heat leaking into a near-critical ethylene system does not produce a clean phase transition; instead, it causes a continuous, difficult-to-manage density decrease.

Frequently Asked Questions

Why does the calculator clamp temperatures outside the 104 K – 282.34 K range instead of extrapolating?

Below 104 K, ethylene exists as a solid, and the density-temperature relationship follows entirely different physics governed by crystalline lattice spacing rather than liquid intermolecular forces. Extrapolating the liquid-phase correlation into the solid region would produce meaningless results.

Above 282.34 K, ethylene becomes a supercritical fluid — a single-phase state where no liquid surface exists. The saturation curve terminates at the critical point, and concepts like "saturated liquid density" lose their physical definition. Clamping to the boundary and displaying a warning is the thermodynamically honest approach.

How accurate is the Antoine equation for ethylene vapor pressure at very low temperatures?

The Antoine equation begins to lose accuracy below approximately 120 K, where the vapor pressure of ethylene drops to fractions of a millibar. At these extremely low pressures, experimental measurement uncertainty is high, and the Antoine constants — regressed from data concentrated around the boiling point — struggle to capture the curvature of the $\ln P$ vs. $1/T$ relationship.

For engineering design work below 120 K, more sophisticated equations of state such as the Peng-Robinson or the Helmholtz-energy-explicit formulation by Jahangiri, Jacobsen, Stewart, and McCarty (1986) should be consulted. However, for the vast majority of practical storage and process applications, which occur between 150 K and 250 K, the Antoine correlation maintains errors well below 2%.

Can this calculator be used for ethylene-ethane mixtures or impure ethylene streams?

No. The density data and vapor-pressure correlation are strictly for pure ethylene along its saturation curve. Real-world ethylene streams from cracking units typically contain trace amounts of ethane, methane, propylene, and acetylene. Even small impurities shift the saturation envelope and alter both the bubble-point temperature and the liquid density.

For mixture calculations, rigorous thermodynamic models such as Soave-Redlich-Kwong (SRK) or Peng-Robinson equations of state with appropriate binary interaction parameters are required. Process simulation software like Aspen HYSYS or PRO/II handles these computations using the full multicomponent phase equilibrium framework.

Professional Conclusion

Accurate knowledge of liquid ethylene density is not a theoretical exercise — it directly determines the safety margin of cryogenic storage vessels, the sizing of transfer pumps and relief valves, and the energy balance of downstream recovery processes. Manual lookup from printed tables introduces interpolation error and consumes engineering time that would be better spent on analysis.

This automated estimator consolidates NIST-validated data with proven correlations (Antoine for vapor pressure, Watson for latent heat) into a single rapid computation. By eliminating the need for manual table lookups and unit conversions, it reduces the risk of transcription errors that can propagate through design calculations with costly consequences.