First, we developed a dynamical simulator of the particle motion around asteroids. Second, we implemented a novel statistical ejecta model. This novel formulation describes the ejecta field as a continuum via probability density functions. The obtained distributions are used to statistically sample the ejecta and they can be integrated analytically to estimate the number of fragments. This feature has been exploited for devising a novel sampling technique to effectively describe the entire ejecta field via “representative fragments.” High-velocity impact on small bodies can create as many as 1014 fragments. It is impossible to generate and propagate such many samples in a timely fashion. The “representative fragments” strategy considerably reduces the number of samples to simulate, allowing a faster propagation and maintaining accuracy (Figure 6, Figure 7). This work was presented at the 2022 AIAA SciTech Forum in January 2022 and published in peer-reviewed journal Icarus (
https://doi.org/10.1016/j.icarus.2023.115432(opens in new window)).
Two possible collection methods have been identified: (1) a collection of particles orbiting the asteroid for a sufficient time after the impact and (2) the collection of those particles passing in the neighbourhood of the Lagrangian point, L2. We studied the effect of the target asteroid's density, size, and equivalent strength; we also considered the impact location. The results of this work include the generation of maps to quickly evaluate a target potential based on its macroscopic properties (i.e. size and density), and the identification of a shortlist of "best targets" for the two collection strategies. The work's outcome was presented at the 72nd International Astronautical Congress (IAC) in October 2021. It is published in the peer-reviewed journal Acta Astronautica (
Figure 1, Figure 2, Figure 3, Figure 4(opens in new window).
We also studied the long-term behaviour of the ejected particles. This is important to assess the availability of samples after several days, allowing easier planning of mission operations. We characterised the number of available particles, their size, and residence time as a function of impact location (Figure 8). This analysis's results were presented at the 73rd International Astronautical Congress in September 2022.
We also performed preliminary analyses of the Hayabusa2 extended mission that still has a projectile among its payloads. We analysed the effects of the projectile impact on the asteroid surface, under the challenging dynamical environment (Figure 9), and the spacecraft's safety during the proximity operations for the projectile release. These analyses have been presented in several meetings held with the astrodynamics and science team of the Hayabusa2 SHARP mission.
A particle density and flux estimation methodology has been developed to assess the number of particles encountered by a collection device as a function of the spacecraft trajectory (Figure 10). Given a specific impact and instrument size, we design optimal collection paths for the spacecraft. Finally, we developed a control algorithm to guide the spacecraft along the optimal path (Figure 11). We investigated several scenarios, proving that in-orbit collection is feasible with the potential to collect hundreds of thousands of samples in the sub-millimetre and millimetre range. These results were presented at the 32nd JAXA Workshop on Flight Mechanics and Astrodynamics in July 2022 and the 73rd International Astronautical Congress in September 2022. They were also presented at the 2nd International Stardust Conference in November 2023.
Finally, we developed a tool for characterising the ejecta plume generated by an impactor. The importance of such an analysis is twofold: a deeper understanding of the impact phenomena for better model developments and knowledge of asteroid composition and formation, and the plume characterisation to identify regions of higher concentrations of particles. We also propose a visualisation tool to generate synthetic images of impact cratering events by exploiting the fragment's sampling, propagation and density estimate procedures (Figure 12).
