For turbulent spray combustion flows a coupled numerical procedure was developed, This method was discretized by using generalized curvilinear coordinates to handle complex geometries. The preconditioning and eigenvalue rescaling techniques were employed to provide efficient convergences over a wide range of subsonic Mach numbers. The accuracy was validated by simulating the laminar cavity flow. The film cooling effect of a liquid rocket engine (KSR-III) were investigated by a spray combustion analysis. The film cooling showed a negative effect on the combustion efficiency. In the combustion chamber wall, the film cooling effect was revealed to be promoted by the production of fuel rich zone.