Reversible limit of processes of heat transfer

Limite reversível de processos de transferência de calor

ARTIGOS GERAIS

Reversible limit of processes of heat transfer

Limite reversível de processos de transferência de calor

Jürgen F. Stilck^I,*; Rafael Mynssem Brum^II

^IInstituto de Física e Instituto Nacional de Ciência e Tecnologia de Sistemas Complexos, Universidade Federal Fluminense, Niterói, RJ, Brasil

^IIInstituto de Física, Universidade Federal Fluminense, Niterói, RJ, Brasil

ABSTRACT

We study a process of heat transfer between a body of heat capacity C(T) and a sequence of N heat reservoirs, with temperatures equally spaced between an initial temperature T₀ and a final temperature T_N. The body and the heat reservoirs are isolated from the rest of the universe, and the body is brought in thermal contact successively with reservoirs of increasing temperature. We determine the change of entropy of the composite thermodynamic system in the total process in which the temperature of the body changes from T₀ to T_N. We find that for large values of N the total change of entropy of the composite process is proportional to (T_N - T₀) / N, but eventually a non-monotonic behavior is found at small values of N.

Keywords: second law of thermodynamics, entropy, reversible and irreversible processes, heat exchange.

RESUMO

Estudamos um processo de transferência de calor entre um corpo de capacidade térmica C(T) e uma sequência de N reservatórios térmicos cujas temperaturas são uniformemente espaçadas entre a temperatura inicial T₀ e a temperatura final T_N. O corpo e os reservatórios térmicos são isolados do resto do universo e o corpo é colocado em contato térmico sucessivamente com reservatórios de temperatura crescente. Determinamos a variação de entropia do sistema termodinâmico composto no processo total, em que a temperatura do corpo varia de T₀ a T_N. Encontramos que para valores grandes de N a variação total de entropia é proporcional a (T_N - T₀) / N, mas eventualmente um comportamento não monotônico surge em valores pequenos de N.

Palavras-chave: segunda lei da termodinâmica, entropia, processos reversíveis e irreversíveis, troca de calor.

1. Introduction

One of the original formulations of the second law of thermodynamics, due to Clausius, states that "no process is possible whose only result is the transfer of heat from a body of lower temperature to a body of higher temperature". As discussed by Fermi in his book [1], since as this statement was made the concept of temperature was empirical, to find out which of two bodies has the higher temperature one had to put them in thermal contact and find out in which sense the heat flows, it will flow spontaneously from the body of higher temperature to the other. So, Clausius's statement could be rephrased into "If heat flows spontaneously from body A to body B, no process is possible whose only result is the transfer of heat from body B to body A". To be more precise, let us isolate the two bodies from the rest of the universe, to assure that they will not exchange heat or other forms of energy with the environment, this is implicit in the term 'only result'. We may notice in this statement that already the idea of irreversibility is present: if heat flows spontaneously from A to B, it will not flow spontaneously from B to A. In other words, the process of transfer of heat between two bodies which are isolated from the environment is irreversible.

The second law off thermodynamics says that not all processes which conserve energy, as required by the first law, occur in nature. In many cases (if not all), if a given process in an isolated system occurs, the inverse process does not. This idea was further formalized by the introduction of the concept of entropy, so that processes in isolated systems which lead to a decrease of the entropy are not allowed. In particular, in the process of heat transfer mentioned above, the entropy of the composite isolated system increases if the heat flows from the body of higher temperature to the one of lower temperature, but would decrease in the inverse process, so that only the further is allowed. It is interesting to recall that in his paper of 1885, Clausius summarized the two first laws of thermodynamics in just two sentences: "The energy of the universe is constant. The entropy of the universe tends to a maximum" [2]. We may notice that the universe is the only composite thermodynamic system which is isolated by definition.

The entropy in equilibrium thermodynamics is a state function, so that it is defined for equilibrium states of the system. When two bodies of different temperatures are brought into thermal contact, the composite system in general will go through a sequence of non-equilibrium states, but if we wait long enough, it will finally reach a new equilibrium state, where both bodies of the composite have the same temperature, and in this state the entropy is larger than it was in the initial state. If we now imagine a process which happens at a very slow rate in time when compared to the relaxation time of the system [3], we may reach a situation close to the one in which all intermediate states of the system are equilibrium states. In this limit the process is called quasi-static, and therefore the whole process may be represented as a path in the space of variables which define the equilibrium states of the system. In the particular case in which the entropy of the final state of the quasi-static process is equal to the entropy of the initial state, and therefore to the entropies of all intermediate states, this process will be reversible.

The process we will study in some detail here is the transfer of heat in isolated systems composed by a heat reservoir (constant temperature) and a body whose heat capacity is described by a function of the temperature C(T). The initial temperature of the body is T₀ and its final temperature is T_N. The change of the temperature of the body is accomplished by putting it in thermal contact with a sequence of N heat reservoirs whose temperatures are equally spaced between T₀ and T_N. For a finite number of reservoirs, each of the sub-processes consists of heat transfer between two bodies at different temperatures and is therefore irreversible, leading to an increase of the entropy of the composite isolated system. In the limit N → ∞ , however, each sub-process will contribute with a negligible increase of entropy and, what may not be obvious initially, the total variation of entropy vanishes as well, so that this is a concrete example of a limit which leads to a reversible, quasi-static process of heat transfer.

This problem, of considerable pedagogical interest, has been studied before in the literature. For a constant heat capacity, such as is found for ideal gases, Calkin and Kiang [4] have shown for a cyclic process (T₀→ T_N → T₀), that the total change of entropy vanishes as 1/N in the limit N → ∞ because, although the number of sub-processes increases, the increase of entropy in each of them is proportional to 1/N² in this limit. A similar situation is proposed in the exercise 4.4-6 in Callen's book [3]. The same process was studied also by Thomsen and Bers [5], again for the case of constant heat capacity and discussing in detail that it is always possible, given any positive value e, to choose a value of N for which the increase of entropy is smaller than e. This argument is then repeated for the inverse process, thus characterizing the process in the limit N → ∞ as reversible.

Here we generalize the problem in two ways: the heat capacity of the system is no longer constant and we develop the increase of entropy as a series in 1/N, so that we find the corrections to the leading term. We also show that some symmetries in the result which are valid in the asymptotic limit close to reversibility are not longer present for finite values of N. In particular, we show that, although in the large N limit the increase of entropy decreases monotonically with N, this may not happen at lower numbers of heat reservoirs. Also, the increase of entropy for the direct process (T₀→ T_N) may be different of the same quantity for the inverse process (T_N → T₀). The determination of the change of entropy is given in section 2 and applications to particular cases and final comments may be found in section 3.

2. Determination of the change of entropy

Let us define the problem in more detail. Our composite system consists of the body, whose heat capacity is C(T) and which will undergo a temperature change ΔT = T_N - T₀, and a set of N heat reservoirs, at equally spaced temperatures, so that the temperature of the reservoir j (j = 1, 2, 3, ..., N) is equal to T_j = T₀ + j ΔT/N. The body is brought into thermal contact with the heat reservoirs in order of increasing values of j until its temperature equals the one of the reservoir and therefore is increased bi ΔT/N. In the first sub-process, for example, the body starts with temperature T₀ and ends with temperature T₀ + ΔT/N. After the N'th process, the body will reach the final temperature T_N. This sequence of heat transfer processes is illustrated in Fig. 1. Actually, one could argue that the time required to reach equilibrium for each sub-process of heat transfer would be infinite, and for this reason, Thomsen and Bers, in their discussion of a similar process [5], have chosen the temperature of the j'th reservoir to be T₀ + (j + 1) ΔT/N, thus assuring that the equilibration time of the system with each reservoir is finite. Here, for simplicity, we will not use this improved version of the process.

If a the body, at temperature T, receives a quantity of heat dQ_b, its temperature will change by dT = dQ_b/C(T). The change of entropy of the body in the whole sequence of heat transfer processes is

The change of entropy of all reservoirs in the sequence of processes will be

where in the last passage, we recall that in each sub-process, the body and the reservoir are isolated, so that the heat received by the reservoir is given by dQ_r = -dQ_b = -C(T) dT. We thus may write the total change of entropy as

If it is possible to perform both integrations, this expression above becomes

where S_b(T) = ∫C(T)/T dT and U(T) = ∫C(T) dT. It is convenient to rewrite the expression for the change of entropy as

where κ_j(T, N) = 1-T/T_j and ϕ(T) = C(T)/T. It is now easy to see that, since stability implies that ϕ(T) > 0, we must have ΔS > 0, as required by the second law of thermodynamics. Let us show this in some detail. Suppose that T_N > T₀, so that the temperature of the body increases in the process. It is then clear that κ_j(T, N) > 0 and as a consequence the entropy of the composite system increases. The same happens in the inverse process, where the temperature of the body decreases from T_N to T₀. In this case the change if entropy is

and since now κ_{j - 1}(T, N) < 0 in the range of integration, again the entropy increases. It may be worth noticing that the difference of the entropy increases between the processes with raising and lowering of the temperature of the body δS(N) = ΔS'(N) - ΔS(N) is

and the sign of this expression is not defined in general.

Another general property of the entropy increase is its dependence of the number of reservoirs N. As we will see below, it vanishes in the limit N → ∞ of a quasi-static process. The question we will address is if the change of the entropy decreases monotonically with N. For convenience, we will define the function

where Θ(T) is the step function. The function κ(T, N) is defined in the whole temperature range [T₀, T_N] and has a sawtooth pattern, with N maxima of decreasing values. The expression for the entropy increase Eq. (4) may then be cast into the following form

The difference of entropy increases for two values of N, (N₂ > N₁) is

It is straightforward to see that if N₂ is a multiple of N₁ the difference of κ functions is non-negative, and therefore ΔS(N₁) - ΔS(nN₁) > 0 for any integer n. For example, we have ΔS(1) > ΔS(2). This, however, is no longer true in general if the ratio N₂/N₁ is not an integer. As an example, in Fig. 2 we show the difference of κ functions for N₁ = 5 and N₂ = 6, and we notice that this difference assumes negative values in part of the temperature domain. We thus reach the conclusion that, depending on the function ϕ(T), the entropy increase for N₂ heat reservoirs may be larger than the one for N₁ reservoirs, if 1 < N₂/N₁ < 2.

We now proceed obtaining the general asymptotic behavior of ΔS(N) for N 〉〉 1. Expanding the last integral in the expression for the increase of entropy, Eq. (2), which corresponds to the amount of heat Q_j exchanged by the body with reservoir j, in a Taylor series around the temperature T_j, we get

where a superscript between parenthesis means a derivative of the corresponding function of the order of the superscript. We may now write the total change of entropy as

We now proceed transforming the sum over the reservoirs into an integral, using Euler-MacLaurin's expansion [6]

where B_2k are the Bernoulli numbers (B₂ = 1/6, B₄ = -1/30, B₆ = 1/42, B₈ = -1/30, .... See Ref. [7] for a list of the numbers). The series in the right hand side is divergent in many cases, since the Bernoulli numbers increase very fast for larger values of k. There are procedures to evaluate the rest [6] and we will discuss this issue in the appendix. The sum we will convert into an integral is

and applying the expansion (12) we find

Now we may substitute this last expression into Eq. (11) and collect the terms in the change of entropy of the body in the sequence of processes in powers of 1/N. We notice, as expected, than the term of order 0 cancels, so that, in the limit N → ∞ the whole process becomes reversible. In the appendix, we proceed calculating all coefficients σ_i of the expansion

but here we we will limit ourselves to the leading term, which is

Integrating by parts, we finally obtain the leading term of the expansion

valid in the limit N 〉〉 1.

Since C(T) > 0, we notice that the leading term vanishes only in the trivial case where C(T) = 0 in the whole range of temperatures and there is no heat transfer, so that the asymptotic behavior as 1/N in the increase of entropy, as was found in the particular case of an ideal gas [4], is universal. We also remark that in the asymptotic regime the entropy increase is invariant if the temperatures are interchanged and also monotonically decreasing with N, properties which are not true in general for small values of the number of reservoirs.

3. Discussions and conclusion

We will now apply the results of the preceding section to some particular cases of physical interest. We start with the ideal gas, which has a constant heat capacity C(T) = C₀. Applying Eq. (2) to this particular case, the increase of the entropy will be

where f = T_N/T₀. The leading term in this case is

In Fig. 3 we show some curves of the entropy increase as a function of 1/N, for different values of the ratio of temperatures f, as well as the dashed straight lines which represent the asymptotic behavior. Notice that, as expected, the same asymptotic behavior is found for f = 2 and f = 1/2. If we consider two processes, one with ratio f₁ > 1 and the other with f₂ = 1/f₁ < 1, that is, related by interchange of the temperatures, both will have the same asymptotic behavior, as may be seen in the particular cases depicted in the figure. Also, the increase of entropy in the second process, where the body is cooled, is always larger than the one in the first process. Another way to state the same property is that, as can be seen in the curves, the increase of entropy is a concave function of 1/N for processes where the body is heated and it is a convex function when the temperature of the body decreases. In the asymptotic regime these properties may be obtained the next term of the Euler-MacLaurin expansion, as will be discussed in the appendix. Also, the increase of entropy for this particular case is a monotonically decreasing function of N.

In most cases the variation of entropy decreases monotonically with N, a property which seems also intuitive, since by decreasing the temperature steps in a certain sense the composite process gets closer to a reversible process. For the ideal gas this is always the case. We have seen above, however, that this may not be generally true. If we look at the Eq. (9) for the difference of entropy increases for different values of N, a non-monotonic result is possible if the function ϕ(T) has one or more peaks at the temperatures where κ(T, N₁) - κ(T, N₂) is negative. One system which could possibly show a non-monotonic behavior is a two-state system [3], since its heat capacity displays a peak, sometimes called "Schottky hump". However, we found that this is not the case, the reason is that the peak is too wide to produce a non-monotonic behavior.

We turn our attention to a situation where the width of the peak in the heat capacity may be changed by adjusting a parameter. The simple expression we will use for the heat capacity is

where κ = T/T_c, so that T_c is the temperature where the maximum of the heat capacity is located, and C is a constant with the dimension of entropy. Although this heat capacity does not correspond to a particular physical system, it is convenient for analytic calculations and consistent with the third law of thermodynamics. Of course, at least for odd values of the parameter m, it is valid only if κ is in the range [0, 2], so we will restrict ourselves to this range of temperatures. As the parameter m grows, the peak in C(T) narrows, and our main interest in this system is that, for sufficiently large values of m, non-monotonic behavior is found in ΔS(N).

In Fig. 4 we show, in the same graph, the function κ(κ, 2) - κ(κ, 3) and the function ϕ(κ) = C(κ)/(Cκ), given by Eq. (20) with m = 16. As is apparent in the graph, with the choices of initial (κ₀ = 0.286) and final (κ_N = 2.0) temperatures we made, the maximum of the peak in C(T) is close to the center of the negative peak in κ(κ, 2) - κ(κ, 3), so that the difference ΔS(2) - ΔS(3) is minimized, since this difference is given by Eq. (9). In Table 1 we list the results for [ΔS(3) - ΔS(2)]/C for different values of m, with the choices above for κ₀ and κ_N, and it is clear that a non-monotonic behavior is found for m > 16. Finally, we found some evidence that a local maximum in the entropy could also be found at N = 4, but if this actually happens for this particular form of heat capacity, it will be at rather large values of m, which lead to numerical problems. We remark that, in order to reduce numerical errors, in all examples we discuss the integrations involved in the calculation of the entropy increase Eq. (2) are performed analytically, as is done in Eq. (3).

As discussed before, we expect, qualitatively, the non-monotonic behavior to be enhanced if the peak in the heat capacity is more pronounced. This leads us, in this final example, to a system which undergoes a continuous phase transition, since in this case a singularity at the critical temperature T_c is found in the heat capacity, of the form C(T) ≈ |1 - T/T_c|^-α, described by the critical exponent α [8]. To the singular behavior of the heat capacity, usually regular contributions should be added, so that we will use the simple expression

where we define the reduced dimensionless temperature κ = T/T_c and which is compatible with the third law of thermodynamics. It should be mentioned that since the exchange of heat when the temperature of the body crosses the critical temperature has to be finite, we should restrict the critical exponent a to values smaller than 1. It is easy to perform the necessary integrations in this case, so that the entropy increase will be given by Eq. (3) with

and

We would expect that the non-monotonic behavior of the entropy should be enhanced as the value of the critical exponent a increases. This is actually true, but at variance with what is seen in the preceding example of a peak in the heat capacity, the local maxima of the entropy increase, for lower values of a, appear at larger values of N. In Fig. 5 we see results of the increase of entropy for two values of the exponent, α = 0.7 and α = 0.5. It is apparent that, for α = 0.7 we already have an increase of the change of entropy between N₁ = 2 and N₂ = 3, the segment with a thicker line in the curve. Many other pairs of successive results with the same property are found for larger values of N. For α = 0.5, the first pair where the change of entropy increases is seen at N₁ = 7, N₂ = 8.

For a given a value of the critical exponent α, we may find N₁, the lowest value of N such that ΔS(N₂) > ΔS(N₁), where N₂ = N₁ + 1. In Fig. 6, the exponent a is plotted as a function of 2/(N₁ + N₂), and we notice that, at least for large values of α, that the function shows steps. The results shown are calculated for κ₀ = 0.2 and κ_N = 2.2. Thus, for example, if α > 0.555, we notice that the first increase of ΔS(N) is seen between N₁ = 2 and N₂ = 3, while if α is in the range [0.555,0.454] this increase is seen at N₁ = 7, N₂ = 8, if a is further decreased the increase shifts to N₁ = 12, N₂ = 13, and so on. This pattern, although qualitatively remaining the same, changes quantitatively for a different temperature range. We notice an interesting behavior for small values of α, which leads to the question if, for any positive α an increase if ΔS(N) is found for α finite value of N. Due to round-off errors, we could not address this question numerically. However, since the coefficients σ_i in the Euler-MacLaurin expansion Eq. (15) diverge in this case for i > 1 and any positive α, the answer to this question should be positive. We discuss this point somewhat more in the appendix.

Finally, we notice that, in opposition to what was found in the first two examples we studied, when the body undergoes a continuous phase transition it is possible that ΔS(N) > ΔS'(N), that is, the increase of entropy in the heating process is larger than the one in the cooling process, a possibility which was anticipated in the discussion after Eq. (6). In Fig. 7 we present results for the difference in entropy increases ΔS'(N) - ΔS(N), for α = 0.7 and α = 0.5, as functions of 1/N. The temperature interval is κ₀ = 0.2, κ_N = 2.2 Although in most cases this difference is positive, for α = 0.7 already for n = 13 we get a negative result, and a set of values of N above this one have the same property. For α = 0.5 the effect is smaller, and the first occurrence of a negative result is for N = 88. Again it would be interesting to study more details of these results at large values of N, but round-off errors prevent this to be done numerically.

A final remark on this example is that, if κ 〈〈 1, C(κ) ≈ Cκ, so that it will show a behavior similar to a gas of non-interacting fermions at low temperature. In this case, if the initial temperature κ₀ approaches 0, the coefficient σ₁ diverges as a logarithm. In other words, in a heating process starting at κ₀ = 0, the tangent to the curve ΔS ×1/N at the origin is vertical. This curve is concave, since σ₂ is negative. The even more unphysical cooling process with vanishing final temperature shows a very uncommon behavior: ΔS(N) diverges for any finite N and vanishes as N → ∞ .

In conclusion, we notice that in general the process of heat transfer between a body and a sequence of heat reservoirs, is a very rich physical situation for understanding of the second law of thermodynamics and the concept of a reversible process. In particular, if the body undergoes a continuous phase transition in the composite process, some curious phenomena appear , which at first seem to defy common sense, but of course none of them violates thermodynamics.

Appendix: Euler-MacLaurin expansion

Here we will develop the terms of the Euler-MacLaurin expansion in some detail. If we substitute Eq. (14) into Eq. (11) and collect the terms in the change of entropy of the body in the sequence of processes in powers of ΔT/N. We notice, as expected, than the term of order 0 cancels, so that, in the limit N → ∞ the whole process becomes reversible. The coefficients we obtain for the change of entropy Eq. (15) are

The two first coefficients of this series, after integrating by parts, are:, are

and

We notice that σ₁> 0, and vanishes only if the temperature interval is zero or if C(T) = 0 in the whole interval. In other words, if there is a non-zero heat exchange between the body and the reservoirs, σ₁ is positive. If this coefficient could be negative, the second law would be violated. It is also invariant with respect to a permutation of the temperatures, as was already noted in the examples discussed above. The sign of the second coefficient σ₂ is not fixed, and, as happens with all coefficients of even order, it switches under permutation of the temperatures. The odd order coefficients do not change their sign under this permutation.

In the particular case of an ideal gas, where C(T) = C₀, the coefficients are

where the last expression is valid for k = 1, 2, 3, .... We notice that, with the exception of the dominant term, only terms with even powers of 1/N appear in the expansion for this particular case. For the heating process T_N > T₀, σ₂ < 0, so that the ΔS is a concave function of 1/N in the limit N → ∞ . For the cooling process T_N < T₀, the function is convex. These properties can be verified in the numerical results presented in Fig. 3.

Very often the expansion in the Euler-MacLaurin summation (12) does not converge. We find that even for N = 1 the truncated asymptotic values are close to the exact ones at rather low orders. However, truncating at higher orders does not necessarily lead to results closer to the exact ones. For example, for N = 1 the result closest to the exact one is found truncating the series at k_max = 6. We notice that higher orders of truncation lead to wrong results, as may be seen in Fig. 8, where the increase of entropy calculated using the expansion Eq. (15), for N = 1 and T_N/T₀ = 2, is shown for several orders of truncation. The results oscillate around the exact value at successive orders. The best result is obtained for k_max= 6 and if we truncate the expansion at larger orders we observe that they start deviating from the exact value.

In the Table 2, we present more results of such calculations. We notice that quite accurate estimates of the entropy change may be obtained if the series is truncated at the most favorable order. The question of determining the order which minimizes the rest R in the Euler-MacLaurin expansion Eq. (12) is rather technical and we refer to the specialized literature for details, but in general it may be assured that, if f(x) in Eq. (12) has always the same sign and if the function and all its derivatives tend monotonically to 0 as x → ∞ , which are valid for f(N) = ΔS(N), then the rest is of the same order and has the same sign of the first neglected term [6].

Finally, it is clear from Eq. (24) for the expansion coefficients, that for the case of a singular heat capacity at a critical temperature T_c (C(T) = Cκ(1 - κ)^α, with κ = T/T_c), all but the first coefficients diverge if the temperature interval includes the critical temperature. This fact may account for the rather peculiar behavior found numerically in the large N limit.

References

[1] E. Fermi, Thermodynamics, (Dover, New York, 1956).

[2] R. Clausius, Annalen der Physik und Chemie, 125, 353 (1865), a translation may be found in William Francis Magie, A Source Book in Physics, (McGraw-Hill, New York, 1935).

[3] H.B. Callen, Thermodynamics and an Introduction to Thermostatistics, (John Wiley & Sons, New York, 1985).

[4] M.G. Calkin and D. Kiang, Am. J. Phys. 51, 78 (1983)

[5] J.S. Thomsen and H.C. Bers, Am. J. Phys. 64, 580 (1996).

[6] K. Knopp, Theorie und Anwendung der unendlichen Reihen, (Springer Verlag, Berlin, 1964).

[7] The denominators and numerators of several Bernoulli numbers may be found in the addresses http://oeis.org/A027642 and http://oeis.org/A000367. respectively. These data are provided by the On Line Encyclopedia of Integer Sequences (http://oeis.org/wiki/Welcome).

[8] J. Yeomans, Statistical Mechanics of Phase Transitions, (Oxford University Press, Oxford, 1992).

Recebido em 23/2/2013;

Aceito em 2/3/2013;

Publicado em 30/10/2013

Publication Dates

• Publication in this collection
  11 Feb 2014
• Date of issue
  Dec 2013

History

• Received
  23 Feb 2013
• Accepted
  02 Mar 2013