Based on control theory,adjoint system for the general problem of turbomachinery aerodynamic optimization was studied and developed in the present paper by using the variation technique in the grid node coordinates combined with Jacobian Matrics of flow fluxes.Then the adjoint system for aerodynamic design optimization of turbine cascade governed by compressible Navier-Stokes equations was derived in detail.With the purpose of saving computation resources,the mathematic method presented in this paper avoids the coordinate system transforming in the traditional derivation process of the adjoint system and makes the adjoint system much more sententious.Given the general expression of objective functions consisting of both boundary integral and field integral,the adjoint equations and their boundary conditions were derived,and the final expression of the objective function gradient including only boundary integrals was formulated to reduce the CPU cost,especially for the complex 3D configurations.The adjoint system was solved numerically by using the finite volume method with an explicit 5-step Runge-Kutta scheme and Riemann approximate solution of Roe's scheme combined with multi-grid technique and local time step to accelerate the convergence procedure.Finally,based on the aerodynamic optimization theory in the present work,2D and 3D inviscid and viscous inverse design programs of axial turbomachinery cascade for both pressure distribution and isentropic Mach number distribution on the blade wall were developed,and several design optimization cases were performed successfully to demonstrate the ability and economy of the present optimization system.