A detailed chemical kinetic mechanism, consisting of 22 species and 104 elementary reactions, has been used in conjunction with the multi-dimensional reactive flow code KIVA-3 to study autoignition of natural gas injected under compression ignition conditions. Calculations for three different blends of natural gas are performed on a three-dimensional computational grid by modeling both the injection and ignition processes. Ignition delay predictions at pressures and temperatures typical of top-dead-center conditions in compression ignition engines compare well with the measurements of Naber et al.  in a combustion bomb. Two different criteria, based on pressure rise and mass of fuel burned, are used to detect the onset of ignition. Parametric studies are conducted to show the effect of additives like ethane and hydrogen peroxide in increasing the fuel consumption rate. Calculations are also performed to study the effect of fuel injection rate and intake temperature on the location of ignition and the fuel burn rate. It is thus shown that apart from accurate predictions of ignition delay, the coupling between multi-dimensional flow and multi-step chemistry is essential to reveal detailed features of the ignition process.