The Navier-Stokes equation coupled with the Darcy equation through interface conditions has attracted scientists’ attention due to its wide range of applications and significant difficulty in the nonlinearity and interface conditions. This presentation discusses a multi-physics domain decomposition method for decoupling the coupled Navier-Stokes-Darcy system with the Beavers-Joseph interface condition. The wellposedness of this system is first showed by using a branch of singular solutions and the existing theoretical results on the Beavers-Joseph interface condition. Then Robin boundary conditions on the interface are constructed based on the physical interface conditions to decouple the Navier-Stokes and Darcy parts of the system. A parallel iterative domain decomposition method is developed according to these Robin boundary conditions and then analyzed for the convergence, especially for the realistic parameters. Numerical examples are presented to illustrate the features of this method and verify the theoretical results.