  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.
