Parallel numerical simulation of impact crater with perfect matched layers

Impact craters are the primary geomorphic features on the surfaces of celestial bodies such as the Moon, and their formation has significant implications for the evolutionary history of the celestial body. The study of the impact crater formation process relies mainly on numerical simulation methods, with two-dimensional simulations capable of reproducing general patterns of impact processes while conserving computational resources. However, to mitigate the artificial reflections of shock waves at numerical boundaries, a common approach involves expanding the computational domain, greatly reducing the efficiency of numerical simulations. In this study, we developed a novel two-dimensional code SALEc-2D that employs the perfect matched layer (PML) method to suppress artificial reflections at numerical boundaries. This method enhances computational efficiency while ensuring reliable results. Additionally, we implemented MPI parallel algorithms in the new code to further improve computational efficiency. Simulations that would take over ten hours using the conventional iSALE-2D code can now be completed in less than half an hour using our code, SALEc-2D, on a standard computer. We anticipate that our code will find widespread application in numerical simulations of impact craters in the future.