“…(5), which expresses the production/destruction rate of the jth specie during the transport, are determined on the basis of the rates of reactions occurring in the DPF. The code is configured to transport 13 gas species (O 2 , N 2 , CO 2 , H 2 O, H 2 , CO, NO, C 3 H 6 , C 3 H 8 ) and the soot, which has been modeled by solid carbon C. By considering the soot as a solid specie transported with the gas, the underlying assumption to neglect any relative motion between the particles and the flow has been taken; this is consistent with the results of the CFD simulations [2], which showed that Stokes numbers of the particles are very low through the channels. In the internal nodes of the channels, the system of conservation equations is a quasi-linear hyperbolic problem, which can be solved by means of shockcapturing numerical methods: in the performed simulations, symmetric finite difference techniques with second order accuracy, such as the MacCormack predictor-corrector method and the twostep Lax-Wendroff method (with the addition of flux-limiting techniques) have been used.…”