Design and Performance Evaluation of a Hydronic Type Compost Heat Exchanger

Abstract While much research has been published on the Compost Heat Recovery Systems (CHRs), little has been documented on the design and performance evaluation of the Hydronic compost heat exchangers using numerical and computational methods, occasionally resulting in compost process inhibition. A CHRs (0.036 m3/7.2 m2) Hydronic-type heat exchanger and 12.43 m2/2.83 m3, compost reactor (CR), was designed and developed with the main objective of evaluating the design and its performance. The numerical design and performance evaluation was achieved by using Kern’s and the effectiveness and Number of Transfer Units methods (ε-NTU), respectively. Empirically, data were captured by using the Polytetrafluoroethylene (PTFE) thermocouples connected to the TC-8 Picolog Data loggers. Data validation (empirical and mathematical), was achieved by modifying a free computer-based software developed by the Chemical Engineering Calculations (CHECAL), into a Hydronic Compost Heat Exchanger design and performance evaluation software (HYDROCOHE). Between the HYDROCOHE and numerical, and between empirical and HYDROCOHE, R2 values of 0.99938–0.9995, and R2 of 0.99269–0.9432 with the effectiveness of 0.4853–0.4848 were achieved with 0.99 kW-empirical and 2.10 kW-HYDROCOHE, respectively. The power disparity may be ascribed to the compost reactor’s insufficient thermal insulation. Counterflow arrangement was more effective (0.4766) than crossflow (0.4622) and parallelflow (0.4430) setups. Parallelflow heat exchanger system, therefore, has the potential to extract heat steadily, minimizing the composting cycle inhibition. Further work on the impact of various flowrates on the direction of flow and heat extraction is recommended.


PUBLIC INTEREST STATEMENT
Although heat utilization from the Compost Heat Recovery Systems (CHRS) has been reported to have started 2000 years ago in China by digging trenches and filling them with biomass materials and covering with a layer of soil for winter crop production, little information has been reported on the design of hydronic heat exchangers to efficiently extract heat and maintain the composting process. With temperatures above 65 °C reported, the composting process if properly utilized and designed could be the source of energy for heating water, pre-heating milk during the sterilizing process, generation of electricity, and other uses to reduce Green House Gases (GhG). The waste products from composting, is a rich soil amendment agent. This paper reports a research that was conducted to improve the heat extraction process during composting to avoid the inhibition of the process. It is believed that coming up with standard composting processing parameters will improve the utilization of energy and farm bio wastes and reduce GhG.

Introduction
Composting is a microorganism, biologically influenced and an exothermic process from which it is possible to capture valuable thermal energy that can be used to heat air or water. It has long been regarded as a high thermal energy generator, although most of it is lost to the surrounding environment and research on utilization and composting heat recovery systems (CHRSs) have been inconsistent (Smith et al., 2017;Walling et al., 2020).
Composting is divided into aerobic, where ample oxygen is present and anaerobic, where decomposition happens in the absence of oxygen (Anaerobic digestion) (Bartocci et al., 2020;Misra et al., 2003). In this study, much attention will be on aerobic composting which comprises of four major phases microbiologically in relations to temperature, namely, mesophilic (25 °C-45 ° C), Thermophilic (45 °C-71 °C), and cooling and maturation. The heat produced is grouped into lowtemperature heat sources and would otherwise be wasted if not used. It can be utilized in heating green and hoop houses, buildings, and hot water (Smith et al., 2017).
Historically, heat utilization from CHRS in China 2000 years ago by digging trenches and filling them with biomass materials, and covering with a layer of soil for winter crop production, has been reported (M. M. M. Smith & Aber, 2018). In France, heat extraction has been reported for winter farming and season extension (Zantedeschi, 2018). In recent years, the Jean Pain type of CHRS has been modified into Biomeiler, thermocompost, bioreactors, or compost reactors (Pelleton, 2014). To bridge up the scarcity of the evidence and data of the performance and heat extraction from the CHRS, experiments have been performed with medium to large-sized compost reactors with full data presentation on the compost feedstock composition, and temperature profiles documented (Zantedeschi, 2018).
Heat recovery methods from the composting process were reviewed and grouped into four recovery systems and challenges highlighted: The CHRS have been used to determine the energy available in the compost reactors for commercial and academic purposes, ranging from a small area (pilot) to large in-vessel composting facilities. The average recovery rates reported are uncertain, and vary widely. The laboratory-scale systems reported average rates of 1.9 MJ hr −1 (1.2 MJkg −1 dm), 20 MJ hr −1 for pilot scale testing (Bajko et al., 2019), and 205 MJ hr −1 commercially (7.1 MJ kg −1 Dry Matter (dm)).
Different researchers have indicated that composting below 60 °C results in quality and highest composting rates and achievable via high airflow rates. Temperatures above 60 °C are not conducive to the thermophilic microorganisms (Gao et al., 2010;Hoitink, 2014;Kulcu & Yaldiz, 2004;Kuter et al., 1985;Liang et al., 2003). Extra heat requires extraction to avoid composting process inhibition and passive or positive aeration can be applicable to achieve that.
In order to achieve an efficient composting process, clear knowledge and understanding of the process in terms of the physical (moisture content, bulk density, etc.), mechanical (porosity, permeability, etc.), and chemical (C/N ratio, nutrients content, etc.) properties of the materials involved, are key. The properties influence the process of aeration and compostability effectiveness. The importance of the physical and biological parameters in the designing process of the composting systems is provided and highlighted by (Keener et al., 1993). The highlighted properties that were key to these studies are; Bulk density (BD wet&dry ), Particle size (diameter), Total Porosity (TP or ), Volumetric water content (V wc ), Free Air Space (FAS), Permeability (K), Superficial velocity (v), Permeability based Reynolds number (Re) (dimensionless), Passibility (η), Pressure difference (ΔP), and Densities of ambient air and air Density in the CR.
However, the efficiency of heat extraction methods depends on the flow rate and temperature of the fluids been used. The higher the flow rate and lower temperatures of the extracting fluid, the higher the rate of extraction hence the higher the chances to inhibit the composting process (Smith et al., 2017).
Major challenges associated with heat recovery without negative impacts on the quality of the compost with the economical value of the energy recovered have been highlighted (R. F. Zhao et al., 2017). This challenges calls for caution, at the design stage because the more heat extracted, the greater the possibilities of deactivating the microbiological process. The effect of heat extraction on the compost process using the pilot scale compost reactors was evacuated and data mined from the tests used in updating the COMSOL Multiphysics TM computational model. Despite the simulated and empirical data been in synchro from the beginning of composting, the data diverged during the thermophilic temperatures (peak) and was attributed to the inappropriate heat transfer boundary conditions in the model (Nwanze & Clark, 2019). The limitations imposed on the design of heat exchangers mostly for compost heat extraction include: unavailability of standardized design procedures (Mason & Milke, 2005), circulating power requirements, requirements for spatial dimensions, unavailability of materials and standards, expertise and technology availability (Rahim & Khaled, 2017). The heat generation of the composting system is site and feedstocks quality specific, the factor that calls for the site and specific composting process designs (J. M. Haug, 1993).
The above challenges associated to the CHRs heat generation and extraction need to be investigated and hydronic heat extractors' design and performance factors, such as the area to volume ratios, flowrates, and feedstock's physical parameters, standardized. These challenges can partly be mitigated by avoiding the try and error heat recovery practices commonly practiced in heat recovery from CHRs by employing the conversional heat exchangers designs. Therefore, this study was based on an opportunity to extract the extra heat generated by the microbial action in the CR. The CR was utilized as a shell side of the heat exchanger, which is a device used for transferring thermal energy between two or more thermal fluids with different temperatures (Sundén & Fu, 2017). Energy is transferred from one thermal fluid to the other across a solid surface.
The most important type of heat exchanger is the recuperator, in which the heat exchanging flowing fluids are on either side of a dividing wall. The second type is a regenerator, in which hot and cold fluids, respectively, move through a space containing a material matrix that provides alternative heat flow means. The third type is the process of evaporation in which a liquid is evaporatively and continually cooled in the same space as the refrigerant. The heat exchange takes place in a direct compact or open heat exchanger by direct mixing of hot and cold fluids and the simultaneous transfer of heat and mass (Adumene et al., 2016).
Nevertheless, unlike the conventional heat exchangers where the fluid in the shell or tube side recirculates thereby gaining or releasing heat, the heat source in this design was the feedstocks in the CR through the aerobic composting process resulting into a 1 shell-pass and multiple tube passes heat exchanger (cooler) arrangement (Alhusseny, 2010;Fateen, 2018). Balancing the shell side (CR) and the Tube side here referred to as Compost Heat Exchanger (COHE), the design was key.
The heat source (CR) and its feedstocks were initially investigated in order to achieve the optimal design of the heat extraction system. Various researchers have associated the attainment of 60 °C and above, by the composting process, to the achievement of optimum aerobic conditions, that resonate with good thermal properties. In their study, H. Ahn et al. (2009) investigated 12 compost-bulking materials for thermal conductivity, thermal diffusivity, and volumetric heat capacity at varying bulk densities, particle sizes, and moisture contents. Linear relationships were established between thermal conductivity and volumetric heat capacity with moisture content and density, while thermal diffusivity was nonlinearly related. Thermal conductivity and volumetric heat capacity levels between 0.12-0.81 W m −1 °C and 1.36-4.08 MJ m −3 °C, respectively, were attained (H. Ahn et al., 2009).
Heat conductivity coefficients of compost from municipal waste are said to depend on the temperature and density. The higher the compost density, the higher the compost thermal conductivity coefficient. A compost temperature of 60 °C, resulted in the thermal conductivity coefficient of 0.309 W m −1 o C at a density of 600 kgm −3 .as compared to the 0.250 W m −1 o C at 442 kgm −3 bulk density at the same temperature. Since the thermal conductivity coefficient decreases with compost age, care should be taken when designing heat extraction from the CR. Aeration is a supply line for both for oxygenation of the microorganisms and extra heat removal (Klejment & Rosiński, 2008).
In this study, techniques and empirical formulae of working out the physical, chemical and mechanical parameters were used as cited and used from the previous research works, and as presented in the methodology and results sections for comparison purposes Keener, 2008;Keener et al., 1991).
Designing of shell and tube heat exchangers is achieved by trial and error calculations. Kern and Bell-Delaware methods are commonly used for heat exchanger designing via performing Thermal Analysis and Hydraulic Analysis. The Kern method is the simplest route but not as accurate as the Bell-Delaware method in terms of factoring in the leakages and other fluid and pressure losses (Dhamodharan, 2018). The key components of the heat exchanger designing procedures are rating (Thermo-hydraulic evaluation) and sizing (Heat transfer rates, fluid flow rates, pressure drop, surface area for heat transfer, and inlet and outlet temperature determination) (NPTEL, 2018; R. K. Shah & Sekulic, 2003).
In this study the Kern method was applied and is systematically described in the methodology section, to mitigate on the design difficulties. Several design methods such as the Log Mean Temperature Difference (LMTD), where the total heat transfer rate is related to the inlet and outlet fluid temperatures and the overall heat transfer coefficient/total heat transfer surface area was applied in this study. The evaluation was achieved by the effectiveness and Number of Transfer Units (ε-NTU), which have been discussed (Theodore, 2011).
In order to use the LMTD method for determination of the heat exchanger size and design parameters, the outlet and inlet temperatures, mass flowrates (Hot and Cold Fluids), should be known. The process to follow is: • Selection of the heat exchanger and its suitability for the application.
• Using the energy balance to work out the unknown inlet and outlet temperatures and the heat transfer rate.
• Calculation of the Mean Temperature Difference (MTD) (T im ) and the correction factor F.
• Overall Heat Transfer Coefficient (U) selection or calculation.
• Finally calculating the area (A) of the heat exchanger transfer surface.
It is recommended that the heat exchanger heat transfer surface area should be equal or larger than (A) (Cengel & Ghajar, 2015).
However, when the outlet Hot and Cold fluid temperatures are not known, the LMTD method requires involving iterations to calculate the parameters. Therefore, effectiveness and NTU, the method formulated by Kays andLondon in 1955, is recommended (Kays &London, 1955). It is a method based on a dimensionless parameter known as the heat transfer effectiveness.
The ε-NTU method enables the determination of the heat transfer rate with no prior knowledge of the outlet temperatures of the fluids involved. The flow arrangement and the geometry of the heat exchanger determines the effectiveness of a heat exchanger. Effectiveness is heat exchanger type specific (Theodore, 2011).
The design and performance assessment of the CHRS heat exchangers is entirely justifiable from the above. Most of the previous research works heighted by (Allen & Chambers, 2009), did little to highlight this aspect leaving a knowledge gap. In his research thesis, Lekic (2005), highlighted the limitation of pipe works in the process of heat recovery from the composting process (R. Zhao et al., 2015).This aspect is key to ensuring that the aerobic-composting route of biomass conversion to energy is competitive, profitable, and act as a source of biological waste disposal method reducing the Greenhouse Gases emissions (GhG). The by-product (Compost Manure) is a good soil amendment product, particularly valuable for reducing cultivation costs and supporting organic farming in rural or small farms and anchoring on number 2, 3, 7, 11, and 13, Sustainable Development Goals (SDGs) (UNDP, 2020).
In this paper, a pilot scale hydronic type, Compost Heat Exchanger (COHE), counterflow (1-shell pass and multiple tube pass) was designed, developed, and evaluated using a pilot-scale static aerated Compost Reactor (CR) whose highest and minimum optimum temperatures were monitored and evaluated prior to heat extraction. The research was set up in Kitale, Trans-Nzioa County in Kenya. The COHE heated the 34 litres water tank.
The primary objectives were; • Designing and developing of the Hydronic Compost Heat Exchanger (COHE) using the Kern's and ε-NTU mathematical models, respectively, as guided by (Shah & Sekulic, 2003) and applied by Irvine.
• Experimental performance evaluation of the developed CHRs Profile.
• Validation of the mathematical and empirical models using the modified free computer-based software developed by the Chemical Engineering Calculations (CHECAL), into a Hydronic Compost Heat Exchanger design and performance evaluation software (HYDROCOHE software) (Chelcal, 2018).
• Evaluating the most effective heat exchanger flow direction suitable for the CHRs (counterflow, parallelflow, and Crossflow (single flow) both fluids unmixed heat exchangers).

Materials and methods
The Compost Reactor built for this study was first analyzed and average temperature profiles empirically captured and analyzed to determine the inlet average temperatures. This information was key and constituted initial parameters such as, the inlet (CR) maximum average temperature (T CR1 ), selection of water as the thermal fluid of the exchanger, and the inlet COHE average temperature (T COHE1 ). Calculation of the average specific heat capacities (C P ) of the CR from the biomass feedstocks composition (C PCR ) and the specific heat capacity of the thermal fluid (water) (Cp COHE ) also adopted to be 4.175 kJ kg −1 °C −1 . These were in line with .
In Figure 1, the methods used in this study are shown, ranging from the Heat source (CR) evaluation to the testing and evaluation of the CR/COHE system. The site weather and geographical conditions were taken into consideration. The methods in the design of the proposed COHE were based on the engineering principles and methodologies outlined by (Abd & Naji, 2017;Bergman et al., 2011;Hayati, 2014;NPTEL, 2019;Shah & Sekulic, 2003).

Site geographical and climatic parameters
This study was carried out in Kitale, Trans-Nzia County in Kenya with a geographical location of Latitude 01° 08ʹ 39 N and Longitude of 034° 59ʹ48E. The study area is located at an average of 1907 m elevation, with mean annual temperature: maximum 26.4 °C and 12.3 °C minimum; average annual mean rainfall of 1259.1 mm with 85.6 sunshine hours per annum.

Experimental measures
The empirical data mining involved the measuring of temperatures, humidity, fluid flowrates, and the relative humidity at various points initially in CR and finally during the experimental trials (heat extraction) with COHE (Tube side) installed inside the CR (shell side).

CR physical properties evaluation
The compost feedstocks mixture and composition were as described . Physical properties of the CR are important in ensuring an efficient aerobic composting system. The following were considered; Temperatures achieved Bulk density (BD wet&dry ), Particle size (diameter), Total Porosity (TP or ), Volumetric water content (V wc ), Free Air Space (FAS), Permeability (K), Superficial velocity (v), Permeability-based Reynolds number (Re) (dimensionless), Passibility (η), Pressure difference (ΔP), and Densities of ambient air and air Density in the CR. (H. M. Keener et al., 1993).

CR temperature profiles
In this research setup, the source of energy was the heat from the CR. The aeration of the CR was statically achieved through the chimney effect created by the spaces left in between the wooden panels used at the ground stage of the CR as shown in Figure 2 which forced the air to move upward, resulting in a counter flow heat exchanger system. The air flowed upward as shown by the label V arrow in Figure 2. The CR was monitored and its temperature profiles analyzed. According to the tendance of natural convection, hot and moist air is pushed into the upper area of the CR; therefore, much heat is concentrated on the upper part (Seng et al., 2018). A condenser-type heat exchanger for CHRs designed and tested by (Bajko et al., 2019), was assembled on top of the compost reactor with consideration of the natural convection effect.
In Figure 2 Tp1 to Tp6 shows the Polytetrafluoroethylene (PTFE) thermocouples connected to the TC-8 TM Picolog Data loggers with a temperature accuracy of ±0.2% of reading and ±0.5 °C were used inside the CR to transmit the data to the PC via a 15 m USB cable. The PTFE were selected due to their robustness in withstanding the harsh conditions within CR due to humidity, pH changes, and microorganisms' actions resulting in a chemical reaction. Y denotes the corn stalk and the mixed compost feedstock in the compost reactor. Bajko et al. (2018), designed a wooden rod that was inserted into the Polypropylene to avoid corrosive environmental effects on the sensors. The aeration in passively aerated compost systems is governed by passive convection, resulting from subjecting a fluid of constant viscosity and not compressible, to a temperature gradient and is described by Darcy's law. The force that drives the airflow is buoyant in nature and calculable by applying Archimedes principles. The buoyant force is said to be equal to the weight of the fluid been displaced (Das & Keener, 1997;Lynch & Cherry, 1996;Yu et al., 2006Yu et al., , 2009).
Average Air velocity; the following equation was used to calculate the average air velocity, Where v = average air velocity (ms −1 ), K = permeability of the compost (m 2 ), ρ Amb = density of ambient air (kg m −3 ), μ = fluid dynamic viscosity of the air (m 2 s −1 ), T Amb and T CR1 = ambient and compost reactor temperatures, respectively, g = acceleration (9.8 m s −2 ), L = the length (m) of the occupancy of the porous medium. The fraction ΔP Δγ = pressure gradient through the compost (Pa m −1 ) (Yu et al., 2006). Equation 1 gives a linear relationship that exists between the steady velocity and the pressure gradient in laminar flow conditions only and ceases under the turbulent and transitional flow. High and above 1 Re, Darcy's equation does not apply and the condition is no longer steady but turbulent (Notton, 2005;Siyu et al., 2005).
The density (ρ Amb ) of the air was worked out using the Psychometric relationships since the temperatures of the ambient air, altitude (1900 m above sea level), and relatively humidity was conveniently measured as explained in the temperature capture section (ASHRAE, 2020). The thermal conductivity coefficient for the composting process is a function of its temperature and density of the feedstocks mixture. The higher the density, the higher the thermal conductivity of the compost reactor (Klejment & Rosiński, 2008). Consequently, the density of the CR is a function of its mass and the space it occupies (volume) (Jain et al., 2019).
The value of the airflow velocity is key in efficient heat generation by supplying oxygen for the reaction of oxidation and by removing heat from the reactor. It is this requirement to remove extra heat that gives an opportunity of thermal energy capture in this study. There is a ceiling value on the aeration flowrate that succumbs to the non-continuous linear relation between the respiration activities and the airflow supply (Mejias et al., 2017).The higher the velocity, the higher the heat loss that can result in inhibition and the lower the velocity the lower the heat generation (Luangwilai et al., 2010).
In order to calculate the average air velocity using Equation 1, figures of the compost substrates such as, permeability (K), air viscosity (μ), Temperatures (T Amb andT CR1 ) in the compost, bulk density, density of ambient temperature, and Free Air Space (FAS), must be worked out. The importance of all these parameters is key in enabling the estimation of the effect of one parameter on another (physical properties) are highlighted by (J. .

Determination of Permeability
Published results were applied to determine the permeability values under the compaction of the compost in the CR, in this study since it was beyond the scope. The permeability values were estimated based on the published values with related values of bulk density, moisture content, and the initial Free Air Space (FAS o ) in respect to the compaction effect (Das & Keener, 1997;Notton, 2005;Yu et al., 2006). This method was applied by (Yu et al., 2006) to develop a practical analytical model of airflow in a passively aerated CR.
The ideal Gas Law and Archimedes principles were applied to compare the buoyant force acting on the air to the temperature difference between the air in the CR bed and ambient air. This was achieved under the assumptions of constant viscosity and negligible compressibility (H. K. Ahn et al., 2008;Siyu et al., 2005;Tiquia, 2005;Yu et al., 2006Yu et al., , 2009).
Darcy's law was used to describe the creeping flow of the fluid that is not compressed and passing through a porous CR and is used under less than 1 Reynolds number (Re) (Yu et al., 2006).

Reynolds number and Pressure drop
The equations applied were: Where Re = Reynolds for flow in a porous medium,ΔP= the pressure (Pa) deference driving force from the point of air entry into the CR until exist,ρ amb and ρ CR = density (kg m −3 ) of ambient and inside compost reactor air, g and H are gravitational acceleration (9.8 m s −2 ). H = the height of the compost reactor after feedstocks loading (m), v = average velocity of air, K = the permeability values, and μ = viscosity of the air.

Density
The density of the air at each point was calculated using: Where ρ amb and ρ CR ambient and CR air density, respectively, P da = is the pressure of dry air (Pa), P pv = water vapour pressure (Pa), T (T amb = 18.2 °C, T CR = 66.8 °C) = air temperature in o C, R d = is the specific gas constant for dry air (287.058 J/(kg o C) and R v = specific gas constant for water vapour equal to 461.495 J/(kg o C) (OC, (Omni Calculator sp. zo.o.), 2020).

Free Air Space (FAS) determination
The equation linking the physical parameters is as stated below as recommended and used by Das & Keener, 1997;McCartney & Chen, 2001;Mejias et al., 2017;Siyu et al., 2005;Yu et al., 2006), was applied in this study.
Where Free Air Space (%) = FAS, Bulk Density (wet basis (wb)) in (kg/m 3 ) = BD, Moisture Content (wb) (%) = MC, the density of water (999 kg/m 3 ) = ρ w , and Particle density (kg/m 3 ) = PD. The equation is applicable when all variables are known. All other particles are easily worked out apart from particle density .

Bulk density (%) (BD)
The bulk density in this study was measured at different pile levels (depth) using the mass per unit volume technique.

Dry and Wet Bulk Density (BD dry ) (kg/m 3 ):
The dry bulk density was calculated using the following equation:

Total Porosity (TP or )
The total porosity was workout using the following equation:

Volumetric water content of the feedstocks
The volumetric water content of the feedstock matrix was calculated using : Given; BD dry and BD wet = dry and wet bulk density, respectively, of the feedstocks used in this experiments (kg m 3 ), V WC = volumetric water content of the compost feedstocks (m 3 m −3 ), TP or = Total porosity, ρ w = density of water (999 kg m −3 ), and ρ p = particle density of the feedstocks (Das & Keener, 1997). Das and Keener (1997) related the FAS, TP or , and V WC using the following equation which was applied in this study:

Relationship between FAS, V WC , andTP or
2.3.11. Moisture Content (%)(MC wb) Water density and Free Air Space relationship with BD Moisture contents were worked out before mixing of the feedstocks to retain a 58% MC wet basis and C:N of 28.8:1 gravimetrically at the Kenya Agriculture and Livestock Research Organization (KALRO), using the methods utilized by  and as highlighted in ). The C/N and pH ratios are as recommended to optimize the microbial action ( Density of water (ρ w ) (kg m −3 ): density of water at standard conditions of 999 kg m −3 was adopted.

Free Air Space (FAS): In this study, the general regression equation which is equivalent to
The following regression equation was used to work out FAS: Equation 11 was used by ) and assumes a constant particle density at all moisture contents under any given compressive force. The relationship between the bulk density and FAS of composting materials is linear (Yu et al., 2009).

Particle density (kg/m 3 ) PD and Permeability (K)
The mean particle density used in this study was 1350 kg m −3 adopted from the ranges as indicated in Table 1 of (Das & Keener, 1997, p. 276).

Permeability (K):
In this study, the relationship based on the physical characteristics of the feedstocks used was utilized to calculate the airflow parameters such as permeability using Ergun (1952)'s work that related the permeability (K) and the possibility (η), that are derived from the Kozeny-Carman model equation (Das & Keener, 1997;Richard et al., 2004;Siyu et al., 2005).
Given, K and η = permeability (m 2 ) and passibility (m), respectively. d p = particle diameter (m) which is an effective particle diameter (weighted average surface to volume ratio) (Macdonald et al., 1979) was measured at KALRO (0.0155 m), A and B = Ergun viscous component constant and inertial component constant, respectively (A = 180 and B = 1.8) as revised by Macdonald et al. (Richard et al., 2004).

Air volumetric and mass flow rate:
Since the volumetric flow of air through a porous media (compost feedstocks inclusive), assumably flows in accordance to Darcy's law, provided the pore flow velocities remain low and maintain the laminar flow, the volumetric flow was calculated in this study using (Lynch & Cherry, 1996;Poulsen & Moldrup, 2007): Where _ V is the volumetric flow rate (m 3 s −1 ), K is permeability (m 2 ), A is cross-section area of the CR (m 2 ), ΔP = pressure drop, L is length over which pressure drop is taking place, μ is viscosity of the air, and ρ is air density.
With all the physical parameters such as temperatures for the CR (T CR1 and T CR2 ) and COHE (water inlet T COHE1 ), density, specific heat capacities, mass flow rates, etc., known, calculating the design parameters followed.

COHE Thermal Fluid Flowrates
The catch and weigh method was used for measuring the flowrates. A two hundred and fifty millilitres (250 ml) bottle was used and the time taken to fill it up measured on 10 runs basis to find the average. This method was used by (Burleson et al., 2020), during the computational modeling and empirical analysis of a biomass-powered drinking water pasteurization technology.
Two Mini DC 12 V 3 M micro submersible brushless water pumps with 5 watt consumption rated at 240 l/hr static flow, were used.
Polyethylene (PE) foam insulation material tightly wrapped around the fittings and connections between the compost reactor and the Water Tank (WT), improved heat insulation. One thermocouple mounted inside the supply pipe from the pile to the WT, monitored the temperature of the water.
The average water temperature was monitored and viscosity and density considered assuming a not more than 45 °C temperature rise.

COHE Surface Area
The surface area of the COHE was to be within the ratio of 1.5:1 to 3:1 (CR/COHE). This was to conform to the data reviewed from the literature were different researchers monitored, utilized the hydronic systems within these ranges (Compost mould/Core ratio) (Biomeiler, 2020; Pitschel & Lowry, 2016;Spade, 2014).

COHE mathematical modeling (design calculations)
The Kern's method (Donald, Q Kern, 1983) and the engineering principles presented by Shah and Sekulic (Shah & Sekulic, 2003) and used by Irvine, were applied to design the COHE used in this study. Further, the step by step design procedures for the heat exchanger outlined and used by different designers were followed (NPTEL, 2019;Shawabkeh, 2015). The process and procedure adopted are shown in Figure 3.
2.5.1. STEP 1: Thermal and hydraulic design 2.5.1.1. Thermal design. Thermal design of a heat exchanger involved the consideration of many interacting design parameters summarized as process and mechanical parameters.

Process parameters.
The process parameters considered are as follows: • Thermal fluid assignments to CR, which was the hot side (shell side), and the COHE which was the cold side.
• Average temperature specifications selection derived from the CR experimentation and water.
• Setting CR (shell side) and COHE (tube side) pressure drop design limits.
• Determination of heat transfer techniques and fouling coefficients for CR (shell side) and COHE (tube side). The properties considered for the design for the water, CR Feedstock (FS), and the transport medium are: temperature flow rate (ṁ) (kg s −1 /kg hr −1 ), density (kg m −3 ) and specific heat capacity C p (kJ kg −1o C −1 ) (for water and the CR) and the thermal conductivity of (k w ) of 0.50 W m −1 °C −1 (ET, (Engineering ToolBox), 2011). The property details for the water and CR average C p , and the HDPE are presented in Table 1 and Table 2 on the results sections, respectively. These were considered at the caloric temperature (Kern, 1983; NPTEL, 2019).

Determination of the Thermal Fluid
2.5.3.1. COHE (Tube side). The thermal fluid used in the COHE was water due to its availability at a lower cost and the good thermal properties, good enough to absorb and store thermal energy (Allen & Chambers, 2009).

STEP 2: Energy Values evaluation (Energy balance)
Since the CR inlet and outlet temperatures, specific heat capacity, and mass flowrates, were known, it was easier to calculate the missing parameters for the COHE (outlet temperature of the water (T COHE2 )). The inlet water average temperature (T COHE1 ) was evaluated prior to the tests.

Heat duty.
Heat duty is the amount of energy the heat exchanger must transfer to the fluid used in the process to heat or cool to the required temperature. In this study, it involved the heat lost by the hot fluid (the CR) denoted by Q CR and the heat gained by the cold fluid (COHE water) denoted by Q COHE . The principals of the law of conservation of energy complimented by the first law of thermodynamics were used to compare the energy generated by the shell side (CR) and assumed that it was gained by the tube side (COHE) (Zohuri, 2018), because energy is neither created nor destroyed. Previous researchers worked out the energy balance in heat exchangers using this method (Bergman et al., 2011;Irvine et al., 2014, May;Sekulic, 2020;M. M. Smith & Aber, 2014) and is expressed as: Where Q CR was the heat energy lost by the CR and Q COHE was the heat energy gained by the COHE.
It should be noted that in this study, no phase change occurred; therefore, the formula used is for sensible heat transfer (Shah & Sekulic, 2003) and the analysis was subject to the following assumed conditions (Theodore, 2011): • The heat exchange is only between the hot and cold thermal fluids (Insulated heat exchanger) • Neglecting the axially conduction along the tubes • Energy changes (potential and kinetic) were negligible • Under constant specific heat capacities of thermal fluids in the CR and COHE • Under constant overall heat transfer coefficient Heat lost by CR Where Q CR is the total heat lost by the CR (kW), _ m CR is the mass flowrate of the CR in kghr −1 heat stream underpinned by the moist air (0.055 kg/s or 198 kg/hr); Cp CR is the average specific heat capacity of the feedstocks used in the CR (3.03 kJ kg −1 °C −1 ). (T CR1 -T CR2 ) is the temperature difference between the initial (T CR1 ) (65.56 °C) average highest temperatures attained during standalone CR tests to the average lowest assumed to avoid CR inhibition (T CR2 ) 53.33 °C.
Heat gained by the COHE: From equation 16, the energy lost by the CR was assumed to be the energy gained by the COHE as shown in Equation 17: Where Q COHE is the total heat gained by the COHE (kW), _ m CR is the mass flowrate of the COHE in kg hr −1 heat stream underpinned by the moist air (0.025 kg s −1 or 90 kg hr −1 ); Cp COHE is the specific heat capacity of the water in the COHE (4.175 kJ kg −1 °C −1 ). (T COHE1 -T COHE2 ) is the temperature difference between the initial (T COHE1 ) (23.36 °C) and the calculated T COHE2 . The Mass flowrate ( _ m COHE ), Cp COHE and the T COHE1 , of the COHE were known and therefore it was easier to work out T COHE2 from Equation 17 (Haug, 1993) and applied by G. Irvine, using the standard heat flow into a material with constant pressure.

STEP 3: Estimation of Overall heat transfer coefficients for the Tube (COHE) and Shell side (CR)
In this study, the range of the initial assumed overall heat transfer coefficient (OHTC) U ass , was 11.34 W m −2 °C (Appendix A), considered the moist air at low pressure from the CR as the hot fluid and water inside the tubes as the cold in a cooler type heat exchanger (Edge, 2020;Sinnott, 2005). The assumed OHTC was in line with the recommendation for low-pressure coolers with fluids at atmospheric pressure, liquid outside or inside and (gas) moist air at low pressure inside or outside the tubes in a cooler type heat exchanger (Edge, 2020;Sinnott, 2005). The assumed OHTC was within the 30% threshold recommended Kern's limitations and used by Abd et al. (Abd & Naji, 2017;Kern, 1983;Towler & Sinnott, 2008). The OHTC assumed was in line with the ranges achieved by previous researchers (Mudhoo & Mohee, 2007;Sylla Boundou et al., 2006). This process was key in helping in calculating the area values using the mean temperature difference (Abd & Naji, 2017). The overall heat transfer coefficient depends on the tube side (COHE) and the shell side (CR) individual heat transfer coefficients and fouling resistances (Dhavle et al., 2018;ET, (Engineering ToolBox), 2003;Kern, 1983).

STEP 4: Determination of the tentative CR (shell)-COHE (Tube passes-n p )
In this step the tentative number of shell (CR) and tube (COHE) passes (n p ) and the Log mean Temperature deference (LMTD) and its correction factor (F) were determined. For the steady and efficient operation of the heat exchanger, the (F), was to be higher than 0.75 (Donald, Q Kern, 1983;NPTEL, 2019). In this study, one shell and multiple tube pass configuration were considered.
2.5.6.1. Logarithmic Mean Temperature Difference (LMTD). The LMTD just like the Effectiveness equation is specific to the heat exchanger fluid flow direction (Bergman et al., 2011). The following equation was used to find the LMTD: For counterflow, the LMTD was given by: In addition, for the Parallelflow, the temperatures were worked out using: In this research, the counterflow arrangement was selected.
2.5.6.2. Log-mean Temperature Difference Correction Factor F. The correction factor (F), is the ratio of the True (effective) Mean Temperature Difference (TMTD) (ΔT m ) to the LMTD and referred to as the log-mean temperature difference correction factor, or exchanger configuration factor. This can also be found by the ratio of the actual heat transfer rate of a heat exchanger to that of a Counterflow heat exchanger with equal UA and fluid terminal temperatures. F is dimensionless and dependents on the effectiveness of the temperature, heat capacity ratio, and the direction of fluid flow of the heat exchanger in consideration (Shah & Sekulic, 2003). The chart for one to two shell pass and two, four, six (multiples of 2) tube passes (Appendix B) was used to configure and compare the mathematically calculated F. To compare, the two dimensionless temperature ratios were calculated as follows: And The above equations for R and S were also applied to calculate the Heat Exchanger Thermal Efficiency (HETE) where the relationship of ε and NTU data was used for conditional flow arrangements. These methods entailed the calculation of the dimensionless parameters S and R, which are determinants of the correction factor (F) by using the following equations.
F was given by: Where TMTD was calculated using: 2.5.6.3. Heat capacity rates. Individual thermal fluid heat capacity rates were worked out using; Where C is the fluid heat capacity, _ m is the mass flow rate of the fluids and C p is the specific heat capacities of the fluids. In this study C min was C COHE standing for the heat capacity of the COHE and C max , was the heat capacity of the CR.
Heat Capacity ratio (C*) was given by; 2.5.7.
Step 5 Heat transfer determination The energy calculated in step 2 using Equation 16, relates with the heat transfer area, Overall Heat Transfer Coefficient selected in step 3 and the T m in Equation 24.
Where Q is the total heat transfer rate, U ass is overall heat transfer Coefficient (Wm −2 K −1 ), A o is the total heat transfer area (m 2 ) and T m is the log-mean Temperature Difference, d o is the outside diameter of the tube and L is the length of the pipe in meters (LMTD) (Abd & Naji, 2017;Ezgi, 2012 Where A o is the provisional total outside surface area of the tubes and A ot is the outside surface area of one turn of the tube.
D b = bundle diameter in m, p t = pitch, d o = tube OD in m, N t = number of turns, K 1 and n 1 are constants from the number of passes and triangular pitch.
2.5.8.6. Selection of the clearance between the shell side and the bundle. In this study, the fixed and U-tube were chosen from the bundle and shell diameter chart in Appendix D (Abd & Naji, 2017). The shell diametrical clearance was selected and applied in the calculation of the clearance of bundle diameter as follows:  , 1995). Mean temperature of the water was worked out; Total flow cross-section area of the tube in this study was worked out using; A at = total flow cross-sectional area.

Water mass velocity (G t );
Water linear velocity (u t );
Water Coefficient: the following equation given from the data of (Eagle & Ferguson, 1930) was used in this study to calculate the water coefficient; Where h i = water coefficient-inside (W m −2o C), t = water temperature ( o C), u t = water linear velocity (m s −1 ), d 1 = inside tube diameter (m), k f = conductivity of the fluid (W m −2o C), Re = Reynolds number (dimensionless). Pr = Prandtl number (dimensionless), μ = viscosity of water (0.0008 Ns m −2 ), j h = heat transfer factor from (Sinnott, 2005, p. 665), and μ w = viscosity of the water at wall temperature (N s m −2 ).

Prandtl number calculations;
Finding the heat transfer factor: The heat factor for this study was worked by dividing the pipe length (L) in m by the internal diameter (ID) (d i ) in m: The tube side heat transfer factor table was used to find the heat transfer factor (Abd & Naji, 2017;Sinnott, 2005, p. 665).

Heat transfer coefficient-Shell side (CR).
Baffle Spacing: the baffle spacing (l b ) used was selected to be as close as possible to give higher heat transfer coefficients (Hayati, 2014;Kern, 1983). Determination of cross-flow area: The cross-flow area (A s ) of the hypothetical row at the middle of the shell (equator), was calculated by: A s = cross-flow area, p t = tube pitch depicted by distance between centers of two tubes (P t = 1.25d o ), d o = tube outside diameter (m), D s = diameter of the shell inside (m), l b = baffle spacing (m).

Determination of mass (G s ) and linear (u t ) velocity of the shell side:
And, W s = fluid flow rate (kgs −1 ) on the shell side, ρ = density of the fluid in the shell side (kgm −3 ).
Finally, the shell side (CR) heat transfer coefficient was calculated by:

STEP 7: Pressure drop
The equation offered for calculating the pressure drop was:

Tube side
Or, Where L = tube length Find j r from the tube side friction factor graph (Hayati, 2014, p. 39) using the Re number found in Equation 41 (Kern, 1983).

Shell side:
Then used the Re in Equation 50 to find the j r from the shell side friction factor segmental baffles graph (Hayati, 2014, p. 40;Kern, 1983).
Pressure drop was given by:

STEP 8: Overall heat transfer coefficient
The overall heat transfer coefficient of the COHE pipe was calculated by adding the stand-alone heat transfer coefficients of the thermal fluids involved using Where o = outside pipe wall (m), i denotes the inside of the pipe wall (m), and U represents the overall heat transfer coefficient (W m −2 o C −1 ).

Experimental setup (Mechanical structure design and Fabrication)
Mechanical Parameters: The mechanical parameters considered were: • The COHE Tubular Exchanger Manufacturers Association (TEMA) layout and the number of passes.
• Lower and upper limits on the length of the tube.
• Design parameters of the CR.
It should be noted that the system was designed to operate at a steady state (Edwards, 2008;Sekulic, 2020;Sölken, 2020).
This section addressed the fabrication of the exchanger based on the mathematical calculations achieved. Fabrications were done according to the design drawings depicting the mathematical calculations and presented in Figure 4. The heat source fluid (air) was influenced by a passively aerated system that moved from the bottom to the top based on the chimney effect as shown in Figure 2. The CR design and functions are based on .
In Figure 4, the Cold fluid (water) flow direction is shown. It entered the COHE from the top and flowed back to the water tank through the bottom. The water, therefore, flowed in a counter direction to the Heat source resulting in a Counterflow heat exchanger (Theodore, 2011).

Performance Analysis Formulation
In this study, the performance indicators of the COHE were performed using the ε-NTU method (Theodore, 2011). To compare the Temperature efficiency between different heat exchangers such as, counterflow, crossflow, and parallelflow the, Log Mean Temperature Difference (LMTD) method was applied. These methods were used by (Adumene et al., 2016). Further step by step heat exchanger design aspects were followed (Edwards, 2008;Shawabkeh, 2015).

Effectiveness of the COHE
The effectiveness of the COHE was worked out using the heat exchanger effectiveness equations that are specific to the heat exchanger fluid flow direction. In this study, the counterflow system was used, based on the flow of the fluids in the CR/COHE system. This method was also adopted by and covered widely.
Where ε is the effectiveness of COHE, Q is actual heat transfer, and Q max is the maximum possible heat transfer. Equation 59 is further identified by the following equations: Where C* is heat capacity ratio expressed in Equation 25.
In addition, Q max is simplified in the following equation: Equation 60 and 61 results in the effectiveness of the counterflow with the equation; The final calculation was to find the COHE water outlet temperature using the following equation:

Number of Transfer Units (NTU) and Temperature Correction Factors
The NTU was given by; NTU ¼

UA C min (63)
Where U ass is the overall heat transfer coefficient, A o is the area of the heat transfer surface, C min minimum heat capacity between the thermal fluids, and NTU is the Number of Transfer Units (Shah & Sekulic, 2003).
The association of the LMTD and ε-NTU was key in carrying out assumptions using P and R values to calculate heat capacity ratio (C*) Equation 25 and effectiveness equation 61. If R ˂ 1, R = heat capacity ratio (C*) Equation 25, and P = ε (equation 61 and equation 62).
The above assumptions were used by (Guimaraes et al., 2015) in the numerical determination of the LMTD correction factor for shell and Heat Exchangers. They are further illustrated by (Shawabkeh, 2015)

Effectiveness and NTU relationship
The ε-NTU relationship for the counterflow (Equation 61), crossflow (single pass) both fluids unmixed (CF(sp)BFU) (Equation 65) and Parallelflow (equation 64) at known NTU and C*, were evaluated using the different effectiveness equations as follows (Theodore, 2011): • Crossflow (single pass) both fluids unmixed The objective was to evaluate the heat exchangers with the better ε under the standard NTU and C*.test conditions used in this study.

Statistical methods
Temperature profiles were captured by means of the Polytetrafluoroethylene (PTFE) thermocouples connected to the TC TM -8 Picolog Data loggers with a Temperature accuracy of ±0.2% of reading and ±0.5 °C, were used to transmit the data to the PC via a 15 m USB cable. The PTFE thermocouples were chosen because of the robustness to withstand harsh conditions within CR due to humidity, pH changes, and microorganisms actions resulting in chemical reactions. Bajko et al. (2018), designed a wooden rod that was inserted into the Polypropylene to avoid corrosive environmental effects on the sensor. Testo 174 H and FreeTec data loggers as were also used to capture data.
The temperature data were then downloaded to the PC and using OriginPro 2018, which uses many nonparametric tests such as Friedman ANOVA and two-sample Kolmogorov-Smirnov test to confirm the normality of the data P values range (P˂0.005).

Computational Model
The Chemical Engineering Calculations Hydronic Compost Heat Exchanger design and performance evaluation software here referred to as HYDROCOHE software was created to be used in this study. A computation model was created in the CheSheets a free web-based spreadsheet program developed by the Chemical Engineering Calculations (CheCal) and distributed on an as-is basis and free software, that provides step by step guidelines on designs of a chemical engineering system. A Google sheet was developed and coded to return output values based on specified input parameters under the Compost reactor situation using the Visual Basic for Applications (VBA) skills (Chelcal, 2018).
The HYDROCOHE input values were based on the Fluid and the Exchanger data. The Fluid data composed of flowrates, Hot-for CR (ṁ CR ) and Cold for COHE (ṁ COHE ), input temperatures (T CR1 and T COHE1 ) and the specific heat capacities (C pCR and C PCOHE ). In the exchanger data, the type of the fluid flow arrangement was considered (Counterflow, parallelflow, and crossflow all fluids unmixed) and the overall heat transfer Coefficients (U).
The outputs were based on the calculations coded using VBA from the mathematical equations discussed in this study, and the outlet temperatures, (T CR2 ) and (T COHE2 ). The calculations achieved were; the heat transfer conductance, the NTU, the Effectiveness (ε), the heat capacity ratio (C*), Heat transfer or energy rate (Q), the mean temperature difference (MTD), the LMTD, and the LMTD Correction factor (F). The outlet temperatures obtainable were for the CR (T CR2 ) and the COHE (T COHE2 ). The ε-NTU curves were done.
The default model was predicating initially on counterflow direction values. During the comparison tests, the input configuration was changed to the crossflow and parallelflow heat exchangers flow sets. The results of each heat exchanger's performance were correlated with each other, with the mathematically obtained values.

Computation and Empirical values Comparisons
After experimental trials of the CR and COHE assembly, the average temperature, flow, and energy values obtained, where simulated into the HYDROCOHE system to output the values and compare. This was important to validate the experiment and the data obtained.

Ethics statement
The collection of the samples was permitted by the owners of the farm in the Malidadi area of Kitale at the geographical location of Latitude 01° 08ʹ 39 N and Longitude of 034° 59ʹ48E and elevation 1907 m. The chicken manure (CM) and the sawdust (SD) were voluntarily provided by a local chicken farmer and sawmill owners, respectively. The National Commission for Science, Technology, and Innovation (NACOSTI), a Government statutory board in Kenya, issued the research license number NACOSTI/P/18/99,643/27,086. Therefore, all experimental procedures conformed to the regulations established by the NACOSTI.

Experimental analysis of the CR
The design and achieved parameters of the CR used in this study are presented in Table 1 and  Table 2 summarizes the calculated physical parameters.

Summary of feedstocks and CR parameters
The feedstocks used are, cow manure, green farm weeds, maize cobs as bulking materials, maize stovers, and sawdust. The C/N and C p calculations are highlighted in . The quantities and physical parameters are as summarized in Table 1.
The composting process is controlled by environmental factors (temperature, pH, moisture content, and aeration) and Feedstocks/substrate natural factors (C/N ratio, nutrient content, and particle size) (Diaz et al., 2007;Makan et al., 2014). Therefore, most importantly, in this study was the balancing of the heat generation from the compost reactor and the efficiency of the heatloving microorganisms responsible for aerobic composting stabilization. Previous studies list heat release as an indication of the aerobically composting process stability and are directly associated to the air supply for the process of aeration (Notton, 2005).
The physical properties considered in this study and highlighted by other researchers are highlighted in Table 2 Notton, 2005;Walling et al., 2020). Air supply to the composting system has threefold functions; Oxygen supply line to the microorganisms, excess moisture removal, and the dissipation of excess heat generated (Diaz et al., 2007;Makan et al., 2014;Notton, 2005). The design and achievement of a good quantity of air supply and its efficiency is linked to a complete understanding of the physical properties, processes, and materials involved.
The temperatures achieved were in the range of 66.8 and 18.4 o C maximum in the CR and ambient, respectively. The ambient air density calculated was 1.20 kg m −3 . Upon entering the CR and at the maximum temperature of 66.8 °C, the air density calculated was 0.96 kg m −3 . This reduction of density as the result of the heat generation from the microbial activities causes the air to rise (weight reduction) to the top leaving a space thereby creating the chimney effect, that aids the inflow of air from the bottom of the CR. Furthermore, the condensate was observed on CR top plastic cover, signaling the condensation of the evaporated water from the compost which further points to air density reduction, thereby exacerbating airflow (Poulsen, 2013).
The drop in density from 1.20 to 0.96 kg m −3 in this study, created a pressure difference of 1.99 Pa. Airflow rate is further affected by the permeability (Barrington et al., 2002;Das & Keener, 1997). Achieved in this study was a permeability of 6265.63 µm 2 , which is in the range of the reported ranges of 2500-25,000 µm 2 and 700-8000 µm 2 for biosolids and cow manure, respectively, over a compressive pressure stress ranging between 0 and 20 kPa. Air permeability is important in that it is indirectly proportional to the stress exerted on the CR (J. M. . The air porosity achieved in this study was 82.38% and was good enough to allow the airflow resistance decrease and permitted the flow and was in the range of the reported and recommended (J. M. . At 82.38% porosity, the CR was provided with enough air and water-filled voids to allow adequate water and air availability to ensure efficient microbial activities at the particle size of 0.0155 m. This resulted in the FAS of 49.56% that is in accordance with the reported values, good enough to command easy air and water flow through the CR feedstocks and had great influence on the heat and mass transport, and microbial kinetics (Alburquerque et al., 2008).
The permeability (K) is dependent on the FAS and the particle size as shown in Equation 12. The superficial velocity also depends on the fluid kinemics viscosity of the air as shown in Equation 1. The superficial velocity achieved in this study was 0.024 m s −1 and a mass flow rate of 0.0055 kg s −1 . Similar figures have been reported (Alkoaik et al., 2019).
The achieved parameters in this study qualified the flow described by Darcy's law since the Reynolds number was less than 1 (0,991) (Poulsen, 2013;Siyu et al., 2005;Yu et al., 2006). The good temperature profiles achieved, were good enough to proceed with heat extraction through the designing of the COHE.
It should be noted that trials on the operation of the Compost reactor were conducted separately before the extraction process was instituted and not reported here. The heat bulb (Exchanger) thermal fluid used was water. The storage water tank was 36 liters capacity and raw water from the Kitale water utility company used. The CR with the above parameters was the heat generator.

Summary of achieved CR and COHE thermal fluids experimental parameters
The initial average temperatures (Input) for the CR (moist air) and COHE fluids achieved were 65.53 °C and 23.36 °C, respectively. In order to prevent the cooling of the CR, the minimum temperature drop expected was 53.33 °C. The assumed overall heat transfer coefficient (U ass ) was 11.34 and was in line with the recommendation in Appendix A (Low-pressure coolers), liquid outside or inside and (gas) moist air at low pressure inside or outside the tubes in a cooler type heat exchanger  (Edge, 2020;Sinnott, 2005). The assumed OHTC was within the 30% threshold recommended Kern's limitations and used by Abd et al (Abd & Naji, 2017;Kern, 1983;Towler & Sinnott, 2008).
Upon having the CR and COHE thermal and flowrates parameters, designing calculations were possible and results are presented below.

Mathematical methods
In Table 4, the cold and hot side physical parameters achieved, are presented.

Energy balance
An energy balance for the two fluids (Shell =CR and Tube = COHE) was achieved using Equation 16. It was assumed that heat loss to the surrounding was negligible, with no potential or kinetic energy changes, no phase, and specific heat capacity changes, and at steady-state conditions (Kuppan, 2000;Mujumdar, 2007). The data and results used to perform the energy/heat transfer and calculated T COHE2 is presented in Table 4.

Design and Performance analysis using Kern and ε-NTU/LMTD methods
The design and performance analysis was achieved by using Kern and ε-NTU methods (Abd & Naji, 2017;Towler & Sinnott, 2013). The inlet temperatures were 65.53 °C for the CR (T CR1 ) and 23.36 °C for the COHE (T COHE1 ).

Mechanical Design: Type of heat exchanger and dimensions.
In this study, a one-shell pass and 56 tube passes counterflow heat exchanger arrangement, was adopted. The LMTD method was applied to calculate the required area of the COHE and detailed results are presented in Table 5.
As indicated, the correction factor (F), was 0.9566, above 0.75 recommended for steady operation of the heat exchanger, thereby validating the efficiency of the configuration applied in this study (Ali & Naji, 2017;Kern, 1983;NPTEL, 2018;Towler & Sinnott, 2013)     The layout used in this study is outlined in Table 5 (0.01905 m-OD and 0.01605 ID). The thermal conductivity of the material used in this study (HDPE) 0.50 W/m/ o C, was lower than the ones in metals such as copper (385 W/m/ o C) or steel (50.2 W/m/ o C) (ET, (Engineering ToolBox), 2011). The higher the thermal conductivity of a material, the better the conduction of heat (Theodore, 2011).

NUMBER OF TUBES
The actual heat transfer area (A o ) was 7.2 m 2 . It was noted in this study that increasing the overall heat transfer coefficient resulted in the reduction of the required total heat surface area of the heat exchanger. Abd and Naji (2017), indicated that the reduced heat transfer area in heat exchanger designing, translates into reduced equipment costs for both manufacturing, fouling, and maintenance costs. In a heat exchanger like this one where the shell side was the heat generator also depending on the microorganisms through composting, care was taken to ensure that the COHE/CR area ratio was not below 1.5 to avoid inhibition.

Tube and Shell side parameters Analysis
The tube side fluid velocity calculated in this study was 1.73 ms −1 at the Re number of 51,609.93, resulting in a turbulent flow. On the other hand, the shell side velocity was 0.18 ms −1 at the Re of 23.8 qualifying it into a laminar flow. For Re ˂ 2100 is laminar flow, 2100˂Re˂10,000 is transitional and Re ˃10,000 is turbulent flow (Theodore, 2011, p. 288).
The baffle size used in this study was 25% cut which is widely used and within 0.2 D s . All dimensionless figures were considered as highlighted in the methodology section and Table 2.
The shell side (CR) heat transfer coefficient achieved in this study was 13.96. Previous studies have reported heat transfer coefficients in the range of 13.55, 17.50, 17.46, and 23,95 W m −2 °C −1 (Sylla Boundou et al., 2006). Mudhoo and Mohee (2007) reported that the overall heat transfer coefficients during composting of organic substrates aerobically are proportional to the temperatures. (The higher the temperature the higher the OHTC.) They found maximum OHTC between 35.5 and 263.9 W m −2 °C −1 , and minimum between 2.44 and 8.15 W m −2 °C −1 . Consequently, the tube and shell side heat transfer coefficients were calculated resulting in the overall transfer coefficient of 13.48, 18.9% of the assumed, and within the recommended 30% (Bergman et al., 2011;Kern, 1983;NPTEL, 2019). The pressure drops were also within the required ranges making the calculations valid for further tests and validation.

Summary of COHE design parameters
The summarized design parameters for the COHE are shown in Table 5 and were used in the fabrication works.

Fabrication and site assembly
The COHE was fabricated following the parameters in Table 5. The HDPE pipe was weaved onto the frame that acted as a baffled stage as shown in Figure 5. It was inserted into the CR as presented in .

Experimental Set-up
The experimental set-up was achieved using the parameters and procedures explained in .

Energy and Temperatures Values achieved
4.2.2.1. Temperature profiles. The temperature profiles achieved by the CR and COHE are shown in Figure 5. On day 1.88, 60 °C was achieved, the highest of 68.9 °C (Point B 1 ) recorded from a starting point of 25.5 °C. This temperature/duration achievement was reported by (Shimizu, 2017) who achieved similar values in 2-4 days. Similarly (Allen & Chambers, 2009) reported 65 °C while composting Horse-based manure from an in-vessel tunnel composting facility in Scotland.
The extraction of heat from the COHE started on day 1.88 when on average, CR was at 65.56 °C, water in the tank (T COHE1 ) at Z was 23.36 °C, T COHE2 at X, was 33.9 °C and the ambient was 20.06 ° C. A to B (T COHE2-T COHE1 ), was at Zero before extraction started. Upon switching on the water circulation pump, the T COHE2 went up to 58.1 °C. This jump in temperature rise is attributed to the reserve water that was in the COHE system (34-liter volume). The water in the tank went up 35.1 °C at this level. However, within 5 hours of extraction commissioning, the temperatures sharply went down to 40.4 °C in the COHE (T COHE2 ), 33.7 °C in the water tank and 61.9°C in the CR, attributing the drop to the hot and cold water mixing. Between day 1.88 and 3.5, the temperatures remained constant between 40.1 o C snf 42 o C in the T COHE2 and 30°C and 38 °C in the T COHE1 .
The temperatures from the COHE went down from 41 to 36.2 o C by 07:28 am on day 3.6. At these time set the T COHE1 temperature was at 34.9 °C, ambient was 16°C and around an average of 51 °C from the compost reactor.
Allen and Chambers (Allen & Chambers, 2009), reported similar temperature drops where an initial activation of the heat extraction using a compost heat extractor with the fluid mixed with water and Glycol was used to heat the water from 5 o C to 48 °C, but dropped to 38 °C after 5 minutes. This signified that more heat has been extracted than what has been generated by the CR. This is the critical problem associated to the heat extraction systems as also reported by (Bajko et al., 2018) and was addressed by staggering the extraction rates and times.
To mitigate on the temperature drop of the in the setup, a 400/750 watts electric water heater procured locally, was introduced to the water tank to heat the system above 55 °C, a temperature recommended for the good composting process. Heating started at point C 1 in Figure 6 and continued for 37 hrs 29 minutes until the water tank attained 70.1 °C. The compost average temperature rose to 69 o C. The maximum COHE IN temperature was 70.1 °C and 68.7 °C for the COHE OUT , minimum of 23.4 and 33.9 °C, respectively. The total energy used by the electrical heater was 28.1 kWh.

CR/COHE and T COHE1 /T COHE2 Temperature differences
In Figure 6, the average temperature differences between the T COHE2 and T COHE1 , and T CR and the T COHE2 , are indicated. Points A and B indicate the period before extraction commenced and the difference between T COHE2 and T COHE1 was Zero. CR was above 65 °C. Areas C and D show the period when the electrical heater was operational. The net temperature was negative because the cold fluid (T COHE1 ), supplied more heat than the T CR . Upon turning off the electrical water heater at point D, the pattern of having the positive value for temperature difference resumed. The average difference temperature throughout the heat extraction period between the T COHE2 and T COHE1 4 ° C and between the T CR and T COHE2 , was 8.5 °C with the maximum and minimum of 32.5 and 31.7 °C and −4.5 and-0.12 °C, respectively.
The ambient temperature was 14.6 °C average. Points 2 1 -5 1 show the effects of solar irradiance on good sun days. The water tank received more heat from the sun resulting in a negative temperature between the T COHE2 and T COHE1 . This was also experienced and highlighted by (Allen & Chambers, 2009).
The analyzed temperature values are in agreement with a section of the approved and existing research work and therefore are credible enough to be applied for the subsequent energy calculations.

Generated and Extracted energy
Extracted energy increased to727 MJ (201.94 kWh) during the extraction trials for 205 hours representing 17.83% of the total generated energy (4,077.4 MJ). The net extracted was 174 kWh (0.99 kW), minus the energy used by the electrical water heater of 28.1 kWh. This resulted into 7.8 kg/kWh or 538 kJ/kg per initial mass compost weight matter and 1273.65 kJ/kg on dry matter and 58,600 kJ/m 2 . Similar extraction levels have been reported by (Seki & Komori, 1995) who extracted 54,550 kJ/m 2 representing 16-22% of 739,827.4 kJ generated using a lab-scale positive aerated cylindrical CR of diameter 0.58 m and height of 0.895 m. In this analysis, a small difference can be attributed to the poor thermal insulation of CR by a thin polypropylene fabric. Average recovery rates of 1159 kJ/kg dm for lab-sized systems, 4302 kJ/kg dm for pilot scale systems, and 7084 kJ/kg dm for commercial systems were presented (M. M. Smith et al., 2017). Therefore, the energy figures achieved in this study are not only reliable but also credible and comparable to other researchers.

Coefficient Of Performance (COP) of COHE
The COP of the COHE during the Extraction period was a maximum of 489.9 and an average of 74.22 and gradually reduced to 71.5 on day 13 as shown in Figure 7 8. Zero COP was registered, during the first 1.74 days the COHE was not activated (No heat extraction). The average energy used by the pump throughout the trial period was 2.7 kWh.

Figure 6. Temperature differences between COHE out & CR and COHE OUT & COHE IN.
The results are higher than (Allen & Chambers, 2009), where a maximum of 12.5 COP was achieved while extracting heat from the compost reactor using a Groundfos UPS 15-50 130ʹ flow pump rated at 55.2 w, 3.4 mhr −3 and 6 m maximum head as compared to the pump used in this setup of 5 w and 48 DCV power requirement. The different power consumption explains the difference in the COP. The less energy the auxiliary system uses (pumping/heating system), the higher the COP. 4.2.5.1. Performance evaluation. The performance analysis results are summarised in Table 3. The cooler fluids exist in counter flow heat exchangers, at the end of the system where the hot fluids join the heat. The hottest cold fluid temperatures greater than the coldest hot fluid temperatures have been reported in the counter flow heat exchangers as opposed to the parallel flow ones (Tawil, 1993).   (Δ 12.23 °C). In this case, where the heat capacity rate of the Hot fluid is higher than the Cold fluid, it is highlighted that the cold fluid would absorb more heat hence larger temperature change (Theodore, 2011).
The effectiveness-NTU method was applied to analyze the heat exchanger in this study because only the inlet temperatures were known initially. This method is recommended by (Theodore, 2011) and (Shah & Sekulic, 2003) and was applied to design and evaluate a waste energy extraction system from the Composting facility located in Scotland, UK and recorded, ε of 0.7370, NTU of 1.5135. This method is simple than using the LMTD that requires prior knowledge of the inlet and outlet temperatures in order to calculate the T IM .
If the heat capacity rate of the hot fluid source was less than that of the cold fluid, the hot fluid would have experienced the larger temperature change and cooled down to the cold fluid temperature (Theodore, 2011). This effect can cause the cooling down of the compost reactor from the thermophilic stage if the temperatures go below 45 °C due to microbial activities interference (Adams, 2005).
The capacity ratio (C*) was 0.6258 and was derived from dividing the smaller to larger heat capacity rate between the CR fluid and the COHE fluid and hence mass flow rate dependant. C*≤1 in this design. A balanced heat exchanger is one with C* =1. The fouling capacity (0.002 W/m 2 /°C) is the standard considered when dealing with a fluid with more than 50 °C temperature (R. K. Shah & Sekulic, 2003).
COHE effectiveness was 0.464. For specific U values and C min , the value of NTU is the surface area configuration. It also represented the ratio of the actual heat transfer rate (2.038 kW) to the maximum possible heat transfer rate (4.404.625 kW) of the COHE. The larger the NTU, the larger the heat exchanger (Mohamad, 2015). The achieved effectiveness in this study is lower than the one reported in of 0.7270 at the actual heat transfer rate of 9889.3897 kW and a maximum possible heat transfer rate of 13,603.883 kW.
The F value found in this study was 0.9566 with a percentage relative error of −4.34% from the recommended unit F for a true counterflow exchanger with maximum temperature potential and with two fluids logically arranged in counterflow. T m was 25.02 °C approaching the T Im (26.15 °C) (R. K. Shah & Sekulic, 2003). Therefore, 4.34% relative error represented the degree of departure of the TMTD from the LMTD and not a signal of high efficiency but an indication in the close operating performance of a heat exchanger with comparable operating conditions of flow rates and temperature inlets of the involved thermal fluids in a counterflow exchanger. Numerical correction factors for the LMTD for shell and tube 1-2 heat exchanger have been highlighted by (Guimaraes et al., 2015).
The LMTD of the Counterflow heat exchanger was 26.15 °C in this study (Equations 18/19). Calculating at the attained input and output temperatures for the CR and the COHE. This would translate to 22.71 °C if a parallelflow heat exchanger was operated instead. This indicates that in the Counterflow heat exchanger, 26.15 °C was the mean difference in the temperatures between the COHE and the CR fluids. This indicates that the Counterflow heat exchangers transfer more heat than the parallelflow heat exchanger when operated under the same conditions. The LMTD correction was 25.02 °C.

Computation Analysis of the COHE
The input values used in the HYDROCOHE software were; inlet temperatures (T CR1 -65.56 °C), the (T COHE1 -23.36 °C), CR flow rate (198 kg/hr), the COHE flowrate (90 kg/hr), the COHE area (7.2 m) and the overall heat transfer coefficient (11.34). The results are shown and discussed in section 3.4.1 where effectiveness comparative analysis of the three types of heat exchangers is explained. Figure 10 shows the input data and results.

Comparison of heat exchangers effectiveness at a given NTU and C*
The achieved C* and NTU in this study was 0.626 and 0.79, respectively. Table 7 shows the comparative parameters of the parallelflow, counterflow, and crossflow heat exchangers at the achieved C* and NTU in this study.
The objective of this section was to compare the heat exchangers in terms of the flow type that would achieve the highest Effectiveness. In Counter flow heat exchangers, the effectiveness is in the range of 0 ≤ε ≤1 (Bergman et al., 2011). The nearer the ε value to 1, the higher the energy, the heat exchanger transfers.
In Figure 11 -13, the modeling results of Effectiveness (ε)-NTU for the Counterflow, Crossflow (single flow) both fluids unmixed, respectively, are shown on the ε-NTU curves. In Figure 11-13 C r = C*. Table 8 shows the modeling results comparing the performance and behavior of the Counterflow; Crossflow (single flow) both fluids unmixed and parallel heat exchangers at the initial input temperatures (T CR1 and T COHE1 ).    Therefore, the following is the analysis and conclusion from the effectiveness (ε), NTU and C* relationship achieved in this study and charts highlighted by (Tawil, 1993): • With small values of NTU up to about 1.5, the value of the effectiveness increase greatly.
Effectiveness ranges from 0 to 1 and reduces with larger values of NTU (larger than 3). For example, at 3.0 NTU in Figure 13, the ε was 0.6, almost equal to ε at 6.0. Operating a heat exchanger with higher NTU (larger sized heat exchanger), is economically not feasible. This is because larger NTU does not translate into good effectiveness but a larger heat exchanger area.
• Comparatively, at a certain NTU and capacity ratio (C*), the highest effectiveness is registered by the Counterflow heat exchangers (0.4766 Figure 11) seconded by the cross-flow heat exchangers (both fluids unmixed) (0.4622 Figure 12). Parallelfow exchangers have the lowest effectiveness values (0.4430 Figure 13). Tables 7 and 8, through the mathematical and computational modeling have demonstrated this principle achieved in this study accurately.
• At NTU values below 0.3, the effectiveness of the heat exchanger is not dependent on the capacity ratio, as shown in Figures 11-13.  • At equal heat capacity rates of two fluids, the heat capacity ratio is 1 (boilers and condensers), the effectiveness is minimum. The maximum effectiveness is achieved at C* = 0 for a given NTU. C*, ranges between 0 and 1.
The larger the NTU number, the larger the size of the heat exchanger (Theodore, 2011).

Model Validation (Mathematical, Experimental, and Computational Analysis)
The Mathematical and computation analysis results between the Counterflow is here presented in Table 9. The average temperatures achieved using the HYDROCOHE compared to the Mathematical, had minor differences with an R 2 value of 0.99938 between the T CR2 outlet temperatures and the R 2 value of 0.9995 between the outlet T COHE2 temperatures. When validating the model, the difference between the experimental and predicted temperatures was considered by (Burleson et al., 2020) and adopted the R 2 value of 0.9992. Hence the range in this study is as what was previously applied, rendering this model valid.
For the mathematical and HYDROCOHE software-computed values, the heat capacity ratio had the R 2 value of 0.9987, NTU registered R 2 of 0.9962, R 2 of 0.9976 for ε, and heat transfer (Q) had 0.9992 R 2 and the R 2 value of 0.9957 for the heat transfer coefficient. In the research to Optimize heat exchangers by a Bio-Inspired Artificial Intelligence Methods, J. Krzywanski used a maximum error lower than ± 3% to validate the approach (Krzywanski, 2019).
In this study, the error was between 0.1 and 0.4%, hence validated. Therefore, the presented procedure can be adopted as a technique in the heat exchanger design and performance analysis in the CHRs utilization systems. It is not only an effective method than the existing approaches, but also an alternative method for the existing trial and error CHRs heat exchanger design and optimization approaches that are complex analytically and numerically, and depends on unstandardized high costing experiments.
In the empirical data, the heat capacity ratio was not affected by the mathematically calculated, because the specific capacities between the CR and COHE fluids did not change and the mass flowrates assumed constant. However, the Overall Heat transfer coefficient dropped to 3.883 W/ m 2 /°C giving a 0.3671 R 2 value. This explains why the energy extracted was giving an R 2 of 0.4848 compared to the mathematically calculated and the HYDROCOHE computed values. This can be attributed to the poor thermal conductivity between the CR feedstock materials (Huet et al., 2012, September) and the scaled & corroded LDPE piping surface due to the caustic environment created in the CR during composting as shown in Figure 12. It is therefore deduced that the higher the overall heat capacity coefficients, the higher the heat transfer in a heat exchanger (Alkoaik et al., 2019). Table 10 shows the comparison in the output temperatures between the T CR2 and the T COHE2 mathematically, computed and experimentally. The R 2 value between the Math and the Comp temperatures was 0.9993 and 0.9995 for the CR (T CR2 ) and the COHE (T COHE2 ), validating the tests. However, between HYDROCOHE and empirical, the R 2 values were 0.9269 & 0.9432, CR (T CR2 ) & (T COHE2 ), respectively. The CR fluid had higher outlet temperatures than the COHE fluid. This might backup the low extraction rate of the COHE experimentally. It can also be attributed to the high fouling rate of the COHE pipes and its lower thermal conductivity. Scientifically, when the transfer surface of the heat exchanger is fouled, the distance between the width of the two thermal fluids may be reduced and is correlated with reduced heat transfer rate and increased pressure drop (Theodore, 2011).
The power extracted experimentally was 0.99 kW as compared to the 2.10 kW computed by the HYDROCOHE and the calculated by the mathematical modeling of 2.042 kW. This is a difference of 51.47-51.52% (R 2 value of 0.4853-0.4848). This divergence between the simulation outcomes and the experimental values was observed by Nwanze and Clark (Nwanze & Clark, 2019), who by using a computer-based model in COMSOL Multiphysics TM , discovered that the empirical data and the Simulation collated in the initial stages of the composting process, but diverged during the thermophilic stage. This was attributed to the heat losses of the compost reactors. Other researchers such as (Mohee et al., 1998), equally confirmed the rapid drop in temperatures during the thermophilic stage of the composting process, hence diverging from the computation setup.
Ventilative heat losses of 70-95% and 36-67% of the total heat flux in full and laboratory-scale compost reactors, respectively, due to conductive, radiation, and convective losses, have been reported by Mason and Milke (Mason & Milke, 2005).  The effects of insulation on the performance of compost reactors with rotary shape have been presented. It is highly recommended that for efficient and quality composting, insulation is effective and necessary for laboratory-scale and small pilot scale CRs with an Area/volume ratio of less than or equal to 6. This is because heat loss is high as A/V increases. The larger the compost piles, the lower the A/V ratio resulting in minimal heat losses (Alkoaik et al., 2019).
In this study, the A/V ratio was 4.4 and only a 3 mm Polypropylene cloth was used as insulation. The power discrepancy of 1,052 kW reported in this research between the empirical data, and the simulated and calculated data can therefore be attributed to the CR's poorly thermal-insulated boundaries as depicted in Figure 12.
Further, low thermal conductivity coefficients of 0.31 W/m/ o C have been reported for a compost density of 600 kg/m 3 at 60 °C and is dependent upon the moisture content, particle size, bulking agent quality and depth of the CR (H. Ahn et al., 2009). The more compacted the CR is, the higher the bulk density hence the higher the thermal conductivity coefficiency (Huet et al., 2012, September).
The CR used in this study had 201.7 kg/m 3 density and within 60 °C, hence poor thermal conductivity coefficient cannot be ruled out, in contributing to the difference in the energy extracted.

Conclusions
The possibility of using the aerobic composting route of biomass wastes conversion to energy has been presented in this study and 175 kWh was extracted. Associated challenges to the designing of heat extraction systems from the CHRs and the mitigation are proposed by using the computerbased HYDROCOHE software. The relationship between the Hydronic compost heat exchanger's mathematical and HYDROCOHE computer-based, design, and experimental performance modeling schemes was discussed in this paper. The mathematical and Computer-based produced similar results with R 2 values between 0.9962 and 0.9992 on the effectiveness, heat capacity ratio, NTU, heat transfer (Q), heat transfer coefficient, and outlet fluids temperatures (Hot and Cold fluids).
The Mathematical and the HYDROCOHE computer model have shown that the Counterflow heat exchanger arrangement, transfers more heat at the same design parameters (heat capacity rates, NTU, heat transfer coefficients, area, and inlet temperatures), than the crossflow and parallelflow heat exchanger arrangements. Since the generation of heat in the CHRs is heat-loving microorganisms-dependent, parallelflow, heat exchanger arrangements are highly recommended because of their ability to maintain the outlet temperatures at lower temperature differences as compared to crossflow and conterflow arrangements. Using the Counter or crossflow arrangements might result in the inhibiting of the compost system. The conclusions drawn from the analysis of the data analyzed are: • The computed effectiveness values from the CHECALC-based software and that realized from the actual empirical model, are almost the same.
• The heat generated by the CR (between 30 and 80 %) was lost via conduction, convection, and radiation.
• The CHECALC-based design software differed in terms of the experimental Outlet temperatures.
• At NTU values below 0.3, the effectiveness of the heat exchanger is not dependent on the capacity ratio, as shown in the Figure 9 to 11.
• At equal heat capacity rates of two fluids, the heat capacity ratio is 1 (boilers and condensers), the effectiveness is minimum. The maximum effectiveness is achieved at C* = 0 for a given NTU. C*, ranges between 0 and 1.
• The larger the NTU number, the larger the size of the heat exchanger.
• Using the LDPE as a piping material affected the heat capture rate because of its lower thermal conduction capacity.
• Despite the energy-extracted differences, the designing and performance, the HYDROCOHE computer-based program proposed in this study is essential for CHRs extraction system designs.
• The flow arrangement and direction have a bigger effect on the amounts of heat extracted, health, and performance of the microorganisms and the temperatures reached during composting.
Further, more work is required to rigorously validate the composting system computation models and transform them into intentionally useful hand tools for accurate CHR design. Accurately designed composting and heat extraction, processes promises to be major contributors to renewable energy sources. Source: (Dhavle et al., 2018;Hayati, 2014;Shawabkeh, 2015)