1. Introduction
In Indonesia, there has been a shift in Coal-Fired Power Plants (PLTU) from using high and medium-grade coal to low-grade coal. This policy change was primarily motivated by the comparatively lower cost of low-rank coal, which accounts for approximately 86.59% of the country's coal reserves [
1]. The largest coal reserves in Indonesia consist of lignite and sub-bituminous coal, characterized by their low calorific value and higher sulfur content. Utilizing low-quality coal will affect all coal-fired power plant's operational reliability [
2,
3]. Unfortunately, the combustion of low-rank coal, which typically has a moisture content ranging from 40 to 50% (wet basis), presents challenges in Coal-Fired Power Plant boilers, leading to environmental issues and significant greenhouse gas emissions [
4]. On the other hand, Indonesia possesses significant biomass potential, including palm oil plantations. The palm oil industry has experienced rapid expansion, with a 7.67% growth rate over the past decade and an 11.09% annual increase in palm oil production [
5]. Unfortunately, the production of palm oil and the establishment of plantations have resulted in various environmental problems, such as greenhouse gas emissions, land conversion, and substantial agricultural waste production [
6]. Such production generated solid waste, including fruit bunch peeling, palm shells, and fiber. Fiber frequently produces electricity and steam for milling [
7]. To ensure the continuous cycling of nutrients and to mitigate soil erosion, it is usual practice to reuse peeled fruit bunches as mulch in the plantation. Compared to other types of solid waste, the mixed biomass of palm fronds (FRD) exhibits superior combustion performance, as indicated by higher Rmax and Tmax values, compared to other biomass mixtures [
8].
Furthermore, using palm frond mixtures (FRD) can mitigate slagging-fouling risks, such as reducing the ash fusion temperature, deposition material, and sticky residue within the combustion chamber [
9]. This process positively impacts the coal particle combustion, ensuring complete burning before entering the furnace exit [
10,
11]. There is a growing focus on utilizing biomass waste as an alternative to fossil fuels. Biomass co-firing is a promising technology that has the potential to substantially decrease CO
2 emissions and improve the utilization of renewable energy sources in the overall energy system [
12]. It is essential to consider the specific characteristics of each biomass when designing and selecting combustion chambers and boiler operating conditions [
13]
. Biomass can be employed in two distinct ways to produce electricity via combustion: dedicated combustion and co-firing. This method offers several advantages when we use co-firing in existing power plants, particularly coal-fueled ones. These advantages encompass enhanced combustion efficiency and reduced capital expenses [
14]. Research on biomass co-firing has progressed rapidly, mainly through experiments and computational analysis. Detailed explanations have been given about the complex modeling of combustion technologies, including the simultaneous burning of coal and biomass [
15].
The process of hydrothermal (HT) processing has been extensively researched as a means of converting biomass with a high moisture content into high calorific solid fuels [
16,
17,
18,
19]. Hydrothermal processing involves hydrolysis, condensation, decarboxylation, and dehydration reactions to convert biomass materials into carbon-rich solids under autogenous pressure at moderate temperatures (150 - 350 °C) [
19]. Numerous studies have demonstrated that because hydrothermal biomass has higher carbon content, a higher calorific value, and is hydrophobic, it is more energy dense [
20,
21]. In addition, the hydrothermal method was deemed necessary for the treatment of palm biomass due to its high conversion efficiency, the abilityto treat biomass under wet conditions, its potential to produce various value-added products, and its ability to remove some inorganic constituents including potassium from the biomass, although it requires specialised equipment and complex process control [
22,
23,
24].
Additionally, a number of studies have used computational fluid dynamics (CFD) with various biomasses to examine how coal and biomass co-firing behave against current traditional methods, such as olive waste [
25], switch grass [
26], pinewood [
27], and the Cynara cardunculus [
28]. The cofiring of various biomass and natural gas blends has also been modeled using CFD [
29]. Co-firing of coal and recovered fuel [
30]. and combustion of coal and biomass for oxy-fuel [
31]. This investigation employs Computational Fluid Dynamics (CFD) as a highly accurate measurement technique to analyze the thermal performance of co-firing combustion [
32]. Unfortunately, to the best of the authors' knowledge, no research has been discovered that discusses and elucidates the combustion characteristics of palm fronds and coal co-firing using numerical simulation on tangentially fired pulverized coal boilers. The main aim of this study is to investigate the process of burning palm fronds (FRD) and coal as a fuel in a 315 MW power plant to analyze the distribution of temperature, gas velocity, and flue gas composition of CO, CO
2, and O
2 in the combustion chamber.
2. Hydrothermal-Based Co-Firing System
Various techniques can be employed to integrate biomass with coal in combustion chambers and injection, including comilling, pregasification, and parallel co-firing. Compared to other strategies, the injection technique is the most efficient due to its lower capital costs and higher co-firing ratios [
33]. Hence, the co-firing injection was selected for this analysis study. In this specific procedure, the palm fronds (FRD) are mixed with powdered coal and inserted into the combustion chamber following individual drying and grinding methods. (
Figure 1) is a schematic diagram that demonstrates the fundamental principle of the FRD co-firing system. The raw FRD is dehydrated at temperatures ranging from 200 to 250 °C [
34,
35], typically taking several minutes to two hours.
This drying process aims to enhance the calorific value and properties of the FRD. Subsequently, the coal fuel is finely ground to the required particle size for efficient combustion within the boiler chamber. Following this, both fuel types are thoroughly mixed and injected into the combustion chamber and air, facilitated by respective heaters [
36,
37]. The superheater absorbs thermal energy from the combustion process, while the economizer generates steam specifically for the turbine. The flue gas, characterized by its comparatively low exergy level, is the fluidization gas in the bed to supply heat to the drying chamber. In this study, (
Figure 2) presents a schematic representation of the combustion chamber employed, illustrating its dimensions, meshing layout, and the burner inlet located in the corner of the boiler. The dimensions of the boiler are 63.7 m in height, 13.7 m in width, and 36.2 m in length. It is utilized in a pre-existing coal-fired power plant with a 315 MWe power generation capacity.
The configuration and area of the boiler's combustion system were established through the commercial software SOLIDWORK 2020. Afterward, the process is proceeded by creating a simulation of the established system domain. In the case of co-firing simulation, the software ANSYS FLUENT 2023 R2 is utilized to calculate velocity distribution, temperature distribution, and exhaust gas composition. The simulation includes the fundamental equations (enthalpy, mass, and momentum), heat transfer radiative, turbulence, and reactions in the gaseous phase and particles. This method is optimal for quantifying fluid flow, heat and mass transfer, chemical reactions, and interactions between solids and fluids. The research examines the co-firing FRD and coal performance using the CFD method [
32,
38,
39]. CFD modeling proves to be significantly more efficient in terms of time and cost than physical experimentation while also being safe and easily scalable. Consequently, it is frequently utilized to validate experimental approaches. The CFD analysis aims to clarify combustion processes, including combustion temperatures and resulting exhaust gas concentrations. This study employed unaltered coal as the standard for combustion. The coal and FRD flow rates are determined based on the specified FRD mass fraction. The study investigates five distinct FRD mass fractions: 0% (representing 100% coal), 5%, 15%, 25%, 35%, and 50%. Observations are made on the temperature distribution and gas concentration generated during co-firing.
2.1. Governing Equations
The Eulerian-Lagrangian approach is frequently used to simulate co-firing as a two-phase reaction flow involving solid and gas phases. This technique uses the Navier-Stokes equations to examine the gas phase, regarding the solid phase as a different entity. Newton's equations of motion are employed to determine the trajectory of each particle, while a spherical technique is employed to simulate the interactions between particles [
40]. The energy and mass transfer calculation for each particle is conducted to determine factors such as temperature and gas concentration. This technique allows for a more accurate simulation by tracking the motion of every individual particle within the system. The interactions between the gaseous phase and solid particles encompass mass, momentum, and energy, performed using the particle source method within a cell. The particle's state undergoes modification as it traverses its trajectory [
41,
42]. Mathematical calculations in this approach are primarily governed by heat transfer and fluid flow. Fluid dynamics is classified as a flow with high viscosity; each governing equation can be expressed in the following fashion:
Momentum equation in the x component
Momentum equation in the y component
Momentum equation in the z component
The fluid flow equations govern mathematical processes in quantitative physics calculations, specifically Computational Fluid Dynamics (CFD). The equations incorporate density, mass conservation, specific mass fraction, temperature, enthalpy, turbulent dissipation rate (ε), and turbulent kinetic energy (k) as essential factor fluids. It is assumed that the representation of viscous flow is achieved through fluid dynamics principles [
43,
44,
45].
2.2. Turbulence
The combustion chamber flow in this study exhibits turbulence due to fluid inertia, complex geometry, and high flow rates. Turbulence plays a crucial role in the co-firing model within the combustion chamber as it significantly influences heat and mass transfer. The k-ε turbulence model is commonly used in computational fluid dynamics (CFD) modeling to calculate turbulent combustion flow. This study utilizes the k-ε model to solve the Reynolds-averaged Navier-Stokes (RANS) equations for co-firing modeling. The k-ε turbulence model is highly regarded for its effectiveness and user-friendly nature [
33,
46], making it a popular choice in numerous industrial applications [
47]. This model's two primary equations govern the turbulent dissipation rate (ε) and turbulent kinetic energy (k). It can be articulated as follows:
Turbulent kinetic energy (
k)
Turbulent dissipation rate (
ε)
The variables
Pb, Ym, Pk, and
μt denote long-term k production, turbulent viscosity, time effect, and the impact of variations in expansion on the overall energy dissipation rate, respectively. The equation determines the values:
Gi, T, and Prt are the gravitational acceleration, temperature, and turbulent Prandtl numbers.
2.3. Radiation
Solid particle combustion, including heating, ignition, and charcoal burning, relies on heat transfer through radiation. Therefore, it is essential to comprehend the rate and amount of volatile substances emitted during the devolatilization phase concerning temperature. The discrete ordinate (DO) radiation model is used in this simulation because it accurately considers the absorption of radiation heat during burning. The absorption coefficient is calculated using the weighted sum of grey gases model (WSGGM) with a fixed value of 0.6 [
48,
49]. To address the issue of radiation heat transfer, the study utilized the P-1 approach, which relies on the expansion of radiation intensity [
46]. The radiation model can be mathematically represented in the following manner:
The variables Ap, σ, V, G, σp, and Ep The terms to be represented include absorption coefficient, incident radiation, Stefan-Boltzmann constant, particulate presence, particle scattering factor, equivalent emission, and volumes all contribute to the equivalent absorption coefficient. Furthermore, fpn, Tpn, εpn, and Apn denote the emissivity, scattering factor, projected area, temperature, and particle n.
2.4. Reaction Mechanisms Particle Phases
The co-firing of coal and HT-FRD is recognised as a unique gas-solid flow system that induces chemical reactions. The model known as Eulerian-Lagrangian incorporates hydrodynamics into the analysis [
50,
51]. Coal particles and FRD chips are modeled separately as a two-phase discrete model. The Rosin-Ramler distribution approach is adopted for the particle distribution model [
52]. Various reactions take place within the particle phase, notably during charcoal combustion. Charcoal undergoes oxidation, producing carbon monoxide (CO), which is then released into the combustion chamber as a large amount of gas. Acknowledging that charcoal produced from biomass generally exhibits greater reactivity and a higher heating rate than charcoal obtained from coal is crucial. This study utilizes a cohesive and all-encompassing response mechanism to assess the process of charcoal combustion in the presence of air. The process combustion of FRD and coal can be described as follows [
36]:
2.5 Reaction Mechanisms Gas Phase
Throughout the devolatilization process, the volatile components from T-FRD and coal are disengaged and subsequently undergo a reaction with oxygen, facilitating additional combustion. The oxidation process of these volatile components is characterized by a two-step overall reaction mechanism in which CO plays a significant role [
53,
54]. The chemical composition and heat of formation of these components are determined through proximate and ultimate analyses. The Westbrook-Dryer mechanism provides insights into the reactions of volatile materials in the gas phase [
55]:
All reactions in the gas or particle phases have their kinetic parameters listed in
Table 1. The kinetic reaction rate coefficients are determined using equation (23), which is derived from the Arrhenius equation that incorporates the activation energy Ei [
40,
56,
57]. It is presumed that the coal in question belongs to the lignite kind.
4. Boundary condition
The simulation involves categorizing boundary conditions into four types: wall, mass flow inlet, pressure outlet, and interior.
Figure 3 shows the simulation domain of co-firing coal and palm fronds in a tangentially pulverized coal boiler. The primary and secondary airflow at the inlet are subjected to mass flow inlet conditions, whereas the boiler outlet is subjected to pressure outlet conditions. The primary air is the fuel feeder gas that moves through the pipes alongside coal particles, and FRD is introduced into the furnace via burner inlets A, B, C, D, and E. However, in this simulation, burner E (standby) is not operated, so burner E is designated as a wall. Additional air gas, required for combustion, enters the furnace through burner inlets AA, AB, BC, CC, DD, DE, EF, and EFF, flowing through the outer pipe.
Each particle is viewed as a solid sphere ranging from 60 to 200 mesh (74-250μm), and the simulated power generator efficiently produces 315 MWe, achieving a fundamental combustion efficiency of 30%.
Table 3 presents the boiler operating performance data. For the combustion simulation, the required amount of air is uniformly introduced into the chamber. The ratios of primary air, secondary air, and overfire air relative ratios remain unchanged compared to the scenario of pure coal combustion. Upon entering the combustion chamber, the particles undergo a series of ongoing events, which include heating, drying, devolatilization, combustion of gases and char, creation of pollutants, and emission of heat radiation [
46,
58]. The Coupled Algorithm, a method suitable for simulating complex geometries and multiple components, is employed to calculate the combined pressure and velocity of the Navier-Stokes equations. An analysis of the gas-solid two-phase flow is conducted using the Eulerian-Lagrangian technique. The Reynolds average Navier-Stokes (RANS) equations are solved in the Eulerian domain to conduct the gas phase modelling.
The coal characteristics and HT-FRD setup information at each burner location are presented in
Table 4. Throughout the simulation, fuel is injected in varying proportions based on the HT-FRD mix ratio of 5-50% across 16 burners (A-D burner areas, while E burner remains on standby). The domain is represented in a 3D combustion chamber model with a mesh count of 737,426 unstructured tetrahedral cells. The CFD simulation tracks the spatial distribution of temperature, velocity, and concentrations of exhaust gases, including CO
2, CO, and O
2.