In this paper, a domain decomposition method for the exterior Helmholtz problem is investigated. The unboundary domain is divided into some non overlapping subdomains. The natural integral operator is used as the artificial boundary conditions on the exterior boundary of the computational domains. The convergence of the algorithm is given in the sense of energy norm. Finally, the discrete problem is discussed and some numerical examples are presented.