def METHOD_NAME(self, inlet_state):
    """
    This is Equation 2.10 in [1]
    """
    time = self.flowsheet().time
    def power_rule(b, t):
        cp_mass = inlet_state[t].cp_mol / inlet_state[t].mw
        power_expr = (
            inlet_state[t].temperature
            * cp_mass
            * inlet_state[t].flow_mass
            * (
                (
                    (
                        pyunits.convert(inlet_state[t].pressure, pyunits.bar)
                        + pyunits.convert(self.boost_pressure[t], pyunits.bar)
                    )
                    / pyunits.convert(inlet_state[t].pressure, pyunits.bar)
                )
                ** b.beta[t]
                - 1.0
            )
        )
        return b.power[t] == pyunits.convert(power_expr, pyunits.kW)
    self.power_eqn = Constraint(time, rule=power_rule)