MARTIN_SIMS

Estimates two-phase convective heat transfer in tube flow using the Martin-Sims correlation. This approach scales a liquid-only heat transfer coefficient by an empirical factor based on phase velocity effects.

The functional structure is:

h_{TP} = h_l\,\phi

where h_l may be supplied directly or computed from provided liquid properties and tube dimensions.

Excel Usage

=MARTIN_SIMS(m, x, D, rhol, rhog, hl, Cpl, kl, mu_b, mu_w, L)

m (float, required): Mass flow rate (kg/s).
x (float, required): Quality at the tube interval (-).
D (float, required): Tube diameter (m).
rhol (float, required): Liquid density (kg/m^3).
rhog (float, required): Gas density (kg/m^3).
hl (float, optional, default: null): Liquid-only heat transfer coefficient (W/m^2/K).
Cpl (float, optional, default: null): Liquid heat capacity at constant pressure (J/kg/K).
kl (float, optional, default: null): Liquid thermal conductivity (W/m/K).
mu_b (float, optional, default: null): Liquid viscosity at bulk conditions (Pa*s).
mu_w (float, optional, default: null): Liquid viscosity at wall temperature (Pa*s).
L (float, optional, default: null): Tube length (m).

Returns (float): Heat transfer coefficient (W/m^2/K).

Example 1: Martin-Sims with provided hl

Inputs:

m	x	D	rhol	rhog	hl
1	0.9	0.3	1000	2.5	141.2

Excel formula:

=MARTIN_SIMS(1, 0.9, 0.3, 1000, 2.5, 141.2)

Expected output:

5563.28

Example 2: Martin-Sims with computed hl

Inputs:

m	x	D	rhol	rhog	Cpl	kl	mu_b	mu_w	L
1	0.9	0.3	1000	2.5	2300	0.6	0.001	0.0012	24

Excel formula:

=MARTIN_SIMS(1, 0.9, 0.3, 1000, 2.5, 2300, 0.6, 0.001, 0.0012, 24)

Expected output:

5977.51

Example 3: Without wall viscosity correction

Inputs:

m	x	D	rhol	rhog	Cpl	kl	mu_b	L
0.7	0.4	0.05	980	2	4180	0.6	0.001	3

Excel formula:

=MARTIN_SIMS(0.7, 0.4, 0.05, 980, 2, 4180, 0.6, 0.001, 3)

Expected output:

20785.9

Example 4: Higher quality flow

Inputs:

m	x	D	rhol	rhog	Cpl	kl	mu_b	mu_w	L
1.4	0.8	0.06	950	2.8	3600	0.52	0.0014	0.0016	4

Excel formula:

=MARTIN_SIMS(1.4, 0.8, 0.06, 950, 2.8, 3600, 0.52, 0.0014, 0.0016, 4)

Expected output:

36660.4

Python Code

Show Code

from ht.conv_two_phase import Martin_Sims as ht_Martin_Sims

def Martin_Sims(m, x, D, rhol, rhog, hl=None, Cpl=None, kl=None, mu_b=None, mu_w=None, L=None):
    """
    Calculate two-phase heat transfer coefficient using the Martin-Sims correlation.

    See: https://ht.readthedocs.io/en/latest/ht.conv_two_phase.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 (-).
        D (float): Tube diameter (m).
        rhol (float): Liquid density (kg/m^3).
        rhog (float): Gas density (kg/m^3).
        hl (float, optional): Liquid-only heat transfer coefficient (W/m^2/K). Default is None.
        Cpl (float, optional): Liquid heat capacity at constant pressure (J/kg/K). Default is None.
        kl (float, optional): Liquid thermal conductivity (W/m/K). Default is None.
        mu_b (float, optional): Liquid viscosity at bulk conditions (Pa*s). Default is None.
        mu_w (float, optional): Liquid viscosity at wall temperature (Pa*s). Default is None.
        L (float, optional): Tube length (m). Default is None.

    Returns:
        float: Heat transfer coefficient (W/m^2/K).
    """
    try:
        return ht_Martin_Sims(
            m=m,
            x=x,
            D=D,
            rhol=rhol,
            rhog=rhog,
            hl=hl,
            Cpl=Cpl,
            kl=kl,
            mu_b=mu_b,
            mu_w=mu_w,
            L=L,
        )
    except Exception as e:
        return f"Error: {str(e)}"

Online Calculator

m *

Mass flow rate (kg/s).

x *

Quality at the tube interval (-).

D *

Tube diameter (m).

rhol *

Liquid density (kg/m^3).

rhog *

Gas density (kg/m^3).

Liquid-only heat transfer coefficient (W/m^2/K).

Cpl

Liquid heat capacity at constant pressure (J/kg/K).

Liquid thermal conductivity (W/m/K).

mu_b

Liquid viscosity at bulk conditions (Pa*s).

mu_w

Liquid viscosity at wall temperature (Pa*s).

Tube length (m).