12 Multiple Char Reaction
12 Multiple Char Reaction
12 Multiple Char Reaction
Introduction
The purpose of this tutorial is to provide guidelines and recommendations to set up and solve
multiple char reactions for coal combustion or gasification using finite-rate/eddy-dissipation
model.
This tutorial demonstrates how to do the following:
Use discrete phase model to set up and solve multiple char reactions for coal combustion.
Activate and set up the finite-rate/eddy-dissipation model for the reactions occurring
during combustion.
Solve the case using appropriate solver settings.
Postprocess the resulting data.
Include the radiation model and study its effect on reaction temperature.
Prerequisites
This tutorial is written with the assumption that you have completed Tutorial 1 from
ANSYS FLUENT 14.5 Tutorial Guide, and that you are familiar with the ANSYS FLUENT
navigation pane and menu structure. Some steps in the setup and solution procedure will
not be shown explicitly. In this tutorial, you will use turbulence and combustion models,
so you should have some experience with them. This tutorial will focus on the application
of these models in coal combustion and will not cover the mechanics of using these models.
Problem Description
The coal combustion system considered in this tutorial is a simple 10 m 1 m twodimensional duct as shown in Figure 1. Only half of the domain width is modeled because
of symmetry. The inlet of the 2D duct is split into two streams. A high-speed stream
near the center of the duct enters at 50 m/s and spans 0.125 m. The other stream enters at 15 m/s and spans 0.375 m. Both streams are air at 1500 K. Coal particles enter
the furnace near the center of the high-speed stream with a mass flow rate of 0.1 kg/s
(total flow rate in the furnace is 0.2 kg/s). The duct wall has a constant temperature of
1200 K. The Reynolds number, based on the inlet dimension and the average inlet velocity,
is approximately 100,000. Thus, the flow is turbulent.
Step 3: Models
1. Enable the Energy Equation.
Models
Energy Edit...
Viscous Edit...
Species Edit...
Step 4: Materials
Materials Create/Edit...
1. Copy the fluids from Fluent Database....
(a) Click Fluent Database... and select fluid from the Material Type drop-down list.
(b) Select carbon-monoxide (co) from the Fluent Fluid Materials list and click Copy.
(c) Similarly copy carbon-solid(c<s>) and hydrogen(h2) from the fluid database.
Note: There are two types of solid material definitions in Fluent.
SolidIt is used for conducting walls and solid bodies where only energy
equation is solved.
Fluid-SolidIn this case, solids like granular materials are defined as
fluids to facilitate solution of flow as well as energy equations.
Further, there are two different fluids available in the database, carbon(c) and
carbon-solid(c<s>). For defining granular carbon, select carbon-solid(c<s>).
(d) Close the Fluent Database Materials dialog box.
2. Ensure that piecewise-polynomial is selected from the Cp drop-down list for co2, co,
c<s>, h2, n2, o2, and h2o species.
Retain the default values in the Piecewise-Polynomial Profile dialog box.
3. Retain Cp as constant for the species coal mv volatiles.
4. Modify the properties of coal-mv-volatiles-air mixture.
(a) Select mixture from the Material Type drop-down list.
(b) Click Edit... for Mixture Species to open the Species dialog box.
C(s) + 0.5O2 CO
(1)
(2)
C(s) + H2 O H2 + CO
(3)
H2 + 0.5O2 H2 O
(4)
CO + 0.5O2 CO2
(5)
2
Particle
Surface
2
c<s>, o2
c<s>=1,
o2=0.5
default
default
1
co
co=1
3
Particle
Surface
2
c<s>,
co2
c<s>=1,
co2=1
default
default
1
co
co=2
default
default
default
default
4
Particle
Surface
2
c<s>,
h2o
c<s>=1,
h2o=1
default
default
2
h2, co
h2=1,
co=1
default
default
5
Volumetric
6
Volumetric
2
h2, o2
2
co, o2
h2=1,
o2=0.5
default
default
1
h2o
h2o=1
co=1,
o2=0.5
default
default
1
co2
co2=1
default
default
default
default
i. Retain the selection of all the reactions from the Reactions list.
ii. Click OK to close the Reaction Mechanisms dialog box.
(e) Retain the selection of incompressible-ideal-gas from the Density drop-down list.
(f) Retain the selection of mixing law from the Cp drop-down list.
(g) Click Change/Create and close the Create/Edit Materials dialog box.
(a) Enable Interaction with Continuous Phase in the Interaction group box.
(b) Enter 40 for Number of Continuous Phase Iterations per DPM Iteration.
(c) Enable Mean Values and RMS Values in the Contour Plots for DPM Variables.
(d) Enter 10000 for Max. Number of Steps in the Tracking Parameters group box.
(e) Enable Specify Length Scale.
(f) Retain default value of 0.01 m for Length Scale.
(g) Click OK to close the Discrete Phase Model dialog box.
Enabling Mean Values and RMS Values leads to addition of extra variables under discrete phase variables in the postprocessing section.
10
Value
0.001
0.03124
10
5
300
0.1
70e-6
200e-6
134e-6
4.52
vii. Enter the same values of X-Position, Y-Position, X-Velocity, Y-Velocity, and
Temperature for the Last Point.
viii. Click the Turbulent Dispersion tab.
A. Enable Discrete Random Walk Model in the Stochastic Tracking group
box.
B. Enter 50 for Number of Tries.
ix. Click OK to close the Set Injection Properties dialog box.
(b) Close the Injections dialog box.
3. Set the properties for the Combusting Particle, coal-mv.
Materials
coal-mv Create/Edit...
Parameter
Density (kg/m3)
Cp (j/kg-K)
Volatile Component Fraction (%)
Binary Diffusivity (m2/s)
Combustible Fraction (%)
Combustion Model
Value
1300
1000
28
5e-4
64
multiple-surface-reactions
11
(a) Click Change/Create and close the Create/Edit Materials dialog box.
A new tab Multiple Reactions will appear in the Set Injection Properties dialog box.
You can click the Multiple Reactions tab and check that the Species Mass Fractions
of c<s> is set to 1.
Step 6: Boundary Conditions
Boundary Conditions
1. Set the boundary conditions for velocity-inlet-2.
Boundary Conditions
velocity-inlet-2 Edit...
Parameter
Velocity Magnitude
Specification Method
Turbulence Intensity
Hydraulic Diameter
Temperature
Species Mass Fractions
Value
15 m/s
Intensity and Hydraulic Diameter
10%
0.75 m
1500 K
o2=0.23
velocity-inlet-8 Edit...
Parameter
Velocity Magnitude
Specification Method
Turbulence Intensity
Hydraulic Diameter
Temperature
Species Mass Fractions
Value
50 m/s
Intensity and Hydraulic Diameter
5%
0.25 m
1500 K
o2=0.23
wall-7 Edit...
(a) In the Thermal tab select Temperature from Thermal Conditions and enter 1200 K
for Temperature.
4. Set the boundary conditions for the pressure-outlet-6.
Boundary Conditions
pressure-outlet-6 Edit...
Parameter
Specification Method
Backflow Turbulence Intensity
Backflow Hydraulic Diameter
Backflow Total Temperature
Species Mass Fractions
12
Value
Intensity and Hydraulic Diameter
5%
1m
2000 K
o2=0.23
Step 7: Solution
1. Use coupled solver scheme.
Solution Methods
(a) Select Coupled from the Scheme drop-down list.
(b) Select PRESTO! from the Pressure drop-down list.
(c) Enable Pseudo Transient.
2. Initialize the flow field.
Solution Initialization Initialize
Hybrid Initialization is the default Initialization Method in ANSYS FLUENT 14.5. Refer
to the section 28.11 Hybrid Initialization, in the ANSYS FLUENT 14.5 Users Guide.
3. Run the calculation for 500 iterations. The residuals are as shown in Figure 3.
The solution converges in approximately 400 iterations.
Run calculation
13
Step 8: Postprocessing
1. Change the view to mirror the display across the symmetry plane.
Graphics and Animations Views...
(a) Select symmetry-5 from the Mirror Planes list and click Apply.
(b) Close the Views dialog box.
2. Change the colormap position.
Graphics and Animations Options...
(a) Select Bottom from the Colormap Alignment drop-down list.
(b) Close the Display Options dialog box.
3. Display the temperature contours.
Graphics and Animations
14
15
16
dpm-summary
The summary for the injections will be displayed in the console window.
Fate
Number
Fate
---Escaped - Zone 6
Fate
%Conv
------------- ---------- ------- ---------- ---------------Escaped - Zone 6
2.965e-05 0.000e+00 100.00
6.777e-05 0.000e+00
100.00
(*) - Multiple Surface Reactions -(*)
Fate
---Escaped - Zone 6
Species
Species Content (kg/s)
Names
Initial
Final
%Conv
------- ---------- ---------- ------c<s>
6.777e-05 6.781e-09
99.99
17
Radiation Edit...
coal-mv-volatiles-air Create/Edit...
18
Results
In this tutorial, an injection group is introduced at the inlet. The coal-mv particles travel
a distance before they start releasing volatiles. The reaction starts at this point and the
temperature increases. The radiation model lowers the peak temperature by taking the
heat away from the reaction zone. In finite-rate/eddy-dissipation coal combustion, coal
particles release volatiles that react with oxygen and produce combustion products. The
stoichiometric coefficients can be calculated once chemical composition of coal volatiles is
known. For information on determining coal volatile composition, see tutorial, Using the
Non-Premixed Combustion Model in the ANSYS FLUENT 14.5 Tutorial Guide.
Summary
Application of multiple char reactions and finite-rate/eddy-dissipation model in a coal combustion case has been demonstrated.
Further Improvements
This tutorial demonstrates a second order solution. You may be able to obtain a more accurate solution by adapting the mesh. Mesh adaption ensures that the solution is independent
of the mesh. In more realistic/complex cases, you can obtain non-reacting solution, reacting flow solution, and then solution with radiation similar to the tutorial, Coal Combustion
with Eddy Break Up (EBU) Model.
19