A comprehensive methodology for efficient full-scale simulation of diesel exhaust aftertreatment systems is presented in this work. Specific algorithms and new numerical solvers were developed to favor easy and efficient mesh generation of DPFs, in order to reduce the computational time without affecting the quality of the results. A new multidimensional parallel FV solver for compressible and reacting flows through thin porous medium layers has been developed in the OpenFOAMR technology for the full scale simulation ofthe hydrodynamics, of the transport and deposition of the soot during the filter loading phase. The proposed techniques were implemented into the Lib-ICE code, which is a set of applications and libraries for multi dimensional engine modeling based on the OpenFOAMR technology. In addition the code has been enhanced to be applied to the simulation of catalytic devices. In particular, the developed libraries have been used to model a NOx abatement system based on the SCR technology, which is usually composed by a urea injector, an hydrolysis + SCR catalyst and an oxidation catalyst.