NU_VCYL_MCADAMS

This function computes natural-convection Nusselt number for a vertical isothermal cylinder using the McAdams-Weiss-Saunders correlation in ht. It accepts optional laminar/turbulent branch forcing.

The formulation is piecewise with Rayleigh-number dependence:

Nu = C\,Ra^n,\qquad Ra = Gr\,Pr

and returns a dimensionless Nusselt number referenced to cylinder height.

Excel Usage

=NU_VCYL_MCADAMS(Pr, Gr, turbulent)

Pr (float, required): Prandtl number of the fluid (dimensionless).
Gr (float, required): Grashof number based on cylinder height (dimensionless).
turbulent (bool, optional, default: null): Whether to force turbulent or laminar regime (dimensionless).

Returns (float): Nusselt number based on cylinder height (dimensionless).

Example 1: McAdams Weiss Saunders example

Inputs:

Pr	Gr
0.7	20000000000

Excel formula:

=NU_VCYL_MCADAMS(0.7, 20000000000)

Expected output:

313.318

Example 2: McAdams Weiss Saunders forced laminar

Inputs:

Pr	Gr	turbulent
0.7	500000000	false

Excel formula:

=NU_VCYL_MCADAMS(0.7, 500000000, FALSE)

Expected output:

80.6992

Example 3: McAdams Weiss Saunders forced turbulent

Inputs:

Pr	Gr	turbulent
0.7	50000000000	true

Excel formula:

=NU_VCYL_MCADAMS(0.7, 50000000000, TRUE)

Expected output:

425.239

Example 4: McAdams Weiss Saunders mid range

Inputs:

Pr	Gr
1	8000000000

Excel formula:

=NU_VCYL_MCADAMS(1, 8000000000)

Expected output:

260

Python Code

Show Code

from ht.conv_free_immersed import Nu_vertical_cylinder_McAdams_Weiss_Saunders as ht_Nu_vertical_cylinder_McAdams_Weiss_Saunders

def Nu_vcyl_McAdams(Pr, Gr, turbulent=None):
    """
    Calculate the Nusselt number for a vertical cylinder using McAdams-Weiss-Saunders.

    See: https://ht.readthedocs.io/en/latest/ht.conv_free_immersed.html

    This example function is provided as-is without any representation of accuracy.

    Args:
        Pr (float): Prandtl number of the fluid (dimensionless).
        Gr (float): Grashof number based on cylinder height (dimensionless).
        turbulent (bool, optional): Whether to force turbulent or laminar regime (dimensionless). Default is None.

    Returns:
        float: Nusselt number based on cylinder height (dimensionless).
    """
    try:
        result = ht_Nu_vertical_cylinder_McAdams_Weiss_Saunders(Pr, Gr, turbulent=turbulent)
        if result is None:
            return "Error: Result is None"
        return float(result)
    except Exception as e:
        return f"Error: {str(e)}"

Online Calculator

Pr *

Prandtl number of the fluid (dimensionless).

Gr *

Grashof number based on cylinder height (dimensionless).

turbulent

Whether to force turbulent or laminar regime (dimensionless).