3D Unsteady Computation of Stall Inception in Axial Compressors
Rotating stall is one of the unsteady phenomena in multistage axial compressors that will damage both of performance and service life of aero engines. Stall inception is a dynamic process including appearance of pre-stall disturbance, evolvement of disturbances into stall cells, and development of stall cells. The main purpose in researching stall inception is to reveal the origins of disturbances and stall cells, investigate the effects of aerodynamic design variations on stall inception, and find the effective ways to prevent engines from turning into rotating stall or surge. Numerical simulation is an economic, reliable and rapid tool to study stall inception. As stall inception is three-dimensional and unsteady, numerical simulation should be capable of describing these aspects. In this paper, a three dimensional unsteady computational model based on the three-dimensional unsteady Euler equations and the three dimensional multi actuator-disks model has been developed. The computational domain can be divided into two kinds. One is blade-free regions, which consist of upstream duct, the axial gaps among blade rows, and downstream duct. The other one is blade rows. The flows in blade-free regions considered inviscid, unsteady, and can be resolved using three-dimensional unsteady Euler equations. The blade rows are replaced by multi actuator-disks with different total-to-static characteristics. By added the correlation functions of inlet and outlet flow angles, we can compute the flow field by combining the Euler equations and the multi actuator-disks model. A two-stage low-speed compressor in NUAA has been investigated, and the predicted results indicates that the second stage comes out stall cell first, and the full developed stall cell rotates at about 40.4% rotor speed, which coincides with the experimental data.