def METHOD_NAME(self, system_frame):
    m = system_frame
    results = solve(m)
    assert_optimal_termination(results)
    # check products
    assert value(
        m.fs.photothermal_water.properties[0].flow_mass_comp["H2O"]
    ) == pytest.approx(199.992, rel=1e-3)
    assert value(
        m.fs.candop_treated.properties[0].flow_mass_comp["H2O"]
    ) == pytest.approx(22.221, rel=1e-3)
    assert value(
        m.fs.candop_treated.properties[0].flow_mass_comp[
            "bioconcentrated_phosphorous"
        ]
    ) == pytest.approx(0.001000222, rel=1e-3)
    assert value(
        m.fs.candop_treated.properties[0].flow_mass_comp["nitrogen"]
    ) == pytest.approx(0.001944444, rel=1e-3)
    assert value(
        m.fs.candop_treated.properties[0].flow_mass_comp["phosphates"]
    ) == pytest.approx(0.0003333333333, rel=1e-3)
    assert value(
        m.fs.candop_byproduct.properties[0].flow_mass_comp["nitrous_oxide"]
    ) == pytest.approx(0.00583356, rel=1e-3)