This thesis contributes to increasing the technology readiness level of hydrogen storage in gas reservoirs, which will be required when hydrogen has become a major energy carrier in the future Dutch energy system. It addresses the mixing processes with resident gases that occur during hydrogen storage operations in gas fields, and analyzes their implementation in a reservoir simulator. The CMG GEM reservoir simulator allows incorporation of mechanical dispersion and effective molecular diffusion into its simulations, while solving the advection dispersion transport equation fully implicitly, as well as gravitational segregation and other macro scale mixing phenomena. Mechanical dispersion is quantified by its main parameter dispersivity, for which a large uncertainty exists in the literature. Dispersivity represents pore scale fluid velocity differences caused by microscale heterogeneities in a porous medium. Due to the lack of adequate hydrogen dispersivity experiments, a range of dispersivity values is collected from literature on (groundwater) dispersivity experiments in sandstones. A sensitivity analysis is conducted, in which the influence of the different mixing processes on the mixing between working gas and cushion gas is analyzed. For this, a conceptual reservoir model (radial and homogeneous), with properties based on Dutch sandstone gas fields was built in CMG GEM. The effect of molecular diffusion on mixing proves to be negligible compared to mechanical dispersion at typical reservoir flow rates. Furthermore, the results of the simulations prove to be significantly influenced by numerical dispersion, which is a calculation error, dependent on grid size and time step. The effect of numerical dispersion compared to mechanical dispersion is quantitatively analyzed, after which various options are introduced to deal with numerical dispersion in a way that the physical processes are most realistically represented. The work in this thesis demonstrates the challenges of realistically implementing hydrogen storage mixing processes in a reservoir simulator, and should be regarded as a foundation for further research.