 Open Access
 Total Downloads : 162
 Authors : G. N. Tiwari, Ibrahim M. AlHelal, Poonam Joshi, Ravi Agrihari
 Paper ID : IJERTV3IS090450
 Volume & Issue : Volume 03, Issue 09 (September 2014)
 Published (First Online): 20092014
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Active Heating of Floating Type Dome Biogas Plant
G. N. Tiwari1, Poonam Joshi1, Ravi Agrihari1
1Centre for Energy Studies,
Indian Institute of Technology Delhi, Hauz Khas, New Delhi 1100 16, India
Ibrahim M. AlHelal2
2Department of Agricultural Engineering College of Food & Agricultural Sciences,
King Saud University,

ox 2460, Riyadh 11451, Saudi Arabia
Abstract – In this paper, an attempt has been made to derive an expression for the slurry temperature of flat plate collectors integrated biogas plant. The flat plate collectors are connected in series and parallel combination to optimise the number of collectors to be in series (N) and parallel (m) for optimum slurry temperature. Effect of mass flow rate ( m f ),
length of heat exchanger (L), heat capacity of slurry (Ms Cs) and number of flat plate collectors (FPCs) on slurry temperature have been studies for climatic condition of Srinagar, India by using antifreeze liquid. It has been observed that the optimum slurry temperature (~ ) has been observed for N=8 (series) and m=5 (parallel) for 40 number of flat plate collectors. Exergy of active biogas plant has also been carried out.
Keywords:Biogas, solar heating, flat plate collectors

INTRODUCTION
Biogas is a mixture of methane (CH4) 5070% and carbon dioxide (CO2) 3050%. It has a calorific value of 2124 MJ/m3 and it is produced at 37 in the absence of oxygen from slurry obtained from 50% dung and 50% water. It is clean energy and environmental friendly which can be used for cooking [15], lighting [56], motor vehicles and small scale industries [710] to meet local energy demand. The temperature of slurry depends on local climatic conditions namely solar intensity and ambient air temperature. The slurry in digester requires to be heated for harsh cold climatic condition to obtain the optimum slurry temperature. Basically, there are two type of biogas plant namely fixed dome and floating dome type biogas plant. In both cases, the slurry in digester is heated either by active method or by hot charging for harsh cold climatic condition. In an active method, a combination of flat plate collectors (FPCs) can be integrated through a coil type heat exchanger placed inside digester of floating/fixed type.
volume 6 ton for hot charging the slurry to ferment the slurry inside digester, Dong and Lu [15]. It has also been found that 2540 m3 biogas is produced with increase of 11.2%. They also observed that there were 14.3% increase of pig manure energy transformation efficiency.
In this paper, we have analysed the slurry temperature of active biogas plant to optimise the number of flat plate collectors to be in series and parallel for a given total number of flat plate collector and heat capacity of slurry. Exergy analysis of active biogas plant has also been carried out.

THERMALMODELLING
In order to write energy balance, the following assumptions have been made:

The floating dome type biogas plant is considered for thermal modelling

The biogas plant is in quasi steady state condition

No stratification along the depth of the slurry in digester due to forced mode and the gas column due to
low heat capacity

Thermal heat capacity (MsCs) of the biogas plant materials are neglected.
Referring to Fig. 1, the energy balance equations during sunshine hours have been formulated as follows:
Sun
Series of
The number of flat plat collectors (FPCs) depends on heat capacity of slurry with local climatic condition. Further, series and parallel combination of collectors should also be optimized accordingly.
Tiwari [1113] have proposed design criteria for active heating of fixed dome biogas plant. Yuan [14] have observed that the heat available from solar flat plate collectors with an effective area of 2m2 and 8 hours operation can met heat demands of a 6m3digester for complete fermentation in slurry. Flat plate collectors with
Collectors
Inlet tank
Inlet pipe
Pump
Partition Wall
Gas outlet
Gas Holder Outlet Tank
Digester
Outlet
Pipe
an effective area of 100.8 m2 is sufficient to heat water of
Fig.1 Schematic view of active floating type biogas plant.
Metallic biogas holder/dome:
2
+ = 1 +
+ + 2 (1) Produced biogas:
1 = 3 (2) Slurry in digester:
= 3 + +
4
+ (3) where,
=
is the rate of evaporation from side of
dome exposed to ambient which can be neglected for simplification of modelling.
and
Fig.2 . Combination of NPVT FPC connected in m rows.
0.016
Referring to Fig. 3, the energy balance for forced
=
circulating water in coil type heat exchanger is given by
= 2 (5)
Equations (13) can be combined into single equation by
1
eliminating Tp and Tg as follows:
= +
x=0 x x+dx
= 1+2
(4) where ,
Tw
= 1+ 2 + 4 +
1
= ;
+
.
Ts
Tw TwdTw
= + ; =
1 3
and = + 2
1 + 3
where,
Fig.3 Elemental length dx of heat exchanger
integrated with digester.
2
= ;
+
1
1 1
= + 1 2 + 1
= ;
=
1
2
1
+
2
+
In order to solve Equation (5), the initial condition namely
is the rate of thermal energy available from NFPC connected in series and m is the number of rows of NFPC connected in series, Fig.2 . Therefore there will be ( Ã—
) FPC in active biogas plant.
=0 = can be used.
Now the solution of Equation (5) is as follows
= 1 21 +
21 (6)Further = = (the outlet of
heat exchanger will be the inlet to the collectors connected in series and parallel), then one gets.
= 1 21 +
21
=
+
or,
= 1 1 + 1 (7)where,
or,
( )
2 r1U L
1 m C
[ ( ) ] f f
The average water temperature over the length inside heat exchanger has been obtained from Equation (6) as
0( )
( )
1 ( )
=
or,
1
0
or,
+ = () (12)
where
= 1 1 1 +
1 1 (8)Further,
1
1
m C
the rate of heat transfer from flowing fluid inside the heat
1 m f f eff
exchanger to the slurry has been obtained as
<> = 21 (9)
= 1 1
() =
+
U L
+ + m C
qab
Following Tiwari (2002), the outlet fluid temperature at the end of N th collector connected in series can be expressed as follows:
and
f f eff
=
+ 1
=
+ m
f C f
eff
Or,
+
The solution of Equation(12) with initial condition i.e. (at t=0) = becomes as
= () 1 exp() + exp() (13)
= + 1 2 + 2
(10)
where,
For a given design and climatic parameters, the hourly slurry temperature () can be obtained from Equation (13).After hourly variation, the maximum slurry temperature (, ) for a particular day will be evaluated.
2 =
Now Equation(9) becomes as
Once, the maximum slurry temperature(, ) is known, the exergy of active biogas system can be obtained from the following equation
= 1 1
(, + 273
=
1 exp 1 1 2
= (, ,
+ 273 , + 273
2
1 1
+


Results and discussion

Design parameters of Table 1 and climatic data of Fig. 4 have been used to evaluate slurry temperature by using
=
(11)
Equation(13). The hourly variation of slurry temperature for different configuration of flat plate collectors (FPCs)
After substituting above equation in Equation(4), one has
has been shown in Fig. 5. It is clear that maximum slurry temperature for all FPCs connected in series (m=1 and
N=40) is about 150C. This is not the optimum temperature (~ 370C) for biogas production. Hence the hourly slurry temperature has been calculated for other configuration as (m=2 and N=20; m=4 and N=10; m=5 and N=8; m=8 and N=5). The results have also been shown in Fig.(69). It is seen from Fig 5 that the optimum slurry temperature is achieved for configuration of m=5 and N=8. This can be possible because of less thermal energy loss for 8 flat plate collectors are connected in series (N=8). In other combination thermal losses are significant.
Table 1:Design parameters of active biogas plant
600
Solar Intensity(W/m2)
500
400
300
200
100
5
4
Ambient Temperature(0C)
3
2
1
0
1
2
3
3 5 8 10 13 15 18 20 23 25 28
Time (hr)
Parameters 
Values 
Parameter s 
Values 
Ac 
2 m2 
Ms 
2500 kg 
Cs 
4190 J/kg K 
K 
204 W/mK 
UL1 
3.56 W/ m2 0C 
FRC 
0.95 
hbf 
100 W/m2 0C 
ULC 
6 W/m2 0C 
Utc,a 
9.5 W/m2 0C 
L 
25 m 
r 1 
0.0125 m 
As 
8.5 m2 
r 2 
0.0175 m 
N 
40 
r 3 
0.625 m 
hc 
58 W/m2 0C 
Av 
10.3 m2 
hrps 
5.2W/m2 0C 
Ah 
8.5 m2 
hsa 
2.8 W/m2 0C 
Av 
4.5 m2 
h 1 
0.66 W/m2 
Ah 
0.9 m2 
h 2 
5 W/m2 
h 3 
1.32 W/m2 
h 4 
0.78 W/m2 
Fig 4. Hourly variation of I(t) and Ta for typical day of Srinagar.
40
Temperature of slurry (0C)
30
20
10
0
0 5 10 15 20 25
Time (hr)
(m=1,N=40) (m=2,N=20) (m=4,N=10) (m=5,N=8)
(m=8,N=5)
Fig.5. Hourly slurry temperature for different configurations.
35
Temperature of slurry (0C)
30
25
20
15
10
0 5 10 15 20 25
Length of pipe(m)
Fig 6. Effect of length on Ts,max.
41
40
Temperature of slurry (0C)
39 m=5,N=8
38
37
36
35
34
33
0.00 0.02 0.04 0.06 0.08 0.10
Mass flow rate (Kg/s)
Fig 7. Effect of mass flow rate on Ts,max .
80
70
Temperature of slurry (0C)
60
50
40
30
20
10
0
0 5000 10000 15000 20000
Mass of slurry (Kg)
Fig 8. Effect of mass on Ts,max.
80
70
60
Exergy (kWh)
50
40
30
20
10
0
0 5000 10000 15000 20000
Mass of slurry (Kg)
Fig 9. Effect of exergy with mass of slurry.
For the optimized number of flat plate collectors in series and parallel (m=5, parallel and N=8, series), parametric studies have been carried out. Fig. 6 shows the effect of heat exchanger length on maximum slurry temperature (Ts,max) and it can be seen that the variation of maximum slurry temperature (Ts,max) becomes insignificant after 25 m length of heat exchanger. It is not economical to have more length of heat exchanger due to copper materials.
Effect of mass flow rate of fluid on maximum slurry temperature (Ts,max) for other optimized parameters namely configuration and length of heat exchanger has been shown in Fig.7. It can be seen that there is not much variation in maximum slurry temperature (Ts,max) after mass flow rate of 0.04 kg/s. Hence the mass flow rate of 0.04 kg/s is optimum for a given other design parameters.
In Fig. 8, the variation of maximum slurry temperature (Ts,max) with different mass of slurry has been shown. The figure indicates that the slurry temperature is maximum at lower value of slurry mass which is not suitable for biogas production. The optimum temperature can be achieved at the mass of 2500 kg.
Equation (13) has been used to evaluate exergy of active biogas plant for different mass of slurry. The results have been shown in Fig. 9. The exergy of active biogas plant at optimum parameters is 40 kWh.
IV. REFERENCES

Bundit Limmeechokchai, and Saichit Chawana, The case of the improved cooking stove and the small biogas digester, Sustainable energy development strategies in the rural Thailand,Renewable and Sustainable Energy Reviews, Volume 11, Issue 5, pp. 818837, June 2007.

Matthew Landi, Benjamin K. Sovacool, and Jay Eidsness, Cooking with gas: Policy lessons from Rwanda's National Domestic Biogas Program (NDBP), Energy for Sustainable Development, Volume 17, Issue 4, pp 34735, August 2013.

Ershad Ullah Khan, Brijesh Mainali, Andrew Martin, and Semida Silveira, Technoeconomic analysis of small scale biogas based polygeneration systems: Bangladesh case study, Sustainable Energy Technologies and Assessments, Volume 7, pp 6878, September 2014.

Karthik Rajendran, Solmaz Aslanzadeh, Fredrik Johansson, and Mohammad J. Taherzadeh,Experimental and economical evaluation of a novel biogas digester,Energy Conversion and Management, Volume 74, Pages 183191, October 2013.

Vianney Tumwesige, David Fulford, and Grant C. Davidson, Biogas appliances in SubSahara Africa, Biomass and Bioenergy, In Press, Corrected Proof, Available online 15 March 2014.

Tushar Jah, and Sujay Basu, Development of a minibiogas digester for lighting in India, Energy, Volume 24, Issue 5, pp 409 411, May 1999.

Li Yingjian, Qiu Qi, He Xiangzhu, Li Jiezhi, andEnergy balance and efficiency analysis for power generation in internal combustion engine sets using biogas,Sustainable Energy Technologies and Assessments, Volume 6, Pages 2533, June 2014.

A. GonzÃ¡lezGonzÃ¡lez, M. CollaresPereira, F. Cuadros, and T. Fartaria, Energy selfsufficiency through hybridization of biogas and photovoltaic solar energy: an application for an Iberian pig slaughterhouse, Journal of Cleaner Production, Volume 65, pp 318323, 15 February 2014.

Abdullah Akbulut, Technoeconomic analysis of electricity and heat generation from farmscale biogas plant: Ã‡iÃ§ekda case study,Energy, Volume 44, Issue 1, pp 381390, August 2012.

C.X. CÃ¡ceres, R.E. CÃ¡ceres, D. Hein, M.G. Molina, and J.M. Pia, Biogas production from grape pomace: Thermodynamic model of the process and dynamic model of the power generation system, International Journal of Hydrogen Energy, Volume 37, Issue 13, pp 1011110117, July 2012.

G.N. Tiwari, A. Chandra, K.K. Singh, S. Sucheta, and Y.P. Yadav, Studies of KVIC biogas system coupled with flat plate collector, Energy Conversion and Management, Volume 29, Issue 4, pp 253 257, 1989.

G.N. Tiwari, S.K. Singh, and Kailash Thakur, Design criteria for an active biogas plant, Energy, Volume 17, Issue 10, pp 955958, October 1992.

R.A. Gupta, S.N. Rai, and G.N. Tiwari, An improved solar assisted biogas plant (fixed dome type): A transient analysisEnergy Conversion and Management, Volume 28, Issue 1, pp 5357, 1988.

Su Yuan, Tian Rui, and Yang Xiao Hong, Research and analysis of Solar Heating Biogas Fermentation Systems, Procedia Environmental Sciences, Volume 11, pp 13861391, 2011.

Feiqing Dong,and Jianbo Lu, Using solar energy to enhance biogas production from livestock residue A case study of the Tongren biogas engineering pig farm in South China, Energy, Vol.57, pp 759765,2013.
Nomenclature A Area (m2)
Horizontal area of the gas holder exposed to solar radiation (m2)
Vertical area of the gas holder which is exposed to solar radiation (m2)
Area of the top (m2)
Slurry vertical area (m2)
Vertical area of the gas holder which is submerged in the slurry (m2)
Area of flat plate collector (m2)
Cf Specific heat capacity of fluid (Antifreeze liquid) (J/kg C)
Cs Specific heat capacity of slurry (J/kg C) dxElemental section
Ex Exergy (W)
FPC Flat plate collector
Flow rate factor (dimensionless)
Heat transfer coefficient (W/m20C)
Radiative heat transfer coefficient (W/m20C)
1Heat transfer coefficient from gas holder plate to gas (W/m20C)
2 Convective heat transfer coefficient from gas holder plate to ambient (W/m2 0C)
3Heat transfer coefficient from gas to slurry (W/m20C) h4Heat transfer coefficient from slurry to ground (W/m20C)
Heat transfer coefficient from gas holder to slurry (W/m20C)
hs Heat transfer coefficient inside the tube from tube to slurry (W/m20C)
hwHeat transfer coefficient inside the tube from water to tube (W/m20C)
Heat transfer coefficient from slurry to air (W/m20C) I (t) Incident solar intensity (W/m2)
K Thermal conductivity (W/m K)
fMass flow rate of flowing fluid (kg/sec) MS Mass of slurry (kg)
N Number of photovoltaic thermal flat plate collector connected in series
Number of sunshine hours (hr)
, Rate of useful thermal energy transfer (kW) r1Inner radii of the tube (m)
r2 Outer radii of the tube (m) t Time (sec)
T Temperature (0C)
Ta Ambient temperature (0C)
TfoN Outlet temperature of fluid of the Nthphotovoltaic thermal flat plate collector (0C)
Tfi Inlet temperature of fluid in the photovoltaic thermal flat plate collector (0C)
Tg Gas holder temperature (0C) Tp Plate temperature (0C)
TsSlurry temperature (0C) TwFluid temperature (0C)
UOverall heat transfer coefficient for the system (W/m20C)
Absorptivity of dome Subscripts
a Ambient air eff Effective
ele Electrical g Glass
s Slurry
w Water Greek letters
Absorptivity of solar cell
Absorptivity of dome
()eff Product of effective absorptivity and Transmissivity Transmissivity
m Module efficiency