Notwithstanding the continuous technological development of spark-ignition engines, the physics of in-cylinder phenomena has not been explicitly clarified to date. Furthermore, the advent of diverse engine technologies and in-cylinder compositions has increased the difficulty of understanding these phenomena using existing analysis methods. Under those circumstances, the modeling based on the fundamental combustion behaviors is becoming critical for analyzing combustion. Built upon the zero-dimensional flame surface density model, therefore, the models proposed herein attempt to reproduce a vast majority of combustion behaviors revealed to date, such as differential diffusion, thermo-diffusive instability, cutoff scales, distributed reaction, and flame quenching. As a result, we successfully predicted the combustion of numerous operating points of four different SI engines, including the extremely lean operating points of gasoline/air mixture up to [Formula: see text]. Furthermore, the extensive thermodynamic ranges are explored for the validation of the models with the aid of continuously variable valve timing modules; the pressure and the calculated temperature of the unburned mixture at ignition timing are varied from 4.5 to 28.6 bar and 559–794 K, respectively. In this study, we suggest an expression for a fully-developed turbulent flame speed well-fitted with the explored ranges of chemical and thermodynamic states by inspecting the influencing factors for flame wrinkling. Additionally, we analyze combustion using the developed models and identify several distinctive features of both stoichiometric and lean-burn points. Lastly, this study reveals the importance of flame thickness (or inner cutoff), which has a comparable or more severe effect on early flame propagation compared to planar laminar flame speed.