Numerical Simulation of Partial Combustion for Biomass Tar Elimination in Two-Stage Gasifier ()
1. Introduction
Biomass gasification is a potential way for the usage of biomass energy. Tar contained in the fuel gas is the main obstacle for the usage of biomass widely. There are quite a lot of methods to reduce tar content [1]. Considering the economy and environment-friendly, Partial oxidative/ (combustion) is a attracting way , for it’s no need for expensive mechanical equipment and can lead to a very high tar reduction. Partial combustion is kind of combustion process that ER is no more than 1. Some scholars have conducted a series of research on this area. Brandt et al. [2] did some work on a 100-kW two-stage gasifier, and the tar content after char bed could be as low as 10 - 40 mg/kg dry wood. X. Guo et al. [3] did research on air-steam gasification of biomass micron fuel (BMF) in a cyclone gasifier. The experimental results showed that the gasification performance was best with ER being 0.37 and S/B being 0.31. Two-stage downdraft gasifier is a kind of low-tar method for biomass conversion. The throat in the gasifier is a crucial part for tar elimination.
The design of throat and choice of reaction condition is crucial. Numerical calculation is a good way to solve this problem, for it’s cheap and high efficient [4]. But there are still very few researches on this area especially on the designing of the throat for partial combustion of biomass tar. This paper gives a detail simulation that will optimize the operation of gasification and the designing of gasifier. A 2D laminar flow reaction was calculated and compared to experiment results to test and verify our model. 3D model of throat combustion area was established to optimize partial oxidation for tar destruction.
2. The Physics and Mathematic Model
Partial combustion is a complex process that contains several parts: fluid flow, mass transfer, heat transfer and chemical reaction. These four parts couple with each other in a particular reactor and it’s nearly impossible to obtain its analytic solution. Numerical calculation is a considerable method to solve it. To use numerical calculation, we first need to establish the physic and mathematic model.
2.1. Flow and Heat Transfer
Fluid flow equations:
(GE-1)
(GE-2)
Heat transfer equation:
(GE-3)
Turbulence k-e equation:
(GE-4)
(GE-5)
Special transport reaction:
(GE-6)
(GE-7)
2.2. Chemical Mechanism
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
Table 1 is the kinetic parameters used in this model. All this parameters are from references.
(17)
Arrhenius reaction rate was calculated from Equation (17).
2.3. Boundary Condition and Solutions
Commercial CFD software FLUENT was used to stimulate this model. Standard model was chosen to calculate turbulent flow. SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm couples velocity and pressure. EDC model can couple turbulent flow and chemical reaction very precisely. Total mesh grid is
Table 1. Kinetic constants used in this model.
about 30,000 and mesh dependence test has been conducted. Calculation was done on a 3.0 GHz computer for 10 hours.
Verification mechanism experiments were conducted on the Laminar flow reactor in Biomass Research Center in SJTU as introduced in our formal works [15].
Figure 1 shows the boundary of the reactor. Boundary conditions: As phenol is the main pyrolysis products of wood, we choose phenol as the tar. Composition and flow rate of the inlet gases: C6H6O (0.45478), CO (0.0303), CO2 (0.14078), CH4 (0.0063), H2 (0.00037), H2O (0.36747); Inlet temperature: 500˚C, Inlet turbulence intensity: 10%, Inlet turbulence length: 0.035 m, Air inlet temperature: 27˚C, Wall temperature: 900˚C.
3. Results and Analysis
3.1. Mechanism Verification Model
This part was to test and verify our model in a Pipe flow reactor. Compared to the experiments in different conditions (ER = 0, 0.029, 0.1, 0.153, 0.2, 0.278, 0.34, 0.4), we can have a clear idea whether this model is reasonable and optimize the ER for our next research.
Figure 2 shows that the temperature increases with the increase of ER. Partial oxidation is an exothermal process. The more oxygen, the more sufficient the reaction will be. But more oxygen will decrease the calorific value of product gas. At the same time, temperature distribution may have an influence on the soot production because high temperature without oxygen may lead PAHs converted into soot.