Study on Powder Coke Combustion and Pollution Emission Characteristics of Fluidized Bed Boilers

: The ﬂuidized reactor is widely used in a number of chemical processes due to its high gas-particle contacting e ﬃ ciency and excellent performance on solid mixing. An improved numerical framework based on the multiphase particle-in-cell (MP-PIC) method has been developed to simulate the processes of gas–solid ﬂow and chemical reactions in a ﬂuidized bed. Experiments have been carried out with a 3-MW circulating ﬂuidized bed with a height of 24.5 m and a cross section of 1 m 2 . In order to obtain the relationship between pollutant discharge and operating conditions and to better guide the operation of the power plant, a series of tests and simulations were carried out. The distributions of temperature and gas concentration along the furnace from simulations achieved good accuracy compared with experimental data, indicating that this numerical framework is suitable for solving complex gas–solid ﬂow and reactions in ﬂuidized bed reactors. Through a series of experiments, the factors a ﬀ ecting the concentration of NO x and SO x emissions during the steady-state combustion of the normal temperature of powder coke were obtained, which provided some future guidance for the operation of a power plant burning the same kind of fuel.


Introduction
The fluidized reactor is widely used in a number of chemical processes (e.g., circulating fluidized bed boiler combustion, biomass pyrolysis, catalytic reaction, and conceptual fluidized bed nuclear reactor [1][2][3][4]) due to its high gas-particle contacting efficiency, wide fuel adaptability, and excellent performance on solid mixing. Powder coke, which is also called semicoke or coke powder, was used as fuel in this fluidized bed boiler research. The powder coke used in the experiment was a high-quality Jurassic coal block produced in China's Shenfu coalfield, which was carbonized and burned in an internally heated vertical retort. Its structure is blocky and its color is light black. It is a fuel and reducing agent for ferroalloys, fertilizers, calcium carbide, blast furnace injection, and other applications. As a new type of fluidized bed boiler fuel, it has a relatively high fixed carbon content; low ash, aluminum, and sulfur contents; and its specific resistance and chemical activity are high. With the slowdown of China's domestic thermal power market growth, higher requirements for energy savings and emissions reductions for new power plants are required. Thus, vigorous development of circulating fluidized bed (CFB) boiler combustion technology for special fuels has important strategic significance for broadening the future circulating fluidized bed market.
In order to study the fuel combustion and emission characteristics of fluidized bed boilers, a large number of experiments have been implemented over the past years to optimize the operating S gi is the sum of the mass changes of each component due to the chemical reaction, and i is the number of gas-phase species. Momentum equation: ∂ ε g ρ g u g ∂t + ∇ · ε g ρ g u g u g = −ε g ∇p + ∇ · ε g τ g + ε g ρ g g + β u s − u g (2) where β is the gas-solid drag coefficient [19].
where Φ represents the viscous dissipation term, .
Q represents the energy source term, and λ g represents gas fluid thermal conductivity.
. q D represents enthalpy diffusion, written as

Solid Phase
Each component solid-phase particle was mapped to a discrete form by statistical weighting, and the m-th component particle having the same particle diameter D m and density ρ m was replaced by N m calculated particles (i.e., parcels). For the M groups of fuel particles, the number of real particles N was calculated by; the spatial distribution of the parcels over time was f X    velocity, diameter, density, and statistical weight of a parcel. For the description of parcel motion, Newtonian mechanics was used for analysis. To simplify the discussion process, the superscript i indicating the i-th parcel was omitted.
The position equation of the parcel: Equation of motion: m p du p (t) dt = m p g + F p,drag (t) + F p,coll (t) (9) where g is gravity, F p,drag (t) is the drag force of the parcel (divided into solid-phase pressure gradient field drag force and gas-solid phase drag force), and F p,coll (t) is the interaction force between the parcels. Since the calculated particles use the same particle size as the actual particles, it was necessary to consider the interaction of W p particles in the calculation of the drag force of the particles. For the i-th calculated particle of the m-th solid particle within the k-th calculation grid, the drag force was where ∇P g,k is the solid-phase pressure gradient within the k-th calculation grid, u g X p is the average velocity of the gas phase at the coordinate X p , V p = W p πD 3 p 6 is the volume of the parcels, and β m is the gas-solid phase energy transfer coefficient.
The gas-solid phase energy exchange of the m-th component in the k-th cell: where V k is the volume of the k-th cell. The frictional stress between the particles is [20,21] where ∇τ is the interparticle pressure. The continuous particle corresponding force model based on Harris and Crighton [22] was selected to calculate the ∇τ: where P s is a pressure constant (usually 1-100 Pa), ε s,cp is the close packing state volume fraction, β is a constant (2 ≤ β ≤ 5 is recommended), and ε is generally 10 −7 , which was used to eliminate the singularity of the numerical calculation [23].

Chemical Reactions
Gas-solid fluidized combustion is a process containing a series of complex chemical reactions. The main components of the fuel in the fluidized bed are fixed carbon, volatiles, moisture, and inert ash. After the fuel enters the furnace, the water evaporates rapidly, the fuel particles are heated and volatilized, and the remaining coke particles react with oxygen or water vapor to burn. Ash stays longer in the furnace than other components, so the calculation of ash can be handled separately. The combustible gases of different compositions generated during the coke combustion and volatilization analysis processes are burned in oxygen. The combustion process of powder coke is shown in Figure 1.
Gas-solid fluidized combustion is a process containing a series of complex chemical reactions. The main components of the fuel in the fluidized bed are fixed carbon, volatiles, moisture, and inert ash. After the fuel enters the furnace, the water evaporates rapidly, the fuel particles are heated and volatilized, and the remaining coke particles react with oxygen or water vapor to burn. Ash stays longer in the furnace than other components, so the calculation of ash can be handled separately. The combustible gases of different compositions generated during the coke combustion and volatilization analysis processes are burned in oxygen. The combustion process of powder coke is shown in Figure  1. In order to calculate the gas-solid phase heterogeneous reaction (such as coke combustion, limestone decomposition, desulfurization, denitration reaction, etc.), the average cell chemistry method was selected. In this method, the parameters of the discrete solid-phase particles are converted into intragrid parameters before the chemical reaction is calculated, and the consumption and generation of each component are calculated separately in each calculation grid.
There are thousands of reactions in a circulating fluidized bed boiler. It is almost impossible to calculate every reaction that occurs in the furnace. In order to simulate those reactions, such as devolatilization, volatile and char combustion, desulfurization, NOx formation, and denitration, only the major reactions have been taken into consideration. The corresponding reaction rates were given in our previous research [24]. In order to calculate the gas-solid phase heterogeneous reaction (such as coke combustion, limestone decomposition, desulfurization, denitration reaction, etc.), the average cell chemistry method was selected. In this method, the parameters of the discrete solid-phase particles are converted into intragrid parameters before the chemical reaction is calculated, and the consumption and generation of each component are calculated separately in each calculation grid.
There are thousands of reactions in a circulating fluidized bed boiler. It is almost impossible to calculate every reaction that occurs in the furnace. In order to simulate those reactions, such as devolatilization, volatile and char combustion, desulfurization, NO x formation, and denitration, only the major reactions have been taken into consideration. The corresponding reaction rates were given in our previous research [24].

Numerical Solution
In the MP-PIC method, the discretization of the solid phase is a major advantage over the two-fluid method. This statistical average and rediscrete processing balance the time consumption of the chemical reaction calculation while maintaining the accuracy of solid-phase flow simulation. In the actual operation process, when selecting the appropriate N m , one can not only maintain the same calculation accuracy as the discrete element method (DEM) but also adapt to large-scale industrial simulation applications. The calculation process is shown in Figure 2 below.
After the calculation starts, the data is first initialized. Initialization includes not only initialization of the gas-phase velocity field, pressure field, density field, and temperature field but also the calculation of particle discretization. The velocity and pressure distribution of the gas phase are calculated by the SIMPLE algorithm [25]. For the calculation of the solid phase, the force and motion are analyzed for every parcel separately. After the gas-solid flow calculation is completed, the information on the discrete particle phase is then averaged inside the grid. A series of heterogeneous chemical reactions, such as particle combustion and gasification, are calculated based on the averaging treatment, and the homogeneous reactions of the gas phase are also calculated. After the chemical reaction calculation is completed, the updated grid information is rediscretized into a discrete distribution of calculated particles, and then the next gas-solid flow calculation is completed. chemical reaction calculation while maintaining the accuracy of solid-phase flow simulation. In the actual operation process, when selecting the appropriate m N , one can not only maintain the same calculation accuracy as the discrete element method (DEM) but also adapt to large-scale industrial simulation applications. The calculation process is shown in Figure 2 below. After the calculation starts, the data is first initialized. Initialization includes not only initialization of the gas-phase velocity field, pressure field, density field, and temperature field but also the calculation of particle discretization. The velocity and pressure distribution of the gas phase are calculated by the SIMPLE algorithm [25]. For the calculation of the solid phase, the force and motion are analyzed for every parcel separately. After the gas-solid flow calculation is completed, the information on the discrete particle phase is then averaged inside the grid. A series of heterogeneous chemical reactions, such as particle combustion and gasification, are calculated based on the averaging treatment, and the homogeneous reactions of the gas phase are also calculated. After the chemical reaction calculation is completed, the updated grid information is rediscretized into a discrete distribution of calculated particles, and then the next gas-solid flow calculation is completed. Figure 3 shows the Dongfang boiler's experimental equipment. Figure 4 illustrates a schematic diagram of the experimental system, which contains a 24.5-m tall furnace, a cyclone, a recirculation bed, powder coke, and limestone feeding system. The main part of the fluidized bed system and the 3D grid of the furnace are shown in Figure 5. Detailed descriptions of the experiments can be found  Figure 3 shows the Dongfang boiler's experimental equipment. Figure 4 illustrates a schematic diagram of the experimental system, which contains a 24.5-m tall furnace, a cyclone, a recirculation bed, powder coke, and limestone feeding system. The main part of the fluidized bed system and the 3D grid of the furnace are shown in Figure 5. Detailed descriptions of the experiments can be found in Table 1. Table 2 is the elemental analysis of the powder coke used in the experiment. The particle size distribution of the powder coke is shown in Figure 6.

Simulation Object and Model Setup
To investigate the impact of the operating parameters on the pollutant discharge characteristics of powder coke, this experiment used the L 9 (3 4 ) type orthogonal test method. The main design and operating parameters closely related to the fluidized bed combustion process (bed temperature, oxygen, and grading combustion) were investigated to provide important scientific guidance for managing and optimizing a powder coke CFB boiler. size distribution of the powder coke is shown in Figure 6.
To investigate the impact of the operating parameters on the pollutant discharge characteristics of powder coke, this experiment used the L9 (3 4 ) type orthogonal test method. The main design and operating parameters closely related to the fluidized bed combustion process (bed temperature, oxygen, and grading combustion) were investigated to provide important scientific guidance for managing and optimizing a powder coke CFB boiler.   size distribution of the powder coke is shown in Figure 6.
To investigate the impact of the operating parameters on the pollutant discharge characteristics of powder coke, this experiment used the L9 (3 4 ) type orthogonal test method. The main design and operating parameters closely related to the fluidized bed combustion process (bed temperature, oxygen, and grading combustion) were investigated to provide important scientific guidance for managing and optimizing a powder coke CFB boiler.

Determination of Fuel Ignition Point by TG-DTG Method
After the fuel particles were ground and subjected to thermogravimetric analysis on the US Perkin Elemer STA6000 instrument, the temperature-weightlessness-weight loss rate curve of the sample was obtained (Figure 7). The ignition point was determined by the thermal gravitydifferential thermal gravity (TG-DTG) joint definition method. It can be seen from Figure 7 that with the increase in temperature, the TG curve below 120 °C had a rapidly declining trend, which was the rapid evaporation process of moisture. When the temperature exceeded 260 °C, the TG curve showed a second downward trend, which corresponded to the precipitation of volatiles, while some lowflammable volatiles began to ignite. The turning point of the TG curve was the ignition temperature, which was about 470 °C. As the temperature rose to 780 °C, the TG curve became a straight line, indicating that the fuel had burned out.

Determination of Fuel Ignition Point by TG-DTG Method
After the fuel particles were ground and subjected to thermogravimetric analysis on the US Perkin Elemer STA6000 instrument, the temperature-weightlessness-weight loss rate curve of the sample was obtained (Figure 7). The ignition point was determined by the thermal gravity-differential thermal gravity (TG-DTG) joint definition method. It can be seen from Figure 7 that with the increase in temperature, the TG curve below 120 • C had a rapidly declining trend, which was the rapid evaporation process of moisture. When the temperature exceeded 260 • C, the TG curve showed a second downward trend, which corresponded to the precipitation of volatiles, while some low-flammable volatiles began to ignite. The turning point of the TG curve was the ignition temperature, which was about 470 • C. As the temperature rose to 780 • C, the TG curve became a straight line, indicating that the fuel had burned out.

Determination of Fuel Ignition Point by TG-DTG Method
After the fuel particles were ground and subjected to thermogravimetric analysis on the US Perkin Elemer STA6000 instrument, the temperature-weightlessness-weight loss rate curve of the sample was obtained (Figure 7). The ignition point was determined by the thermal gravitydifferential thermal gravity (TG-DTG) joint definition method. It can be seen from Figure 7 that with the increase in temperature, the TG curve below 120 °C had a rapidly declining trend, which was the rapid evaporation process of moisture. When the temperature exceeded 260 °C, the TG curve showed a second downward trend, which corresponded to the precipitation of volatiles, while some lowflammable volatiles began to ignite. The turning point of the TG curve was the ignition temperature, which was about 470 °C. As the temperature rose to 780 °C, the TG curve became a straight line, indicating that the fuel had burned out.

Simulation Results of the First Set of Experimental Conditions
For the sake of verifying the accuracy of the 3D computational fluid dynamics (CFD) model based on the MP-PIC method, a set of stable operating conditions were selected for simulation. The main operating parameters of this working conditions are shown in Table 3. The calculated results were compared with the experimental results to verify the accuracy of the model.  10 show the main gas-solid parameters obtained by simulation. It can be seen from Figure 8 that there was a general tendency to move downward due to the injection of fuel particles and the backmixing of the large particle size from the return port at the bottom of the furnace. As the height increased, the fuel particles underwent a drying and crushing action, and the particle size gradually decreased. The movement of the particles in the 2.5-m dense phase region was complicated.
When the height reached about 5 m, particles in the core area had an overall upward trend, but the particles near the wall moved downward along the water wall, showing a typical "core-annulus" shape. The particles in the core area flowed upward under the action of drag. Some of the particles continuously migrated to the vicinity of the wall. When the height exceeded 10 m, reaching the dilute phase region of the furnace, the backmixing between the core rings was gradually weakened, and the small particles in the dilute phase region tended to move upward. Figure 9 shows the distribution of the solid-phase volume fraction along the height. It can also be seen that the region below about 5 m was a dense phase region, and the solid-phase fraction was gradually reduced at 5 m or more.

Simulation Results of the First Set of Experimental Conditions
For the sake of verifying the accuracy of the 3D computational fluid dynamics (CFD) model based on the MP-PIC method, a set of stable operating conditions were selected for simulation. The main operating parameters of this working conditions are shown in Table 3. The calculated results were compared with the experimental results to verify the accuracy of the model.  10 show the main gas-solid parameters obtained by simulation. It can be seen from Figure 8 that there was a general tendency to move downward due to the injection of fuel particles and the backmixing of the large particle size from the return port at the bottom of the furnace. As the height increased, the fuel particles underwent a drying and crushing action, and the particle size gradually decreased. The movement of the particles in the 2.5-m dense phase region was complicated.
When the height reached about 5 m, particles in the core area had an overall upward trend, but the particles near the wall moved downward along the water wall, showing a typical "core-annulus" shape. The particles in the core area flowed upward under the action of drag. Some of the particles continuously migrated to the vicinity of the wall. When the height exceeded 10 m, reaching the dilute phase region of the furnace, the backmixing between the core rings was gradually weakened, and the small particles in the dilute phase region tended to move upward. Figure 9 shows the distribution of the solid-phase volume fraction along the height. It can also be seen that the region below about 5 m was a dense phase region, and the solid-phase fraction was gradually reduced at 5 m or more.   The curves A-E in Figure 10 are distributions of the particle size distributions of the solid particles at heights of 0.3, 5, 15 and 20 m, respectively, in the furnace. As the height of the furnace increased, the proportion of small-sized particles became larger and larger due to the influence of solid particle precipitation and backmixing, while the proportion of large-sized particles became smaller and smaller. After some of the particles entered the cyclone with the flue gas, the particles that had undergone cyclone separation were returned to the furnace through the return sealing device. Large-sized particles were mostly concentrated in the dense phase region at the bottom of the furnace. Figure 11 shows the comparison of temperature distribution under stable operating conditions between experiments and simulation. The temperature distributions obtained by simulation and experiment were quite close. After the coke was rapidly separated into volatiles after entering the furnace, the remaining coke was burned in the lower layer of the furnace to release a large amount of heat, resulting in the temperature of the lower layer of the furnace being relatively high. With the increase in furnace height, the oxygen concentration decreased with the consumption in the area below the secondary tuyere, and the temperature of the upper dense phase zone decreased with the addition of low-temperature fuel and limestone particles. At about 5 m, the fresh air brought from the secondary tuyere caused an increase of oxygen concentration, and the volatile combustion increased, forming a locally higher temperature. The temperature gradually decreased due to the cooling of the furnace wall above the secondary air vent.   The curves A-E in Figure 10 are distributions of the particle size distributions of the solid particles at heights of 0.3, 5, 15 and 20 m, respectively, in the furnace. As the height of the furnace increased, the proportion of small-sized particles became larger and larger due to the influence of solid particle precipitation and backmixing, while the proportion of large-sized particles became smaller and smaller. After some of the particles entered the cyclone with the flue gas, the particles that had undergone cyclone separation were returned to the furnace through the return sealing device. Large-sized particles were mostly concentrated in the dense phase region at the bottom of the furnace. Figure 11 shows the comparison of temperature distribution under stable operating conditions between experiments and simulation. The temperature distributions obtained by simulation and experiment were quite close. After the coke was rapidly separated into volatiles after entering the furnace, the remaining coke was burned in the lower layer of the furnace to release a large amount of heat, resulting in the temperature of the lower layer of the furnace being relatively high. With the increase in furnace height, the oxygen concentration decreased with the consumption in the area below the secondary tuyere, and the temperature of the upper dense phase zone decreased with the addition of low-temperature fuel and limestone particles. At about 5 m, the fresh air brought from the secondary tuyere caused an increase of oxygen concentration, and the volatile combustion increased, forming a locally higher temperature. The temperature gradually decreased due to the cooling of the furnace wall above the secondary air vent. The curves A-E in Figure 10 are distributions of the particle size distributions of the solid particles at heights of 0.3, 5, 15 and 20 m, respectively, in the furnace. As the height of the furnace increased, the proportion of small-sized particles became larger and larger due to the influence of solid particle precipitation and backmixing, while the proportion of large-sized particles became smaller and smaller. After some of the particles entered the cyclone with the flue gas, the particles that had undergone cyclone separation were returned to the furnace through the return sealing device. Large-sized particles were mostly concentrated in the dense phase region at the bottom of the furnace. Figure 11 shows the comparison of temperature distribution under stable operating conditions between experiments and simulation. The temperature distributions obtained by simulation and experiment were quite close. After the coke was rapidly separated into volatiles after entering the furnace, the remaining coke was burned in the lower layer of the furnace to release a large amount of heat, resulting in the temperature of the lower layer of the furnace being relatively high. With the increase in furnace height, the oxygen concentration decreased with the consumption in the area below the secondary tuyere, and the temperature of the upper dense phase zone decreased with the addition of low-temperature fuel and limestone particles. At about 5 m, the fresh air brought from the secondary tuyere caused an increase of oxygen concentration, and the volatile combustion increased, forming a locally higher temperature. The temperature gradually decreased due to the cooling of the furnace wall above the secondary air vent. As can be seen from Figure 12, the concentration of oxygen contained in the primary air at the inlet of the furnace gradually decreased after the gas-solid phase heterogeneous reaction occurred with the fuel particles in the dense phase region. The oxygen brought about by the secondary air caused a transient rise in the oxygen concentration near the secondary tuyere, and then the oxygen concentration gradually decreased as the gas-phase reaction proceeded. The changing trend of CO2 concentration was just opposite to that of O2, and it had been increasing after a short decline near the secondary air outlet. From the trend of CO concentration, it can be seen that in the dense phase region at the bottom of the furnace, the rate of CO formation was much larger than the consumption rate. Due to the insufficient oxygen concentration and the high concentration of combustibles in the dense phase region, the combustion was inadequate and a large amount of CO was generated. In the dilute upper phase of the furnace, the concentration of CO was gradually reduced, mainly because the consumption of combustibles was exhausted, and the remaining O2 in the flue gas reacted with CO to form CO2. Generally, sulfur in coal can be divided into two parts: organic sulfur and pyrite sulfur, in addition to a small amount of sulfate sulfur. Since the fluidized bed combustion temperature was generally controlled below 1000 °C, the sulfur present in the fuel as sulfate was not SO2, so the amount As can be seen from Figure 12, the concentration of oxygen contained in the primary air at the inlet of the furnace gradually decreased after the gas-solid phase heterogeneous reaction occurred with the fuel particles in the dense phase region. The oxygen brought about by the secondary air caused a transient rise in the oxygen concentration near the secondary tuyere, and then the oxygen concentration gradually decreased as the gas-phase reaction proceeded. The changing trend of CO 2 concentration was just opposite to that of O 2 , and it had been increasing after a short decline near the secondary air outlet. From the trend of CO concentration, it can be seen that in the dense phase region at the bottom of the furnace, the rate of CO formation was much larger than the consumption rate. Due to the insufficient oxygen concentration and the high concentration of combustibles in the dense phase region, the combustion was inadequate and a large amount of CO was generated. In the dilute upper phase of the furnace, the concentration of CO was gradually reduced, mainly because the consumption of combustibles was exhausted, and the remaining O 2 in the flue gas reacted with CO to form CO 2 . As can be seen from Figure 12, the concentration of oxygen contained in the primary air at the inlet of the furnace gradually decreased after the gas-solid phase heterogeneous reaction occurred with the fuel particles in the dense phase region. The oxygen brought about by the secondary air caused a transient rise in the oxygen concentration near the secondary tuyere, and then the oxygen concentration gradually decreased as the gas-phase reaction proceeded. The changing trend of CO2 concentration was just opposite to that of O2, and it had been increasing after a short decline near the secondary air outlet. From the trend of CO concentration, it can be seen that in the dense phase region at the bottom of the furnace, the rate of CO formation was much larger than the consumption rate. Due to the insufficient oxygen concentration and the high concentration of combustibles in the dense phase region, the combustion was inadequate and a large amount of CO was generated. In the dilute upper phase of the furnace, the concentration of CO was gradually reduced, mainly because the consumption of combustibles was exhausted, and the remaining O2 in the flue gas reacted with CO to form CO2. Generally, sulfur in coal can be divided into two parts: organic sulfur and pyrite sulfur, in addition to a small amount of sulfate sulfur. Since the fluidized bed combustion temperature was generally controlled below 1000 °C, the sulfur present in the fuel as sulfate was not SO2, so the amount Generally, sulfur in coal can be divided into two parts: organic sulfur and pyrite sulfur, in addition to a small amount of sulfate sulfur. Since the fluidized bed combustion temperature was generally controlled below 1000 • C, the sulfur present in the fuel as sulfate was not SO 2 , so the amount of SO 2 produced by fluidized bed combustion was lower than that of a pulverized coal furnace. At the same time, the coal minerals in the fluidized bed boiler could absorb the SO 2 released by the combustion process. This characteristic of the fluidized bed combustion process is called the inherent Energies 2019, 12, 1424 13 of 18 desulfurization capacity of coal (i.e., the self-desulfurization capacity). Due to the low-temperature combustion characteristics of the fluidized bed combustion process, the addition of a desulfurizing agent (e.g., an alkaline earth metal oxide such as CaO) in the furnace absorbed SO 2 to form a stable sulfate compound, thereby realizing the possibility of efficient desulfurization. The main factors affecting the desulfurization efficiency of a CFB boiler include the Ca/S molar ratio, operating bed temperature, fuel particle size, operating oxygen, staged combustion (primary air rate and secondary air ratio), SO 2 residence time in the furnace, circulation material quantity, and load change [26]. It can be seen from Figure 12 that the concentration of SO 2 gradually increased with the increase of the height in the dense phase region, mainly because the sulfur element in the fuel in the bottom region of the furnace was gradually released, and a large amount of SO 2 was formed. Since powdered coke fuel has a molar ratio of calcium to sulfur of 5.66, the powder coke had strong self-desulfurization performance. As the height increased, a large amount of SO 2 was absorbed by the fuel's own CaO to form CaSO 4. When no additional desulfurizer was added, the SO 2 concentration at the exit of the furnace was only 45 mg/Nm 3 .
There are two main NO x generation pathways. One is that the nitrogen element in the air reacts with oxygen at a high temperature to form NO x , that is, thermal NO x . The other is that the nitrogen element in the fuel is thermally decomposed and oxidized to form NO x , that is, fuel-type NO x . The NO formed by fuel nitrogen accounts for more than 95% of the total NO x emission of the fluidized bed combustion mode [27]. It was assumed that the main component of NO x is NO, and the effect of N 2 O on NO x emissions was ignored in these simulations. When the coal was burned, part of the combustion was analyzed by volatilization (i.e., volatile nitrogen), and the other part remained on coke (i.e., coke nitrogen). Studies have shown that the main factors affecting the NO x emissions of circulating fluidized bed boilers are operating temperature, excess air coefficient, staged combustion, desulfurizer, fuel properties, cycle rate, selective noncatalytic reduction (SNCR), and ammonia injection [26]. It also can be seen from Figure 12 that after entering the furnace, the fuel particles rapidly generated a large amount of fuel-type NO x in a high-temperature environment, and then the NO x concentration gradually decreased with the reduction of CO and the dilution of the secondary air. The NO x concentration at the exit of the furnace was only about 90 mg/Nm 3 .
The highest concentration of pollutants occurred near the fuel feed point. When there was no limestone desulfurization and SCR/SNCR denitrification, the total discharge of pollutants from the furnace outlet was only about 135 mg/Nm 3 , which shows that the CFB boiler has a relatively low rate of pollutant emissions.

Factors Affecting Pollutant Emissions
Based on the above experiments and simulation results, it can be found that the flow combustion parameters' distribution in the test bench was very similar to that of the actual power plant fluidized bed. Therefore, a series of studies on the factors affecting the actual power plant emissions were carried out on this test bench. Table 4 shows the powder coke burning test conditions. Figure 13a-c shows three factors affecting the discharge of pollutants, which are operating bed temperature, operating oxygen, and primary air rate. Because the temperature, pressure, and moisture in the flue gas of different fuels and boilers are different, for the convenience of horizontal comparison, the concentrations of pollutants were converted to 6% O 2 , and dry flue gas concentration under standard conditions was used (1 atm and 0 • C). It can be seen in Figure 13a that the increase of bed temperature will significantly affect the emission of SO2. As the bed temperature increases, the self-desulfurization ability of the fuel decreases, resulting in a sharp rise in SO2 emission concentration. The relationships between SO2 It can be seen in Figure 13a that the increase of bed temperature will significantly affect the emission of SO 2 . As the bed temperature increases, the self-desulfurization ability of the fuel decreases, resulting in a sharp rise in SO 2 emission concentration. The relationships between SO 2 emission, operating oxygen, and primary air rate are shown in Figure 13b,c. The change trends of these two curves are quite similar. When the operating oxygen is 3.0% and the primary air rate is 40%, the SO 2 emission is the highest. With the increase of operating oxygen and primary air rate, SO 2 emissions are gradually reduced. After the operating oxygen exceeds 3.6% and the primary air rate exceeds 46%, the changes of the two have little effect on SO 2 emissions.
It also can be seen from Figure 13 that the factors affecting the NO x emission concentration during the normal bed temperature and stable combustion of the powdered coke fuel are from large to small in the order of operating oxygen > operating bed temperature > primary air rate. After the bed temperature exceeds 916 • C, the generation of NO x does not change substantially. Similarly, the amount of NO x generated during the increase of primary air rate from 40% to 46% is also small, but the change of NO x emission and operating oxygen quantity is linear. The NO x emission concentration value of the powdered coke fuel under normal bed temperature and stable combustion conditions increases with the increase of oxygen content and primary air rate. Considering the influence of bed temperature on SO 2 emissions, high temperatures (such as 965 • C) and low temperatures (such as 890 • C) are not conducive to reduce pollutant discharge in fluidized bed powder coke combustion; the optimal operating bed temperature is about 916 • C. After the L 9 (3 4 ) type orthogonal table test data processing, it was found that the optimal test conditions for making the NO x emission concentration the lowest are a bed temperature of 916 • C, an oxygen content of 3.0%, and a primary air rate of 40%.
The comparison between the simulation results of these nine test conditions and the experimental results shows that the grid-based average chemical reaction calculation and the MP-PIC-based gas-solid flow simulation have relatively high precision for powder coke fluidized bed combustion and pollutant discharge. However, when the primary air rate was 51%, the NO x emission in the simulation results showed a certain error compared with the experimental results, but the overall change trend was basically correct, which indicates that the calculation model still has room for improvement.

Conclusions
As a new type of fluidized bed fuel, powder coke, which is carbonized by coal, has been studied in this paper for its fluidized combustion and pollutant emission characteristics. The factors affecting the stability and pollutant emission of powder coke in the furnace combustion of a fluidized bed boiler have been studied through a series of experiments. At the same time, a 3D CFD model of the furnace based on the MP-PIC method has been developed. The particle size distribution of the particles was taken into consideration in this model. The particle packing form was used to describe the solid phase to reduce the quantity of calculations. The discrete gas-solid heterogeneous reactions were averaged in the grid cell. The simulated results were compared with the experimental results. The calculated temperature distribution in the furnace was similar to the experimental results, the maximum error was less than 25 • C, and the porosity distribution and solid-phase velocity distribution in the furnace were correct. The distribution of gases was similar to that of other fluidized bed boilers [24]. The results show that the model has fairly high precision and can be used for further research.
The NO x pollutant emission value of the normal bed temperature stable combustion test conditions ranged from 70.5 to 507.9 mg/Nm 3 (dry flue gas, 6% oxygen). The analysis of NO x emission characteristics of multiple test conditions showed that CFB was used. The NO x emission value was affected by various factors, such as operating bed temperature, oxygen volume, and primary air rate, during combustion. The self-desulfurization characteristics in the CFB furnace were mainly affected by fuel type, self-calcium ratio of calcium to sulfur, and boiler operation mode. For powder coke fuel fluidized combustion, operating oxygen (≤3.0%) and primary air rate (≤40%) were low and the average operating bed temperature was high (≥960 • C), which seriously affected the self-desulfurization capacity of the powder coke fuel and led to an increase in SO 2 emission concentration in the operation of oxygen from high to low (such as from 4.2% to 3.0% or less.) Since the SO 2 emission concentration is highly likely to rise sharply, it is therefore recommended to avoid the above operating conditions during CFB boiler operation and to design a process to improve the powder coke fuel furnace. The self-desulfurization ability in the internal combustion process enables the implementation of desulfurization agents such as limestone and an SO 2 emission value that meets the current strict SO 2 emission standards. This article mainly introduced the combustion and pollutant emission test of powder coke in the fluidized bed boiler test bench. At the same time, we developed a 3D CFD model based on the MP-PIC method for a fluidized bed boiler fluidization reaction system to investigate the internal gas-solid flow and reactions. As a special fluidized bed boiler fuel, powder coke has a different linearity of composition and pollutant emissions than traditional coal. In future research, we will try to extend this method to the application of a large-scale commercial 600-MW circulating fluidized bed. In the near future, we will combine the research results of the core-annulus cell model to simulate the Baima Demonstration Power Station, further verify the accuracy and stability of the MP-PIC method, and provide technical guidance for future plant design, operation, and operations based on the calculation results.