A method for optimizing operational settings for a plurality of energy devices includes representing each of the plurality of energy devices in terms of a set of decision variables and operational parameters. The decision variables and operational parameters are constrained based on operational conditions and interrelationship within the plurality of energy devices. A two-tiered model of the plurality of energy devices is generated wherein a top tier of the model represents interaction of various sub-models and a bottom tier of the model includes a set of the sub-models that form the top tier, each sub-model representing detailed operation of the plurality of energy devices. The two-tiered model is optimized to provide either a schedule of operation for the plurality of energy devices or real-time control for the plurality of energy devices.