In this paper, by using diagonal fourth order central difference method
and TLNS equations, the numerical solution of the steady viscous supersonic
axisymmetric flow is implemented over blunt cone with shock-fitting method.
Because of using high order terms of Taylor series in discretization of derivation,
this method has high accuracy and low numerical error (dispersion error) with
respect to low order method. Boundary-closure scheme has an important role in
stability of this method. By Using coarse grid in this method, the results of
numerical solution is very close to fine grid of implicit second order (Beam &
Warming) method. High accuracy of this method is more visible than second
order method when the grid is being fined. The convergence rate of this method
is higher than second order method. This method can be converged to the
machine accuracy.
