CIP method and higher-order finite differernce method were introduced to simulate combuting multiphase flow. The both methods had advantages and shortcomings which were oriented in their characteristics. We revealed their characteristics in computational fluid dynamics and concluded that it was the best to combine the both methods. CCUP method was also adapted to overcome the instability in the simulation of compressible fluid flows and to be free from CFL regulation for speed of sound. Through the methods, we could simulate various combusting flow fields stably with a large time interval. We finally combined the method with large eddy simulation to simulate larger scale systems and to cooperate with experimental studies of multiphase combution.