docs: def Serth_HEDH(D, sigma, Hvap, rhol, rhog):
      r'''Calculates critical heat flux for nucleic boiling of a tube bundle
      according to [2]_, citing [3]_, and using [1]_ as the original form.
  
      .. math::
          q_c = KH_{vap} \rho_g^{0.5}\left[\sigma g (\rho_L-\rho_g)\right]^{0.25}
  
          K = 0.123 (R^*)^{-0.25} \text{ for 0.12 < R* < 1.17}
  
          K = 0.118
  
          R^* = \frac{D}{2} \left[\frac{g(\rho_L-\rho_G)}{\sigma}\right]^{0.5}
  
      Parameters
      ----------
      D : float
          Diameter of tubes [m]
      sigma : float
          Surface tension of liquid [N/m]
      Hvap : float
          Heat of vaporization of the fluid at T, [J/kg]
      rhol : float
          Density of the liquid [kg/m^3]
      rhog : float
          Density of the produced gas [kg/m^3]
  
      Returns
      -------
      q: float
          Critical heat flux [W/m^2]
  
      Notes
      -----
      A further source for this would be nice.
  
      Examples
      --------
      >>> Serth_HEDH(D=0.0127, sigma=8.2E-3, Hvap=272E3, rhol=567, rhog=18.09)
      351867.46522901946
  
      References
      ----------
      .. [1] Zuber N. "On the stability of boiling heat transfer". Trans ASME
         1958 80:711-20.
      .. [2] Serth, R. W., Process Heat Transfer: Principles,
         Applications and Rules of Thumb. 2E. Amsterdam: Academic Press, 2014.
      .. [3] Schlünder, Ernst U, and International Center for Heat and Mass
         Transfer. Heat Exchanger Design Handbook. Washington:
         Hemisphere Pub. Corp., 1987.
      '''
      R = D/2*(g*(rhol-rhog)/sigma)**0.5
      if 0.12 <= R  <= 1.17:
          K = 0.125*R**-0.25
      else:
          K = 0.118
      return K*Hvap*rhog**0.5*(g*sigma*(rhol-rhog))**0.25