In the arbitrary curvilinear coordinates,intermediate and final production formations of the reverse-flow annular combustor with swirl cup were numerically investigated by the four-step reaction mechanism. The renormalization group (RNG) k-ε model,eddy-dissipation model (EDM) and discrete-ordinate method were used to simulate the turbulent viscosity,reaction rate and radiative heat transfer,respectively. The stochastic particle trajectory model was applied to simulation of the droplet trajectories,size and temperature history. Discretization equations were solved with pressure implicit split-operator (PISO) algorithm. The three-dimensional two-phase reacting flow fields and profiles of combustion production were predicted on a non-staggered gird system. Calculations are in reasonable agreement with experimental data. It is shown that the numerical methods and models adopted are reliable and can be used to provide technology support for the optimum design and reducing pollutant emission in the reverse-flow combustor.