code: def tyn_calus(target, VA, VB, sigma_A, sigma_B, temperature='pore.temperature',
                viscosity='pore.viscosity'):
      r
      T = target[temperature]
      mu = target[viscosity]
      A = 8.93e-8*(VB*1e6)**0.267/(VA*1e6)**0.433*T
      B = (sigma_B/sigma_A)**0.15/(mu*1e3)
      value = A*B
      return value
docs: r"""
      Uses Tyn_Calus model to estimate diffusion coefficient in a dilute liquid
      solution of A in B from first principles at conditions of interest
  
      Parameters
      ----------
      target : OpenPNM Object
          The object for which these values are being calculated.  This
          controls the length of the calculated array, and also provides
          access to other necessary thermofluid properties.
  
      VA : float, array_like
          Molar volume of component A at boiling temperature (m3/mol)
  
      VB : float, array_like
          Molar volume of component B at boiling temperature (m3/mol)
  
      sigmaA:  float, array_like
          Surface tension of component A at boiling temperature (N/m)
  
      sigmaB:  float, array_like
          Surface tension of component B at boiling temperature (N/m)
  
      pressure : string
          The dictionary key containing the pressure values in Pascals (Pa)
  
      temperature : string
          The dictionary key containing the temperature values in Kelvin (K)