CHEN_EDELSTEIN
This function estimates the two-phase flow boiling heat transfer coefficient using the Chen-Edelstein correlation for saturated tube flow. It blends nucleate-boiling and convective mechanisms using empirical enhancement and suppression factors.
The model structure is:
h_{tp} = S h_{nb} + F h_{sp,l}
where the single-phase liquid term and nucleate-boiling term are combined to represent mixed boiling behavior across a range of flow conditions.
Excel Usage
=CHEN_EDELSTEIN(m, x, D, rhol, rhog, mul, mug, kl, Cpl, Hvap, sigma, dPsat, Te)
m(float, required): Mass flow rate (kg/s).x(float, required): Quality at the tube interval (dimensionless).D(float, required): Tube diameter (m).rhol(float, required): Liquid density (kg/m^3).rhog(float, required): Gas density (kg/m^3).mul(float, required): Liquid viscosity (Pa*s).mug(float, required): Gas viscosity (Pa*s).kl(float, required): Liquid thermal conductivity (W/m/K).Cpl(float, required): Liquid heat capacity (J/kg/K).Hvap(float, required): Heat of vaporization (J/kg).sigma(float, required): Surface tension (N/m).dPsat(float, required): Saturation pressure difference (Pa).Te(float, required): Excess wall temperature (K).
Returns (float): Heat transfer coefficient (W/m^2/K), or an error message if invalid.
Example 1: Example from reference
Inputs:
| m | x | D | rhol | rhog | mul | mug | kl | Cpl | Hvap | sigma | dPsat | Te |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0.106 | 0.2 | 0.0212 | 567 | 18.09 | 0.000156 | 0.00000711 | 0.086 | 2730 | 200000 | 0.02 | 100000 | 3 |
Excel formula:
=CHEN_EDELSTEIN(0.106, 0.2, 0.0212, 567, 18.09, 0.000156, 0.00000711, 0.086, 2730, 200000, 0.02, 100000, 3)
Expected output:
3289.06
Example 2: Higher quality with larger diameter
Inputs:
| m | x | D | rhol | rhog | mul | mug | kl | Cpl | Hvap | sigma | dPsat | Te |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0.18 | 0.45 | 0.03 | 820 | 12 | 0.00021 | 0.000011 | 0.11 | 3000 | 210000 | 0.022 | 90000 | 6 |
Excel formula:
=CHEN_EDELSTEIN(0.18, 0.45, 0.03, 820, 12, 0.00021, 0.000011, 0.11, 3000, 210000, 0.022, 90000, 6)
Expected output:
4517.15
Example 3: Low quality at small diameter
Inputs:
| m | x | D | rhol | rhog | mul | mug | kl | Cpl | Hvap | sigma | dPsat | Te |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0.08 | 0.15 | 0.018 | 930 | 7.5 | 0.00025 | 0.000013 | 0.13 | 3400 | 230000 | 0.028 | 70000 | 2.5 |
Excel formula:
=CHEN_EDELSTEIN(0.08, 0.15, 0.018, 930, 7.5, 0.00025, 0.000013, 0.13, 3400, 230000, 0.028, 70000, 2.5)
Expected output:
5167.53
Example 4: Mid-range inputs
Inputs:
| m | x | D | rhol | rhog | mul | mug | kl | Cpl | Hvap | sigma | dPsat | Te |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0.12 | 0.3 | 0.024 | 740 | 10 | 0.00019 | 0.000009 | 0.095 | 2600 | 195000 | 0.019 | 85000 | 4.5 |
Excel formula:
=CHEN_EDELSTEIN(0.12, 0.3, 0.024, 740, 10, 0.00019, 0.000009, 0.095, 2600, 195000, 0.019, 85000, 4.5)
Expected output:
3944.98
Python Code
Show Code
from ht.boiling_flow import Chen_Edelstein as ht_Chen_Edelstein
def Chen_Edelstein(m, x, D, rhol, rhog, mul, mug, kl, Cpl, Hvap, sigma, dPsat, Te):
"""
Compute the Chen-Edelstein boiling heat transfer coefficient.
See: https://ht.readthedocs.io/en/latest/ht.boiling_flow.html
This example function is provided as-is without any representation of accuracy.
Args:
m (float): Mass flow rate (kg/s).
x (float): Quality at the tube interval (dimensionless).
D (float): Tube diameter (m).
rhol (float): Liquid density (kg/m^3).
rhog (float): Gas density (kg/m^3).
mul (float): Liquid viscosity (Pa*s).
mug (float): Gas viscosity (Pa*s).
kl (float): Liquid thermal conductivity (W/m/K).
Cpl (float): Liquid heat capacity (J/kg/K).
Hvap (float): Heat of vaporization (J/kg).
sigma (float): Surface tension (N/m).
dPsat (float): Saturation pressure difference (Pa).
Te (float): Excess wall temperature (K).
Returns:
float: Heat transfer coefficient (W/m^2/K), or an error message if invalid.
"""
try:
return ht_Chen_Edelstein(m=m, x=x, D=D, rhol=rhol, rhog=rhog, mul=mul, mug=mug,
kl=kl, Cpl=Cpl, Hvap=Hvap, sigma=sigma, dPsat=dPsat, Te=Te)
except Exception as e:
return f"Error: {str(e)}"Online Calculator
Mass flow rate (kg/s).
Quality at the tube interval (dimensionless).
Tube diameter (m).
Liquid density (kg/m^3).
Gas density (kg/m^3).
Liquid viscosity (Pa*s).
Gas viscosity (Pa*s).
Liquid thermal conductivity (W/m/K).
Liquid heat capacity (J/kg/K).
Heat of vaporization (J/kg).
Surface tension (N/m).
Saturation pressure difference (Pa).
Excess wall temperature (K).