Use of spline approximations for higher-order accurate solutions of Navier-Stokes equations in primitive variables