An Eulerian-Lagrangian two phase large-eddy simulation (LES)is described for calculating the instantaneous gas-liquid two-phase turbulentreacting flow in an annular combustor by using different sub-grid scalecombustion models. The k-equation sub-grid scale model is used to simulatethe sub-grid eddy viscosity. The six-heat flux model is applied topredictions of the radiant heat transfer. The EBU combustion sub-grid scalemodel(E-A model), the turbulent combustion second-order moment algebraicmodel(SOM-A Model) and the second-order moment transport model (SOMmodel) are,respectively, employed to determine the chemical reaction rate.Three-dimensional body fitted grid is generated by an elliptic gridgeneration procedure. The gas phase is solved with SIMPLE algorithm andhybrid scheme in the non-staggered grid system. The liquid phase equations aresolved in a Lagrangian frame of reference by PSIC algorithm. From simulationresults, the exchange of mass, momentum and energy between gas and fueldroplet can be obtained. Experimental and predicted results are shown to bein qualitative agreement, and the comparison shows that the large-eddysimulation method in an arbitrary curvilinear coordinates may be used toinvestigate two-phase reacting flows and the combustion process by usingdifferent sub-grid scale combustion models. SOM model is better than E-AModel and SOM-A Model in the simulations.