# Soot

## Soot Equations

Soot formation, growth, and oxidation is modeled using the hybrid-method of moments (HMOM) model developed by Mueller et al. [1]. This approach combines the numerical ease of the method of moments with interpolative closure (MOMIC) with the ability to capture the bimodal nature of the soot number density function (NDF) provided by the direct quadrature method of moments (DQMOM). \(M_{x,y}\) is the moment of the soot NDF, where \(x\) is the order for volume and \(y\) for surface area; these terms are modeled according to

where \(R\) is the order of the polynomial interpolation, \(a_{r,k}\) are the interpolation coefficients, \(N_0\) is the weight of the delta function, and \(V_0\) and \(S_0\) are the volume and surface area of the nucleated spherical soot particles, respectively. The location of the delta function is fixed at coordinates \(V_0\) and \(S_0\) and assumed to be equal to the nucleated particle size. The nucleated particle volume and surface area are fixed according to \(V_0 = 2 W_C C_{\rm{dimer}} / \rho_{\rm{soot}}\) and \(S_0 = (36 \pi)^{1/3} V_0^{2/3}\), where \(W_C\) is the molar mass of carbon, \(C_{\rm{dimer}}\) is the average number of carbon atoms per dimer, and \(\rho_{\rm{soot}}\) is the density of soot (\(\rho_{\rm{soot}} = 1800 {\text{ kg/m}}^3\)). If the first-order polynomial interpolation of the moments (\(R=1\)) is used, the above equation reduces to

where \(V_L\) and \(S_L\) are the mean volume and surface area of the second mode (large particles). For first-order polynomial interpolation, four transport equations are solved: \(M_{0,0}\) (number density), \(M_{1,0}\) (volume fraction, also denoted as \(f_v\)), \(M_{0,1}\), and \(N_0\). The governing equations for the soot moments are [2]

where \(\boldsymbol{J}_{M}\) is the soot mass flux and \(\dot{M}_{x,y}\) is the soot source term. The current formulation ignores the molecular diffusion and thermophoretic effects of the soot, making the first term of the right-hand side zero.

For more details regarding the HMOM model, users are encouraged to consult the references cited.

## Soot Flags and Inputs

In the

`GNUmakefile`

, specify`USE_SOOT = TRUE`

and`NUM_SOOT_MOMENTS = N`

where`N`

is the number of moments to be used in the solution, either 3 or 6.Depending on the gas phase solver, soot solving functionality can be turned on in the input file using

`pelec.add_soot_src = 1`

or`peleLM.do_soot_solve = 1`

.The chemistry model specified with

`Chemistry_model =`

in`GNUmakefile`

must contain the PAH inception species, as well as`H2`

,`H`

,`OH`

,`H2O`

,`CO`

,`C2H2`

, and`O2`

.A PAH inception species must be provided in the input file using

`soot.incept_pah =`

. Currently, only one of three inputs are accepted:`A2`

,`A3`

, or`A4`

; which correspond to naphthalene (C10H8), phenathrene (C14H10), or pyrene (C16H10).If the inception species is named something other than

`A#`

in the chemistry model, a different name can be specified using`soot.pah_name =`

. However,`soot.incept_pah`

must be set to`A2`

,`A3`

, or`A4`

.