Volterra series is one of the powerful system identification methods for representing the nonlinear dynamic system behavior. The methods of step response and impulse response are commonly applied to a discrete aerodynamic Computational Fluid Dynamic (CFD) to identify the first- and second-order Volterra kernels. A critical problem, however, is the difficulty of identifying the second-order Volterra kernels correctly in CFD-based method. In this paper the second-order Volterra kernel function is expanded in terms of Chebyshev functions to reduce the size of the problem and the accuracy of the identification is also improved based on a third-order reduced model of Volterra series