def METHOD_NAME(self):
    mods = [
        chemicals.volume.Yen_Woods_saturation,
        chemicals.volume.Rackett,
        chemicals.volume.Yamada_Gunn,
        chemicals.volume.Townsend_Hales,
        chemicals.volume.Bhirud_normal,
        chemicals.volume.COSTALD,
        # chemicals.volume.Campbell_Thodos,  # numba version not working
        # chemicals.volume.SNM0,  # numba version not working
        # chemicals.volume.CRC_inorganic,  # requires rho
        # chemicals.volume.COSTALD_compressed,  # requires Psat
    ]
    h2o = op.phase.Species(network=self.net, species='water')
    Vm = []
    for f in mods:
        Vm.append(op.models.phase.chemicals_wrapper(h2o, f=f).mean())
    assert_allclose(Vm, 1.88e-5, rtol=0.2)