GitHub repository for SPA impact crater ejecta analysis using iSALE3D hydrocode.
The main objective of this project is to determine the distribution of impact ejecta during the excavation stage of the impact. The model parameters for the analysis are listed in the table below. Two evaluation scenarios are analysed: 50 K/km (hotter) with 10 km/s impact velocity and 10 K/km (colder) with 15 km/s impact velocity. Data from the first 400 steps of the simulations is used for the analysis.
The filtered tracer data outputs from the full workflow is available here and can be used to generate final visualizations. Check access with iSALE3D developers for simulation data.
This notebook filters the simulation data file and generates an output csv file for all ejecta tracers that land on the moon. This csv file contains all relevant ballistic parameters for analysis / plotting. Note: for large simulation data files, the compute required by this notebook is large.
Takes in the mirror csv file and generates the ejecta distribution on mollweide projection. Tracers are filtered by ejection velocity and target and impactor material.
This notebook creates 2x2 degree equal surface area bins. The layer thickness is determined by computing the histogram of landed ejecta volume and dividing by the surface area of the bins. The plots are filtered by depth of origin and shock pressure of ejecta particles.
This notebook creates ejecta thickness as a function of depth of origin along a 45 degree radial line from the SPA basin center with inverted stratigraphy. This notebook filters tracers with azimuth between 40-50 and bins these by radial distance. The histogram is divided by surface area as before to produce the layer thickness data.
This notebook calculates the great circle distance and bearing between any two latitude longitude coordinates on a sphere using the haversine formula.