A new gradient projection active set approach for reliable and economic operation of combined energy systems