Heat and Mass Transfer and Thermo Physical Analysis of Water Distillation Using Passive Solar Still ()
1. Introduction
In many areas of the world especially in developing countries, the shortage of drinking water for livelihood is an important problem. An immediate attention is required as the drinking water resources are being exploited and polluted rapidly by human beings [1]. The potable water is very scarce and the establishment of a human habitat in these areas strongly depends on how such water can be made available in arid regions. The importance of supplying potable water can hardly be over stressed [2]. Many techniques are available for purification of saline water among them solar still is simple and economical [3] [4]. Two major advantages favor the use of solar stills—clean and free energy and friendly to environment [5] [6]. Also the availability of clean water is a necessity for reducing the spread of diseases in these countries [7] [8]. Solar still is the widely used solar desalination device. Single slope solar still is the basic model of solar still and it is shown in Figure 1. The solar still at its lower part consists of basin which is insulated with insulating material to reduce the thermal losses to the surroundings. An inlet is provided for the access of saline water. Condensing covers are located on the top of the still. An outlet is provided to collect the distilled water [9].
2. Parameters of Solar Still Influencing the Performance
The solar still productivity and efficiency depends on the various parameters like solar radiation intensity, temperature, basin water depth, glass cover material thickness, and its inclination, wind velocity and the heat capacity of the still [9]. The parameters such as solar radiation intensity, temperature and wind velocity are uncontrolled as they are metrological parameters. Other parameters such as collector area, aperture ratio, water depth, glass thickness and insulation are controlled since they are managed by humans. Omar et al. [10] revealed that the productivity and efficiency of the solar still improves with increase in intensity. Also due to higher solar intensity the overall heat loss coefficient increases until it reaches the maximum. Kalidasa Murugavel et al. [11] suggested that the temperature difference between the water and glass is high, the productivity increases.
Figure 1. Skeleton of a single basin solar.
Ahmed et al. [12] conducted experimental analysis on double slope solar still based on temperature and also verified its numerical results. The productivity is enhanced by increasing the temperature of water usingwater heater in the base of solar still [13] [14]. Velmurugan et al. [15] inferred that an increase in collector area results in more productivity. Bilal et al. [16] concluded that the productivity of the solar still decreases linearly with increase in depth of water at the rate of 0.4 lit./cm. and found that the optimum tilt angle for more production as 35˚ during the month of May in Jordan, where the latitude angle is 31˚. Singh and Tiwari [17] and Velmurugan et al. [18] found that the annual yield of the solar still is maximum when the condensing glass coverinclination is equal to the latitude of the place. Solar still with glass cover plate with 3 mm thickness gives 15.5% more production than the cover of 6 mm glass plate thickness [19].
3. Experimental Setup
A single basin solar still made up of Copper sheet was designed for the basin size of 900 × 300 × 50 mm and 2 mm thick as shown in Figure 2. The Copper has higher thermal conductivity and the rate of heat transfer to water in the still is more. The basin is enclosed in a wooden box of inner cross section 1050 × 350 × 430 mm and thickness 15 mm. Plywood is used as outer cover to keep the still basin inside. The gap of 25 mm between the sides of the tray and the wooden box is filled with saw dust to prevent heat loss. Copper sheet was made in to a rectangular tray by sheet metal work of bending and cutting. The effective area of the solar still for saline water is 0.27 m2. The height of the solar still was taken 15 cm for lower end of glass. The glass of size 1175 × 320 mm was used as the roof for the still. The top of the basin is covered with a 5 mm thick transparent glass of mm. The cover is sealed tightly using silicon sealant to reduce the vapor leakage and the basin becomes air tight. The still is positioned in such a way that the sloping sides always face the sun. The whole experimental setup was kept in the North-South direction, with the inclination of 9˚. The solar stills were properly oriented, and directly exposed to the solar radiation. The measured parameters and quantities were: the solar radiation intensity (I), the glass temperature (Tg), the basin water temperature (Tw), and the ambient air temperature (Ta).
Figure 2. Solar still made up of copper.
4. Experimental Procedure
The measurements of the solar intensity, wind velocity, ambient, glass inside, glass outside, vapor, water, basin temperature and distillate yield are taken hourly to study the effect of each parameter on the still productivity. The water level in the solar basin is maintained to a level of 1 cm for constant basin water mass of 2.7 kg. Readings are taken at an interval of one hour till the sun sets in the evening. The solar intensity, wind velocity, wet bulb temperature, dry bulb temperature, basin plate temperature, basin water temperature, glass temperature, vapor temperature and distilled water productivity are measured every 1 hour. The same measurement process is repeated by varying the water level inside the solar still to 3 cm and 5 cm.
5. Heat Transfer
The energy balance equations at various portions of the solar still [20] [21] are described as follows:
The rate of energy gained by the glass and the rate of energy gained from the water surface to the glass by radiation, convection and evaporation are equal to the rate of energy lost to air.
5.1. Outside Surface of Glass Cover
Rate of energy received from inner glass cover by conduction = Rate of energy lost to the ambient by convection and radiation
(1)
Equation (1) for the energy balance equation of the outside surface glass cover becomes,
(2)
By simplifying Equation (2), the glass outside temperature
is written as,
(3)
5.2. Inner Surface of Glass Cover
Rate of energy absorbed from solar radiation + Rate of energy received from water mass by convection, evaporation and radiation = Rate of energy transferred to glass outside surface
(4)
Equation (3)-(4) for the energy balance equation of the inside surface of the glass becomes,
(5)
By substituting Equation (3) into Equation (5), the temperature of the inside surface of the glass is written as,
(6)
By rearranging Equation (6),
(7)
where
By simplifying the above equation, we get,
(8)
5.3. Basin Liner
Rate of energy absorbed from solar radiation + Rate of energy lost to water mass by convection = Rate of energy lost to the ambient by conduction and convection
(9)
Equation (9) becomes
(10)
After simplifying Equation (10), the basin temperature of the solar still is,
(11)
where
5.4. Water Mass
Rate of energy absorbed from solar radiation + Rate of energy received from basin liner by convection + Rate of energy received from external devices = Rate of energy stored + Rate of energy lost to the glass inner surface by convection, evaporation and radiation
(12)
The energy balance equation of the water mass in the solar still becomes,
(13)
By substituting values of Tgi and Tb from Equations (8) and (11) into Equation (13),
(14)
(15)
where
,
and
In the case of a passive solar still, additional thermal energy becomes zero,
.
The above equations can be expressed in simplified form as follows:
(16)
where,
and
The following assumptions have been made to find the approximate analytical solution:
1) The time interval ∆t (0 < t < ∆t) is small;
2) The value of “a” is constant during the time interval ∆t;
3) The function “ƒ(t)” is constant for the time interval between 0 and t;
4) i.e.
.
By using the boundary condition,
, the solution for above first order differential equation is
(17)
The hourly yield is given as follows,
(18)
The total daily yield is given as follows,
(19)
5.5. Overall Thermal Efficiency of the Olar Still
The instantaneous thermal efficiency is calculated as
(20)
The overall thermal efficiency of passive solar still is
(21)
where, mew = mass of distillate collected in kg/h; A = Area of the basin in m2; I(t) = Solar radiation with respect to time W/m2; and L = Latent heat of vaporization in kJ/kg.
6. Results and Discussion
The experimental study on single basin single slope passive solar still with black coated copper sheet is conducted at Nazareth, Tamil Nadu, India (Latitude: 9˚N; Longitude: 78˚E and an altitude of 9 m above sea level). Experiments are carried out from 7:00 h to 17:00 h during clear sunny days. During the experiments, hourly variation of solar intensity, wind velocity, ambient, glass inside, glass outside, vapor, water, basin temperature and distillate yield are recorded. The average solar radiation is 755 W/m2. It has been observed that solar radiation varies parabolically with time, with its maximum value at 13.00 h. In the morning, the temperature of water is low; therefore it needs high energy to change its phase from saturated liquid to saturated vapor phase. The ambient temperature reached the minimum value in the early morning and it is gradually increased after sun rise. The results show that temperature and required heat are inversely proportional. In the early afternoon the temperature of water reaches the maximum so it needs less heat to vaporize and vice versa in the late afternoon. The average wind velocity on test day was 3 m/sec. The wind velocity is an important parameter which affects the glass outside temperature of the solar still due to which condensation at higher rate is attained and ultimately it is expected to give higher yield.
The hourly variation of measured parameters such as ambient temperature, basin, water temperature, vapour temperature, glass inside temperature and glass outside temperature are shown in Figure 3. The maximum water temperature of 77˚C and glass inside temperature of 70˚C was observed at 14:00 h.
The hourly variation of yield is shown in Figure 4. Copper has higher thermal conductivity and coating of black paint absorbs all the incident radiation so that the basin water temperature is increased and hence the increased in the evaporation rate. It was found that the higher distilled water yield was obtained with higher basin water temperature. It revealed that there is proportionality increase in distilled water production with respect to the basin water temperature.
Figure 5 shows the yield of solar stills in summer and winter for the PSS-Cu, B. It is observed that the maximum yields of 0.19 kg for the PSS-Cu, B was obtained at 13:00 h during the summer. This was due to the higher water temperatures in the basin of the solar still which increases the evaporation rate and in turn increases yield. It also shows that maximum yield of 0.13 kg for the PSS-Cu, B at 12:00-13:00 h during the winter. Daily climatic conditions have resulted in slight variations in water evaporation rates under low radiation conditions. By comparing the operation of the solar stills in the summer and winter, lower yield was observed for the solar still during the winter as a result of lower radiation and ambient temperatures. However, the percentage increase in yield was as high as it was during the summer. This can be attributed to less thermal losses under low temperature operating conditions resulting in an increase in yield.
Figure 6 shows the comparison of convective heat transfer coefficient in summer and winter. The convective heat transfer coefficient is not varied drastically and it is dependent on wind velocity. The maximum convective heat transfer coefficients of 2.91 W/m2 K and 2.41 W/m2 K for the PSS-Cu, B were obtained for the summer and winter respectively. The convective heat transfer coefficient is PSS-Cu, B is lower due to the vacuum present inside the solar still. The temperature difference between (ΔT) evaporative surface (water) and condensing surface (glass inside surface), which plays significant role in raising the yield and in evaporative heat transfer coefficient as well.
Figure 3. Hourly variation of different temperature.
Figure 4. Hourly variation of hourly yield.
Figure 5. Comparison of hourly yield in summer and winter.
Figure 6. Convective heat transfer coefficient in summer and winter.
Figure 7 shows the comparison of evaporative heat transfer coefficient in summer and winter. The evaporative heat transfer coefficient increases with time and reaches a maximum at 14:00 h during the summer and 15:00 h during the winter prior to decreasing with time. The maximum value of the evaporative heat transfer coefficient measured during the summer and the winter was 75.13 W/m2 K and 42.52 W/m2 K respectively. The evaporative heat transfer coefficient is much higher than the convective and radiative heat transfer coefficients. The evaporative heat transfer coefficient is temperature dependent and is higher in the summer due to higher water temperature in the basin and lower in the winter. This is due to the higher thermal conductivity of copper sheet.
Figure 8 shows the comparison of radiative heat transfer coefficient in summer and winter. The radiative heat transfer coefficient did not vary much compared to the evaporative heat transfer coefficient and it was independent of temperature under normal operating conditions (<80˚C). The maximum radiative heat transfer coefficients of 7.96 W/m2 K and 6.91 W/m2 K were measured during the summer and winter respectively. Further it is also noticed that the convective and radiative heat transfer coefficients remain independent of the operating temperature.
Figure 9 shows the comparison of total heat transfer coefficient insummer and winter. The maximum total heat transfer coefficients were 110.62 W/m2 K and 65.43 W/m2 K in the summer and winter respectively, which can be attributed to the higher operating water temperature in the basin. The total heat transfer coefficient was much lower. The maximum values of total heat transfer coefficient 85.01 W/m2K and 52.10 W/m2K in the summer and winter respectively.
Figure 10 shows the hourly variation of the convective, radiative, evaporative and total heat transfer coefficients. It also shows that the evaporative and total heat transfer coefficients were much higher than the convective and radiative heat transfer coefficients. The evaporative heat transfer coefficient increased with time until 14:00 h and then its value decreased with time. The maximum values of the total and evaporative heat transfer coefficients were 58.6 W/m2K and 50.4 W/m2 K respectively. The higher value of the total and evaporative heat transfer coefficients are due to the increase in the basin water temperature. The convective heat transfer coefficient did not vary much and mainly depended on the wind velocity.
The maximum thermal efficiencies achieved by passive solar still with black coated copper solar still is 40% as shown in Figure 11 and 22% for the summer and winter respectively. The yield of the still was higher when compared to the still made up of galvanised iron sheet in both the summer and winter.
Figure 7. Evaporative heat transfer coefficient in summer and winter.
Figure 8. Radiative heat transfer coefficient in summer and winter.
Figure 9. Total heat transfer coefficient in summer and winter.
Figure 10. Hourly variation of heat transfer coefficients.
Figure 11. Hourly variation of the thermal efficiency.
7. Conclusion
The performance of single basin passive solar still is analysed on normal sunny days. The radiative heat transfer (Qr), convective heat transfer (Qc) and evaporative heat transfer (Qe) for internal and external heat transfer rates inside the still are predicted. The instantaneous efficiency and overall efficiency are calculated for all the studies. Heat transfer coefficients predicted under internal and external heat transfer modes for still performance. The thermal conductivity and dynamic viscosity increase with respect to time and possess almost the same trend. The variation of density of water is found for the still under these experiments and it is observed as 1.094133 to 1.1025812 kg∙m−3. It is observed that the density decreases with respect to increase in water temperature and it starts to increase with the decrease in water temperature. The amount of daily distillate yield collected is 1.09 kg/day. The instantaneous efficiency is increased according to the time.
Acknowledgements
We would like to acknowledge our sincere thanks to the management of our colleges, our Colleagues and our family members who have supported and helped us different stages of this Research work.
Nomenclature
A: area, m2
Cp: specific heat, J/kgK
I(t): solar intensity on an inclined surface, W/m2
Id: diffuse radiation on a horizontal surface, W/m2
Ig: global radiation on a horizontal surface, W/m2
p: partial pressure, N/m2
Q: heat transfer rate, W
q: heat flux, W/m2
T: temperature, ˚C
dt: time interval, s
dT: temperature difference, ˚C
h: heat transfer co-efficient, W/mK
hfg: enthalpy of evaporation at Tw, J/kg
mc: condensate, kg/hr
m: mass, kg/s
U: side heat loss co-efficient from basin to ambient, W/m2 K
V: wind velocity, m/s
K: thermal conductivity, W/mK
L: latent heat of vaporization, kJ/kg
hcb: basin liner convection heat transfer coefficient, W/m2 K
hb: basin liner overall heat transfer coefficient, W/m2 K
hcg: glass cover convection heat transfer coefficient, W/m2 K
hcw convective heat loss coefficient from water surface, W/m2 K
hew: heat loss coefficient by evaporation from water, W/m2 K
hrb: basin liner radiative heat transfer coefficient, W/m2 K
hrg: glass cover radiative heat transfer coefficient, W/m2 K
hrw: basin water radiative heat transfer coefficient, W/m2 K
htg: total glass heat transfer loss coefficient, W/m2K
hw: convective heat transfer coefficient from basin to water W/m2 K
htw: total water surface heat transfer loss coefficient, W/m2 K
kins: insulation thermal conductivity, W/mK
Lins: insulation thickness, m
M: total mass productivity/day, kg/day
MCw: water heat capacity rate of water per unit area, J/m2K
Pg: glass saturated partial pressure, N/m2
Pw: water saturated partial pressure, N/m2
qg: rate of total energy from the glass cover, W/m2
Subscripts
a: ambient
b: basin
c: convective
e: evaporative
g: glass
r: radiative
w: water
eq: equivalent loss side loss
d: daily
t: time
Symbols
ε: emissivity
α: absorptivity
ρ: density, kg/m3
σ: Stefan-Boltzmann constant, 5.6697 × 10−8 W/m2 K4
η: efficiency, %
β: inclination angle
θi: incident angle on an inclined surface
θh: incidentangle on a horizontal surface