docs: def Domanski_Didion(x, rhol, rhog, mul, mug):
      r'''Calculates void fraction in two-phase flow according to the model of 
      [1]_ also given in [2]_ and [3]_.
  
      if Xtt < 10:
      
      .. math::
          \alpha = (1 + X_{tt}^{0.8})^{-0.378}
          
      Otherwise:
      
      .. math::
          \alpha = 0.823- 0.157\ln(X_{tt})
      
      Parameters
      ----------
      x : float
          Quality at the specific tube interval []
      rhol : float
          Density of the liquid [kg/m^3]
      rhog : float
          Density of the gas [kg/m^3]
      mul : float
          Viscosity of liquid [Pa*s]
      mug : float
          Viscosity of gas [Pa*s]
  
      Returns
      -------
      alpha : float
          Void fraction (area of gas / total area of channel), [-]
  
      Notes
      -----
      [1]_ has been reviewed. [2]_ gives an exponent of -0.38 instead of -0.378
      as is in [1]_. [3]_ describes only the novel half of the correlation.
      The portion for Xtt > 10 is novel; the other is said to be from their 31st 
      reference, Wallis.
      
      There is a discontinuity at Xtt = 10.
          
      Examples
      --------
      >>> Domanski_Didion(.4, 800, 2.5, 1E-3, 1E-5)
      0.9355795597059169
  
      References
      ----------
      .. [1] Domanski, Piotr, and David A. Didion. "Computer Modeling of the 
         Vapor Compression Cycle with Constant Flow Area Expansion Device." 
         Report. UNT Digital Library, May 1983. 
      .. [2] Xu, Yu, and Xiande Fang. "Correlations of Void Fraction for Two-
         Phase Refrigerant Flow in Pipes." Applied Thermal Engineering 64, no. 
         1-2 (March 2014): 242–51. doi:10.1016/j.applthermaleng.2013.12.032. 
      .. [3] Dalkilic, A. S., S. Laohalertdecha, and S. Wongwises. "Effect of 
         Void Fraction Models on the Two-Phase Friction Factor of R134a during 
         Condensation in Vertical Downward Flow in a Smooth Tube." International 
         Communications in Heat and Mass Transfer 35, no. 8 (October 2008): 
         921-27. doi:10.1016/j.icheatmasstransfer.2008.04.001.
      '''
      Xtt = Lockhart_Martinelli_Xtt(x, rhol, rhog, mul, mug)
      if Xtt < 10:
          return (1 + Xtt**0.8)**-0.378
      else:
          return 0.823 - 0.157*log(Xtt)