def METHOD_NAME(self, system_frame):
    m = system_frame
    results = solve(m)
    assert_optimal_termination(results)
    # check products - treated water
    assert pytest.approx(0.011574, rel=1e-3) == value(
        m.fs.filtered_water.flow_mass_comp[0, "H2O"]
    )
    assert pytest.approx(4.62963e-11, rel=1e-3) == value(
        m.fs.filtered_water.flow_mass_comp[0, "arsenic"]
    )
    assert pytest.approx(6.94444e-11, rel=1e-3) == value(
        m.fs.filtered_water.flow_mass_comp[0, "uranium"]
    )
    assert pytest.approx(2.893519e-8, rel=1e-3) == value(
        m.fs.filtered_water.flow_mass_comp[0, "nitrate"]
    )
    assert pytest.approx(1.157407e-10, rel=1e-3) == value(
        m.fs.filtered_water.flow_mass_comp[0, "phosphates"]
    )
    assert value(m.fs.filtered_water.flow_mass_comp[0, "iron"]) < 1e-12
    assert value(m.fs.filtered_water.flow_mass_comp[0, "filtration_media"]) < 1e-12
    # check products - byproducts
    assert value(m.fs.byproduct.flow_mass_comp[0, "H2O"]) < 1e-12
    assert pytest.approx(4.16667e-10, rel=1e-3) == value(
        m.fs.byproduct.flow_mass_comp[0, "arsenic"]
    )
    assert pytest.approx(6.25e-10, rel=1e-3) == value(
        m.fs.byproduct.flow_mass_comp[0, "uranium"]
    )
    assert pytest.approx(8.680556e-8, rel=1e-3) == value(
        m.fs.byproduct.flow_mass_comp[0, "nitrate"]
    )
    assert pytest.approx(1.041667e-9, rel=1e-3) == value(
        m.fs.byproduct.flow_mass_comp[0, "phosphates"]
    )
    assert value(m.fs.byproduct.flow_mass_comp[0, "iron"]) < 1e-12
    assert pytest.approx(8.680556e-9, rel=1e-3) == value(
        m.fs.byproduct.flow_mass_comp[0, "filtration_media"]
    )