Flux difference splitting for three-dimensional steady incompressible navier-stokes equations in curvilinear coordinates