Exploration of icy moons in the search for extra-terrestrial life is becoming a major focus in the NASA community. As such, the Exobiology Extant Life Surveyor (EELS) robot has been proposed to survey Saturn's Moon, Enceladus. EELS is a snake-like robot that will use helically grousered wheels to propel itself forward through the complex terrains of Enceladus. This moon's surface is composed of a mixture of snow and ice. Mobility research in these types of terrains is still under-explored, but must be done for the EELS robot to function. As such, this thesis will focus on the methodologies required to effectively simulate wheel interaction with cohesive media from a computational perspective. Three simulation tools will be briefly discussed: COMSOL Multiphysics, EDEM-ADAMS, and projectChrono. Next, the contact models used in projectChrono will be discussed and the methodology used to implement a custom Johnson Kendall Roberts (JKR) collision model will be explained. Finally, initial results from a cone penetrometer test in projectChrono will be shown. Qualitatively, the final simulations look correct, and further work is being done to quantitatively validate them as well as simulate more complex screw geometries.