© 2017, American Institute of Aeronautics and Astronautics Inc, AIAA. All rights reserved.The gradient based sensitivities required by design optimization are obtained by three different methods based on three dimensional Euler equations. Finite difference, Direct and Adjoint methods are used to compute objective sensitivities. A cell centered, upwind based finite volume method is implemented to discretize the Euler equations. The flow solution is obtained by preconditioned matrix-free Newton-GMRES algorithm. The required derivatives for Adjoint and Direct methods are obtained by analytical derivation of discrete flow equations. Resulted linear systems are solved by PARDISO solver. Design variables are selected as Non-Uniform Rational B-Spline (NURBS) curve and surface control points while Radial Basis Function (RBF) is used in volume mesh deformation.