Authors: Anatoly M. Tsirlin (ivsukin@gmail.com), Alexander I. Balunov (ivsukin@gmail.com; 2Department of Cybernetics, Yaroslavl State Technical University, 150023 Yaroslavl, Russia), Ivan A. Sukin (ivsukin@gmail.com)
Categories: Article, irreversible thermodynamics, optimization, problem statements, thermodynamic balances, Carathéodory theorem, averaging, entropy production
Source: Entropy
Doi: 10.3390/e27060649
Authors: Anatoly M. Tsirlin, Alexander I. Balunov, Ivan A. Sukin
In this manuscript, the typical problems of “finite-time thermodynamics”, their general methodology, and the general features of their solutions are considered. We also consider the role of minimal dissipation processes, the properties of the irreversibility index, and the consequences of its existence. A generalization of the Carathéodory theorem for averaged optimization problems corresponding to cyclic processes and the properties of optimal solutions following from it are given. The existence of the irreversibility index for economic macrosystems and their analogies to and differences from thermodynamic systems are proven.
The development of thermodynamics, starting with the work of S. Carnot, is closely connected with extremal problems on the ultimate possibilities of thermodynamic systems. Carnot would hardly have obtained a solution to the problem of the ultimate efficiency of a heat engine if he had posed it in a mathematically rigorous manner, such as the find a law T(t) of change in the temperature of the working fluid of a heat engine receiving heat from a source with temperature T+ and giving the heat to a source with temperature T− such that the ratio ρ of the work produced to the heat taken from the hot source is maximal. After all, one of the variables sought for in this problem is the duration of the cycle τ. The set of admissible values of this variable is limited only by the condition of non-negativity and therefore is not closed. Consequently, according to the Weierstrass theorem, which Carnot could not have known, the problem may not have a solution. Indeed, the maximum of ρ does not exist. The upper limit of the efficiency supremum is achieved in the limit as τ tends to infinity.
This feature was discovered to be characteristic of other extreme problems in thermodynamics (on the minimum work for separation of mixtures, on the efficiency of cycles for systems with sources of finite capacity, etc.). The solutions to these problems are reversible processes in which the exchange flows are arbitrarily close to zero; therefore, an arbitrarily long time is required to exchange a finite amount of substance or energy.
Note that a reversible Carnot heat engine can also have finite power but only in cases where the heat exchange coefficients between the heat sources and the working fluid of the heat engine (the dimensions of the engine) are arbitrarily large. In this case, the Carnot efficiency is the ratio of the “reversible power” to the heat flow taken from the hot source.
Apparently, the first problem in optimization thermodynamics was the “maximum power problem”, about a form of the heat engine cycle receiving heat from a source of infinite capacity with temperature T+ and giving off heat to a source with temperature T−, for which the power of the heat engine would be maximum. This was considered and solved soon after Carnot’s work; however, this result was not noticed. Much later, the problem of the maximum power cycle was solved in [1,2]. At that time, nuclear power began to develop. The cost of building nuclear power plants is very high, and the cost of fuel is much less than for thermal power plants. Under these conditions, obtaining maximum power is much more important than achieving maximum efficiency.
In the problem of maximum power, instead of the quantities of heat and work, heat and work flows appear within the kinetics of heat transfer when the dimensions of the heat engine are arbitrarily large. The authors of the above-mentioned and numerous other studies often solved this problem independently of each other but according to the same Heat exchange flows in contact with each of the sources were assumed to be proportional to the difference in temperature between the source and the working fluid (Newtonian kinetics).The desired cycle was assumed a priori to consist of two isotherms and two adiabats, similar to the Carnot cycle, and the temperatures of the working fluid in contact with the sources were selected based on the condition of maximum power, taking into account the law of conservation of energy and the fact that the entropy flow coming from a hot source must be equal to the entropy flow given to a cold source. The latter meant that the processes inside the working fluid were assumed to be reversible.
It was revealed that the maximum power of a heat engine in contact with reservoirs having temperatures T+ and T− is equal to(1)Nmax=α1α2α1+α2T+−T−2. where α1 and α2 are heat exchange coefficients in contact with sources. The efficiency of the maximum power cycle does not depend on the heat exchange coefficients and is equal to(2)η=1−T−T+.
These works did not answer the following For what heat transfer kinetics does the maximum power cycle consist of two isotherms and two adiabats?Is the efficiency of the maximum power cycle always independent of the heat transfer coefficients?What is the shape of the heat engine cycle that has the maximum efficiency at a fixed power?
Answers to these questions were given in the works [3,4], in which the methods of averaged optimization developed in [5] were used. It was found that the following statements are For any heat exchange kinetics satisfying the natural condition that the direction of heat flow coincides with the sign of the temperature difference of the contacting bodies, the maximum power cycle consists of two isotherms and two adiabats.The efficiency corresponding to this cycle for kinetics that differ from Newtonian depends on the kinetic coefficients.The cycle of a heat engine with maximum efficiency at a fixed power for arbitrary kinetics can consist of three isotherms and three adiabats. There, a condition was obtained under which the number of isotherms is equal to two.
The development of optimization thermodynamics from the ≪maximum power problem≫ to its current state began after the publication of the work of Curzon and Ahlborn [2]. This happened largely because the talented researchers P. Salamon, B. Andresen and K. H. Hoffmann ended up at the Department of Chemistry at the University of Chicago, headed by R.S. Berry. They realized that in real thermodynamic processes, thermal, mass transfer, chemical, etc., the factor of limited duration plays a crucial role. It is important to evaluate the capabilities of systems with an additional restriction on the duration of processes, i.e., in the class of irreversible processes. Thus, a new field of thermodynamics arose, which they called ≪Finite-time thermodynamics≫ (see [6,7,8,9] and many other publications).
The name ≪Finite-Time Thermodynamics≫ (FTT) has taken hold and is hard to abandon. However, this name cannot be considered successful. The overwhelming majority of problems on optimal irreversible processes do not contain time at all. The processes occur in stationary or cyclically operating systems. In the latter case, the indicators are averaged over the cycle time and time is not included in the problem after averaging.
Three classes of problems are Optimal periodic processes of limited durationOptimal cyclic processesOptimal processes in open stationary systems.
At the same time, for any class of problems, including problems on periodic processes, their conditions lead to non-zero flows due to their limited duration.
The correct name for this area of thermodynamics would be ≪Thermodynamics with non-zero flows≫. It emphasizes the features of all problems arising in it.
In order to be able to use the apparatus of equilibrium thermodynamics when solving problems that arise, the system is divided into several equilibrium subsystems, and irreversibility is determined only by the processes that occur when equilibrium subsystems come into contact with each other.
Here, thermodynamic systems should be understood as all systems consisting of a large number of interacting elements, each of which cannot be controlled. Control of such systems is possible only at the macro level and this control imposes the same conditions on all microelements. Such systems are called macrosystems. Examples of macrosystems may be economic systems, where microelements are sellers and buyers, individual ≪households.≫; social systems, where microelements are taxpayers, transport passengers, and others. An important class of macrosystems are segregated systems [10], whose elements (households, units in fluidized bed apparatuses, biomass in biochemical reactors, etc.) evolve under the influence of the environment. The parameters of the environment are determined by external influences and the average effect of the units.
All macrosystems are characterized by the division of the variables characterizing them into intensive and extensive. That implies the existence of processes that proceed spontaneously in one direction, without the expenditure of external resources, and in the other requiring the involvement of external resources for their implementation, i.e., the processes in such systems are irreversible. The irreversibility index for any processes in closed macrosystems increases. However, the form and the fact of the existence of such an index require proof for each type of macrosystem. For thermodynamic systems, such an indicator is entropy, for economic systems, the ≪utility function≫, the proof of the existence of which is given below. The irreversibility indicator is always the product of the flow and the driving force causing this flow. The sign of the flow coincides with the sign of the driving force, so their product is non-negative.
In macrosystems, there is a special subsystem, an intermediary. The intermediary (a reseller in economics, a working fluid in thermodynamics) reduces the growth rate of the irreversibility indicator and, due to this, extracts the target resource from the system (capital in economics, work in thermodynamics). For the intermediary to operate, a difference in the values of intensive variables in the subsystems is necessary. In economics, a difference in prices, in thermodynamics, a difference in temperatures or chemical potentials.
In the following sections, we will present the formulations of typical FTT problems, then the general methodology, and the results of their solution. Let us consider the mathematical features of solving problems of resource exchange and optimization of cyclic processes.
The present work is a review of some results obtained by the authors in solving FTT problems and in developing mathematical methods necessary for such a solution. Our goal is to show by examples that FTT has developed as a complete section of thermodynamics with its own methodology and class of problems, that it has allowed formulating and solving problems that could not even be posed within the framework of classical thermodynamics, that the FTT methodology is also applicable to other types of macrosystems, such as microeconomics, that FTT required the development of averaged optimization methods, and that there are unsolved problems in FTT.
Let us list these Processes of minimum dissipation.Construction of attainability sets of thermodynamic systems.Extraction of maximum work from a non-uniform system in a limited time.Limit capabilities of heat exchange systems.Microeconomic analogies.Optimality conditions for averaged problems and Lyapunov-type equations.
Each non-isolated thermodynamic system exchanges energy and material flows with its environment. In steady-state conditions, these flows can be stationary or periodic. In the latter case, the flow intensity is understood as its average intensity over a period. The mechanism of system operation (kinetics of heat and mass transfer, chemical reactions, etc.) establishes a connection between the incoming and outgoing flows. From the outgoing flows, a target flow can be identified or formed. Its intensity is the system’s productivity. The incoming flows form a cost flow. For a heat engine, the target flow is power. The cost flow is heat taken from a hot source.
The most studied types of macrosystems are thermodynamic. In the first section, we will present the formulations, methodology, and results of solving typical problems for such systems. In the second section we will focus on systems of an economic nature, and in the third we will prove the conditions that determine the type of optimal solution to problems with averaging, characteristic of macrosystems.
The problem of the maximum possible productivity of a thermodynamic system of arbitrary nature under certain constraints and the value of the cost flow that ensures this maximum productivity is a direct generalization of the maximum power problem.
This raises the following For which systems is productivity limited from above, and for which systems, by increasing the cost flow, can productivity be made arbitrarily large, and therefore the problem of maximum productivity has no solution?There is a system of two or more thermodynamic reservoirs and a working fluid that contacts each of them in a steady state or alternately and produces a target flow. How should the contacts of the working fluid be organized to obtain the maximum value of the target flow? What should be considered the efficiency of such a system?What will change if in problem 2, instead of reservoirs, there are sources of finite capacity? In particular, what is the maximum work that can be extracted in a closed thermodynamic system with a fixed duration? This problem coincides with the problem of calculating the exergy of a system in the particular case where the duration of the process is unlimited.How to organize thermodynamic processes so that for a given average flow intensity increase in entropy of the system is minimal (minimum dissipation processes)?In particular, what criterion should be used to evaluate the heat exchange process? How to organize the heat exchange process of two vector flows so that for a given heat load and total heat exchange coefficient the entropy production is minimal?How to construct the region of realizable modes of a thermodynamic system with limited kinetic coefficients in a space along the axes of which the flow intensities are plotted.
The above formulations, of course, do not exhaust the problems in FTT, but they allow us to judge the nature and applied focus of the problems that arise. Note that in most of them the duration of the process does not figure.
Thermodynamic systems are characterized by two types of intensive and extensive. The first type do not change when combining subsystems if before combining they were equal in each subsystem (temperature, pressure, concentration …), while the second type are are added when combining such subsystems (volume, number of moles, internal energy …).
In this case, subsystems can be either passive, in which case intensive variables are determined by extensive variables and the equation of state, or active (the working fluid of a heat engine, or an absorption-desorption cycle), in which case intensive variables (they are controls) are selected to achieve one or another goal.
Let us first characterize the general scheme by which these and similar problems can be solved. The first step in studying the ultimate possibilities of thermodynamic systems is to compile balance relations for matter, energy, and entropy. The last of these relations will include a term characterizing the irreversibility of entropy production rate σ. It is equal to the rate of growth of the entropy of the system. This term is equal to zero if all processes in the system are reversible, and greater than zero for irreversible processes. The non-negativity of dissipation, due to the equations of thermodynamic balances, determines a certain set of realizability in the parameter space of input and output flows.If additional conditions of finite duration of processes, given average flow intensity, and limited kinetic coefficients are imposed on the system, then the dissipation value that is minimally possible under these restrictions is found. In any real system σ≥σmin, which narrows the set of realizability. Now this set takes into account the kinetics of the processes as well as the dimensions of the installation through the heat and mass transfer coefficients.The third step is to obtain from the balance equations the relationship between one or another indicator of the system’s efficiency and the dissipation σ. As a rule, natural efficiency indicators monotonically deteriorate with the growth of σ and reach their limit values in a reversible process, which leads to estimates similar to the Carnot efficiency for processes of very different nature.Since in a complex system dissipation additively depends on dissipation in each of the elementary processes, an important stage of the study is to identify the conditions of minimal dissipation. Optimal organization of processes in a complex system comes down to coordinating individual processes of minimal dissipation with each other.
Thermodynamic balances are a system of equations of material, energy, and entropy balances [11,12]. For simplicity, we will consider them for an open system. Some flows enter the system from outside, others are generated in the system. One of such flows generated in the system is entropy production, which is non-negative. This circumstance turns the entropy balance equation into an inequality. Together with the other equations, it identifies a region of feasibility in the space of flows, the boundary of which corresponds to reversible processes. If, however, in one way or another it was possible to solve the problem of the minimum possible entropy production σmin under the restrictions imposed on the process, then the equations of thermodynamic balances with the condition σ≥σmin identify a region whose boundary is determined by processes of minimal dissipation. This region lies inside the region bounded by reversible processes.
Open system. Thermodynamic balances establish a connection between the flows of each substance, energy, and entropy that the system exchanges with the environment, as well as the occurrence of these quantities in the system and the rate of change of their quantity. We will sum up all flows further, considering the incoming flows as positive, and the outgoing flows as negative. We will divide the flows into convective and diffusive, marking the latter with the index d. Unlike the convective flow, the diffusion flow depends on the difference between the intensive variables of the system under study at the point where it enters or leaves and the intensive variables of the environment. In addition, we will use the following j—flow index, ej, vj—internal energy and volume of one mole of the corresponding flow, and Pj—its pressure, hj=ej+Pjvj—molar enthalpy, hdj—enthalpy in the flow entering by diffusion, qj—j-th heat flow, Na—power produced by the system.
Let us present the general form of balance (3)dEdt=∑jgjhj+∑jqdj+∑jqj−Na,(4)dNidt=∑jgjxij+∑jgdjxdj+∑ναiνWν,(5)dSdt=∑jgjsj∑jqdj−∑igdjμdijTdj+∑iνμiνniνTν+∑jqjTj+σ. Here niν=−αiνWν is the intensity of formation of the i-th substance in the ν-th reaction, Tν is the temperature of the ν-th reaction. If there are no diffusion flows, then(6)dEdt=∑jgjhj+∑jqj−Na,(7)dNidt=∑jgjxij+∑ναiνWν,(8)dSdt=∑jgjsj+∑jqjTj+∑iνμiνniνTν+σ, where the number of heat flows includes heat flows released or absorbed during chemical reactions, which depend on the reaction rate.
If we consider a stationary mode of the process, i.e., when dE/dt=dNi/dt=dS/dt=0, then the written equations from differential equations are transformed into final relations. When considering a cyclic process, balances can be written not for each moment in time, but on average for the cycle of the installation. Since the state of the system is the same at the beginning and end of the cycle, the total change in energy, amount of substance and entropy of the working fluid for the cycle are zero. In this case, balances are also reduced to a system of relations linking the average values of the terms on the right-hand sides of the equations for the cycle.
For closed systems consisting of several equilibrium subsystems, thermodynamic balances have the formE˙0=∑iE˙i,N˙0=∑iN˙i,S˙0=∑iS˙i, where i is the subsystem number, and index 0 refers to the system as a whole. In turn, E˙i, N˙i and S˙i are determined by the relations (3)–(5).
The requirement of minimum entropy production in processes of a given average intensity leads to the problem of such an organization of the process, in which the entropy production associated with it will be minimal [13]. Let us present a scheme for solving this problem for the scalar case.
Let two bodies contact, characterized by their extensive variables Y (volume, internal energy, entropy, number of moles of substance, …) and intensive u (temperature, composition, pressure, …). The difference in intensive variables leads to the appearance of an exchange flux J. The entropy increase over time L is equal to the integral of the product of the flux J and the driving force X, which, like J, depends on the intensive variables. Moreover, J and X always have the same sign. Dissipation is equal to the average rate of entropy increase. In this case, the intensive variables themselves depend on the extensive ones due to the equation of state. We will assume that the intensive variable u2 of the second of the two contacting systems can be changed in an optimal way, and u1 changes due to the influence of the exchange flow on it. The average value of the flow is given. The problem will take the (9)σ¯=1L∫0LJ(u1,u2)X(u1,u2)dl→minu2∈V under conditions (10)du1dl=φ(u1,u2),u1(0)=u10,(11)1L∫0LJ(u1,u2)dl=J¯. In these equations l is the duration of contact.
The value of L can be either fixed or subject to optimal selection.
The optimality conditions for this problem, as well as the optimality conditions for vector flows characterized by Onsager kinetics, are given in [14,15]. For the most important case, when (12)φ(u1,u2)=c(u1)J(u1,u2) we obtain (13)J2(u1,u2)=λ2∂J(u1,u2)∂u2:∂X(u1,u2)∂u2.
Let the organization of processes in the system correspond to the conditions of minimal dissipation and the corresponding production of entropy σmin, being substituted into the entropy balance equation, distinguishes the realizability region in the space of flows. The following statement is true [15]:
If the system productivity monotonically depends on the flow of mechanical, electrical power, or separation power generated by the system, and the flow of costs monotonically depends on the heat flow, then the region of realized modes in the productivity-costs plane is limited from above (Figure 1a). If productivity is determined by the heat flow, and costs are determined by mechanical or electrical power, then productivity monotonically increases with increasing costs (Figure 1b).
The first type includes heat engines, rectification processes, absorption refrigerators, etc. The second type includes compression refrigerators and heat pumps, membrane separation systems, electric heaters, etc.
Let there be a non-uniform thermodynamic system consisting of several subsystems (reservoirs, subsystems of finite capacity), the intensive variables of which at the initial moment of time differ from each other. In addition, the system has a working body that can establish contact with each of the subsystems. The contact function U(t) takes the value of one in the presence and zero in the absence of contact. The coefficients of heat and mass transfer for such contact are limited. It is required to organize the contacts of the working body with the subsystems in such a way as to extract the maximum work from the system in a given time τ. The extracted work is equal to the change in the internal energy of the system, and since the initial state is given, the goal of the solution is the maximum decrease in the internal energy in a limited time. In this case, the state of the working body at the beginning and end of the process is usually assumed to be the same.
The result of solving the problem are two statements [16]:
Statement In a thermodynamic system consisting of reservoirs and a working fluid with a given initial state, for any laws of heat and mass transfer, the maximum work extracted during time τ corresponds to a process for which
the vector of intensive variables u and contact functions U on the interval (0,τ) is piecewise constant, and the number of values it takes does not exceed r+m+2, where r is the number of conditions imposed on the final state of the subsystems and m is the dimension of the concentration vector;
at the beginning and end of the process, the intensive variables of the working fluid change abruptly to some optimal values corresponding to optimal pressures;
the entropy of the system grows on the interval (0,τ) as a piecewise linear function.
Statement 1 follows from the properties of solutions of averaged problems presented in the previous section. Depending on the given boundary conditions, the maximum work can be greater or less than zero. In the latter case, it corresponds to the minimum of the expended work. We emphasize that such a structure of the optimal process is characteristic of any kinetics of heat and mass transfer.
Corollary: When no restrictions are imposed on the composition and entropy of the working fluid and the increments of the extensive variables of the reservoirs at t=τ, the entropy of the system in the optimal process for any laws of heat and mass transfer increases at a constant rate, and the working fluid contacts the same reservoirs throughout the process.
In the presence of subsystems of finite capacity, the problem of maximum work turns out to be an optimal control problem with integer variables U(t). In this case, it is true that
Statement On each interval of constancy of the contact function between the working fluid and the subsystem of finite capacity, the law of change of the vector u(t) of intensive variables of the working fluid in the optimal process must satisfy the conditions of minimal dissipation.
In multi-flow heat exchange systems with given input temperatures and heat capacities of hot flows, total heat exchange coefficient and heat load for Newton’s law of heat exchange, the minimum of entropy production corresponds to such an organization in which at each point of contact the ratio m of absolute temperatures of the heated and heating flows is the same, as is the temperature T¯ of cold flows at the system outlet [17].
The following relationships are T¯=∑iTi0Wi−q¯∑iWi,q*(Ti0)=Wi(Ti0−T¯),α*(Ti0)=α¯Wi(lnTi0−lnT¯)∑iWi(lnTi0−lnT¯),m=1−1α¯∑iWi(lnTi0−lnT¯),σ¯=α¯(1−m)2m,α(Ti0)=q*(Ti0)=Wi=0,Ti0≤T¯.
Here Ti0, Wi0 are the temperature and heat capacity of the i-th hot stream, q is the heat flux, σ is the entropy production rate, and α¯ is the total heat transfer coefficient for the system. The last equality indicates that hot streams whose temperature does not exceed T¯ should be excluded from the system.
In the economy, economic agents (EA) function at the micro level, each of which is uncontrollable. Together, they form a macrosystem that is largely analogous to the thermodynamic one [18].
The analogy between thermodynamic and microeconomic systems is obvious and has attracted the attention of many researchers [19,20,21,22]. The reason for this analogy is that both systems consist of a large number of individually uncontrolled elements, economic agents (EA), i.e., they are macrosystems. In order to make such an analogy strict, it is necessary to introduce a variable for microeconomic systems that is an indicator of the irreversibility of the processes occurring in them. The proof of the existence of such an indicator is given in [16]. Here it is given with minor abbreviations.
In exchange processes, the EA buys and sells resources, simultaneously changing the stocks of resources and capital. Let us introduce the function U, the differential of which(14)dU=dM+∑ipidNi. We will call this function capitalization or total capital of the EA, since its change takes into account the change in both the basic resource M and the tied capital F=∑ipiNi. The vector p depends on the stocks of resources. The values pi have the dimension of the cost of the i-th resource in units of the basic resource.
Let some firm carry out an equilibrium exchange with the EA, exchanging some types of resources for others through purchases and sales. The exchange occurs reversibly and in such a way that the initial and final states of the EA in the space with coordinates Ni coincide. If in such a process the firm could extract a certain amount of the basic resource, this would mean that the possibilities of its extraction are not limited, since it can be obtained only from one EA, without causing any changes in its state or in the state of the environment. From the impossibility of such an economic “perpetual motion machine” it follows that r=const∮∑ipi(N,M)dNi=0. From this, in turn, it follows that there exists a function Q(N,r), whose partial derivatives with respect to Ni are equal to pi, and the differential has the formdQ=∑ipidNi+∂Q∂rdr. Here r is the valuation of the EA capital.
The expression (14) can be rewritten asdU=dM+dQ−∂Q∂rdr=d(M+Q)−∂Q∂rdr. Denoting M+Q=Y, and −∂Q∂r=γ, we obtain that(15)dU=dY+γdr.
Expression (15) depends on three Y,γ and r. These variables are related to each other. If this were not so, then on the plane with coordinates Y,r it would be possible to move from a given initial state to any other with constant capitalization by choosing γ. In reality, the possibilities of such a transition are limited.
Indeed, let point 1 correspond to the initial state on the plane Y,r, and point 2 to the final state, such that with reversible exchange with the environment with a constant value of r, the capitalization of the EA decreases (dU<0). If it were possible to get to the same point at dU=0, then the EA could organize a cyclic exchange process in which its state would change from point 1 to point 2 along the trajectory with dU=0, and would return from point 2 to point 1 with a constant value of r and an increase in capitalization (dU>0). So it could increase its capitalization by any amount in a reversible exchange, which is impossible. Thus, γ is a function of Y and r. The differential of the total capital is thus a Pfaffian form of two variables, which always has an integrating factor.
Recall that a Pfaffian form is a first-order differential form, i.e., the sum of products of functions of several variables by the differentials of these variables,dK=∑i=1nFi(x)dxi. If n=2 and the functions Fi are differentiable, then there will always be an integrating factor f(x) such that dS=f(x)dK is a total differential, i.e., S depends on x, and ∮dS=0. The factor f(x) is usually not unique.
Before finding the integrating factor p0(Y,r), let us clarify the dependence of γ on these variables. To do this, consider a system consisting of two EAs with identical estimates of the basic resource (r1=r2=r0). Capitalization U and variable Y and the basic resource reserves for the system are additive, so that U=U1+U2, Y=Y1+Y2. So that Q=Q1+Q2. Accordingly, γ=∂Q∂r=γ1+γ2. Thus, the dependence of γ on Y satisfies the superposition γ(Y1+Y2,r)=γ(Y1,r)+γ(Y2,r) and has the form(16)γ(Y,r)=G(r)Y+l(r),
Thus, for the capitalization differential we obtain the expressiondU=dY+(G(r)Y+l(r))dr. The integrating factor can be any function p0 such thatdS=p0dU=p0dY+p0(G(r)Y+l(r))dr is a total differential. Since the mixed derivative of the function S does not depend on the order of differentiation, the function p0 must satisfy the condition(17)∂S∂r∂Y=∂S∂Y∂r→∂p0∂r=∂[p0(G(r)Y+l(r))]∂Y.
We will look for p0 as a function of r, then the condition (17) leads to the equation(18)∂p0∂r=p0G(r). One of the solutions to this equation is the function p0(r)=exp[G(r)]. It uniquely depends on r and can also serve as an estimate of the base resource.
Thus, there is some function of state (extensive variables) S(N,M) such that its differential has the form(19)∮dS=0. In a process with non-zero flows, when prices do not coincide with estimates, the function S (the irreversibility index) always increases.
The estimates of resources can be expressed in terms of the function S as(20)p0=∂S∂M,pi=∂S∂Ni/∂S∂M,i=1,2,… In this case, the assessment of the basic resource p0 is positive for all economic agents, and pi can be negative if the i-th resource requires, for example, utilization or storage costs. Cash can be considered one of the resources with a reserve of Nj, the assessment of which pj=1/p0. The function S(N,M) is called the welfare function or, more briefly, welfare.
When exchanging resources between economic agents, the conditions of voluntariness must be observed, which consist in the fact that none of the welfare functions Sν decreases (the exception is associated exchange, charity). The conditions of voluntariness make it impossible to directly exchange one type of resource if its assessments by the EAs in contact with each other are of the same sign. Such an exchange becomes possible only if there is an intermediary.
If the welfare function is measured in national currency, then the value p0 characterizes the value for the EA of the international currency and has the dimension [unit of national currency/unit of the basic resource]. The assessment p0 of the basic resource in the processes of exchanging money by an economic agent on the foreign exchange market plays the same role as the assessment p in the exchange of resources.
The description of economic systems becomes formally close to the relations of thermodynamics if we introduce an ‘economic temperature’,T=1p0. This notation was introduced in [23], the interpretation of this quantity, called liquidity, and its properties were given considerable attention in [24].
When the scale of the EA changes without changing its properties, the stocks of resources and capital are proportional to the scale. It is natural to assume that the welfare function changes in the same way, i.e., it is an extensive variable like N and M. In this case, the S function is homogeneous of the first order, and its derivatives with respect to N and M (estimates) are homogeneous functions of the zeroth order. According to Euler’s theorem for homogeneous functions, it can be written in the form(21)S(N,M)=p0(M,N)∑ipi(M,N)Ni+M. The dependence p(N,M) can be found by formulas (20) and experimentally by the behavior of the EA in exchange processes.
If the existence of a welfare function S is postulated, then the estimates are determined through an extremal problem in which they enter as (22)S(N,M)→max/∑ipiNi+M=U, where the value of U is fixed.
The Lagrange function of the problem (22) isL=S(N,M)−p0∑ipiNi+M, where p0 is an indefinite Lagrange multiplier. The stationarity conditions for L lead to the relations(23)pi(N,M)=∂S∂Ni/∂S∂M, which coincides with the expression (20). The function S is often assumed to be twice continuously differentiable, monotonically increasing in each of the arguments, strictly convex upward and equal to zero at the origin. Its partial derivatives tend to infinity when the corresponding variable tends to zero.
As a consequence of these assumptions, the solution to the problem (22) is unique and positive, and pi decreases with an increase in Ni. With this description, the EA is similar to a thermodynamic subsystem of finite capacity. For an economic reservoir, the resource estimates pi and p0 are constant, and the welfare function is linear.
Note that in the general case, the welfare function of a system is not equal to the sum of the welfare functions of its constituent subsystems. Moreover, for each subsystem, the welfare function can have its own dimension. In contrast, the capitalization of each of the subsystems, the volume of the basic resource and the associated capital have the same dimension and their sum makes sense for the system as a whole under the conditions of internal equilibrium of each of the EA included in it.
Let’s write the differential of the function S,(24)dS=p0dM+∑i=1npidNi=p0dU. Solving (24) with respect to dM, we obtain(25)dM=dSp0−∑i=1npidNi. From equality (21) it follows that(26)M=Sp0−∑i=1npiNi,(27)dM=dSp0+Sd1p0−∑i=1n(pidNi+Nidpi).
Comparing equalities (27) and (25), we obtain a relationship linking the estimates of resources and capital(28)Sd1p0−∑i=1nNidpi=0. Similarly, comparing the differential S found from equality (21) with the expression (24), we obtain(29)Mdp0+∑i=1nNid(p0pi)=0. Conditions (28) and (29) follow from the existence of the function S and its homogeneity. They are economic analogues of the Gibbs-Duhem equations. For example, it follows from them that if the state of the system changes so that the resource estimates are constant, then the capital estimate p0 is also unchanged.
Due to the symmetry of the matrix of second derivatives for a twice differentiable function, the sensitivity of estimates to changes in resource stocks and capital are related by the equalities(30)∂(p0pi)∂Nj=∂(p0pj)∂Ni=∂2S∂Ni∂Nj,(31)∂p0∂Nj=∂(p0pj)∂M=∂2S∂M∂Nj. From conditions (30) and (31) it follows that(32)∂pi∂Nj+pi∂pj∂M=∂pj∂Ni+pj∂pi∂M,i,j=1,…,n. Equalities (30) and (31) are the economic analogue of Maxwell’s relations.
Let us return to the cyclical process of interaction of a firm with one EA and require that the average intensity of exchange processes be fixed. Then the firm when purchasing a resource will be forced to raise prices compared to the estimates pi, and sell at prices that are lower than these estimates. The capitalization of the EA will increase, since(33)ΔU=∮∑(pi(N,M)−ci)dNi≥0, and the firm will suffer losses in comparison with the reversible process in the amount of ΔU.
We will call the intensity of the firm’s losses due to irreversibility,(34)σ(t)=∑ini(pi,ci)(pi−ci)≥0, the dissipation of capital associated with the irreversibility of the resource exchange process. In this case, it has the meaning of trading costs.
The condition (33) of non-decreasing capitalization (and therefore welfare) in economic exchange is an analogue of the Clausius integral. Similarly, the law according to which at the contact of two EA, the resource passes from the EA for which its valuation is lower to the EA for which its valuation is higher, and at the same time the total value of the tied capital does not decrease (Δ(F1+F2)≥0), is an analogue of the second law of thermodynamics. This allows us to construct an irreversible microeconomics, in many ways similar to finite-time thermodynamics.
The laws of conservation of matter and energy in microeconomics are the laws of conservation of resources. Here we will focus on the economic analogy of the second law of thermodynamics.
There are several formulations of the second law of thermodynamics, each of which can be considered a consequence of the others. Let us discuss the analogs of some of these formulations in microeconomics.
Among the numerous formulations of the second law, we will highlight the formulation of Clausius with Planck’s ≪Heat by itself cannot pass from a cold body to a hotter body without leaving other changes≫, and also the formulation of Leontovich: ≪It is impossible to construct a device, as a result of the action of which positive work would be produced only due to the cooling of one body without any other changes≫.
In microeconomics, the following statements correspond to the above *A scalar resource flow cannot pass from an EA with a higher valuation to an EA with a lower valuation without leaving other changes.*It is impossible to extract capital by exchanging resources with one EA without any other changes.
M. Planck formulated the following statement as a consequence of the second law of ≪In an irreversible process in a closed thermodynamic system, entropy can only increase and the exergy of the system can only decrease. The equilibrium state of such a system corresponds to the maximum entropy and the minimum performance under conditions corresponding to the imposed restrictions≫. Similarly, for economic systems resource exchange processes in isolated microeconomic systems proceed in such a direction that the total tied capital of economic agents increases and reaches a maximum, and the potential ability of extracting the basic resource (profitability) decreases and reaches a minimum, under the restrictions imposed on the system, which include the conditions of voluntariness.
The notations adopted here and listed in Table 1 T− and T—temperatures of the reservoir and the system in contact with it, p−—valuation of the resource on the market, c—price of the resource set by the monopolist firm, N—stock of the resource, U—internal energy of the system and total capital, q and n—flows of heat and resource, M and F—basic resource and tied capital.
Let us write the balance equations for a heterogeneous economic system exchanging resource and capital flows with its environment. In this case, we assign index i to the i-th type of resource, and index j to the j-th subsystem. External resource and capital flows entering the system will be considered positive, and those leaving it will be considered negative. These flows are divided into two those forced into the system and changing their intensity over time under the influence of external factors, and those dependent on prices set by external sellers and buyers, and on resource valuations in the subsystem under consideration. The first of these, by analogy with thermodynamics, will be called convective and marked with index k, and the second diffusive (index d). In addition, the subsystem can transform some types of resources into others.
The balance equation for the i-th resource is(35)Ni˙=∑jnijk(t)+nijd(pj,cj)+Wj(pj)αij,i=1,2,… Here the summation is carried out over all subsystems; Wj(pj)—the intensity of the resource transformation process in the j-th subsystem, and the coefficients αij>0 if the i-th resource arises in the j-th subsystem and αij<0 if this resource is spent, they determine how much of the i-th resource arises (is spent) per unit of time; cj—the price vector during the exchange of the j-th subsystem with the environment.
The balance equation for the basic resource is(36)M˙=∑jmjk(t)−∑icijnijd(pj,cj) and the balance sheet equation for tied capital is(37)F˙=∑i,jpij(Nj,Mj)(nijk(t)+nijd(pj,cj))+σ, where capital dissipation σ is equal to(38)σ=12∑j∑νnjν(pj,pν)(pj−pν)+∑jWj(pj)Aj. Here pj and pν are vectors of resource estimates for the j-th and ν-th subsystems with components pij and piν, Aj=∑iαijpij, and njν=−nνj is a vector function of resource flow.
The first term in (38) is a change in associated capital due to resource exchange, and the second is due to the transformation of resources. pj and pν are vectors of resource estimates for contacting subsystems with components pji and pνi. Finally, njν=−nνj—vector function of the exchange flow.
The value σ(p)≥0, so that in a heterogeneous open system in the absence of convective flows, the bound capital in the outgoing flows is not less than in the flows, entering the system. The equal sign corresponds to a homogeneous system.
As in thermodynamics, the condition σ(p1,p2)≥0 together with the balance relations (35)–(37) distinguishes the boundary of the feasibility region of an economic system in the class of reversible processes. The conditions imposed on the intensity of a particular flow allow us to determine the value σmin>0, achievable under these conditions. In this case, the area of feasibility narrows, since instead of the inequality σ≥0, the inequality σ≥σmin.
If Δpjν=pj−pν is small, and the kinetic function niν is differentiable with respect to the set of arguments, then σ is a positive-definite quadratic form of the variables Δpiνj. In the steady state, the right-hand sides of Equations (35)–(37) are equal to zero. In cyclic mode, whenN(0)=N(τ),M(0)=M(τ),F(0)=F(τ) the integrals of the right-hand sides of these equations are equal to zero.
The condition of voluntary exchange imposes restrictions on the set of possible states of an open system. To find these restrictions, we write down the balance relations for the resource and capital for j-th subsystems;(39)Nij˙=nijk(t)+nijd(pj,cj)+∑νnijν(pj,cjnu)+Wj(pj)αij,i=1,2,…(40)Mj˙=mjk(t)−∑i[cijnijd(pj,cj)+∑νcijnunijν(pj,cjν)]. Here cjν is a vector of intermediate prices with components cjνi, which is determined from the continuity condition of flows,(41)n˜νji(pj,cjν)=−n˜jνi(cjν,pν)=nνji(pj,pν),i=1,…,m. For any subsystem due to the conditions of voluntariness,Uj˙=Sj˙/p0j=∑i(pijNij˙+Mj˙)≥0, which, taking into account (39) and (40), leads to restrictions on estimates and flows in the (42)∑i[pij[nijk(t)+Wj(pj)αij]+∑ν[nijν(pj,cjν)(pij−cijν)]++nijd(pj,cj)(pij−cij)]+mjk(t)≥0,j=1,2,….
For an isolated system, there are no external flows and the balance relations (35)–(37) take the form(43)Ni˙=∑jWj(pj)αij,Ni(0)=∑jNij(0),(44)M˙=∑jMj˙=0,(45)F˙=∑jFj˙=σ≥0. In equilibrium, tied up capital is maximal, and the flows nij and the rates of resource transformation Wj(pj) are equal to zero. In this case, the equilibrium distribution of the basic resource M¯ depends on the kinetics of resource exchange between subsystems.
In thermodynamics, the equilibrium conditions can be obtained from the problem about the maximum of the total entropy of the system for a given total value of one or another extensive variable of the subsystems. For example, if the total volume of the subsystems is given, and the remaining extensive variables are unchanged, then the maximum of the total entropy corresponds to such a distribution of volumes, in which the derivatives of the entropy of each system with respect to its volume are the same, and this leads to equality of pressures. Similarly, the distribution of the total amount of thermal energy in equilibrium leads to equality of temperatures, etc.
In microeconomics, the conditions for the equilibrium distribution of resources can also be derived from the requirement of the maximum total tied capital of the system, given the total stock of resources. This leads to a distribution in which the derivatives of the tied capital of each subsystem with respect to the stock of its allocated resource (resource valuation) are the same.
For each j-th subsystems of the isolated system(46)Nji˙=∑νnνji(pj,pν)+Wj(pj)αij,(47)Mj˙=−∑ν,in˜νji(pj,cjν)cjνiν,j=1,…,n. Here cjν are intermediate prices determined by conditions (41). Thus, the price vector, and hence the right-hand side of the Equation (47) depends on the type of kinetics n˜νj of demand and supply of resources.
After determining from (41) the intermediate price cjν(pj,pν) and substituting it into n˜jν and n˜νj, these functions turn out to be the same and equal to the kinetic function with components njν(pj,pν), appearing in expressions (38) and (46).
As an example, let us consider the kinetics of the exchange of two EAs linear with respect to the difference between the price and the valuation and calculate the dissipation of capital σ in this (48)n1(p1,c)=a1(p1−c),(49)n2(p2,c)=a2(p2−c). Let us find c(p1,p2) using the condition −n1=n2=n:a1(p1−c)+a2(p2−c)=0, where(50)c=a1p1+a2p2a1+a2,(51)n(p1,p2)=−n1(p1,c(p1,p2))=a¯(p2−p1),(52)a¯=a1a2a1+a2. Then the resource exchange dissipation is(53)σ(p1,p2)=(p2−p1)a¯(p2−p1)=a¯(p2−p1)2=n2(p1,p2)a¯.
The equilibrium distribution of the basic resource M¯ is determined by the exchange kinetics, and the equilibrium stocks of resources N¯ also depend on M¯, since they satisfy the equality(54)pj(Mj¯,Nj¯)=pν(Mj¯,Nj¯)=λ,∀j,ν.
In a number of cases, the set Q of values of the vector M¯ that can be reached from a given initial state of the system for various functions n˜(p,c) of demand and supply is of interest.
For each j-th subsystem, the minimum capital increase ΔMj=Mj¯−Mj0 is achieved when the exchange with other subsystems is conducted at prices cjν arbitrarily close to pj, i.e., reversibly, so that the value Mj¯min can be found by the condition Sj(Nj¯,Mj¯)=Sj(Nj0,Mj0).
The maximum value of Mj¯max corresponds to an exchange at which cjν are arbitrarily close to pν. Thus, in the space with coordinates Mj, it is possible to construct a parallelepiped with boundaries Mj¯min≤Mj¯≤Mj¯max. The section of this parallelepiped by the plane(55)∑jMj¯=∑jMj0 selects the set Q of all possible equilibrium distributions of the basic resource between subsystems for different kinetic functions.
Exchange flows in macrosystems depend on intensive variables of subsystems, and the rate of change of extensive variables is determined by these flows. This leads to problems in which extensive variables are not included in the right-hand sides of differential equations that determine their change. For example, the change in the internal energy of two contacting subsystems is characterized by the differential equation(56)dU1dt=q(T1,T2).
The internal energy itself is not included in the right-hand side of this equation. If the values of the internal energy are subject to constraints U1(0)=u10,U1(τ)=u1τ, then these constraints correspond to the condition imposed for a given value of τ on the average value of the current1τ∫0τq(T1,T2)=1τΔU1.
The Equation (56) are called Lyapunov-type equations by L.I. Rozonoer [3]. When solving problems, these equations can be discarded, replacing them with conditions imposed on the average rate of change of state variables, which leads to problems of averaged optimization. The same is typical for cyclically occurring processes in systems with intermediaries. In such systems, their efficiency indicators are averaged over a time equal to the duration of the cycle, which also leads to averaged problems. Below we present a proof of the theorem that defines the structure and type of optimal solutions to such problems.
The mean value of a continuous scalar function f(x(t)),t∈[0,τ],x∈V⊂Rn can be calculated over time as(57)ft(x)¯=1τ∫0τf(x(t))dt or over a set as(58)fp(x)¯=∫Vf(x)p(x)dx. The function p(x) is called the distribution density. In the case where x(t) is a random function, p(x) is the distribution density of a random variable. It is non-negative and its integral on V is equal to one. In our case, x(t) is a deterministic function. Let us dwell in more detail on the properties of the function p(x), such that the results of averaging by Formulas (57) and (58) are the same.
We will assume that the variable x is scalar, the set V here and below is bounded and closed, and introduce the function θ(x0),x0∈V, equal to the total duration of those time intervals t for which x(t)≤x0. Obviously, this function does not exceed τ. By P(x0) we denote the ratio θ(x0)τ, i.e., the fraction of the interval [0,τ] for which x(t)≤x0. This function increases monotonically with x0, changing from zero to one. It is similar to the distribution function of a random variable.
The distribution density is(59)p(x0)=dP(x0)dx0=1τdθ(x0)dx0=1τ1∑νdxνdtxν=x0. Here the interval θ increases with x0 for any sign of the derivative for those values xν of the function x(t) for which it is equal to x0.
If for some value x0 the function x(t) is constant during a fraction γ of the interval [0,τ], then the function P(x0) experiences a jump of magnitude γ at this point, and the distribution density at it is equal to γδ(x−x0).
Example 1—Linear functions. Let x(t)=htτ. Then by Formula (59) we obtain p(x)=1h=const. The same distribution density corresponds to all triangles with base [0,τ] and height h.
Example 2—Piecewise constant functions. These functions take discrete values xi, each during a fraction γi of the interval [0,τ]. Any such function, according to Formula (59), corresponds to the distribution density(60)p(x0)=∑iγiδ(x−xi),γi>0,∑iγi=1. The order in which a piecewise constant function assumes one or another of the possible values does not matter.
From these examples it is clear that to each function x(t) there corresponds a distribution density of its values p(x), defined on V, and to each distribution density there corresponds an arbitrary number of functions x(t), for which fp(x)¯=ft(x)¯. An exception is the distribution density of the form p(x)=δ(x−x1). In this case, the corresponding function x(t)=x1=const on the entire interval [0,τ] and it is unique.
Let us consider the case when the function f depends on several, for simplicity on two, variables x1(t) and x2(t). In this case, the distribution function P(x0) of the values of the vector x is a fraction of the interval [0,τ] for which two inequalities are x1(t)≤x10 and x2(t)≤x20. This function monotonically increases with the growth of each of the arguments, when the first component of the vector x0 is maximum (p1(x1)=1), it is equal to P(x20), and its derivative is equal to the distribution density p(x1max,x2)=p2(x2); similarly in the case when x2=x2max,p(x2max,x1)=p1(x1). The functions x1(t) and x2(t) are independent of each other so that p(x1,x2)=p1(x1)p2(x2).
The desired solution to the averaged optimization problem is the distribution density p*(x) of the vector x on the set V of its admissible values. To implement this solution in time, one must find one of the possible functions x(t) that have the distribution p*(x). The solution to this last problem is significantly simplified by the features of the optimal solutions p*(x), proven in the next section.
We will denote the averaging operation by a bar drawn over the averaged function or vector. Thus,x¯=∫Vxp(x)dx,f(x)¯=∫Vf(x)p(x)dx.
The simplest problem of averaged optimization is the problem of maximizing the average value of a scalar function f(x) for a given average value of its (61)f(x)¯→max/x¯=x0,x∈V⊂Rn or in more detailed notation(62)∫Vf(x)p(x)dx→max/∫Vxp(x)dx=x0,p(x)≥0,∫Vp(x)dx=1. The unknown in this problem is p(x) (the distribution density of the vector of unknown variables). This function is non-negative and its integral on the set V is equal to one.
Carathéodory’s theorem [25] on convex hulls of sets states that any element of the convex hull CoD of a compact set D in a Euclidean space of dimension n can be represented as an element of a simplex that has at most n+1 vertices (base points), each of which belongs to D.
In particular, a subgraph of a function f(x) can be a set D. The convex hull of a function is called the convex hull of a subgraph. A function that depends on n variables is the boundary of a set in the space Rn+1 of dimension n. The base points are certainly on this boundary, which means that their number does not exceed n+1. Below we will call Carathéodory’s theorem the theorem on convex hulls of functions.
The ordinate of the convex hull of the function f0(x) at the point x0, which belongs to the convex hull of the set of definition of the function, is the value of the (63)f0(x)¯→maxp(x)/xi¯=xi0,i=1,n¯,x∈V⊂Rn where V is compact.
According to Carathéodory’s theorem, the optimal solution to this problem isp*(x)=∑j=0nγjδ(x−xj),γj≥0,∑i=0nγj=1. That is, the optimal distribution is concentrated at no more than the (n+1)-th base point.
This fact allows us to rewrite the problem as a nonlinear programming (NP) (64)∑j=0nγjf0(xj)→max/∑j=0nγjxj=x0,xj∈V⊂Rm,∑j=0nγj=1,γj≥0. The variables are the base vectors xj and the vector of weight coefficients γ. We use the Kuhn-Tucker theorem [25] to solve If y* is a solution to a nonlinear programming problem(65)f(y)→max/φi(y)≤0,yj≥0,i=1,2,…,m,j=1,…,n, then there exists a nonzero vector of factorsλ=λ0,…,λm(λ0is equal to0or1,λi≤0fori>0), such that for the Lagrange functionR=λ0f(y)+∑i=1mλiφi(y) the following conditions (66)∂R∂yjy=y*=0ifyj*>0;∂R∂yjy=y*≤0ifyj*=0, and(67)λi=0ifφi(y*)<0;λi≤0ifφi(y*)=0.
For problem (64), the Lagrange function has the form(68)R=∑j=0nγjf0(xj)+∑i=1nλixij−Λ, where Λ is the Lagrange multiplier corresponding to the condition that the sum of the weight coefficients is equal to one.
The Kuhn–Tucker local unimprovability conditions with respect to weight coefficients lead to the requirements(69)R0(xj,λ)=f0(xj)+∑i=1nλixij<Λifγj=0,R0(xj,λ)=f0(xj)+∑i=1nλixij=Λifγj>0,j=0,…,n+1. Here R0 is the Lagrange function of problem (64) in the absence of averaging. Such a problem will be called original. Thus, we have
In the NP problem with averaging of functions defining relationships between variables, it is required to achieve the maximum of the mean value of the function f0(x) on the set V of admissible values of x, provided that the mean value of the vector function f(x)=(f1(x),…fi(x),…fm(x)) is zero. Formally(70)f0(x)¯→max/fi(x)¯=0,i=1,…,m,x∈V∈Rn, yielding.
Note that here and below, the conditions in the form of the maximum principle do not require the functions defining the averaged problem to be smooth in x; the set V may not be simply connected.
Example 1. Consider a system consisting of an electric motor, a pump it rotates, and a tank. The motor consumes power n on which the pump performance g depends. The dependence g(n) is shown in Figure 2. It is required to find the mode for which, for a given average expended power n¯, the average productivity g¯ is maximum.
This is a problem about the ordinate of the convex hull of the function g(n) at the point n¯. The number of base points is equal to two, one of them is the origin, and the second n1, is determined by the condition of the maximum of the Lagrange function R=g(n)+λn at it, the same as at n=0. Eliminating λ from the conditions for the maximum of the Lagrange function and the requirement that this maximum be zero, we arrive at the equation for n1:g(n)n=dg(n)dn. There are as many optimal implementations of this solution in time as you like. For each of them the pump power takes the value zero for one fraction of the period, and n1 for the remaining fraction, and the fraction of the interval τ for which n=n1 is equal to 1−n¯n1. The maximum value of the interval τ is determined by the value of the capacity G, it is equal toτmax=2Gg(n1).
The value of the problem isg¯*=g(n1)1−n¯n1. It does not depend on G. When the capacity tends to zero, the sliding mode becomes the optimal solution.
In averaged problems there can be two types of variables, random and deterministic. There is no averaging over the variables of the second type. Let us consider the NP problem, in which there is no averaging over some variables.
A problem with averaging over some variables will take the form(76)f0(x,y)¯→max/fj(x,y)¯=0,x∈V⊂Rn,y∈Vy⊂RK,j=1,…,m. The functions f0,…, fm are continuous and continuously differentiable with respect to the set of arguments, the bar corresponds to averaging over x∈V, and the sets V and Vy are closed and bounded.
For any y, this problem is an averaged nonlinear programming problem (70), and therefore, by Theorem 2, the optimal density function of x is concentrated at most at the (m+1)-th base point, so that p*(x)=∑0mγjδ(x−xj), and there exists a nonzero vector λ such that at each of these points the Lagrange function of the original problem(77)R=∑j=0mλjfj(x,y),x∈V⊂Rn,y∈Vy⊂RK is maximal with respect to x.
The Lagrange function of problem (76), in which the density function x is p*(x), has the form(78)R*¯=∑j=0mλj∑i=0mγifj(xi,y),xi∈V⊂Rn,y∈Vy⊂RK.
For any density function p(x) of randomized variables, problem (76) is a nonlinear programming problem and, according to the Kuhn-Tucker theorem, there exists a nonzero vector λ with components λ0=(0;1),λj,j=1,…,m, such that for the function (77) the optimal solution satisfies the conditions of local non-improvability with respect to y, in particular, the stationarity conditions(79)∂R*¯∂ylδyl≤0,l=1,…K. Here δyl is the admissible variation of yl.
In FTT the specificity of problems on the optimal exchange of resources and problems with cyclically changing variables of the working fluid or mediator leads to averaged problems of nonlinear programming. Their optimal solutions are of a switching nature between the basic values of the variables. The maximum number of basic values is one more than the number of averaged conditions, and in each of them the Lagrange function of the non-averaged problem reaches a maximum.
The presented results of solving the problems of “Thermodynamics at finite time” illustrate the nature of its applications for systems of various natures, but do not pretend to be an exhaustive presentation. Not all the problems of FTT have been solved. In particular, the problem of the conditions of minimal dissipation for vector flows whose kinetics differs from Onsager’s awaits its solution. It is not clear what to consider the efficiency of a machine receiving energy from several sources. It is not clear what the irreversibility indicator looks like for macrosystems of a social nature, etc.