This thesis outlines an iterative approach for determining the scattered wave for two dimensional multiple acoustic scattering problems using high order local absorbing boundary conditions and second order finite difference. We seek to approximate the total wave as it is scattered off of multiple arbitrarily shaped obstacles. This is done by decomposing the scattered wave into the superposition of single scattered waves. We then repeatedly solve the single scattering system for each obstacle, while updating the boundary conditions based off the incident wave and the scattered wave off the other obstacles. We solve each single scattering by enclosing the obstacle in a circular artificial boundary and generating a curvilinear coordinate system for the computational region between the obstacle and the artificial boundary. We impose an absorbing boundary condition, specifically Karp's Farfield Expansion ABC, on the artificial boundary. We use a finite difference method to discretize the governing equations and to discretize the absorbing boundary conditions. This will create a linear system whose solution will approximate the single scattered wave. The forcing vector of the linear system is determined from the total influence on the obstacle boundary from the incident wave and the scattered waves from the other obstacles. In each iteration, we solve the singular acoustic scattering problem for each obstacle by using the scattered wave approximations from the other obstacles obtained from the previous iteration. The iterations continue until the solutions converge. This iterative method scales well to multiple scattering configurations with many obstacles, and achieves errors on the order of 1E-5 in less than five minutes. This is due to using LU factorization to solve the linear systems, paired with parallelization. I will include numerical results which demonstrate the accuracy and advantages of this iterative technique.