The problem of modeling acoustic waves scattered by an object withNeumann boundary condition is considered. The boundary condition istaken into account by means of the fictitious domain method, yieldinga first order in time mixed variational formulation for theproblem. The resulting system is discretized with two families of mixed finite elements that are compatible withmass lumping. We present numerical results illustrating that the Neumann boundary condition on the object is not always correctly taken into account when the first family of mixed finite elements is used. We, therefore, introduce the second family of mixed finite elements for which atheoretical convergence analysis is presented and error estimates areobtained. A numerical study of the convergence is also considered fora particular object geometry which shows that our theoreticalerror estimates are optimal.