Numerical simulations of galaxy formation provide a powerful technique for calculating the non-linear evolution of cosmic structure formation. In fact, they have played an instrumental role in establishing the current standard cosmological model known as LCDM. However, unlocking the predictive power of current petaflop and future exaflop computing platforms requires a targeted effort in developing ...