Numerical modelling and fast-solver calculation of approximately normal shocks