The paper describes the 1-D fluid dynamic modeling of unsteady flows with chemical specie tracking in the ducts of a four-cylinder s.i. automotive engine, to predict the composition of the exhaust gas reaching the catalyst inlet. A comprehensive simulation model, based on classical and innovative numerical techniques for the solution of the governing equations, has been developed. The non-traditional shock-capturing CE-SE (Conservation Element-Solution Element) method has been extended to deal with the propagation of chemical species. A comparison of the MacCormack method plus FCT or TVD algorithms with the CE-SE method has pointed out the superiority of the latter scheme in the propagation of contact discontinuities. A realistic composition of the exhaust products in the cylinder, evaluated by a two-zone combustion model including emission sub-models, has been imposed at the opening of the exhaust valve, considering the effect of short-circuit of air during valve overlap. The fluid dynamic code has been applied to investigate the transport of the resulting chemical species along the exhaust manifold, towards the catalyst. A detailed 1-D model of the catalytic converter has been adopted, considering the distributed and concentrated pressure losses in the system.