Study of a Flame Kernel Evolution in a Turbulent Mixing Layer Using LES with a Laminar Chemistry Model

Temporal evolution of an ignition kernel in a spark ignited turbulent hydrogen–air mixing layer is studied using the large eddy simulation approach with an implicit treatment of the reaction source terms. The applied numerical code is based on a high-order compact difference approximation combined with a weighted essentially non-oscillatory scheme, which provide an accurate resolution of the small scale phenomena. The spark is modelled by an energy deposition model coupled with the enthalpy equation. Since the fuel and oxidiser streams move in the opposite directions the flow is dominated by shear stresses and vortical structures. The ignition follows three different scenarios depending on the spark location. An interaction between the developing flame kernel and large turbulent structure is analysed starting from the energy transfer up to a fully reacting state. The size and shape of the flames are correlated with the initial ignition scenario. A 3D visualisation of the transient position of the flame kernel shows its different spatio-temporal behaviour. It is observed that the flame can stay close to the initial position and spread equally in all directions or it can move far from the initial location and follow the evolving flow field. In the latter scenario two separate sub-stages are convincingly identified. Concerning the ignition probability (Pign\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{ign}$$\end{document}), when the spark is located in the immediate vicinity of the mixing layer, the Pign\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{ign}$$\end{document} on the fuel-rich side is higher. On the other hand, when the ignition occurs far from the mixing layer, the Pign\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{ign}$$\end{document} is significantly larger on the lean side. Unlike the local composition the impact of the strain rate of the velocity field was found to have very limited impact on the Pign\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{ign}$$\end{document}.


Introduction
The complete understanding of the turbulent flame ignition is extremely important both from the fundamental as well as practical point of view. Safety and efficiency of technical combustion devices has been getting attention for last decades. A large number of industrial burners, engines and aircraft combustors are based on non-premixed systems 1 3 and combustion initiation through an externally introduced energy source (e.g. an electric spark, laser, plasma). Further development of these devices is a challenging task and it requires more experimental and numerical works on the ignition process and flame kernel evolution. A considerable amount of research in this field has been done to date but no univocal conclusions have been found. Most of these studies focus on spark ignition of jets (Birch et al. 1981;Smith et al. 1988;Ahmed and Mastorakos 2006;McCraw et al. 2007;Lacaze et al. 2009;Chen et al. 2017;Zhang et al. 2019), counter-flows (Ahmed et al. 2007) and mixing layers Chakraborty and Mastorakos 2008;Ahmed and Mastorakos 2016;EidiAttarZade et al. 2017;Tyliszczak 2018, 2019). Due to simple configurations, the analysis of these problems offers a deep insight into physics of the flame/turbulence interactions. Strongly unsteady phenomena occurring in the igniting flows including the flame kernel formation, its propagation and stabilisation are of current interest for Computational Fluid Dynamics (CFD) community. The fundamental findings on the ignition phenomenon, discussion of its probabilistic nature, and modelling aspects for gaseous and two-phase flows were presented in review papers by Mastorakos (2009Mastorakos ( , 2017. In the following, we focus only on some selected problems that are crucial for the present work. Spark ignition in axisymmetric turbulent jets was studied experimentally by Birch et al. (1981) and Smith et al. (1988) in the 1980s. These early experiments demonstrated that fluctuations of the mixture fraction are critical for the flame kernel generation process when the temperature rises and chemical reactions proceed. The results showed that the probability of successful flame kernel generation ( P ker ) and probability of finding flammable mixture ( P f ) at the centreline of the jets were approximately equal, P ker ≈ P f . Generation of the flame kernel does not necessarily mean the successful ignition, which is considered as a four-step process, i.e. energy deposition→flame kernel formation→propagation→stabilization. It was observed that the probability of ignition ( P ign ) can be lower than P ker (Smith et al. 1988). In the jet type flows the kernel formation is regarded as the process, that happens from the time moment of spark initiation up to the time when the downstream edge of the flame starts propagation with the flow and along the radial directions. Then, two different scenarios are possible: (1) a whole flame moves with the flow direction (blow-off); (2) its upstream edge stays for a short time at the initial location and begins the propagation against the incoming reactants-the successful ignition with the flame stabilized as lifted one or attached.
More recently, the differences between P f , P ker and P ign were reported by Ahmed et al. (2007) in a counter-flow configuration. They showed that P ker < P f resulted from the fact that in a region of strong mixing, the flame kernel could have been quenched due to local strain effects or large turbulent fluctuations. Hence, even if the flame was initiated at favourable mixture conditions, the probability of ignition was lower, i.e. P ign < P ker < P f . One should assume that the same situations can be observed for jet type flames. It is worth to note that Ahmed et al. (2007) observed ignition in locations of unfavourable mixture composition, i.e. in the places with P f = 0 . Their results showed that even when the spark was initiated outside of the ignitable region (e.g. at very lean/rich mixtures) the heat and radicals could be convected towards flammable regions and initiate the flame there. In such cases, the turbulent structures seem to be the most important, as on the one hand their presence may lead to enhanced mixing promoting the ignition process, while on the other, if the turbulence level is too high, the flame kernel can be quenched. In fact, depending on the spark location with respect to a vortical structure, three significantly different ignition scenarios (two successful and one failed) can be identified (Wawrzak and Tyliszczak 2019). It was found that in the case of successful ignition the initial scenario directly affects further 1 3 flame development phase. Worth noting is that not only the size and energy of the sparks can have impact on the success of ignition but also the size and strength of the turbulent structures present in the spark vicinity (Wawrzak and Tyliszczak 2018).
The spark ignition process in a methane jet with particular attention to a spatial distribution of the ignition probability was analysed experimentally by Ahmed and Mastorakos (2006). They considered varying spark parameters and jet velocities and found that: (1) downstream the axial location where the mean mixture fraction drops below the stoichiometric value ⟨ ST ⟩ the P ign decreases significantly; (2) the P ign is always smaller than one, even along ⟨ ST ⟩ isolines; (3) the region of high P ign decreases with an increasing jet velocity; (4) the P ign increases for a longer spark duration. The failure of the ignition process was attributed to the lack of the flame kernel or its blow-out at downstream locations, similarly as in Smith et al. (1988). In the work of Ahmed et al. (2007) the failed ignition events associated with the latter situation were observed for the spark locations beyond some critical distance above which the flame could not propagate against the incoming reactants. These findings were partially confirmed by McCraw et al. (2007) in a jet configuration with a pilot flame as the ignition source.
In numerical simulations of the spark ignition process, the spark is modelled one of two ways, i.e. by an energy (or temperature) source term added to the energy equation or by a volume of burning mixture, which, when placed in a flow domain, mimics the flame kernel. The former approach is often coupled both with direct numerical simulations (DNS) Chakraborty and Mastorakos 2008) as well as LES (Lacaze et al. 2009;Richardson et al. 2007;Jones and Tyliszczak 2010) with various combustion models. The latter method has been mostly used with LES (Rosiak and Tyliszczak 2016;Tyliszczak et al. 2014;Wawrzak and Tyliszczak 2017). The above mentioned experiment of Ahmed and Mastorakos (2006) was reproduced by Lacaze et al. (2009) using LES combined with the thickened flame combustion model (Colin et al. 2000;Jones and Prasad 2011) with the help of LES and the Eulerian stochastic field method (Valiño 1998), and more recently by Chen et al. (2017) using LES combined with the joint sub-grid PDF approach for partially premixed combustion and Zhang et al. (2019) with the help of LES and the Conditional Moment Closure model. In the first two studies, the spark was defined as the energy source term according to an energy deposition model (Lacaze et al. 2009), while in the latter, it was mimicked by a burning mixture introduced to the flow domains. In all above mentioned papers, the measured size and position (lift-off height) of the flame, as well as its propagation speed were predicted accurately. The study of Chen et al. (2017) revealed that the propagating flame follows a spiral-shaped trajectory. Spatial ignition probability was studied only by Lacaze et al. (2009). It was in good agreement with the measurements in a mixing layer close to the inlet, whereas further downstream from the inlet, the probability was significantly over-predicted.
Recent experimental and numerical research on the spark ignition in a shear-less mixing layer (methane and air) performed in Mastorakos (2016), EidiAttarZade et al. (2017) focused on a V-shaped region of the ignition probability. It was found that the highest P ign was located around ST . In lean mixtures, it was consistent with the lower flammability limit, whereas in the rich regions, the values of P ign were reduced to zero before reaching the location of the upper flammability limit.
To summarize, the above findings show that placing the spark in the region where P f = 0 does not preclude generation of the flame kernel and its subsequent growth. On the other hand, the location of the spark at the point with P f > 0 , where P ker > 0 , does not guarantee successful ignition. The established flame kernel may grow or may be too weak to propagate against the vortices in a highly turbulent flow. Definitely, the ignition process has a stochastic nature which is largely conditioned by flow regimes.
The present paper extends the research performed in Tyliszczak (2018, 2019). The objective of the study is to analyse the spark ignition in flow conditions dominated by shear stresses and the presence of strong vortical structures. We investigate the mechanism of flame development in a temporally evolving mixing layer formed between streams of fuel ( H 2 ∕N 2 ) and air flowing in opposite directions, similarly as in Hawkes et al. (2007) for a case of a temporally evolving jet. Usually, the successful spark ignition is considered as a multi-step process in which the main parts are (1) the initial kernel generation, (2) the flame propagation, and (3) stabilisation. In the present configuration, the stable or steady state conditions with a well-defined flame position cannot exist by definition, and therefore, we treat the ignition event as successful when the flame starts to propagate across the domain. We apply the LES method with a detailed chemistry for hydrogen oxidation and 'no model' approach for reaction terms closure, where the chemical source terms are computed directly using the filtered (resolved) temperature and species mass fractions. Such a 'no model' approach is sometimes labelled as 'implicit LES (ILES)' (Duwig et al. 2011), which should not be confused with an implicit approach for sub-filter terms resulting from the filtering operation applied to the convective parts of the governing equations. The two major problems analysed in this paper are: (1) the mechanism of flame kernel formation in the presence of convective forces and strong shear stresses; (2) the impact of turbulent structures on the initial phase of the flame evolution. We focus on two different successful ignition scenarios identified in Wawrzak and Tyliszczak (2019). Furthermore, we put attention to the ignition probability and try to assess its distribution depending on a spark location with respect to an initial shear layer position and a size and strength of the vortical structures formed in its vicinity.

Mathematical Approach
We consider a low Mach number reacting flow for which the continuity, Navier-Stokes and scalar transport equations, complemented with the equation of state, are defined within the framework of the LES method as: where the bar and tilde symbols denote filtered quantities (Geurts 2004). The variables appearing in the above equations are the velocity vector , the density , the hydrodynamic pressure p, the species mass fractions Y and enthalpy h. The term Q represents the energy source, whose role is to mimic a real spark. Its definition is discussed later. The subscript represents the index of the species = 1, … , N-species, the symbols p 0 and R are the thermodynamic pressure and the gas constant, respectively. The quantities and D , D are the viscous stress tensor and mass and heat diffusivities. The sub-filter tensor is given by SGS = 2 t , where is the rate of strain tensor of the resolved velocity field and t is the sub-filter viscosity modelled as in Vreman (2004). The sub-filter diffusivities in Eqs.

'No Model' Approach for the Reaction Terms
The main difference between LES with the combustion models (e.g. flamelet, thickened flame, conditional moment closure, etc.) and a 'no model' approach relies on treatment of the reaction source terms ẇ . In the latter approach, the sub-filter scales resulting from filtering of the convective terms are modelled as usually (in this work, the Vreman's model (Vreman 2004) is applied), but the chemical source terms ẇ are computed directly using the filtered variables: where the impact of the sub-filter unresolved scales is neglected. The 'no model' assumptions would certainly fail in the Reynolds Averaged Navier-Stokes (RANS) framework as the fluctuations in RANS models are large in general. On the other hand, in the laminar flow simulation and in DNS, where all turbulent flow scales are resolved, the 'no model' approach is perfectly valid. Hence, one may assume that for sufficiently dense computational meshes, when the sizes of grid cells are comparable with the Kolmogorov length scale, the use of this method is appropriate, as shown for instance in Duwig et al. (2011), Prasad et al. (2014. The correctness of the 'no model' approach for a flow configuration similar to the one studied in the present work was verified in Wawrzak and Tyliszczak (2018), and here it will be only partially repeated for completeness of the presentation.

Numerical Method
We apply an in-house numerical code (SAILOR) for low Mach number flows with variable density and temperature. The solution algorithm is based on the projection method for pressure-velocity coupling (Tyliszczak 2014(Tyliszczak , 2016. The time integration is performed by means of an operator splitting approach where the transport in physical space and chemical terms are solved separately. The convective and diffusive parts of the Eqs. (2)-(4) are advanced in time using a predictor-corrector technique with the 2nd order Adams-Bashforth/Adams-Moulton methods. The chemical reactions are computed using CHEMKIN interpreter (http://www.react ionde sign.com/produ cts/open/chemk in.html) and the detailed mechanism of Mueller et al. (1999) involving 9 species and 21 reactions for the hydrogen combustion. The reaction terms are stiff and therefore they are integrated in time applying the VODPK (Brown and Hindmarsh 1989) solver that is well suited for stiff systems. The spatial discretization is performed on half-staggered meshes applying the 6th order compact difference approximation for the Navier-Stokes and continuity equations (Santhanam et al. 2003;Laizet and Lamballais 2009;Tyliszczak 2014). The 5th order WENO scheme (Shu 2003) is used for discretization of the convective terms in the equations for 1 3 the species mass fraction and enthalpy. The applied code has been used and thoroughly verified in previous studies devoted to gaseous flows and flames (Tyliszczak 2015a, b;Rosiak and Tyliszczak 2016;Wawrzak and Tyliszczak 2018).

Spark Model
The spark is modelled by the energy deposition model, where the source term Q in Eq. (4) denotes the energy added to the flow through the enthalpy equation. Such a spark does not changes the mixture composition directly. This happens only when chemical source terms activate due to increased temperature level. The source term Q has the Gaussian distribution in time and space and is defined as Lacaze et al. 2009): where r is the radial distance to the centre of the spark r 0 , t denotes the time and t 0 is the time when the spark energy reaches the maximum. Three control parameters are the energy deposited , s = s ∕a and t = t ∕a , where s and t are the characteristic size and duration of the spark, respectively. The coefficient a = 4 √ ln(10) ensures that 98% of the deposited energy is present in the domain 3 s × t (Lacaze et al. 2009). The spark defined by Eq. (7) has a spherical shape, which does not change in time. The flow surrounding the spark and being inside it does not affect its size. However, when the spark is switched on, its intensity (power) changes in space and time according to Eq. (7) and reaches the maximum at t = t 0 . As will be shown later, before this occurs, an initially spherical flame kernel arising around the spark is significantly torn by the neighbouring vortical structures.

Configuration and Simulation Details
A sketch of the computational geometry is shown in Fig. 1. It is a rectangular box with the dimensions L x = 0.025 m, L z = 0.01 m and L y = 0.04 m in the streamwise (x), spanwise (z) and transversal (y) directions, respectively. The flow is periodic in the x and z-directions while for y = ±L y ∕2 the moving walls are assumed. The fuel flows in the upper area and air flows in the lower one, whereas both streams have the same initial temperature of 300 K. The fuel is a mixture of hydrogen and nitrogen with the mass fractions Y H 2 = 0.1 and Y N 2 = 0.9 for which the stoichiometric value of the mixture fraction is ST = 0.225 . In this work, following the Bilger's definition, the mixture fraction is computed as: where Y's are the elemental mass fractions, W's are the atomic weights and the subscripts, f and o, correspond to the fuel and air, respectively.
The initial velocity field is defined by a hyperbolic tangent function as u(y) = U ∞ tanh(2y∕ ) where U ∞ is the free stream velocity and = 2U ∞ ∕|du∕dy| max is the initial vorticity thickness assumed equal to = 0.5 mm. The Mach number based on the free stream velocity is approximately equal to 0.05, which justifies the use of low Mach number approximation. The shape factor of the velocity profile H = * ∕ = 2.25 , where dy is the momentum thickness, indicates that the initial velocity profile is typical for a boundary layer in transitional regimes. The Reynolds number is Re = U ∞ ∕ air = 600 , which corresponds to jet flows at moderate speeds. The assumed tangent hyperbolic profile reflects very well the velocity profile, which in a jet forms 1-2 diameters downstream a nozzle exit. In the temporally evolving mixing layer the tangent hyperbolic profile was used by Moser and Rogers (1993), among others. In the region of the mixing layer ±10 we introduce an isotropic initial disturbance defined through the Passot-Pouquet kinetic energy spectrum, which in terms of a wave number, k, is given as (Passot and Pouquet 1987): where k 0 is an adjustable wave number allowing generation of a turbulent flow field with the required integral length scale ≈ k −1 0 and Taylor length scale = Here, the symbol u ′ is the initial velocity fluctuation. The Passot-Pouquet turbulent kinetic energy (TKE) spectrum is known to skip the inertial range and is widely used to represent the flow scales having dissipative character (Martinez et al. 1997). The energy contained in is approximately 20% larger than in the range k ∈ [0, k 0 ] and therefore the initial conditions largely depend on the small eddies of the Kolmogorov length scale . The Passot-Pouquet spectrum was used to generate the initial conditions proposed by Doom and Mahesh (2009) for auto-ignition studies in a shear-less mixing layer of hydrogen-air mixtures and by Kolla et al. (2016) in research focused on behaviour of a temporally evolving jet in varying Damköhler number regimes. The generator of the initial conditions requires a random set of numbers (seeds), which along with the assumed values of k 0 and u ′ , lead to the velocity fields that follow the same TKE but differ in the spatial distribution of the velocity fluctuations. In this study we consider 16 initial

Fig. 1 Schematic view of the computational configuration (not in scale)
velocity fields, S0-S15, whose common global parameters are listed in the Table 1. For an assumed turbulence intensity u � ∕U ∞ equal to 1%, the Reynolds number based on the Taylor length scale of the initial fluctuations is equal to Re = u � ∕ air = 1.2 and the Kolmogorov length scale is ≈ ∕10 . When the computations start and time passes, the flow develops and larger vortices arise. At the ignition moments, the values of Re vary from approximately 30-60, which is of the same order as those used in Hawkes et al. (2007), Doom and Mahesh (2009), Chakraborty and Cant (2011).

Verification of the Computational Settings
The computational mesh is uniform in the periodic directions and compacted towards the centre of the domain in the y-direction. The compactness of the grid nodes is defined by a hyperbolic tangent function in such a way that in the vicinity of the mixing layer ( −2 < y < 2 ) the nodes are distributed almost uniformly. Figure 2 shows the computational domain with sample results illustrating the size of the vortical structure (2D contours and 3D iso-surfaces of the vorticity magnitude) and the flame kernel (2D contours of temperature) with respect to the initial shear layer thickness. These results were obtained for the successful ignition event and correspond to the time moment when the spark has been just switched-off and the flame has begun to freely propagate in space. The mesh used in all the simulations consisted of 512 × 128 × 256 nodes in the x, z and y-direction, respectively. Compared to the Kolmogorov length scale the cell sizes were equal to: ∕ x = 0.98 , ∕ z = 0.61 , ∕ y min = 1.72 (at centreline). The mesh resolution slightly changes in the region of the spark placement ( −4 ≤ y ≤ 4 ) and at the locations at y = ±4 we have ∕ y = 0.96 , which from the point of view of the LES ensures very good representation of the smallest resolved scales. Note that in the Vreman's sub-grid model (Vreman 2004) the sub-grid viscosity vanishes when a mesh meets DNS requirements. The mesh lines shown in Fig. 2 are plotted every four mesh points in all directions. However, it can be seen that the smallest sub-filter curvatures of the vorticity or temperature isolines are smooth and comparable with the cell sizes. A grid independence study for a flow configuration similar to the one analysed in the present work were performed in Wawrzak and Tyliszczak (2018). There it was shown that the grid with 512 × 128 × 256 nodes is adequate to well predict the phenomena crucial to ignition, i.e. vortices formation and their evolution at the later time. A convincing argument supporting this statement is provided in Fig. 3 showing the ratio of the sub-filter velocity to the free stream velocity and the ratio of the turbulent viscosity to the molecular one. The presented results have been extracted along the x-direction in the central cross-section plane and are conditioned by the stoichiometric mixture fraction, i.e. only these values of u SGS and T for which = ST are presented. The sub-filter scale velocity fluctuation was computed based on the sub-grid viscosity as is the local LES filter width and the model constant is C = 0.1 (Pope 2000). The results presented in Fig. 3 correspond to two specific time instants, t = 1.00 ms and t = 1.25 ms, which are the times of formation of large scale vortices and their roll-up process, respectively. At these specific time moments the spark energies will reach maximum values. At t = 1.00 ms and t = 1.25 ms the Reynolds numbers computed based on the Taylor length scale in the homogeneous z-direction Table 1 List of parameters of the initial velocity fields S0-S15 S0-S15 0.5 600 0.01 1.2 /10 in the region y∕ = ±10 are equal to Re z = 29 and Re z = 58 . This means that the flow significantly evolved and the vortices grew noticeably compared to their initial sizes. Nevertheless, all flow structures are accurately represented by the computational mesh. This can be deduced from the fact that the sub-filter velocity fluctuations are very small and the turbulent viscosity is almost two orders of magnitude smaller than the molecular viscosity, as can be seen in Fig. 3. In the whole domain and during the entire simulation time, the ratio of the subfilter kinetic energy estimated as k SGS = 2 T ( C ) 2 to the total turbulent kinetic energy is less than 0.1 in all computations performed in this paper. Hence, according to Pope's criterion (Pope 2004) the used mesh provides proper resolution from the point of view of LES needs. It holds not only taking into account the initially assumed flow field but also during its later evolution. Nevertheless, to show the level of dependence of the results on the mesh density, the solutions obtained using two additional meshes consisting of different number of grid points in the ' x − y ' plane, i.e. 256 × 192 and 512 × 384 were performed. Concerning the spanwise (z) direction the use of the mesh with 192 nodes virtually does not change the solution. The results obtained from comparative simulations will be presented in Sects. 3.2 and 4.
As discussed above, the computational mesh applied is sufficiently dense for accurate and well resolved LES studies. However, its correctness for the 'no model' approach requires a bit more attention. Figure 4 shows a regime diagram for LES of turbulent combustion modelling proposed by Pitsch and Duchamp De Lageneste (2002). In this plot, particular regions are defined with respect to a local sub-grid Karlovitz number ( Ka SGS ) and the ratio of the LES filter width and the laminar flame thickness ( l ). The subgrid Karlovitz number combines the sub-grid quantities with the chemical scales, and hence, it may serve to verify whether the 'no model' approach method is valid at given mesh resolution. Following Pitsch and Duchamp De Lageneste (2002), the Karlovitz number can be estimated based on the laminar flame thickness, burning velocity ( S l ) and sub-grid fluctuations as: The flame thickness is roughly estimated as l = D∕S l = 2.01 × 10 −5 m, where D = 2.512 × 10 −5 m 2 ∕s is the heat diffusion coefficient for the stoichiometric mixture and the burning velocity of diluted hydrogen equals to S l = 1.3 m∕s (Li et al. 2018). Taking into account that x = 4.88 × 10 −5 m, z = 7.81 × 10 −5 m, y min = 2.79 × 10 −5 m and the (10) Fig. 4 Combustion regime diagram for LES according to Pitsch and Duchamp De Lageneste (2002) minimum filter width is min = 4.73 × 10 −5 m, it is clear that the flame front cannot be correctly resolved on the present mesh, even assuming that the above definition of l is too rigorous and the real flame thickness can be even a few times larger (Poinsot and Veynante 2001). The results shown in Fig. 4 correspond to the ones presented previously in Fig. 3. It can be seen that all points lie in the laminar flamelet regimes, where the turbulent velocity fluctuations are small and accurately captured by the computational mesh ( ≈ ), and where the flame thickness is smaller than the LES filter width ( > l ). According to Duwig et al. (2011) the 'no model' approach works correctly when the flow inside the LES cells can be treated as in perfectly stirred reactors, which means that it has to be homogeneous and the fluctuations must be small to have no significant impact on chemical reactions. A simple method allowing the assessment of the mixedness level of the reactants seems to be offered by using the mixture fraction gradients ( |∇ | ), whose values close to zero indicate the regions of homogeneous mixture or the regions of a pure fuel or oxidiser. Figure 5 presents the contours of |∇ |∕|∇ | MAX in the central cross-section plane at the time moments t = 1.00 ms and t = 1.25 ms. It can be seen that large values occur only in the central part of the domain, where the fuel and oxidiser streams mix and the vortical structures develop. The enlarged part of the figures show the places in which |∇ |∕|∇ | MAX changes from zero to one across a thin layer. One can see that this layer is covered by at least 5-6 grid cells, which is certainly sufficient to correctly represent its spatial evolution, especially with the use of high-order discretization method. However, the assumption that within these particular cells the mixture is homogeneous, as would be advisable for the 'no model' approach, is probably slightly exaggerated. Nevertheless, possible modelling errors have the impact only on a very small region of the flow, not necessarily corresponding to the spark location. Thus, neglecting the sub-filter scales in the reaction terms for the simulation of combustion in the present flow configuration should not falsify the underlying physics.

Spark Parameters
We assume the spark with an overall duration time t = 0.5 ms, total energy = 1 mJ and diameter s = 1.5 mm = 3 . These parameters are comparable with those used in Ahmed and Mastorakos (2006) for the ignition of a fuel jet and in our previous study (Wawrzak and Tyliszczak 2018) focussed on the importance of the spark parameters. The sparks are placed at x∕ = 10 , z∕ = L z ∕2 = 10 and at various vertical locations: y∕ = 0, ±1, ±2, ±4 . Note that the spark size is three times larger than that could suggest that its relatively small movement will not result in a significantly different ignition scenario. However, as will be presented later, this supposition turned out to be false. Figure 6 shows the spark locations (the red stars) along with the vorticity magnitude contours and isoline of the stoichiometric mixture (black line) at the selected time instants.
Characteristic phenomena for the temporally evolving mixing layer include continuous growth of the initial disturbances, formation of the vortical structures and their roll-up. Shortly after the simulation starts, at t = 0.5 ms, the fuel and oxidiser streams are still clearly separated and a boundary between them is only slightly waved. The first vortices, which are the results of growth of the most amplified disturbance mode, start to be seen at t = 0.75 ms (Fig. 6b). Evolution of particular disturbance modes can be analysed using a globally defined integral parameter (Lesieur et al. 1988;Moser and Rogers 1993): where the variable û x ( ) is the amplitude of the Fourier mode of the x-velocity component and û * x ( ) is its conjugate value. Figure 7 shows the time evolution of A of the most amplified mode ( A 1 ) and its subharmonics ( A 1∕2 ). The maximum of A 1 occurs at Fig. 6 Contours of instantaneous vorticity magnitude and ST isoline at various time instants prior to the spark initiation. The results are shown for the case S0 t A 1 MAX ≈ 1.00 ms and it corresponds to the time instant when the small vortices, visible along the shear layer at t = 0.75 ms, have fully developed and rolled-up. As can be seen in Fig. 6, at the later time they grow and roll-up when A 1∕2 attains the maximum value at t A 1∕2 MAX ≈ 1.25 ms. This scenario agrees very well with the findings of Moser and Rogers (1993). The results presented in Fig. 7 were obtained on three different meshes labelled as M1, M2 and M3. The differences between the solutions obtained on two densest meshes, M2 and M3, are negligible during almost the entire simulation time and are noticeable only from t = 1.75 ms. As will be discussed later, at this time instance, the spark will be already switched-off and the result of the ignition event (successful or failed) and a state of the flame (a size and direction of propagation) will be fully known. Hence, the mesh M2, which is used for all performed simulations is sufficiently dense to correctly predict the mechanism of flame kernel evolution. Concerning the mesh M1, which is four times coarser than M3, it can be seen that its use also ensures a qualitatively good solution. The observed differences are only quantitative. The time at which A 1 and A 1∕2 achieve maxima is predicted very well on all meshes and only their values on the mesh M1 are slightly under-predicted. This, however, can be expected as the content of the resolved energy on the coarser mesh is smaller by definition. Intuition suggests that for two different spark initiation times, the sparks will be surrounded by vortical structures of different sizes and this will lead to different behaviour of the flame kernel, i.e. to the success or failure of the ignition. Hence, we divided the test cases into two groups depending on spark initiation time. First we considered the situation when the spark was turned on at t sp = 1.0 ms ≈ t A 1 MAX and reached the maximum when the small vortices rolled-up and larger vortical structures were present in the flow at t 0 = t A 1∕2 MAX (see Fig. 6c, d). In the second case, the spark was turned on at the earlier time, t sp = 0.75 ms, when the mixing layer was only wrinkled (see Fig. 6b). Most of the attention was devoted to the first situation, while the second one was considered only for the ignition probability analysis. To study this issue we generated sixteen initial flow fields (S0-S15) having the same global characteristics ( Re , Re , T i , ) but resulting from a different random seed used to compute the velocity fluctuations, as discussed in Sect. 3.

Flow Characteristics Prior to Spark Ignition
The ignition process depends on the mixture composition, its gradients and the velocity field. In the following sections we focus on two cases corresponding to the initial conditions labelled as S1 and S8, which were selected from all simulations as the representative for highlighting differences in the ignition process. The flow conditions, i.e. the mixture fraction ( ), the scalar dissipation rate (SDR, = 2D∇ ⋅ ∇ ) and module of the strain rate tensor ( �S� = √ 2S ij S ij ) representing the mixture fraction and velocity gradients, are presented in Table 2. It shows the values for the cases S1 and S8 and also the mean values calculated based on all 16 cases at all spark locations ( x∕ = z∕ = 10 , −4 ≤ y∕ ≤ 4 ) at the time t sp = 1.0 ms. The underlined values correspond to the cases which are analysed in detail in the rest of the paper. In Table 2 the mixture fraction is normalized by its stoichiometric value and the strain rate tensor by its maximum value at the beginning of the simulation ( |S| INI MAX ). As can be seen the instantaneous vorticity thickness ( (t) ) is approximately three times larger then the initial one. In the case S1 the values of |S| and (t) are slightly higher than for the case S8, however, they are not far from the mean values. Significant changes are found when the spark location changes, especially in the region −2 < y∕ < 2 . There, the strain rate and the SDR for the case S1 are higher than in the case S8. Note that mean values of |S| are almost symmetric across the mixing layer while the SDRs are not. This is because in some realizations the spark is placed almost exactly on the interface between the fuel and the oxidizer stream, where the mixture fraction gradients are very large. Such events increase the mean value to the level that cannot be compensated by other events due to small number of samples. Worth noticing is that most of the sparks is placed in the fuel rich locations, but it should be reminded that the local conditions in the spark region changes dynamically as the vortical structures continuously evolve in the flow (see Fig. 8). Hence, pockets of rich mixture are not constantly present at the ignition location. Table 2 List of conditions prior to the ignition at the spark centre ( x∕ = z∕ = 10 , y∕ ) for the initial velocity field S1, S8 and the mean values from the cases S0-S15 The underline indicates the cases, which are in detail analysed in the paper Analysis of the time sequences of the obtained results shows that in the region where the strain rate is high ( −2 < y∕ < 2 ) the sparks "experience" significantly different mixture compositions.

Spark Ignition Scenario
Once the spark is turned on, a spherical region of increased energy appears in the flow domain. This is shown in Fig. 8 presenting the location and size of the spark (black circle), vorticity magnitude contours and temperature iso-surfaces during the ignition process at four selected time instants. Neither the shape of the spark nor its intensity, represented by varying colour inside the black circle, is affected by the flow field in the neighbouring vortices. It can be seen that along with the increasing spark energy, the region in its vicinity becomes hotter and high temperature values appear. Unlike the spark, the temperature iso-surfaces, which initially also had spherical shape, are quickly deformed by the surrounding vortices and it seems as if they were washed away from the spark surface. The temperature rises quickly and at t = 1.25 ms the flame engulfs the whole spark region and spreads in its vicinity, see Fig. 8c. It can be seen that the temperature iso-surfaces are much 'smoother' than the vorticity contours. Furthermore, the vorticity inside the hot region has significantly smaller values than on the boundary of the flame. This can be attributed to the fact that for a higher temperature, the molecular viscosity increases, and thus, the local Reynolds number decreases and turns out to have a stabilising effect. If the chemical reactions are activated before the spark is switched-off, the flame kernel becomes self-propagating by spreading depending on the local mixture composition and being torn and stretched by the neighbouring vortices. Similar situation was considered in previous study (Wawrzak and Tyliszczak 2019) as successful ignition scenario No. (1). In the present paper we will stay with this notation. When the spark is placed outside the flammable region ( P f = 0 ), the heat produced by the spark decays without initiating the chemical reactions. It may happen that even if the spark overlaps with flammable mixture regions, the increase of the temperature is insufficient to initiate the reactions. Moreover, one should note that the failed ignition can be also observed when the spark is placed close or inside the region overlapping the stoichiometric mixture. Such a location of the spark can be only seemingly favourable for ignition as the presence of strong shear stresses or vortical structure may either destroy the flame kernel or move it to inflammable very rich/ lean flow regions. According to Wawrzak and Tyliszczak (2019), such a failed ignition scenario is termed as No. (2). However, the opposite situation is also possible, i.e. the hot fluid created in such regions can be brought by the vortices to the region with favourable conditions where the flame can develop. This is the second successful ignition scenario (No. (3)) that leads to a stable flame (Wawrzak and Tyliszczak 2019).

Ignition Kernel Evolution
To further investigate two different successful ignition scenarios, the flame kernel behaviour after the spark initiation had to be analysed with scrutiny. Its spatio-temporal evolution is illustrated in Figs. 9 and 10 using the high temperature isosurface coloured by the heat release for the cases with spark ignition scenarios No. (1) and No. (3), respectively. Based on all ignition attempts we found that the ignition scenario No. (1) prevailed on the rich side of the mixing layer (100% for y∕ = 1 , y∕ = 2 ; 87.5% at y∕ = 0 ) and the scenario No.
(3) occurred significantly more often on the lean side (83.3% at y∕ = −2 ). The results presented in Figs. 9 and 10 correspond to the initial conditions labelled as S1 and S8, as discussed in Sect. 3.2. The temperature value of 1200 K is 'a priori' chosen to represent the shape and volume of the flames. The isolines of the vorticity magnitude at central ' x − y ' and ' z − y ' cross-sections are also shown to visualise the turbulent structures in the neighbourhood of the spark location (black circle). It can be seen that at the time instant t = t 0 , at which the sparks reach the maximum energy, the turbulent structures already exist in the flow. Although it seems that in the cases S1 and S8, the conditions close to the spark are quite similar (see Table 2), the flame evolutions at the later stages are substantially different. In the case S1, the flame develops near its initial position and virtually does not move. In the case S8, the situation is different. In the first phase of the ignition, called later the sub-stage I, the flame kernel approaches the favourable local mixture conditions where its further expansion occurs. This process is accompanied by a rolling up vortical structure and its movement towards the left boundary. From now on, the flame development and propagation is dominant (sub-stage II). Well mixed fuel-air pockets ignite along the vortex chain which give the effect of flame stretching. This causes that the flame deforms and its shape becomes highly irregular and mostly stretched in the streamwise direction and spanwise direction. The colours representing the heat release shows in which direction the flame propagates. The temperature distribution shows that when the spark duration ends, the flame S1 extends over 15 in x-direction whereas the S8 spreads over the distance smaller than 10 . In contrast, the flame S8 is wider in the z-direction compared to the case S1. It seems obvious that the size of the flame in case S8 is smaller, as it was ignited later, however, it will be shown momentarily that the speeds of the flame volume growth are also different in the cases S1 and S8.
To give more light on each ignition scenario Fig. 11 shows the transient evolution of the maximum heat release ( Q MAX ) found at T = 1200 K iso-surface and the mixture fraction corresponding to these specific locations for ignition events S1 and S8. It is seen that the ignition occurs earlier for the case S1, as the sudden increase in Q MAX is observed around t = 1.27 ms whereas for the case S8 the Q MAX occurs at t = 1.31 ms. At the beginning, when the kernel expansion overwhelms the convective motion, the fluctuations of are large compared to the later phase of the flame development. Nevertheless, the locations of the Q MAX suggest that the flames propagate along the lean and nearly stoichiometric mixture for both ignition scenarios.
Two sub-stages identified for the ignition scenario No.
(3) are further analysed in Fig. 12 showing the time evolution of the maximum temperature ( T MAX ) inside the computational domain, the normalized strain rate ( S = |S|∕|S| INI MAX ) at the locations of the mesh points where T MAX is found, and the rate of change of a flame volume ( dV FLAME ∕dt ). The flame volume is defined as the sum of the volumes of computational cells in which the temperature is larger than 1200 K ( V FLAME = ∑ V T≥1200 K cell ). Figure 12 also shows the solutions obtained applying the mesh M3 (case S8-M3). As can be seen, the profiles almost overlap, which to some extent justifies the use of 'no model' approach as the sub-filter activity in Fig. 8 Flame kernel evolution when the spark is turned on at y∕ = 1 . The black circle shows the spark location and the colour inside the circle corresponds to the spark intensity Q . The contours show the vorticity magnitude in the ' x − y ' and ' y − z ' cross-section planes crossing the spark. The iso-surfaces show the temperature distribution in the vicinity of the spark. The values of temperature iso-surfaces are given in the subfigure captions 1 3 the vicinity of very intense reactions is at very low level. The T MAX profile shown for the case S1 (scenario No. (1)) characterizes a steep increase during the sub-stage I, followed by a small temperature drop. In the sub-stage II the temperature stabilizes and the flame volumes grows rapidly. A significantly different behaviour is reflected through the two-steps rise of T MAX observed for the case S8 (scenario No. (3)). Here, the temperature initially increases to a moderate level (sub-stage I), it remains there for a short time or slightly decreases, and then suddenly rises (sub-stage II) and stabilizes similarly as in the case S1. Such a two-step flame evolution is attributed to unfavorable flow conditions in the place/ time of the spark generation (sub-stage I) and the movement of hot pockets of gas towards a region where the chemical reactions can proceed (sub-stage II). Comparison of the S evolution (Fig. 12b) shows that delay of the ignition in scenario No. (3) in the case S8, is caused by large velocity gradients in the region of the flame kernel formation. Although at the Fig. 9 Instantaneous temperature iso-surfaces (1200 K) coloured by heat release and vorticity magnitude isolines in the central cross sections at subsequent time instants following the successful ignition scenario with the spark located at y∕ = −1 . Results for the case S1 Fig. 10 Instantaneous temperature iso-surfaces (1200 K) coloured by heat release and vorticity magnitude isolines in the central cross sections at subsequent time instants following the successful ignition scenario with the spark located at y∕ = −1 . Results for the case S8 initial moment of the spark generation the strain rate for the case S8 is approximately two times smaller than for S1 (see Table 2), later, in the sub-stage I it significantly grows and its maximum is 1.3 times larger than the maximum of S observed for the case S1. Large velocity gradients that contribute to S are partially caused by density differences due to the spark generation and partially by the movement of strong vortical structure close to it (see Fig. 10). Regarding the difference in the flame volumes in the cases S1 and S8, the V FLAME is approximately two times smaller for S8 during the entire simulation time. Figure 12c shows the speed of the flame volume growth dV FLAME ∕dt during the ignition and at the later phase. The red lines represent the linear fit of the results and they indicate that before the T MAX stabilizes the flame in the case S1 expands slightly faster, and on the other hand, when the T MAX becomes equal for both the cases then the flame growth is more intense for the case S8, though the difference is not very large. Worth noticing is that when the T MAX is stable the dV FLAME ∕dt is approximately two times bigger than during the transient phase.
In the following the results obtained from the ignition trials for different initial solutions (S1, S4, S13 and S14) at the spark locations y∕ = 0 and y∕ = −2 are compared. Figure 13 shows the ratios T (1) MAX ∕T (3) MAX and V (1) FLAME ∕V (3) FLAME , where the superscript denotes the ignition scenario. The results at three time instances are presented, t = 1.35 ms, t = 1.50 ms and t = 1.65 ms. They are denoted by right triangle, left filled triangle and diamond, respectively. The first time moment corresponds to the time shortly after the spark reached the maximum energy, the second one matches the end of the spark duration, and the last one is the time when the maximum temperatures stabilize. Initially, at t = 1.35 ms the T MAX and V FLAME ratios are large, which means that the flame volumes differ substantially. At the later time, as the temperature maxima stabilize and the ratios T (1) MAX ∕T (3) MAX approach to one, the ratios V (1) FLAME ∕V (3) FLAME decrease. Nevertheless, in all analysed cases, the flame volumes in the ignition events following the scenario No. (1) were larger and the ratios V (1) FLAME ∕V (3) FLAME were in the range 2-4.

Flame in Mixture Fraction Space
Structure of the flame in mixture fraction space is presented in Fig. 14. The scatter plots show the temperature coloured by the mass fraction of OH radical for three cases, which differ by the spark location ( y∕ = 0 , y∕ = −1 ) and randomness of the initial flow field (S0, S1). It can be seen that both these factors have a strong impact on the solution. The ignition was successful for both the cases but the initial flame development scenarios were significantly different. This can be concluded from the fact that the points with high temperature values are located in different parts of -space. It can be seen that at t = t 0 the flame begins to develop only for the case S1 with the spark at y∕ = 0 . The increased OH mass fraction values indicate the proceeding chemical reactions. The most interesting seems to be the case S0, with the spark located at y∕ = 0 , for which at t = t 0 the majority of the 'hot points' lies in a very rich region, even close to = 1 . As can be seen in Fig. 14b when the time passes they all move towards the lean fuel side where they fit under an envelope typical for the flamelet type solution. A sample flamelet profile for the case with Y H 2 = 0.1 , Y N 2 = 0.9 , T = 300 K for = 1.0 and Y O 2 = 0.1 , Y N 2 = 0.9 , T = 1200 K for = 0.0 is represented in Fig. 14b by the down oriented black triangles.
Concerning the results at particular time moments, one should realize that the points with high temperature on the rich side at t = t 0 (open circles) cannot be related to the heat release and chemical reactions as the flame cannot exist above the flammability limits. The OH concentration is at a very low level that also indicates that the flame does not 14 Scatter plots of the temperature coloured by the OH mass fraction for the cases S0 and S1 at t = t 0 (a) and for the case S0 for three time instants (b) evolve in this range of . Hence, it seems that this situation results from diffusive heat transfer and movement of hot gases from the burning region towards the upper fuel-rich side of the domain. The latter effect is dominant and is caused by strong mixing, which is manifested by large values of the scalar dissipation rate (SDR, = 2D∇ ⋅ ∇ ). Its conditional values ( | T>1200 K ) at the selected time instants are shown in Fig. 15. Additionally, the conditional mass fraction of the OH and heat release values are presented in Figs. 16 and 17. It can be seen that at t = t 0 , none of these two quantities exist on the rich side, whereas the SDR reaches high levels mainly for > 0.5 . This unequivocally means the presence of intense mixing (large mixture fraction gradients → large ) but no chemical reactions. At the later times, when the flame has partially developed, the maximum values of the conditional SDR concentrate only in the vicinity of ST = 0.225 and they decrease with time. Based on the scatter plots showing OH radicals and heat release, one can deduce that after the spark is switched-off, the flame expands in physical space mainly along the lean and near stoichiometric regions, which is consistent with results presented in Fig. 11 and findings of Chakraborty et al. (2007).

Spark Ignition Probability
Finally, based on 16 different initial fields we focus on the probability of the flame kernel generation ( P ker ) and ignition ( P ign ) across the mixing layer. For each of 224 ignition attempts (16 trials × 7 locations × 2 ignition times) the maximum temperature and the mass fraction of OH species were monitored. The flame kernel formation was classified as successful if the temperature increased high enough to initiate the chemical reactions and the Y OH raised over 2 × 10 −4 . This condition was chosen in analogy to the one used in monitoring the occurrence of auto-ignition spots. If after the spark was switched-off the high temperature persisted and the flame volume continuously grew such an event was treated as successful ignition. In contrast, if the temperature and OH content dropped down shortly after the spark deposition, the ignition was treated as failed. There were also ignition events in which the flame kernel has not been generated at all, i.e. the OH was not produced. In the cases in which the above depicted signs of ignition were observed the simulations were continued for the time period 3 t , otherwise they lasted only for 1.5 t from the spark initiation. We compare the results analysed previously for the ignitions at the time t 0 = t A 1∕2 MAX and also the solutions for the cases with the sparks turned-on at the earlier time with t 0 = t A 1 MAX , i.e. before the small vortices rolled up and transformed to the large structures.  Figure 18 shows the profiles of mean mixture fraction ( (y) MEAN ) variations across the mixing layer, which were calculated by spatial averaging in the ' x − z ' planes. The length of the horizontal solid line corresponds to the spark size s = 3 . It shows that the regions of two neighbouring sparks overlap. Although the initial fields are characterised by the same global parameters ( Re , Re , T i , , mean velocities), it can be seen that the profiles of (y) MEAN are very similar but not exactly the same and differences increase with the time. They would be identical only in the case of infinite domain in the x and z-directions. The fact that depending on the initial velocity disturbances the evolved vortical structures had slightly different statistical properties allowed us to presume that at the ignition times the flow fields in the vicinities of the sparks differed from each other. It was verified by careful checking of every single flow field. Hence, in each case, the sparks experienced different flow conditions, even if they were placed in the same spatio-temporal locations. This led to failed or successful ignition with different flame development phases, as was discussed in Sect. 4.2. Note that the relatively small number of trials slightly diminishes credibility of the probability analysis, however, it gives an idea on how sensitive the ignition event is on the location of the spark close to the mixing layer. The monitored ignition event is binary ('no', 'yes'), so its variance is given by c(1 − c) , where c is the number of 'yes' divided by the total number of tests N. The statistical error is given by √ (c(1 − c))∕ √ (N) , so, the statistical uncertainty is the highest for P ign = 0.5 and equals to 0.125. It is smaller at lower and high probabilities.
The computed probabilities are presented in Fig. 19 where the bars denote the statistical error at particular probability values. It can be seen that P ker as well as P ign distributions manifest the same trend, they are not symmetric in y-direction and vary from one point to another even if the spark regions overlap. Comparing directly their values it can be seen that there are only a few points in which P ker < P ign , which means that most of the flame kernels successfully evolved in time and space. Substantial sensitivity of the P ker and P ign on the spark localization is in line with Mastorakos (2009) and confirms that this process is very much dependent on local flow regimes. The highest P ign occurs in the centre of the domain where the mean velocity is close to zero and it decreases drastically when the spark is moved outward. It falls to zero when sparks are placed at y∕ ≥ ±4 , where the mean mixture compositions are beyond the flammability limits. This is consistent with the findings of Ahmed and Mastorakos (2006) for a jet type flame. They observed that in a region close to a nozzle the flow can be ignited only if the spark is placed in a thin region of the mixing layer in which a fuel is mixed with the oxidiser. In the present study, we found that ignition probability increases for the cases with the spark initiated later ( t 0 = t A 1∕2 MAX ) as in this time moment the fuel and oxidiser are better mixed and the shear stresses existing in (a) (b) Fig. 18 Mean mixture fraction profiles across the mixing layer for the cases S0-S15 (a) (b) Fig. 19 Probability of flame kernel formation (a) and ignition (b) across the mixing layer 1 3 the initial flow field become weaker. Knowing that larger vortices are present at the later time it seems that their presence does not have a significant impact on the P ign distribution as it shows the same trends for t 0 = t A 1∕2 MAX and t 0 = t A 1 MAX . Furthermore, we observed that: 1. When the ignition occurs in the region ( −y∕ , y∕ ), the kernels formed at leaner mixtures are less resistant than these at the fuel-rich side. This is somewhat surprising as the mean strain rates are similar on both sides (see Table 2) while the mean mixture composition in (−y∕ , 0) is closer to the stoichiometric regimes-see Fig. 18 and Table 2. 2. When the ignition occurs at y∕ ± 2 then the P ign drastically drops down on the rich side, whereas for the sparks located on the leaner side, the P ign stays at a comparable level as in the centre of the domain. 3. Regarding the ignition probability in view of the mean values presented in Table 2 it can be said that the strain rate and the scalar dissipation rate play less important role than the mixture composition. The P ker and P ign are the definitely higher on the fuel rich side (2.5 ST < < 3.3 ST ).

Conclusions
The paper presented LES studies of spark ignition process in the temporally evolving mixing layer formed between the hydrogen/nitrogen mixture and air. The performed test computations and analyses of the sub-filter scale quantities showed that the applied mesh was sufficiently dense to apply LES with 'no model' approach for the reactive terms. The aim of the study was to analyse the flame kernel behaviour and its topology during the spark initiation and soon after it is switched-off. We focused on the flames following two different ignition scenarios (No. (1) and No. (2)) that depend on the localization of the spark with respect to the vortex structures rolling-up in the mixing layer. It was shown that some sparks led to a sudden and high temperature rise resulting in a quick development of the flame whose initial shape was torn and stretched by the vortical structures. In such a case, the flame kernel was formed near the spark position. Different ignition scenario was observed mainly for the spark locations shifted from an easily ignitable region. If this was the case, the flame kernel formation was significantly delayed and two sub-stages of the ignition process were identified: (1) the temperature raised only to some intermediate level at which the chemical reactions could not proceed; (2) the heat from the volume of the hot gas diffused and in the same time this volume was torn by the vortices towards the flow parts with ignitable mixture composition where the flame could rapidly develop. In such a situation, the flame shape was very irregular and strongly dependent on the flow conditions around it. It was found that in the cases when the developing flames were moved from the place of initiation, their volumes were significantly smaller compared to the volume of the flames growing near the spark localization. It was found that the scenario No.
(1) prevailed when the spark were placed if the fuel rich locations while the scenario No.
(3) dominated on the lean side. In all successful ignition events the flames expanded mainly along the lean and near stoichiometric regions. The analysis of the solutions in mixture fraction space showed that this happens even if the spark was placed in the fuel rich side. In such events the nascent flame first propagated quickly towards the leaner conditions and then developed.
Concerning the ignition probability, the present study confirmed the known fact that the ignition is hardly possible when the spark is placed far from the mixing layer region. The unsuccessful ignition trials were observed also for the sparks placed in the regions of ignitable mixture but in unfavourable flow conditions where surrounding vortices either destroyed the nascent flame or moved it towards non-flammable flow parts. It was found that when the spark was placed in the mixing layer and close to the initial boundary between the fuel and oxidiser, the ignition probability was asymmetric and larger on fuel richer side. The presence of larger vortices seems to have a large impact on the flame evolution but not on the ignition probability. Unlike the local composition the impact of the strain rate of the velocity field was also very limited in view of the ignition success. This aspect certainly needs more research with varying spark parameters (size, duration, energy) and substantially larger number of samples, in particular, as only then the statistical analysis can be fully credible. Nevertheless, taking into account that the sparks used in the present study were three times larger than the initial mixing layer thickness and with the energy significantly exceeding the minimum ignition energy, the presented distribution of ignition probability showed how sensitive the ignition process is to the location of the spark and local instantaneous flow regimes.