An Experimental and Simulation Study of Early Flame Development in a Homogeneouscharge SparkIgnition Engine
^{1}
Department of Mechanical & Nuclear Engineering, The Pennsylvania State University,
University Park,
PA  USA
^{2}
Department of Engineering “Enzo Ferrari”, University of Modena and Reggio Emilia,
Modena  Italy
^{3}
Department of Mechanical Engineering, University of Michigan,
Ann Arbor,
MI  USA
^{*} Corresponding author email: dch12@psu.edu
Received:
15
November
2016
Accepted:
7
August
2017
An integrated experimental and LargeEddy Simulation (LES) study is presented for homogeneous premixed combustion in a sparkignition engine. The engine is a singlecylinder twovalve optical research engine with transparent liner and piston: the Transparent Combustion Chamber (TCC) engine. This is a relatively simple, open engine configuration that can be used for LES model development and validation by other research groups. Pressurebased combustion analysis, optical diagnostics and LES have been combined to generate new physical insight into the early stages of combustion. The emphasis has been on developing strategies for making quantitative comparisons between highspeed/highresolution optical diagnostics and LES using common metrics for both the experiments and the simulations, and focusing on the important early flame development period. Results from two different LES turbulent combustion models are presented, using the same numerical methods and computational mesh. Both models yield CycletoCycle Variations (CCV) in combustion that are higher than what is observed in the experiments. The results reveal strengths and limitations of the experimental diagnostics and the LES models, and suggest directions for future diagnostic and simulation efforts. In particular, it has been observed that flame development between the times corresponding to the laminartoturbulent transition and 1% massburned fraction are especially important in establishing the subsequent combustion event for each cycle. This suggests a range of temporal and spatial scales over which future experimental and simulation efforts should focus.
© Y. Shekhawat et al., published by IFP Energies nouvelles, 2017
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Introduction
CycletoCycle Variations (CCV) of flow and combustion in sparkignition engines have been the subject of extensive research over the last few decades [1]. Highspeed/highresolution optical diagnostics and LargeEddy Simulation (LES) have been brought to bear to understand the root causes of CCV, and significant progress has been made in that regard. Multiple sources of variability have been identified. These include local variations in the flow, temperature and mixture composition, and in the spark discharge, as well as global variations in equivalence ratio, dilution and trapped mass [2–5]. While these contributing factors have been recognized for some time, their interactions and influences on combustion remained unclear due to a lack of multidiagnostic data and the difficulty of performing accurate LES simulations over a sufficiently large number of cycles. Lacour and Pera [5] and Baum et al. [6] performed multidiagnostic experiments that allow a deeper understanding of the coupled physics in engines. And recently these interactions have been investigated using multiparameter experimental [6,7] and simulation [6,8] approaches.
Earlier experimental and simulation studies (e.g., [9,10]) have shown that what happens during the earliest stages of combustion can determine, to a large extent, the subsequent combustion process for that cycle. Under realistic engine operating conditions, the velocity magnitude [7,8] and velocity gradient parameters [7] in the vicinity of the spark plug at the time of ignition have been found to correlate with the Indicated Mean Effective Pressure (IMEP) for the cycle. Motivated by the difficulties of making accurate pressurebased combustion measurements during the early flame development period, some LES studies have focused on earlyflamekernel development, toward developing an understanding of the governing factors that result in cyclic variation of flame growth. Granet et al. [11] demonstrated that initial flame convection within the spark plug gap is a probable reason for incomplete combustion cycles due to local quenching, while Goryntsev et al. [12,13] emphasized the superposition of flow variation and mixture quality for directinjection engines, whose combined effect on combustion CCV was elucidated by LES.
For these reasons, the focus here is on the early burn: ignition through fully developed turbulent flame. A specific goal is to assess the predictive capability of two different ignition and turbulent flame propagation models, using the same CFD code and computational mesh. The approach is to compare the simulated and measured combustion using a single set of metrics that are applicable to both the simulated and measured data, focusing on the early burn. The metrics are derived from incylinder pressure and highspeed optical measurements, and are applied to perform equivalent analysis of the simulation data. The study was conducted in a research engine that exhibits high cycletocycle variability. As will be shown later, fast and slow burning cycles are established during the Early Flame Development (EFD), defined here as the period between the Start Of Ignition (SOIgn) and a fully developed turbulent flame. Variation of the EFD duration (either between the experiments and simulations, or between individual cycles in either the experiments or the simulations) changes the combustion phasing, and thus is effectively the same as variability in the ignition timing. Thus, the turbulent flame propagation after the EFD encounters earlier or later thermodynamic states and turbulence properties from one cycle to another. In turn, the resulting combustionphase variability affects the production of work (IMEP) and the time for completion of combustion (combustion efficiency), as revealed by experimental pressurebased combustion analysis [14,15]. Thus, improper simulation of EFD will result in improper pressure or mass burned fraction evolution, even if the turbulent flame propagation model is perfect. This motivates the focus on EFD for this study.
Historically, the EFD period has been characterized using the MassBurned Fraction (MBF) of the combustion progress, computed from an Apparent Heat Release (AHR) analysis based on the measured incylinder pressure trace. Typically, a crank angle corresponding to between one and ten percent MBF (CA01–CA10) is the earliest value that is sufficiently resolved experimentally, and that has been used as the metric to characterize the EFD variability. Although values as low as CA01 can be calculated from the pressure data, singlecycle analysis has high uncertainty due to recording noise and the large dynamic range requirements for the pressure transducer and recording system (here 40–2000 kPa). As a second EFD metric, the transition between slow and rapid flame growth has been measured optically. Beretta et al. [16] observed the earlyflamegrowth period between the start of ignition and the transition to a turbulent flame, which they defined in terms of the transition from a laminar to a turbulent characteristic time scale. This concept was introduced for modeling ignition in RANS simulations by Abraham et al. [17]. The duration of this laminartoturbulent transition period, τ_{lamturb}, is measured here from the twodimensional area projection of the threedimensional flame volume, as computed from both the optical flame measurements and the LES. The laminartoturbulent transition is quantified as the timescale transition revealed in the log(A_{Flame}) versus log(t) plots (where A_{Flame} is the timevarying flame area), as suggested by Arpaci et al. [18]. It is demonstrated that τ_{lamturb} occurs at a MBF much less than one percent in the experiments.
An emphasis in this paper is to apply equivalent metrics to compare results from the experimental measurements with results from two different LES ignition and turbulent combustion models, focusing on the EFD period. As will be shown later, both models were found to produce higher variability than the experiment. Both the measured and simulated cycles show that the combustion phasing of early and lateburning cycles is established by CA10. The projected flamearea growth rate between the start of the spark discharge (SOIgn) and CA10 is compared to τ_{lamturb}, to CA01, and to the sparkplasma duration to estimate the time and length scales of the flame kernel during these periods. To assess the causes of the reasons for the disparity between the experiments and simulations, both the thermodynamic and the flow properties at SOIgn are compared.
Contributions of this paper are: (1) the first comparisons between combustion LES and experiment for a configuration that is fully accessible to other research groups; (2) comparisons of results from two different LES turbulent combustion models using the same mesh and numerics; and (3) combining pressurebased combustion analysis, optical diagnostics and LES to generate new physical insight into the early stages of combustion and the limitations of current modeling approaches. Specific questions that are addressed include:

Are there specific local or global conditions that result in a fast cycle versus a slow cycle?

Can useful criteria/metrics be established for the transition between ignition and fully developed turbulent flame propagation?

What controls the heat release rate during the early burn?

How well do current LES models represent this phase of combustion?

What experiments/models are needed to make further progress in understanding and modeling the early burn?
1 Engine Experiments
The University of Michigan transparent combustion chamber (TCCIII) optical engine is illustrated in Figure 1 [19–21]. This is a singlecylinder, sparkignited, twovalve, fourstrokecycle, pancakeshapedcombustionchamber engine with a transparent liner and piston. The geometric compression ratio is 10:1, and bore × stroke are 92 × 86 mm. This is a research engine that has been developed specifically to facilitate optical access for ParticleImage Velocimetry (PIV) and other highfidelity laserbased diagnostics, as well as simple mesh generation for CFD. This geometrically simple engine configuration is known to exhibit highly variable incylinder flow. The combination of a pancakeshaped combustion chamber, large pistontovalve diameter ratio (92 mm/30 mm), and undirected ports produces large CCV of the largescale, motored flow structures [22,23]. Details of the geometric configuration and operating conditions, and extensive motored engine PIV data, are available at [21]. Several groups have been using the motored TCC engine data to test and validate physical models and numerical algorithms for incylinder LES [19,20,23–25].
The measured combustion data are from the fired TCCIII engine experiments conducted by Schiffmann [7]. There, the measured data for a range of operating conditions (lean, rich, N_{2} dilution, methane fuel, propane fuel) were analyzed to explore the sources of CCV. This first comparison between simulated and measured combustion data is conducted at a robust combustion condition: 40 kPa Manifold Absolute Pressure (MAP), 1300 rpm, stoichiometric premixed propaneair reactants, and no N_{2} dilution. Ignition was at 342 crank angle degrees after TopDeadCenter exhaust (°aTDCe). The 40 kPa/101.5 kPa intake/exhaust manifold absolute pressures result in backflow during exhaustvalve opening and valve overlap periods. Details of the engine and optical methods can be found in [7]. Engine and operating parameters key to this paper are provided in Table 1, and details are available at [21].
Incylinder pressure measurements were recorded every 0.5 CrankAngle Degrees (CAD) of rotation. The AHR analysis was performed on both the measured and simulated pressure traces using the same crankangle resolution and similar values of the polytropic coefficients n; for the experiments, individualcycle values of n were used [7], whereas averaged values were used for the LES. The instantaneous heat release for crankangle interval i is given by (1)
The cumulative heat release (Q) is then computed by integrating equation (1), and the AHRbased massburned fraction is given by MBF = Q/Q_{total}, where Q_{total} = Q_{max}.
Imaging of OH* chemiluminescence (Fig. 1) was used to measure the area of the twodimensional projection of early flame growth with 0.5 × 0.5 mm spatial resolution. Highspeed imaging was used to acquire images every 2 CAD from SOIgn (342 °aTDCe) to topdeadcenter combustion (corresponding to an average MBF of approximately CA10), at which time the flames had grown outside of the field of view. The OH* signal was filtered by a combination of optical filters (UG11 and a 330 nm shortpass filter) that is centered on 307 nm with a 25 nm wide 10% cutoff to isolate the 306.4 nm band from emissions of other major combustion products over 330 nm [26]. To account for spatial differences in intensifier efficiencies, all acquired images were first whitefield corrected and a spatial calibration was applied. Then a binarization threshold was determined by the average value of the maximum values in each pixel column. The parameterization of this threshold improves capturing the burnedgasarea region and allows a higher degree of automation. The flamekernel growth rate can be quantified in this manner prior to the time when acceptable pressure measurements are available, and the measured growth rates then can be compared with the simulated flamekernel growth rates.
Highspeed PIV was used to measure the twodimensional, twocomponent velocity in the plane y = −4.5 mm near the spark plug (Fig. 1). The present analysis focused on data from a 6 × 12 mm region (x:z = −2 → −8 mm: −8 → +4 mm), as shown with the white box in Figure 1. The spatial resolution was 1 mm, and was determined by the 1 mm thick laser sheet and the 1 × 1 mm (32 × 32 pixel) interrogation spots on a 0.5 × 0.5 mm grid (50% interrogationspot overlap). The velocity dynamic range was 0.6–25 m/s, as determined by the 0.2to8pixel PIV resolution. The velocity gradients were computed using central differencing on the grid interior, with onesided differencing at the boundaries. The measurement plane was offset by 4.5 mm from the spark gap (y = −4.5 mm, Fig. 1) to avoid scattering from the sparkplug elements. Details of the PIV measurement procedures are provided in [7].
Figure 1 Schematic illustration of the TCCIII engine showing the coordinate system, PIV measurement plane (y = −4.5 mm) and the flame imaging fieldsofview. The white box in the PIV fieldofview image shows a 12 × 6 mm region near the spark plug (x = 4 to −8 mm, z = −2 to −8 mm) in which flow parameters have been averaged in the subsequent analysis. 
TCCIII engine geometry, valve timing and operating condition.
2 Numerical Simulations
LES has been performed using STARCD v4 [27]. The computational domain includes the incylinder region, intake and exhaust ports, and intake and exhaust plenums. The average incylinder mesh size is approximately 0.75 mm, with finer meshes near the sparkplug tip and in the valvecurtain regions. Timedependent pressure and temperature boundary conditions from a onedimensional model are prescribed at the plenum inlet/outlet, with no CCV [19]. A standard Smagorinsky subfilterscale turbulence model is used. Results from two different ignition and turbulent flame propagation models (designated as Model 1 and Model 2) are presented, with all other numerical, flow and heat transfer modeling aspects being the same. For each combustion model, 60 engine cycles have been simulated. For Model 1, three blocks of 20 consecutive cycles were simulated, whereas for Model 2, four blocks of 15 consecutive cycles were simulated. The same initial conditions were used for both models, which were obtained from different cycles of a preliminary LES study where a different subfilterscale turbulence model and turbulent combustion model were used. No patterns of correlation between different blocks of cycles were evident, as shown in Figure 2.
Figure 2 Maximum incylinder pressure versus cycle number for Model 1 (top: 3 sets of 20 consecutive cycles) and Model 2 (bottom: 4 sets of 15 consecutive cycles). 
2.1 Model 1
Model 1 uses a simple energy deposition model for ignition, and a modified Thickened Flame Model (TFM) [28] implemented by the authors for turbulent flame propagation. The energy deposition model creates a spherical temperature source with a radius of 0.75 mm, a temperature of 3000 K, and a duration of 4.0 crank angle degree. This hightemperature source initiates a flame based on a twostep global propaneair chemical mechanism [11], which has been tuned to give the correct laminar flame speed and adiabatic flame temperature for the engine conditions of interest. After the flame initiation, the thickened flame model takes over to propagate the flame. The essential features of the thickened flame model can be understood in terms of classical thermal theory for laminar premixed flames, according to which the laminar flame speed is proportional to the square root of the product of the molecular diffusivity D and the reaction rate ω, and the laminar flame thickness is proportional to the square root of the quotient of D and ω: (2)
The first essential idea in TFM is to multiply the diffusivity D and to divide the chemical source term ω in the species and enthalpy equations by the same factor F. From equation (2), it follows that this increases the flame thickness by the factor F (for F > 1), while the laminar flame speed (or equivalently, the mass rate of reactant consumption) remains unchanged. This suggests an approach by which one can recover the correct burn rate in an “underresolved” numerical simulation of a flame where the underlying flame structure is not fully resolved: choose F to be the ratio of the actual mesh size to the mesh size that is needed to fully resolve the underlying flame structure (typically ∼0.1 ). In LES of turbulent flames, an additional factor is needed to account for the wrinkling of the flame (thereby enhancing the local burn rate) by unresolved turbulence scales. The second essential element in TFM, then, is to maintain the correct overall burn rate (or turbulent flame speed) in a numerical simulation of a turbulent flame, even when all turbulence scales are not fully resolved. This effect is implemented through a multiplicative factor Ξ that is applied to both the diffusive term and the chemical source term in the species and enthalpy equations. Multiplying both D and ω in equation (2) by the same factor Ξ increases the flame speed (or mass rate of reactant consumption) by the factor Ξ (for Ξ > 1), without changing the flame thickness. In contrast to F, Ξ has a rather complicated functional form that depends on local turbulence parameters, and that is informed by theoretical analysis and by direct numerical simulations. The resulting species mass fraction (Y_{m}) equation is then: (3) and similarly for enthalpy. The factors F and Ξ are applied only within the turbulent flame brush. Several models for Ξ have been formulated. Charlette et al. [29] proposed the following powerlaw: (4)
Here is the combustion filter size, β is a model constant, the function Γ represents the integration of turbulent strain rate induced by scales smaller than Δ (which is fitted by a function [29]), and is the subfilterscale velocity fluctuation. An expression for the latter was proposed by Colin and Truffin [30]: (5) where is a model constant. The resolved velocity gradients are felt by the turbulent combustion model through the Smagorinsky subfilterscale turbulence model and through . Two important modifications for applications of TFM in engines are to limit the turbulence length scale in vicinity of walls (the length scale cannot be larger than the distance to nearest wall), and to disable the unrealistically high rate of wall heat transfer that is otherwise implied by the model when the flame reaches a wall. Details can be found in [31].
2.2 Model 2
Model 2 uses the Imposed Stretch SparkIgnition Model (ISSIMLES) [30], implemented in STARCD by the authors and coupled with the ECFMLES turbulent combustion model [32]. The model is based on the theoretical ability of the flame surface density (FSD) approach to account for flame kernel growth to a fully developed flame, from the first appearance of a minimum selfsustained flame nucleus. The simplified model of the coil ignition system from Duclos and Colin [33] is adopted, accounting for electric energy discharge after the breakdown electrode voltage has been exceeded. The electric circuit of the TCCIII engine delivers 36.5 mJ per cycle, and this value is used in the calculations. The model reproduces the sparktogas energy transfer rate in the glow phase as a function of voltage and electric intensity, so that a finiteduration (noninstantaneous) ignition process is simulated. Energy losses to the electrodes lead to an energy transfer efficiency of approximately 60%, following Verhoeven's measurements [34]. Whereas the model for the electric circuit is able to account for restriking events, no information on this is available from the experiments. The spark plug is fully meshed in the 3D model, and the energy at breakdown is deposited at the midpoint of the electrode gap (1.15 mm). A corresponding amount of burned gas is initialized, equal to the mass of a cylinder with radius 2δ_{l}, the electrode gap as the initial length and unburnt density ρ_{u}. The modeling of a reduced size flame kernel is a challenging task for the FSDLES equation. For this purpose, a modified FSDLES equation was proposed in [32]. With r_{b} denoting the burnedgas (flame) radius and the flame surface density, the term accounts for subfilterscale stretch during ignition. As stated in [30], the hypothesis that the wrinkling factor Ξ = 1 is made, given the low turbulence intensity at ignition time. The modified FSDLES equation is: (6)
The transition from subfilterscaleonly to subfilterscale and resolvedscale contributions to flame area corresponds to the continuous transition from ignition to a fully propagating turbulent flame. The use of a unique transport equation for FSD from the early stages allows the simulation of ignition events occurring in regions of high flow velocity and/or velocity gradients. The transition is governed in the model by an algebraic α function, which depends on the filtered flame radius and is expressed as in equation (7). (7) The flame radius is calculated from the equivalent enflamed spherical volume V_{b}. This is defined as shown in equation (8); it is function of the massbased progress variable , where and are the Favreaveraged concentrations of fuel and fuel tracer: (8)
At ignition (a ≅ 0), the ignition source term is active and accounts for flame surface production due to subfilterscale stretch. As the enflamed volume grows and α increases from 0 to 1, the ignition source term is gradually suppressed, while the C_{sgs}, , C_{res} and S_{res} terms are activated. These are the subfilterscale and resolvedscale curvature, the flame surface density production and the resolved strain contributions to the FSD equation. These terms restore the standard flamesurfacedensitybased LES turbulent flame propagation equation. The use of a transport equation for the entire flame history starting from ignition allows for full interaction between the local flow conditions (resolved velocity and velocity gradients) and the initial flame growth. The gradient of the resolved velocity affects the resolved strain , while the subfilterscale velocity fluctuation u′ is accounted for through the efficiency function that was proposed in [28], and which is included in the subfilterscale strain term S_{sgs}. In equation (9), Δ is the filter length, and s_{l} and δ_{l} are the laminar flame speed and thickness, respectively: (9)
Finally, a regionspecific wall temperature of 1000 K is imposed for the spark plug metal electrodes, given the impact of heat transfer from the initially small enflamed kernel to the surrounding solid parts.
3 Results
The comparisons between experiments and simulations are organized into four subsections: analysis of thermodynamic conditions at time of ignition, analysis based on apparent heat release and related quantities, analysis based on experimental flame imaging and analogous quantities extracted from the simulations, and analysis of the velocity fields in the vicinity of the sparkplug gap at the time of ignition.
3.1 Thermodynamic Analysis
The first comparison between experiments and simulations is the thermodynamic state of the mixture at SOIgn, 342 °aTDCe. The log P−log V plots of Figure 3 show 60 cycles of experimental data and 60 cycles from each model through compression, combustion and expansion. This provides a comparison of the volumeaveraged thermodynamic states experienced in the experiments versus those experienced by the two models. Looking first at the compression stroke, it can be seen that there is imperceptible CCV in both the measurements and the simulations, and that the measured data suggest a larger value for the polytropic coefficient n compared to the models: n_{Meas} = 1.33, n_{Mod1} = 1.31, n_{Mod2} = 1.31. There is very little droop in the expected straight lines, as indicated by the very high values of the measured and simulated R^{2}. The lower polytropic coefficients for the LES can only be caused by different trapped mass or by heat loss, since there is no blowby in the LES and the air/fuel composition is well known. The unknown temperature of the quartz cylinder in the experiments and uncertain wallheattransfer model in the LES are the most likely reasons for the disagreement between the measured and simulated polytropic coefficient values. It can be estimated from the polytropic compression that this difference in the values of n would result in a 6% lower bulk temperature for the LES compared to experiment (out of approximately 600 K) at SOIgn = 342 °aTDCe.
The pressure CCV near TDC are also revealed in Figure 3. From these, it can be seen that some of the LES cycles burn quite late, and that in turn yields higher CCV in the expansion stroke in the simulations compared to the experiments.
Figure 3 log P versus log V diagrams of the compression, combustion and expansion strokes for the measured data (top), Model 1 (center) and Model 2 (bottom). 
3.2 Apparent Heat Release Analysis
The combustion variability is more evident in Figure 4, where measured and computed Apparent Heat Release Rate (AHRR) traces are compared. There, 60 measured cycles are shown, so that the sample size is the same as in the simulations. It will be shown later that these 60 measured cycles are representative of the 1500 total available experimental cycles (500 cycles each from three separate tests). The measured results in Figure 4 show that combustion has completed by 30 CAD after TDC for all cycles. This was suggested by Ayala and Heywood [4] to be a necessary condition for stable engine operation. Both LES models exhibit longer combustion durations and higher CCV compared to what is observed in the experiments. In particular, the simulations (especially Model 1) show some extremely slow cycles that are not seen experimentally for this operating condition.
It is customary to use the MBF obtained from the pressurebased cumulative AHR, to quantify the 10, 50 and 90percent MBF periods (CA10, CA50 and CA90). Figure 5a compares the individualcycle MBF based on the cumulative AHR with the MBF based on the actual fuel consumption for LES Model 1, and Figure 5b shows the corresponding CA10, CA50 and CA90 values. Results for Model 2 are similar (not shown). Together, Figure 5a and b demonstrate that the pressurebased AHRR analysis of equation (1), which has typically been applied to experimentally measured data, is a precise metric that accurately reflects the actual fuel burn rate. It is therefore useful as a single metric that can be applied to both the measured and simulated data.
Figure 6 shows the range of CA10 CCV, and the correlations between CA50, CA90, and the Location of Peak Pressure (LPP) with CA10. Linear correlations of CA50, CA90 and LPP with CA10 are evident in both the experiments and the LES. The experiments exhibit a consistently smaller range of values for CA10 and CA90 (lower CCV) compared to the simulations. However, in all cases the burning after CA50 becomes more variable, as indicated by the lower correlation of CA10 with CA90 compared to CA50. Although the ranges of the values from the simulations overlap with the measured values, it is clear that both models exhibit some cycles with considerably slower burns compared to the measurements. More importantly for the present purposes, Figure 6 demonstrates that the CCV of combustion is highly correlated with the early flame kernel development time, even for the verylateburning simulated cycles (evident from the R^{2} values). Therefore, it is essential to correctly capture the early flame development period in the simulations, and to better understand what is happening during this critical period. The subsequent comparisons between experiments and simulations focus on the early flame development.
While it is possible to compute the AHRR from measured pressured data prior to CA10, experimental values of AHRR at very early times can be unreliable because of inaccuracies and low signaltonoise ratios. Figure 7a shows this breakdown for the very early burn (CA005). Figure 7b shows that CA01 correlates linearly with CA10 in both the experiments and the simulations. As already noted, the experiments are difficult for very early times. And in the LES models used here, factors such as discharge repeatability and heat transfer to electrodes, are not fully taken into account; therefore, a correlation between CA01 and CA10 might be expected in the simulations, to some extent. The arbitrary imposition of a cycleindependent ignition energy might reduce the CA01/CA10 correlation, but possibly for the wrong reason. With these caveats, the results suggest that the phasing of CA10 (and from what has been shown earlier, the phasing of the entire combustion event) has largely been established by CA01 for this operating condition. Thus, it is of interest to explore additional measures to dissect the ignition and flame propagation processes prior to CA01. For this purpose, flame imaging was used to capture the earlyflame growth, and to estimate the notional laminartoturbulent transition time, τ_{lamturb}, suggested by Abraham et al. [17], as described in the following subsection.
Figure 4 Measured and computed individualcycle AHRR versus crank angle. Top: measured. Center: Model 1. Bottom: Model 2. 
Figure 5 Comparisons of MBFs based on the cumulative AHR versus on the fuel consumption rate, for Model 1. 
Figure 6 Scatter plots of CA50, CA90 and Location of Peak Pressure (LPP) versus CA10. Left: Measured. Center: Model 1. Right: Model 2. 
Figure 7 (a) Measured MBF at crank angles prior to CA10, demonstrating limitations of pressurebased AHRR analysis prior to CA01. (b) Comparisons of CA01 and CA10 distributions between simulated (Model 1 and Model 2) and measured (labeled “UM”) cycles. 
3.3 Flame Imaging Analysis
The flame growth was quantified here as the projected 2D area of the 3D burnedgas volume, as measured from the OH* images and the simulations (Fig. 8a) every 2 CAD. For the projected flame area in the simulations, the 3D flame surface was defined as the location of 50% local mass burned fraction, and the 2D projected area was computed as in the experiments. This is a somewhat lessthanperfect comparison between experiments and simulations, because of the missing outofplane flamegrowth information. Nevertheless, it is of interest to see whether or not this relatively simple approach reveals anything of interest. Applying the concept of Arpaci et al. [18], the transition time τ_{lamturb} is defined here as the time between SOIgn and max(d^{2}A_{flm}/dt^{2}), as shown in example of Figure 8b; therefore the rolloff of the flame growth after 1 ms is due to the limited experimental field of view. Figure 8c shows that this transition occurs well before CA01.
It is recognized that the single view that is available experimentally is limited in its ability to capture outofplane growth, for example. However, this metric is intended only to identify the overall growth rate, not the absolute volume. This metric was found to correlate with CA10 and to discriminate between the stoichiometric and rich/lean/dilute mixtures of the experimental measurements of Schiffmann [7]. Figure 9 shows that there is a weak linear correlation between τ_{lamturb} and CA10 (the bimodal distribution of growth rates is believed to be an artifact of the low sampling rate (2 CAD) compared to the growth rate). This may be due to multivariable influences as well as the single view of the measurement. Nonetheless, τ_{lamturb} identifies the transition from one growth rate to another, and in principle can be applied equally to both the measured and simulated data. However, a sharp transition in growth rate was not observed in the simulations for either model (not shown).
To further dissect the EFD period, the measured cycletocycle sparkduration distribution is compared to the measured duration of CA10 (CA00–CA10) in Figure 9. The plasma duration distribution is presented here as a reminder that both electrical and chemical energies contribute to the early flame kernel during this period. Figure 9 demonstrates that the glowdischarge ignition system used here creates a plasma whose duration is longer than τ_{lamturb}, and that it lasts until nearly CA10. This suggests that the flame growth measured by the OH* emissions have commenced well before the plasma has ended. Furthermore, neither spark duration nor τ_{lamturb} correlate strongly with CA10. This figure also illustrates that the 60 measured cycles that were selected to compare with the 60 computed cycles from each model are representative of the 1500 available measured cycles.
A length scale that is expected to be relevant for the early flame development is computed as the diameter of a circle that has the same projected 2D area as that of the flame kernel. Figure 10 shows the ensembleaveraged equivalent diameter of the flame as a function of time after start of ignition, for the experiments and the two LES models. The range of the CCV in this quantity (dotted lines) is quantified by the minimum and maximum areas observed in the three data sets. Also shown are the average values of τ_{lamturb}, CA01, End of Spark (EoSpk), and CA10. From this figure, estimates of the relevant flamekernel length and timescales during the period between SOIgn and CA10 can be deduced. Since τ_{lamturb} has been shown to be poorly correlated with CA10, while CA01 is highly correlated with CA10, it can be deduced that the correlation with CA10 develops in the period between 0.4 ms and 1.5 ms (CA01), where the flame diameter is between 5 and 15 mm. This can be compared to the estimated integral time scales of 1–1.3 ms and integral length scales of L_{ux}, L_{uz}, L_{wz} and L_{wx} equal to 10–20 mm, 8 mm, 3 mm and 5 mm, respectively [7], where L_{i}_{,j} denotes the correlation length for the ith velocity component in the jth direction. This then suggests that simulated and measured flow properties at these scales could be compared, to develop further insight into the early flame development.
Figure 8 (a) Examples of measured and simulated 2D projected flame kernel area. A singlecycle example of the burnedgas area evolution for a measured cycle is used to quantify the laminarturbulent transition time, τ_{lamturb}, shown by the reddiamond in (b), which occurs much earlier than the 1% MBF from the AHR analysis (c). 
Figure 9 Measured sparkdischarge duration and τ_{lamturb} distributions from 1500 cycles (500 cycles from each of three independent tests − in black), and the 60 sampled cycles (in red) that are used to provide a sample size that is the same as that available from the simulations. 
Figure 10 Ensembleaveraged (solid lines) and cycletocycle range (dotted lines) of the flame's equivalent diameter for all measured (1500) and simulated (60 each) cycles. The vertical lines mark crank angles of interest, from which the temporal and spatial scales of the early flamekernel growth can be deduced. 
3.4 Local Flow at Spark Gap Analysis
Comparisons between computed and measured mean and rms velocity fields for motored and firing conditions are reported in [31], and the reader is referred there for detailed quantitative comparisons. An example of mean velocity fields on two cutting planes that pass close to the spark gap at an instant shortly before SOIgn are shown in Figure 11, to give a general impression of the largescale flow structure. There are quantitative differences between the two LES models, between LES and PIV, and between the motored and firing preignition flow fields, but the general flow patterns are similar.
Multiple earlier experimental and simulation studies have shown that local flow conditions in the vicinity of the spark gap at the time of ignition affect the early flame kernel development and subsequent combustion process [7–13]. In particular, flame plasmas and kernels can be advected away from the spark plug, which affects the combustion phasing. Thus, it is prudent to compare not only the thermodynamic state between the experiments and simulations at the time of ignition (Fig. 3), but also the velocity and velocity gradients. Scatter plots of velocity magnitude, V, and the direction of the velocity in the 2D measurement plane (Fig. 1), u and w, are shown in Figure 12. Both the measured and simulated values have been spatially averaged over the 12 mm × 6 mm region illustrated in Figure 1. Therefore, these represent lowpass spatially filtered values (over a spatial scale of ∼10 mm) in the vicinity of the spark plug. As noted earlier, the measurement plane does not pass directly through the sparkplug gap (it is offset by 4.5 mm). Nevertheless, the measured values of V were found to correlate with CA10 in the statistical analysis of Schiffmann [7]. This offset does not detract from the point of the comparison here, which is to determine whether or not the simulations are predicting similar velocities near the spark plug. Figure 12 shows that the range of V values from the simulations tend to lie within the lower half of the measured range of values. Furthermore, the direction of the largescale flow (as indicated by the trend in the w versus u plot) is consistent between the experiment and simulations, although again, the overall lower velocity magnitudes in the simulations are evident.
Since Figure 12 demonstrates that there are differences between the measured and simulated largescale velocity magnitudes near the spark plug, it is reasonable to assess if the measured and simulated resolvedscale velocity gradients are similar or significantly different. The resolvedscale velocity gradients play an important role in subfilterscale turbulence and turbulent combustion modeling, as discussed earlier in the modeling section. In Figure 13, measured and computed distributions of the twocomponent resolvedscale velocity gradient components over the 12 mm × 6 mm measurement region displaced 4.5 mm from the sparkplug gap (the same region for which areaaveraged velocities are shown in Fig. 12) are plotted. These distributions were sampled for all 60 simulated cycles and for the 60 representative measured cycles. The measured (PIV) velocity gradients were computed directly from 1 mmresolution interrogations spots on a 0.5 mm grid, whereas the simulation distributions were sampled on the computationalgrid and subsequently mapped onto the PIV grid. The most striking feature of Figure 13 is the significantly narrower distribution of measured velocity gradient values compared to the simulated velocity gradient values.
To assess if the larger gradients from the simulations are the result of higher computationalgrid resolution compared to the experimental spatial resolution, a linear correlation was performed between the gradients mapped from the computational grid to the PIV grid and gradients computed from the remapped computed velocities. The linear correlation demonstrated that the gradients computed from the remapped velocity were only 8% lower than the higherresolution simulation gradients, with a linear correlation coefficient of 0.98. Thus, it was concluded that the narrower distributions of the measured gradients are not due to the remapping. Rather, Figure 13 reveals the limited ability of PIV measurements to resolve steep gradients. This may be the result of three factors: the error between the particles' motion and the local gasphase flow velocity; the limited dynamic range of the instantaneous velocity measurements (0.2 pixels out of 8 = 40×); and the difficulty of PIVcorrelationbased peakdetection in regions of steep velocity gradients (vector drop out). Thus, based on these measurements, it is not possible to conclude if the simulations accurately predict the resolvedscale gradients in the flow. Nevertheless, since the two LES models yielded different combustion phasing (e.g., Fig. 4), it is of value to compare the computed gradient distributions between the two models.
To assess the velocitygradient distributions that the LES models present to the ignition/combustion models, all nine velocitygradienttensor components were sampled and compared from the two simulations in a threedimensional volume around the spark plug at the time of ignition. The distributions were sampled from all 60 cycles, on a 0.5 mm grid in a rectangular prism volume centered at the sparkplug gap: x = −8 mm → 4 mm, y = −4.5 mm → 4.5 mm, z = −8 mm → −2 mm. In addition, the distributions from the 10 fastest burning and 10 slowest burning cycles are plotted in Figure 14. These distributions are the fraction of the total sample in the histograms (bin size of 250 1/s) normalized by the number of samples. Thus, the fraction of the total sample (370k samples) can be compared to the 37k samples of the 10 fastest and slowest cycles, respectively. The distributions displayed in Figure 14 show that Model 2 has more lowmagnitude gradients (−1000/s to +1000/s) compared to Model 1; Model 1 yields broader distributions with larger excursions from the mean. Figure 14 also shows slight variations in skewness between the two simulations. As demonstrated by Petersen and Ghandhi [35], the turbulent dissipation in engine flows is expected to occur in relatively localized highshear regions between larger structures. Thus, the occurrence of the largest gradients from Eulerian measurements is expected to be infrequent compared to the smaller magnitude gradients at the interior of the larger structures. This is consistent with the distributions in Figure 14.
To better reveal the large dynamic range and relative frequency of occurrence of highmagnitude gradients, it is of value to replot the fractional frequency in Figure 14 using a logarithmic scale (Fig. 15). Figure 15 more clearly reveals the comparatively smaller range of velocity gradients presented to the ignition and turbulent combustion models by Model 2 compared to Model 1. The reason for this difference is unclear, at this point. It is perhaps more revealing to compare the fastest and slowest cycles from each model. Figures 14 and 15 show that for both models, the slowest cycles correspond to a narrower distribution of velocity gradient magnitudes (fewer samples with high velocity gradient magnitudes) compared to the mean distribution, while the fastest cycles correspond to a broader distribution of gradient magnitudes (more samples with high velocity gradient magnitudes) compared to the mean. Higher velocity gradients may lead to accelerated flame growth, and exposing the early kernel to a broader range of fluctuating velocity gradients may also contribute. Conversely, lower velocity gradient fluctuations and the absence of high velocity gradient magnitudes result in slower early flame growth.
Figure 11 Comparisons between measured and simulated mean velocity fields on two cutting planes at 340 °aTDCe (2° before SOIgn). 
Figure 12 Comparisons between measured and simulated velocity magnitude V distributions with the duration of CA10 (left), and 2D velocity components u and w showing distribution of the velocity direction (right). In all cases, the velocities have been spatially averaged over a 12 mm × 6 mm region in the vicinity of the sparkplug gap (y = −4.5 mm, Fig. 1) at SOIgn (342 °aTDCe). For the experiments, all 1500 cycles and 60 representative cycles are shown, the latter to provide a clearer visual comparison with the 60cycle LES results. 
Figure 13 Measured and computed twocomponent velocity gradient distributions in a measurement region in the vicinity of the spark plug (Fig. 1) at SOIgn (342 °aTDCe). Distributions of all four 2D velocity gradient components are shown. 
Figure 14 Distributions of nine components of the velocity gradient tensor in the volume centered at the spark plug at SOIgn (342 °aTDCe) for two different LES models. The distributions are sampled at all volume nodes in the region for the 60 simulated cycles (total sample size of 370k). Three distributions are shown for each component: all 60 cycles, the 10 fastestburning cycles, and the 10 slowestburning cycles. 
Figure 15 Distributions of nine components of the velocity gradient tensor in the volume centered at the spark plug at SOIgn (342 °aTDCe) for two different LES models. The distributions are sampled at all volume nodes in the region for the 60 simulated cycles (total sample size of 370k). Three distributions are shown for each component: all 60 cycles, the 10 fastestburning cycles, and the 10 slowestburning cycles. This is the same as Figure 14, but plotted using a logarithmic scale. 
Conclusion
Pressurebased combustion analysis, optical diagnostics and LES have been used to generate new physical insight into the early stages of combustion in homogeneouscharge sparkignition engines. The emphasis has been on developing strategies for making quantitative comparisons between highspeed/highresolution optical diagnostics and LES using common metrics for both the experiments and the simulations, focusing on the important early flame development period.
These are the first comparisons between combustion LES and experiments for the TCC engine, an engine configuration that is freely accessible to other research groups. Results from two different LES ignition turbulent flame propagation models have been presented, using the same numerical methods and computational mesh. The results reveal strengths and limitations of the experimental diagnostics and LES models. Specific findings are as follows:

Local flow conditions in the vicinity of the spark gap at the time of ignition play a significant role in determining the subsequent combustion event for the cycle. This largely confirms results from multiple earlier experimental and simulation studies.

Two widely used LES models yield rather different combustion behavior compared to experiments. Both models show somewhat lower velocity magnitudes in the vicinity of the spark gap at the time of ignition compared to experiment, and both models show higher combustion CCV compared to experiment. In particular, the models give some very slowburning cycles that are not observed in the experiments. Nevertheless, the comparisons between models, and between models and experiments, provide useful insight into the early flame development and into the limitations of the models and the experiments.

In the simulations, fastburning cycles correspond to cycles where there is a broader distribution of velocity gradient magnitudes (more samples with high velocity gradient magnitudes) in the vicinity of the sparkplug gap at the time of ignition compared to the average distribution. Conversely, slowburning cycles correspond to cycles where there is a narrower distribution of velocity gradient magnitudes in the vicinity of the sparkplug gap at the time of ignition compared to the average distribution. Higher velocity gradient magnitudes may lead to accelerated flame growth, and exposing the early kernel to a broader range of fluctuating gradients may also contribute. To explore this experimentally, it will be necessary to modify the PIV technique that has been used here to better resolve steep local velocity gradients.

The fate of the combustion event for each cycle is largely determined by what happens between the laminartoturbulent transition time and CA01. This suggests a range of temporal and spatial scales over which subsequent comparisons between experiments and simulations should be performed.
Acknowledgments
Coauthors from the Pennsylvania State University and the University of Michigan gratefully acknowledge funding from GM Global R&D. Coauthors from the Pennsylvania State University and the University of Modena and Reggio Emilia gratefully acknowledge Siemens PLM Software for providing access to, and support for, their STARCD CFD code.
References
 Ozdor N., Dulger M., Sher E. (1994) Cyclic variability in spark ignition engines: a literature survey, SAE Technical Paper 940987. [Google Scholar]
 Johansson B. (1996) Cycle to cycle variations in S.I. engines − the effects of fluid flow and gas composition in the vicinity of the spark plug on early combustion, SAE Technical Paper 962084. [Google Scholar]
 Pajot O. (2000) Etude expérimentale de l'influence de l'aérodynamique sur le comportement et la structure du front de flamme dans les conditions d'un moteur à allumage commandé, Thèse, University of Orléans, Orléans, France. [Google Scholar]
 Ayala F.A., Heywood J.B. (2007) Lean SI engines: the role of combustion variability in defining lean limits, SAE Technical Paper 2007240030. [Google Scholar]
 Lacour C., Pera C. (2011) An experimental database dedicated to the study and modelling of cyclic variability in sparkignition engines with LES, SAE Technical Paper 2011011282. [Google Scholar]
 Baum E., Peterson B., Böhm B., Dreizler A. (2014) On the validation of LES applied to internal combustion engine flows, Part 1: comprehensive experimental database, Flow Turb. Combust. 92, 1–2, 269–297. [Google Scholar]
 Schiffmann P. (2016) Root causes of cycletocycle combustion variations in spark ignited engines, PhD thesis, University of Michigan, Ann Arbor, MI, USA. [Google Scholar]
 Truffin K., Angelberger C., Richard S., Pera C. (2015) Using largeeddy simulation and multivariate analysis to understand the sources of combustion cyclic variability in a sparkignition engine, Combust. Flame 162, 12, 4371–4390. [CrossRef] [Google Scholar]
 Bates S.C. (1989) Flame imaging studies of cyclebycycle combustion variation in a SI fourstroke engine, SAE Technical Paper 892086. [Google Scholar]
 Pera C., Knop V., Reveillon J. (2015) Influence of flow and ignition fluctuations on cycletocycle variations in early flame kernel growth, Proc. Combust. Inst. 35, 2897–2905. [CrossRef] [Google Scholar]
 Granet V., Vermorel O., Lacour C., Enaux B., Dugué V., Poinsot T. (2012) Largeeddy simulation and experimental study of cycletocycle variations of stable and unstable operating points in a spark ignition engine, Combust. Flame 159, 1562–1575. [CrossRef] [Google Scholar]
 Goryntsev D., Sadiki A., Klein M., Janicka J. (2009) Large eddy simulation based analysis of the effects of cycletocycle variations on airfuel mixing in realistic DISI ICengines, Proc. Combust. Inst. 32, 2759–2766. [CrossRef] [Google Scholar]
 Goryntsev D., Sadiki A., Janicka J. (2013) Analysis of misfire processes in realistic direct injection spark ignition engine using multicycle large eddy simulation, Proc. Combust. Inst. 34, 2969–2976. [CrossRef] [Google Scholar]
 Heywood J.B. (1988) Internal combustion engine fundamentals, McGraw Hill, New York. [Google Scholar]
 Matekunas F.A. (1983) Modes and measures of cyclic combustion variability, SAE Technical Paper 830337. [Google Scholar]
 Beretta G.P., Rashidi M., Keck J.C. (1983) Turbulent flame propagation and combustion in spark ignition engines, Combust. Flame 51, 217–245. [CrossRef] [Google Scholar]
 Abraham J., Bracco F.V., Reitz R.D. (1985) Comparisons of computed and measured premixed charge engine combustion, Combust. Flame 60, 309–322. [CrossRef] [Google Scholar]
 Arpaci V.S., Ko Y., Lim M.T., Lee H.S. (2003) Spark kernel development in constant volume combustion, Combust. Flame 135, 315–322. [CrossRef] [Google Scholar]
 Liu K., Haworth D.C., Yang X.S., Gopalakrishnan V. (2013) Largeeddy simulation of motored flow in a twovalve piston engine: POD analysis and cycletocycle variations, Flow Turb. Combust. 91, 373–403. [Google Scholar]
 Schiffmann P., Gupta S., Reuss D., Sick V., Yang X., Kuo T.W. (2016) TCCIII engine benchmark for largeeddy simulation of IC engine flows, Oil Gas Sci. Technol. – Rev. IFP 71, 1, 1–27. [Google Scholar]
 Sick V., Reuss, D.L., Yang X., Kuo T.W. (2014) https://deepblue.lib.umich.edu/handle/2027.42/108382. [Google Scholar]
 Reuss D.L. (2000) Cyclic variability of largescale turbulent structures in directed and undirected IC engine flows, SAE Technical Paper 2000010246. [Google Scholar]
 Abraham P.S., Yang X., Gupta S., Kuo T.W., Reuss D.L., Sick V. (2016) Flowpattern switching in a motored spark ignition engine, SAE Int. J. Eng. Res. 16, 323–339. [Google Scholar]
 Montorfano A., Piscaglia F., Onorati A. (2015) An extension of the dynamic mesh handling with topological changes for LES of ICE in OpenFOAM, SAE Technical Paper 2015010384. [Google Scholar]
 Ameen M.M., Yang X., Kuo T.W., Xue Q., Som S. (2015) LES for simulating the gas exchange process in a spark ignition engine, in: ICEF20151002, Proceedings of the ASME 2015 Internal Combustion Engine Division Fall Technical Conference, Houston, TX, 8–11 November, 2015. [Google Scholar]
 Zhao H., Ladommatos N. (2001) Engine Combustion Instrumentation and Diagnostics, SAE International, Warrendale, PA, p. 842. [CrossRef] [Google Scholar]
 Siemens PLM Software (2016) http://mdx.plm.automation.siemens.com/starcd. [Google Scholar]
 Colin O., Ducros F., Veynante D., Poinsot T. (2000) A thickened flame model for large eddy simulations of turbulent premixed combustion, Phys. Fluids 12, 1843–1863. [CrossRef] [Google Scholar]
 Charlette F., Meneveau C., Veynante D. (2002) A powerlaw flame wrinkling model for LES of premixed turbulent combustion. Part I: nondynamic formulation and initial tests, Combust. Flame 131, 159–180. [CrossRef] [Google Scholar]
 Colin O., Truffin K. (2011) A spark ignition model for large eddy simulation based on an FSD transport equation (ISSIMLES), Proc. Combust. Inst. 33, 3097–3104. [Google Scholar]
 Shekhawat Y. (2017) Largeeddy simulations of motored flow and combustion in a homogeneouscharge sparkignition engine, Ph.D. thesis, The Pennsylvania State University, University Park, PA, USA. [Google Scholar]
 d'Adamo A., Breda S., Fontanesi S., Cantore G. (2015) LES modelling of sparkignition cycletocycle variability on a highly downsized DISI engine, SAE Int. J. Eng. 8, 2029–2041. [Google Scholar]
 Duclos J.M., Colin O. (2001) Arc and kernel tracking ignition model for 3D sparkignition engine calculations, COMODIA, 343–350. [Google Scholar]
 Verhoeven D. (1995) Spark heat transfer measurements in flowing gases, SAE Technical Paper 952450. [Google Scholar]
 Petersen B.R., Ghandhi J.B. (2010) High resolution scalar dissipation and turbulence length scale measurements in an internal combustion engine, SAE Technical Paper 2010010185. [Google Scholar]
All Tables
All Figures
Figure 1 Schematic illustration of the TCCIII engine showing the coordinate system, PIV measurement plane (y = −4.5 mm) and the flame imaging fieldsofview. The white box in the PIV fieldofview image shows a 12 × 6 mm region near the spark plug (x = 4 to −8 mm, z = −2 to −8 mm) in which flow parameters have been averaged in the subsequent analysis. 

In the text 
Figure 2 Maximum incylinder pressure versus cycle number for Model 1 (top: 3 sets of 20 consecutive cycles) and Model 2 (bottom: 4 sets of 15 consecutive cycles). 

In the text 
Figure 3 log P versus log V diagrams of the compression, combustion and expansion strokes for the measured data (top), Model 1 (center) and Model 2 (bottom). 

In the text 
Figure 4 Measured and computed individualcycle AHRR versus crank angle. Top: measured. Center: Model 1. Bottom: Model 2. 

In the text 
Figure 5 Comparisons of MBFs based on the cumulative AHR versus on the fuel consumption rate, for Model 1. 

In the text 
Figure 6 Scatter plots of CA50, CA90 and Location of Peak Pressure (LPP) versus CA10. Left: Measured. Center: Model 1. Right: Model 2. 

In the text 
Figure 7 (a) Measured MBF at crank angles prior to CA10, demonstrating limitations of pressurebased AHRR analysis prior to CA01. (b) Comparisons of CA01 and CA10 distributions between simulated (Model 1 and Model 2) and measured (labeled “UM”) cycles. 

In the text 
Figure 8 (a) Examples of measured and simulated 2D projected flame kernel area. A singlecycle example of the burnedgas area evolution for a measured cycle is used to quantify the laminarturbulent transition time, τ_{lamturb}, shown by the reddiamond in (b), which occurs much earlier than the 1% MBF from the AHR analysis (c). 

In the text 
Figure 9 Measured sparkdischarge duration and τ_{lamturb} distributions from 1500 cycles (500 cycles from each of three independent tests − in black), and the 60 sampled cycles (in red) that are used to provide a sample size that is the same as that available from the simulations. 

In the text 
Figure 10 Ensembleaveraged (solid lines) and cycletocycle range (dotted lines) of the flame's equivalent diameter for all measured (1500) and simulated (60 each) cycles. The vertical lines mark crank angles of interest, from which the temporal and spatial scales of the early flamekernel growth can be deduced. 

In the text 
Figure 11 Comparisons between measured and simulated mean velocity fields on two cutting planes at 340 °aTDCe (2° before SOIgn). 

In the text 
Figure 12 Comparisons between measured and simulated velocity magnitude V distributions with the duration of CA10 (left), and 2D velocity components u and w showing distribution of the velocity direction (right). In all cases, the velocities have been spatially averaged over a 12 mm × 6 mm region in the vicinity of the sparkplug gap (y = −4.5 mm, Fig. 1) at SOIgn (342 °aTDCe). For the experiments, all 1500 cycles and 60 representative cycles are shown, the latter to provide a clearer visual comparison with the 60cycle LES results. 

In the text 
Figure 13 Measured and computed twocomponent velocity gradient distributions in a measurement region in the vicinity of the spark plug (Fig. 1) at SOIgn (342 °aTDCe). Distributions of all four 2D velocity gradient components are shown. 

In the text 
Figure 14 Distributions of nine components of the velocity gradient tensor in the volume centered at the spark plug at SOIgn (342 °aTDCe) for two different LES models. The distributions are sampled at all volume nodes in the region for the 60 simulated cycles (total sample size of 370k). Three distributions are shown for each component: all 60 cycles, the 10 fastestburning cycles, and the 10 slowestburning cycles. 

In the text 
Figure 15 Distributions of nine components of the velocity gradient tensor in the volume centered at the spark plug at SOIgn (342 °aTDCe) for two different LES models. The distributions are sampled at all volume nodes in the region for the 60 simulated cycles (total sample size of 370k). Three distributions are shown for each component: all 60 cycles, the 10 fastestburning cycles, and the 10 slowestburning cycles. This is the same as Figure 14, but plotted using a logarithmic scale. 

In the text 