Developing a MILP model to study heat pump flexibility in power systems