A multi-zone model was used to predict the operating range of homogeneous charge compression ignition (HCCI) engine, the boundaries of the operating range were determined by knock (presented by ring intensity), misfire (presented by sensitivity of indicated mean effective pressure to the initial temperature). A HCCI engine fueled with Di-Methyl Ether (DME) was simulated under different initial temperature and equivalence ratios, and the operating range was well produced by the model. Furthermore, the model was applied to develop the operating range for thermal stratification in the preceding condition of initial temperature and equivalence ratios. The computations were conducted using Senkin application of the CHEMKINII kinetics rate code.