Agnes L Bodor

Alberto Cuoci

Tiziano Faravelli

Benedetta Franzelli

A forward approach for the validation of soot size models using LII

While validating the numerical modeling of the primary particle size distribution (PPSD) in sooting flames, a common practice is to compare the numerical results to the corresponding experimental data obtained with the Time-Resolved Laser-Induced Incandescence (TiRe-LII) technique.

Since the PPSD is not directly measured by TiRe-LII, but derived with a post-processing procedure, various uncertainties and errors can potentially affect the consistency of such comparison requiring the estimation of many input parameters. On the contrary, nowadays detailed numerical simulations provide access to a more complete set of data, which can be used

to reconstruct the incandescence signal. In this work, a forward approach for the generic validation of numerical models for the PPSD is performed.

It is based on the numerical reconstruction of the temporal evolution of the incandescence from the numerical results and its comparison with the measured signal. First, two indexes are proposed to quantify the agreement between the numerically synthesized and the measured signals. Then, the effectiveness of the proposed approach is demonstrated a priori by quantifying the potential errors that can be avoided with this new strategy compared to the classical approach. Finally, the feasibility of the proposed procedure is proven by comparing synthesized signals to the experimental ones available in literature for a laminar premixed flame. It is shown that the proposed approach can be used to strengthening the analysis on numerical model performances, additionally to the classical approach.

Introduction

Among the different pollutants produced by combustion, soot is getting more and more attention from researchers, societal actors and industry.

From environmental and health perspectives, it is important to characterize the emitted particulate matter, not only in terms of the total mass of the particles, but also of their surface [START_REF] Hansen | Soot climate forcing via snow and ice albedos[END_REF][START_REF] Shiraiwa | Hazardous components and health effects of atmospheric aerosol particles: Reactive oxygen species, soot, polycyclic aromatic compounds and allergenic proteins[END_REF].

Unfortunately, the numerical prediction of soot particles' morphology represents a difficult challenge since it has been shown that the biggest particles, namely aggregates, have fractal-like nature [START_REF] Sorensen | The Morphology of Macroscopic Soot[END_REF]. In specific, soot aggregates are constituted of small quasi-spherical particles, called primary particles. The aggregates surface is then classically retrieved once the number and the size of the primary particles are known. For this reason, it is essential to dispose of experimental [START_REF] Kholghy | The evolution of soot morphology in a laminar coflow diffusion flame of a surrogate for Jet A-1[END_REF][START_REF] Kempema | Combined optical and TEM investigations for a detailed characterization of soot aggregate properties in a laminar coflow diffusion flame[END_REF][START_REF] De Iuliis | Three-angle scattering/extinction versus TEM measurements on soot in premixed ethylene/air flame[END_REF][START_REF] Abid | On evolution of particle size distribution functions of incipient soot in premixed ethylene-oxygen-argon flames[END_REF][START_REF] Zhao | A comparative study of nanoparticles in premixed flames by scanning mobility particle sizer, small angle neutron scattering, and transmission electron microscopy[END_REF][START_REF] De Iuliis | Effect of hydrogen addition on soot formation in an ethylene/air premixed flame[END_REF][START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] and numerical [START_REF] Park | An aerosol model to predict size and structure of soot particles[END_REF][START_REF] Patterson | Models for the aggregate structure of soot particles[END_REF][START_REF] Zhang | Implementation of an advanced fixed sectional aerosol dynamics model with soot aggregate formation in a laminar methane/air coflow diffusion flame[END_REF] tools allowing the characterization of the primary particle size distribution (PPSD).

Concerning the experimental determination of the primary particle size, the Transmission Electron Microscope (TEM) analysis of thermophoretically sampled soot has been extensively applied to determine the mean primary particle diameter d pp and the primary particle size distribution [START_REF] Kempema | Combined optical and TEM investigations for a detailed characterization of soot aggregate properties in a laminar coflow diffusion flame[END_REF][START_REF] De Iuliis | Three-angle scattering/extinction versus TEM measurements on soot in premixed ethylene/air flame[END_REF][START_REF] Abid | On evolution of particle size distribution functions of incipient soot in premixed ethylene-oxygen-argon flames[END_REF][START_REF] Zhao | A comparative study of nanoparticles in premixed flames by scanning mobility particle sizer, small angle neutron scattering, and transmission electron microscopy[END_REF][START_REF] De Iuliis | Effect of hydrogen addition on soot formation in an ethylene/air premixed flame[END_REF].

When using the TEM technique, the d pp is directly measured from pictures taken by microscope. Such technique brings information about the shape of the particles, such as the transition from spherical particles to aggregates.

The TEM is one of the best techniques to visualize soot particles down to the size of ∼10 nm [START_REF] Zhao | A comparative study of nanoparticles in premixed flames by scanning mobility particle sizer, small angle neutron scattering, and transmission electron microscopy[END_REF]. In addition to this technique, the Helium-Ion Microscopy (HIM) was developed for microscopy investigations at the nanoscales [START_REF] Ward | Helium ion microscope: A new tool for nanoscale microscopy and metrology[END_REF].

The retrieved information is similar to TEM, but HIM is classically used to characterize soot particles as small as 2 nm [START_REF] Schenk | Morphology of nascent soot in ethylene flames[END_REF][START_REF] Schenk | Imaging Nanocarbon Materials: Soot Particles in Flames are Not Structurally Homogeneous[END_REF]. Despite their relevant advantages, these techniques present some similar drawbacks. The measurements are performed intrusively and point-wise, limiting the analysis to few points in the flame. In addition, results depend on the algorithm and/or on the researcher criteria used to visually determine the number and the size of the primary particles in an aggregate.

Another option is the use of the Scanning Mobility Particle Sizer (SMPS), which provides information about the particle mobility diameter [START_REF] Sorensen | The mobility of fractal aggregates: A review[END_REF][START_REF] Camacho | Mobility size and mass of nascent soot particles in a benchmark premixed ethylene flame[END_REF]. In case of spherical particles, the primary particle size and the particle size coincide so that this technique can be considered for the determination of the size of the primary particles. On the contrary, the extraction of d pp from SMPS distributions remains very challenging for large aggregates.

Alternatively to these aforementioned intrusive techniques, optical diagnostics represent a very accurate experimental approach suitable to perform investigations even in the harsh conditions of combustion. Among them, the Time-Resolved Laser-Induced Incandescence (TiRe-LII) is an in-situ minimally invasive diagnostic technique that is more and more often used to retrieve information on d pp and PPSD [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Chen | Determination of soot particle size using time-gated laser-induced incandescence images[END_REF]. In addition, the Planar Laser-Induced Incandescence can be used to obtain information on the spatial distribution of the PPSD in the whole laser plane.

Due to their advantageous properties, the TiRe-LII technique and laser diagnostics in general are powerful tools to investigate soot formation. However, when using laser diagnostics, the fundamental quantities of interest are usually not directly measured, but they are the results of a post-processing procedure applied to the detected signal (inverse method). Most of the time, such post-processing requires the quantification of various quantities simultaneously, leading to additional measurements, post-processing and/or assumptions. Unfortunately, each of them may be considered as a new source of uncertainty. Concerning the TiRe-LII approach, a long list of quantities are required to obtain the d pp and the PPSD from the measured incan-descence signal. Beside the LII signal decay, information on the gas phase, such as its temperature or pressure, and on the soot population, such as the number of particles per aggregate, are needed to correctly post-process the signal. Disposing of all these quantities simultaneously is more than challenging, since most of them are not even directly measured. This greatly complexifies the validation of the numerical models for the PPSD, which usually relies on comparing numerical results to experimental data.

On the contrary, recent numerical simulations can provide access to a wide range of information. For example, Salenbauch et al. [START_REF] Salenbauch | Detailed modeling of soot particle formation and comparison to optical diagnostics and size distribution measurements in premixed flames using a method of moments[END_REF] proposed to directly compare the numerical spatial evolution of the mass concentration of small and large particles to the spatial evolution of the Laser induced Fluorescence and LII signals, respectively. Unfortunately, such comparison is not quantitative since it is not based on the same quantities (numerical particles concentration on the one side, experimental signals on the other size).

Alternatively, the numerical data can be utilized to reconstruct the experimental signal. A direct comparison between the numerically synthesized and the experimental signals can then be performed (forward method). Examples of a forward method for the validation of numerical simulations can be found in literature [START_REF] Boyce | Computational fluid dynamics validation using multiple interferometric views of a hypersonic flowfield[END_REF][START_REF] Danehy | Numerical Simulation of Laser-Induced Fluorescence Imaging in Shock-Layer Flows[END_REF][START_REF] Connelly | A paradigm shift in the interaction of experiments and computations in combustion research[END_REF]. However, to the author's knowledge, no previous effort on comparison of synthesized and experimental TiRe-LII signals has been done in literature for the validation of the numerical models in terms of PPSD.

The objective of this work is therefore to investigate the potential of a forward method for the validation of the numerically modeled PPSD, by comparing numerically synthetised signals to TiRe-LII experimental data.

In this context, it is worth to be mentioned that the reconstruction of the PPSD from the experimental TiRe-LII signal is strongly affected by the post-processing and the experimental uncertainties. The aim of this paper is not to improve the quality of such reconstruction but to propose a new post-processing technique allowing for a consistent comparison of experimental and numerical results to strengthen the validation of the numerical modeling.

The structure of the paper is as follows. The uncertainties that are potentially introduced by the post-processing of the TiRe-LII signal to obtain the PPSD are quantified in Section 2 by reviewing results from the literature. Two indexes are suggested to characterize the LII signal in Section 3. An uncertainty evaluation on these quantities is carried out to a priori identify the possible advantages of the validation through synthesized TiRe-LII signals. Then, both the inverse and the forward validation methods are applied to a premixed ethylene flame in Section 4. The input numerical results are obtained from detailed simulations already validated in [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF] and the experimental LII signals are extracted from [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF]. Differences between measured and synthesized signals are highlighted and conclusions on the potential of the proposed approach for validating numerical models for the PPSD prediction are finally reported.

Literature overview on the current accuracy of d pp determination by TiRe-LII

The TiRe-LII technique relies on the heating up of soot particles via a nanosecond laser pulse and the recording of the subsequent induced incandescence signal [START_REF] Michelsen | Understanding and predicting the temporal response of laser-induced incandescence from carbonaceous particles[END_REF]. As large particles cool more slowly than small ones, the temporal evolution of the incandescence signal can be used to determine the primary particle size.

In order to obtain the PPSD from the signal decay, a post-processing inverse procedure is needed. The LII signal is reconstructed by solving the mass and the energy balance equations for theoretical PPSDs. Usually, the problem is recasted as a minimization problem, which allows to identify the distribution parameters minimizing the deviation between the measured signal and the theoretical one [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF][START_REF] Betrancourt | Investigation of the size of the incandescent incipient soot particles in premixed sooting and nucleation flames of n-butane using LII, HIM, and 1 nm-SMPS[END_REF][START_REF] Dankers | Determination of primary particle size distributions from time-resolved laser-induced incandescence measurements[END_REF][START_REF] Franzelli | Multi-diagnostic soot measurements in a laminar diffusion flame to assess the ISF database consistency[END_REF][START_REF] Foo | Experimental investigation of acoustic forcing on temperature, soot volume fraction and primary particle diameter in non-premixed laminar flames[END_REF]. Such post-processing can be affected by two main sources of uncertainties. First, the reconstruction of the LII signal is known to be strongly dependent on the model and the parameters (such as soot density) selected for the mass and the energy balance equations as exhaustively discussed by Schulz et al. [START_REF] Schulz | Laser-induced incandescence: Recent trends and current questions[END_REF] and Michelsen et al. [START_REF] Michelsen | Modeling laser-induced incandescence of soot: A summary and comparison of LII models[END_REF]. Second, most of the models rely on various input quantities, summarized in Table 1, which are not necessarily known. The role of each parameter in the signal modeling is not described here in details, but can be found in [START_REF] Schulz | Laser-induced incandescence: Recent trends and current questions[END_REF][START_REF] Michelsen | Modeling laser-induced incandescence of soot: A summary and comparison of LII models[END_REF]. Many of these variables are properties of the experimental equipments ("Equipment"), some of them are properties of the gaseous field around the investigated soot particles ("Ambient") and numerous of them are characteristics of the soot particles themselves ("Soot"). The errors introduced by the parameters of the equipment are supposed to be limited by accurately designing the experimental setup and are not discussed here. 1 Some of the remaining variables can be measured ("M"), but this may possibly introduce additional errors and requires additional devices and/or post-processing. Some other parameters are generally assumed ("As").

With the forward approach proposed in Section 3, uncertainties related to the LII model cannot still be avoided but some errors due to the input parameters can be reduced by using the values available in the numerical simulations. The concerned parameters are marked by "N". As, M PPSD shape As, M, N Number of primary particles in the aggregate (n p ) As, M, N Thermal accommodation coefficient (α T )

Equipment

As, M Table 1 Main input parameters needed for the modeling of the temporal evolution of the TiRe-LII signal according to [START_REF] Hofmann | A web-based interface for modeling laser-induced incandescence (LIISim)[END_REF].

1 As an example, a non-homogenous profile of the laser beam can potentially affect the post-processing results [START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF]. In the proposed forward technique, it is possible to account for it, reducing the uncertainties on the comparison between experimental and numerical results.

Additionally to the five variables that are known from the numerical simulations (pressure p, temperature T , molecular weight of gas phase M g , PPSD shape and number of primary particles in the aggregate n p ), the thermal accommodation coefficient and the soot absorption function are today not well known but they are expected to depend on soot properties. Therefore, these two last quantities may be potentially recovered from numerical simulations once the relations with other particles characteristics are identified. In the following, the effect of the uncertainties introduced by these parameters on d pp are recalled from earlier studies.

2.1 Assumed or measured parameters 2.1.1 Soot absorption coefficient Various values of the soot absorption function E(m), constant or wavelength-dependent, can be found in literature [START_REF] Michelsen | Modeling laser-induced incandescence of soot: A summary and comparison of LII models[END_REF][START_REF] Migliorini | Investigation of optical properties of aging soot[END_REF][START_REF] Snelling | Determination of the soot absorption function and thermal accommodation coefficient using low-fluence LII in a laminar coflow ethylene diffusion flame[END_REF][START_REF] Krishnan | Extinction and Scattering Properties of Soot Emitted From Buoyant Turbulent Diffusion Flames[END_REF]. In particular, it has been found that E(m) varies between ∼0.2 to ∼0.4 even for a fixed laser wavelength of 1064 nm [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF]. However, such high uncertainty does not greatly affect PPSD results since many of the LII models [START_REF] Dankers | Determination of primary particle size distributions from time-resolved laser-induced incandescence measurements[END_REF][START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF] are insensitive to E(m) [START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF]. More specifically, it has been shown in [START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF] that a variation of ±0.1 from the reference value (E(m) = 0.3) did not lead to significant modifications of the PPSD, except when using the model of Lehre et al. [START_REF] Lehre | Size distributions of nanoscaled particles and gas temperatures from time-resolved laser-inducedincandescence measurements[END_REF].

Mass accommodation coefficient

The mass accommodation coefficient is a key parameter for the description of particle sublimation. In the literature, this quantity is characterized by a high uncertainty: its value varies between 0.5 and 1 [START_REF] Michelsen | Modeling laser-induced incandescence of soot: A summary and comparison of LII models[END_REF]. To avoid the possible errors due to an inaccurate estimation of the mass accommodation coefficient, recent experimental studies suggest to work with a low laser fluence so that the sublimation phenomenon is negligible compared to the other processes occurring after the interaction of the laser with the particles [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Michelsen | Understanding and predicting the temporal response of laser-induced incandescence from carbonaceous particles[END_REF][START_REF] Michelsen | Modeling laser-induced incandescence of soot: A summary and comparison of LII models[END_REF][START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF][START_REF] Liu | Effects of primary particle diameter and aggregate size distribution on the temperature of soot particles heated by pulsed lasers[END_REF]. Uncertainties on β can then be neglected.

Thermal accommodation coefficient

In the literature the thermal accommodation coefficient α T is characterized by a large variability [START_REF] Michelsen | Laser-induced incandescence: Particulate diagnostics for combustion, atmospheric, and industrial applications[END_REF]. Values from 0.23 to 0.9 are classically used [START_REF] Michelsen | Modeling laser-induced incandescence of soot: A summary and comparison of LII models[END_REF]. The wrong assumption can lead to significant inaccuracies in the estimation of d pp , since a variation of just ±0.05 can result in significantly different PPSDs [START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF]. As long as the connection between α T and other soot properties is not established the uncertainty related to this parameter is not avoidable. However, once the relation between α T and other particles characteristics will be known, the LII signal reconstruction from numerical results can potentially overcome this issue.

2.2 Numerically available parameters 2.2.1 Gas temperature Several researches concern the impact of the uncertainty of the gas phase temperature T g on the experimental estimation of d pp from the Tire-LII signal. The exact deviation depends on the used models and only some examples are reported in the following.

Concerning a monodisperse distribution, Sun et al. [START_REF] Sun | Singleshot, Time-Resolved planar Laser-Induced Incandescence (TiRe-LII) for soot primary particle sizing in flames[END_REF] demonstrated that a variation of ±200 K from 1700 K (±12%) leads to errors in ±3 nm for d pp = 15 nm (±20%) and ±5 nm for d pp = 40 nm (±13%). The significant impact of the flame temperature on the derived primary particle size was presented also by Will et al. [START_REF] Will | Performance characteristics of soot primary particle size measurements by time-resolved laser-induced incandescence[END_REF]. For a variation of ±180 K from 1800 K (±10%) the relative error of the diameter was about ±6 nm (±30%) and ±8 nm (±20%) for a particle of 20 and 40 nm, respectively. Cenker et al. [START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF] observed a deviation up to ±20% in the primary particle diameter for a temperature variation of ±15%.

Concerning lognormal distributions, Daun et al. [START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF] quantified the differences between two PPSDs derived from the same LII signal when perturbing the temperature around its exact value. The measure of difference was the Cramér-von Mises (CVM) goodness-of-fit parameter, which was defined as the area contained between the recovered and the exact particle size cumulative distribution functions, i.e. a zero CVM value means that the two PPSDs are identical. The CVM parameter for a temperature perturbation of ±30 K from the nominal 1700 K varied between ∼0.1 to ∼0.43 2 , depending on the model used for the LII signal calculation [START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF].

Experimental uncertainties on temperature depend both on the used experimental technique as well as the investigated configuration with an error that can reach ±200 K [START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Franzelli | Multi-diagnostic soot measurements in a laminar diffusion flame to assess the ISF database consistency[END_REF][START_REF] Liu | Soot temperature and volume fraction retrieval from spectrally resolved flame emission measurement in laminar axisymmetric coflow diffusion flames: Effect of self-absorption[END_REF][START_REF] Cignoli | Two-dimensional twowavelength emission technique for soot diagnostics[END_REF][START_REF] Ma | Absolute light calibration using S-type thermocouples[END_REF][START_REF] Guo | Digital camera measurements of soot temperature and soot volume fraction in axisymmetric flames[END_REF][START_REF] Stasio | Influence of the soot property uncertainties in temperature and volume-fraction measurements by two-colour pyrometry[END_REF][START_REF] Kuhn | Soot and thin-filament pyrometry using a color digital camera[END_REF][START_REF] Jenkins | Soot pyrometry using modulated absorption/emission[END_REF]. Even worse, the temperature value is often assumed when extracting PPSD information from TiRe-LII signal [START_REF] Franzelli | Multi-diagnostic soot measurements in a laminar diffusion flame to assess the ISF database consistency[END_REF][START_REF] Wu | Particle size distribution of soot from a laminar/diffusion flame[END_REF][START_REF] Stirn | Comparison of particle size measurements with laser-induced incandescence, mass spectroscopy, and scanning mobility particle sizing in a laminar premixed ethylene-air flame[END_REF] greatly increasing the error bars on d pp . On the contrary, the numerical simulations provide the gas phase temperature in the whole domain, so that the uncertainty on this quantity can be reduced by numerically reconstructing the LII signal and comparing it to the experimental one (supposing that the numerical model has been previously validated for the gaseous phase description).

PPSD shape

The presumed PPSD shape is one of the most relevant uncertainties. In many cases a monodisperse distribution is assumed for the sake of simplicity [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF]. However, the primary particle distribution is often observed to be a lognormal function [START_REF] Kempema | Combined optical and TEM investigations for a detailed characterization of soot aggregate properties in a laminar coflow diffusion flame[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Lehre | Size distributions of nanoscaled particles and gas temperatures from time-resolved laser-inducedincandescence measurements[END_REF], characterized by a geometric mean diameter d g and standard deviation σ g :

f (d) = 1 √ 2πd ln σ g exp - (ln d -ln d g ) 2 2(ln σ g ) 2 . (1)
When assuming a monodisperse distribution, all particles have the same diameter d mono . In case of a lognormal distribution, the relation between a monodisperse equivalent mean particle diameter d eq mono (obtained by preserving the total volume and primary particle number but assuming a monodisperse distribution) and d g is related to σ g , as depicted in Figure 1. Depending on σ g , the ratio of the two characteristic diameters can vary significantly.

An exhaustive analysis was performed by Cenker et al. [START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF] on the pos-

1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 (d g -d mono )/ d g [%] σ g [-]
Fig. 1 The relation between the monodisperse equivalent mean particle diameter d eq mono and the geometric mean diameter dg for lognormal PPSD distribution as function of the standard deviation σg.

distribution with a wide range of {d g ,σ g } pairs (at 1750 K, 1 bar and n p =1).

In this work, the significant difference between d mono and d g extracted from the same signal was presented. As an example, for {d g ,σ g }={30 nm,1.4 nm} the resultant d mono was 40 nm. A similar study was carried out by Chen et al. [START_REF] Chen | Determination of soot particle size using time-gated laser-induced incandescence images[END_REF] with similar results. Sun et al. [START_REF] Sun | Singleshot, Time-Resolved planar Laser-Induced Incandescence (TiRe-LII) for soot primary particle sizing in flames[END_REF] also reevaluated a typical TiRe-LII signal with both monodisperse and lognormal distribution assumptions (at 1700 K, 1 bar and n p =100). The retrieved characteristic diameters were d mono = 19 nm and d g = 15 nm with σ g = 1.25 nm. However, it is important to mention that in both investigations [START_REF] Sun | Singleshot, Time-Resolved planar Laser-Induced Incandescence (TiRe-LII) for soot primary particle sizing in flames[END_REF][START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF] On the one side, all these works enlighten the fact that accounting for polydispersity may drastically change the reconstructed mean d pp value and that a consistent definition between the experimental and the numerical mean d pp diameter is necessary. On the other side, Daun et al. [START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF] showed that when the lognormal PPSD shape is valid entirely different {d g ,σ g } pairs minimize the objective function of the iterative process. This means that various PPSDs can be considered as a good match to the experimental signal taking into account the experimental uncertainties, greatly weakening the validation of the numerical prediction in terms of the PPSD.

Finally, it has been shown that a lognormal PPSD shape cannot be considered as a generality [START_REF] Kholghy | The evolution of soot morphology in a laminar coflow diffusion flame of a surrogate for Jet A-1[END_REF][START_REF] Kempema | Combined optical and TEM investigations for a detailed characterization of soot aggregate properties in a laminar coflow diffusion flame[END_REF][START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF][START_REF] Vargas | Pressure dependence of primary soot particle size determined using thermophoretic sampling in laminar methane-air diffusion flames[END_REF][START_REF] Liu | Effects of gas and soot radiation on soot formation in a coflow laminar ethylene diffusion flame[END_REF], so that no shape could be presumed with confidence even if it is necessary when post-processing TiRe-LII data.

By reconstructing the synthesized signal from the numerical results, uncertainties on the PPSD shape can be strongly limited improving the comparison accuracy.

Number of primary particles in the aggregate

In most of the cases, single spherical particles are assumed for the TiRe-LII signal evaluation [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Franzelli | Multi-diagnostic soot measurements in a laminar diffusion flame to assess the ISF database consistency[END_REF][START_REF] Sun | Singleshot, Time-Resolved planar Laser-Induced Incandescence (TiRe-LII) for soot primary particle sizing in flames[END_REF], even if the cooling process strongly depends on the number of primary particles per aggregate (n p ) [START_REF] Chen | Determination of soot particle size using time-gated laser-induced incandescence images[END_REF][START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF][START_REF] Michelsen | Laser-induced incandescence: Particulate diagnostics for combustion, atmospheric, and industrial applications[END_REF][START_REF] Bladh | On the dependence of the laserinduced incandescence (LII) signal on soot volume fraction for variations in particle size[END_REF]. Liu et al. [START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF][START_REF] Liu | Effects of primary particle diameter and aggregate size distribution on the temperature of soot particles heated by pulsed lasers[END_REF] studied the relation between n p and the signal decay, providing a method to account for the shielding effect. 3Accounting for n p may have a strong impact on the results [START_REF] Crosland | Instantaneous in-flame measurement of soot volume fraction, primary particle diameter, and aggregate radius of gyration via auto-compensation laser-induced incandescence and two-angle elastic light scattering[END_REF][START_REF] Bauer | Can soot primary particle size distributions be determined using laser-induced incandescence?[END_REF]. As an example, the reevaluation of the measured TiRe-LII signal with n p =100 instead of n p =1 by Bladh et al. [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] led to the drop of the derived mean d pp by 30%. The same exercise was performed by Sun et al. [START_REF] Sun | Singleshot, Time-Resolved planar Laser-Induced Incandescence (TiRe-LII) for soot primary particle sizing in flames[END_REF], where the difference between the two derived diameters is more than 30% (29 nm for n p = 1 and 19 nm for n p = 100).

Unfortunately, as long as the LII experiment is not accompanied by TEM measurements, the aggregate structure is unknown, leading to uncertainties on n p and consequently on d pp . On the contrary, the numerical prediction provides information about the number of primary particles in the aggregate, so that the LII signal reconstruction can account for the phenomenon, allowing for a pertinent comparison of experiments and numerical predictions. The novel validation approach is then expected to reduce the n p related experimental uncertainties.

Smallest detectable particle

Once the mean d pp is obtained by postprocessing the TiRe-LII signal, the comparison with the numerical results may be affected by an additional parameter: the smallest particle d det,min that can emit a detectable LII signal. While earlier studies indicated that the smallest soot particles detected by TiRe-LII are in the range of 6-10 nm [START_REF] Abid | Size distribution and morphology of nascent soot in premixed ethylene flames with and without benzene doping[END_REF][START_REF] Conturso | Effect of alkylated aromatics on particle formation in diffusion flames: An experimental study[END_REF][START_REF] Bladh | Probing the smallest soot particles in lowsooting premixed flames using laser-induced incandescence[END_REF], a recent study of Betrancourt et al. [START_REF] Betrancourt | Investigation of the size of the incandescent incipient soot particles in premixed sooting and nucleation flames of n-butane using LII, HIM, and 1 nm-SMPS[END_REF] pointed out that soot particles as small as 2 nm can absorb laser energy and emit an LII signal.

However, the lower detection limit may depend on the sensitivity of the system, as it was pointed out by Sirignano et al. [START_REF] Sirignano | Detection of nanostructures and soot in laminar premixed flames[END_REF]. Even if the value of d det,min only slightly affects the soot volume fraction determination since the contribution of small particles is negligible [START_REF] Schulz | Laser-induced incandescence: Recent trends and current questions[END_REF], it can generally affect the value of d pp calculated from the numerical PPSD [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF]. Unfortunately, d det,min is not available for all experiments, potentially affecting the accuracy of the inverse validation method. On the contrary, small particles do not influence significantly the LII signal so that the comparison of synthesized LII signals to experimental ones is expected to be less affected by

d det,min .
3

The forward method and its potential improvement

In the forward method, the temporal evolution of the LII signal is first calculated for the numerically predicted PPSD extracting the quantities needed by the LII modeling from the numerical simulations. Then, the synthesized signal is compared with the experimental data. For this, different moments M k of the LII signal I LII (t) may be considered:

M k = t k I LII (t)dt (2)
with t denoting the time starting from the peak of LII signal (t = 0 ns when I LII is maximum). In the following, two indexes have been kept: the expected value

E = M 1 M 0 (3)
and the non-dimensional standard deviation

σ = σ E (4)
with the standard deviation

σ = M 2 /M 0 -(M 1 /M 0 ) 2 .
The skewness and the kurtosis indexes have also been considered but they did not provide a better description of the LII signal compared to E and σ , so that they have not been used in the following analysis. The LII signal is reconstructed by an in-house code utilizing the signal modeling description provided by Hofmann et al. [START_REF] Hofmann | A web-based interface for modeling laser-induced incandescence (LIISim)[END_REF] and the conduction model developed by Liu et al. [START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF][START_REF] Liu | Effects of primary particle diameter and aggregate size distribution on the temperature of soot particles heated by pulsed lasers[END_REF].

Effect of exponential fitting of the signal

Cenker et al. [START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF] investigated the error introduced by approximating the LII signal by a single exponential fit. They found that it can be considered as a good approximation in case of a monodisperse distribution, so that the characteristic decay time provides sufficient information to determine the diameter. However, as already said, the PPSD is rarely monodisperse and from the theoretical point of view the decay rate is not a single exponential function. For polydisperse distributions, the problem becomes more complex as the signal is the superposition of several quasi-single exponential signals with various decay times. To overcome this issue, Dankers and Leipertz [START_REF] Dankers | Determination of primary particle size distributions from time-resolved laser-induced incandescence measurements[END_REF] suggested to derive two characteristic LII signal decay times, τ 1 (∆t 1 ) and τ 2 (∆t 2 ) for two subsequent measurement time intervals (∆t 1 and ∆t 2 ), resulting in a double exponential function. They demonstrated that for a lognormal distribution the correlation between {τ 1 ,τ 2 } and {d g , σ g } is unambiguous when the signal noise is not considered. This approach has then been used in recent studies [START_REF] Chen | Determination of soot particle size using time-gated laser-induced incandescence images[END_REF][START_REF] Franzelli | Multi-diagnostic soot measurements in a laminar diffusion flame to assess the ISF database consistency[END_REF][START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF][START_REF] Wu | Particle size distribution of soot from a laminar/diffusion flame[END_REF].

However, both fits are approximations of the original signal and may lead to the loss of information. In the following, the effect of the approximation is discussed by comparing the E and σ indexes calculated from the original signal (subscript LII ), the single-(subscript SExp) and the double-(subscript DExp) exponential fitted signals for both monodisperse and lognormal distributions. This discussion will also clarify how to interpret the E and σ indexes and will quantify the errors introduced by the fitting procedure. Fig. 2 Expected value E of the signal (LII) and the curves obtained by the singleexponential (SExp) and double-exponential (DExp) fit. The decay time τ of the exponential fit is also provided.

Monodisperse distribution

10 to 30%, which confirms that for non-monodisperse distributions the single exponential fit cannot properly reproduce the LII signal.

In Figure 3(a-c), the time interval was chosen to avoid the effect of a limited detection time (18000 ns). However, in practical applications, the detection time is finite. To characterize the potential error, the evaluation was also performed by considering the signal only up to 800 ns, corresponding to the detection time of the experimental results [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] that will be used in Sec. 4. The resultant E LII is presented in Figure 3(d). It can be observed that for large diameters the values are significantly lower than at the full length of the signal, indicating that the choice of the detection time can affect the results. However, the numerical simulation can imitate the detection length and its temporal resolution, so that the validation is not affected by this deviation when using the forward approach. As the same E LII may originate from various {d g , σ g } pairs, a second index should be considered. In Figure 4, σ is provided for the original signal and for the double-exponentially fitted curve σ DExp . For a singleexponentially fitted curve, the value is σ SExp = 1 for all {d g , σ g }. Results for the original signal show that σ LII is almost insensitive to d g for d g > 10 nm, whereas it strongly varies with σ g . Therefore, the σ LII value may be used in the inverse approach to identify the σ g and, subsequently, detecting d g from E LII . However, it should be reminded that such procedure is not valid for small d g values (below 10 nm) since the σ LII is also strongly d g -dependent (Figure 4 imental signals can be directly compared without any need to retrieve d g and σ g . In addition, it is important to point out that a limited detection time provides a significantly different σ LII profile for large diameters, as observed in 4(d). In particular, the dependence on σ g is inverted since for increasing σ g the σ LII decreases. It is therefore clear that in most pratical situations, where a limited detection time is used for the experiments, it would not be possible to recover d g and σ g from E LII and σ LII , proving once again the limits of the inverse approach. Alternatevely, a longer detection time could be considered by using multiple detectors as suggested by [START_REF] Mansmann | Sequential signal detection for high dynamic range time-resolved laser-induced incandescence[END_REF]. On the contrary, using the forward approach, it would be possible to compare the experimental and numerical synthetized signals to quantify the inaccuracy of the soot model. The double-exponential fit shows a similar behaviour, but does not recover the same values of the original signal. In particular, for high σ g , it significantly under-predicts the index highlighting that the double-exponential fitting is not adequate for highly polydisperse distributions. Therefore, the d g and σ g values obtained in the inverse strategy may be affected by the double-exponential fitting procedure.

In conclusion, the use of a forward comparison of synthesized and experimental signals, characterized by their E LII and σ LII , is expected to reduce the uncertainties on the validation procedure. Such post-processing approach can be applied to any PPSD shape and the signals can be reconstructed to meet the experimental resolution and detection time. In addition, it controls the errors related to the experimental fitting of the signals, even if additional indexes may be needed to improve the characterization of the signals. In this sense, the skewness and kurtosis indexes have also been tested but it was found that they were not of interest for the considered signals. In the following, the evolution of E LII and σ LII with variations of the numerically retrievable input parameters will be discussed to quantify a priori the additional potential improvements of the forward approach compared to an inverse comparison.

Sensitivity study

In order to quantify the potential errors avoided when using a forward procedure, the impact of the variations of the input parameters on the reconstruction of the LII signal is quantified here. For this, a reference signal is produced with a set of base parameters Φ0 : Φ0 = (n p,0 , T 0 , M g,0 , p 0 )

for four different d pp values (5, 15, 25 and 35 nm) with a monodisperse distribution (where not stated differently). The reference parameters are:

n p,0 = 1, T 0 = 1700 K, M g,0 = 28 kg/kmol, and p 0 = 1 bar. Then, the parameters of Φ = (n p , T, M g , p) are perturbed one-by-one 4 and the change in the indexes is evaluated as:

∆E LII E LII,0 = E LII ( Φ) -E LII ( Φ0 ) E LII ( Φ0 ) (6)
The signals are calculated based on the physical model described by Hofmann et al. [START_REF] Hofmann | A web-based interface for modeling laser-induced incandescence (LIISim)[END_REF] and the conduction model by Liu et al. [START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF][START_REF] Liu | Effects of primary particle diameter and aggregate size distribution on the temperature of soot particles heated by pulsed lasers[END_REF]. The parameters not recoverable from the numerical simulation (α T and E(m)) were kept constant, with α T = 0.37 and E(m) = 0.32.

Temperature

The variation of E LII , introduced by a temperature perturbation, is shown in Figure 5 along with the relative errors. The impact of the temperature modification is large. The relative error of E LII is in the range of ±20% for the large particle diameters, and it is even larger for 4 The sensitivity studies of the following parameters were carried out by varying simultaneously the parameter in focus and the temperature to examine if they affect each other. The results showed that the error was always the superposition of the two individual errors. Therefore, for the clarity of the plots, only the single parameter perturbation results are presented. the 5 nm diameter particle (almost 40%, when the temperature is modified by 10%). From results in Fig. 5, it can be observed that by underestimating/overestimating the gas temperature, significantly higher/lower diameter are derived. The use of a proper temperature is therefore essential to minimize errors on E LII . As discussed before, the inverse approach requires the experimental measurements of temperature. Alternatevely, the forward approach allows the use of the numerical temperature to perform the validation of the numerical results. Of course, the accuracy of the numerically synthetized E LII depends on the quality of the numerical prediction for both the gaseous and the solid phases. In this sense, the forward approach could be used to quantify the inaccuracy of the whole numerical strategy and not of the single soot phase modeling.

Pressure

A sensitivity analysis is performed by applying a perturbation of ±2.5% and ±5% with respect to 1 bar. In Figure 6, E LII and its relative error are shown for the pressure perturbation. It can be observed in Figure 6 that E LII does not vary significantly with the pressure for small d pp . On the contrary, for higher values of d pp , the E LII becomes more sensitive to the pressure and the deviation is not negligible anymore. However, when looking at the relative error it can be deduced that the perturbation effect on E LII is almost d pp independent and reaches up to 5% for a pressure variation of 5%. Supposing that the dependence between E LII and d mono is linear, similarly to the results of Fig. 2, this means an error in the diameter up to 5%. Such error can be avoided with the novel forward approach, since information on pressure is classically available in the numerical simulations.

Molecular weight

To experimentally evaluate the gas phase molecular weight M g , information on the local composition is required. Many species can be measured with high accuracy. However, extra experimental equipment and post-processing are needed to access these quantities. Even when using ethylene as fuel (whose molecular weight is similar to nitrogen), the gas phase molecular weight can differ significantly from the one of air in the sooty region since soot is accompanied by several high molecular weight gas phase species, like heavy PAHs, but also lighter species as hydrogen, water vapor or acetylene. Therefore, the generally used approach, i.e. assuming air properties when retrieving the PPSD from the measured LII signal [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] Franzelli | Multi-diagnostic soot measurements in a laminar diffusion flame to assess the ISF database consistency[END_REF], is inaccurate. As an example, the M g value has been extracted from the numerical simulations of different sooting flames at atmospheric pressure [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF]. For a premixed ethylene laminar flame with an equivalence ratio of 2.1 and co-flow ethylene laminar flames with various ethylene content in the fuel stream (32%, 40%, 60% and 80 % volumetric ratio), the M g value can vary between 25kg/mol to 30 kg/kmol in the sooty region, whereas the general assumption is M g = 28 kg/kmol. The sensitivity study presented in Figure 7 reveals that, though other parameters (like n p or temperature) cause higher variations in E LII , the error introduced by an inaccurate M g estimation is not negligible. The error of E LII for the variation of M g with ±10% can reach up to ±5%. The relative error due to M g perturbation does not show strong dependency on the d pp in this range. Since the gas phase molecular weight is one of the base variables of numerical simulations, the novel forward approach can remove the uncertainty related to M g without any additional cost compared to the classical validation strategy.

Number of primary particles in the aggregate

As already discussed in Sec. 2.2, accounting for the morphology of the soot particles, i.e. the number of primary particles, drastically affects the d pp value obtained from the experimental LII signal. To additionally quantify the sensitivity of E LII to the number of primary particles in the aggregate, n p was varied from 1 to 100 (1, 5, 30, 60 and 100). In Figure 8, the E LII and its relative error from the case with n p = 1 is presented.

For all the considered diameters, E LII drastically changes moving from a single particle assumption to an aggregate structure, even when it is char- acterized by a small number of primary particles. The relative error reaches 25% already at n p = 10. The difference in E LII between n p = 60 and n p = 100 is not so relevant as the difference between n p = 1 and n p = 10. When the aggregate nature of the particle is not considered, a smaller This large source of errors can be avoided by reconstructing the LII signal from a numerical simulation that provides information about the aggregate nature of the particle, strengthening the validation of numerical modeling compared to the traditional inverse method.

PPSD shape

As mentioned before, the PPSD is not always lognormal. In order to see the effect on the indexes of a different PPSD shape, several tests are carried out. The first test corresponds to a 'bimodal' distribution similar to the one experimentally observed by Abid et al. [START_REF] Abid | On evolution of particle size distribution functions of incipient soot in premixed ethylene-oxygen-argon flames[END_REF][START_REF] Abid | Size distribution and morphology of nascent soot in premixed ethylene flames with and without benzene doping[END_REF].

It is imitated by a superposition of a first mode f 1 with a large number of very small particles below 5 nm and a second "moving mode" f 2 with a lognormal shape, characterized by a varying d 2 g =3 -20 nm and a constant As discussed previously, two indexes are used to characterize the LII signal emitted by a polydisperse distribution. In Figure 10, the E LII and σ LII are presented for the resultant bimodal distribution f 1 + f 2 and for the sole lognormal distribution f 2 . By comparing the two cases, the contribution of the first mode compared to a classical lognormal assumption can be evidenced. In both cases, the E LII increases with d 2 g , as expected from Fig. 3. The values related to the bimodal and the single lognormal PPSD are very similar, the deviance between the two values being below 1%. This means that the presence of a first mode in the PPSD only slightly affects the E LII results compared to a classical lognormal assumption as observed in [START_REF] Cenker | Determination of small soot particles in the presence of large ones from time-resolved laser-induced incandescence[END_REF]. Concerning σ LII , it can be observed that it decreases for increasing To investigate additionally the effect of the PPSD shape on the LII signal, a second test is performed. This time, the PPSD is obtained by superimposing two standard lognormal distributions: one (f 1 ) with fixed characteristics {d 1 g , σ 1 g } = {25 nm, 1.2 nm} and another (f 2 ) with a fixed

f 1 (d) =        0 : d < d m d α m d -α-1 : d ≥ d m (
σ 2 g = 1.
2 and a d 2 g varying from 10 to 60 nm ('moving mode'). Some examples of the distributions are presented in Figure 12. Results on E LII and σ LII for the combined PPSD f 1 + f 2 and for the 'moving mode' f 2 are plotted in Figure 13.

Concerning E LII , a quasi-linear behaviour is observed with d 2 g when considering only the 'moving mode' f 2 . When looking for the combined PPSD, it can be noticed that E LII is initially constant until d 2 g < 25 nm. Therefore, the obtained value is equal to the E LII value for a single lognormal function d g =25 nm. This is because the fixed lognormal mode f 1 has a constant d 1 g =25 nm and its contribution to the combined LII signal f 1 + f 2 is higher compared to the second moving mode

f 2 when d 2 g < d 1 g . Once d 2 g > d 1 g , it can be observed that E LII for f 1 + f 2 is increasing with d 2
g and almost collapsing with the results for f 2 . Concerning σ LII , the variations in Figure 13 are so small that they can be considered as negligible. This is consistent with the fact that both σ 1 g and σ 2 g are kept constant in this test case. From E LII and σ LII , it can be deduced that the LII signal is governed by the lognormal distribution with the higher d g values and that the presence of a second mode with a smaller d g cannot be easily identified from the LII signal.

Finally, a PPSD is obtained by superimposing two standard lognormal distributions: one (f 1 ) with fixed characteristics ({d 1 g , σ 1 g } = {40 nm, 1.2 nm}) and another (f 2 ) with a fixed d 2 g = 25 nm and σ 2 g varying from 1.05 to 1.6 nm ('moving mode'). Results are presented in Fig. 14. It can be observed that E LII for f 1 + f 2 is almost constant for σ 2 < 1.45 and higher than the result obtained when considering only f 2 . This proves that E LII for a combined PPSD is mainly governed by the mode with the highest d g (here f 1 ) for small σ 2 g . Once σ 2 g >> σ 1 g , the E LII is also affected by the presence of the first mode. Similarly, σ for f 1 + f 2 is mainly governed by the mode with the highest d g , i.e. f 1 , except for σ 2 g >> σ 1 g .

In conclusion, the mode with the highest d g mainly governs the LII signal at least as long as the other modes are not characterized by largely greater standard deviations, meaning by a large amount of big particles.

On the one side, this implies that it is extremely hard to correctly obtain the d pp from the experimental LII signal unless the PPSD is known. Therefore, the accuracy of the validation of the numerical modeling based on d pp , i.e. the inverse method, is strongly affected by the assumption of the shape of the PPSD. With the forward method, information on the PPSD are extracted from the simulation. On the other side, it means that by applying the forward method it is not possible to assess with certainty that the predicted PPSD agrees with the experimental ones, even if the synthesized and the experimental signals are in good agreement, since similar LII signals can be generated from different PPSDs. However, a disagreement on E LII and σ LII indicates that the predicted PPSD is not correct.

Discussion

All the above discussed uncertainties decrease the quality of the traditional validation method of the numerical PPSD by TiRe-LII. Recent numerical simulations provide access to most of the needed parameters. Therefore, the LII signal can be synthetized with a reduced number of assumptions. While temperature, pressure and molecular weight are basic output quantities of CFD simulations, depending on the soot model, the "Soot" quantities in Table 1 are available on a different level. For detailed modeling, information on PPSD and n p are classically available. Furthermore, once the relation of α T to other soot properties is known, this parameter can be also determined from the numerical calculation and related uncertainties can be avoided by the new approach. Utilizing these numerically derived quantities and the synthesized signal comparison approach, a more consistent validation of the whole modeling strategy can be performed. However, it has to be mentioned that as several PPSDs can lead to the same LII signal, the forward approach can not definitely confirm the validity of the PPSD, but it can be used to detect incorrect PPSDs and avoid the exclusion of a correct PPSD. Such issue is however intrinsic to TiRe-LII signal and affects the inverse method as well. In addition, general relations of α T and E(m) with other particles characteristics are today not available. Finally, it has to be reminded that the numerical synthesized reconstruction of the LII signal also depends of the quality of the gaseous phase description. However, also with the inverse technique, a correct prediction of the gas phase is a necessary step before validating the solid phase. In the following, the new validation approach is applied to the study of a premixed laminar sooting flame to prove the feasibility of the technique in the context of a real flame investigation.

Application to a premixed flame

In the following, both the traditional and the synthesized LII signal comparisons are performed for a target flame of the soot research community.

The objective here is not the validation of the numerical model used for the simulation since it has been already discussed in [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF] but to prove the feasibility of the forward approach. First, the mean diameters are calculated from the numerical PPSD obtained in [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF] and compared to measurement results [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF]. Then, the forward validation procedure is applied by reconstructing the synthesized LII signals from the numerical results at 4 different height above burners (HABs) and by comparing the E LII and σ LII indexes to the experimental data.

Flame configuration

The third laminar premixed ethylene target flame from the 4 th International Sooting Flame (ISF) Workshop [START_REF]International Sooting Flame (ISF) Workshop[END_REF] has been used here due to the availability of experimental LII signals in literature [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF]. In these works, a so-called McKenna burner is used to produce a premixed flat ethylene/air flame with an equivalence ratio of Φ = 2.1. The total flow rate is 10 l/min (at 0 • C and 1 atm) equivalent to 6.44 cm/s cold flow inlet velocity.

The LII measurements carried out by Bladh [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] were performed with a 1064 nm beam produced by a Q-switched Nd:YAG laser. The detection system was equipped with a bandpass filter at 575 nm with a 32 nm full width at half maximum (FWHM). The laser fluence was set to 0.13 J/cm 2 and the temporal profile of the laser pulse had a FWHM of ∼11 ns. This information, needed for the synthetization of the LII signal, has been extracted from [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF].

Numerical setup and chemical mechanism

The associated numerical simulations have already been performed and validated in [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF]. They were performed with the OpenSMOKE++ framework [START_REF] Cuoci | OpenSMOKE++: An object-oriented framework for the numerical modeling of reactive systems with detailed kinetic mechanisms[END_REF]. ing the equations for the primary particle number density of the sections with the method described in [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF]. Due to the simplicity of the flame, the errors introduced by an improper modeling of the flow field of the flow-flame interaction are expected to be negligible.

The inlet velocity and the species mass fractions were prescribed according to the experimental setup. The temperature profile was imposed as suggested by the ISF Workshop [START_REF]International Sooting Flame (ISF) Workshop[END_REF] in the reference calculation. This temperature profile differs from the one experimentally obtained by Bladh et al. [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] so that a numerical simulation was also performed using the experimental gas temperature in the [0.7,1.8] cm region. No relevant difference in terms of PPSD has been evidenced between the two calculations (not shown) so that in the following the validation is performed only with the results from the ISF temperature profile.

Classical inverse strategy: comparison of the mean primary particle diameter

The classical inverse validation of the model was performed by Bodor et al. [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF]. Only the main results from [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF] are summarised here to facilitate the reading. The LII measurements reported in Fig. 15 are based on an assumption of monodisperse population. Therefore, the monodisperse equivalent size recovered by Bladh et al. [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] is considered from the TEM measurements. For the same reason, the numerical PPSD results have been post-processed to obtain an equivalent mean diameter d eq mono [START_REF] Bodor | A post processing technique to predict primary particle size of sooting flames based on a chemical discrete sectional model: application to diluted coflow flames[END_REF]. The significant mismatch between the experimental results can be related to the different porous plug materials [START_REF] Migliorini | How "flat" is the rich premixed flame produced by your McKenna burner?[END_REF] used in the studies (a stainless steel porous plug was used [START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Stirn | Comparison of particle size measurements with laser-induced incandescence, mass spectroscopy, and scanning mobility particle sizing in a laminar premixed ethylene-air flame[END_REF], a bronze plug was used in [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF]).

In Figure 15, the comparison of the numerical and the experimental results [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Stirn | Comparison of particle size measurements with laser-induced incandescence, mass spectroscopy, and scanning mobility particle sizing in a laminar premixed ethylene-air flame[END_REF] is presented. Concerning the numerical results, d eq mono has been calculated with d det,min = 2 nm and d det,min = 5 nm (bottom and top blue dashed lines respectively) providing a range of possible solutions (brown area) since the diameter of the smallest experimentally captured particle d det,min is not indicated for all databases.

At low heights above burners (HABs), the numerical simulation tends to overpredict the d mono and a slower diameter increase is obtained compared to all experimental data.

It can be noticed that depending on the used experimental database, it is possible to conclude that the numerical model overpredicts the mean d pp in the early nucleating region or in the later surface growth region, leading then to different modifications of the model. The reevaluation with n p = 100 indicates that the difference may be also caused by the fact that the shielding effect was neglected when post-processing the experimental data.

Depending on the smallest particle d det,min considered when calculating the mean diameter from the numerical PPSD, the conclusions change again. In conclusion, the large error bars related to the uncertainty in the retrieved d pp prevents to determine the strategy for the numerical model development. In order to draw further conclusions on the validation and/or to identify potential paths of improvements, a more reliable evaluation is needed. In particular, to identify the reactions that may need an update, it is crucial to accurately locate the region of the d pp mismatch. Therefore, in the following section, the forward approach is tested on this configuration in order to understand if additional indications on the accuracy of the model can be obtained. The LII signal was reconstructed by an in-house code using the model equations reported in [START_REF] Hofmann | A web-based interface for modeling laser-induced incandescence (LIISim)[END_REF][START_REF] Liu | Influence of polydisperse distributions of both primary particle and aggregate sizes on soot temperature in low-fluence laser-induced incandescence[END_REF][START_REF] Liu | Effects of primary particle diameter and aggregate size distribution on the temperature of soot particles heated by pulsed lasers[END_REF]. The differential equation was solved with 4 th order Runge-Kutta method using the characteristics of the laser in the measurements of Bladh et al. [START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF]: laser wavelength = 1064 nm, fluence = 0.13 J/cm 2 , shot duration = 11 ns, bandpass filter = 575±16 nm. Vaporization is neglected since a low fluence has been considered.

While in the first study of Bladh et al. [START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF] a varying α T (within the range of 0.5-0.61) has been used, the second study [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] was performed with a constant value of α T (0.37). Most of the works on TiRe-LII signal evaluations are based on a constant value of α T [START_REF] Michelsen | Understanding and predicting the temporal response of laser-induced incandescence from carbonaceous particles[END_REF][START_REF] Schulz | Laser-induced incandescence: Recent trends and current questions[END_REF], though α T was shown to depend on the soot particle characteristics [START_REF] Bladh | Influence of spatial laser energy distribution on evaluated soot particle sizes using two-colour laser-induced incandescence in a flat premixed ethylene/air flame[END_REF][START_REF] López-Yglesias | Soot maturity and absorption cross sections[END_REF]. In the following, the constant value indicated by Bladh et al. [START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] is used here since no reliable estimation of α T evolution along the flame was possible.

The indexes of the synthesized and experimental LII signals are compared in Figure 16. The signal was reconstructed from the numerical PPSD experimental database, comparing synthesized and experimental LII signals is a more consistent approach to validate the numerical models.

Q f Z Q A r p G L v t p 9 8 n M h j V o + G z t R y w = " > A A A C 4 n i c j V F d S x t B F D 1 u 1 U b b a l o f R V g M h b 4 0 b B q p + h b 0 p Q E f I h g T S N I w u 5 m k Q / a L m d l i C H n y r W / F V / 9 A X / X H S P 9 B / R f e G V e x B G n v s r t n z r 3 n z N y 5 f h o K p T 3 v 9 4 L z Y n F p + W V h Z f X V 6 z d r 6 8 W 3 7 0 5 V k s m A N 4 M k T G T b Z 4 q H I u Z N L X T I 2 6 n k L P J D 3 v L H h y b f + s 6 l E k l 8 o i c p 7 0 V s F I u h C J g m q l / c 6 i o x i t j X r t J M 9 q d d z c / 0 9 K h e n 8 3 c z s d e v 1 j y y p 4 N d x 5 U c l B C H o 2 k e I M u B k g Q I E M E j h i a c A g G R U 8 H F X h I i e t h S p w k J G y e Y 4 Z V 0 m Z U x a m C E T u m 7 4 h W n Z y N a W 0 8 l V U H t E t I r y S l i / e k S a h O E j a 7 u T a f W W f D P u c 9 t Z 7 m b B P 6 + 7 l X R K z G N 2 L / p X u o / F + d 6 U V j i D 3 b g 6 C e U s u Y 7 o L c J b O 3 Y k 7 u P u l K k 0 N K n M E D y k v C g V U + 3 L N r N c r 2 b u 6 W 2 f w f W 2 l Y s w 7 y 2 g y 3 5 p R 2 w P s m P j + O c x 6 c f i p X q u X q 8 U 6 p d p C P u o B N b O M D z X M X N X x B A 0 3 y P s c v X O H a G T g / n J / O x X 2 p s 5 B r N v B X O J d 3 Y D a b u Q = = < / l a t
For the interpretation of the results, it should be reminded that only the first 800 ns have been considered to calculate the indexes. Therefore, the tendencies should be interpreted by looking at Figs. 3(d) and 4(d). Results on d mono are also summarized in Table 2. The tendencies of the indexes are qualitatively similar between the experimental and the numerical results. The index E LII increases up to 13 mm. This increase slows down for HAB>13 mm. Concerning σ LII , it decreases with the HAB almost reaching a plateau for HAB ≥ 13mm. This behaviour indicates that d g and/or σ g increase with the height above the burner. When compared to the experimental data, the numerical simulation over-predicts E LII (and underestimates σ LII ) for HAB=7 mm, whereas it under-predicts E LII (and overestimates σ LII ) for HAB=13 mm, with a quite good agreement for HAB=9 mm. This indicates that d g and/or σ g are initially overestimated by the simulation, whereas they are under-predicted downstream the flame. These results obtained with the forward approach confirm the results on d mono obtained with the inverse technique in Table 2. Therefore, the forward analysis can be used to confirm the conclusions obtained with the inverse approach for d det,min = 2nm.

However, an additional indication can be deduced by looking at the results with the forward approach. In Table 2 on d mono is around ±30% for the different available points. On the contrary, when looking to σ LII a significant disagreement between the experimental and the numerical data is observed at HAB=7 mm. This indicates that the error on the PPSD prediction, presented in Fig. 17, is probably more significant at this height compared to the other positions. One possible scenario to explain the fact that the σ LII value is strongly underestimated is that σ g is highly overestimated by the numerical results. This can possibly indicate that in the numerical PPSD of Fig. 17 the appearance of the second mode (composed of big particles largely contributing to the LII signal) is predicted too quickly. The forward procedure seems to indicate that the transition between unimodal to bimodal PPSD, often obtained as a results of collisional processes, is predicted too quickly. Therefore, a possible path of improvements of the model may consist in an accurate reduction of the collisional processes intensity, leading to a slower transition from a uni-modal to bimodal PPSD, i.e. to higher σ g values, and the verification of its effect on the synthesized LII signal. Additionally, it can be concluded that a longer detection time for experimental LII signals should be used to identify more easily if the discrepancies are related to d g or σ g misprediction using results

from Figs. 3(a) and 4(a).

Conclusions

In this work, a novel forward approach has been proposed to validate the numerical models for the primary particle size distribution prediction with experimental results obtained with the TiRe-LII technique. It is based on the numerical reconstruction of the temporal evolution of the LII signals from the simulation information. The difficulties of the traditional inverse strategy of validation were pointed out. They are mainly related to the uncertainties of the parameters governing the LII signal such as the gas temperature, the pressure, the molecular weight, the PPSD shape and the number of particles in the aggregate.

On the contrary, comparing directly the numerically synthesized and the experimental LII signals is very attractive since many of these quantities are generally available from detailed numerical simulations. In order to compare the experimental and the synthesized LII signals, two characteristic indexes have been suggested, the expected value E LII and the non-dimensional standard deviation σ LII . Thanks to these indexes, the error possibly introduced by the exponential fitting procedure frequently used for the TiRe-LII signal evaluations can be avoided.

The feasibility of the forward method was tested on a premixed ethylene ISF target flame [START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF] together with the inverse approach. The traditional inverse comparison is characterized by a significant experimental error bars and a strong effect of the definitions of the numerical mean d pp and the experimental d det,min so that no conclusion can be definitely assessed. On the contrary, the forward method showed that larger primary particles and/or a larger standard deviation of the PPSD should be predicted in the upper flame region. These conclusions are consistent with the results on d mono obtained with the inverse method but they are less affected by the discussed uncertainties. Therefore, the novel forward approach can be considered as a consistent validation strategy complementary to the inverse method. Even if the forward strategy allows to reduce the uncertainties, it should be reminded that similar LII signal decay corresponds to various possible PPSDs, so that the agreement between synthesized and measured signal does not ensure the validity of the PPSD. This issue is not due to the post-processing validation strategy but it is intrinsic to TiRe-LII measurements. However, when a mismatch is observed in such a comparison, it can be used to state that the predicted PPSD is incorrect and add further indications for the modeling development. In addition, it has to be reminded that the numerical synthesized signal depends on the model's accuracy for both the gas and the solid phases, allowing to investigate the quality of the whole numerical strategy and not only of the solid phase simulation.

The potential improvement of the new approach is related to the thermal accommodation coefficient, which is one of the sources of uncertainty even in the novel approach. The dependency of α T on other soot properties is not clearly established yet. However, once the relation will be found, the reconstruction of the thermal accommodation coefficient could be done providing accurate predictions even in the presence of a nascent and mature soot particle mixture.

sible d mono values retrieved from LII signals constructed from a lognormal

  the diameter d mono derived by assuming monodispersity from the LII signal does not coincide with the monodisperse equivalent mean particle diameter, d eq mono of the lognormal distribution. Franzelli et al. [31] evaluated the measured LII signals of a series of co-flow laminar flames with both monodisperse and lognormal PPSD assumption. In this work, beside d mono , the arithmetic mean primary particle diameter defined as d amean = exp ( ln(d g ) + ( ln(σ g ) 2 ) / 2 ) was presented. While d amean varied in the range of 10 to 65 nm, d mono reached up to ∼ 90 nm, meaning a difference as high as 30% for the heavily sooting flame.

EFig. 3

 3 Fig. 3 Expected value E as function of {dg, σg} pairs. Results for (a) the original signal ELII , (b) the single-exponential fitted curve ESExp, (c) the doubleexponential fitted curve EDExp and (d) the original signal considering only up to 800 ns.

Fig. 4

 4 Fig. 4 Non-dimensional standard deviation as function of {dg, σg} pairs. Results for (a) the original signal σ LII , (b) the double-exponential fitted curve σ DExp , (c) the original signal zoomed on small dg values and (d) the original signal considering only up to 800 ns.

Fig. 5

 5 Fig. 5 ELII (triangles) and its relative error (circles) as a function of the temperature for four dpp=5, 15, 25 and 35 nm (red, green, blue and yellow).

Fig. 6

 6 Fig. 6 ELII (triangles) and its relative error (circles) as a function of pressure for four dpp=5, 15, 25 and 35 nm (red, green, blue and yellow).

  weight [kg/kmol]

Fig. 7

 7 Fig. 7 ELII (triangles) and its relative error (circles) as a function of the molecular weight for four dpp=5, 15, 25 and 35 nm (red, green, blue and yellow).

Fig. 8

 8 Fig. 8 ELII (triangles) and its relative error (circles) as a function of np for four dpp=5, 15, 25 and 35 nm (red, green, blue and yellow).

E

  LII is assigned to the same d pp , in agreement with the lower decay time observed for the similar tests performed by Liu et al.[START_REF] Liu | Effects of primary particle diameter and aggregate size distribution on the temperature of soot particles heated by pulsed lasers[END_REF]. This implies an over-prediction of the d pp value since for E LII (for example E LII = 160 ns) a larger diameter (d pp =35 nm) is estimated for a spherical particle (n p = 1) compared to big aggregates (for n p = 100, d pp = 25 nm).

Fig. 9

 9 Fig. 9 First mode f1 and examples of "moving mode" f2 with various d 2 g (indicated in the label)

7 )

 7 with α = 3 and d m = 2.1mm. The first mode and a few examples of the second mode are presented in Figure 9. The distribution is split up to 0.5 nm wide sections.

Fig. 10 ELII

 10 Fig. 10 ELII and σ LII of the LII signal originating from the bimodal distributions f1 + f2 (black circle) and from the moving lognormal mode f2 only (white square). The moving mode has varying d 2 g and constant σ 2 g .

d 2 g

 2 values in analogy with results from Figure 4. The difference between the two signals remains small, being below 1% for this second index, despite the very different PPSDs, indicating an almost negligible effect of the first mode on the LII signal, meaning that very similar LII signals are obtained with the two distributions f 2 and f 1 + f 2 . In Fig. 11, results of the mean d pp for the bimodal and the lognormal distributions are presented. Since the LII signals are similar, it can be deduced that when investigating a bimodal distribution by assuming a log-normal shape, the d pp value extracted from the LII signal overstimates the real d pp . On the contrary, when comparing the synthetized and the experimental signal, conclusions are not biased by the assumption of the PPSD shape.

Fig. 11 Fig. 12

 1112 Fig. 11 Primary particle diameter dpp for the bimodal distributions f1 + f2 (circle) and for the moving lognormal mode f2 only (line) as function of d 2 g .

Fig. 13 ELII

 13 Fig.13ELII and σ LII of the LII signal originating from two superimposed lognormal distributions f1 +f2 (solid line) and from the moving mode f2 only (dashed line). The moving mode has varying d 2 g and constant σ 2 g .

Fig. 14 ELII

 14 Fig. 14 ELII and σ LII of the LII signal originating from two superimposed lognormal distributions f1 +f2 (solid line) and from the moving mode f1 only (dashed line). The moving mode has constant d 1 g and varying σ 1 g .

Fig. 15

 15 Fig.15Experimentally[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF][START_REF] Axelsson | Laser-induced incandescence for soot particle size measurements in premixed flat flames[END_REF][START_REF] Stirn | Comparison of particle size measurements with laser-induced incandescence, mass spectroscopy, and scanning mobility particle sizing in a laminar premixed ethylene-air flame[END_REF] and numerically obtained primary particle diameter with monodisperse assumption. "np = 100" marks the reevaluation performed at HAB = 1.7 cm by Bladh et al.[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] with aggregate structure assumption

4. 4

 4 Novel forward strategy: comparison of TiRe-LII signalsIn this section, the measured and the synthesized LII signals are compared in terms of E LII and σ LII to reduce the uncertainties of the post-processing technique and potentially strengthen the validation procedure. By doing so, some potential weaknesses of the numerical model can be investigated additionally to the classical inverse strategy. Four HABs (7, 9, 13 and 17 mm) were selected for the LII signal reconstruction and to perform a comparison with the experimental data[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF].

Fig. 16

 16 Fig.[START_REF] Schenk | Imaging Nanocarbon Materials: Soot Particles in Flames are Not Structurally Homogeneous[END_REF] Comparison of synthesized and measured[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] LII signal indexes at four HABs(7, 9, 13 and 17 mm)

Fig. 17

 17 Fig.17Numerical PPSD for the premixed laminar flame experimentally investigated by[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF].

  Fig.17Numerical PPSD for the premixed laminar flame experimentally investigated by[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF].

  SExp of the single exponential fit, which is also shown in Fig.2. It can be observed that the decay time τ SExp and the expected value E SExp of the curve constructed from the single-exponential fit overlap, as expected.However, even for the monodisperse case, the single-exponentially fitted curve does not capture the expected value of the original curve E LII for large diameters, showing the inaccuracy of the single-fit procedure. On the contrary, the expected value of the double-exponential fit E DExp shows a good agreement with E LII at least for a monodisperse distribution. For a monodisperse variation, a quasi-linear relation between E LII and d mono is observed meaning that an error observed on E LII directly implies the same error on d mono . This provides a first indication to interpret the agreement

		300
			E LII
	Characteristic value [ns]	250 200 150 100	E DExp E SExp τ SExp
		50
		10	20	30	40	50	60
				d mono [nm]
							6, the error
	remains below 5%. The error for the single-exponential fit varies between

Results for E LII , E SExp and E DExp for a monodisperse distribution are presented for different values of d mono in Fig. 2. It can be observed that E LII depends almost linearly on d mono , as

shown by the dotted line in Fig.

2

. Therefore, a single index is sufficient to characterize the signal in case of a monodisperse distribution. Supposing a perfectly exponential decay, E LII is expected to be equal to the decay time τ between experimental and numerically synthetized signals. For example, assuming a monodispersed population a disagreement of 20% on E LII between the experimental and numerical signals implies an error of 20% on the predicted d mono . For a polydisperse population, this interpretation is more complex requiring the analysis at least of 2 indexes as discussed in the following.

3.1.2 Polydisperse distribution

In Figure

3

, E LII , E SExp , and E DExp are presented for a lognormal PPSD with different {d g , σ g } pairs. A similar behaviour of E is observed for the three cases. First, it can be observed that a given value of E corresponds to multiple {d g , σ g } couples so that a second index is needed to discriminate the PPSD, as expected. The recovered curves both from single and double-exponential fits show a similar behavior. The double-exponential fit provides a good agreement for σ g < 1.

  The kinetic mechanism includes 189 gaseous species and a chemical discrete sectional soot model (CDSM) with 20 sections (BINs), which are divided into subsections based on the H/C ratio. While heavy polycyclic aromatic hydrocarbons (PAHs) (from BIN 1 to BIN 4 ) and small soot particles (from BIN 5 to BIN 10 ) have 3 H/C subsections, large soot particles (BIN 11 and BIN 12 ) and soot aggregates (from BIN 13 to BIN 20 ) are divided into 2 H/C subsections. Information on the PPSD is obtained by transport-

  , it can be observed that the error HAB [mm] Bladh et al.[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] Numerical data Error[-]

	7	5.8	8.23	+30%
	9	12.6	11.4	-10%
	13	23	15.1	-34%
	17	28	17.1	-39%

Table 2

 2 dmono [nm] at different HABs derived from TiRe-LII by Bladh et al.[START_REF] Bladh | Optical soot characterization using two-color laser-induced incandescence (2C-LII) in the soot growth region of a premixed flat flame[END_REF] and from a numerical simulation with dpp,min = 2 nm. The corresponding errors are indicated in the last columns.

As an example, the CVT parameter is 0.658 when comparing a lognormal distribution characterized by a geometrical mean dg and a standard deviation σg {dg,σg} = {30 nm, 1.25 nm} and a lognormal distribution with {dg,σg} = {26 nm, 1.23 nm}[START_REF] Daun | Determining aerosol particle size distributions using time-resolved laser-induced incandescence[END_REF].

In an aggregate, the inner particles are thermally shielded by the outer particles from the cold ambient gas molecules so that they cool down more slowly.

Acknowledgements

This project has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska-Curie Grant Agreement No 643134. B. Franzelli acknowledges the support of the European Research Council (ERC) under the European Unions' Horizon 2020 research and innovation program (Grant Agreement No. 757912).