The ideas in this paper are motivated by an increased need for systematic data-enabled resource management of large-scale electric energy systems. The basic control objective is to manage uncertain disturbances, power imbalances in particular, by optimizing available power resources. To that end, we start with a centralized optimal control problem formulation of system-level performance objective subject to complex interconnection constraints and constraints representing highly heterogeneous internal dynamics of system components. To manage spatial complexity, an inherent multi-layered structure is utilized by expressing interconnection constraints in terms of unified power variables and their dynamics. Similarly, the internal dynamics of components and sub-systems (modules), including their primary automated feedback control, is modeled so that their input-output characterization is also expressed in terms of power variables. This representation is shown to be key to managing the multi-spatial complexity of the problem. In this unifying energy/power state space, the system constraints are all fundamentally convex, resulting in the convex dynamic optimization problem, for typically utilized quadratic cost functions. Based on this, an interactive multi-layered modeling and control method is introduced. While the approach is fundamentally based on the primal-dual decomposition of the centralized problem, it is proposed for the first time to utilize sensitivity functions of distributed agents for solving the primal distributed problem. Iterative communication complexity typically required for convergence of point-wise information exchange is replaced by