An efficient non-linear harmonic methodology has been developed for predicting unsteady blade row interaction effects in multistage axial flow compressors. Flow variables are decomposed into time averaged variables and unsteady perturbations, resulting in time averaged equations with deterministic stress terms depending on the unsteady perturbation. The non-linear interaction between the time averaged flow field and the unsteady perturbations are included by a simultaneous pseudotime integration approach, leading to a strongly coupled solution. The stator/rotor interface treatment follows a flux averaged characteristic based mixing plane approach and includes the deterministic stress terms due to upstream running potential disturbances and downstream running wakes, resulting in the continuous nature of all parameters across the interface. The basic computational methodology is applied to the three-dimensional Navier-Stokes equations and validated against several cases. Results show that this method is much more efficient than the non-linear time-marching methods while still modeling the nonlinear unsteady blade row interaction effects.Copyright © 2000 by ASME