As human presence in space, both biological and technical, becomes more and more prevalent, a good understanding of how the sun-Earth system behaves is of great importance. Magnetospheric substorms are a key part of the transfer of energy from the solar wind to the Earth and they affect the dynamical nature of the Earth’s magnetosphere on a daily basis. They consist of a loading of energy from the solar wind via magnetic reconnection and then a rapid release of that energy. One of the results of this energy release is the intense brightening of the polar aurora observed at Earth. At some point the energy input reaches a peak, the system recovers back to a quiet time condition and the cycle begins again. This thesis seeks to address one of outstanding questions in regards to substorms, specifically regarding the substorm recovery phase, defined by the retreat of the magnetic reconnection site, or neutral line, in the magnetotail, which is ‘why does the substorm recovery phase start when it does?’ One result presented by Baumjohann et al. [9] is that the expansion of the dipolarization of the magnetic field in the inner magnetosphere causes the neutral line to move tailward. This thesis uses the Lyon-Fedder-Mobarry magnetohydrodynamic code to model seven substorms to find out what causes the neutral line retreat in the simulation. The result presented here is that the simulation reproduces the substorms with the loading and unloading of energy, but retreat of the neutral line is a directly-driven process that is determined by when the solar wind driver is turned off.