In this work, cyclic combustion simulations of a spark ignition engine were performed using the Large Eddy Simulation techniques. The KIVA-4 RANS code was modified to incorporate the LES capability. The flame surface density approach was implemented to model the combustion process. Ignition and flame kernel models were also developed to simulate the early stage of flame propagation. A dynamic procedure was formulated where all model coefficients were locally evaluated using the resolved and test filtered flow properties during the fully developed phase of combustion. A test filtering technique was adopted to use in wall bounded systems. The developed methodology was then applied to simulate the combustion and associated unsteady effects in a spark ignition engine. The implementation was validated using the experimental data taken from the same engine.Results show that, even with relatively coarser meshes used in this work, present LES implementation has been able to resolve the evolution of a large number of in-cylinder flow structures, which are more influential for engine performance. Predicted combustion rate and pressure rise is also in good agreement with the measurements. The limits of cyclic variations are well within the experimentally observed range. It has also been able to demonstrate the limits of cyclic fluctuations to a reasonable degree even with a fewer number of simulation cycles. A significant variation of flame propagation has also been predicted by the simulations.