© 2014 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.The main purpose of this study is to analyze hypersonic flow field around the blunt bodies and to design of that bodies in order to obtain minimum pressure drag. Modeling of non-equilibrium must be done properly. In this study, non-equilibriums of thermal and chemical modes are considered. Translational and rotational energy modes are assumed that energy exchange between these modes is so fast. Vibrational and electronic energy terms are neglected. Therefore, one temperature is used to model thermo-chemical non-equilibrium. Flow field is assumed as inviscid and continuum region. Moreover, there is not diffusion. Thus, to model chemical non-equilibrium, finite rate chemistry can be used. For the thermal- nonequilibrium, enthalpy, entropy and specific heat constants are obtained from curve fitting methods. The coupled flow field equations are solved by using Newton’s methods. To solve Newton’s method, Jacobian matrices evaluation is required. In terms of convergence, Jacobian matrices are obtained by using analytical methods. At the design part, sensitivities are obtained by using adjoint design methods. The aim of design part is finding a hypersonic blunt geometry with minimum pressure drag while keeping the maximum temperature smaller than the baseline value.