Key points are not available for this paper at this time.
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 and providing high resolution while conserving computational resources. In previous cratering simulations, the code iSALE-2D has been widely used. However, to suppress the influence of artificial reflections from the numerical boundaries, a much larger model is usually necessary in iSALE-2D simulations (a model domain of at least 10 times of the basin diameter). This undoubtedly greatly reduces computational efficiency. To overcome this limitation, in this study, we developed a novel two-dimensional code SALEc-2D based on SALE. This new hydrocode 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. To validate the reliability of our code, we compare the results from SALEc-2D and iSALE-2D with same input parameters. The results from both codes were nearly identical when PML is disabled. 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 with 16 cores. We anticipate that our code will find widespread application in numerical simulations of impact craters in the future.
Li et al. (Fri,) studied this question.