Thermodynamic modeling, energy equipartition, and nonconservation of entropy for discrete-time dynamical systems

In this paper we develop thermodynamic models for discrete-time large-scale dynamical systems. Specifically, using compartmental dynamical system theory, we develop energy flow models possessing energy conservation, energy equipartition, temperature equipartition, and entropy nonconservation principles for discrete-time, large-scale dynamical systems. Furthermore, we introduce a new and dual notion to entropy, namely, ectropy, as a measure of the tendency of a dynamical system to do useful work and grow more organized, and show that conservation of energy in an isolated thermodynamic system necessarily leads to nonconservation of ectropy and entropy. In addition, using the system ectropy as a Lyapunov function candidate we show that our discrete-time, large-scale thermodynamic energy flow model has convergent trajectories to Lyapunov stable equilibria determined by the system initial subsystem energies.


I. INTRODUCTION
Thermodynamic principles have been repeatedly used in continuous-time dynamical system theory as well as information theory for developing models that capture the exchange of nonnegative quantities (e.g., mass and energy) between coupled subsystems (see [1] and the numerous references therein). In particular, conservation laws (e.g., mass and energy) are used to capture the exchange of material between coupled macroscopic subsystems known as compartments. Each compartment is assumed to be kinetically homogeneous, that is, any material entering the compartment is instantaneously mixed with the material in the compartment. These models are known as compartmental models and are widespread in engineering systems as well as biological and ecological sciences [2]- [4]. Even though the compartmental models developed in the literature are based on the first law of thermodynamics involving conservation of energy principles, they do not tell us whether any particular process can actually occur; that is, they do not address the second law of thermodynamics involving entropy notions in the energy flow between subsystems.
The goal of the present paper is directed toward developing nonlinear discrete-time compartmental models that are consistent with thermodynamic principles. Specifically, since thermodynamic models are concerned with energy flow among subsystems, we develop a nonlinear compartmental dynamical system model that is characterized by energy conservation laws capturing the exchange of energy between coupled macroscopic subsystems. Furthermore, using graph theoretic notions we state three thermodynamic axioms consistent with the zeroth and second laws of This research was supported in part by AFOSR under Grant F49620-03-1-0178 and NSF under Grant ECS-0133038. thermodynamics that ensure that our large-scale dynamical system model gives rise to a thermodynamically consistent energy flow model. Specifically, using a large-scale dynamical systems theory perspective, we show that our compartmental dynamical system model leads to a precise formulation of the equivalence between work energy and heat in a large-scale dynamical system.
Next, we give a deterministic definition of entropy for a large-scale dynamical system that is consistent with the classical thermodynamic definition of entropy and show that it satisfies a Clausius-type inequality leading to the law of entropy nonconservation. Furthermore, we introduce a new and dual notion to entropy, namely, ectropy, as a measure of the tendency of a large-scale dynamical system to do useful work and grow more organized, and show that conservation of energy in an isolated thermodynamically consistent system necessarily leads to nonconservation of ectropy and entropy. Then, using the system ectropy as a Lyapunov function candidate we show that our thermodynamically consistent large-scale nonlinear dynamical system model possesses a continuum of equilibria and is semistable, that is, it has convergent subsystem energies to Lyapunov stable energy equilibria determined by the large-scale system initial subsystem energies. In addition, we show that the steady-state distribution of the largescale system energies is uniform leading to system energy equipartitioning corresponding to a minimum ectropy and a maximum entropy equilibrium state. Finally, we note that the proofs of the results in this paper are similar to the proofs given in [1] and hence are omitted. For details of the proofs see [5].

II. MATHEMATICAL PRELIMINARIES
In this section we introduce notation, several definitions, and some key results needed for developing the main results of this paper. Let R denote the set of real numbers, Z + denote the set of nonnegative integers, R n denote the set of n×1 column vectors, (·) T denote transpose, and let I n or I denote the n×n identity matrix. For v ∈ R q we write v ≥≥ 0 (respectively, v >> 0) to indicate that every component of v is nonnegative (respectively, positive). In this case we say that v is nonnegative or positive, respectively.
Let R q + and R q + denote the nonnegative and positive orthants of R q , that is, if v ∈ R q , then v ∈ R q + and v ∈ R q + are equivalent, respectively, to v ≥≥ 0 and v >> 0. Finally, we write ∆V (x(k)) for V (x(k + 1)) − V (x(k)).
The following definition introduces the notion of nonnegative functions [6].
if and only if w : V → R q is nonnegative.
The following definition introduces several types of stability for the discrete-time nonnegative dynamical system (1).
Definition 2.2: The equilibrium solution r(k) ≡ r e of (1) is Lyapunov stable if, for every ε > 0, there exists (1) is semistable if it is Lyapunov stable and there exists δ > 0 such that if r 0 ∈ B δ (r e ) ∩ R q + , then lim k→∞ r(k) exists and converges to a Lyapunov stable equilibrium point. The equilibrium solution r(k) ≡ r e of (1) is asymptotically stable if it is Lyapunov stable and there exists δ > 0 such that if r 0 ∈ B δ (r e ) ∩ R q + , then lim k→∞ r(k) = r e . Finally, the equilibrium solution r(k) ≡ r e of (1) is globally asymptotically stable if the previous statement holds for all r 0 ∈ R q + .

III. THERMODYNAMIC MODELING FOR DISCRETE-TIME DYNAMICAL SYSTEMS
The fundamental and unifying concept in the analysis of complex (large-scale) dynamical systems is the concept of energy. The energy of a state of a dynamical system is the measure of its ability to produce changes (motion) in its own system state as well as changes in the system states of its surroundings. These changes occur as a direct consequence of the energy flow between different subsystems within the dynamical system. Since heat (energy) is a fundamental concept of thermodynamics involving the capacity of hot bodies (more energetic subsystems) to produce work, thermodynamics is a theory of large-scale dynamical systems [1]. As in thermodynamic systems, dynamical systems can exhibit energy (due to friction) that becomes unavailable to do useful work. This is in turn contributes to an increase in system entropy; a measure of the tendency of a system to lose the ability to do useful work.
To develop discrete-time compartmental models that are consistent with thermodynamic principles, consider a discrete-time large-scale dynamical system G involving q interconnected subsystems. Let E i : Z + → R + denote the energy (and hence a nonnegative quantity) of the ith subsystem, let S i : Z + → R denote the external energy supplied to (or extracted from) the ith subsystem, let σ ij : R q + → R + , i = j, i, j = 1, . . . , q, denote the exchange of energy from the jth subsystem to the ith subsystem, and let σ ii : R q + → R + , i = 1, . . . , q, denote the energy loss from the ith subsystem. An energy balance equation for the ith subsystem yields or, equivalently, in vector form, where Equation (2) yields a conservation of energy equation and implies that the change of energy stored in the ith subsystem is equal to the external energy supplied to (or extracted from) the ith subsystem plus the energy gained by the ith subsystem from all other subsystems due to subsystem coupling minus the energy dissipated from the ith subsystem. Note that (3) or, equivalently, (2) is a statement reminiscent of the first law of thermodynamics for each of the subsystems, with E i (·), S i (·), σ ij (·), i = j, and σ ii (·), i = 1, . . . , q, playing the role of the ith subsystem internal energy, energy supplied to (or extracted from) the ith subsystem, the energy exchange between subsystems due to coupling, and the energy dissipated to the environment, respectively.
To further elucidate that (3) is essentially the statement of the principle of the conservation of energy let the total energy in the discrete-time large-scale dynamical system G be given by U , and let the energy received by the discrete-time largescale dynamical system G (in forms other than work) over the discrete-time interval . Then, premultiplying (3) by e T and using the fact that e T w(E) ≡ e T E, it follows that where ∆U U (k 2 ) − U (k 1 ) denotes the variation in the total energy of the discrete-time large-scale dynamical system G over the discrete-time interval {k 1 , . . . , k 2 }. This is a statement of the first law of thermodynamics for the discrete-time large-scale dynamical system G and gives a precise formulation of the equivalence between variation in system internal energy and heat.
It is important to note that our discrete-time large-scale dynamical system model does not consider work done by the system on the environment nor work done by the environment on the system. Hence, Q can be interpreted physically as the amount of energy that is received by the system in forms other than work. The extension of addressing work performed by and on the system can be easily handeled by including an additional state equation, coupled to the energy balance equation (3), involving volume states for each subsystem [1]. Since this slight extension does not alter any of the results of the paper, it is not considered here for simplicity of exposition.
For our large-scale dynamical system model G, we assume that σ ij (E) = 0, E ∈ R q + , whenever E j = 0, i, j = 1, . . . , q. This constraint implies that if the energy of the jth subsystem of G is zero, then this subsystem cannot supply any energy to its surroundings nor dissipate energy to the environment. Furthermore, for the remainder of this paper we assume that This constraint implies that the energy that can be dissipated, extracted, or exchanged by the ith subsystem cannot exceed the current energy in the subsystem. Note that this assumption implies that E(k) ≥≥ 0 for all k ≥ k 0 .
Next, premultiplying (3) by e T and using the fact that e T w(E) ≡ e T E, it follows that Now, for the discrete-time large-scale dynamical system G define the input u(k) S(k) and the output y(k) d(E(k)). Hence, it follows from (6) that the discrete-time large-scale dynamical system G is lossless [7] with respect to the energy supply rate r(u, y) = e T u − e T y and with the energy storage function U (E) e T E, E ∈ R q + . This implies that (see [7] for details) and is the maximum amount of stored energy which can be extracted from the discrete-time large-scale dynamical system G at any discrete-time instant K, and U r (E 0 ) is the minimum amount of energy which can be delivered to the discrete-time largescale dynamical system G to transfer it from a state of minimum potential E(−K) = 0 to a given state E(k 0 ) = E 0 , it follows from (7) that the discrete-time large-scale dynamical system G can deliver to its surroundings all of its stored subsystem energies and can store all of the work done to all of its subsystems. In the case where S(k) ≡ 0, it follows from (6) and the fact that σ ii (E) ≥ 0, E ∈ R q + , i = 1, . . . , q, that the zero solution E(k) ≡ 0 of the discrete-time large-scale dynamical system G with the energy balance equation (3) is Lyapunov stable with Lyapunov function U (E) corresponding to the total energy in the system.
The next result shows that the large-scale dynamical system G is locally controllable.
Proposition 3.1: Consider the discrete-time large-scale dynamical system G with energy balance equation (3). Then for every equilibrium state E e ∈ R q + and every ε > 0 and It follows from Proposition 3.1 that the discrete-time large-scale dynamical system G with the energy balance equation (3) is reachable from and controllable to the origin in R q + . Recall that the discrete-time large-scale dynamical system G with the energy balance equation (3) is reachable from the origin in R q there exists a finite time k f ≥ k 0 and an input S(k) defined on {k 0 , . . . , k f } such that the state E(k), k ≥ k 0 , can be driven from E(k 0 ) = E 0 to E(k f ) = 0. We let U r denote the set of all admissible bounded energy inputs to the discrete-time large-scale dynamical system G such that for any K ≥ −k 0 , the system energy state can be driven from E(−K) = 0 to E(k 0 ) = E 0 ∈ R q + by S(·) ∈ U r , and we let U c denote the set of all admissible bounded energy inputs to the discrete-time large-scale dynamical system G such that for any K ≥ k 0 , the system energy state can be driven from E(k 0 ) = E 0 ∈ R q + to E(K) = 0 by S(·) ∈ U c . Furthermore, let U be an input space that is a subset of bounded continuous R q -valued functions on Z. The spaces U r , U c , and U are assumed to be closed under the shift operator, that is, if S(·) ∈ U (respectively, U c or U r ), then the function S K defined by S K (k) = S(k+K) is contained in U (respectively, U c or U r ) for all K ≥ 0.
The nonlinear energy balance equation (3) can exhibit a full range of nonlinear behavior including bifurcations, limit cycles, and even chaos. However, a thermodynamically consistent energy flow model should ensure that the evolution of the system energy is diffusive (parabolic) in character with convergent subsystem energies. Hence, to ensure a thermodynamically consistent energy flow model we require the following axioms. For the statement of these axioms we first recall the following graph theoretic notions.

Definition 3.1 ( [1]):
A directed graph G(C) associated with the connectivity matrix C ∈ R q×q has vertices {1, 2, . . . , q} and an arc from vertex i to vertex j, i = j, if and only if C (j,i) = 0. A graph G(C) associated with the connectivity matrix C ∈ R q×q is a directed graph for which the arc set is symmetric; that is, C = C T . We say that G(C) is strongly connected if for any ordered pair of vertices (i, j), i = j, there exists a path (i.e., sequence of arcs) leading from i to j.
Recall that C ∈ R q×q is irreducible, that is, there does not exist a permutation matrix such that C is cogredient to a lower-block triangular matrix, if and only if G(C) is strongly connected (see Theorem 2.7 of [8]). Let φ ij (E) σ ij (E) − σ ji (E), E ∈ R q + , denote the net energy exchange between subsystems G i and G j of the discrete-time largescale dynamical system G.
Axiom i): For the connectivity matrix C ∈ R q×q associated with the large-scale dynamical system G defined by and The fact that φ ij (E) = 0 if and only if E i = E j , i = j, implies that subsystems G i and G j of G are connected; alternatively, φ ij (E) ≡ 0 implies that G i and G j are disconnected. Axiom i) implies that if the energies in the connected subsystems G i and G j are equal, then energy exchange between these subsystems is not possible. This is a statement consistent with the zeroth law of thermodynamics which postulates that temperature equality is a necessary and sufficient condition for thermal equilibrium. Furthermore, it follows from the fact that C = C T and rank C = q− 1 that the connectivity matrix C is irreducible which implies that for any pair of subsystems G i and G j , i = j, of G there exists a sequence of connected subsystems of G that connect G i and G j . Axiom ii) implies that energy is exchanged from more energetic subsystems to less energetic subsystems and is consistent with the second law of thermodynamics which states that heat (energy) must flow in the direction of lower temperatures. Furthermore, note that φ ij (E) = −φ ji (E), E ∈ R q + , i = j, i, j = 1, . . . , q, which implies conservation of energy between lossless subsystems. With S(k) ≡ 0, Axioms i) and ii) along with the fact that φ ij (E) = −φ ji (E), E ∈ R q + , i = j, i, j = 1, . . . , q, imply that at a given instant of time energy can only be transported, stored, or dissipated but not created and the maximum amount of energy that can be transported and/or dissipated from a subsystem cannot exceed the energy in the subsystem. Finally, Axiom iii) implies that for any pair of connected subsystems G i and G j , i = j, the energy difference between consecutive time instants is monotonic, that is, Next, we establish a Clausius-type inequality for our thermodynamically consistent energy flow model. Proposition 3.2: Consider the discrete-time large-scale dynamical system G with energy balance equation (3) and assume that Axioms i), ii), and iii) hold. Then for all E 0 ∈ R q + , k f ≥ k 0 , and S(·) ∈ U such that E(k f ) = E(k 0 ) = E 0 , where c > 0, Q i (k) S i (k) − σ ii (E(k)), i = 1, . . . , q, is the amount of net energy (heat) received by the ith subsystem at the kth instant, and E(k), k ≥ k 0 , is the solution to (3) with initial condition E(k 0 ) = E 0 . Furthermore, equality holds in (12) if and only if ∆E i (k) = 0, i = 1, . . . , q, and E i (k) = E j (k), i, j = 1, . . . , q, i = j, k ∈ {k 0 , . . . , k f − 1}.
Inequality (12) is analogous to Clausius' equality and inequality for reversible and irreversible thermodynamics as applied to discrete-time large-scale dynamical systems. It follows from Axiom i) and (3) that for the isolated discretetime large-scale dynamical system G; that is, S(k) ≡ 0 and d(E(k)) ≡ 0, the energy states given by E e = αe, α ≥ 0, correspond to the equilibrium energy states of G. Thus, we can define an equilibrium process as a process where the trajectory of the discrete-time large-scale dynamical system G stays at the equilibrium point of the isolated system G. The input that can generate such a trajectory can be given by S(k) = d(E(k)), k ≥ k 0 . Alternatively, a nonequilibrium process is a process that is not an equilibrium one. Hence, it follows from Axiom i) that for an equilibrium process φ ij (E(k)) ≡ 0, k ≥ k 0 , i = j, i, j = 1, . . . , q, and thus, by Proposition 3.2 and ∆E i = 0, i = 1, . . . , q, inequality (12) is satisfied as an equality. Alternatively, for a nonequilibrium process it follows from Axioms i) − iii) that (12) is satisfied as a strict inequality.
Next, we give a deterministic definition of entropy for the discrete-time large-scale dynamical system G that is consistent with the classical thermodynamic definition of entropy.
Definition 3.2: For the discrete-time large-scale dynamical system G with energy balance equation (3), a function S : for any k 2 ≥ k 1 ≥ k 0 and S(·) ∈ U, is called the entropy of G.
Next, we show that (12) guarantees the existence of an entropy function for G. For this result define the available entropy of the large-scale dynamical system G by where E(k 0 ) = E 0 ∈ R q + and E(K) = 0, and define the required entropy supply of the large-scale dynamical system G by where E(−K) = 0 and E(k 0 ) = E 0 ∈ R q + . Note that the available entropy S a (E 0 ) is the minimum amount of scaled heat (entropy) that can be extracted from the large-scale dynamical system G in order to transfer it from an initial state E(k 0 ) = E 0 to E(K) = 0. Alternatively, the required entropy supply S r (E 0 ) is the maximum amount of scaled heat (entropy) that can be delivered to G to transfer it from the origin to a given initial state E(k 0 ) = E 0 . Theorem 3.1: Consider the discrete-time large-scale dynamical system G with energy balance equation (3) and assume that Axioms ii) and iii) hold. Then there exists an entropy function for G. Moreover, S a (E), E ∈ R q + , and S r (E), E ∈ R q + , are possible entropy functions for G with S a (0) = S r (0) = 0. Finally, all entropy functions (16) Remark 3.1: It is important to note that inequality (12) is equivalent to the existence of an entropy function for G. Sufficiency is simply a statement of Theorem 3.1 while necessity follows from (13) with E(k 2 ) = E(k 1 ). For nonequilibrium process with energy balance equation (3), Definition 3.2 does not provide enough information to define the entropy uniquely. This difficulty has long been pointed out in [9] for thermodynamic systems. A similar remark holds for the definition of ectropy introduced below.
The next proposition gives a closed-form expression for the entropy of G.

Proposition 3.3:
Consider the discrete-time large-scale dynamical system G with energy balance equation (3) and assume that Axioms ii) and iii) hold. Then the function S : R q + → R given by where log e (ce + E) denotes the vector natural logarithm given by [log e (c + E 1 ), . . . , log e (c + E q )] T and c > 0, is an entropy function of G. Remark 3.2: Note that the entropy function given by (17) satisfies (13) as an equality for an equilibrium process and as a strict inequality for a nonequilibrium process.
The entropy expression given by (17) is identical in form to the Boltzmann entropy for statistical thermodynamics. Due to the fact that the entropy is indeterminate to the extent of an additive constant, we can place the constant q log e c to zero by taking c = 1. Since S(E) given by (17) achieves a maximum when all the subsystem energies E i , i = 1, . . . , q, are equal, entropy can be thought of as a measure of the tendency of a system to lose the ability to do useful work, lose order, and settle to a more homogenous state.
Next, we introduce a new and dual notion to entropy, namely, ectropy, describing the status quo of the discretetime large-scale dynamical system G. First, however, we present a dual inequality to inequality (12) that holds for our thermodynamically consistent energy flow model. Proposition 3.4: Consider the discrete-time large-scale dynamical system G with energy balance equation (3) and assume that Axioms i), ii), and iii) hold. Then for all E 0 ∈ R q + , k f ≥ k 0 , and S(·) ∈ U such that E(k f ) = E(k 0 ) = E 0 , where E(k), k ≥ k 0 , is the solution to (3) with initial condition E(k 0 ) = E 0 . Furthermore, equality holds in (18) if and only if ∆E i = 0 and E i = E j , i, j = 1, . . . , q, i = j. Note that inequality (18) is satisfied as an equality for an equilibrium process and as a strict inequality for a nonequilibrium process. Next, we present the definition of ectropy for the discrete-time large-scale dynamical system G.
for any k 2 ≥ k 1 ≥ k 0 and S(·) ∈ U, is called the ectropy of G.
For the next result define the available ectropy of the large-scale dynamical system G by where E(k 0 ) = E 0 ∈ R q + and E(K) = 0, and the required ectropy supply of the large-scale dynamical system G by where E(−K) = 0 and E(k 0 ) = E 0 ∈ R q + . Note that the available ectropy E a (E 0 ) is the maximum amount of scaled heat (ectropy) that can be extracted from the large-scale dynamical system G in order to transfer it from an initial state E(k 0 ) = E 0 to E(K) = 0. Alternatively, the required ectropy supply E r (E 0 ) is the minimum amount of scaled heat (ectropy) that can be delivered to G to transfer it from an initial state E(−K) = 0 to a given state E(k 0 ) = E 0 .
Theorem 3.2: Consider the discrete-time large-scale dynamical system G with energy balance equation (3) and assume that Axioms ii) and iii) hold. Then there exists an ectropy function for G. Moreover, E a (E), E ∈ R q + , and E r (E), E ∈ R q + , are possible ectropy functions for G with E a (0) = E r (0) = 0. Finally, all ectropy functions (22) The next proposition gives a closed-form expression for the ectropy of G.
Proposition 3.5: Consider the discrete-time large-scale dynamical system G with energy balance equation (3) and assume that Axioms ii) and iii) hold. Then the function E : R q + → R given by is an ectropy function of G. Remark 3.3: Note that the ectropy function given by (23) satisfies (19) as an equality for an equilibrium process and as a strict inequality for a nonequilibrium process.
It follows from (23) that ectropy is a measure of the extent to which the system energy deviates from a homogeneous state. Thus, ectropy is the dual of entropy and is a measure of the tendency of the discrete-time large-scale dynamical system G to do useful work and grow more organized.
Inequality (13) is analogous to Clausius' inequality for equilibrium and nonequilibrium thermodynamics as applied to discrete-time large-scale dynamical systems, while inequality (19) is an anti-Clausius' inequality. Moreover, for the ectropy function defined by (23), a thermodynamically consistent discrete-time large-scale dynamical system is dissipative [7] with with respect to the supply rate E T S and with storage function corresponding to the system ectropy E(E). For the entropy function given by (17) note that S(0) = 0, or, equivalently, lim E→0 S(E) = 0, which is consistent with the third law of thermodynamics (Nernst's theorem) which states that the entropy of every system at absolute zero can always be taken to be equal to zero.
(24) Inequality (24) implies that, for any dynamical change in an isolated (i.e., S(k) ≡ 0 and d(E(k)) ≡ 0) discrete-time large-scale system, the entropy of the final state can never be less than the entropy of the initial state. It is important to stress that this result holds for an isolated dynamical system. It is however possible with energy supplied from an external dynamical system (e.g., a controller) to reduce the entropy of the discrete-time large-scale dynamical system. The entropy of both systems taken together however cannot decrease. The above observations imply that when an isolated discrete-time large-scale dynamical system with thermodynamically consistent energy flow characteristics (i.e., Axioms i)−iii) hold) is at a state of maximum entropy consistent with its energy, it cannot be subject to any further dynamical change since any such change would result in a decrease of entropy. This of course implies that the state of maximum entropy is the stable state of an isolated system and this state has to be semistable. Analogously, it follows from (19) that for an isolated discrete-time large-scale dynamical system G the fundamental inequality is satisfied, which implies that the ectropy of the final state of G is always less than or equal to the ectropy of the initial state of G. Hence, for the isolated large-scale dynamical system G the entropy increases if and only if the ectropy decreases. Thus, the state of minimum ectropy is the stable state of an isolated system and this equilibrium state has to be semistable. The next theorem concretizes the above observations. Theorem 3.3: Consider the discrete-time large-scale dynamical system G with energy balance equation (3) with S(k) ≡ 0 and d(E) ≡ 0, and assume that Axioms i) − iii) hold. Then for every α ≥ 0, αe is a Lyapunov equilibrium state of (3). Furthermore, E(k) → 1 q ee T E(k 0 ) as k → ∞ and 1 q ee T E(k 0 ) is a semistable equilibrium state. Finally, if for some m ∈ {1, . . . , q}, σ mm (E) ≥ 0, E ∈ R q + , and σ mm (E) = 0 if and only if E m = 0 1 , then the zero solution E(k) ≡ 0 to (3) is a globally asymptotically stable equilibrium state of (3). Theorem 3.3 implies that the steady-state value of the energy in each subsystem G i of the isolated large-scale dynamical system G is equal, that is, the steady-state energy of the isolated discrete-time large-scale dynamical system G given by E ∞ = 1 q ee T E(k 0 ) = 1 q q i=1 E i (k 0 ) e is uniformly distributed over all subsystems of G. This phenomenon is known as equipartition of energy [1] and is an emergent behavior in thermodynamic systems. The next proposition shows that among all possible energy distributions in the discrete-time large-scale dynamical system G, energy equipartition corresponds to the minimum value of the system's ectropy and the maximum value of the system's entropy.
Proposition 3.6: Consider the discrete-time large-scale dynamical system G with energy balance equation (3), let 1 The assumption σmm(E) ≥ 0, E ∈ R Furthermore, E min E(E * ) = 1 2 β 2 q and S max S(E * ) = q log e (c + β q ) − q log e c. It follows from (24), (25), and Proposition 3.6 that conservation of energy necessarily implies nonconservation of ectropy and entropy. Hence, in an isolated discretetime large-scale dynamical system G all the energy, though always conserved, will eventually be degraded (diluted) to the point where it cannot produce any useful work. Hence, all motion would cease and the dynamical system would be fated to a state of eternal rest (semistability) wherein all subsystems will posses identical energies (energy equipartition). Ectropy would be a minimum and entropy would be a maximum giving rise to a state of absolute disorder. This is precisely what is known in theoretical physics as the heat death of the universe [1].

IV. CONCLUSION
Motivated by energy flow modeling of large-scale interconnected systems, in this paper we develop discrete-time nonlinear compartmental models that are consistent with thermodynamic principles. Specifically, using a discretetime large-scale systems perspective, we develop some of the key properties of thermodynamic systems involving conservation of energy and nonconservation of entropy and ectropy using dynamical systems theory. The concept of entropy for a large-scale dynamical system is defined and shown to be consistent with the classical thermodynamic definition of entropy.