def METHOD_NAME(b, t):
    return (
        b.feed_side.properties_out[t].pressure_osm_phase["Liq"]
        - b.permeate_side[t, 1.0].pressure_osm_phase["Liq"]
    ) / b.feed_side.properties_out[t].pressure