A numerical simulation model has been developed to predict the attenuation features of silencers in the frequency domain. The model is fully nonlinear, solving the conservation equations for one-dimensional, unsteady, compressible flow in duct systems. The numerical techniques adopted are the MacCormack and the Lax-Wendroff shock-capturing finite difference schemes, and the mesh method of characteristics. The numerical code is able to simulate the gasdynamic response of a muffling pipe system to a harmonic pressure perturbation, enabling the evaluation of the transfer function of the equivalent acoustic filter. The simulation has been extended to several silencers, resorting to appropriate boundary conditions and acoustically equivalent duct systems. Helmholtz resonators, expansion chambers with internal orifices, extended inlets and outlets, and other more complex silencing duct systems have been modeled, achieving a satisfactory agreement between predicted results and experimental data, measured in a semi-anechoic room. Furthermore, an investigation of perforated pipes has been carried out, following a detailed approach to model plane-wave propagation in the perforate and in the cavity liner. The method suggested takes into account every single hole with its diameter, length, and friction coefficient, allowing for the correct modeling of the reactive and dissipative behavior of the silencer. ©