Numerical solution of transonic normal shock wave-boundary layer interaction using the Bohning-Zierep model