Coupled Analysis of Leachate Recirculation and Heat Transfer for a Landfill with Buried Waste Tyres

Huijie Zhang*, Zequan Li, Haixia Xu, Xinqian Guo, Huiling Wu, Junyi Sun
Faculty of Civil Engineering and Mechanics, Jiangsu University, No. 301 Xuefu Road Zhenjiang, Jiangsu, 212013, China

Article Metrics

CrossRef Citations:
Total Statistics:

Full-Text HTML Views: 176
Abstract HTML Views: 147
PDF Downloads: 92
ePub Downloads: 50
Total Views/Downloads: 465
Unique Statistics:

Full-Text HTML Views: 133
Abstract HTML Views: 75
PDF Downloads: 63
ePub Downloads: 40
Total Views/Downloads: 311

© 2018 Zhang et al.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: ( This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author at the Faculty of Civil Engineering and Mechanics, Jiangsu University, No. 301 Xuefu Road Zhenjiang, Jiangsu, 212013, China; Tel: 0086-0511-88797036; E-mail:



In bio-reactor landfills, both moisture and temperature have a significant impact on the bio-degradation processes. In order to speed up the leachate recirculation, a new drainage system is proposed to build a space network of bonded whole waste tyres in landfills.


In this study, a coupled dual-permeability flow and the heat transfer model was constructed to study the interaction between leachate migration and temperature evolution.

Results and Conclusion:

The established model was applied to a simplified bio-reactor landfill and the simulation results demonstrated that the presence of waste tyres could significantly speed up the heat transfer. And the intermittent injection and the varying injection pressure might affect the distribution of temperature as well. Additionally, cooling the leachate injected could be another measure in avoiding the occurrence of hot spots in landfills.

Keywords: Bioreactor landfill, Leachate recirculation, Bonded waste tyres, Temperature distribution, Cooling, Bio-degradation processes.


In recent years, bioreactor landfills have become an innovative and efficacious means for the safe disposal of Municipal Solid Waste (MSW) and biogas harvesting [1, 2]. The biogas production is influenced by water content, temperature, pH, and other factors, e.g., availability of oxygen, alkalinity, nutrients and inhibitors [3]. In bioreactor landfills, the collected leachate is recycled into municipal solid waste, eventually increasing the water content and accelerating the biodegradation of municipal solid waste [4-9]. In the presence of recirculated leachate, pH tends to remain stable probably because the injected leachate serves as a buffer solution; whereas, pH of the landfill without leachate recirculation declines too much below neutrality which inhibits biodegradation after a certain period [10]. The field data has proved the more rapid degradation and distinctly higher proportion of methane of the landfill with leachate recirculation compared with those without recirculation [11]. Furthermore, leachate recirculation can maximize the efficiency and waste volume reduction rate of landfill sites [12]. However, the leachate should be recirculated to all part of landfills by the employment of a recirculation system, because non-uniform distribution could result in short-circuiting which could reduce gas production rate [13].

The leachate recirculation system in bioreactor landfills includes vertical wells (VWs) [14-17], Horizontal Trenches (HTs) [18-21], and Drainage Blankets (DBs) [22, 23]. Waste tyre chips are commonly used to form the leachate drainage layer, but there is still controversy over whether spontaneous combustion will occur in the leachate trench [24].

There is a consensus that reducing the thickness of the tire-filled layer would help to avoid potential combustion of the tire. On the other hand, laboratory experiences demonstrated that increases in temperature from 20°C to 40°C raise considerably the CH4 generation rate, while further temperature elevation causes cessation of methanogenesis [25]. Therefore, the heat needs to be controlled for the sake of both efficiency and safety of bio-reactor landfills.

The recirculation of leachate is extensively operated under pressure in order to distribute the liquid within a limited time [26]. However, the high injection pressure may lead to a excessive pore water near the slope, which potentially results in a slope failure. Therefore, the manner of injection, namely, continuous or intermittent operation has been investigated [27] to mitigate the crisis due to pressured injection.

In the authors’ preceding work [28], a spacial network of bonded whole waste tyres is assumed to be buried in municipal waste during landfill. Interfaces between tyres and waste are anticipated to make preferential flow paths for the leachate, and the connection of tyres would contribute to the enforcement of waste when subjected to any potentials of slope slide. Meanwhile, the new pattern of dispersed tyres results in a thinner drainage layer compared with the conventional tire grooves, thereby reducing the risk of underground fire. A dual-permeability mathematical model combined with the momentum balance equation is employed to simulate the flow of leachate in the bio-reactor landfill and local safety factor for slope under the condition of a pressured leachate injection. The results have proved that the spatial tire network can obviously accelerate the uniform distribution of leachate and enhance the slope stability as well.

Based on the theory of heat transfer and fluid mechanics, this study continues to establish a coupled mathematical model to investigate the interaction between leachate recirculation and heat distribution in a landfill with waste tyres embedded. The numerical simulation of both fluid and the temperature field for a simplified landfill aims to provide fundamental knowledge for the future practice of this newly proposed drainage system.


2.1. Leachate Flow Model

The Richards equation is used to describe the flow of landfill leachate without tyres, whilst the dual- permeability model combines two Richard’s Equations to represent the flow within the domain of waste matrix and the priority flow at the tire-waste interface separately. The two domains influence and act on each other. When using w, denoting the volumetric ratio of the priority domain or matrix domain with respect to the total soil volume (-), the multiple bulk soil properties are constituted by domain-specific matrix and fracture pore system properties as shown in Eq (1a-1e).


The two flow domains are considered to be overlapped and interacted [29] (Eq. 2 and 3)


where the subscript f indicates the preferential flow domain and the subscript m indicates the matrix domain,C is the differential water capacity (dθ/dh) (L-1), Se is the effective saturation (–), h is the pressure head (L), t is the time (T), K is the hydraulic conductivity (LT-1), Ss is the specific storage (L-1), and Γw denotes the water exchange term (T-1) between the two domains.

The van-Genuchten function is used to describe the hydraulic properties of both the matrix and preferential flow domains [30]. The water exchange rate is assumed to be proportional to the pressure head difference between the fracture and matrix pore systems as follows [31] (Eq. 4).


in which αw is a first-order mass transfer coefficient (L-1T-1) for water defined as Eq. (5).


Where β is a dimensionless geometry coefficient for describing the shape of the matrix structures, a is the characteristic length of the matrix structure (L), γw is a dimensionless scaling coefficient for water absorption, and Ka is the relative hydraulic conductivity. A single average value of 0.4 for γw was found to be applicable for different hydraulic properties and initial conditions [32] and the coefficient β is equal to 3 [33]. The relative hydraulic conductivity Ka is calculated by averaging the hydraulic conductivities of two pore domains as Ka = (Km+Kf)/2 [34, 35].

2.2. Heat Transfer Model

According to Fourier’s law and the law of conservation of energy, we can get (Eq. 6).




In the above formulas (Eq. 7), the subscript s denotes the skeletal of waste, and l is the pore liquid, namely the leachate in this study.T is the temperature (K), t represents time (s), cs is the heat capacity of solid portion of waste (kJ/(m3K)), cl is the heat capacity of pore liquid, λ is the body of the heat conduction coefficient (J/(ms)·K), is the Darcy flow rate of pore liquid (m/s), Φ is the heat source, Se is the effective saturation (-).

In the landfills, the temperature rise is mainly due to waste heat caused by degradation. The instantaneous heat released by the degradation of organic waste at any time t is approximated as [36, 37] Eq. (8).


and the coefficient k is a function of temperature:

2.3. Coupled Dual-Permeability and Thermal Model

For the temperature-dependent dynamic viscosity coefficient, it can be calculated according to empirical formula [8, 11] (Eq. 9).


Where T denotes the temperature in a unit of degree (Celsius), and the resultant µ in cm2/s. Taking into account the relationship between conductivity and permeability coefficient, the conductivity at current temperature can be obtained provided that the conductivity is known as listed in Table 1.

Table 1. Summary of parameters used (at normal temperature).
Symbol Parameter Value Unit
γa Unit weight 5.5 kN/m3
α Inverse of air-entry pressure 1.4 1/kPa
θr Residual moisture content 0.14 -
θs Saturated moisture content 0.58 -
n Van Genuchten fitting parameter 1.6 -
a Van Genuchten fitting parameter 0.37 -
b Van Genuchten fitting parameter 0.5 -
c Van Genuchten fitting parameter 0.5 -
kv Hydraulic conductivity 1×10-7 m/s
Kh/kv Hydraulic conductivity anisotropy of MSW - -
λbs Thermal conductivity of waste solid 0.216 W • m-1 • K-1
cl Heat capacity of Leachate 4600.1 kJ/(m3 • K)
cs Heat capacity of waste solid 1200.5 kJ/(m3 • K)
a The van Genuchten model parameters of MSW are after [18] b Thermal analysis parameters refers to [22]

Two Richards equations are coupled by the water exchange term, and provide the pore liquid pressure and the associated liquid velocity for the calculation of the convective flux of heat transfer. The schematic diagram of coupling approach is illustrated in Fig. (1). The detailed computational method and the verification of dual-permeability modelling refer to the authors’ previous work [28].

Fig. (1). Coupling method of dual-permeability flow and heat transfer models.


3.1. Application to a Landfill Profile

The configuration of the landfill considered is analogous to the previous studies [26] as shown in Fig. (2). The leachate was injected through the horizontal trench both at a constant pressure and in a one-week on-and-off pattern, respectively. The Leachate Collection and Removal System (LCRS) is located at the bottom of the landfill and consists of sand and other substances, which has a thickness of 0.3 m. The pore water pressure at the bottom is fixed to zero because of the drainage layer. Since the bottom is in contact with the leachate for long-term, the temperature is assumed to be 50°C. For both leachate migration and heat transfer, the top and side slope are no-flux boundaries while the left-side is a symmetric boundary. At the top, there is no flux of leachate and the calculated annual average temperature is specified (22.5°C in the present study) to describe its contact with the atmosphere.

Fig. (2). Conceptual model for leachate recirculation using horizontal trench system.

Since the present study focuses on the effects of the presence of bonded tyres on the leachate flows and heat transfer, the influences of site-specific situations, such as the location of the landfill, abd the composition of waste, are represented by appropriately choosing the parameters used in the following simulations, which are listed in Table 1. In this study, the unsaturated hydraulic properties of the MSW are not varied with respect to the depth because: (1) very little published information is available on the evolution of unsaturated hydraulic properties of MSW as a function of overburden pressure; and (2) the unsaturated hydraulic properties have a relatively small impact on the key design parameters at steady flow state conditions [22].

3.2. Effect of Heat Transfer on the Rate of Leachate Recirculation

A base scenario without thermal effect is simulated using dual-permeability model and the flow rate, namely, the integration of leachate flux around injection trench is taken as an index to measure the effectiveness of injection. Fig. (3) demonstrates that the coupled models include the influence of varying temperature on dynamic viscosity of leachate and lead to a much higher unit flow rate. This implies the importance of considering the influence of temperature in predicting the leachate flow in landfills.

Fig. (3). Effect of temperature field on the flow rate: Red and blue colors are for injection pressure of 49 kPa and 98 kPa, respectively.

3.3. Effect of Tyres on the Distribution of Temperature

Four representative points are used to probe the temperature distribution in landfill. They are located at point 1 (15,10), point 2 (15,25), point 3 (15,40), point 4 (55,25) and the coordinates are in the unit of meter. Fig. (4) illustrates that the temperature is higher in deeper location and the leachate recirculation speeds up the release of internal heat. The indicates that the convective heat flux due to the leachate flow significantly decreases the subsurface temperature, which consequently avoids too hot temperature where consequently the microbial activity would otherwise be inhibited.

Fig. (4). Effect of leachate recirculation on the temperature field.

3.4. Effects of Injection Method on the Distribution of Temperature

3.4.1. Influence of Injection Pattern

Intermittent injection is commonly utilized to avoid excessive pore pressure and the slope failure. Here, Fig. (5) shows the temperature evolution with time for point 1, 2 and 4. The temperature at shallower location experiences larger fluctuation than deeper points. The reason is probably because of the location of the trench. Whatever the depth is, the average temperature for intermittent injection is slightly higher than that for a continuous pattern. Further, the same trend is also confirmed in the contour illustrations as shown in Fig. (6).

Fig. (5). Effect of intermittent leachate recirculation on the temperature field (injection pressure = 98 kPa): solid curves are for intermittent injection, and dotted curves are for continuous injection.

Fig. (6). Evolution of the temperature field for both continuous (left) and intermittent (right) injection: The temperature is in unit of °C.

3.4.2. Influence of Injection Pressure

Fig. (7) demonstrates that raising the injection pressure results in a lower temperature due to the quicker leachate flow and the related larger convective heat flux. However, the influence of higher injection pressure on the heat transfer is not as significant as the contribution of tyres. Therefore, the importance of waste tyres as a new leachate drainage system has been approved.

Fig. (7). Effect of injection pressure on the temperature field: solid curves are for pressure of 98 kPa, and dotted curves are for injection pressure of 49 kPa; and red, blue and black curves are for point1, point2 and point3, respectively.

3.5. Effect of Cooling on the Distribution of Temperature

To account for the possible extra heat generation because of the presence of tyres, a preliminary discussion is made here in the absence of an accurate description due to the scarcity of field data. The initial temperature for the cases with tyres is designated to be 10°C above the normal temperature. To meet the challenge of rising temperature (which would exceed the favorite temperature range for biodegradation), cooling by adding ice to the recirculated leachate is attempted by specifying the temperature of leachate injected at 10°C below the normal temperature (Fig. 8). It has been found that cooling makes a considerable difference in the upper part of the landfill, but much smaller at a deeper section. Therefore, the cooling well should be placed near the hot spots. In addition, further investigations of in-situ temperature monitoring need to be performed to improve the accuracy of heat generation and the bottom temperature boundary condition for numerical simulation.

Fig. (8). Effect of cooling on the temperature field when enhanced heat generation due to tyres is considered: solid curves are for cases with cooling, and dotted curves without cooling; and red, blue and black curves are for point1, point2 and point4, respectively.

3.6. Feasibility Discussion

When compared with the conventional practices for a landfill which does not accept waste tyres, the additional costs of the proposed drainage system in this study arise from waste tyres baling and cooling of leachate. These two spendings are both routine for a waste tyres monofill or an inert waste landfill which assign certain cells for tyres exclusively. Regularly, a landfill operator requires that the waste tyres should be baled prior to delivery because it will be easier to handle, help save the valuable landfill space, and will have a greater potential for recovery in the future [38]. Especially for the passenger tyres, landfilling will be a better choice than recycling. Since the large capital investment required for a relatively low return, due to the high percentage of fibre and lower volumes of rubber compared with other tyres, makes it unlikely to be recycled in the short to medium term. Therefore, receiving the waste passenger tyres to construct drainage system would be an excellent choice and the fee paid by tyres collectors would help cover the expense related to tyres baling and possible cooling of leachate.

Further, the high temperature in bio-reactor landfill appears to cause methanogenic mesophile and thermophile bacteria become dormant or perish. Meanwhile, hyperthermophilic bacteria will thrive and the hot zone expands until the entire landfill is converted to a high-temperature environment, which presents many operational challenges [39]. It would be more cost-efficient if the cooling facility is integrated with a geothermal retrieval system. A case of small-scale landfill in New Hampshire, USA, has demonstrated that the payback period of the cooling system was around 7 years [39]. However, the method proposed here in this study could help alleviate the emergency for a landfill without looping in a simple manner.

Finally, floating of tyres [40] to break through the daily cover soil will help mitigate the slope stability problem, which is caused by the low conductivity of the cover soil. The cover soil was demonstrated to bring risk in the slope stability because it intercepts the leachate and causes a subsequent lateral transmission of leachate to the side slope [4]. Therefore, the improvement in safety by virtue of this new drainage system is an extra benefit besides the economic consideration.


This study coupled dual-permeability with heat transfer model to describe the non-isothermal leachate migration in a bio-reactor landfill with waste tyres embedded. The effectiveness of leachate recirculation was found to be improved by the decreasing dynamic viscosity of leachate due to the increasing temperature, which implies the importance of including the interaction between leachate seepage and heat transfer. Further, the injection method, namely, the pattern of intermittent or continuous injection and the injection pressure would influence the distribution of the internal heat in landfill. However, the presence of waste tyres net leads to a more significant lowering of temperature. Therefore, it would help release the concern of occurring of subsurface hot spots when applying this new drainage system to landfills. Additionally, cooling the leachate before injecting back at a trench near a potential hot spot could be another effective measure in avoiding the excess internal temperature. Finally, the limitations remain in this study because of the assumption of the heat generation in the presence of waste tyres and the temperature boundary condition at the bottom, which need to be resolved in further field monitoring. However, the preliminary work aims to shed light on the interaction of leachate seepage and heat transfer, and draw more attention to this alternatively new drainage system in bio-reactor landfills.


Not applicable.


The authors declare no conflict of interest, financial or otherwise.


This work was supported by the National Natural Science Foundation of China (Grant No. 51308259); the Foundation for Advanced Talent of Jiangsu University, PRC (Grant No. 13JDG060 and No. 14JDG023); and the sponsorship of Jiangsu Overseas Visiting Scholar Program for University Prominent Young & Middle-aged Teachers and Presidents.


[1] Reinhart DR, Townsend TG. Landfill bioreactor design and operation 1997.
[2] Sharma HD, Reddy KR. Geoenvironmental engineering: site remediation, waste containment, and emerging waste management technologies 2004.
[3] Manzur SR. Effect of leachate recirculation on methane generation of a bioreactor landfill M.S. Thesis, University of Texas at Arlington, Arlington, TX, USA, 2010.
[4] Barlaz MA, Ham RK, Schaefer DM. Microbial, chemical and methane production characteristics of anaerobically decomposed refuse with and without leachate recycling. Waste Manag Res 1992; 10(3): 257-67.
[5] Chugh S, Clarke W, Pullammanappallil P, Rudolph V. Effect of recirculated leachate volume on msw degradation. Waste Manag Res 1998; 16: 564-73.
[6] Reinhart DR, McCreanor PT, Townsend T. The bioreactor landfill: Its status and future. Waste Manag Res 2002; 20(2): 172-86.
[7] Reddy KR. “Geotechnical aspects of bioreactor landfills”, in Indian Geotechnical Conference 2006-Geotechnical Engineering-Indian Experience. 2006, pp.79-94.
[8] Townsend TG, Powell J, Jain P, Xu Q, Tolaymat T, Reinhart D. Sustainable practices for landfill design and operation 2016.
[9] Ko JH, Yang F, Xu Q. The impact of compaction and leachate recirculation on waste degradation in simulated landfills. Bioresour Technol 2016; 211: 72-9.
[10] Morris JWF, Vasuki NC, Baker JA, Pendleton CH. Findings from long-term monitoring studies at MSW landfill facilities with leachate recirculation. Waste Manag 2003; 23(7): 653-66.
[11] Sponza DT, Ağdağ ON. Impact of leachate recirculation and recirculation volume on stabilization of municipal solid wastes in simulated anaerobic bioreactors. Process Biochem 2004; 39: 2157-65.
[12] Giri RK, Reddy KR. Slope stability of bioreactor landfills during leachate injection: Effects of unsaturated hydraulic properties of municipal solid waste. Int J Geomech 2013; 8: 144-56.
[13] Sawatdeenarunat C. Effects of leachate recirculation on anaerobic treatment of municipal solid waste 3rd International conference on geoinformation technology for natural disaster management and rehabilation vol.4: 448-60.2010;
[14] Feng SJ, Cao BY, Zhang X, Xie HJ. Modeling of leachate recirculation using vertical wells in bioreactor landfills. Environ Sci Pollut Res Int 2015; 22(12): 9067-79.
[15] Feng SJ, Cao BY, Xie HJ. Modeling of leachate recirculation using spraying-vertical well systems in bioreactor landfills. Int J Geomech 2017; 17: 04017012.
[16] Khire MV, Haydar MM, Haydar MM. Leachate recirculation using geocomposite drainage layer in engineered msw landfills Geo-Frontiers Congress 2005; 1-11.
[17] Jain P, Townsend TG, Tolaymat TM. Steady-state design of horizontal systems for liquids addition at bioreactor landfills. Waste Manag 2010; 30(12): 2560-9.
[18] Reddy KR, Giri RK, Kulkarni HS. Design of horizontal trenches for leachate recirculation in bioreactor landfills using two-phase modeling. Int J Environ Waste Manag 2015; 15: 347-76.
[19] Chan GYS, Chu LM, Wong MH. Effects of leachate recirculation on biogas production from landfill co-disposal of municipal solid waste, sewage sludge and marine sediment. Environ Pollut 2002; 118(3): 393-9.
[20] Reddy KR, Giri RK, Kulkarni HS. Validation of two-phase flow model for leachate recirculation in bioreactor landfills. Waste Manag 2014; 24: 1-24.
[21] Haydar MM, Khire MV. Leachate recirculation using horizontal trenches in bioreactor landfills. J Geotech Geoenviron 2005; 131: 837-47.
[22] Reddy KR, Giri RK, Kulkarni HS. Two-phase modeling of leachate recirculation using drainage blankets in bioreactor landfills. Environ Model Assess 2015; 20: 475-90.
[23] Haydar MM, Khire MV. Leachate recirculation using permeable blankets in engineered landfills. J Geotech Geoenviron 2007; 133: 360-71.
[24] Hall D, Drury D, Keeble R, Morgans A, Wyles R. Review and investigation of deep-seated fires within landfill sites Environment Agency, Bristol, Science Report no. SC010066, 2007.
[25] Kasali GB, Senior E. Effects of temperature and moisture on the anaerobic digestion of refuse. J Chem Technol Biotechnol 1989; 44: 31-41.
[26] Xu Q, Tolaymat T, Townsend TG. Impact of pressurized liquids addition on landfill slope stability. J Geotech Geoenviron 2012; 138: 472-80.
[27] Sang NN, Soda S, Inoue D, Sei K, Ike M. Effects of intermittent and continuous aeration on accelerative stabilization and microbial population dynamics in landfill bioreactors. J Biosci Bioeng 2009; 108(4): 336-43.
[28] Zhang H, Jie H, Jeng D, Xu H. Effects of bonded tyres on leachate recirculation of bioreactor landfills. Open Civ Eng J 2017; 11: 552-62.
[29] Gerke HH, Van Genuchten MT. A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media. Water Resour Res 1993; 29: 305-19.
[30] Van GMT, Wierenga PJ, O’Connor GA. Mass transfer studies in sorbing porous media, 3: Experimental evaluation with 2,4,5-t trichlorophenoxyacetic acid, movement of herbicide through soil. Soil Sci Soc Am J 1977; 41: 1-10.
[31] Ray C, Ellsworth TR, Valocchi AJ, Boast CW. An improved dual porosity model for chemical transport in macroporous soils. J Hydrol (Amst) 1997; 193: 270-92.
[32] Gerke HH, Van Genuchten MT. Evaluation of a first-order water transfer term for variably saturated dual-porosity flow models. Water Resour Res 1993; 29: 1225-38.
[33] Gerke HH, Van Genuchten MT. Macroscopic representation of structural geometry for simulating water and solute movement in dual-porosity media. Adv Water Resour 1996; 19: 343-57.
[34] Arora B, Mohanty BP, McGuire JT. Inverse estimation of parameters for multidomain flow models in soil columns with different macropore densities. Water Resour Res 2011; 47(4): WR009451.
[35] Laine-Kaulio H, Backnäs S, Karvonen T, Koivusalo H, McDonnell JJ. Lateral subsurface stormflow and solute transport in a forested hillslope: A combined measurement and modeling approach. Water Resour Res 2014; 50: 8159-78.
[36] Yang J. Study of temperature-chemical coupling effect in landfill disposal of household rubbish in urban area PhD thesis, Southwest Jiaotong University, Chengdu, China, 2007.
[37] Han Z, Liu D, Li Q. Influence of leachate recirculation on producing methane in anaerobic bioreactor. Environ Sci Technol 2008; 31: 160-3.
[38] Matthews M. Investigation into the environmental, social and economic impacts of a potential banning of used tyres to landfill and of mandatory recycling of used tyres Available from: http:// media/files/ documents/ impacts_tyre_policies.pdf
[39] Grillo RJ. Energy recycling–landfill waste heat generation and recovery. Current Sustainable/Renewable Energy Reports 2014; 1: 150-6.
[40] Solid Waste Agency “The Problem with Waste tyres”, Available from: Accessed: 2018- 07-17.