def METHOD_NAME(self):
    mods = [
        # The following 3 require constants
        # 'Viswanath_Natarajan_3',
        # 'Viswanath_Natarajan_2',
        # 'Viswanath_Natarajan_2_exponential',
        'Letsou_Stiel',
        # 'Przedziecki_Sridhar',  # Requires molar volume at temperature
        # 'Lucas',  # Requires saturation pressure at temperature
    ]
    pn = op.network.Demo()
    h2o = op.phase.Species(network=pn, species='water')
    mu = []
    for m in mods:
        f = getattr(chemicals.viscosity, m)
        mu.append(op.models.phase.chemicals_wrapper(phase=h2o, f=f).mean())
    assert_allclose(mu, 0.000621, rtol=0.02)