def METHOD_NAME(b, t):
    return b.heat_consumption[t] == pyunits.convert(
        b.properties_in[t].flow_mass_comp["nonbiodegradable_cod"]
        * b.heat_required_per_vfa_mass[t],
        to_units=pyunits.kJ / pyunits.s,
    )