Within the context of finite-time thermodynamics (FTTs) some models of convective atmospheric cells have been proposed to calculate the efficiency of the conversion of solar energy into wind energy and also for calculating the surface temperature of the planets of the solar system. One of these models is the Gordon and Zarmi (GZ) model, which consists in taking the sun-earth-wind system as a FTT-cyclic heat engine where the heat input is solar radiation, the working fluid is the earth's atmosphere and the energy in the winds is the work produced. The cold reservoir to which the engine rejects heat is the 3 K surrounding universe. In the present work we apply the GZ-model to investigate some features of the convective zone of the sun by means of a possible structure of successive convective cells along the well-established convective region of the sun. That is, from 0.714 RS up to RS being RS the radius of the sun. Besides, we estimate the number of cells of the model, the possible size of the cells, their thermal efficiency, and also their average power output. Our calculations were made by means of two FTT regimes of performance: the maximum power regime and the maximum ecological function regime. Our results are in reasonable agreement with others reported in the literature.
1. Introduction
The problem of thermal balance between the planets of the solar system and the sun under a finite-time thermodynamics approach has been treated by several authors [1–8]. In some of these articles the question of the conversion of solar energy into wind energy is also treated. In particular, De Vos [3] demonstrated that cosmic radiation, starlight, and moonlight can be neglected for the thermal balance of any of the planets of the solar system and only the following quantities have an influence: the incident solar influx or solar constant Isc, the planet’s albedo ρ, and the greenhouse effect of the planet’s atmosphere crudely evaluated by means of a coefficient γ. This coefficient can be taken as the normalized greenhouse effect introduced by Raval and Ramanathan in [9]. When only the global thermal balance between the sun and a planet is considered, one can roughly obtain the planet’s surface temperature assumed as a uniform temperature Tp. If the conversion of solar energy into wind energy is to be modeled, it is necessary to involve at least two representative atmospheric temperatures for making the creation of work possible; that is, to take the planet’s atmosphere as a working fluid that converts heat into mechanical work. In 1989, Gordon and Zarmi [1] introduced a FTT-model taking the sun-earth-wind system as a FTT-cyclic heat engine where the heat input is solar radiation, the working fluid is the earth’s atmosphere, and the energy in the winds is the work produced; the cold reservoir to which the engine rejects heat is the 3K surrounding universe. By means of this simplified model, Gordon and Zarmi were able to obtain reasonable values for the annual average power in the earth’s winds and for the average maximum and minimum temperatures of the atmosphere, without resorting to detailed dynamic models of the earth’s atmosphere, and without considering any other effect (such as earth’s rotation, earth’s orbital motion around the Sun, and ocean currents). Later, De Vos and Flater [2] extended the GZ model to take into account the wind energy dissipation by means of a maximum power criterion. This model was extended by De Vos and van der Wel [4, 5] by constructing a model based in convective Hadley cells. All the models used in [1–5] are endoreversible ones in the sense of FTT [10], that is, all irreversibilities are located in the exchanges between the engine and the external world. The GZ model was later studied under a nonendoreversible approach and by using the so-called ecological optimization criterion [6, 7]. This approach [11] consists of maximizing a function E that represents a good compromise between high-power output and low-entropy production. The function E is given byE=P-TextΔSu,
where P is the power output of the cycle, ΔSu the total entropy production (system plus surroundings) per cycle, and Text is the temperature of the cold reservoir. This optimization criterion for the case of the so-called Curzon-Ahlborn cycle [12], for instance, leads to a cycle configuration such that for maximum E it produces around 75% of the maximum power and only about 25% of the entropy produced in the maximum power regime [13]. By means of employing this criterion in a nonendoreversible GZ model, the authors of [6] also found reasonable values for the annual average power of the winds and for the extreme temperatures of the earth’s troposphere. Later, the non-endoreversible GZ model was applied to calculate the surface temperature of planets of the solar system [8], considering two regimes of performance: maximum power regime and maximum ecological function regime. In this work, we apply the GZ model to the convective zone of the sun which is located between 0.714 RS and RS [14]. Our FTT approach leads to a possible structure of the convective region of the sun consisting in approximately sixteen coupled cells. It is important to remark that these sixteen convective Carnotian cells are only a kind of idealized cells, thermodynamically equivalent to the complex structure of the actual convective zone of the sun. The paper is organized as follows: in Section 2, we present a brief review of the GZ model for the convective cells under both the maximum power and the ecological function regimes. In Section 3, we applied the GZ model to the convective zone of the sun and finally in Section 4 we present some concluding remarks.
2. Endoreversible GZ Model for Atmospheric Convection
The endoreversible GZ model is based on annual average quantities and thus it does not represent actual convective cells but a kind of annual virtual cell that takes into account the global thermodynamic restrictions over the convection as a dominant energy transfer mechanism in the air (which has a large Rayleigh number). Besides, this kind of model must only be taken as one that producing better upper bounds than those calculated by means of classical equilibrium thermodynamics, which is one of the main purposes of FTT.
2.1. Maximum Power Regime
In Figure 1, a schematic view of a simplified sun-earth-winds system as a heat engine cycle is depicted. This cycle consists of four branches: (1) two isothermal branches, one in which the atmosphere absorbs solar radiation at low altitudes and one in which the atmosphere rejects heat at high altitudes to the universe and (2) two intermediate instantaneous adiabats [10] with rising and falling currents. In [15], it was shown that a Curzon-Ahlborn FTT cycle in the endoreversible limit with instantaneous adiabats is reached for large compression ratios. In the GZ virtual cells, it is feasible to consider that this condition is fulfilled. According to GZ, this oversimplified Carnot-like engine corresponds very approximately to the global scale motion of wind in convective cells. Below, we use all of GZ model’s assumptions.
Scheme of a simplified solar-driven heat engine (taken from [1]).
For example, the work performed by the working fluid in one cycle W, the internal energy of the working fluid U, and the yearly average solar radiation flux qs are expressed per unit area of the earth’s surface. The temperatures of the four-branch cycle are taken as follows: T1 is the working fluid temperature in the isothermal branch at the lowest altitude, where the working fluid absorbs solar radiation for half of the cycle. During the second half of the cycle, heat is rejected via black-body radiation from the working fluid at temperature T2 (highest altitude of the cell) to the cold reservoir at temperature Text (the surrounding 3K universe). In the GZ model, the objective is to maximize the work per cycle (average power) subjected to the endoreversibility constraint [10], that is,ΔSint=∫0to{qs(t)-σ[T4(t)-Text4(t)]T(t)}dt=0,
where ΔSint is the change of entropy per unit area, t0 is the time of one cycle, σ is the Stefan-Boltzmann constant (5.67×10-8W/m2K4), and qs, and T are functions of time t, taken as [1]
T(t)={T1;if0≤t≤t02,T2;ift02≤t≤t0,qs(t)={0;ift02≤t≤t0,Isc(1-ρ)2;if0≤t≤t02,
in the same way, Text=3K for 0≤t≤t0, with Isc the yearly average solar constant (1373 W/m2) and ρ=0.35 [2], the effective average albedo of the earth’s atmosphere. The GZ model maximizes the work per cycle W, taken from the first law of thermodynamics:ΔU=-W+∫0t0{qs(t)-σ[T4(t)-Text4(t)]}dt=0,
by denoting average values as,T¯=T1+T22,Tn¯=T1n+T2n2,qs¯=Isc(1-ρ)4,
where n is an integer with values n=3 or 4. The factor of 1/4 arises from a factor of 1/2 to account for the day/night difference and a geometric factor of 1/2 to account for the earth’s cross section, which is intercepted by solar radiation, as opposed to the corresponding hemispherical surface area of the earth. From (4) and (5) and taking into account the constraint given by (2), GZ construct the following Lagrangian L;
L=T4(t)+λ[qs(t)T(t)-σT3(t)],
where λ is a Lagrange multiplier. The Euler-Lagrange formalism will be used, by using ∂L(t)/∂T(t)=0, GZ found the following values for the earth’s atmosphere: T1=277K, T2=192K, and Pmax=Wmax/t0=17.1W/m2. These numerical values are not so far from “actual” values, which are P≈7W/m2 [16], T1=290K (at ground level), and T2≈195K (at an altitude of around 75–90 Km). However, as GZ assert, their power calculation must be taken as an upper bound due to several idealizations in their model. In [6], another endoreversible case was analyzed in which the tropopause layer with Text=200K was used as cold reservoir. In this case, the following Lagrangian was used:L(t)=qs+σText4¯-σT4¯-α[qs¯T1-σ(T13+T23)2-σText4(1T1+1T2)],
with α a Lagrange multiplier. By numerically solving ∂L(t)/∂T(t)=0, they obtained T1=293.387K and T2=239.267K, which are excellent values for convective cells restricted to the troposphere. If these temperature values are substituted in the expression for the average power (see [6])P=qs+σText4¯-σT4¯,
a value of P=10.758W/m2 is obtained, which is a good value for the wind power [16].
2.2. Ecological Function Regime
As De Vos and Flater [2] state, no mechanism guarantees that the atmosphere maximizes the wind power. In fact, some authors [17–19] have recognized that the earth’s atmosphere operates at nearly its maximum efficiency; thus, from an FTT point of view, an ecological-type criterion seems feasible. This is due to the properties of the E function, which at its maximum value represents an austere compromise between power and entropy production, additionally leading to a high efficiency [11, 13]. This ecological criterion, as previously occurred with the concepts of power output and efficiency [20], has also been used in the context of irreversible thermodynamics [21–23]. In particular, in [7] the so-called ecological criterion was applied to the GZ model. This criterion consists in maximizing equation (1). By means of the second law of thermodynamics, first, we calculate ΔSu, the total entropy change per cycle (system plus surroundings),ΔSu=∫0t0{-qs(t)+σ[T4(t)-Text4(t)]T(t)}dt.
From (3), we obtainΔSu=∫0t0/2{-qs(t)T1+σ(T13-Text4T1)}dt-∫t0/2t0{σ(T24-Text4Text)}dt.
Thus, the total entropy production is given by [7, 8],Σ=ΔSut0≈qs¯T1+σ2(T13+T24Text),
here, we have used the approximation qs¯≫σText4(223W/m2≫4.59×10-6W/m2) with Text=3K. So, the ecological function E for this case isE=qs¯-σT4¯+Textqs¯T1-σText2(T13+T24Text).
By using (12) and the constraint given by (2), we proposed the following Lagrangian function LE:LE=qs¯-σT4¯+Textqs¯T1-σText2×(T13+T24Text)-α[qs¯T1-σT3¯],
with α being the Lagrange multiplier. By substituting the values of qs¯, σ, and Text and numerically solving ∂L(t)/∂T(t)=0, we find T1=294.08K,T2=109.54K and P=6.89W/m2, which are reasonable values for T1 and P, but not for T2. However, if we use as a cold reservoir, the tropopause layer with Text=200K, we can now use the Lagrangian function: [24],LE=qs+σText4¯-σT4¯+(qs¯+σText42)TextT1-σText2(T13+T24Text)-σText42-β[qs¯T1-σ(T13+T23)2+σText4(1T1+1T2)],
with β a Lagrange multiplier. By using again the Euler-Lagrange formalism, we numerically obtain T1=303K,T2=219K, and P=7W/m2 which are very good values, for T1, T2, and P. Besides, these values are restricted to typical values in the troposphere, where the climatic phenomena occurs. It is important to note that the power values (6.89W/m2 and 7W/m2), which were calculated by the means of the ecological function, were deduced without considering the greenhouse effect (γ coefficient). When the later is taken into account, the values of P are bigger than 7W/m2 [7, 8]. These scenarios lead to larger upper bounds for the wind’s power permitting an energy excess for other relevant dissipative processes such as ocean currents and biological structuring.
3. The GZ Model Applied to the Convective Zone of the Sun
The core of the sun goes from 0 to 0.2 RS, where RS (6.96×108m[14]) is the radius of the sun. The radiative zone embraces the region between 0.2 RS and 0.714 RS and beyond that lies the convective zone. The later is estimated to have a width of approximately 0.286 RS [14]. In (8) and (11) the input data were qs and Text, the thermal energy and the temperature of the surrounding cold thermal bath for the earth’s atmospheric cells, respectively. In the case of the convective zone of the sun, first we will use the maximum power criterion. In Figure 2, we show the heat fluxes balance for the convective zone of the sun. Then, by using (2), (3), (4), and (5) we obtain the following Lagrangian functional:
L(T1,T2,λ)=qs2+σText42-σ2(T14+T24)-λ[qs2T1-σ2(T13+T23)+σText42(1T1+1T2)],
where λ is a Lagrange multiplier, Text=3K, T1=2.18×106K [14] is the temperature of the spherical layer at 0.714 Rs and qs=σT14 the input thermal energy at the lower layer of the convective zone. The energy transport through the sun can be considered as a “sandwich”, that is, there are two regions in which radiation transports the energy separated by a region where convection transports it [25]. Strictly speaking, qs should be calculated by means of a diffusive model based on kinetic theory of gases [25]. However, for simplicity, in our thermodynamic model we take the 0.714 Rs layer at T1≈2×106K as a blackbody radiant system (qs=σT14, see Figure 2). The radiation emitted by this layer is rapidly absorbed by the gases at the bottom of the convective zone. For the definition of T2, see Figure 2. By using the Euler-Lagrange formalism over the Lagrangian of (15), that is, ∂L(t)/∂T(t)=0, we obtain the following equations:T15-λσ(qs4+3σ4T14+σText44)=0,T25-λ(34T24+Text44)=0.
By eliminating λ from these equations and by using the restriction given by (2), we obtain
T24(4T13T2+T24-3T14)-Text4(T14+T24)=0.
In this equation, the only unknown variable is T2. Then, we numerically solve (17) to obtain T2, the temperature of the upper bound for the first convective cell starting from T1=2.187761×106K [14]. Our next step is to take the obtained T2 value of the first cell as the temperature of the lower layer of the following successive cell. This new T2 value is taken as T1 in (17) and then we calculate a new T2 for the second cell. For the following successive cells we use the same recursive procedure until to reach a final T2 coinciding approximately with the well-known value of the average surface temperature of the sun, which is TS≈5780K [14].
Schematic diagram of the energy fluxes present in the first internal convective cell. T1 at 0.714 Rs is taken as the temperature of the first isothermal layer, Text=3K is taken as the cold reservoir temperature, and T2 is taken as the upper shell temperature of the first convective cell. Short arrows indicate that emitted radiation is rapidly absorbed by opaque gases.
In Table 1 we show that after 16 successive Carnotian convective cells we reach a final T2≈6000K. Table 1 shows our results for Carnotian cells performing in the maximum power regime. As one can see in this Table (third column), the widths of the cells are decreasing toward the outer regions. The total width is around 0.280 Rs which is not so far of the value 0.286 Rs given by other sun models [26]. If we take as the mode of thermodynamic performance of the sun’s convective cells the so-called maximum ecological regime [11], in a similar way as (15), then we obtain the following Lagrangian functional:LE(T1,T2,λ)=qs2(1-TextT1)+σText4-σ(T14+T24)+σText2(T13+T23)-σText2(T14T2+T24T1)-σText52(1T1+1T2)-λ[qs2T1-σ2(T13+T23)-σText42(1T1+1T2)].
By using the Euler-Lagrange formalism over the Lagrangian of (18), that is, ∂LE(t)/∂T(t)=0 and following a similar procedure as in the case of (17), we obtain8Text4T1T2(T14+T24)+3Text5(T15+T25)+TextT2(4T18+13T15T23-16T13T25-7T28)-T1T25(32T13T2+8T24-24T14)=0.
Similarly to (15), the only unknown variable in this equation is T2. Following a similar numerical procedure as in the case of maximum power conditions, we can calculate a convective cell structure. In Table 2 we present the numerical results for the maximum ecological function. We can see in Table 2 that with 16 successive Carnotian convective cells we can reach a final T2≈6000K.
Maximum power regime case: First column shows the normalized radial position of the hot layers corresponding to the sixteen virtual convective cells. The following columns give, respectively, second, the cell’s widths; third, the hot isotherms; fourth, the cold isotherms; fifth, the average power output; sixth, thermal efficiency.
No.
r(T1)/Rs
Δr=r(T2)-r(T1)(Km)
T1×106(K)
T2×106(K)
W¯(erg/cm2·s)
η=η(T1,T2)
1
0.714
61370.9
2.187761
1.51504
6.49519×1017
0.307495
2
0.802177
42499.6
1.51504
1.04917
4.65789×1017
0.307495
3
0.863239
29431.2
1.04917
0.726555
1.07122×1017
0.307495
4
0.905525
20381.3
0.726555
0.503143
2.46361×1016
0.307495
5
0.934809
14114.1
0.503143
0.348429
5.66582×1015
0.307495
6
0.955088
9774.09
0.348429
0.241289
1.30303×1015
0.307495
7
0.969131
6768.61
0.241289
0.167094
2.99671×1014
0.307495
8
0.978856
4687.29
0.167094
0.115713
6.89186×1013
0.307495
9
0.985591
3245.97
0.115713
0.0801319
1.585×1013
0.307495
10
0.990254
2247.85
0.0801319
0.0554917
3.64521×1012
0.307495
11
0.993484
1556.65
0.0554917
0.0384283
8.38353×1011
0.307495
12
0.995721
1077.99
0.0384283
0.0266118
1.92829×1011
0.307495
13
0.99727
746.511
0.0266118
0.0184288
4.43714×1010
0.307495
14
0.998342
516.963
0.0184288
0.012762
1.02289×1010
0.307495
15
0.999085
357.999
0.012762
0.00883776
2.37682×109
0.307495
16
0.999599
247.916
0.00883776
0.00612019
5.7099×108
0.307495
Maximum ecological regime case: first column shows the normalized radial position of the hot layers corresponding to the sixteen virtual convective cells. The following columns give, respectively, second, the cell’s widths; third, the hot isotherms; fourth, the cold isotherms; fifth, the average power output; sixth, thermal efficiency.
No.
r(T1)/Rs
Δr=r(T2)-r(T1)(Km)
T1×106(K)
T2×106(K)
W¯(erg/cm2·s)
η=η(T1,T2)
1
0.714
61370.8
2.187761
1.51504
6.49519×1017
0.307494
2
0.802176
42499.5
1.51504
1.04917
4.65788×1017
0.307494
3
0.863239
29431.1
1.04917
0.726558
1.07123×1017
0.307494
4
0.905525
20381.2
0.726558
0.503147
2.46362×1016
0.307493
5
0.934808
14114.1
0.503147
0.348433
5.66589×1015
0.307492
6
0.955087
9774.07
0.348433
0.241293
1.30306×1015
0.307491
7
0.96913
6768.59
0.241293
0.167098
2.99682×1014
0.307489
8
0.978855
4687.28
0.167098
0.115718
6.89224×1013
0.307486
9
0.98559
3245.97
0.115718
0.0801366
1.58513×1013
0.307482
10
0.990254
2247.85
0.0801366
0.0554965
3.64566×1012
0.307476
11
0.993483
1556.65
0.0554965
0.0384331
8.38504×1011
0.307468
12
0.99572
1077.98
0.0384331
0.0266166
1.9288×1011
0.307456
13
0.997269
746.509
0.0266166
0.0184336
4.43885×1010
0.307439
14
0.998341
516.961
0.0184336
0.0127669
1.02346×1010
0.307414
15
0.999084
357.998
0.0127669
0.00884262
2.37873×109
0.307378
16
0.999599
247.916
0.00884262
0.00612506
5.71627×108
0.307326
Our results in Table 2 again show that the width of the cells decrease with increasing radius. The total width in this case is around 0.2859 Rs which is practically the value 0.286 Rs given by other sun models [26].
A remarkable fact observed in Tables 1 and 2 (third column) is that between the cell number 10 and 16, the vertical linear sizes are between 2247 Km and 247 Km, respectively. These are values near to those reported for the linear sizes of granules in [25], which are typically around 900–1000 Km, reaching their largest values up to 2000 Km in diameter. On the other hand, in the highest convective cell of our model, the average power has a value of 5.7×109erg/cm2s, which is of the order of the power reported in [25] for convection in the photosphere (which is 7×109erg/cm2s). Our highest cell overlaps with photosphere. This result is also of the order of the power reported for a mixing length theory of convection in [25], which is 10–20×109erg/cm2s. Clearly, our oversimplified model coincides with those reported in [25] in that the energy transported by convection must increase rapidly as we go below the surface region of the convective zone. Finally, it is very interesting that all 16 cells in Tables 1 and 2 have practically the same thermal efficiency, η≈0.307.
4. Concluding Remarks
In the present work we have used a simplified finite-time thermodynamic method to describe the global thermal properties of the convective zone of the sun. This method was previously used by Gordon and Zarmi to describe convective motions of the air in the earth’s atmosphere. These authors assert that this FTT-approach corresponds very approximately to the global scale motion of the wind in convective cells. However, it is necessary to remark that convective cells of this kind of FTT-models are only virtual cells performing by unit area and yearly averages. Thus, they only represent the global thermodynamic properties stemming from the first and second laws of thermodynamics; that is, kind of thermodynamically equivalent cells that only captures global average quantities and discards any other dynamical detail. Nevertheless, all these simplifications permit to obtain reasonable values for some thermal quantities associated to the convective zone of the sun. Our simplification is mainly based in taking several spherical virtual layers as black-body radiant surfaces, whose emitted radiation is rapidly absorbed by the opaque gases of the convective zone. This radiant energy is taken as the driver energy of convective cells.
Acknowledgments
This work was supported in part by CONACYT, COFAA, and EDI-IPN-México.
GordonM.ZarmiY.Wind energy as a solar-driven heat engine: a thermodynamic approach198957995998De VosA.FlaterG.The maximum efficiency of the conversion of solar energy into wind energy199159751754De Vos A.1992Oxford, UKOxford University PressDe VosA.van del WelP.Endoreversible models of the conversion of solar energy into wind energy1992177789De VosA.van der WelP.The efficiency of the conversion of solar energy into wind energy by means of Hadley cells199346419320210.1007/BF00865706Barranco-JiménezM. A.Angulo-BrownF.A nonendoreversible model for wind energy as a solar-driven heat engine199680948724876Barranco-JiménezM. A.Angulo-BrownF.A simple model on the influence of the greenhouse effect on the efficiency of solar-to-wind energy conversion2003265235Barranco-JiménezM. A.Chimal-EguíaJ. C.Angulo-BrownF.The Gordon and Zarmi model for convective atmospheric cells under the ecological criterion applied to the planets of the solar system20065232052122-s2.0-33747451300RavalA.RamanathanV.Observational determination of the greenhouse effect198934262517587612-s2.0-002485586110.1038/342758a0RubinM. H.Optimal configuration of a class of irreversible heat engines—I19791931272127610.1103/PhysRevA.19.1272Angulo-BrownF.An ecological optimization criterion for finite-time heat engines199169117465746910.1063/1.347562CurzonF.AhlbornB.Efficiency of a Carnot engine at maximum power output197543122Arias-HernándezL. A.Angulo-BrownF.A general property of endoreversible thermal engines1997817297329792-s2.0-0001624884CarrolB. W.OstlieD. A.2007San Francisco, Calif, USAPearson-Adisson WesleyGutkowicz-KrusinD.ProcacciaI.RossJ.On the efficiency of rate processes. Power and efficiency of heat engines1978699389839062-s2.0-0006943760GustavsonM. R.Limits to wind power utilization1979204438813172-s2.0-0018330220CurryA.WesterP. J.1999Academic PressInternational Geiphysics SeriesLorentzE.196086New York, NY, USAPergamon PressSchulmanL.A theoretical study of the efficiency of the general circulation197734559580CaplanS. R.EssigA.1983Cambridge, Mass, USAHarvard University PressAngulo-BrownF.SantillánM.Calleja-QuevedoE.Thermodynamic optimality in some biochemical reactions1995171879010.1007/BF02451604Angulo-BrownF.Arias-HernándezL. A.SantillánM.On some connections between first order irreversible thermodynamics and finite-time thermodynamics200248182192SantillánM.Arias-HernándezL. A.Angulo-BrownF.Some optimization criteria for biological systems in linear irreversible thermodynamics1997191991092-s2.0-336450773981991215thWilliam Benton PublisherMullanD. J.2010Boca Raton, Fla, USACRC PressSturrockP. A.HolzerT. E.MihalasD. M.UlrichR. K.McCormacB. M.Physics of the Sun, Vol. I: The solar interior1986Calif, USALoockheed Palo Alto Research Laboratory