Abstract: A local discontinuous Petrov-Galerkin method for numerical simulation of two kinds of air pollution models is constructed. Firstly, air pollution model equations are transformed into equivalent first-order differential equations with variable substitution. Secondly, discontinuous Petrov-Galerkin method is used to solve the differential equations. The method can choose different test function and trial function space, and maintains advantages of the intermittent Petrov-Galerkin method. Compared with local discontinuous finite element method, calculation formula of the method is simpler. Numerical examples show that the method has third-order accuracy and less error than the finite volume method. The algorithm provides a practical tool for numerical simulation of air pollution models.