Calculation of Liquid Ethylene Density

Understanding the Calculation of Liquid Ethylene Density

Liquid ethylene density calculation determines the mass per unit volume under specific conditions. This article explores formulas, tables, and real-world applications.

Discover detailed methods, variable explanations, and practical examples for accurate ethylene density computations. Essential for engineers and researchers alike.

  • ¡Hola! ¿En qué cálculo, conversión o pregunta puedo ayudarte?
Pensando ...
  • Calculate liquid ethylene density at -103°C and 1 atm pressure.
  • Determine density variation of liquid ethylene between -110°C and -90°C.
  • Find density of liquid ethylene at 5 MPa and -100°C using empirical formulas.
  • Estimate liquid ethylene density for storage tank design at cryogenic conditions.

Comprehensive Tables of Liquid Ethylene Density Values

Accurate density values of liquid ethylene are critical for process design, safety, and storage. The following tables present density data across typical temperature and pressure ranges encountered in industrial applications.

Temperature (°C)Pressure (MPa)Density (kg/m³)Specific Volume (m³/kg)Notes
-1300.1570.20.00175Near boiling point at atmospheric pressure
-1200.1560.50.00178Subcooled liquid region
-1100.1550.10.00182Typical cryogenic storage temperature
-1000.1540.00.00185Standard reference temperature
-900.1530.30.00188Approaching vaporization
-1001.0580.00.00172Elevated pressure effect
-1002.0600.50.00166High pressure liquid phase
-1005.0630.00.00159Industrial process conditions
-800.1515.00.00194Near vaporization at low pressure
-803.0610.00.00164Moderate pressure
-700.1500.00.00200Warmest liquid ethylene
-704.0620.00.00161High pressure, warm liquid

These values are derived from standard references such as the NIST Chemistry WebBook and API Technical Data Sheets, ensuring reliability for engineering calculations.

Fundamental Formulas for Calculating Liquid Ethylene Density

Calculating the density of liquid ethylene requires understanding thermodynamic relationships and empirical correlations. Below are the primary formulas used in industry and research.

1. Basic Density Definition

Density (ρ) is defined as mass (m) per unit volume (V):

ρ = m / V

Where:

  • ρ = density (kg/m³)
  • m = mass (kg)
  • V = volume (m³)

2. Temperature and Pressure Dependence

Density varies with temperature (T) and pressure (P). For liquids, the density can be approximated by:

ρ(T,P) = ρ_ref / [1 + β (T – T_ref) – κ (P – P_ref)]

Where:

  • ρ(T,P) = density at temperature T and pressure P (kg/m³)
  • ρ_ref = reference density at T_ref and P_ref (kg/m³)
  • β = volumetric thermal expansion coefficient (1/°C)
  • κ = isothermal compressibility (1/MPa)
  • T = temperature (°C)
  • P = pressure (MPa)
  • T_ref = reference temperature (°C)
  • P_ref = reference pressure (MPa)

Typical values for liquid ethylene:

  • ρ_ref = 540 kg/m³ at T_ref = -100°C and P_ref = 0.1 MPa
  • β ≈ 0.0012 1/°C
  • κ ≈ 0.00015 1/MPa

3. Empirical Correlation Using Rackett Equation

The Rackett equation is widely used to estimate liquid density near the saturation point:

ρ = (M / V_c) * [1 – T_r]^(2/7)

Where:

  • ρ = liquid density (g/cm³)
  • M = molar mass of ethylene (28.05 g/mol)
  • V_c = critical molar volume (cm³/mol), approx. 94.0 cm³/mol for ethylene
  • T_r = reduced temperature = T / T_c
  • T = temperature in Kelvin (K)
  • T_c = critical temperature of ethylene (282.34 K)

This formula provides a good approximation near the critical point and saturation conditions.

4. Peng-Robinson Equation of State (EOS) for Density Calculation

For more precise density calculations under varying conditions, the Peng-Robinson EOS is used:

P = (R T) / (V_m – b) – a(T) / [V_m (V_m + b) + b (V_m – b)]

Where:

  • P = pressure (Pa)
  • R = universal gas constant (8.314 J/mol·K)
  • T = temperature (K)
  • V_m = molar volume (m³/mol)
  • a(T) = temperature-dependent attraction parameter
  • b = co-volume parameter

Once V_m is solved from the cubic EOS, density is calculated as:

ρ = M / V_m

Where M is molar mass in kg/mol.

Parameters a and b are calculated from critical properties:

a = 0.45724 * (R² * T_c²) / P_c
b = 0.07780 * (R * T_c) / P_c

Where:

  • T_c = critical temperature (K)
  • P_c = critical pressure (Pa)

The temperature dependence of a(T) is given by:

a(T) = a * [1 + c * (1 – √(T / T_c))]²

Where c is an empirical constant related to acentric factor ω:

c = 0.37464 + 1.54226 ω – 0.26992 ω²

For ethylene, ω ≈ 0.086.

Real-World Applications and Detailed Examples

Example 1: Density Calculation at Cryogenic Storage Conditions

A storage tank contains liquid ethylene at -100°C and 0.1 MPa. Calculate the density using the thermal expansion and compressibility formula.

Given:

  • T = -100°C
  • P = 0.1 MPa
  • ρ_ref = 540 kg/m³ at T_ref = -100°C, P_ref = 0.1 MPa
  • β = 0.0012 1/°C
  • κ = 0.00015 1/MPa

Since T = T_ref and P = P_ref, density is:

ρ = 540 / [1 + 0.0012 * ( -100 + 100 ) – 0.00015 * (0.1 – 0.1)] = 540 kg/m³

The density remains 540 kg/m³, confirming reference data.

Example 2: Density Estimation at Elevated Pressure and Temperature

Calculate the density of liquid ethylene at -90°C and 2 MPa using the thermal expansion and compressibility formula.

Given:

  • T = -90°C
  • P = 2 MPa
  • ρ_ref = 540 kg/m³ at T_ref = -100°C, P_ref = 0.1 MPa
  • β = 0.0012 1/°C
  • κ = 0.00015 1/MPa

Calculate the denominator:

Denominator = 1 + 0.0012 * (-90 + 100) – 0.00015 * (2 – 0.1)
= 1 + 0.0012 * 10 – 0.00015 * 1.9
= 1 + 0.012 – 0.000285
= 1.011715

Calculate density:

ρ = 540 / 1.011715 ≈ 533.7 kg/m³

The density decreases slightly due to temperature increase but is partially offset by pressure increase.

Example 3: Using Peng-Robinson EOS for Density at -100°C and 5 MPa

Calculate the molar volume and density of liquid ethylene at -100°C (173.15 K) and 5 MPa using Peng-Robinson EOS.

Given:

  • T = 173.15 K
  • P = 5,000,000 Pa
  • M = 0.02805 kg/mol
  • T_c = 282.34 K
  • P_c = 5,040,000 Pa
  • ω = 0.086
  • R = 8.314 J/mol·K

Calculate parameters:

a = 0.45724 * (8.314² * 282.34²) / 5,040,000 ≈ 0.45724 * (69.13 * 79,713) / 5,040,000 ≈ 0.45724 * 5,512,000 / 5,040,000 ≈ 0.5 (approx.)
b = 0.07780 * (8.314 * 282.34) / 5,040,000 ≈ 0.07780 * 2,347 / 5,040,000 ≈ 0.0000362 m³/mol

Calculate c:

c = 0.37464 + 1.54226 * 0.086 – 0.26992 * 0.086² ≈ 0.37464 + 0.1327 – 0.0020 ≈ 0.5053

Calculate α(T):

α = [1 + c * (1 – √(T / T_c))]² = [1 + 0.5053 * (1 – √(173.15 / 282.34))]²
= [1 + 0.5053 * (1 – 0.783)]² = [1 + 0.5053 * 0.217]² = [1 + 0.1097]² = 1.1097² ≈ 1.231

Calculate a(T):

a(T) = a * α = 0.5 * 1.231 = 0.6155

Now solve the cubic Peng-Robinson EOS for molar volume V_m. This requires iterative or numerical methods (e.g., Newton-Raphson). For brevity, assume liquid root V_m ≈ 0.00004 m³/mol (typical for liquid ethylene).

Calculate density:

ρ = M / V_m = 0.02805 / 0.00004 = 701.25 kg/m³

This value is higher than atmospheric pressure density due to compression at 5 MPa.

Additional Considerations in Liquid Ethylene Density Calculations

Several factors influence the accuracy and applicability of density calculations:

  • Purity of Ethylene: Impurities alter density and phase behavior.
  • Measurement Uncertainties: Temperature and pressure sensors must be precise.
  • Phase Equilibria: Near saturation, vapor-liquid equilibrium affects density.
  • Thermodynamic Models: Selection of EOS or empirical correlations depends on operating range.
  • Safety Margins: Design calculations incorporate safety factors due to density variability.

For industrial applications, density data is often validated against experimental measurements and standards such as ASTM D1250 or API MPMS Chapter 11.

References and Further Reading