Neutron stars are one of the most exciting nuclear physics laboratories in the Universe. With interior densities well above nuclear saturation density they allow us to probe conditions impossible to replicate on Earth. In addition the thermal energy of the star is negligible compared to the Fermi energy, and neutrons in the interior will be superfluid.Superfluidity affects the dynamics of the star, as now neutrons can flow relative to the ‘normal’ components of the star with little viscosity. A direct probe of such an effect is thought to come from pulsar ‘glitches’, sudden jumps in frequency observed in otherwise spinning down radio pulsars. Most theories of glitches are based on the idea that a large scale superfluid component of the star is decoupled from the spin-down of the ‘normal’ component, and its sudden re-coupling leads to a glitch.On theoretical grounds we expect this effect as a superfluid rotates by forming an array of quantised vortices, and these vortices are strongly attracted, or ‘pinned’, by ions in the neutron star crust (or superconducting flux tubes in the core). If the superfluid cannot expel ‘pinned’ vortices it cannot spin-down and builds up a lag with respect to the normal component, until hydrodynamical lift forces become strong enough to break the pinning.Despite the success of this picture in interpreting glitches, only recently has progress been made in quantitatively describing glitches with large scale hydrodynamical simulations, and statistics throughout the pulsar population with small scale quantum-mechanical simulations of vortex motion.This proposal aims to bridge the gap between these two scales by using inputs from quantum mechanical simulations to describe vortex unpinning in hydrodynamical simulations, which will include state of the art crustal physics and thermal conduction. We will thus quantitatively describe the response of the star to different kinds of glitches and obtain, for the first time, robust statistics.