On the Impact of the Dynamics of Heat Transfer of the Thermal Machine Working Fluid and Heat Sources on the Shape of the Boundary of the Set of Realizable Regimes

From the point of view of finite time thermodynamics, the performance boundaries of thermal machines are considered, taking into account the irreversibility of the heat exchange processes of the working fluid with hot and cold sources. We show how the dynamics of heat exchange affects the shape of the optimal cycle of a heat engine and its performance, in particular, energy conversion efficiency in the maximum power mode. This energy conversion efficiency can depend only on the ratio of the heat transfer coefficients to the sources, or not depend on them at all. A class of dynamic functions corresponding to “natural” requirements is introduced and it is shown that, for any dynamics from this class, the optimal cycle consists of two isotherms and two adiabats, not only for the maximum power problem, but also for the problem of maximum energy conversion efficiency at a given power. Examples are given for calculating the parameters of the optimal cycle for the cases when the heat transfer coefficient to the cold source is arbitrarily large, and for dynamics in the form of a linear phenomenological (Fourier heat transfer) law.


Introduction
Systems converting heat into work have been the main object of thermodynamic research since the time of S. Carnot, who gave an upper boundary for the efficiency of such a transformation.He considered a heat engine as a system consisting of a working fluid which is in contact with two heat sources, a hot and a cold one.Carnot efficiency (the ratio of the work produced to the heat received from the hot source) was achieved in the limiting reversible process when both the power generated by the machine and the heat flows tend to zero, or when the heat transfer coefficients between the working fluid and the heat sources tend to very large numbers.The friction effects and inhomogeneities of the temperature field inside the working fluid leading to irreversibilities were negligible.In real heat engines, the power is fixed, the coefficients of heat exchange is limited and their efficiency is significantly less than the efficiency of Carnot.We, as is customary in finite time thermodynamics, consider each element of the system under consideration of internal equilibrium, and we consider the associaton with the entropy production processes, concentrated on the boundaries of the subsystems during their contacts with each other.
With the development of nuclear power, when the part of fuel cost in generated energy cost per kilowatt is much less, and the part of capital costs is much more than for non-nuclear thermal power plants (see [1]), the problem of the maximum power of thermal machines with limited coefficients of heat exchange with heat sources arose.I.I.Novikov assumed in his publication [2] that the cycle of maximum power, like the Carnot cycle, consists of two isotherms and two adiabats.He showed that the maximum power is limited, found the corresponding value, depending on the coefficient of heat transfer with sources and the corresponding efficiency, that turned out to be related to the Carnot efficiency η c by the expression Novikov assumed that the dynamics of heat exchange (the dependence of heat flow on the temperature of the contacting bodies) is linear with fixed heat transfer coefficients.The heat flux q i from the i-th source is proportional to the temperature difference between the source T i and the working fluid T: where i denotes the index of an arbitrary source and α i is the transfer coefficient.This kind of dynamics is often called Newtonian.With a limited coefficient α i , the heat exchange is irreversible.Note that η N as the Carnot efficiency does not depend on the heat transfer coefficients with the sources, and that is not obvious.Novikov's conclusions was later independently confirmed by Curzon and Ahlborn [3], assuming that the heat transport to the cold heat bath is irreversible, too.
Later, the problems of the maximum power of the thermal machine and its maximum energy conversion efficiency at a given power were considered in [4][5][6] and it was shown that for any heat transfer dynamics in the maximum power case, the optimal cycle consists of two isotherms and two adiabats.In the case of maximum energy conversion efficiency at a fixed power, or for the maximum power at a fixed heat consumption, the optimal cycle can consist of no more than three isotherms and three adiabats.An example of such a cycle was given in [7] for a somewhat exotic form of dynamic dependencies.These results were based on the methods of averaged optimization [8,9].The problems of thermodynamic analysis of heat conversion to the work of separation are also considered in [10][11][12][13][14][15][16][17] as well.There are a number of papers in which similar problems are considered, for example, [18], where the relationship between efficiency at maximum power and Carnot efficiency for a particular type of kinetic equations is considered.Another problem is given in [19], where reversible cycles with a finite heat capacity of the working fluid are considered.The duration of the cycle in this task does not appear.
There are the following problems in connection with these studies: 1. Find the dynamic laws of heat transfer with sources for which the efficiency of the machine in the maximum power mode does not depend on the dynamic coefficients 2. Determine the laws of dynamics in which the cycle corresponding to the maximum efficiency at a given heat consumption consists of two and three isotherms.3.If the optimal cycle consists of three isotherms, then sources two of them are realized?This work is devoted to solving these problems.

Averaged Optimization and Convex Hulls of a Function
The working fluid of a heat engine changes its parameters cyclically in time or in space.In the latter case, it circulates in a closed loop, receiving and giving heat to sources and doing work, therefore, all parameters of flows (power, heat) are necessary to average per cycle.This leads to the fact that it is convenient to use the methods of average optimization for the cycle analyzing.Now we will give certain information about these methods.

Convex Hulls of a Function
The convex hull a function f (x) on the set x V is the maximum of the mean value of the function with a given mean value of its argument.The formal description of this problem: Here, the angle brackets correspond to averaging, P(x), the function of measure.This function is non-negative and has an area equal to one.Exactly the measure function is the required solution to problem (3).The Caratheodory theorem [20] reduces this problem to the common task of nonlinear programming.It approves that the measure functions have a view at the optimal solution.The n is the dimension of vector x.The values of x ν of the vector x are called basic.Thus, the optimal measure is concentrated no more than at the (n + 1)-th base point on V, and γ ν are the weight coefficients of the corresponding values of x and f (x) After the substitution of the measure in the form (4) into the problem (3) it takes the form: This is a finite-dimensional problem of nonlinear programming.For functions of one variable, there are no more than two base values.The convex hull either coincides with the graph of the function f (x) (for one base value) or passes above the graph and represents a segment of a straight line tangent to the function f (x) in two base points.Similarly, for functions of two variables, a convex hull, if it does not coincide with f (x), it may either be a line segment (for two points), or a region of the plane (for three base points).

Averaged Problem of Conditional Optimization
In [8,9] and very briefly in Appendix A, the averaged problem of conditional optimization was considered: As proven for this problem, from the theorem Caratheodory it follows that the measure function on the optimal solution have a form where the m is the dimension of the vector function f (x) (the number of conditions).The values of x ν of the vector x are also called base values here.Their number, in contrast to the problem of a convex hull of function, does not depend on the dimension of x and no more than one more than the dimension of f (x).Thus, in the problem with a single average condition, the number of base values is not more than two, and in a problem with two conditions, not more than three values.As will be shown below, the problem of the maximum power of the heat machine has one, and the problem of maximum efficiency for a given power has two averaged conditions.

"Natural" Dynamic Functions and Optimal Cycles Structure
We will call natural dynamic dependencies relating the heat flow q with temperatures T i , T source and the working fluid, if they satisfy the following conditions: 1. q i (T i , T) are continuous and differentiable with respect to their arguments everywhere except for > 0 and ∂q i ∂T < 0 hold true (these are the monotonicity conditions).The same applies to the signs of the derivatives T i , T on q i .2. The sign of the heat flux coincides with the sign of the temperature difference T i − T (the second law of thermodynamics).q i (T i , T) = 0 at T i = T holds true. 3. The dynamic dependence can be written in the form q i (T i , T) = α i f i (T i , T), where the heat transfer coefficient α i > 0 is proportional to the contact surface of the working fluid with the source.For example, the dependency of type q i (T i , T) = α i f i (T i , T) + β i r i (T i , T) does not match this requirement.This type of dynamic was investigated in [7].The optimal cycle consisted of three isotherms and three adiabats in this work.
Cyclicity in the state of the working fluid causes what can distinguish two half-cycles, heating and cooling; the first of which, the working fluid, is contacted with hot and secondly with a cold source.

Statement of the Problem and Solutions Structure Corresponding to the Boundary of Realizability Region
The area of realizability of a thermal machine is a region in a plane whose axis of ordinates is the machine's power, and whose abscissa is the heat flux taken away from the hot source.This region is for specific forms of dynamic dependence and has been built in many papers, see [21] and others.We shall consider it peculiar properties, not specifying the form of dynamic functions q + , q − , which determine the heat transfer from hot and cold sources to the working fluid.
We assume a working fluid heat engine uniform temperature at each time or, in each section.Its extensive variables, such as internal energy and entropy, vary cyclically.
Let q( T) = α f ( T) be the heat flux of exchange of working fluid with sources, which changes the sign and shape of the dynamic function depending on the sign ∆ i = T i − T.
Each point of the boundary of the realizability region in the working load range is the solution of the following problem: Given a power p of the heat machine, find such a temperature vector T = {T i , T} of sources and the working fluid so that the energy conversion efficiency of the machine η = p q + is maximal.This problem is equivalent to problems about the maximum power of the machine at a given value of the average heat flux q + , taken from the hot source.
The formal optimization takes the form under conditions q T = 0; Here, q + > 0 is the heat supplied to the working fluid from the hot source, and q − < 0 is the heat drawn from it to the cold source.The union of these two heat fluxes is the heat flux q between the working fluid and the sources.So, q + is equal to q, when q is positive, which can mathematically be expressed as q + = 0.5(q + |q|).Analogously, q − is equal to q, when q is negative, which can be expressed as q − = 0.5(q − |q|).
The condition (9) follows from the entropy balance equation (the entropy change rate of the working fluid is q T , and its entropy does not change at a cycle).
The set of possible values of the vector T is defined by the fact that the range of variation of working fluid temperature T lies inside the range of variation of sources temperature T i .The range of variation of T i takes the values T + and T − .
The problem (8)−( 10) contains the averaged constraint ( 9) and the restriction on the average difference value between the input and output heat fluxes (machine power).From the results of the average optimization in [8,9] it follows that: 1.The optimal solution of this problem can take no more than three (base) values of the vector T = Tν ; i.e., no more than three isothermal sections with adiabatic transitions from one isotherm to another, where ν is the base solution index.2. The duration of each isothermal section is a certain part γ ν of the cycle time, and the adiabatic transitions occurs instantaneously.
In the problem of limiting power, the condition ( 10) is absent, and the average problem has no more than two base points.This means that for any form of dynamic, the optimal cycle of change in the state of the working fluid may consist of only two isotherms and two adiabatic (the decision having one base point does not satisfy the condition ( 9)).
We show that for any heat transfer dynamic that satisfies the "natural" requirements listed in the introduction and the mentioned additional inequality, the problem ( 8)-( 10) has two base solutions, i.e., the corresponding optimal cycle, like the Carnot and Novikov cycles, consists of two isotherms and two adiabats.

Graphic Interpretation of the Solution
Below we give a graphic interpretation of the solution of the provided problems and show under what condition is the optimal solution of the problem of the boundary of the set of realizable regimes, and hence of the maximum efficiency at a given power, consists of two isotherms and two adiabats.If the additional inequality (11) not met then the optimal cycle may contains three isothermal sections, the two of them can be only in the half-cycle of working fluid heating.
We use the index "+", when something is related to the source with a high temperature and an index "−", when something is related to the source with a low temperature.
First, we deduce the system of equations and inequalities that define the boundary of the working region for an arbitrary form of dynamic dependencies.
We divide the set of admissible values of the vector T into two disjoint subsets.The first of which T i = T + > T, and the second is T i = T − < T. The line T i = T does not contain the desired solution, since q = 0 on it.The optimal solution can not be unique, since in this case, for T i = T the heat flux retains its sign and the condition (9) can not be satisfied.
We introduce alongside heat fluxes two functions characterizing fluxes of entropy supplied and removed from the working fluid.Monotonicity conditions ensure that the dependencies T 1 (q + ) and T 2 (q − ) are unambiguous and monotonic; they are zero when T 1 = T + , T 2 = T − .Having calculated the derivatives of curves s + (q + ) and s − (q − ), we obtain that for any heat transfer laws their slopes at the origin are 1/T + and 1/T − , respectively (see Figure 1).This means that in the neighborhood of the origin condition s + (q + ) + s − (q − ) = 0 corresponds to the equality and the efficiency of a machine having arbitrarily low power for any dynamics is The problem of the maximum power in the absence of a restriction (10) on the heat flux is reduced to the construction of a convex hull of the set bounded below by the curves s + (q + ) and s − (q − ).The boundary of the convex hull has only two common points with these curves, in which it is a, tangent.At the point of intersection of this boundary with the abscissa axis, the average entropy flux is zero, and the power (average heat flux value) is maximum.The corresponding value of q = p * .Values of weight coefficients easily find so as their sum of equal is one, and the attitude γ 1 /γ 2 is equal to the ratio of the segments of the convex shell lying above and below the axis of the abscissa.
The average value of the flow is q + fixed the problem of maximum efficiency.If the dependence s + (q + ) convex downward, its convex hull throughout it coincides with the average value of the heat flow corresponding to the point on this curve with ordinate s + ( q + ).If from this point we were to conduct tangent to the curve s − (q − ), then the intersection of this tangent with the abscissa corresponds to the maximum average value of heat flow.The point of this intersection is at the maximum distance from the origin if only the drawn line is a tangent.
Note that the tangent can be drawn in the case when the curve s − (q − ) is not smooth and not convex.In any case, the temperature T 2 corresponding to the touch point is the only one and, consequently, the optimal half-cycle cooling the working fluid is always composed of a single isotherm.
If the average heat flux q + is given and the curve s + (q + ) is not convex, then the point corresponding to the average heat flux can not be on the curve s + (q + ), but on its convex hull.This means that the corresponding flow can be realized by alternately switching the temperatures corresponding to the tangency points of the convex hull with the curve s + (q + ).
Figure 2 shows the case when this dependence is not convex and the given average value of the input heat flux q 0 + corresponds to a point on its convex hull, such that the tangent drawn from it to the curve s − (q − ) intersects with the abscissa axis more to the right than the tangent drawn from the point s + (q 0 + ).The optimal heating half-cycle in this case contains two isotherms, and the cooling half-cycle always contains one.Note that when changing the set average value of q + , the temperature of these isotherms is not changed, but only a fraction of the cycle length that corresponds to them is changed.Under the conditions of the monotonicity, the slope of the curve s + (q + ) = q + T 1 (q + ) almost always increases with increasing q + (the numerator increases linearly, and the denominator decreases).The condition of nonnegativity of the second derivative leads to inequality: that guarantees the convexity of this dependence.

The Optimal Solution
The fact that the averaged problem ( 8)-( 10) has an optimal solution in the form of a cycle containing two isotherms and two adiobates reduces its solution to a finite-dimensional problem with three variables: T 1 , T 2 and γ, where γ is the fraction of the cycle during which the working fluid contacts the hot source.The fraction of contact time with a cold source is 1 − γ.
The value of γ can be expressed through the coordinates of the points on the curves s − (q − ) and s + (q + ) through which the tangent to the function s − (q − ) passes.You can also find the optimal temperature of the working fluid when it is heated and cooled T 1 and T 2 .
The point on the heating curve s 0 + = s + (q 0 + ), q 0 + is given and the dependence T 1 (q + ) determines the value T 0 1 = T 1 (q 0 + ).The tangent equation has the form: Here and further the subscript "0" shows that the derivative is calculated in point q = q 0 − .To determine the coordinates of this point, we have a tangent condition: From this equation find q 0 − , and then s 0 − and T 2 .The point of intersection of the tangent with the abscissa axis determines the maximum power achievable at a spending of heat q 0 + : The condition of constancy of working fluid entropy γs 0 + + (1 − γ)s 0 − = 0 defines time share his contact with the sources: The Relation between the Power and the Entropy Production in the System The production of entropy in the system is caused by the irreversible process of heat exchange between the sources and the working fluid.It is equal to: With increasing flux intensity, the driving forces enclosed in the round parentheses, grow, so the production of entropy with increasing intensity of fluxes grows faster than linear dependence.So, for Newtonian dynamics: Eliminating q − which is q + − p we obtain: from the condition ( 16).This expression is an analogue of the well-known Stodola formula (see [22]), written down for fluxes of heat and work.Maximum power at the fixed value of q + corresponds to the minimum production entropy for any form of dynamic dependencies.
The maximum power for a fixed value of q + corresponds to minimum entropy production for any form of dynamic dependencies.Since the entropy production grows faster than the first term with increasing q + , then the power at q + = 0 is zero, then increases to certain maximum, and then decreases to zero.The machine's efficiency is associated with the production of entropy as: We note that for any heat transfer laws, the production of entropy decreases monotonically with increasing heat transfer coefficients α i .The Lagrange function in the problem of maximum power (minimum dissipation) for a given sum of heat transfer coefficients has the form: The conditions of its stationarity for the sought-for variables leads to the requirement of equality: Solving the problem for the case when the half-cycle of heating can consist of two isotherms, does in two stages.At first solves the problem of constructing convex hull of s + (q + ).Its solution determines the temperature of two isotherms and the proportion of the duration of heating half cycles for each of the dependent given average heat flow q + .The second stage solves the problem of determination cooling isotherm and corresponds to the tangency line drawn from the calculated at the first stage point of intersection of the vertical line with abscissa of q + with the convex hull s + (q + ), the dependence of s − (q − ).At this stage, we determine the shares of the cooling cycle and heating from the overall cycle time.

Characteristic Points of the Boundary
The boundary of the realizability region has three characteristic points (see Figure 3): 1.The origin of coordinates where q + and p tend to zero.The energy and conversion efficiency for any heat exchange law at this point is η C = 1 − T − T + .2. The point at which the power of the machine is maximal.This point is denoted by p * .
The corresponding energy conversion efficiency and heat flux will be denoted by η * and q * + .3. The point at which the irreversibility of the processes is so great that the power of the machine turns out to be zero.The value of the heat flux corresponding to this point is denoted by q max + .
We define the energy conversion efficiency of converter as and find the relationship between the power p, the energy conversion efficiency and the heat consumption q + .From the condition of entropy invariance of the working fluid, it follows that

Boundary Construction
The boundary of the realizability region is determined by the conditions: q + (T + , T 1 ) which at fixed heat exchange coefficients α + , α − determine the values of T 1 (p) and T 2 (p).The latter allows us to find the machine's energy conversion efficiency as For 0 < p < p * the solution of the Equations ( 22) and ( 23), and hence, the value of η(p), is not unique, see Figure 3.The area lying to the left of q * + corresponds to a higher energy conversion efficiency value.For p = p * the energy conversion efficiency value is unique and equal to η * .The uniqueness condition for the real root of p = p * allows us to find p * and η * .For p > p * , the Equations ( 22) and ( 23) have no real roots.
These values can also be found from the solution of the problem under the condition q + (T + , T 1 ) which determines the temperatures of the working fluid for the machine at maximum power.At the root q max + the energy conversion efficiency is zero, therefor T 1 = T 2 = T 0 , see (21).The temperature T 0 is determined from the equation The left side of this equation decreases monotonically, while the right increases monotonically with increasing T 0 , so it has a unique solution, which lies in the range (T − , T + ).For "natural" dynamic dependencies the value of T 0 depends only on the ratio of heat transfer coefficients.

Conditions of Independence Efficiency of the Machine in the Mode of Maximum Power from Coefficients of Heat Exchange with Sources
Let us show in which case the ratio of the working fluid temperatures T 2 T 1 in the solution of the problem ( 25) and ( 26), and hence the efficiency of the heat engine does not depend on the heat transfer coefficients.The Lagrange function for this the problem, has the form The conditions for its stationarity on T 1 and T 2 lead after excluding λ to the equation: Here the sign " " marks the partial derivatives of the dynamic functions with respect to T 1 and T 2 , respectively.After substituting each of the functions in the equality (28) as Thus, the temperature ratio for p = p * , and hence the energy conversion efficiency (21), depends on the heat transfer coefficients only through T 1 , T 2 , which can be found by the solution of the equation system (30) and (26).In (26), as follows from (29), only the ratio of heat transfer coefficients appears.
If the dynamic dependences are such that the right-hand side of the optimality condition (30) does not depend on the temperatures of the working fluid (for Newtonian dynamics it is equal to T − /T + ), then the energy conversion efficiency in the maximum power mode does not depend on the value of the heat transfer coefficients.In the general case, it can depend only on their ratio.So, if the heat transfer areas of the machine change proportionally, then its energy conversion efficiency η * remains unchanged.

Novikov Case
In a number of cases, the problem of calculating the boundaries of the realizability region is facilitated.This applies to the situation when the coefficients fulfill the relation α − >> α + .For thermal or solar power plants, for example, the heat transfer coefficient of condensing steam with the cold source (cooling water) is usually much greater than the heat transfer coefficient of evaporating water with the fuel gases.
In this case, we can assume that the temperature T 2 = T − .The relationship between power and heat consumption then takes the form: and the efficiency is The maximum power condition ∂p ∂q + = 0 in this case leads to the equation for T 1 : The positive root of this equation, lying between T − and T + , determines the energy conversion efficiency in the maximum power mode.Since the ratio of the heat flux and its derivative with respect to the temperature T 1 does not depend on the heat transfer coefficient, the solution of the Equation (32), and therefore, η * does not depend on the heat transfer coefficients for any kind of "natural" dynamic function.

The Boundary of the Realizability Region for Two Types of Heat Transfer Dynamics
In the following we will consider the realizability region for two well-known engines: the Curzon-Ahlborn engine with 1. Newtonian heat transfer and 2. Fourier heat transfer.
The temperatures: T 2 = T 2 (q − , T − ), T 1 = T 1 (q + , T + ).Given these dependencies and conditions (20) you can find the relationship of p with η and q + for specific laws of heat transfer.

Newtonian Heat Transfer Laws
In this case we have Hence, it follows We can conclude Thus, the power and energy conversion efficiency for Newtonian laws of heat transfer are related by equality: Obviously, for the case when η is equal to the Carnot energy conversion efficiency of an reversible heat engine the power p(η k ) = 0, as for η = 0. Let us find the value of η, for which the power of the heat engine is maximal, and the value p * .It is not difficult to show that the function p(η) is convex up, so the condition for its stationarity determines the maximum whence taking into account the fact that p > 0, we obtain which is the well-known Curzon-Ahlborn efficiency.Consequently, the maximum power is As follows from the expression (27), the maximum heat flux is at which the power of the machine becomes zero.

The Linear Phenomenological Law
When we have The linear phenomenological law (Fourier heat transfer), and so , and The maximum possible value of the heat flux is reached at η = 0 and is equal to where ᾱ is again defined as in the Newtonian case.The power vanishes again at η = 0 and at η = η C , for some η = η * it reaches its maximum.Let us find this value by the condition of stationarity of the right-hand side of the Equation (38).We obtain a quadratic equation with respect to η: After dividing the left side of the equation by α − T − and denoting In the case when the heat transfer coefficient at contact with the cold source is much greater than with the hot one, γ goes to 0, and the root of the Equation ( 41) is equal to: In the general case, the machine's energy conversion efficiency in the maximum power mode is equal to the positive real root of Equation (41): In deriving the dependencies between the power and the incoming heat, and between the power and energy conversion efficiency, we used the hypothesis of internal equilibrium of the working fluid, considering the production of entropy in it is equal to zero.If there is entropy production in the system that is not related with heat exchange between the working fluid and the sources, the value of σ will be large, and therefore, for the same heat flux q + values η and the power p decrease.So the dependence p(q + ), constructed in Figure 3, is the boundary of the realizability region for the irreversible conversion of thermal energy into work in a system with two reservoirs.No real heat engine with non-zero power can work in the region lying above the curve p(q + ), and have an energy conversion efficiency greater than that corresponding to the boundary points of this region realizability.The reasonable boundary of real engines is the left branch of this dependence corresponding to the growth of p with the growth of q + .We should note one important fact: the obtained dependencies do not include equations of the state of the working fluid.They require only knowledge of the heat transfer dynamics.

Conclusions
In this article, we show that the cycle heat engine in the maximum power regime for any heat dynamic consists of two isotherms and two adiabatic; we also demostrate the conditions under which its efficiency does not depend on heat transfer coefficient with the source.
Irreversible cycle of a machine with a maximum efficiency can have three isotherms and three adiabats, and two of them can only be in half-cycle heating.
Therefore, all stationary points of the Lagrange function x 0 (λ) are searched for, and the Equation (A1) is solved not for all points of maximum, but for all points of stationarity.In this case, a non-degenerate solution of the actual value of λ exists.
Whether a solution is degenerate is not known in advance; therefore, it is usually assumed that it is not degenerate, and only when there is a doubt, are looking for extreme points of a system of constraints.For a scalar function f , check if it has extremes in which its value is zero.The Lagrange function is written in the more general form considering that the factor λ 0 can take two values: one in the nondegenerate and zero in the degenerate case, for each case define stationary points L, one of which is the desired solution.
Thus, if the point x * is the optimal solution of the conditional maximum problem and the functions f 0 (x) and f (x) are smooth, then there are Lagrange multipliers for which the Lagrange function L is stationary at this point.
Note that at the boundary points of the set D and in the case when the restriction has the form of an inequality of the form f (x) ≥ 0 as for the unconditional maximum problem, the stationarity conditions are replaced with the conditions of local unimprovability of the Lagrange function.
If it were possible a priori to guarantee the existence of a real root of the Equation (A1), in which x * (λ) is a maximum point of L, then one would refuse the requirement of smoothness of functions f 0 , f ,snd find a solution using the Lagrange procedure and in the case when x can take, for example, integer values.This is very tempting, but it is possible only in the transition from the problem of conditional optimization to its averaged formulation.
Averaged optimization and maximum principle.In some problems, the maximum is not searched for the function f 0 (x), but its average value with averaged constraints of the form The set V represents all the values of x that satisfy the constraints imposed on the components of this vector autonomously.It is wider than the set D, which was admissible in problems of conditional maximum.
The solution of the averaged problem (A3) and (A4) of conditional optimization is not the vector x, but the measure function P(x), which can be understood as the distribution density of the values x.If this vector takes only one value x * , then the function P(x) is equal to the Dirac function δ(x − x * ) concentrated at the point x * .
The Dirac function is zero everywhere except for the point at which it is concentrated, the integral of it is equal to unity, and for any continuous function W So, if it is considered a priori that the solution of the averaged problem is sought in the class of δ-functions, then the averaged problem turns into a conditional optimization problem.
In fact, the optimal form of the measure function may be different, so the set of feasible solutions in the averaged problem is wider than in the problem without averaging, and therefore the maximum of the average value of the criterion is not less than the conditional maximum f 0 (x).If it is strictly greater, they say that the transition to the averaged problem is effective/ The simplest averaged problem is the problem of the maximum of the average value of the function f 0 (x) provided that the average value of the argument is given In the problem without averaging, such a statement does not make sense, you just need to substitute x 0 in f 0 (x).
In the averaged problem, for some values of x 0 the average value of the function f 0 (x) may be greater than f 0 (x 0 ).For example, the function f 0 (x) = 1atx = 0 and x = 1, and for other values of the argument is zero.The value of x 0 = 0.5.Choosing P(x) = 0.5(δ(x + δ(x − 1)) in the averaged problem, we obtain that the average value of the function at the point x 0 = 0.5 is equal to not zero, but one.
The dependence of the performance of the energy costs for the majority of the pumps is that high performance can be achieved with less energy if the time the pump operates with a nominal power and the rest of the time off.So the transition to the averaged problem is effective.It is because of these savings build water towers or use other accumulate unit.
The dependence of the maximum average value in the averaged task from a predetermined average value of the argument is called textit convex hull of function.For the function of one variable convex hull coincides with the elastic thread stretched over the top of the graph of the function.At some point it coincides with the graph, while in others it is higher.In the above example, the convex hull coincides with the function only at the points of zero and one, and for every point between the averaged transition to the problem is effective.
Central to the averaged optimization is the Caratheodory theorem: The convex hull of a function f 0 (x) depending on n variables for any value of x 0 can be calculated as the average of at most (n + 1) the ordinates f 0 (x ν ), ν = 0, 1, ..., n.The values of x ν of the vector x are called basic.Thus, for a function of one variable, depending on a given average value of the argument, either one or two basic values in the averaged problem can be.In the first case, the transition to averaging is not effective; in the second, it is effective.
Theorem Caratheodory essentially states that in the problem of the maximum of the average values of n variables for a given average argument value it is possible to use as a measure P(x) more narrower class of functions defined in the form where γ ν is a weight coefficients.Problem (A6) for the case where x is a vector of dimension n given (A5) takes the form We came to a conditional optimization problem, in which the variables are the basic values of the vector x and corresponding to each of these weights.We can write the Lagrangian for this problem and seek a solution among the points that satisfy the conditions of stationarity of this function with respect to x ν and conditions of local unimprovability with respect to γ ν .The latter has a form: This takes into account that the weight factors γ ν ≥ 0, which means their allowable variation is non-negative, so for γ ν > 0 ∂L 0 ∂γ ν , and for γ ν = 0 ∂L 0 ∂γ ν ≤ 0. Returning to more general setting, when the restrictions are imposed not on the mean value of x, but on the average value of one or more functions of x.This formulation can be reduced to a problem (A6), if we introduce into consideration reachability function: Here, f (x) may be of a vector dimension m as vector c.
The following statement is true: The solution of the averaged problem (A3) and (A4) corresponds to the maximum of the average value of the reachability function, provided that the average value of c = 0. Indeed, one way of solving the averaged problem may be as follows: first, for each value of c solve the problem of conditional maximum and memorize x * (c) and f * (c) = f 0 (x * (c)); then solve averaged problem (A8) of the maximum of the average value of the reachability function at zero, i.e., about the value of its convex hull at the origin.
We write the formulation of this problem, taking into account the fact that c ν = f (x ν ) and f * (c ν ) = f 0 (x ν ): (A12) Note that the number of base points is determined by the dimension of the argument of the reachability function, and not the dimension of x.The function f can be a vector even if x is a number, and it can be a scalar function of many variables.
The Lagrange function for the problem (A12) has the form: where L(x, λ) is the Lagrange function for the non-averaged, initial, conditional optimization problem.The local unimprovability conditions (A10) on γ ν lead to the requirement: Therefore, for the base values of x ν (they are included in the solution of the problem with positive weights), there are Lagrange multipliers for which the Lagrange function L of the original problem reaches a maximum in x.According to Caratheodory's theorem, the quantity of such maxima of the Lagrange function can be no more than m + 1.The values of L at these points are the same and equal to the minimum in λ of the maximum L(x, λ) in x.Thus, in the averaged problem, the functions f 0 , f may not be smooth, and the set V may include isolated values of x.

Figure 1 .
Figure 1.The dependence of the entropy flux on the heat flux supplied and removed from the working fluid.

Figure 3 .
Figure 3.The pattern of the dependence of converter power on the cost of heat.