Sparse quasi-Newton method for Navier-Stokes solution

A sparse finite difference Newton method and a sparse quasi-Newton method have been applied to the Navier-Stokes solution. Much faster convergence to the steady state has been achieved compared to the conventional time marching method. For multidimensional applications, a block line Gauss-Seidel iterative method has been used for the solution of the resulting linear system. The methods have been demonstrated for hypersonic flow solution around a sharp cone using Osher's flux difference splitting scheme for spatial discretization.