def Nu_Shitsman(Re, Pr_b, Pr_w):
    """Calculates internal convection Nusselt number for turbulent vertical
    upward flow in a pipe under supercritical conditions according to [1]_ and
    [2] as shown in both [3]_ and [4]_.
        
    .. math::
        Nu_b = 0.023 Re_b^{0.8}(min(Pr_b, Pr_w))^{0.8}

    Parameters
    ----------
    Re : float
        Reynolds number with bulk fluid properties, [-]
    Pr_b : float
        Prandtl number with bulk fluid properties, [-]
    Pr_w : float
        Prandtl number with wall fluid properties, [-]

    Returns
    -------
    Nu : float
        Nusselt number with bulk fluid properties, [-]

    Notes
    -----
    [3]_ states this correlation was developed with D = 7.8 and 8.2 mm and with
    a `Pr` approximately 1. [3]_ ranked it third in the enhanced heat transfer
    category, with a MAD as 11.5%
    
    [4]_ cites a [1]_ as the source of the correlation. Neither have been
    reviewed, and both are in Russian. [4]_ lists this as third most accurate
    of the 14 correlations reviewed from a database of all regimes.

    Examples
    --------
    >>> Nu_Shitsman(1E5, 1.2, 1.6)
    266.1171311047253

    References
    ----------
    .. [1] M. E Shitsman, Impairment of the heat transmission at supercritical 
       pressures, High. Temperature, 1963, 1(2): 237-244
    .. [2] Miropol’skiy ZL, Shitsman ME (1957). Heat transfer to water and 
       steam at variable specific heat. J Tech Phys XXVII(10): 2359-2372
    .. [3] Chen, Weiwei, Xiande Fang, Yu Xu, and Xianghui Su. "An Assessment of
       Correlations of Forced Convection Heat Transfer to Water at 
       Supercritical Pressure." Annals of Nuclear Energy 76 (February 2015): 
       451-60. doi:10.1016/j.anucene.2014.10.027.
    .. [4] Yu, Jiyang, Baoshan Jia, Dan Wu, and Daling Wang. "Optimization of 
       Heat Transfer Coefficient Correlation at Supercritical Pressure Using 
       Genetic Algorithms." Heat and Mass Transfer 45, no. 6 (January 8, 2009): 
       757-66. doi:10.1007/s00231-008-0475-4.
    """
    return 0.023 * Re ** 0.8 * min(Pr_b, Pr_w) ** 0.8
