The Hubble Space Telescope (HST) has been operational for almost 30 years and throughout that time it has been bombarded by high energy charged particles colloquially referred to as cosmic rays. In this paper, we present a comprehensive study of more than 1.2 billion cosmic rays observed with HST using a custom written python package, \texttt{HSTcosmicrays}, that is available to the astronomical community. We analyzed $75,908$ dark calibration files taken as part of routine calibration programs for five different CCD imagers with operational coverage of Solar Cycle 23 and 24. We observe the expected modulation of galactic cosmic rays by solar activity. For the three imagers with the largest non-uniformity in thickness, we independently confirm the overall structure produced by fringing analyses by analyzing cosmic ray strikes across the detector field of view. We analyze STIS/CCD observations taken as HST crosses over the South Atlantic Anomaly and find a peak cosmic ray flux of $\sim1100$ $CR/s/cm^2$. We find strong evidence for two spatially confined regions over North America and Australia that exhibit increased cosmic ray fluxes at the $5\sigma$ level.

Granular materials of different sizes are present on the surface of several atmosphere-less Solar System bodies. The phenomena related to granular materials have been studied in the framework of the discipline called Granular Physics; that has been studied experimentally in the laboratory and, in the last decades, by performing numerical simulations. The Discrete Element Method simulates the mechanical behavior of a media formed by a set of particles which interact through their contact points. The difficulty in reproducing vacuum and low-gravity environments makes numerical simulations the most promising technique in the study of granular media under these conditions. In this work, relevant processes in minor bodies of the Solar System are studied using the Discrete Element Method. Results of simulations of size segregation in low-gravity environments in the cases of the asteroids Eros and Itokawa are presented. The segregation of particles with different densities was analysed, in particular, the case of comet P/Hartley 2. The surface shaking in these different gravity environments could produce the ejection of particles from the surface at very low relative velocities. The shaking causing the above processes is due to: impacts, explosions like the release of energy by the liberation of internal stresses or the re accommodation of material. Simulations of the passage of impact-induced seismic waves through a granular medium were also performed. We present several applications of the Discrete Element Methods for the study of the physical evolution of agglomerates of rocks under low-gravity environments.