The present work is focused on the detailed investigation of internal combustion engines in terms of flow and combustion analysis using an advanced numerical simulation method in the framework of computational fluid dynamics. The numerical approach used in this work is based on the Large Eddy Simulation technique implemented in STAR-CD, which has the intrinsic capability to accurately observe highly unsteady phenomena of engine flow characteristics. An extended version of the Coherent Flame Model, called 3-Zones Extended Coherent Flame Model (ECFM-3Z), which provides a universal combustion modelling approach available for all modes of technical reacting flows, is applied within the LES context. The ECFM-3Z is a model specially developed, to describe the combustion process taking place in modern internal combustion engines, in which depending on the mixture different modes of combustion exists. This study provides a comprehensive validation of the LES combined with the ECFM-3Z combustion model and reveals a detailed insight into incylinder physical phenomena.After an introduction to the fluid dynamics background in the form of a mathematical and physical description of turbulent reactive flows and their numerical implementation, the ECFM-3Z LES approach is first used for turbulent reactive and non-reactive flows of technical combustion systems of varying complexity to deal with different combustion modes. The results obtained for these flows provide a model performance assessment or predicting experimental measurements for the flow and scalar distribution.The main part of this work is focused on a detailed analysis of a single cylinder spark-ignition direct-injection DISI engine. To guarantee a proper definition of boundary conditions, the intake and exhaust ports are included in the numerical domain until the positions of pressure measurements. A multi-cycle LES was carried out on two meshes with different spatial resolution and provides an insight into the physical processes that shape the engine flow structure under motored conditions. Based on the comparison of first and second order velocity moments in experiments and LES, the accuracy of the numerical method employed was determined. Additionally the ability to reliably predict the flow field and to identify the nature of cyclic variations of IC engines was assessed.The ability to predict the combustion process within IC engines is analysed based on a comparison between LES and experiments in terms of flame presence probability and its temporal evaluation. Causes of combustion variability are presented in the form of large scale flow motion behavior in time.