def METHOD_NAME():
    m = construct_dummy_model(
        {
            "sigma": 4.443,
            "epsilon": 215.7,
            "mw": 30.070e-3,
            "omega": 0.099,
        }
    )
    ChapmanEnskogLennardJones.visc_d_phase_comp.build_parameters(m.params, "Vap")
    expr = ChapmanEnskogLennardJones.visc_d_phase_comp.return_expression(
        m.props[1], m.params, "Vap", m.props[1].temperature
    )
    expr_micropoise = pyunits.convert(expr, pyunits.micropoise)
    # Pulled from Table 9.2, Properties of Gases and Liquids
    T_list = [47, 117, 247]  # Gas temperature in C
    visc_list = [100, 120, 156]  # Experimental viscosities in millipoise
    for i in range(3):
        m.props[1].temperature.value = T_list[i] + 273.15
        assert pyo.value(expr_micropoise) == pytest.approx(visc_list[i], rel=0.03)
    assert_units_equivalent(expr, pyunits.Pa * pyunits.s)