This paper presents a new method to analyse the unsteady gas ow in both intake and exhaust systems of internal combustion engines. Such a method is based on the simultaneous use of a one-dimensional model applied to describe the phenomena in ducts, together with a lumped parameter scheme to investigate the cylinder or other volume behaviour, coupled with a three-dimensional model, able to guarantee detailed information on ow behaviour in complex geometry, retaining the advantages of all methods, accuracy as well as fast processing and high ow pattern resolution. The description of the one-dimensional model developed with an example of its application is presented. The integrated approach with the coupling procedure is then described. Finally the results of a multicylinder exhaust system simulation are illustrated.