Soot Model
Soot: Theoretical Foundations
Overview
Professor, why is soot important in combustion CFD?
Soot is carbonaceous fine particles (particle size 10-100 nm) generated by incomplete combustion, and it is important for three reasons. (1) It is subject to exhaust gas regulations (PM: Particulate Matter), (2) It contributes significantly to radiative heat transfer (soot dominates flame radiation), (3) Health hazards (carcinogenicity). Soot prediction in diesel engines, aircraft engines, and industrial furnaces is an essential task.
Please explain the soot formation mechanism.
Soot formation progresses through a 4-stage process.
1. Nucleation: The initial soot nuclei are formed by the polymerization of PAHs (Polycyclic Aromatic Hydrocarbons). C2H2 (acetylene) is the main precursor for PAH growth.
2. Surface Growth: Carbon deposits on the soot particle surface via the HACA mechanism (H-Abstraction-C2H2-Addition).
3. Coagulation: Particles collide and coalesce, increasing in size.
4. Oxidation: Soot is combusted and eliminated by O2 or OH.
Moss-Brooke Model
Please explain the governing equations for the soot model used in CFD.
The Moss-Brooke 2-variable model solves two transport equations: for soot mass fraction $Y_s$ and soot number density $N$ (number of particles/kg).
The soot volume fraction $f_v$ is obtained by the following equation.
Here, $\rho_s \approx 1800$ kg/m$^3$ is the density of soot, and $d_p$ is the average particle diameter.
Key Parameters for Soot Formation
Under what conditions is more soot produced?
Let's summarize the main factors for soot formation.
| Factor | Direction for Increased Soot | Reason |
|---|---|---|
| Equivalence Ratio | Rich ($\phi > 1$) | Insufficient oxygen leads to incomplete combustion |
| Temperature | 1500-1800 K | Optimal temperature range for nucleation |
| Pressure | High Pressure | Increased collision frequency |
| Fuel Structure | Aromatic > Straight-chain | Easier generation of PAH precursors |
| Residence Time | Long | Ensures time for soot growth |
So 1500-1800 K is the temperature window for soot formation.
Correct. Below this temperature, nucleation rate is slow; above it, OH oxidation becomes dominant and soot burns out. This "soot formation window" is visualized as a $\phi$-T map. The $\phi$-T map for diesel combustion (Dec diagram) forms the basis for simultaneous reduction strategies of soot and NOx.
So the soot model is a fusion of chemical kinetics and particle dynamics.
Yes. Because it requires accurate description of both gas-phase PAH chemistry and soot particle dynamics, it is one of the most challenging fields in combustion modeling.
Soot is a Nano-sized Carbon Particle โ The Challenge of Describing the Formation of Objects Below 1nm in Diameter with Equations
Soot formation is a process where carbon particles with diameters of 1โ100 nm form, grow, and coagulate in just a few milliseconds during combustion. To represent this in CFD, the four processes โ "nucleation โ surface growth โ coagulation โ oxidation" โ must be formulated mathematically. There are different schools of thought, such as the Fenimore-Jones two-equation model and Frenklach's detailed soot model, but describing the particle size distribution (PSD) is particularly difficult. Particles of 1 nm and 100 nm behave completely differently, yet CFD grids are on the order of millimeters โ meaning we are dealing with a problem of different dimensions where the object size is one-millionth of the computational cell.
Computational Methods for Soot
Details of Numerical Methods
Please explain the numerical implementation of the soot model.
Soot models used in CFD are broadly classified into three categories.
| Model | Accuracy | Computational Cost | Features |
|---|---|---|---|
| Empirical 2-Variable Model | Low-Medium | Low | Moss-Brooke, transports $Y_s$ and $N$ |
| Method of Moments (MoM) | Medium-High | Medium | Transports moments of particle size distribution |
| Sectional Method | High | High | Resolves particle size distribution with discrete sections |
Empirical 2-Variable Model
Please start with the simplest model.
The Moss-Brooke model (standard in Fluent) describes the particle population with two variables: $Y_s$ (soot mass fraction) and $N$ (number density). It uses Arrhenius-type rate expressions for each process: nucleation, surface growth, coagulation, and oxidation. It's simple, but particle size distribution information is limited to average values.
Method of Moments (MOMIC)
What is the Moment Method?
It is a method that solves transport equations for the moments $M_r = \int_0^\infty v^r n(v) dv$ of the particle size distribution function $n(v,t)$ ($v$ is particle volume). $M_0$ corresponds to number density, $M_1$ to volume fraction. MOMIC (Method of Moments with Interpolative Closure) proposed by Frenklach & Harris is available in Fluent 2020 and later.
Sectional Method
What are the advantages of the Sectional Method?
The particle size range is divided into discrete sections (bins), and the number density in each section is transported individually. It can represent arbitrary shapes of particle size distributions with the highest accuracy, but requires transport of 20-30 additional scalar variables, resulting in high computational cost. Available in STAR-CCM+ and CONVERGE.
Fluent Setup
Please explain the soot model setup procedure in Fluent.
1. Models > Species > Species Transport (assuming combustion model is already set)
2. Models > Soot > Moss-Brooke (simple) or MOMIC (recommended)
3. Setup of PAH precursor species (C2H2, C6H6, etc.) โ must be included in the reaction mechanism
4. Coupling with radiation model โ pass soot absorption coefficient to the radiation model
An important note: the soot model requires a reaction mechanism that includes PAH precursors (at least C2H2). Global one-step mechanisms do not include C2H2, so soot calculation is not possible. Use mechanisms like DRM-19 or higher.
Coupling of Soot and Radiation
How are soot and radiation coupled?
Soot particles emit and absorb radiation across a continuous spectrum. The absorption coefficient of soot is approximated by the following equation.
Here, $C_0$ and $C_2$ are optical constants. The continuous radiation from soot is added to gas radiation (band radiation from CO2, H2O), so in flames with high soot content, radiative loss increases significantly.
So the soot model is a triple coupling of combustion model + particle model + radiation model.
Yes. Balancing model complexity and computational cost is crucial. A practical approach is to first grasp trends with the Moss-Brooke 2-variable model, then move to MOMIC or the Sectional method as needed.
Numerical Methods for Soot Models โ Choosing Between Transport Equation Methods and Sectional Methods
For soot analysis in CFD, there are mainly "2-equation models (Moss-Brookes: soot number density N + soot volume fraction f)" and "Sectional methods (tracking particle size distribution with discrete bins)." The 2-equation model has low computational cost and is standard in Fluent, but cannot provide detailed particle size distribution. The sectional method can resolve particle nucleation, coagulation, and surface growth by particle size, offering high prediction accuracy for light scattering/absorption properties. Recently, MoM (Method of Moments) is gaining attention as an intermediate approach, offering a good balance between computational cost and modeling accuracy. Practical implementation of sectional methods is accelerating for compliance with NVPM (Non-Volatile Particulate Matter) regulations for aircraft engines.
Experience the theory firsthand with the interactive simulator for this field
All Simulators