Operation Adaptation of Moving Bed Biomass Combustors under Various Waste Fuel Conditions

: This paper analyzes a moving grate biomass boiler operating with three alternative waste fuels, including biomass pellets, wood waste, and refuse-derived fuel (RDF) from a combination of thermal, economic, and environmental perspectives. The focus of this paper is on system functionality adaptation to retroﬁt the current systems operational conditions. A one-dimensional numerical bed model integrated with a black-box overbed model was developed to carefully investigate the fuel bed’s thermal characteristics, as well as the boiler’s output. According to the results, the system operates more e ﬃ ciently under the biomass pellets feeding and annually generates 548 GJ heat, while it drops signiﬁcantly in other scenarios. The system was economically evaluated based on a 25-year life cycle cost analysis. Subsequently, an internal rate of return (IRR) of 36% was calculated for biomass pellets, while the value reduced by 50% and 27% regarding wood waste and RDF, respectively. The fuel cost was identiﬁed as the main contributor to the total life cycle cost of the heating system, regardless of which feeding fuel was utilized. A long-term environmental impacts assessment of the boiler operation emerged, to show how plant-based fuels can signiﬁcantly decrease the impacts of climate change that have originated from fossil fuel usage. The current study concludes that all the proposed scenarios are feasible to di ﬀ erent degrees, and can extensively beneﬁt a diverse set of energy sectors.


Introduction
Thanks to advancements in technology, in recent decades, renewable energy has been deemed a very practical alternative to using conventional fossil fuels to supply cleaner sustainable energy systems [1]. Biomass, a carbon-based renewable resource, is capable to replace fossil fuels with the least modification in conventional facilities compared to other renewable resources. Among all possible biomass conversion processes, yet direct combustion is the most commercialized one [2,3]. Other advancing technologies, such as anaerobic digestion (AD), hydrothermal carbonization (HTC), as well as hydrothermal gasification (HTG) are generally less commercialized, owing to poor infrastructure of raw fuel supply, high operational and capital cost, and mostly are in prototype-scale so far [2][3][4][5].
Moving bed, fluidized bed, and pulverized fuel combustors are identified as the commercialized combustion technologies to burning biomass fuels. Based on a certain application, each of the aforementioned systems can be a more fitting choice, as discussed elsewhere [6]. In the course of small-scale purposes, moving grate furnaces are one of the most highly utilized approaches to burning raw and treated biomass waste. While, in addition, within the small to medium scale applications, the moving grate furnaces require less fuel preparation.
Currently, research has been conducted in order to modify the combustion system efficiency, as well as the contaminant emission issue that are associated with these systems [7][8][9][10][11][12]. Throughout the literature, the majority of the developed biomass combustion models were exercised in micro-scale system analyses in order to understand the biomass combustion properties more comprehensively [13][14][15]. These models were rarely governed as a system analysis tool in order to evaluate the functional system outputs and required operation adaptation for more efficient combustion. For example, a numerical model was developed by Zhou et al. [8] to study the effect of variant excess air ratio on combustion efficiency. The authors concluded that, under O 2 enriching combustion, which is defined as an increase in air flow rate from the stoichiometric condition to a higher fraction; the optimum combustion improvement is achieved when the O 2 volume fraction raises from 21% to 25%. Meanwhile, the maximum bed temperature is restricted to 1400 K. Similarly, bed and overbed zones were integrated on the same 3D grid, while the bed zone modeled through a porous medium [16]. The model then characterized a flow pattern and a specific temperature, as well as, pollutant profiles in the combustion chamber. Flame temperature, species concentration, and ignition front velocity for a packed-bed biomass combustor under CO 2 -rich condition was analyzed using an integrated 1D bed model into a 3D overbed model via Karim et al. [17]. The CO 2 -rich combustion indicates recycling flue gas as a fraction of primary air flowing into the fuel bed. The ignition front rate sped up to almost double the rate with CO 2 -rich air feed conditions. A connected 1D bed model and 3D overbed model was developed to cope with how different quality fuel influenced operations [18]. Larger high-temperature combustion zones were observed as a result of higher quality fuel, and, thus, consequently higher system efficiency. The effect of biomass fuel compositions variability on combustion properties in a moving grate combustor was studied through an independent bed model coupled with a simplified overbed combustion model [19]. The results demonstrated a major impact of composition variability on ignition front rate, conversion time, and output heat flux, as well as a minor impact on the flame temperature. Existing models help engineers to design for a better biomass combustion, which is quite necessary for technology advancement. Meanwhile, system analysis using black-box model without paying attention to operational detail can be beneficial only for general decision making [20]. However, for many of the domestic and industrial facility's users, it is impractical to take advantage of such modeling applications. In other words, for example, it is crucially important for them to know how system operation should be readjusted when switching from variant feeding fuels for example from specific biomass pellets to wood chips. In addition, what the economic and ecosystem consequences of this fuel switching action are.
In this paper, functional thermal properties of a moving bed biomass boiler fed with three specific waste-based fuels: biomass pellets, wood waste, and refuse-derived fuel (RDF), are carefully determined through a routine of one-dimensional transient numerical bed model integrated with an overbed black-box model. In this study, instead of focusing on biomass combustion diagnosis and micro-scale combustion analyses, we propose the required adjustments to facilitate system operation settings when switching between the proposed fuels. Furthermore, to illustrate the long-term economic and environmental impacts of system operation with each fuel, life cycle costing (LCC) and life cycle analysis (LCA) are also conducted. It should be noted that moving grate boilers can concurrently burn different fuel types, such as those that are applied in this study. In consequence, the results of this article are highly informative for system users.

Modeling Procedure
The process of biomass combustion inside the packed bed deals with thermochemical reactions consisting of drying, devolatilization, and char combustion. To some degree, this is a complicated process that still has many unknown variables. The existing 1D, 2D, and 3D empirical models can predict the actual combustion process from an engineering standpoint.
Drying as an endothermic reaction has a significant impact on combustion. Many studies have proved that higher moisture content of biomass particles can diminish thermal performance and increase contaminant emissions [15,[21][22][23]. By analyzing the literature, where different drying models can be identified, Arrhenius and heat sink models have been applied in recent studies. The influence of moisture content of different particle shapes on drying behavior was examined using the Arrhenius and heat sink models by Khodaei et al. [24]. Through a lack of experimental data regarding particle de-moisturization, the heat sink model exhibited more compatible results that was opted for in this work. Pyrolysis occurs after the evaporation of the initial moisture content. This includes several interrelated parallel reactions that initiate at disparate temperatures and begin at a higher rate than drying. Here, the fuel mass of the bed profoundly diminishes while volumetric fuel bed reduction is insignificant. Varied pyrolysis models have been listed in Reference [25], although only a few of them have been recently applied in the biomass conversion models due to their relative complexity and inaccuracy. Among all, the model with three parallel reaction-terms has been dominantly utilized, as it has a trade-off between computational cost and precision, and, therefore, simulates the devolatilization process in this article. The most complex reactions occur during char combustion, in which both diffusion and thermal mechanisms compete to overtake the reaction. It is imagined that char reaction consists of four overlapping reactions [25], while, in most of the research, the fourth reaction, the reaction of carbon with hydrogen, is overlooked, due to its negligible effect on the whole process. For a detailed description of each subprocess and the corresponding models, Reference [26] is recommended.
Fuel bed simulation has been conducted through varying techniques in the publications. They can be categorized as empirical bed models, separate bed models, discrete particle models (DPM), and porous medium bed models. A more thorough explanation of each technique can be found through further research beyond the scope of this article [27]. The separation of fuel bed and overbed zone modeling is the most preferred strategy in biomass furnace modeling, in which it facilitates the sensitive and in-depth analysis of solid particle conversion. Instead of 3D modeling of fuel beds, 2D models (or 1D transient models) are usually governed, because fields gradient in the z-direction of bed zones are insignificant. According to the literature, this strategy can be implemented by utilizing variant approaches, including the single-particle model, the continuous medium model, and the particle-resolved model. The single-particle model delivers proper results for a fluidized bed system, because the proposed particles have uniform boundary conditions, while in the moving grate combustors, particles have non-uniform boundary conditions. The particle-resolved modeling has modified the problems with single-particle molding, however, the computational cost of this is quite high. In this article, the continuous medium model is governed to simulate fuel conversion, where it has a fair trade-off between accuracy and computational cost. It should be mentioned that this method could raise dispersion in results when large-size particles are deployed. Figure 1 visualizes the fuel particles arrangement under the actual condition and continuous medium modeling approach along with heat transfer mechanisms inside the fuel bed. From this figure, the main limitation of the continuous medium model is to neglect radiation between particles inside the fuel bed. The furnaces side walls are assumed to be insulated so that no heat leakage can occur from the fuel bed zone to the outside of the walls. Reactions happen in the vertical direction inside the fuel bed, whereby the horizontal gradient is insignificant. The one-dimensional transient model is executed, in order to eliminate the complexity of two-dimensional modeling. This method, the so-called walking column, tracks changes in a tiny column of the fuel bed from the opening of the grate until all fuel in the column is consumed [28]. The column is split into small cells, as each cell consists of the solid and gas phase and is updated during the conversion process. The solid phase initially includes a predefined amount of moisture, volatile, and char, which are gradually depleted during the process until all solid parts convert to gaseous elements. A one-dimensional walking column strategy which models a 2D steady-state fuel bed geometry, coupled with overbed black-box modeling, is visualized in Figure 2. Heat flux is then obtained from oxidation of combustible gases in the overbed zone.  The solid and gas phase governing equations used in this work are demonstrated in Table 1.
Other information related to the modeling, such as kinetic formula of subprocesses, coefficients of specific heat, conductivity, source terms of governing equations, and so on, are presented in Table 2. Gas pressure drop for laminar flow in the packed fuel bed is also introduced in Table 2, where L is the fuel bed height, g c is the gravitational constant, ε is effective porosity, µ air is absolute viscosity, ν is air velocity, d p is effective particle diameter, and ρ is air density [29]. Table 1. Conservation equations regarding the solid and gas phase of fuel packed bed conversion. Table 2. Modeling coefficients and fuel properties.

Parameter
Value REF. Solid heat capacity. J/kg K C p,s k Y k C p,k , k = drying, pyrolysis, char burning [30] Effective solid conductivity, W/m λ e f f k Y k λ k , k = drying, pyrolysis, char burning [30] Bed void fraction φ 0.42 [30] Emissivity ω 0.85 [30] Gas pressure drop in packed bed, Pa ρν 2 dp [29] During the conversion of fuels, volumetric shrinkage of single biomass particles causes the fuel bed to decay, which is referred to as fuel bed porosity [36]. An explicit approach is used in this article, so that the particle mass consumption inside the grid is continuously observed until the all grid mass is depleted. In each iteration, each grid is filled by its upper grid mass. In this way, bed shrinkage is fairly traceable while biomass conversion proceeds [37]. a 1 , a 2 , and a 3 imply a shrinkage factor of drying, pyrolysis and char oxidation, at which values 0.1, 0.15, and 0.75 are assigned, respectively.
Outflowing gas species from solid fuel conversion are presumed to be CO, CO 2 , CH 4 , H 2 , H 2 O, and C 6 H 6 , whereby some react with oxygen mainly in the overbed zone. To quantify the fraction of volatile gases during combustion, the experimental results of another publication were employed [38]. Homogeneous gas-phase reactions [39], together with the corresponding enthalpy of reaction, are presented in Equations (2)- (5). The overbed gas-phase reactions along with the in-bed char combustion eventually form the source of heat generation in the combustion chamber. The accumulation of all the previously mentioned reactions are then accounted for in heat flux calculation.

Solution and Validation
The finite difference method (FDM) is employed to discretize the conservation equations based on the given solution algorithm in Figure 3. FDM is a well-documented method suits the simple geometries. A fully implicit scheme backward difference in time and the central difference in space (BTCS) is opted for, in order to resolve the energy equation for the solid phase (heat equation), an upwind scheme for the gas-phase mass and momentum equations (advection equation), by using a BTCS scheme. Finally, an upwind scheme for the gas-phase energy and species equations (advection-diffusion equation) using the backward difference in time, the forward difference in first-order derivative term, and central difference for the second-order derivative term. Although the selected schemes are unconditionally stable, the accuracy of the solution will be protected by choosing the proper step size. Varied time step (10 −2 to 10 −3 s) has opted for the solution along with 100 grid cells for the walking column. As most elements of the time-space matrix were zero, a sparse matrix method was used to speed up the solution procedure. The Python software was applied to code the solution algorithm. Bed temperature and mass loss profiles concerning the developed model were validated versus two separate experiments implemented on packed-bed biomass combustors [40,41]. From Figure 4, there is a satisfactory agreement between the model and experiments, unless there is a reasonable deviation, caused by inherent modeling error, model assumption, measuring error, etc. The dispersion between the temperature profile in terms of the experiments and model predictions within the drying stage can be interpreted in this way: the fuel moisture resided inside the particle pores absorbs a high proportion of the heat flowing from the overbed zone. In other words, most of the absorbed heat by particle in beginning step of conversion is used to increase moisture temperature, therefore, the particle surface temperature change is much slower. In doing so, the inside particle temperature rises, while the outer surface temperature, which the thermocouple measures, does not rise at the same rate as inside. This is deemed a limitation of the fuel bed experiment. Additionally, the predicted fuel bed conversion obtained by current simulation is in line with an earlier experiment, based on a study by Yin et al. [42], showing the contribution of drying, devolatilization, and char burnout in different lengths of a moving grate.

Results and Discussion
Before demonstrating and discussing the results, the key assumptions deployed through this study are highlighted here. For thermal analysis: 1. fuel bed is uniformly propagated in depth, and reaction rate is negligible in horizontal direction; 2. no air infiltration has been taken in calculation; 3. no heat loss from the wall occurs; 4. primary and secondary air inject to the system at 373 K. In terms of economic analysis: 1. cost of technical maintenance has accumulated in capital cost; 2. operational electricity cost of the feeder, conveyer, and air fans are equal for all scenarios; 3. no government incentive has affected in cash flow diagram. Finally, for LCA: 1. the analysis boundary has been cut in a way that no manufacturing of equipment, ash landfill, and recycling process are considered; 2. the same supply location of raw fuel is selected for all scenarios.

System Adaptation Options
In order to obtain required information for boiler adaptation when switching between fuels, combustion of three variant biomass waste fuels including biomass pellets, wood waste, and RDF were examined in a moving grate biomass boiler. Fuel composition characteristics together with corresponding stoichiometric air ratio and particle density are represented in Table 3 [18,43,44]. A biomass boiler with a bed height of 20 cm, a grate length of 2 m, and an initial overbed temperature of 1250 K were applied in this study. This overbed temperature is an initial estimation of overbed temperature and would be updated during the calculation. It would also converge to the actual value as combustion proceeds. The primary to secondary air ratio was set on an appropriate value of 40:60 according to the literature [45], along with the primary air temperature on 100°C to accelerate drying. The fuel properties are used within the developed model in order to compare the combustor's properties under different feeding fuel conditions. The temperature contour of the fuel bed conversion for each fuel is shown in Figure 5. The reaction front in moving grate boilers with cross-current flow regime propagates from bed surface downward the grate (the first phase), and then inversely from the grate toward the surface of the bed, until the fuels would be completely burned (the second phase). The first phase of reaction front for the system fed with biomass pellets finished earlier, by 11 min, compared to other fuel options by 15 and 18 min for wood waste and RDF, respectively. This can be attributed to the fact that biomass pellets contain less moisture content than the other two fuels. In the second phase, however, fuel with less fixed carbon content experiences a shorter conversion time. Furthermore, it can be concluded that the maximum temperature inside the fuel bed has a direct correlation with the fixed carbon content of the fuel. From the system adaptation view, since the complete combustion of biomass pellets, wood waste, and RDF requires 42.8 min, 35.23 min, and 34.56 min correspondingly, for a grate length of 2 m, the grate velocity must be regulated by 0.77 mm/s, 0.94 mm/s, and 0.96 mm/s, respectively. Moreover, as drying, pyrolysis, and char oxidation happen in different lengths of the grate, the primary air distribution should be updated by the wind-box mechanism under the grate, in order to speed up the conversion process inside the fuel bed. According to Yin et al. [42,46], about 40% of primary air is injected to the fuel bed from the grate opening until complete volatile decomposition, and the rest of air is distributed between middle and last zone with a ratio of 2:1. In the utilized biomass combustor, three equal wind-boxes are equipped under the grate in which the primary air is distributed for an efficient conversion. Here, Table 4 is proposed in order to adapt primary air distribution of biomass bed regarding each waste fuel.  To elaborate on the mass loss of each proposed fuel during the combustion process, dashed-line profiles in Figure 6 are referred to. From the figure, it is clear that, despite the lack of similarity between the bed decay of wood waste and RDF, the mass loss behavior is quite similar, while the biomass pellet shows a different trend. The mass-loss rate of wood waste and RDF in drying stage outperform biomass pellet because of less fuel particle density; nevertheless, as the moisture content of biomass pellet is far less than others, its devolatilization stage is triggered earlier and consequently, its mass loss dominates two other ones for a few minutes. In the final stage that char starts to burn out, again, wood waste and RDF mass loss surpass the biomass pellet, owing to less density and fixed carbon content.
Solid particle temperature profiles of bed surface regarding each biomass fuel are shown in Figure 6, in which the decelerating effect of moisture content on temperature evolution is demonstrated. In the char combustion step, the temperature becomes converged, as there is no difference between all particles. Due to the fact that biomass pellets include less moisture and more char content, the temperature rise takes place in the shorter grate length and results in higher heat generation per grate length for pellets, which will be quantified later on. Volumetric concentration in terms of mole per hour for all emitting gas from the moving fuel bed is illustrated as Figure 7. Higher water vapor of the wood waste and RDF than the pellets was expected, due to more moisture content in the fuel. The carbon monoxide outflowing rate for RDF dominates one for wood waste, which can be interpreted by higher volatile matter, whereas the carbon dioxide rate originating from fixed carbon content was inversely higher for wood waste. Among the emitting gas species from fuel bed conversion, CO 2 and H 2 O remain neutral and ones that were listed through Equations (2)-(5), which react with oxygen together with heat flux from char oxidation and form the heat source in the system. Useful output heat of the combustor can be calculated using Equation (6), where N f is mole of fuel, h s, f sensible enthalpy of fuel, h s,i sensible enthalpy of each products, h f g latent heat of water vapor, and q loss is heat leakage [47]. It is assumed that flue gases leave the stack at the temperature 400 K. According to the outcomes, the useful heat flux per kg of the wet biomass pellets was obtained at 15.282 MJ, while it lowered by 32% and 30% for wood waste and RDF, respectively. It must be noted that water vapor is not condensed within the combustor, which wastes a notable amount of heat.
Taking the conversion rate of each fuel into the account given in Table 5, annual heat generation from the proposed system was attained at 548 GJ for biomass pellets, contrasted with 48% and 35% lower in terms of wood waste and RDF. Here, the point worth making is that operating time has assumed being on full-load whole year. It was similarly concluded that for complete combustion of biomass pellets, a residence time of 42.8 min must be accounted for, and it must be updated when wood waste or RDF are deployed. Gas species concentration emitted from the bed reactions are taken to find the flame temperature of the overbed zone via a black-box method. Due to the high air infiltrations rate in moving bed biomass boilers, the flame temperature is much lower than stoichiometric condition in practice. According to the literature [48,49], excess air coefficient usually changes between 1.2-3 for such systems, and the total excess air ratio of 2 is utilized here. This unwanted excess air causes significant heat loss through the combustion chamber. Although applying a screw feeder can mitigate this, it still poses a great challenge against moving bed combustor. Lower flame temperature of pellets to other fuels can be attributed to higher CO/H 2 ratio, because the hydrogen flame temperature is almost 200 • C higher.

Economic Analysis
In this section, the proposed biomass combustor under each feeding fuel condition are economically investigated for the hot water supply applications. The biomass combustor is considered in a hot water supply system equipped with hot water tanks, piping and instruments, and water circulating pumps. The approaches such as net present value (NPV), internal rate of return (IRR), annual cost (AC), and so on, have been used within the existing literature [50,51]. In this article, in addition to NPV and AC analysis of the operating system, IRR will be also addressed, in order to make the analysis more comprehensive. The AC indicator serves up a powerful metric for contrasting the costs of various equipment, especially for those with low capital and high operational costs and vice versa.
The cost of producing energy from the proposed fuels depends on the system's capital cost (C c ), installation cost (C in ), fuel storage cost (C st ), fuel cost (C f ), electricity cost (C el ), maintenance cost (C m ), unseen cost (C un ), and salvation cost (C s ). The salvation cost is the income at the end of system lifetime associated with equipment disposal. Table 6 displays the relevant values regarding the aforementioned parameters for the 25 kW biomass boiler, which can be updated via Equation (7) for variant system size [50][51][52][53]. In this research, the capacity of the proposed biomass boiler is 30 kW, therefore, the capital, installation, and electricity cost are updated in the calculation.
where C and C 0 are the capital cost of the new system and base system, and N and N 0 are system capacity corresponding alongside factor g, ranging from 0.4-0.8 for process equipment. Hence, the life cycle cost (C LCC ) of energy generation during the system's lifespan, of 25 years will be as follows: where c and i are escalation rate 4% and interest rate 2%, respectively. In this way, the annualized cost (C a ) can be obtained from Equation (9) which is a function of capital recover factor (e), achieved from interest rate (i) and equipment lifespan (n), as pointed out via Equation (10). Eventually, the cost of heat (COH) using annual generated heat Q ann and the annualized cost is calculated by Equation (11). Additionally, the IRR, which is the interest rate at which the NPV of every cash flow during the lifespan of a project becomes zero, is counted via Equation (12) where CF means annual cash flow.
The results of the economic analysis are presented in Table 7. Levelized COH reveals that RDF has the lowest heat price by 7.63 ¢/kWh, compared to 8.23 and 9.80 ¢/kWh for pellets and wood waste correspondingly, which can be attributed to the lower fuel cost of RDF. Although the annualized cost of heat from the pellet-fueled system is about 64% more than other ones, because of the corresponding higher annual heat production, the related COH factor has a marginal difference with two other options. Taking the wood waste heat price into account, which is the highest reported value between all options in Table 7, it still outperforms the US national average residential rate of 11.88 ¢/kWh, regardless of potential government incentives supporting the energy generation from cleaner resources [54]. The LCC breakdown for each particular scenario, as demonstrated in Figure 8, reveals the sway of fuel cost over the total life cycle cost of the system followed by system capital cost invested in the first year. Here, in Figure 8, the capital and installation costs are aggregated in capital cost term, which is roughly 14% for biomass pellet and 19% and 20% for wood waste and RDF, consecutively. Finally, higher IRR of pellets to RDF, in spite of lower COH of RDF, is explained by the fact that the boiler generates higher annual heat rate under the pellet's feeding conditions. cleaning and non-technical maintenance 0.01* (C c + C in + C st ) unseen cost, C un system design and loss risks, etc. 0.01* (C c + C in + C st ) salvation cost, C s system disposal 0.05*C c

Environmental Analysis
Ordinarily, the production, collection, and exploitation of the proposed fuels are associated with a variety of combinations of environmental profits and losses. From the literature, the boundary of LCA could vary from the very beginning of the process cycle to the very end of that, corresponding to the subjects of the study [55,56]. Negative impacts of the collection and treatment of raw fuels are compensated and even surpassed during the replacement of fossil fuels with these alternatives in combustion systems, leading to emission discount. An LCA approach is applied in this section, in order to compare the environmental impacts of each alternative fuel. The LCA boundary from fuel preparation until the heat generation for proposed fuels is depicted in Figure 9. Considering the aim is to compare the impact of each fuel utilization on the ecosystem, the same heat production is assumed for all scenarios. Table 8 lists the main input data employed in the LCA. The system operating site is located in Alma, Quebec, Canada where the closest fuel supplier is in distance of 160 km [57]. The boiler electricity consumption rate was assumed to be equal for all waste fuel scenarios. The quantity of emitting pollutants from each fuel combustion was accounted via the built-in SimaPro database. SimaPro software and relevant database for North America (TRACI 2.1) was employed for pollutant indices calculations. It must be mentioned that fuels originated from wood resources are deemed neutral CO 2 , and, therefore, no carbon dioxide emission is counted during the combustion process in LCA calculation, while RDF as a mixture of biogenic and inorganic materials, so it does emit CO 2 in the combustion process.  The results of the environmental analysis are demonstrated in Table 9, based on different impacts. To acquire a certain impact indicator, different factors must be combined into one prevailing factor. For instance, the value of CH 4 is converted to the equivalent of the CO 2 amount using the relevant conversion factor, 21, and is then added to the basic CO 2 amount to represent the total equivalent CO 2 for global warming impact [58]. On the same heat generation, RDF environmentally becomes more destructive in global warming, smog, and acidification factors than others, while the biomass pellets have the worst impact on ozone depletion, eutrophication, carcinogenic, and ecotoxicity, but have the most positive effect on fossil fuel depletion. Environmentally speaking, wood waste is seen as the least harmful fuel, mainly due to avoiding high fuel preprocessing. Global warming impacts are probably the most concerning factor among all sorted items in Table 9. From the results, the RDF has a negative impact in terms of global warming three times greater than that of biomass pellets, mainly caused by emitted CO 2 in combustion. In order to elaborate on the contribution of each particular subprocess on emitting greenhouse gases, Figure 10, is displayed and only shows the top three substantial affecting contributors. By switching between alternative fuels from biomass pellets to wood waste and RDF, the key contributor to the greenhouse emission is the pelletizing process, electricity use by the boiler, and pollutant emission from RDF combustion combined.

Conclusions
In this article, a thermo-economic and environmental analysis was conducted for a moving grate biomass boiler to compare the utilization of three variant waste fuels, including biomass pellets, wood waste, and RDF in the proposed system. Fuel conversion inside the fuel bed was delicately modeled using an in-house routine numerical model. The useful heat flux from the boiler was calculated through an overbed black-box model based upon each feeding fuel, and the required time of complete fuel combustion was determined. It was concluded that biomass pellets deliver the highest system thermal efficiency so that the pellet-fueled system provides 548 GJ heat annually, opposite to 282 GJ utilizing wood waste, and 355 GJ utilizing RDF. Furthermore, the maximum solid bed temperature for biomass pellets is about 11% higher than the estimation in the RDF scenario. This is mainly because of the higher fixed-carbon content of biomass pellets. Moreover, a separate comprehensive economic analysis was practiced for the proposed biomass boiler with respect to each alternative fuel, and, again, the results revealed the privilege of using biomass pellets among all other waste fuels. Although the COH of RDF scenario was lowest by 7.63 ¢/kWh, which was 7% and 22% less than biomass pellets and wood waste respectively, the IRR of the pellet-fueled system was highest by 36%, because of its higher heating capacity compared to all others. Regardless of the scenario, the fuel cost reserved the most contribution to the total life cycle cost of heat production. In the last section of the article, a comparative LCA ratified the importance of employing fuels with no carbon footprint. The CO 2 emission from the RDF combustion was obtained far higher than other plant-based fuels, which were counted as CO 2 neutral emission.
Biomass boilers should be considered to be promising facilities for heating supply in remote areas [59]. In future studies, it is recommended to investigate the deployment of combustible waste from cargo ships for heating purposes from thermo-economic and environmental viewpoints. Carrying waste is problematic for shipping companies. On the other hand, the randomness of fuel arrangements inside the fuel bed, along with fuel constituents, and size uncertainty, are some challenges encountering biomass combustor operation and modeling in the present day. Quantifying the relevant effects on the system operation is recommended. Lastly, the LCA system boundary of this study was limited to some degree, which can be extended up to the very beginning and very end for future works.