Next Article in Journal
System Design and Analysis of a Directly Air-Assisted Turbocharged SI Engine with Camshaft Driven Valves
Previous Article in Journal
The Use of Energy Storage Systems for Supporting the Voltage Needs of Urban and Suburban Railway Contact Lines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical and Experimental Study of Recycling Baffled Double-Pass Solar Air Heaters with Attached Fins

Energy and Opto-Electronic Materials Research Center, Department of Chemical and Materials Engineering, Tamkang University, Tamsui, New Taipei City 25137, Taiwan
*
Author to whom correspondence should be addressed.
Energies 2013, 6(4), 1821-1842; https://0-doi-org.brum.beds.ac.uk/10.3390/en6041821
Submission received: 25 January 2013 / Revised: 14 March 2013 / Accepted: 14 March 2013 / Published: 28 March 2013

Abstract

:
The study of the heat transfer of solar air heaters with a new design using an absorbing plate with fins and baffles, which facilitate the recycling of flowing air, is reported. The mathematical formulation and analytical analysis for such a recyclic baffled double-pass solar air heater were developed theoretically. The performance of the device was studied experimentally as well. The theoretical predicted and experimental results were compared with another design, i.e., a downward-type single-pass solar air heater without recycle and double-pass operations reported in our previous work. Significant improvement in heat-transfer efficiency is achieved with the baffle and fin design due to the recycling heating and the extended heat transfer area. The effects of mass flow rate and recycle ratio on the heat-transfer efficiency enhancement as well as on the power consumption increment are also discussed.

1. Introduction

The flat-plate solar air heater [1,2,3,4,5] is considered to be a simple device consisting of one or more glass (or transparent) covers situated above an absorbing plate with the air flowing over [4,5], under [1,2,3] or simultaneously over and under [6] the absorbing plate. The conventional flat-plate solar air heater has been investigated for heat-transfer efficiency improvement by introducing free convection [7] forced convection [8,9] extended heat-transfer area [10,11] and increased air turbulence [12,13]. Moreover, the adoption of the recycle-effect concept can effectively enhance the heat transfer rate and has been applied in many separation processes and reactor designs, as confirmed by the previous works [14,15,16,17].Extending from our previous work [18] this study focuses on a new design featuring recycling as well as turbulence and heat transfer area enhancement, as shown in Figure 1.
Figure 1. Double-pass solar air heater with internal baffles and fins.
Figure 1. Double-pass solar air heater with internal baffles and fins.
Energies 06 01821 g001
This work presents the theoretical mathematical formulations and the experimental verification for the new solar air heater design. The effects of the recycle ratio and air mass flow rate on the heat transfer rate, heat-transfer efficiency improvement, and the power consumption increment are also delineated.

2. Mathematical Formulation

The device proposed is conceptually depicted in Figure 1. An absorbing plate inserted into the air flow channel creates the upper and lower subchannels, which facilitates the double-pass air flow. On both sides of the absorbing plate, baffles and fins guide the air flow and provide turbulence and extra heat transfer area. Like in conventional solar air heaters, two layers of glass with a gap in between are used on top of the upper subchannel to reduce the heat loss to the environment and other sides of the device are thermally insulated.
Figure 2 shows the configuration of the device with recycle. The absorbing plate is inserted in parallel into a parallel-plate channel with width of W and length of L, and the open conduit is divided into two parts both with height of H. Before entering the lower subchannel, the inlet air with mass flow rate and inlet temperature is premixed with part of the air leaving the lower subchannel with mass flow rate and outlet temperature. The steady-state one-dimensional mathematical model was developed with the assumptions of: (1) the temperature of the absorbing plate, bottom plate and bulk fluids are functions of the flow direction only; (2) both glass covers and flowing air do not absorb radiant energy; (3) except the glass cover, all outside surfaces of the solar air collector are well insulated thermally.
Figure 2. Configuration of a double-pass solar air heater with recycle.
Figure 2. Configuration of a double-pass solar air heater with recycle.
Energies 06 01821 g002

2.1. Effect of Attaching Fins with Baffles

The total heat transfer from the absorbing plate with attached fins and baffles to a passing airflow can be written as the following via modifications from Bannett and Myers [19]:
q = h i ( A t + A f η f + A b η b ) ( T p T i ( z ) ) ,   i = a , b
where the fin efficiency is defined as:
η f = tanh  m f h f m f h f
in which:
m f = 2 h i ( L + t f ) k f A t f
The baffle efficiency η b is affected by the density of the baffles, i.e., W/De and l/L. Following empirical correlations [20,21] the baffle efficiency can be obtained for the fins attached to the lower and upper subchannels as:
η b = 26.361 ( W B D e ) 0.454 ( l L ) 0.634  for lower subchannel
and:
η b = 15.583 ( W B D e ) 0.0518 ( l L ) 0.227  for upper subchannel

2.2. Temperature Distributions for the Flowing Air in the Flow Direction

For the steady-state one-dimensional mathematical model, the energy balances for the flowing air streams, glass cover 1, absorbing plate and bottom plate, are given in Equations (6–10).
For the airflow in upper subchannel:
m ˙ R C p W d T b ( z ) d z = h b ϕ b ( T p T b ( z ) ) h b ( T b ( z ) T c 1 )
For the airflow in lower subchannel:
[ ( R + 1 ) m ˙ C p W ] d T a ( z ) d z = h a ϕ a ( T p T a ( z ) ) h a ( T a ( z ) T R )
For the glass cover 1 (inner cover):
h r , p c 1 ( T p T c 1 ) + h b ( T b ( z ) T c 1 ) = U c 1 s ( T c 1 T s )
For the absorbing plate:
I 0 α p τ g 2 = U T ( T p T s ) U B ( T p T s ) h b ϕ b ( T p T b ( z ) ) h a ϕ a ( T p T a ( z ) )
For the bottom plate
h a ( T a ( z ) T R ) + h r , p R ( T p T R ) = U B s ( T R T s )
The heat transfer coefficients between airflow and the conduit walls were assumed to be equal, in which:
ϕ a = 1 + A f , a A t η f , a + A b , a A t η b , a
ϕ b = 1 + A f , b A t η f , b + A b , b A t η b , b
Rearranging Equations (8–10) and solving for (TpTb(z)), (Tb(z)−Tc1), (TpTa(z)) and (Ta(z)−TR), and then, substituting the results into Equations (6,7), we have the temperature distributions with the dimensionless coordinate in the flow direction as follows:
d ( T b ( ξ ) T s ) d ξ = B 1 ( T b ( ξ ) T s ) + B 2 ( T a ( ξ ) T s ) + B 3
d ( T a ( ξ ) T s ) d ξ = B 4 ( T b ( ξ ) T s ) + B 5 ( T a ( ξ ) T s ) + B 6
ξ = z / L .
The boundary conditions for solving Equations (13,14) are:
ξ = 0 ,   T a , i 0 = T a , i + R T a , L ( 1 + R )
ξ = 1 ,   T a ( 1 ) = T a , L ,   T b ( 1 ) = T b , L
By following the similar mathematical treatment performed in our previous work [18], the analytical solutions to this device were obtained after transforming Equations (13,14) into ( D B 1 ) ( T b T s ) B 2 ( T a T s ) = B 3 and B 4 ( T b T s ) + ( D B 5 ) ( T a T s ) = B 6 with the use of the differential operator notation, in which D = d / d ξ . Operate on the second of these with ( 1 / B 4 ) ( D B 1 ) and add to the first equation to obtain Equation (18) with using boundary conditions in Equations (16),(17), and hence, Equation (19) was thus obtained. The temperature distributions of the flowing air for both lower and upper subchannels are thus obtained, respectively:
T b ( ξ ) = Y 1 B 5 B 4 C 1 e Y 1 ξ + Y 2 B 5 B 4 C 2 e Y 2 ξ B 5 ( B 3 B 4 B 1 B 6 ) B 4 ( B 1 B 5 B 2 B 4 ) B 6 B 4 + T s
T a ( ξ ) = C 1 e Y 1 ξ + C 2 e Y 2 ξ + B 3 B 4 B 1 B 6 B 1 B 5 B 2 B 4 + T s
All the coefficients Bi, Gi, Mi, Yi, Ci, Fi and Ii are referred to the Appendix. Once the inlet and outlet temperatures are known, the useful gain may be estimated from the relation:
Q u = m ˙ ( 1 + R ) C p ( T a , L T a , i 0 ) + m ˙ C p ( T b , o T b , L )
or:
Q u = m ˙ ( 1 + R ) C p ( T a , L T a , i ) + m ˙ C p ( T b , o T b , L ) = m ˙ ( T b , o T a , i )
Finally, the collector efficiency η D of a baffled double-pass solar air heater with external recycle and attached internal fins, relates the actual useful energy gain to the incident solar radiation, was defined as follows Ho et al. [18]:
η D = Q u ( U s e f u l   g a i n   o f   e n e r g y   c a r r i e d   a w a y   b y   a i r ) I 0 A c ( T o t a l   s o l a r   r a d i a t i o n   i n c i d e n t )
= m ˙ C p ( T b , o T a , i ) I 0 A c = τ g 2 α p U L ( T p , m T s ) / I 0
and the airflow outlet temperature is readily obtained from Equation (22), i.e.:
T b , o = T a , i + ( η D I 0 A c / m ˙ C p )
Equating Equations (22,23), we have:
T p , m = T s + ( I 0 τ g 2 α p / U L ) m ˙ C p ( T b , o T a , i ) A c U L
Substitution of Equation (24) into Equation (25) to eliminate T b , o results in:
T p , m = T s + ( I 0 / U L ) ( τ g 2 α p η D )
Equation (26) is the expression of T p , m in terms of η D .

2.3. Heat Transfer Coefficients

The total heat loss rate from the solar air heater to the surrounding is the summation of the top, edge and bottom loss rates, i.e.:
Q loss = U L A c ( T p , m T s ) = U T A c ( T p , m T s ) + U B A c ( T R , m T s ) + U E A E ( T a , m + T b , m 2 T s )
where the average temperatures T a , m , T b , m , T p , m and T R , m can be obtained by integration and rearrangement of Equations (10,18,19,26), and U T , U B and U E are estimated by Equations (28,32). An empirical equation for U T was developed [22] following the basic procedure in [23] as:
U T = { N C a i r T p , m [ ( T p , m T s ) N + f ] e + 1 h w } 1 + σ ( T p , m + T s ) ( T p , m 2 + T s 2 ) { ( ε p + 0.00591 N h w ) 1 + [ 2 N + f 1 + 0.133 ε p ] ε g N }
where:
f = ( 1 + 0.089 h w 0.1166 h w ε p ) ( 1 + 0.07866 N )
C a i r = 520 ( 1 0.00005 β 2 )  for  0 < β < 70 ;   for  β > 70  use  β = 70
e = 0.43 ( 1 100 / T p , m )
The edge and bottom loss coefficients depend mainly on the insulation thickness:
U B = U E = k s / l s
Moreover, all heat transfer coefficients in Equations (6,7) are estimated by the following empirical and correction equations. The forced convective heat-transfer coefficient for air flowing in the ducts between two flat plates may be estimated by Kay’s data [24,25]:
N u , i = h i D e , i / k = 0.0158 Re i 0.8 ,    i = a , b
The modification of McAdams [25] for turbulent flow in a short conduit:
N u , i = h i D e , i / k = 0.0158 Re i 0.8 [ 1 + ( D e , i / L ) 0.7 ] ,     i = a , b
The correlation equation for laminar flow was expressed by Heaton et al. [26]:
N u , i = 4.4 + 0.00398 ( 0.7 Re j D e , i / L ) 1.66 1 + 0.0114 ( 0.7 Re i D e , i / L ) 1.12 ,   i = a , b
in which the characteristic length is the equivalent diameter of the air duct:
D e , 0 = 4 ( 2 H W ) 2 ( 2 H + W ) ,   D e , a = D e , b = 4 ( H W 5 H f t f ) 2 ( H + W ) + 10 H f
for downward-type single-pass and double-pass devices, respectively. The average velocities of the downward-type single-pass device and the lower and upper subchannels of the double-pass device are:
v 0 ¯ = m 2 H W ρ , v a ¯ = m ( 1 + R ) ( H W 5 h f t f ) ρ ,   v b ¯ = m ( H W 5 h f t f ) ρ
Thus; from Equations (36,37); one obtains the Reynolds numbers of downward-type single-pass and double-pass devices; respectively; for the rectangular ducts:
Re 0 = ρ v 0 ¯ D e , 0 μ , Re a = ρ v a ¯ D e , a μ ,   Re b = ρ v b ¯ D e , b μ
The heat transfer coefficients and radiation coefficients in Equations (8–10) were calculated by the following empirical equations. The thermal resistance from glass cover 1 through glass cover 2 to the ambient air, may be expressed as:
1 U c 1 s = 1 U T 1 h r , p c 1 + h 1
or:
1 U c 1 s = 1 h w + h r , c 2 s + 1 h c 1 c 2 + h r , c 1 c 2
where the heat-transfer coefficient for free convection of air between two glass covers may be estimated using Hottel’s empirical Equation [23] as:
h c 1 c 2 = 1.25 ( T c 1 , m T c 2 , m ) 0.25
and the convective heat-transfer coefficient for air flowing over the outside surface of glass cover 2, can be calculated using the following empirical equation given by McAdams [25]:
h w = 2.8 + 3.0   V
The radiation coefficients between the two air-duct surfaces may be estimated by assuming a mean radiant temperature equal to the mean fluid temperature, viz.:
h r , p c 1 4 σ T b , m 3 / [ ( 1 / ε p ) + ( 1 / ε c 1 ) 1 ]
h r , p R 4 σ T a , m 3 / [ ( 1 / ε P ) + ( 1 / ε R ) 1 ]
while those between two glass covers and from glass cover 2 to the ambient are, respectively:
h r , c 1 c 2 = σ ( T c 1 , m 2 + T c 2 , m 2 ) ( T c 1 , m + T c 2 , m ) ( 1 / ε g ) + ( 1 / ε g ) 1
h r , c 2 s = ε g σ ( T c 2 , m 2 + T s 2 ) ( T c 2 , m + T s )

2.4. Theoretical Calculation Method for Collector Efficiency

The calculation procedures for theoretical predictions of collector efficiency and collector efficiency improvement were described as follows. Firstly, with the given collector geometries ( A c = 0.09 m2, H = 0.05 m, hB = hf = 0.05 m, L = 0.3 m, W = 0.3 m, l = 0.12 m, lB = ls =0.06 m, tf = 0.002 m, WB = 0.03 m) and the system properties [ α p = 0.96, ε g = 0.94, ε p = 0.8, ε R = 0.94, τ g = 0.875, kB = kf = 14.9 W/(m K), ks = 0.033 W/(m K)], as well as the operating conditions [ m ˙ = 38.52, 57.96 and 77.04 kg/hr, Tin = (30 ± 0.1) °C, Ts = (30 ± 0.1) °C, I0 = (830 ± 40), and (1100 ± 40) W/m2, V = 1.0 m/s], the initial guesses for the average temperatures Tp,m, TR,m, Tc1,m, Tc2,m are estimated to obtain the heat transfer coefficients and to calculate the temperature distributions of both upper and lower subchannels from Equations (13–17).Secondly, the values of Tp,m, TR,m, Tc1,m, Tc2,m are then checked by using Equations (10–12,26), respectively, and new iteration values of Tp,m, TR,m, Tc1,m, Tc2,m may be obtained. If the new calculated values of Tp,m, TR,m, Tc1,m, Tc2,m are different from the previous iteration values, re-calculated with the proper equations by iteration procedures until the error tolerance is small enough (say less than 10−3), and thus the collector efficiencies ηD can be calculated by using Equation (23).The application of recycle-effect concept to a solar air heater actually increases the convective heat-transfer coefficient by producing a preheating effect with the hot outgoing air. Results show that this advantage influence can partially or completely offset the decrement of teat transfer driving force (temperature gradient).The collector efficiency improvement ID with comparing the collector efficiency between downward-type single-pass device and the present design of solar air heaters is:
I D = η D η S η S × 100 %
in which ηD and ηS, respectively, denote the collector efficiency in baffled double-pass type with fins attached and downward-type single-pass solar air heaters.

3. Apparatus and Experimental Methods

A baffled double-pass solar air heater with external recycle (stainless steel of 7 × 10−5 m thickness) and with five fins attached on both above and under the absorbing plate for experimental studies, are shown in Figure 1 and Figure 2. The collector (0.3 m long, 0.3 m wide, 0.05 m in height of each subchannel) consists of two glass covers (2.75 × 10−3 m thickness and 0.05 m space), a black absorbing plate, baffles welded to the fins, 42 electrical energy supplies (110 V, 125 W) in the part B in Figure 3 and thermal insulation using Styrofoam (0.06 m thickness), as shown in Figure 3 for artificial simulation. The incident solar radiationI0(arrived irradiation)were measured and recorded with an Epply Laboratory Inc. pyranometer and an Instantaneous Solar Radiation meter Model No. 455 also from Epply Laboratory Inc. Before air entered the lower subchannel, a pre-heater was used to ensure a specified temperature [i.e., Tin = (30 ± 0.1) °C] at the inlet to mix the recycle airflow from the lower subchannel regulated by an adjustable valve. There are two air boxes (0.3 m × 0.3 m × 0.05 m) with 14 holes (0.015 m diameter) were used for air distribution at inlets of both upper and lower subchannels, respectively. Except for the outer glass cover, all parts of the solar air heater were well insulated thermally to minimize the energy loss.
Figure 3. Schematic diagram of a double-pass solar air heater with artificial simulation.
Figure 3. Schematic diagram of a double-pass solar air heater with artificial simulation.
Energies 06 01821 g003
Airflow was steadily supplied by a blower (Model BL model 552, Redmond Co., Owosso, MI, USA) and the mass flow rate controlled by a transformer was measured with Model 24-611 anemomaster from Kanomax Japan Inc. (Osaka, Japan). The air temperatures of absorbing plate and at the inlet and outlet of the solar collector were measured with thermocouples. The ambient temperature [Ts = (30 ± 0.1) °C] measured at a position 0.15 m above the outer glass cover with a thermocouple was monitored by using an air conditioner with a fan to provide wind flow of 1.0 m/s (as indicated by a Chicago Cenco anemometer, TSI Inc., St. Paul, MN, USA) in the experimental apparatus.

4. Results and Discussion

The essential physical properties employed in this study are presented in Table 1, and the operating conditions are Ta,I = 303 K, and m ˙ = 38.52, 57.96 and 77.04 kg/h. By substituting the specified values into the appropriate equations, theoretical predictions for collector efficiency and airflow outlet temperature with incident solar radiation, mass flow rate and recycle ratio as parameters were obtained. The theoretical predictions and experimental results are plotted in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9.
Table 1. Physical properties of air at 1atm [8].
Table 1. Physical properties of air at 1atm [8].
T (K)ρ(kg/m3)Cp [J/(kg K)]k [W/(m K)]μ [kg/(m s)]
2731.29210060.02421.72 × 10−5
2931.20410060.02571.81 × 10−5
3131.12710070.02721.90 × 10−5
3331.05910080.02871.99 × 10−5
3530.99910100.03022.09 × 10−5
Figure 4. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 830 W/m2, m ˙ = 38.52 kg/h.
Figure 4. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 830 W/m2, m ˙ = 38.52 kg/h.
Energies 06 01821 g004
Figure 5. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 830 W/m2, m ˙ = 57.96 kg/h.
Figure 5. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 830 W/m2, m ˙ = 57.96 kg/h.
Energies 06 01821 g005
Figure 6. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 830 W/m2, m ˙ = 77.04 kg/h.
Figure 6. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 830 W/m2, m ˙ = 77.04 kg/h.
Energies 06 01821 g006
Figure 7. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 1100 W/m2, m ˙ = 38.52 kg/h.
Figure 7. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 1100 W/m2, m ˙ = 38.52 kg/h.
Energies 06 01821 g007
Figure 8. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 1100 W/m2, m ˙ = 57.96 kg/h.
Figure 8. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 1100 W/m2, m ˙ = 57.96 kg/h.
Energies 06 01821 g008
Figure 9. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 1100 W/m2, m ˙ = 77.04 kg/h.
Figure 9. Effect of recycle ratio on collector efficiency in a double-pass flat-plate solar air heater with recycle; I0 = 1100 W/m2, m ˙ = 77.04 kg/h.
Energies 06 01821 g009
This experimental study for recyclic baffled solar air heaters with attaching internal fins was conducted under artificial simulation, as shown in Figure 3. The agreement between the experimental results and theoretical predictions is fairly good, and the operation with recycle-effect device substantially improves the collector efficiency, as confirmed by Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9. In order to validate the accuracy of the experimental results as compared to theoretical predictions, the following definition of error analysis was made:
E = 1 N exp i = 1 N exp | η t h e o , i η exp , i | η t h e o , i
where Nexp, ηtheo,i and ηexp,i indicate the experimental measurements, theoretical predictions and experimental data of collector efficiencies, respectively. The error analysis may be calculated with the results in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 by using Equation (48). The error analysis of the experimental measurements is of 3.19 × 10−2E ≦ 9.86 × 10−2 and presented in Table 2.
The theoretical predictions of I D for the system of the present interest were also calculated by Equation (47), and the results are also listed in Table 3 and Table 4. It is shown that the considerable improvement in collector efficiency is obtained if the operation is carried out with fins and baffles attached under external recycle.
Table 2. The accuracy analysis of the experimental results.
Table 2. The accuracy analysis of the experimental results.
I0 (W/m2) m ˙ (kg/h)RecycleRecycle + FinsRecycle + Fins + Baffles
83038.528.79× 10−26.80× 10−23.20× 10−2
57.963.83× 10−25.24× 10−23.44× 10−2
77.044.41× 10−23.40× 10−23.19× 10−2
110038.529.86× 10−28.53× 10−23.96× 10−2
57.964.57× 10−28.50× 10−23.41× 10−2
77.044.41× 10−25.10× 10−24.18× 10−2
Table 3. Theoretical predictions of heat-transfer efficiency improvement ID for I0 = 830 W/m2.
Table 3. Theoretical predictions of heat-transfer efficiency improvement ID for I0 = 830 W/m2.
RRecycleRecycle + FinsRecycle + Fins + Baffles
m ˙ (kg/h)
-38.5257.9677.0438.5257.9677.0438.5257.9677.04
0.0029.139.195.8030.4240.5415.7088.3579.1965.46
0.2550.1641.3535.2752.7552.9736.47111.9782.1667.15
0.5067.6454.8646.3869.5860.5447.10117.4884.0569.81
0.7578.9662.9752.6680.5865.6853.38120.0685.1471.01
1.0086.7368.3856.7688.3568.9257.25121.6885.6871.50
1.2592.8872.4359.6693.8572.1659.90122.9886.4971.98
1.5097.4175.1461.5998.3874.0562.08123.6286.7672.22
1.75100.9777.5763.29101.6274.8663.53123.9587.0372.46
Table 4. Theoretical predictions of heat-transfer efficiency improvement ID for I0 = 1100W/m2.
Table 4. Theoretical predictions of heat-transfer efficiency improvement ID for I0 = 1100W/m2.
RRecycleRecycle + FinsRecycle + Fins + Baffles
m ˙ (kg/h)
-38.5257.9677.0438.5257.9677.0438.5257.9677.04
0.0029.229.217.1130.1940.6517.1684.4279.4067.16
0.2550.3241.4637.2552.9253.1238.48111.6982.3869.36
0.5067.8655.0148.2869.8160.7049.02117.5384.2872.06
0.7578.9063.1454.6680.5265.8555.15120.1385.3773.28
1.0087.0168.5658.8288.3169.1159.31122.0885.9173.77
1.2592.8672.3661.7693.8372.3662.01123.0586.7274.26
1.5097.4075.3463.7398.3874.2564.22123.7086.9974.51
1.75100.9777.5165.44101.9575.0765.69124.0387.2674.75
2.00103.9079.4066.67104.5575.6166.91124.3587.5375.00
The collector efficiency improvement increases with increasing incident solar radiation and recycle ratio as well as mass flow rate owing to turbulence increment and forced convective heat-transfer coefficient enhancement, as shown in Table 3 and Table 4, especially for operating at larger air flow rate m ˙ . The improvement in collector efficiency by recycle operation is obtained more than 64% as compared to downward-type single-pass device. In addition, the larger incident solar radiation is associated with the higher absorbing plate temperature of the solar air collector, resulting in the increment of the heat transfer efficiency. For instances, the heat-transfer efficiency improvements with fins and baffles for m ˙ = 77.04 kg/h and R = 2, are 72.71% and 75.0% for I 0 = 830 and 1100 W/m2, respectively. Accordingly, the order of heat transfer efficiencies of flat plate collectors are: collector with fins and baffles> collector with fins > collector without fins and baffles, as indicated in Table 3 and Table 4.
The power consumptions of the solar air heater with fins attached are calculated by the correlation equation. The power consumptions for downward-type single-pass and double-pass solar air heaters without attached baffles, respectively, are as follows:
P S = m w f , s = m 2 f F v 0 ¯ 2 L D e
P D = m w f , a + m ( 1 + R ) w f , b = m 2 f F v a ¯ 2 L D e , a + m ( 1 + R ) 2 f F v b ¯ 2 L D e , b
in which:
f F = 24 Re ,  for laminar flow
f F = 0.0790 Re 0.25 ,  for turbulent flow ,  Re < 10 5
The power consumption increment, I P , may be defined as:
I P = P D P S P S
The hydraulic dissipated energy in a single-pass device without recycle are P S = 7.68 × 10 6 , 2.36 × 10 5 and 5.23 × 10 5 W, respectively, for m ˙ = 38.52 , 57.96 and 77.04 kg/h for mass flow rate. The suitable selections of the operating parameters (R and m ˙ ) on considering of both the heat-transfer efficiency improvement, I D , and the power consumption increment, I P , are presented in Figure 10.
Figure 10. ID/IP vs. recycle ratio R with mass flow rate as a parameter.
Figure 10. ID/IP vs. recycle ratio R with mass flow rate as a parameter.
Energies 06 01821 g010
The economically optimal design must take into account the costs of energy and devices. However, for simplicity, only the power consumption is considered in this study. The index of I D / I P is adopted for comparing different designs. Operating double-pass solar air heaters with recycle at R o p t i m a l 0.5 gives the maximum value of I D / I P , as seen from Figure 10, while the values of I D / I P decrease as the recycle ratio moves away from 0.5. A remark must be stated that even if only energy criterion is to be considered, because thermal energy and electric power are different in terms of the quality of energy, employing exergy analysis can provide different results from those shown in Figure 10.
The present work is actually the extension previous work [18] except the configuration of recycle. Figure 11 illustrates the graphical representation for comparisons with some experimental results and theoretical predictions obtained in [18] under the same design and operating parameters. With those comparisons, the advantage this study is evident for all mass flow rates and recycle ratios.
Figure 11. Comparisons of collector efficiency among the present and previous works [18].
Figure 11. Comparisons of collector efficiency among the present and previous works [18].
Energies 06 01821 g011

5. Conclusions

As demonstrated by theoretical and experimental studies, the introduction of external recycle has positive effects on the heat transfer efficiency of solar air heaters. Moreover, the advantages of baffled solar air heaters with internal fins attached provide higher turbulence and extending the heat transfer area, and thus, the enhancement of heat-transfer efficiency. It is believed that the availability of such a mathematical formulation as developed here for a recycling baffled solar air heater with internal fins attached is an important contribution to the design and analysis of double-pass devices with internal or external recycle. Although the lager mass flow rate and recycle ratio are beneficial to collector efficiency improvement, higher energy dissipation occurs at the same time when operating at a higher mass flow rate and recycle ratio. The value of I D / I P vs. recycle ratio R was represented with the mass flow rate as a parameter in Figure 10 to recommend an optimal design. After all, both the heat-transfer efficiency improvement and power consumption increment are significant factors to be taken into account simultaneously in designing recyclic double-pass solar air heaters.

Acknowledgments

The authors wish to thank the National Science Council of the Republic of China for its financial support.

Nomenclature

A c
surface area of the collector = LW (m2)
A b
total surface area of attached baffles (m2)
A E
surface area of the edge of collector = 2LH (m2)
A f
total surface area of attached fins (m2)
A t
the area of absorbing plate less the welded area of fins (m2)
A t f
cross section area of fins (m2)
Bi
coefficients defined in Equations (A16–A21)
Ci
coefficients defined in Equations (A26,A27)
C p
specific heat of air at constant pressure (J/kg K)
D e , 0
equivalent diameter of downward-type single-pass device defined in Equation (36) (m)
D e , a
equivalent diameter of lower subchannel of double-pass device defined in Equation (36) (m)
D e , b
equivalent diameter of upper subchannel of double-pass device defined in Equation (36) (m)
Fi
coefficients defined in Equations (A28–A30)
f F
Fanning friction factor
G i
coefficients defined in Equations (A5–A8) and (A13–A15)
H
height of both upper and lower channels (m)
Hf
height of fin (m)
HB
height of baffles (m)
Hf
height of fins (m)
ha
heat-transfer coefficient of air in the lower subchannel (W/m2 K)
hb
heat-transfer coefficient of air in the upper subchannel (W/m2 K)
h r
radiative heat-transfer coefficient between two parallel plates (W/m2 K)
h c 1 c 2
heat transfer coefficient for free convection of air between two glass covers, defined in Equation (44)(W/m2 K)
h r , c 1 c 2
radiation heat transfer coefficient between two covers, defined in Equation (45) (W/m2 K)
h r , c 2 s
radiation heat transfer coefficient from cover 2 to the ambient, defined in Equation (46)(W/m2 K)
h r , p c 1
radiation heat transfer coefficient between cover 1 and absorbing plate, defined in Equation (43) (W/m2 K)
h r , p R
radiation heat transfer coefficient between absorbing plate andbottom plate, defined in Equation (44) (W/m2 K)
h w
convective heat-transfer coefficient for air flowing over the outside surface of glass cover (W/m2 K)
Ii
coefficients defined in Equations (A31,A32)
I 0
incident solar radiation (W/m2)
I D
percentage of collector efficiency improvement, defined in Equation (47)
I P
the power consumption increment,, defined in Equation (53)
k
thermal conductivity of air (W/m K)
k B
thermal conductivity of bottom plate (W/m K)
k f
thermal conductivity of fins (W/m K)
k s
thermal conductivity of insulator (W/m K)
L
channel length (m)
l
distance between baffles (m)
l B
thickness of bottom plate (m)
l s
thickness of insulator (m)
w f , a
lower subchannel friction loss of double-pass solar air heater (J/kg)
w f , b
upper subchannel friction loss of double-pass solar air heater (J/kg)
w f , s
friction loss of downward-type single-pass solar air heater (J/kg)
m ˙
total air mass flow rate (kg/h)
m f
coefficient defined in Equation (3)
M a
coefficient defined in Equation (A23) (W/m2 K)
M b
coefficient defined in Equation (A22) (W/m2 K)
N
number of glass cover
N exp
the number of experimental measurements
N u
Nusselt number
PB
power consumption of baffled solar air heaters (W)
PD
power consumption of baffled double-pass solar air heaters (W)
PS
power consumption of downward-type solar air heaters (W)
Q u
useful energy gain carried away by air (W)
Qloss
total heat loss rate from the air collector to the surrounding area(W)
q
total heat flow rate of solar air collector with fins attached (W)
q
total heat flow rate of solar air collector with baffles attached (W)
R
recycle ratio
Re 0
Reynolds number in downward-type single-pass solar air heater
Re a
Reynolds number in lower subchannel of double-pass solar air heater
Re b
Reynolds number in upper subchannel of double-pass solar air heater
T i n
inlet air temperature (K)
T a ( z )
axial fluid temperature distribution in lower subchannel (K)
T b ( z )
axial fluid temperature distribution in upper subchannel (K)
T c
temperature of glass cover (K)
T a , m
mean fluid temperature in lower subchannel (K)
T b , m
mean fluid temperature in upper subchannel (K)
T c , m
mean temperature of glass cover (K)
T p
temperature of absorbing plate (K)
T p , m
mean temperature of absorbing plate (K)
T R
temperature of bottom plate (K)
T R , m
mean temperature of bottom plate (K)
T s
ambient temperature (K)
t f
thickness of fins (m)
U B
loss coefficient from the bottom plate of solar air heater to the ambient environment (W/m2 K)
U B s
loss coefficient from the surfaces of edges and the bottom plate of the solar collector to the ambient environment (W/m2 K)
U c 1 s
loss coefficient from glass cover 1 to the ambient environment (W/m2 K)
U E
loss coefficient from the edge of solar air heater to the ambient environment (W/m2 K)
U L
overall loss coefficient (W/m2 K)
U T
loss coefficient from the top of solar air heater to the ambient environment (W/m2 K)
W
width of both upper and lower subchannels (m)
WB
baffle width (m)
v 0 ¯
mean air velocity in the downward-type single-pass solar air heater (m/s)
v a ¯
mean air velocity in lower subchannel of double-pass solar air heater (m/s)
v b ¯
mean air velocity in upper subchannel of double-pass solar air heater (m/s)
V
wind velocity (m/s)
Yi
coefficients defined in Equations (A24,A25)
z
axial coordinate (m)

Greek Letters 

α p
absorptivity of absorbing plate
η b
baffle efficiencies of upper and lower channels
η D
collector efficiency of downward-type single-pass solar air heater
η exp , i
experimental data of collector efficiency
η f
fin efficiency, defined in Equation (2)
η t h e o , i
theoretical prediction of collector efficiency
τ g
transmittance of glass cover
ε g
emissivity of glass cover
ε R
emissivity of bottom plate
ε p
emissivity of absorbing plate
ρ
air density (kg/m3)
μ
air viscosity (kg/s m)
ξ
dimensionless channel length, defined in Equation (15)
σ
Stefan-Boltzmann constant = 5.67 × 10−8 (W/m2 K4)

Appendix

Rearranging Equations (8,9), one obtains, respectively:
h r , P c 1 [ ( T P T b ( z ) ) + ( T b ( z ) T c 1 ) ] + h b ( T b ( z ) T c 1 ) U c 1 s [ ( T b ( z ) T s ) ( T b ( z ) T c 1 ) = 0
S 0 α p τ g 2 U T [ ( T P T b ( z ) ) + ( T b ( z ) T s ) ] U B [ ( T P T a ( z ) ) + ( T a ( z ) T s ) ] h b ( T P T b ( z ) ) h a [ ( T P T b ( z ) ) + ( T b ( z ) T s ) ( T a ( z ) T s ) ] = 0
Solving Equations (A1),(A2) simultaneously for ( T p T b ( z ) ) and ( T b ( z ) T c 1 ) , we have:
T P T b ( z ) = G 1 ( T b ( z ) T s ) + G 2 ( T b ( z ) T s ) + G 3
T b ( z ) T c 1 = G 4 [ ( U c 1 s + G 1 h r , P c 1 ) ( T b ( z ) T s ) G 2 h r , P c 1 ( T a ( z ) T s ) h r , P c 1 G 3 ]
where:
G 1 = ( h a ϕ a + U T + U B ) / ( U T + U B + h a ϕ a + h b ϕ b ) ,
G 2 = h a ϕ a / ( U T + U B + h a + h b ) ,
G 3 = S 0 a p τ g 2 / ( U T + U B + h a ϕ a + h b ϕ b )
G 4 = ( h r , p c 1 + h b + U c 1 s ) 1 .
Again, rearranging Equations (9,10), one obtains, respectively:
S 0 α p τ g 2 U T [ ( T P T a ( z ) ) + ( T a ( z ) T s ) ] U B [ ( T P T a ( z ) ) + ( T a ( z ) T s ) ]
h b ϕ b [ ( T P T a ( z ) ) + ( T a ( z ) T s ) ( T b ( z ) T s ) ] h a ϕ a ( T P T a ( z ) ) = 0
h r , P R [ ( T P T a ( z ) ) + ( T a ( z ) T R ) ] + h 1 ( T a ( z ) T R ) U B s [ ( T a ( z ) T s ) ( T a ( z ) T R ) = 0
Solving Equations (A9,A10) simultaneously for ( T p T a ( z ) ) and ( T a ( z ) T R ) , we have:
T P T a ( z ) = G 5 ( T a ( z ) T s ) + G 6 ( T b ( z ) T s ) + G 3
T a ( z ) T R = G 7 [ ( U B s + G 5 h r , P R ) ( T a ( z ) T s ) G 6 h r , P R ( T b ( z ) T s ) h r , P R G 3 ]
where:
G 5 = ( h b ϕ b + U T + U B ) / ( U T + U B + h a ϕ a + h b ϕ b )
G 6 = h b 1 / ( U T + U B + h a ϕ a + h b ϕ b )
G 7 = ( h r , P R + h a + U B s ) 1
Substitution of Equations (A3,A4) into Equation (6) and substitution of Equations (A11,A12) into Equation (7) result in Equations (13,14), respectively, in which:
B 1 = ( h r , p c 1 h 1 G 3 G 4 h 1 ϕ 1 G 3 U c 1 s h 1 G 4 ) / M b
B 2 = ( h 1 ϕ 1 G 2 + h r , p c 1 h 1 G 2 G 4 ) / M b
B 3 = ( h 1 ϕ 1 G 3 + h r , p c 1 h 1 G 3 G 4 ) / M b
B 4 = ( h 2 ϕ 2 G 6 + h r , p R h 2 G 6 G 7 ) / M a
B 5 = ( h r , p R h 2 G 5 G 7 h 2 ϕ 2 G 5 U B s h 2 G 7 ) / M a
B 6 = ( h 2 ϕ 2 G 3 + h r , p R h 2 G 3 G 7 ) / M a
M b = [ m ˙ C P W L ] = m ˙ C p / A c
M a = [ m ˙ ( 1 + R ) C P W L ] = m ˙ ( 1 + R ) C p / A c
Y 1 = ( B 1 + B 5 ) + ( B 1 B 5 ) 2 + 4 B 2 B 4 2
Y 2 = ( B 1 + B 5 ) ( B 1 B 5 ) 2 + 4 B 2 B 4 2
C 1 = [ ( 1 + R R   e Y 2 ) ( F 2 + F 3 ) / F 1 ] / ( e Y 2 I 2 B 4 ) F 2 F 1 { 1 + R R   e Y 1 [ e Y 1 I 1 B 4 ( 1 + R R   e Y 2 ) ] } / ( e Y 2 I 2 B 4 )
C 2 = [ ( 1 + R R   e Y 1 ) ( F 2 + F 3 ) / F 1 ] / ( e Y 1 I 1 B 4 ) F 2 F 1 { 1 + R R   e Y 2 [ e Y 2 I 2 B 4 ( 1 + R R   e Y 1 ) ] } / ( e Y 1 I 1 B 4 )
F 1 = B 1 B 5 B 2 B 4
F 2 = B 3 B 4 B 1 B 6
F 3 = B 3 B 5 B 2 B 6
I 1 = Y 1 B 4 B 5
I 2 = Y 2 B 4 B 5

References

  1. Close, D.J.; Dunkle, R.V. Behaviour of adsorbent energy storage beds. Sol. Energy 1976, 18, 287–292. [Google Scholar] [CrossRef]
  2. Liu, C.H.; Sparrow, E.M. Convective-radiative interaction a parallel plate channel-application to air-operated solar collectors. Int. J. Heat Mass Transf. 1980, 23, 1137–1146. [Google Scholar] [CrossRef]
  3. Seluck, M.K. Solar Air Heaters and Their Applications; Sayigh, A.A.M., Ed.; Academic Press: New York, NY, USA, 1977. [Google Scholar]
  4. Tan, H.M.; Charters, W.W.S. Experimental investigation of forced-convective heat transfer for fully-developed turbulent flow in a rectangular duct with asymmetric heating. Sol. Energy 1970, 13, 121–125. [Google Scholar] [CrossRef]
  5. Whillier, A. Plastic covers for solar collectors. Sol. Energy 1963, 7, 148–154. [Google Scholar] [CrossRef]
  6. Yeh, H.M.; Ho, C.D.; Hou, J.Z. Collector efficiency of double-flow solar air heaters with fins attached. Energy 2002, 27, 715–727. [Google Scholar] [CrossRef]
  7. Kreith, F.; Kreider, J.F. Principles of Solar Engineering, 2nd ed.; McGraw-Hill: New York, NY, USA, 1978. [Google Scholar]
  8. Duffie, J.A.; Beckman, W.A. Solar Engineering of Thermal Processes, 3rd ed.; Wiley: New York, NY, USA, 1980. [Google Scholar]
  9. Tonui, J.K.; Tripanagnostopoulos, Y. Improved PV/T solar collectors with heat extraction by forced or natural air circulation. Renew. Energy 2007, 32, 623–637. [Google Scholar] [CrossRef]
  10. Gao, W.; Lin, W.; Liu, T.; Xia, C. Analytical and experimental studies on the thermal performance of cross-corrugated and flat-plate solar air heaters. Appl. Energy 2007, 84, 425–441. [Google Scholar] [CrossRef]
  11. Mohamad, A.A. High efficiency solar air heater. Sol. Energy 1997, 60, 71–76. [Google Scholar] [CrossRef]
  12. Verma, S.K.; Prasad, B.N. Investigation for the optimal thermohydraulic performance of artificially roughened solar air heaters. Renew. Energy 2000, 20, 19–36. [Google Scholar] [CrossRef]
  13. Yeh, H.M. Theory of baffled solar air heaters. Energy 1992, 17, 697–702. [Google Scholar] [CrossRef]
  14. Garg, H.P.; Sharma, V.K.; Bhargava, A.K. Theory of multiple-pass solar air heaters. Energy 1985, 10, 589–599. [Google Scholar] [CrossRef]
  15. Ho, C.D.; Yang, W.Y. An analytical study of heat-transfer efficiency in laminar counterflow concentric circular tubes with external refluxes. Chem. Eng. Sci. 2003, 58, 1235–1250. [Google Scholar] [CrossRef]
  16. Ho, C.D.; Yeh, C.W.; Hsieh, S.M. Improvement in device performance of multi-pass flat-plate solar air heaters with external recycle. Renew. Energy 2005, 30, 1601–1621. [Google Scholar] [CrossRef]
  17. Ho, C.D.; Yeh, H.M.; Chiang, S.C. Mass-Transfer enhancement in double-pass mass exchangers with external refluxes. Ind. Eng. Chem. Res. 2001, 40, 5839–5846. [Google Scholar] [CrossRef]
  18. Ho, C.D.; Yeh, H.M.; Wang, R.C. Heat-transfer enhancement in double-pass flat-plate solar air heaters with recycle. Energy 2005, 30, 2796–2817. [Google Scholar]
  19. Bannett, C.O.; Myers, J.E. Momentum, Heat, and Mass Transfer, 3rd ed.; McGraw-Hill: New York, NY, USA, 1982. [Google Scholar]
  20. Chou, W.H. A Study in Collector Efficiency of Baffled Solar Air Heaters. M.S. Thesis, Cheng Kung University, Tainan, Taiwan, 1987. [Google Scholar]
  21. Yeh, H.M. Energy balances for upward-type baffled solar air heaters. Energy 1994, 19, 919–924. [Google Scholar] [CrossRef]
  22. Klein, S.A. Calculation of monthly average transmittance-absorptance product. Sol. Energy 1979, 23, 547–551. [Google Scholar] [CrossRef]
  23. Hottel, H.C.; Woertz, B.B. Performance of flat-plate solar-heat collectors. Trans. ASME 1942, 64, 91–104. [Google Scholar]
  24. Kays, W.M. Convective Heat and Mass Transfer, 2nd ed.; McGraw-Hill: New York, NY, USA, 1980. [Google Scholar]
  25. McAdams, W.H. Heat Transmission, 3rd ed.; McGraw-Hill: New York, NY, USA, 1954. [Google Scholar]
  26. Heaton, H.S.; Reynolds, W.C.; Kays, W.M. Heat transfer in annular passages. Simultaneous development of velocity and temperature fields in laminar flow. Int. J. Heat Mass Transf. 1964, 7, 763–781. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Ho, C.D.; Chang, H.; Wang, R.C.; Lin, C.S. Analytical and Experimental Study of Recycling Baffled Double-Pass Solar Air Heaters with Attached Fins. Energies 2013, 6, 1821-1842. https://0-doi-org.brum.beds.ac.uk/10.3390/en6041821

AMA Style

Ho CD, Chang H, Wang RC, Lin CS. Analytical and Experimental Study of Recycling Baffled Double-Pass Solar Air Heaters with Attached Fins. Energies. 2013; 6(4):1821-1842. https://0-doi-org.brum.beds.ac.uk/10.3390/en6041821

Chicago/Turabian Style

Ho, Chii Dong, Hsuan Chang, Rei Chi Wang, and Chun Sheng Lin. 2013. "Analytical and Experimental Study of Recycling Baffled Double-Pass Solar Air Heaters with Attached Fins" Energies 6, no. 4: 1821-1842. https://0-doi-org.brum.beds.ac.uk/10.3390/en6041821

Article Metrics

Back to TopTop