1

Modeling and Experimental Study for Winter Performance of an Earth to Air

Heat Exchanger: An Alternative Energy Source for Greenhouse
1
M.K.Ghosal, Sujata Nayak2, G.N.Tiwari3 and N.Sahoo4
1, 4
College of Agricultural Engineering and Technology,
Orissa University of Agriculture and Technology,
Bhubaneswar-751003, India
2, 3
Centre for Energy Studies, Indian Institute of Technology Delhi,
Hauz Khas, New Delhi-11 00 16, India.
1
Corresponding author email: mkghosal1@rediffmail.com / mkghosal@yahoo.com

ABSTRACT

A complete numerical model has been developed to investigate the potential of using the stored
thermal energy of ground for space heating with the help of an earth to air heat exchanger
(EAHE) system integrated with the greenhouse located in the premises of IIT, Delhi, India. The
analysis was based on quasi-steady state condition. Experiments were conducted extensively
during winter period from November 2002 to March 2003, but the model, developed, was
validated against the clear and sunny days. The performance of the system was evaluated in
terms of total heating potential obtained from EAHE, coefficient of performance (COP) and
thermal load leveling. The heating potential of the system has also been standardized by the
characteristic curve for greenhouse similar to that of flat plate collector. Temperatures of
greenhouse air were found to be on an average 7-8 0C more than the same greenhouse when
operating without earth air heat exchanger. The temperature fluctuations of greenhouse air were
also less when operated with EAHE as compared to without EAHE. Predicted and measured
values of greenhouse air temperature in the model developed, exhibited fair agreement.

Keywords: Greenhouse, solar energy, earth air heat exchanger, coefficient of performance
(COP), thermal load leveling, modeling, India.

1. INTRODUCTION

Heating of greenhouse is one of the most energy consuming activities during winter periods.
Lack of heating has adverse effects on the yield, cultivation time, quality and quantity of the
products in the greenhouse (Santamouris et al. 1994a). But studies on greenhouse heating
strategies have shown that the cost of heating even exceeds 30% of the overall operational cost
of the greenhouse (Coffin, 1985). Due to high relative cost of energy, only a small number of
greenhouse owners can afford to the use of auxiliary heating systems. The use of low-cost and
alternative heating system is therefore of primary importance for a greenhouse to provide
optimum indoor conditions during winter months.
Efforts to decrease energy consumption have directed the researchers to use alternative energy
sources for heating of greenhouse. Several types of passive solar systems and techniques have
_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 2

been proposed and used (Santamouris et al. 1996, Barral et al. 1999, Bargach et al. 2000, Kumari
et al. 2007) for the substitution of conventional fuels with solar energy as a low cost technology.
As solar energy is the ideal source of energy for space heating particularly in the northern
hemisphere where it is available sufficiently. Also it is practically inexhaustible and its use does
not result in pollution. Efforts to harness solar energy have been accelerated during the last
decade as world demand for energy grows (Saravia et al. 1997). In all the heat collection
systems, the basic strategy is to reduce the heat losses and at the same time to store the surplus
energy for use during the energy shortage period. Continuous research in this area and several
successful demonstration projects has resulted in rapid advancements and commercialization of
these technologies with satisfactory results.
The greenhouses, which utilize solar energy for heat purposes, are equipped with heat collection
systems integrated into the cell of the greenhouse. The important heat storage mediums are
namely water, latent heat materials, rock bed, buried pipes etc. Out of these storage mediums,
buried pipes systems have gained increasing acceptance for better and easy exploitation of
thermal energy from the ground (Santamouris et al. 1994b).
In buried pipe systems, the nearly constant and stored thermal energy of earth at a certain depth
is usually extracted with the help of an arrangement called earth air heat exchanger (EAHE). The
stored thermal energy and thereby the earth’s surface and sub surface temperatures at any given
location, are determined by the balance between the solar energy absorbed at the surface and heat
losses by the outgoing long wave radiations and convective heat exchange with ambient air mass
(Sodha et al. 1981). An earth air heat exchanger system herein is defined (Puri, 1987) as the
study of heat transfer between soil, tubes and air flowing through the tube when the tubes are
placed below the ground surface at a certain depth where temperature of soil remains nearly
constant throughout the year. As air travels the length of the tube, it gets heated in the winter
period and gets cooled during the summer period resulting in the space conditioning due to its
entry into the enclosed space. Earth air heat exchanger system has therefore the potential of
being used throughout the year. Hence considering the importance of EAHE as a simple,
inexpensive and alternative source of energy, the system has been incorporated in IIT, model
greenhouse, New Delhi, during the whole winter period with a view to study its thermal
performance for heating of the greenhouse through a model, developed, in the composite climate
of India.
2. EXPERIMENTAL PROCEDURES

The EAHE under study was used in the greenhouse located in the premises of IIT, Delhi, India.
The climate of the place is composite i.e., it remains hot dry for five months, warm and humid
for three months, moderate for one month and cold for three months. The absolute minimum
temperature of ambient air during winter period is close to 4 0C while mean minimum is close to
9 0C. The greenhouse combined with EAHE was of even span type of greenhouse with floor area
6m x 4m. and was oriented from east to west direction. The EAHE was installed outside in west
side of the greenhouse. Total length and diameter of buried pipes used were 39m and 0.06m
respectively. EAHE also consisted of PVC pipes buried under bare surface at the depth of 1m in
a serpentine manner with 8 nos. of turns. The blower was attached in the suction end of the
EAHE. The suction and delivery ends of EAHE were placed in the southwest and northwest
corners of the greenhouse for allowing uniform mixing of air. The isometric view of
_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 3

experimental greenhouse integrated with EAHE is shown in Fig. 1a. Experiments were

N
Ground
surface W E
Ground S

North wall

Outlet of burried pipe

Burried pipe
(dia=6 cm)
South wall
Blower
(inlet to burried pipe)

Exhaust fan

All Dimensions in cm

Figure 1a. Isometric view of even span greenhouse integrated with EAHE arrangement

conducted continuously for two days in a week in clear and sunny days from Nov’2002 to
March’2003 with 1st day without any heating arrangement and 2nd day with EAHE system.
However the experimental validation was done for typical date (clear sunny day) of observations
i.e., on 23-01-03 for greenhouse with EAHE, since January is the coldest month for Delhi.
Hourly observations of solar radiation and temperatures of air for ambient condition, greenhouse
enclosure, suction end and delivery end were recorded during the experimentation with the help
of calibrated solarimeter and mercury thermometer, respectively.

3. THERMAL ANALYSIS

The energy balance equations for various components of greenhouse combined with earth to air
heat exchanger can be written on the basis of following assumptions:
(i) Analysis is based on quasi steady state conditions,
(ii) There is no radiative heat exchange between the walls and roofs of greenhouse, due to
negligible temperature differences,
(iii) Flow of air is uniform along the length of buried pipes,
(iv) Heat flow is one-dimensional.
_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 4

Energy balance equations for north wall, floor and room air of greenhouse are as follows:
a) North wall
α n (1 − rn ) Fn (1 − r ){∑ Ai I iτ i } = hnr (T y = 0 − Tr ) An + hna (T y = 0 − Ta ) An (1)
b) Floor
α g (1 − rg )(1 − Fn )(1 − r ){∑ Ai I iτ i } = hgr (T x =0
− Tr ) Ag − hg∞ (T x =0
− T∞ ) Ag (2)
At larger depths, the temperature of ground is assumed to be equal to ambient air temperature,
T∞ = Ta , then Eq. (2) becomes
α g (1 − rg )(1 − Fn )(1 − r ){∑ Ai I iτ i } = hgr (T x =0
− Tr ) Ag + hg∞ (T x =0
− Ta ) Ag (3)
c) Greenhouse air
hnr (T y = 0 − Tr ) An + hgr (T x=0
− Tr ) Ag + Q& u = ∑ AiU i (Tr − Ta )
dTr (4)
+ 0.33NV (Tr − Ta ) + M aCa
dt
&
The term i.e., Qu in Eq. (4) is the useful thermal energy obtained from EAHE and is expressed
by the equation, Q& u = FR m& a C a (T0 − T fi ) (5)
2 π r1 hgf
− L′
where FR = 1 − e m& a C a
. Now eliminating T y =0
from Eq. (1) and after rearrangement,
I effN
hnr (T y =0
− Tr ) = F1 − U n (Tr − Ta ) (6)
An
hnr (h )(h )
where I effN = α n (1 − rn ) Fn (1 − r )(∑ Ai I iτ i ) , F1 = and U n = nr na
hnr + hna (hnr + hna )
Similarly eliminating T x =0
from Eq. (3) and after rearrangement,
I effF
hgr (T x=0
− Tr ) = F2 − U g (Tr − Ta ) (7)
Ag
hgr (hgr )(hg∞ )
where I effF = α g (1 − rg )(1 − Fn )(1 − r )(∑ Ai I iτ i ) ; F2 = ; Ug = and T fi = Tr .
hgr + hg∞ (hgr + hg∞ )
Now substituting Eqs. (6) and (7) in Eq. (4) and simplifying, Eq. (4) can be written in the
dTr
following first order differential equation, + aTr = B (t ) (8)
dt
F (t ) + (UA) eff Ta a1
where B (t ) = and a = ; F (t ) = F1I effN + F2 I effF + FR m& aCaTo
M aCa M a Ca
(UA)eff = U n An + U g Ag + 0.33NV + (∑ AiU i ) ; a1 = U n An + U g Ag + 0.33 NV + ( ∑ AiU i ) + FR m& aCa
( ∑ Ai I iτ i ) = ( Ae I eτ e + Aww I wwτ ww + Asr I srτ sr + Anr I nrτ nr + As I sτ s )
(∑ AiU i ) = ( AeU e + AwwU ww + AsrU sr + AnrU nr + AsU s ) ; U e = U ww = U sr = U nr = U s = U

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 5

Ln 1 Lg 1 L 1 1 1 −1 1 1
hna = [ + ]−1 , hg∞ = [ ]−1 , U n = [ + n + ]−1 , U g = [ + ] ; U = [ + ]−1 ,
K n h0 Kg hi K n h0 hgr hg∞ hi h0
hnr = hgr = hgf = hi . The analytical solution of Eq. (8) can be written as
B (t )
Tr = (1 − e −at ) + Tro e −at (9)
a
where, Tro is the greenhouse air temperature at t = 0 and B (t ) is the average of B (t ) for the
time interval 0 and t , and a is constant during the time. From Eq. (9), the temperature of air
inside greenhouse, combined with earth air heat exchanger can be determined for analysis.

3.1 Instantaneous Thermal Efficiency (η i ) Characteristic Curves for Greenhouse

The instantaneous thermal efficiency (η i ) is defined as the ratio of thermal energy used in raising
temperature of greenhouse air from Tro to Tr to input energy and is expressed as
M a C a (Tr − Tro )
ηi = (10)
t (∑ Ai I i )
Putting the expressions of Tr , B (t ) , F (t ) , (UA) eff , a and a1 in Eq. (10) and simplifying, the
Tro − Teff
instantaneous thermal efficiency becomes, η i = F ′ [(ατ ) eff − U eff ] (11)
I
1 − e − at F1 I effN + F2 I effF a
Where F ′ = ; (ατ ) eff = ; U eff = 1 ; I = ∑ I i ; at is dimensionless;
at ∑ Ai I i ∑ Ai
FR m& a C a To + (UA) eff Ta
Teff = and F ′ is the efficiency factor for greenhouse (dimensionless). Eq.
a1
(11) is the function of design and climatic parameters and is similar to the characteristic equation
for flat plate collector (Duffie and Beckman, 1991). This equation is helpful for comparison and
standardization of various heating methods inside the greenhouse.

3.2 Computational Procedure and Input Parameters

The energy balance equations derived for greenhouse with EAHE have been solved with the help
of a computer program based on Matlab software. The design and operating parameters given in
Table-1 have been used as input parameters for the mathematical model developed. The
closeness of predicted and experimental values has been presented with coefficient of correlation
(cr) and root mean square of percent deviation (er). Solar radiation falling on different walls and
roofs of greenhouse was calculated with the help of Liu and Jordan (1962) formula by using the
beam and diffuse components of solar radiation incident on the horizontal surface. The heat
removal factor for EAHE has been calculated from steady state energy mechanism as shown in
Fig. 1b and as per Eq. (5). The mass flow rate of the circulating air was kept constant with 100
kg/hour.
_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 T T+dT 6
L

ma To
Inlet Tfi Outlet Tfo

Ground

X
dx
T(x) T(x)+dT(x)

Figure 1b. Energy exchange between ground and flowing air in elementary segment of the buried pipe

The performance of EAHE has been evaluated in terms of thermal load leveling (TLL) (Singh
and Tiwari, 2000) heating potential and COP as per the following
Tr ,max − Tr ,min
expressions: TLL = ; Qh = ∑ m& a C a (Td − Tsc ) Δt and
Tr ,max + Tr ,min
output energy
COP = . Thermal load leveling gives an idea about the
Energy spent to get output energy
fluctuations of air temperature inside the greenhouse. The less the fluctuations, the better is the
environment for plants inside the greenhouse. In winter, TLL should have lower values by
incorporating heating method due to the increase of (Tr , max + Tr , min ) as well as decrease of
(Tr , max − Tr , min ) as compared to TLL without heating arrangement. The temperatures of ground
i.e., To were recorded with the help of data logger through the thermocouples located at the depth
of 1.0m under EAHE arrangement.

Table 1. Input parameters used for computations

Parameters Values Parameters Values Parameters Values
2 20
Ae 8.3 m hi 2.8 W/m C N 1-300
Af 24.0 m2 ho 5.7 W/m2 0C r1 0.03m
An 12.0 m2 hna 1.9 W/m2 0C U 1.8 W/m2 0C
As 12.0 m2 hgr 5.7 W/m2 0C v 0.5-1.5 m/s
2 20
Anr 13.8 m hnr 5.7 W/m C V 60 m3
Asr 13.8 m2 Kn 0.84 W/m0C rg 0.2
Aww 10.0 m2 Kg 0.52 W/m0C rn 0.2
Ca 1012 J/kg 0C L′ 39m αg 0.4
Fn 0.09-0.15 Lg 1m αn 0.6
FR 0.64 m& a 0.02 kg/s τ 0.5
20
hgf 2.8 W/m C Ma 72 kg

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 7

4. RESULTS AND DISCUSSION

The hourly variations of temperature for ambient air, greenhouse air when operating with EAHE
for typical winter day (23-01-2003) and without EAHE (22-01-2003) have been presented in Fig.
2. From the figure, it is seen that the minimum as well as maximum temperatures for ambient air,
greenhouse air with EAHE and without EAHE varied between 5-20 0C, 13-24 0C and 5-28 0C
respectively indicating the less fluctuations of air temperature inside the greenhouse with EAHE
as compared to ambient air and greenhouse air without EAHE. This is due to the addition of
supplementary thermal energy to the greenhouse by EAHE arrangement. The temperature of
ground on the above day at the depth (1m) in which the EAHE system was installed was
recorded to be about 23 0C. The predicted values of greenhouse air have been validated with
experimental values for the above typical day (23-01-2003) and they showed fair agreement with
(cr ) as 0.98 and (er ) as 6.93. By examining closely the daily temperature profiles of
greenhouse air from Fig. 3, it was found that the delivery temperatures of EAHE were 7-9 0C
higher than the suction temperatures from 7 pm to 9 am (heating of greenhouse air). In between 9
am and 10 am both delivery and suction temperatures became almost same (with zero heating
potential obtained from EAHE) and from 10 am to 6 pm, the delivery temperatures were 3-6 0C
less than the suction temperatures (cooling of greenhouse air). As a result of which, the
temperatures of greenhouse air were maintained in the range of 13-24 0C for creating healthy
environment for the growth of plants during winter period. Also it was observed that the
temperatures of greenhouse air were around 7-8 0C higher than those of greenhouse air without
EAHE system for thermal heating during nighttime.
After knowing the suction and delivery temperatures of EAHE as well as mass flow rate, the
diurnal variations of total heating potential obtained from the system for the typical day in the
winter months were calculated and have been shown in Fig. 4. From the figure, it is evident that
the air in the greenhouse was heated during off-sunshine hours and cooled during sunshine hours
(10 am to 6 pm) causing the reduction of its undesirable rise of temperature in the latter period.
Similarly the total heating potentials obtained from EAHE for a typical day in each winter
months have been computed and presented in Fig. 5. From the results, it is seen that the heating
potentials obtained from EAHE were higher in the month of January followed by December,
February, November and March. The higher value of heating potential in January (coldest
month) is due to the more differences of temperature in suction and delivery ends. The
coefficient of performance determined for typical day in each month has also been discussed in
Fig. 6 to know the applicability of the system. The values of coefficient of performance were
highest in the month of January (2.8), followed by December (2.45), February (2.1) and
November (1.94). However, in the month of March its value was below the dashed line (value
less than 1) indicating the system needs to be discontinued during this month. The values of
thermal load leveling achieved for typical days in each month have been calculated and
presented in Fig. 7 in order to know the efficacy of system during the study. From the computed
results, it is seen that the values of TLL in each month for greenhouse with EAHE were lower
than those without EAHE proving the former to be more effective for reducing the daily swings
of temperature of air in greenhouse. After computation of room air temperature, instantaneous
thermal efficiency characteristic curve was evaluated from Eq. (11) for the greenhouse with and

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 8

without EAHE. The equation represents the equation of straight line between efficiency in Y-
axis and (Tro − Teff ) / I in X-axis. The intercept F ′ (ατ )eff is the gain term where as slope of
gradient − ( F ′U eff ) is the loss factor. Instantaneous thermal efficiency curves standardize and
compare the heating potential of different heating methods. The slope of gradient (m) in the
curve represents the magnitude of various thermal losses from enclosed air to ambient where as
the intercept (gain term) refers to the thermal energy rise of the enclosure (greenhouse) air
particularly by incoming solar radiation and by other auxiliary heating sources. For heating of an
enclosure, the loss factor should be as minimum as possible and the gain term should be as
maximum as practicable. From Fig. 8, it is evident that gain term was more in case of greenhouse
with EAHE as compared to without EAHE. On the contrary, the loss factor was less in
greenhouse without EAHE. The reason for the higher values of loss factor in case of greenhouse
with EAHE is due to the more temperature differences between greenhouse air and ambient air
resulting in more heat transfer in that condition. However, overall heating effect was more in
greenhouse with EAHE as compared to without EAHE due to additional thermal energy gain in
the former by the EAHE arrangement. Also theoretical and experimental characteristic curves
showed good agreement in both the conditions of experiment.

30
Date:23-01-2003 Tr
27 Tr=greenhouse temp. with EAHE Tr1
Tr1=predicted greenhouse temp. Tr2
with EAHE Ta
24 Tr2=greenhouse temp. without EAHE
Ta=ambient temp.
21
cr=0.98
Temperature 0C

18 er=6.93

15

12

0
1 3 5 7 9 11 13 15 17 19 21 23
Time of day in hour

Figure 2. Hourly variations of greenhouse air temperature for typical winter day by EAHE

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 30
Date: 23-01-2003 Td
27 Td=delivery temp. Ts
Ts=suction temp.
Tr=greenhouse air temp.
Tr
Tr1
9
24 (with EAHE 23-01-03) Ta
Tr1=greenhouse air temp.
21 (without EAHE 22-01-03)
Ta=ambient air temp.

Temperature 0C
18

15

12

0
1 3 5 7 9 11 13 15 17 19 21 23
Time of day in hour

Figure 3. Hourly variations of suction, delivery, greenhouse air (with and without EAHE)
and ambient air temperatures during experimentation

5
x 10
12
11 Date:23-01-2003
10
9
8
7
6
5
Thermal energy (J)

4
3
2
1
0
-1
-2
-3
-4
-5
-6
-7
-8
1 3 5 7 9 11 13 15 17 19 21 23
Time of day in hour

Figure 4. Hourly variations of heating potential by EAHE for a typical winter day

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 10

140 Total heating potential

Thermal energy (MJ)

120
100
80
60
40
20
0
Nov Dec Jan Fab March
Winter months

Figure 5. Monthly variations of total heating potential obtained from EAHE during experimentations

3 COP

2.5
2
COP

1.5
1
0.5
0
Nov Dec Jan Feb March
Winter months

Figure 6. Monthly variations of coefficient of performance (COP) during experimentations

0.7 With EAHE

0.6 Without EAHE
0.5
0.4
TLL

0.3
0.2
0.1
0
Nov Dec Jan Feb March
Winter months

Figure 7. Monthly variations of thermal load leveling (TLL) during experimentations

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 11

0.5

22-01-03 (without EAHE)

23-01-03(with EAHE) exptl ηi(without EAHE)
0.45
predictedηi(without EAHE)

Instantaneous thermal efficiency(ηi) (decimal)

exptl ηi(with EAHE)
(Ueff)=9.46, (ατ)eff=0.38 predictedηi(with EAHE)
0.4

0.35

0.3 (Ueff)=5.67, (ατ)eff=0.24

0.25

0.2

0.15

0.1
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
(Tro-Teff) / I (oC/w/m2)

Figure 8. Instantaneous thermal efficiency characteristic curve for greenhouse with and without EAHE

5. CONCLUSIONS

From the above results, the main conclusions for the present study are as follows:
(i) There occurs 7 –8 0C rise of temperatures for greenhouse air during winter period due to
the incorporation of EAHE as compared to without EAHE,
(ii) Relative fluctuations of temperature for greenhouse air are less in EAHE arrangement
than without that system
(iii)The predicted and experimental temperatures of greenhouse air in the model
developed, with EAHE arrangement exhibit fair agreement.
(iv)The computed and experimental thermal efficiency characteristic curves for
greenhouse compare well with each other.

6. ACKNOWLEDGEMENTS

The authors are grateful to Indian Council of Agricultural Research, Krishi Bhawan, New Delhi,
India for providing financial support to carry out the above experiment.

7. REFERENCES

Barral, J.R., P.D. Galimberti, A. Barone and A.L. Miguel. 1999. Integrated thermal
improvements for greenhouse cultivation in the central part of Argentina. Solar Energy,
67(1-3): 111-118.

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 12

Bargach, M.N., R. Tadili, A.S. Dahman and M. Boukallouch. 2000. Survey of thermal
performances of a solar system used for the heating of agricultural greenhouses in
Morocco. Renewable Energy, 20: 415-433.
Coffin W.1985. Design and testing of a cold climate greenhouse. ISES Congress, Montreal.
p.611.
Duffe, J.A., W.A. Beckman. 1991. Solar Engineering of Thermal Processes, 2nd edition, New
York: John Wiley and Sons.
Kumari, N., G. N.Tiwari, and M.Sodha. 2007. Performance evaluation of greenhouse having
passive or active heating in different climatic zones of India. Agricultural Engineering
International: the CIGR Ejournal. Manuscript EE 06 011. Vol. IX.
Liu B.Y.H. and R.C. Jordan 1962. Daily insolation on surfaces tilted towards the equator.
ASHRAE Journal. 3(10): 53.
Puri V.M. 1987. Earth tube heat exchanger performance correlation using boundary
element method. Transactions of the ASAE. 30(2): 514-520.
Santamouris M., C.A. Balaras, E. Dascalaki and M. Vallindras. 1994a. Passive solar agricultural
greenhouse: a worldwide classification evaluation of technologies and systems used for
heating purposes. Solar Energy. 53(5): 411-426.
Santamouris M., A. Argiriou and M. Vallindras. 1994b. Design and operation of a low energy
consumption passive solar agricultural greenhouse. Solar Energy. 52 (5): 371- 378.
Santamouris M., G. Mihalakakou, C.A. Balaras, J.O. Lewis, M. Vallindras and A. Argirious.
1996. Energy conservation in greenhouses with buried pipes. Energy. 21(5): 353-360.
Saravia L, R. Echazu, C. Cadena, M. Condori, C. Cabanillas, A. Iriarte and S. Bistoni.1997.
Greenhouse solar heating in the Argentinian northwest. Renewable Energy. 11(1): 119-
128.
Singh R.D. and G.N. Tiwari. 2000. Thermal heating of controlled environment greenhouse: a
transient analysis. Energy Conversion and Management. 41: 505-522.
Sodha M.S., N.K. Bansal and A.K. Seth. 1981. Variation of ground temperature distribution.
Applied Energy. 8: 245.

8. NOMENCLATURE

A - Area, m2
Ca - Specific heat of air, J/kg 0C
Fn - Fraction of solar radiation falling on north wall, dimensionless, decimal
FR - Heat removal factor for EAHE from underground earth’s surface
hi - Heat transfer coefficient from greenhouse cover to inside
greenhouse air, W/m2 0C, ( 2.8 + 3.0 v ), (Duffie and Beckman, 1991)
ho - Heat transfer coefficient from greenhouse cover to ambient,
W/m2 0C, ( 5.7 + 3.8 v ),(Duffie and Beckman, 1991)
hgf - Convective heat transfer coefficient from underground earth’s surface to
flowing air inside the buried pipes, W/m2 0C

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 13

hg∞ - Heat transfer coefficient from floor to larger depth of ground, W/m2 0C
hna - Heat transfer coefficient from north brick wall to ambient, W/m2 0C
hnr - Heat transfer coefficient from north wall to greenhouse air, W/m2 0C
hgr - Heat transfer coefficient from floor to greenhouse air, W/m2 0C
I - Solar radiation falling on inclined surface or greenhouse cover, W/m2
K - Thermal conductivity, W/m 0C
Kg - Thermal conductivity of ground, W/m 0C
L - Thickness, m
L′ - Total length of buried pipes (EAHE), m
m& a - Mass flow rate of air entering into the buried pipes, kg/s
Ma - Total mass of air in greenhouse enclosure, kg
N - Number of air changes per hour
Qh - Heating potential offered by EAHE for greenhouse air, J
Q& u - Useful thermal energy obtained from EAHE for greenhouse air, W
r - Reflectivity from greenhouse cover, dimensionless, decimal
rg - Reflectivity from greenhouse floor, dimensionless, decimal
rn - Reflectivity from north wall, dimensionless, decimal
r1 - Radius of buried pipe in EAHE, m
t - Time in second
Δt - Time interval in hour
T - Temperature, 0C
Td - Delivery temperature, 0C
To - Temperature of ground in which pipes are spread in EAHE, 0C
T fi - Temperature of inlet fluid or temperature at suction point, 0C for EAHE
Tsc - Suction temperature, 0C
U - Overall heat transfer coefficient for greenhouse cover, W/m2 0C
Ug - Overall heat transfer coefficient from greenhouse air to floor, W/m2 0C
(UA) - Overall heat loss from greenhouse, W/ 0C
v - Velocity of air, m/s
V - Volume of greenhouse, m3
Greek letters
α - Absorptivity, dimensionless
τ - Transmissivity, dimensionless
∞ - Infinity (at larger depth)
(ατ )eff -Effective transmittance-absorptance product for greenhouse
Subscript
a - Ambient
_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.
 14

e - East wall of greenhouse

g - Floor of greenhouse
i - Different walls and roofs of greenhouse
n - North wall
r - Greenhouse room
s - South wall
nr - North roof
sr - South roof
ww - West wall
eff - Effective

_____________________________________________________________________________
M.K.Ghosal, Sujata Nayak, G.N.Tiwari, N.Sahoo. “Modeling and Experimental Study for
Winter Performance of an Earth to Air Heat Exchanger”. Agricultural Engineering International:
the CIGR Ejournal. Manuscript EE 07 012. Vol. X. January, 2008.