def METHOD_NAME(self, model):
    assert model.props[1].flow_mol_comp["D"].value == pytest.approx(
        7.0849, rel=1e-4
    )
    assert model.props[1].mole_frac_phase_comp["Liq", "A"].value == pytest.approx(
        0.09916, rel=1e-4
    )
    assert model.props[1].mole_frac_phase_comp["Liq", "B"].value == pytest.approx(
        0.09916, rel=1e-4
    )
    assert model.props[1].mole_frac_phase_comp["Liq", "C"].value == pytest.approx(
        0.09916, rel=1e-4
    )
    assert model.props[1].mole_frac_phase_comp["Liq", "D"].value == pytest.approx(
        0.70253, rel=1e-4
    )
    assert value(model.rxns[1].k_rxn["R1"]) == pytest.approx(0.72121, rel=1e-4)
    assert value(model.rxns[1].reaction_rate["R1"]) == pytest.approx(
        value(
            model.rxns[1].k_rxn["R1"]
            * model.props[1].mole_frac_phase_comp["Liq", "A"] ** 1
            * model.props[1].mole_frac_phase_comp["Liq", "B"] ** 1
        ),
        rel=1e-6,
    )
    assert value(model.rxns[1].k_eq["R2"]) == pytest.approx(71.4505, rel=1e-4)
    assert value(model.rxns[1].k_eq["R2"]) == pytest.approx(
        value(
            model.props[1].mole_frac_phase_comp["Liq", "D"] ** 1
            * model.props[1].mole_frac_phase_comp["Liq", "B"] ** -1
            * model.props[1].mole_frac_phase_comp["Liq", "C"] ** -1
        ),
        rel=1e-6,
    )