A moving-grid, shock-fitting, finite element method has been implemented that can achieve high-order accuracy for flow simulations with shocks. In this approach, element edges in the computational mesh are fitted to the shock front and moved with the shock throughout the simulation. The Euler or Navier-Stokes equations are solved on the moving mesh in an arbitrary Lagrangian-Eulerian framework. The method is implemented in two-dimensions in the context of a streamwise upwind Petrov-Galerkin finite element discretization with unstructured triangular meshes and mesh adaptation. It is shown that the shock interface motion equation has a wave nature, and disturbances can propagate along the shock interface. A SUPG stabilization term is introduced to the interface motion equation that is critical for ensuring that interface disturbances do not lead to non-convergent solution behavior. The formal order of accuracy of the scheme is verified, and the performance of the proposed scheme is assessed for both inviscid and viscous problems. It was found that the present scheme predicts smooth and noise-free surface heating for hypersonic flow over a cylinder with purely irregular triangular elements.