Experimental assessment of presumed filtered density function models

Measured filtered density functions (FDF) as well as assumed beta distribution model of mixture fraction and ’subgrid’ scale (SGS) scalar variance z2, used typically in Large Eddy Simulations, were studied by analysing experimental data, obtained from two-dimensional planar, laser induced fluorescence measurements in isothermal swirling turbulent flows at a constant Reynolds number of 29000 for different swirl numbers (0.3, 0.58 and 1.07). Two-dimensional spatial filtering, by using a box filter, was performed in order to obtain the filtered variables, namely resolved mean and ’subgrid’ scale scalar variance. These were used as inputs for assumed beta distribution of mixture fraction and top-hat FDF shape estimates. The presumed beta distribution model, top-hat FDF and the measured filtered density functions were used to integrate a laminar flamelet solution in order to calculate the corresponding resolved temperature. The experimentally measured FDFs varied with the flow swirl number and both axial and radial positions in the flow. The FDFs were unimodal at flow regions with low SGS scalar variance, z2 <0.01, and bimodal at regions with high SGS variance, z2 >0.02. Bimodal FDF could be observed for a filter size of approximately 1.5-2 times the Batchelor scale. Unimodal FDF could be observed for a filter size as large as four times the Batchelor scale under well-mixed conditions. In addition, two common computational models (a gradient assumption and a scale similarity model) for the SGS scalar variance were used with the aim to evaluate their validity through comparison with the experimental data. It was found that the gradient assumption model performed generally better than the scale similarity one.

In addition, two common computational models (a gradient assumption and a scale similarity model) for the SGS scalar variance were used with the aim to evaluate their validity through comparison with the experimental data. It was found that the gradient assumption model performed generally better than the scale similarity one.

I. INTRODUCTION
In Large Eddy Simulations (LES) of turbulent reacting flows, spatially-averaged versions of the equations describing fluid motion are usually solved on a grid, which is 'coarse' relative to the smallest scales of fluid motion. In LES, the dynamics of scales that are larger than the filter width (the resolved scales) are captured explicitly while the effects of the smallest or 'subgrid' scales that are not resolved by the filter need to be modelled. The filtered quantity f (x, t)(velocity, mixture fraction etc.) in physical space is simply defined as a convolution of the non-filtered field f (x , t) with a function or convolution kernel, G ∆ (x − x ) as follows 1 : The low-pass spatial filter function G ∆ has characteristic filter width ∆, satisfying G ∆ (x − x )dx = 1 and is independent of spatial location. The instantaneous value of any physical quantity f in a turbulent flow can be split into a spatially filtered component (the resolved field) and the fluctuations around the filtered component (the residual field) 2 f =f + f . The spatially filtered component refers to the resolved mixture fraction fieldz, which is computed as a convolution (moving average) of measured mixture fraction fields z with a box filter of various sizes. The effect of the filter size on the resolved mean will be discussed in section III.
A laminar flamelet approach can be used in LES to model reacting flows due to its simplicity and flexibility. This approach is based on a concept of thin, laminar, locally onedimensional structures (called laminar flamelets) embedded into the turbulent flow field.
The main idea behind this approach is to separate the detailed and complex chemistry calculations and the flow fluid dynamics. The laminar flamelet approach reduces the complex chemistry to only two variables, which are the mixture fraction and the scalar dissipation rate. All the dependent variables can then be linked to a single parameter, e.g. the mixture fraction, or to both the mixture fraction and the scalar dissipation rate. The dependent variable is the variable that is computed as a function of the mixture fraction and the scalar dissipation rate.
In order to embed locally one-dimensional laminar flamelets into the turbulent flow, a concept of probability density function in the context of Reynolds-Averaged Navier Stokes However, the most simplified and practical modelling approach is based on presumed filtered density functions (FDFs), which are usually approximated by beta distribution model or top-hat FDF, parameterized by the first two statistical moments of the mixture fraction, namely the filtered mixture fraction (the resolved mean) and the 'subgrid' scale mixture fraction variance. The validity and applicability of the β−function approximation were investigated by using DNS data of non-premixed reacting flows 7 . The results confirm that the β−function provides a good estimate for the FDF of the mixture fraction. However, it was shown 8,9 that the FDF could substantially deviate from assumed beta distribution model Parameters a and b of the beta distribution model are related to the resolved mixture fraction and the 'subgrid' scale mixture fraction variance as follows: Parameters z a and z b of the top-hat function are also related to the resolved mixture fraction and its 'subgrid' scale scalar variance as follows: where z 2 is resolved 'subgrid' scale scalar variance,z is resolved mean computed from measured mixture fraction field The filtered mixture fraction is determined by the solution of the transport equations and its variance is usually modelled 11,12 . If the FDFs are known, the resolved quantities from the laminar flamelet solutions can be computed directly; for instance, if chemistry is assumed to be infinitely fast, local equilibrium is reached and species mass fractions, reaction rates or red temperature are functions of mixture fraction z only. The mixture fraction z is defined as the ratio of mass of fluid having its origin in the fuel stream to the mass of the stream, in two stream problems. The mixture fraction varies between zero and one, corresponding to no fuel (0) and pure fuel (1).
where Y f uel is fuel mass fraction, Y ox is oxidizer mass fractions, M W is molar mass, v is number of moles of oxidizer needed to completely oxidise one mole of fuel.
The spatially averaged quantities are directly obtained by an integration with given filtered density function as follows: where φ(x, t) denotes any spatially filtered variable and φ(z) is any variable from the laminar flamelet solution (e.g. temperature, species mass fractions etc.).
The filtered density function f of a scalar φ is written as 8 .
whereφ is the sample-space variable for φ, δ φ (x , t) −φ is the Dirac delta function The filtered density function is in fact a probability density function of a subfilter state, which characterizes the distribution of the scalar within the confines of the filter. The FDF can be constructed from experimentally measured mixture fraction fields or any measured scalars as follows: at a given spatial location, a box filter of known size is applied and the scalar values within the filter confines are extracted. The values are distributed into a number of bins, which are converted into the FDF by counting frequencies in each bin and dividing by a relative bin width (bin width times the number of elements in the dataset).
where φ i is number of mixture fraction values in i histogram bin, ∆ bin is bin width used in histogram, n φ total number of mixture fraction values within the confines of the filter.
The shape of the beta density or the top-hat FDF is defined by the first two statistical moments of the mixture fraction. In LES, the resolved mean is available directly from the solution of the filtered transport equations, while a model is needed to approximate the 'subgrid' scale scalar variance. The SGS scalar variance transport equation can in principle be solved, which does not require modelling. However, the SGS scalar variance transport equation, which is regarded as the most accurate in the prediction of the SGS scalar variance, can in some cases result in a non realistic solution. It was showed that the variance transport equation tended to underpredict the SGS scalar variance due to numerical errors that were related to the numerical approximation of the SGS scalar transport equation 13 . The SGS scalar variance is important in the scalar dissipation rate modelling and hence an accurate prediction of the SGS scalar variance is paramount.
In this work, we use the 'true' measured SGS scalar variance to parametrize the presumed filtered density functions and calculate a modelled SGS scalar variance to compare with the experimentally measured one. We refer to the measured mixture fraction scalar variance within the box filter as the 'true' SGS scalar variance. The SGS scalar variance (z 2 ) can be computed directly from experimentally measured mixture fraction as follows: where ∆ is box filter width and z i is i−value of mixture fraction within the confines of the box filter The computation of SGS variance z 2 is similar to a low-pass filtering procedure and can be implemented directly. However, computation according to Eq. 11 can be time consuming and non-efficient. In order to decrease computation time it is convenient to rewrite the SGS scalar variance as. In order to simplify the implementation of the algorithm, Taking into account that the filter function is normalised, i.e. G ∆ (x − x )dx = 1, we obtain as, where the dependence on x and t is omitted for clarity.
We use equation 13 to compute the SGS scalar variance, as the measured squared instantaneous mixture fraction field convolved with a box filter minus the squared resolved mixture fraction field. The two SGS scalar variance models are also briefly introduced here for convenience. The first model is based on a gradient assumption 14 . According to this model, the SGS scalar variance is directly related to the gradient of the resolved scalar field via a modelling constant C z and the filter width ∆: The second SGS scalar variance model is based on the similarity between large and small scales 15 . The relationship between large and small scales is often approximated by the selfsimilarity assumption. Under this assumption the SGS scalar variance can be computed from larger scales by applying the so called test-level filtering. The test-level filtering operation is applied to the resolved field and the SGS scalar variance is computed according to the following formula, where ∧ denotes a test filter, which is often taken to be two times larger than the LES filter 15 .
The scalar similarity coefficient C ss can be prescribed or defined by using a dynamic procedure, which requires another test filter that must be larger than the test filter used in the model. However, the assumption of self-similarity is not correct for arbitrary LES filters.
It was shown that the base filter (LES filter) and the test filter had to be similar 16 . This means that if, for instance, a box filter is used in the LES simulation, the test filter must be also a box filter. The similarity concept and the similarity assumption model have been evaluated both numerically and experimentally 17 . It was generally found that this model performed better than gradient assumption models.
In addition, we quantify the probability of having bimodal FDF distribution at several locations in the flow. We compute the SGS scalar variance from 1500 instantaneous realisations for a fixed spatial location, fixed filter size and report the results in terms of the p.d.f. of the fluctuations of the logarithm of the 'subgrid' scale scalar variance. The p.d.f.
was computed from the logarithm of the instantaneous values of the 'subgrid' scale scalar variance minus the mean and normalizing by its standard deviation. If the resulting shape of the p.d.f of the SGS scalar variance resembles the Gaussian distribution, then the SGS scalar variance will be log-normally distributed.
Also, it is generally expected that if the LES filter size is chosen close to the Batchelor scale, the resulting FDF will be unimodal. However, this statement is not always true and depends on flow dynamics. We demonstrate that the FDF can be bimodal even when the filter size is comparable to the Batchelor scale λ β . In turn, the FDF can be unimodal when the filter size is 4-5 times larger than the Batchelor length scale.
where Pz 2 is probability density function of SGS scalar variance and z 2 is resolved subgrid scale scalar variance Even though several experimental attempts to measure filtered density function of mixture fraction in both reacting and non-reacting flows exist 8,18,19 , no measurements of FDF in swirling flows with low or high swirl number exist to the best of our knowledge. In this paper we analyse the properties of presumed beta density, top-hat FDF, measured FDF as well as the p.d.f. of the SGS scalar variance. We evaluate associated dependent variable integration errors that arise from the discrepancy between the measured FDF values and the theoretical presumed FDF output, especially for cases of bimodality. The remaining paper is structured as follows. The next section describes the swirling air flow geometry and the optical instrumentation used to measure the mixture fraction, which is then used to evaluate the presumed FDF modelling methods and the p.d.f. of the SGS scalar variance.
The last section describes the results and discusses the findings in terms of the differences between the measurements and the results from the presumed FDFs. The paper ends with a summary of the main conclusions.

II. EXPERIMENTAL SETUP
The burner comprised two concentric pipes with the annulus supplying air and the central pipe delivering air seeded with acetone vapour, which was the measured scalar quantity 20,21 .
The central pipe had an inner diameter D f = 15 mm and an outer diameter of 18 mm, was 0.75 m long and was located concentrically in the outer pipe of inner diameter, D = 50.8 mm and centred within it by three locking screws placed 25 mm upstream of the burner exit.
The annular flow development section was 0.264 m long. The annular air stream was split into two separately metered streams. One ('tangential air') passed through a static swirler containing six tangential slots to impart angular momentum. The static swirler was located in a plenum chamber in which the swirling air was combined with the second stream, which delivered 'axial air'. Metal plates were installed in the axial and tangential air sections of the plenum chamber to ensure that the axial and tangential air streams were distributed uniformly upstream of the inlets into the annular air supply stream of the burner, where they were combined to control the strength of swirl at the burner exit. The tangential and axial air flowrates were metered by flowmeters after correction to atmospheric pressure and temperature. The axial direction, parallel to the flow propagation, was denoted as 'y' and the radial direction, perpendicular to the main direction of the flow, as 'x' (Fig. 1). The swirl number is defined as follows 21 : where S is swirl number, R is radius of the outer pipe, D is diameter of the outer pipe, r i is radius of the inner pipe, ρ is density, U is axial velocity component, W is tangential velocity  A double frame CCD camera (LaVision Inc., Imager Intense) was used and was equipped with a Nikkor lens with focal length of 50 mm and f-number of 1.4; a 13 mm extension ring was also used after the camera lens for tighter focusing. A BG3 bandpass Schott filter was used in front of the camera lens in order to block any remaining 532 nm light that might interfere with the acquired signal. Any scattered light at 266 nm UV wavelength that might come from the test section was blocked by the glass material of the camera lens. The array of the CCD camera was 1376×1040 pixels and the pixel size was 6.45 µm. The optical magnification was determined by using a calibration target plate and was found to be 0.0263 mm/pixel. The laser sheet thickness was measured as 0.16 mm.
The measured FDFs were computed from the instantaneous images for two box filter sizes of 0.3 (≈ λ β ) and 1.0 mm (≈ 3.3λ β ). The smaller filter size was chosen to test an assumption that the LES results should be in better agreement with experiments if the LES filter size ∆ is chosen to be as close as possible to the Batchelor length scale. Filter sizes larger than 1 mm were not chosen, because bimodality could already be observed.
We measured the SGS scalar variance, the FDF and the probability density function of the SGS scalar variance from nine different spatial locations (Fig. 2) at every laser sheet downstream position (Fig. 1). In order to quantify the probability of having bimodal FDF distributions at specific spatial locations in the flow, the p.d.f of the SGS scalar variance was computed from 1500 images at all 9 locations.
We have chosen the dependent variable to be a temperature obtained from a laminar flamelet solution as a function of mixture fraction. We used Chemkin-Pro to compute the flame temperature of a methane/air diffusion flame at atmospheric pressure and temperature as a function of mixture fraction. The San Diego reaction mechanism, thermodynamic database and the transport database from UC San Diego were chosen for this calculation.
The laminar flamelet solution, obtained from the planar diffusion flame simulation, was used with the measured FDF, presumed top-hat and beta density and the associated integration errors were computed as follows: where ∆ β and ∆ l are laminar flamelet solution integration errors for presumed beta density and top-hat FDF, T fz , T β , T Π are resolved mean temperatures computed by integration of laminar flamelet solution using measured FDF, presumed beta density and top-hat FDF.

A. Effect of experimental noise on FDFs
Filtered density functions of mixture fraction computed from raw mixture fraction images do not usually correspond to the so-called 'true' filtered density functions. This is due to digital noise that is always present in the experimental data. The resulting shape of the corresponding filtered function is a combination of 'real' mixture fraction distribution within the confines of a filter and the digital noise. Depending on the amount of noise in the experimental data, the finally computed FDF may not reflect the actual mixture fraction distribution. It is, therefore, worth considering both raw and denoised mixture fraction images. The signal-to-noise ratio in the raw images is ≈100 and the spatial resolution is η M T F @2% ≈ 300 µm. We quantified the spatial resolution by measuring the Modulation Transfer Function (MTF) 22 using a slanted knife edge method 23 and, as a measure of the resolution, we used the width where the MTF is at 2% of its peak value. We corrected the effects of experimental noise and optical system blurring by using a Wiener filter, which relies on measured signal and noise spectra and on a model for the high-wavenumber decay of the 'true' signal spectrum 24,25 . The result of this signal reconstruction process can be bounded within computed limits and we calculated that in our measurements, depending on the swirl number, the accuracy of the scalar variance measurement is better than 2% and, typically, better than 1%. These errors are, also, within the bounds calculated by modelling the filtering procedure using model spectra and the present turbulent Reynolds number (based on the Taylor length scale) and MTF 26 . When we present our results, we also show the effect of the de-noising process on the measured small scale distributions.

III. RESULTS AND DISCUSSIONS
An example of an instantaneous mixture fraction distribution (before any spatial filtering) and the corresponding resolved means, for various filter widths, are shown in Figure 3, at axial distance y/D f = 1 and S=0.3. The resolved mean was computed for three different filter sizes (from top to bottom ∆ =1.0, 0.7, 0.5 mm; ≈ 3.3λ β , 2.3λ β , 1.7λ β ) as a convolution of an original measured mixture fraction field with a box filter. The convolution can be performed, e.g. by using the principle of discrete convolution (or moving average) or in the Fourier domain; we used discrete convolution in the spatial domain. Figure 3 shows that the filtering procedure is equivalent to a low-pass filtering, which results in smoother mixture fraction fields. The filter, in a broad sense, removes smaller scales resulting in a coarser representation of the mixture fraction field.
The subgrid scale scalar variance was computed for four different filter sizes using the bottom image in Figure 3. The measured SGS variance was computed by using Eq. 12 and is shown in Figure 4, which shows that larger box filter corresponds to larger SGS scalar variance. It is also expected that larger SGS scalar variance will change the shape of the measured filtered density functions. In addition, large SGS scalar variance is found in regions where large spatial gradients of the mixture fraction are found. Therefore, it is suggested that the SGS scalar variance (and hence the shape of the FDF) is a function of local mixing regime.   Figure 3, which is the measured reference image.

A. Measured filtered density function
In this section, we present measured filtered density functions of the mixture fraction and presumed filtered density functions and we study the effect of the box filter size on the resulting shape of the FDF. An example of filtered density function computation from an experimentally measured mixture fraction field for two different filter sizes is shown in Figure   5. The smallest filter, which was comparable to the Batchelor scale (≈300 µm) contained almost uniformly distributed mixture fraction values within the confines of the filter. In this case, the FDF would be expected to be unimodal, as it is illustrated in Figure 6. On the other hand, if the filter contains highly segregated mixture fraction values, i.e. with larger fluctuations around the mean, the resulting FDF would be expected to be bimodal ( Fig. 6 bottom left).   bimodal distribution is observable for the largest filter size and an approximately unimodal distribution for the smallest filter size. It should be noted that for fixed filter size, the shape of the FDF will be different as the filter moves along the spatial domain. This is also illustrated in Figure 4, where for fixed filter size, the SGS scalar variance is different and hence the FDF will also be different.
In addition, we studied the shape of the FDF as a function of the SGS scalar variance.
By using trial and error procedures we experimentally found a threshold SGS scalar variance value above which the FDF would be bimodal. The threshold was computed to be approximately z 2 =0.02. This value was found by plotting the measured FDF as a function of the SGS scalar variance and examining the resulting shape. We, also, define the term 'high' and 'low' SGS scalar variance as values that are higher or lower than a threshold value, namely z 2 > 0.02 and z 2 < 0.02. An interesting observation can be seen in Figure 8, which shows measured FDF computed using a filter size of ∆/λ β =1.5; the presumed beta density and top-hat FDF are also shown, which were calculated using the measured resolved mean and measured SGS scalar variance from the same filter. The figure shows that for this small filter, which is comparable to the Batchelor scale, the bimodal FDF distribution can still be observed. In fact, Figure 8 demonstrates the so-called 'filter overlapped state', which is schematically illustrated in the cartoon form in Figure 6 (left side).
Calculated presumed beta density and top-hat FDF are also plotted in Figure 8 and in Figure 9. For large SGS variance neither the beta density nor the top-hat FDF are able to correctly represent the true measured FDF. However, for small SGS variance, the beta density and the top-hat FDF coincide with the measured FDF becoming delta functions and almost indistinguishable.
In addition, the unimodal FDF distribution can be found, e.g. in the jet potential core, where no mixing takes place and hence for any arbitrary filter size, the distribution of the mixture fraction values within the confines of the filter is nearly homogeneous. Similarly, under well-mixed conditions, e.g. at far away distances from the nozzle exit or highly swirling flows, where rapid mixing takes place, the shape of the FDF is likely to be unimodal. Figure 9 shows an example of measured FDFs from two spatial locations for large filter sizes (∆/λ β = 4) and well-mixed regime at y/D f = 7 for S=0.58. As it can be clearly seen, even for large filter sizes, the shape of the FDF can be unimodal.  The effect of the swirl number cannot be readily assessed from the shape of the FDF because the FDF is by definition computed from the instantaneous realisations. General remarks, which will be supported later by examining the p.d.f. of the SGS scalar variance, are given here as follows. An increase in the level of mixing increases the homogeneity of the flow field, which in turn means that the scalar is locally well-mixed. A well-mixed regime indicates that the SGS scalar variance tends to be low, which results in a unimodal FDF distribution. It was generally found that unimodal distributions were observed at all locations in the flow at y/D f = 7 for S=1.07 and for all filter sizes. This is also consistent with Wang & Tong 9 where it was pointed out that the conditional FDF was found to be unimodal, and not far from Gaussian, when the SGS variance was small.
The SGS scalar variance is, therefore, the key parameter and is the measure of the level of mixing within one LES cell. It is convenient to relate the measured SGS scalar variance to specific positions in the flow or to operating conditions by considering the probability density function of the SGS scalar variance fluctuations. We report the p.d.f. of the SGS scalar variance, as an example, for one position in the flow (window 3 of Figure 2), according to Eq. 16. A box filter of given size, e.g. 1.0 mm, is fixed at a specific position in the flow and the SGS scalar variance is computed using the filter at every instantaneous realisation.
The SGS p.d.f. convergence was checked and we found that the p.d.f. of the SGS scalar variance converged for 1000 realisations rather than the full 1500 realisations. An example of a raw and the corresponding denoised mixture fraction image is shown in Figure 11 for S=1.07 at y/D f = 1, as well as the corresponding FDFs. It is seen that the level of denoising affects the resulting shape of filtered density functions. However, similar trends are generally observed for both raw and denoised FDFs, i.e. bimodality or unimodality. For the smallest filter size, which are comparable with the Batchelor scale, the FDF computed using the denoised mixture fraction distribution can result in unimodal distribution, while the FDF computed using the same filter at the same spatial location at the raw image can result a bimodal distribution. It should be noted that the experimentally measured FDF in Figure 11 corresponds to the 'true' filtered density function, which was corrected for both digital noise and optical system blurring.
Finally, we comment on the integration errors between the temperature computed using the presumed FDF and the corresponding measured FDF. For simplicity and for the sake of brevity, we consider the dependent variable to be the temperature obtained from the laminar flamelet solution. The resolved temperature was computed only for realisations with high SGS scalar variance z 2 >0.02, i.e. only when a bimodal FDF distribution was observed. The integration errors were below 1% and a direct correlation between the shape of the presumed beta density and the associated integration error was not found. The weak dependence between the temperatures computed using the presumed FDF and the measured FDF suggests that the SGS scalar variance models red, and not the presumed beta density, are perhaps the main source of error in the LES simulations that are based on presumed FDF methods. This is also consistent with Floyd 7 and Cook & Riley 15 .
In addition, we present several integrations of the laminar flamelet solution with the measured FDF, the beta density and the top-hat FDF in a table form. Table I  in resolved variables. Based on these facts, which are supported by Table I, it has been, therefore, suggested that the actual source of error in resolved variables (temperature, species fractions etc.) could be possibly found in the 'subgrid' scale models rather than in presumed FDF methods.

B. Assessment of SGS scalar variance models
We examine here the two common models for the SGS scalar variance: a) a gradient assumption model with prescribed constant (Eq. 14); b) a scale similarity approach with a prescribed constant (Eq. 15). Figure 12 shows scatter plots of the SGS scalar variance

IV. SUMMARY
These experiments focused on the filtered density functions (FDF) as well as on the 'subgrid' scale (SGS) scalar variance. The FDF were studied by analysing experimental data, obtained from two-dimensional planar, laser induced fluorescence scalar measurements in isothermal swirling turbulent flows. For a filter size, which is comparable to the Batchelor scale, the filtered density functions demonstrated unimodal distributions that were similar to delta functions. The bimodal distributions were observed for filter sizes as small as two times the Batchelor scale and even for ∆/λ β ≈1.5. Higher swirl numbers lead to well-mixed regime, resulting in unimodal FDFs, even for the filter size as large as four times the Batchelor scale. It was also demonstrated that the bimodal distribution of the measured FDF did not automatically introduce errors when integrating the laminar flamelet data. A denoising procedure was also used in order to assess the effect of digital noise on experimentally measured FDFs. It was shown that, even though the shape of the resulting FDF can be quite different, bimodality and unimodality were still observed. Therefore, the above mentioned conclusions are still applicable to raw (non-denoised) data. However, this statement needs to be checked against other denoising techniques as well as for different spatial resolutions and experimental configurations.
The 'subgrid' scale scalar (SGS) variance was also measured in this work for different filter sizes and swirl numbers. It was demonstrated that the SGS scalar variance increases with box filter size. The p.d.f. of the SGS scalar variance was also presented for a fixed filter size and a spatial location in the flow. It was demonstrated that the probability of having bimodal FDF distribution (high SGS variance) shifted to very low values for the higher swirl numbers of 0.58 and 1.07, which indicated that the distribution was generally unimodal. It was also shown that the p.d.f. of the SGS scalar variance was close to a log-normal distribution. A direct comparison between the measured SGS scalar variance and a modelled variance was provided. Two SGS variance models were chosen, namely a gradient assumption model and a scale similarity approach, for their simplicity and common use. The discrepancy between the measured SGS scalar variance and the model outputs was negligible with filter width close to the Batchelor scale, for both gradient assumption and scale similarity models. For larger filter size the gradient assumption model performed better than the scale similarity model for a single value of the modelling constant. Based on the fact that the bimodal FDFs, which were measured in this work, and the corresponding presumed FDF were not correlated in terms of associated laminar flamelet solution integration errors, it has been suggested that the actual source of error in resolved variables (temperature, species fractions etc.) could possibly be found in the 'subgrid' scale models rather than in presumed FDF methods.