Is ortho -Cresol a Viable Lignocellulosic Blendstock? A Kinetic Study of Its Co-Oxidation within a Surrogate Fuel

: The viability of the use of ortho -cresol as a bio-blendstock or antiknock additive from lignocellulosic biomass is assessed; Ignition delays of ortho -cresol within blends with iso -octane are measured with the ULille rapid compression machine, and compared with results from the literature; It is shown that ortho -cresol has a strong inhibiting effect on the reactivity towards ignition, most notably in the Negative Temperature Coefﬁcient region; This effect is found to originate from competition with iso -octane on the OH radicals, where the reactivity of ortho -cresol with these radicals does not lead to radical chain-branching.


Introduction
The quest for viable and affordable bio-blendstocks has led researchers to explore the potential of lignin-based compounds for the production of alternative fuels. Such alternative fuels have strong requirements for compliance with existing spark-ignition engine technologies, among which sufficient resistance to engine knock. Knock is caused by unwanted auto-ignition of the fuel/air mixture and is commonly evaluated by octane number measurements. Lignin is a cheap source of species in the structure of which the presence of an aromatic ring is often viewed as a guarantee of strong anti-knock behavior. Among such species, ortho-cresol has been identified as a promising additive in a kinetic modeling study aiming at comparing the anti-knock tendencies of various oxygenated aromatics [1]. This effect was discussed and suggested to originate from the relative stability of 2-methylene, cyclohexa-3,5-dien-1-one, a species that can be formed easily from the two most expected radicals in the case of ortho-cresol. These two radicals, later referred to as R1 and R2, are presented in Figure 1 along with the structure of ortho-cresol. Both of them are stabilized by resonance, and can therefore be long-lived in reactive systems, and difficultly undergo chemical chain-branching. It is also possible that such radicals recombine with the HO 2 radical, thereby forming an unstable hydroperoxide, as observed for alkenes [2,3] or the benzyl radical [4]. These pathways have however not been considered in [1]. Blending experiments within a diesel fuel have also demonstrated a reduction in NO x and soot emissions from a compression-ignition engine when ortho-cresol is added [5].
Energies 2021, 14, x FOR PEER REVIEW 2 of 10 in order to include all the necessary pathways into kinetic models, as demonstrated in the past [4,[10][11][12]. To serve that purpose, because of its two-stage ignition features and relevance to spark-ignition engine fuels, iso-octane was used as a co-oxidant in this study. Therefore, the present study aims at evaluating the potential of ortho-cresol as a bio-blendstock, and gaining insight into the low-temperature oxidation pathways of ortho-cresol at high pressures, through an experimental study of the ignition of isooctane/ortho-cresol blends inside a rapid compression machine (RCM). The obtained results ortho-cresol is also a very common product of the oxidation of aromatics [4,6] in combustion applications, as well as of the pyrolysis of biomass [7][8][9]. There is however a lack of available literature on its combustion kinetics. Moreover, in practical conditions, ortho-cresol can be in presence of a radical pool that may not exist if it had been on its own. Especially, in the case of automotive fuels, species that are often considered unreactive can undergo co-oxidation because of the presence of a more reactive compound inside the fuel blend. Studying co-oxidation conditions is therefore of crucial importance in order to include all the necessary pathways into kinetic models, as demonstrated in the past [4,[10][11][12]. To serve that purpose, because of its two-stage ignition features and relevance to sparkignition engine fuels, iso-octane was used as a co-oxidant in this study.
Therefore, the present study aims at evaluating the potential of ortho-cresol as a bioblendstock, and gaining insight into the low-temperature oxidation pathways of ortho-cresol at high pressures, through an experimental study of the ignition of isooctane/ortho-cresol blends inside a rapid compression machine (RCM). The obtained results will be compared with those from a model available in the literature, after which this model will be used to gain insight into the results and the effect of ortho-cresol on the ignition delays of iso-octane.
To the authors' knowledge, this constitutes the first experimental study of the combustion kinetics of ortho-cresol.

Experimental Methods
The ULille RCM was used to measure the ignition delays of stoichiometric isooctane/ortho-cresol blends at several blending ratios between 0 and 40% mol. ortho-cresol in the fuel blend, and at pressures between 14 and 20 bar, and temperatures between 688 and 850 K. Because this experimental setup has been described extensively in the literature before [11][12][13], a schematic being given in [14], only the features relevant to this work will be described here. Mixtures of iso-octane, ortho-cresol, dioxygen, and inert gases were prepared using the partial pressure method inside a heated mixture preparation facility, whose temperature was fixed as superior to the one necessary for a vapor pressure of all compounds superior to twice the required one. All lines leading to the RCM were heated to a temperature superior by 5 K to this temperature to prevent any condensation. The temperature at top dead center was calculated from the adiabatic core assumption and varied by using a mixture of diluents including CO 2 , N 2 , and Ar. The composition of the mixtures used in this study is given in the Supplementary Material.
First-Stage (FSIDT) and total Ignition Delay Times (IDT) are defined as the time lapse between the Top Dead Center (TDC) and their respective maximum in pressure derivative. A creviced piston and a reproducible compression time of 45 ms were used to prevent the formation of a piston corner vortex during compression. During the modeling of the results, volume profiles processed from experiments where O 2 was replaced by N 2 were used as input to model the compression phase and the heat losses after compression. These volume profiles are also provided as Supplementary Material. An uncertainty of ±5 K and ±0.1 bar is assumed on the reported temperatures and pressures, respectively. All of the measured IDTs are reported in the figures, demonstrating the reproducibility of the measurement.

Kinetic Modeling
All simulations were performed using the Cantera solver [15]. The model used in this work is a kinetic model for the oxidation of anisole-isooctane blends [11]. The submechanism describing the oxidation of iso-octane originates from the work by Fang et al. [16]. The anisole submechanism was updated from a model by Wagnon et al. [17] on the basis of new ab initio calculations. It contains a submechanism for the combustion kinetics of cresol sharing large similarities with the submechanism proposed by Bounaceur et al. [18]. To the authors' knowledge, detailed kinetic models dealing with the oxidation of cresols are scarce in the literature and have not been the subject of dedicated studies. In the early 1990s, Emdee et al. [19] observed the formation of cresols during the oxidation of toluene in a flow reactor. They explained their kinetics of formation by an electrophilic addition of an O atom on the aromatic ring by toluene + O = cresoxy + H or the addition of a methyl radical on phenol giving cresol and an H-atom. The cresoxy radical can further recombine with an H atom to form cresol. Lindstedt et al. [20] refined this model by considering the consumption of cresols via H loss and H-atom abstraction by H and OH at the hydroxyl position followed by the decomposition of the cresoxy radical yielding benzene, CO, and an H-atom. The mole fraction profiles of toluene oxidation products were measured in a JSR (Jet-Stirred Reactor) and the ignition delay times of toluene/oxygen/argon mixtures were measured in a shock tube [18]. The authors developed a detailed kinetic model reproducing their experimental data. This model contains primary and secondary mechanisms of the oxidation of toluene, and among them, a sub-mechanism for the oxidation of cresols and derived radicals. The H-atom abstraction reactions forming two different radicals were considered. The reaction rate parameters for the formation of the methylphenoxyl radical (I < 3500 wonOC 6 H 4 CH 3 ) were taken equal to phenol [21] and for the formation of hydroxylbenzyl radical (HOC 6 H 4 CH 2 ) were taken equal to toluene [18]. The consumption of methylphenoxyl radicals happens via CO elimination to form the methycyclopentadienyl radical, when the consumption of hydroxylbenzyl radicals occurs by additions to molecular oxygen or to the hydroperoxyl radical. These two last reactions mainly lead to the formation of methylbenzaldehyde. It can be noted that none of these models consider the ortho-, meta-, or parastructure of the cresol isomers, which are lumped into one species.

Ignition Delays
The effect of the addition of ortho-cresol on the ignition delays of iso-octane/orthocresol blends is displayed in Figure 2 at two different compressed temperatures T C = 686 and 850 K. At T C = 686 K, the mixture displays two-stage ignition, representative of the transition from low-temperature to intermediate-temperature combustion. As the orthocresol fraction increases, both first-stage and total ignition delay increase monotonically. Both ignition delays are simulated with an acceptable agreement in this condition, with a maximum deviation of about 30%. This discrepancy is nearly constant, and also observable in the pure iso-octane case. This suggests that in these conditions, the effect of ortho-cresol addition is well captured by the model. However, at T C = 850 K, the experiments show little to no effect of ortho-cresol addition on the ignition delays, while the model predicts a large increase in the IDT. Note that no first-stage ignition was observed at this higher temperature, in agreement with similar studies with iso-octane [14]. This is consistent with the reduced performance of the model in this temperature range in the case of iso-octane/anisole blends, as observed in [11].
These conclusions can be generalized over a wide range of pressures, as demonstrated in Figure 3. The discrepancy between experimental and simulated ignition delays is consistent all over the investigated pressure range for pure iso-octane mixtures, with an underestimation of the IDT of about 30% at T C = 850 K, and a similar overestimation of both IDTs at T C = 686 K. While this remains true at all investigated ortho-cresol mole fractions for the lowest temperature, where two-stage ignition is observed, the model again overestimates the effect of ortho-cresol addition on the IDT at all pressures. Experimentally, the addition of ortho-cresol hardly modifies the IDT at T C = 850 K.
To gain further insight, an 80/20% mol. iso-octane/ortho-cresol mixture was further studied. A comparison of the evolution of its IDTs as a function of temperature with a pure iso-octane case is pictured in Figure 4. Both cases demonstrate Negative Temperature Coefficient (NTC) behavior, but the temperature domain at which this behavior is observed is however shorter for the mixture case in comparison to pure iso-octane. Ortho-cresol addition results in an increase in the IDT at all investigated temperatures, this effect being maximal in the NTC zone. Around 750 K, the IDTs are multiplied by about a factor of 6 by blending 20% of ortho-cresol into iso-octane. The temperature range at which first-stage ignition could be observed is also shorter than in the case of pure iso-octane, for which first-stage ignition is still observed for T C above 750 K. However, at the highest investigated temperature T C = 850 K, the effect of ortho-cresol addition is nearly negligible, as discussed above. The model captures the increase in the IDT and reproduces the experimental first-stage and total IDTs with a fair agreement for the pure iso-octane and blend cases, especially considering that it was not validated on experimental data on ortho-cresol before, and was not modified as part of this study. Additional work on the accurate reproduction of the ignition chemistry of iso-octane in the NTC temperature range could help improve these results. These conclusions can be generalized over a wide range of pressures, as demo strated in Figure 3. The discrepancy between experimental and simulated ignition delay is consistent all over the investigated pressure range for pure iso-octane mixtures, with a underestimation of the IDT of about 30% at TC = 850 K, and a similar overestimation both IDTs at TC = 686 K. While this remains true at all investigated ortho-cresol mo fractions for the lowest temperature, where two-stage ignition is observed, the mod again overestimates the effect of ortho-cresol addition on the IDT at all pressures. E perimentally, the addition of ortho-cresol hardly modifies the IDT at TC = 850 K.   These conclusions can be generalized over a wide range of pressures, as demonstrated in Figure 3. The discrepancy between experimental and simulated ignition delays is consistent all over the investigated pressure range for pure iso-octane mixtures, with an underestimation of the IDT of about 30% at TC = 850 K, and a similar overestimation of both IDTs at TC = 686 K. While this remains true at all investigated ortho-cresol mole fractions for the lowest temperature, where two-stage ignition is observed, the model again overestimates the effect of ortho-cresol addition on the IDT at all pressures. Experimentally, the addition of ortho-cresol hardly modifies the IDT at TC = 850 K.  To gain further insight, an 80/20% mol. iso-octane/ortho-cresol mixture was further studied. A comparison of the evolution of its IDTs as a function of temperature with a reproduces the experimental first-stage and total IDTs with a fair agreement for the pu iso-octane and blend cases, especially considering that it was not validated on expe mental data on ortho-cresol before, and was not modified as part of this study. Additio work on the accurate reproduction of the ignition chemistry of iso-octane in the N temperature range could help improve these results.

Kinetic Analysis
Because this study focuses on the co-oxidation of ortho-cresol in conditions wher will not show auto-ignition, the following discussion will examine further the TC = 686 condition, where two-stage ignition takes place. A brute-force sensitivity analysis on first-stage and total ignition delays has been performed and is shown in Figure 5. Becau of the size of the mechanism and the focus on explaining the effect of ortho-cresol, t sensitivity analysis was performed only taking into account the submechanism for t species, as well as radicals R1 and R2, phenol, and the phenoxy radical. First, one c note that the two hydrogen abstraction reactions by OH leading to R1 and R2 have highest sensitivity coefficient for both ignition delays. The fact that both have posit sensitivity coefficients illustrates that ortho-cresol acts mainly as a radical scavenger OH, thereby reducing the reactivity in the low-temperature zone. This is also true for and total (full symbols, solid lines) ignition delays for stoichiometric pure isooctane/'air' (black) and 80/20% mol. iso-octane-ortho-cresol/'air' (red) mixtures with the core gas temperature T C . P C = 20 bar.

Kinetic Analysis
Because this study focuses on the co-oxidation of ortho-cresol in conditions where it will not show auto-ignition, the following discussion will examine further the T C = 686 K condition, where two-stage ignition takes place. A brute-force sensitivity analysis on the first-stage and total ignition delays has been performed and is shown in Figure 5. Because of the size of the mechanism and the focus on explaining the effect of ortho-cresol, this sensitivity analysis was performed only taking into account the submechanism for this species, as well as radicals R1 and R2, phenol, and the phenoxy radical. First, one can note that the two hydrogen abstraction reactions by OH leading to R1 and R2 have the highest sensitivity coefficient for both ignition delays. The fact that both have positive sensitivity coefficients illustrates that ortho-cresol acts mainly as a radical scavenger of OH, thereby reducing the reactivity in the low-temperature zone. This is also true for the hydrogen abstraction reaction by OH on ortho-hydroxybenzaldehyde, demonstrating the importance of this secondary species, which will be confirmed below. Moreover, the formation of R2 inhibits the reactivity more than the corresponding reaction forming R1. This can be explained by the fact that according to the simulations, R2 is mainly recycled into cresol, when R1 can undergo further reactivity with HO 2 and O 2 . This creates a reactivity loop where R1 then reacts with HO 2 to yield back ortho-cresol, and this reaction therefore shows a positive sensitivity coefficient on the first-stage ignition delay. This conclusion can however be toned down by the fact that further reaction pathways for radical R2, such as the ones yielding cyclohexa-3,5-dien-1-one described in [1], are not present in the kinetic model. One can also notice the negative sensitivity coefficient of the HO 2 addition reaction on radical R1. This results in the formation of a hydroperoxide which can readily decompose into a hydroxybenzoxy radical and an OH radical, therefore facilitating ignition by converting a relatively unreactive HO 2 radical into more reactive radicals. The hydrogen abstraction reaction by HO 2 that yields R1 however has a negative sensitivity coefficient on the total ignition delay while its sensitivity coefficient for the first-stage ignition delay is nearly equal to zero, showing that during the second stage of ignition (or in the intermediate temperature range) ortho-cresol can reduce the ignition Energies 2021, 14, 7105 6 of 9 delay by facilitating the conversion of HO 2 into H 2 O 2 , whose decomposition into two OH radicals is commonly associated with the final ignition event.
which can readily decompose into a hydroxybenzoxy radical and an OH radical, therefore facilitating ignition by converting a relatively unreactive HO2 radical into more reactive radicals. The hydrogen abstraction reaction by HO2 that yields R1 however has a negative sensitivity coefficient on the total ignition delay while its sensitivity coefficient for the first-stage ignition delay is nearly equal to zero, showing that during the second stage of ignition (or in the intermediate temperature range) ortho-cresol can reduce the ignition delay by facilitating the conversion of HO2 into H2O2, whose decomposition into two OH radicals is commonly associated with the final ignition event. A reaction pathway analysis was also performed to shed light on the main conversion channels of ortho-cresol and its products at two stages of the ignition delay of the 80/20 iso-octane/ortho-cresol blend at TC = 686 K. The first one corresponds to 5% conversion of ortho-cresol and the first-stage ignition chemistry, while the second one corresponds to 75% conversion of ortho-cresol and the chemistry corresponding to the final ignition stage. The results of this analysis are shown in Figure 6. They show that during the first stage of ignition, R2 is formed preferentially by hydrogen abstraction by OH. However, since none of its conversion pathways are active in these conditions, it easily converts back to ortho-cresol, mostly by converting an HO2 radical into O2. This loop is effective in reducing the reactivity of the mixture towards first-stage ignition, by converting two radicals into stable species, but could result from missing reaction pathways, Figure 5. Brute-force sensitivity analysis on the first-stage and total ignition delays of a stoichiometric 80/20 iso-octane/ortho-cresol mixture. T C = 686 K, P C = 20 bar.
A reaction pathway analysis was also performed to shed light on the main conversion channels of ortho-cresol and its products at two stages of the ignition delay of the 80/20 iso-octane/ortho-cresol blend at T C = 686 K. The first one corresponds to 5% conversion of ortho-cresol and the first-stage ignition chemistry, while the second one corresponds to 75% conversion of ortho-cresol and the chemistry corresponding to the final ignition stage. The results of this analysis are shown in Figure 6. They show that during the first stage of ignition, R2 is formed preferentially by hydrogen abstraction by OH. However, since none of its conversion pathways are active in these conditions, it easily converts back to ortho-cresol, mostly by converting an HO 2 radical into O 2 . This loop is effective in reducing the reactivity of the mixture towards first-stage ignition, by converting two radicals into stable species, but could result from missing reaction pathways, as mentioned above. When R1 is formed, it mainly reacts by recombining with HO 2 , thereby forming a hydroperoxide that decomposes into the ortho-hydroxyl-benzoxy radical, which will in turn form ortho-hydroxybenzaldehyde. This compound is also formed minorly by the addition of R1 to O 2 .
During the ignition stage, the conversion of R1 still leads mostly to hydroxy-benzaldehyde. However, the balance between the direct channel of the addition of R1 to HO 2 and the stabilization of the hydroperoxide is shifted towards the former pathway. At all temperatures, ortho-hydroxybenzaldehyde reacts through hydrogen abstraction on the carbonylic carbon atom site, and the thereby formed radical undergoes the elimination of carbon monoxide to yield the phenoxy radical. The fate of the phenoxy radical lies in losing a hydrogen atom to form cyclohexadienone isomers, or to abstract a hydrogen atom to yield phenol. At higher temperatures, an alternate pathway leading to CO and the cyclopentadienyl radical is also possible.
Energies 2021, 14, 7105 7 of 9 R1 to HO2 and the stabilization of the hydroperoxide is shifted towards the former pathway. At all temperatures, ortho-hydroxybenzaldehyde reacts through hydrogen abstraction on the carbonylic carbon atom site, and the thereby formed radical undergoes the elimination of carbon monoxide to yield the phenoxy radical. The fate of the phenoxy radical lies in losing a hydrogen atom to form cyclohexadienone isomers, or to abstract a hydrogen atom to yield phenol. At higher temperatures, an alternate pathway leading to CO and the cyclopentadienyl radical is also possible. Figure 6. Rate of production analysis for the submechanism of ortho-cresol in the case of a stoichiometric 80/20 iso-octane/ortho-cresol blend at TC = 686 K. PC = 20 bar.Italic numbers correspond to 5% ortho-cresol conversion (i.e., first-stage ignition), bold numbers to 75% ortho-cresol conversion (i.e., final ignition).

Discussion
Insight on the inhibiting effect of ortho-cresol can be gained from the kinetic analysis. The inhibiting effect originates from a radical scavenging effect of ortho-cresol, which is mainly due to the implication of OH radicals into hydrogen abstraction reactions. OH radicals are typically formed in this case by the active low-temperature branching originating from iso-octane which induces first-stage ignition. This scavenging effect results in reduced heat release in the first-stage ignition, thereby limiting the temperature increase and increasing the duration of the second stage of ignition. This results in comparable first-stage ignition delays for blends with up to 30% ortho-cresol, as demonstrated in Figure 2, while the total ignition delay increases monotonically as the ortho-cresol content increases. When R2 is formed in the first stage of ignition, it does not lead to radical chain-branching since its most important conversion channels require higher temperatures. R1 can undergo further reactivity, but then yields phenol, cyclohexadienones, or ortho-hydroxybenzaldehyde, i.e., species which do not contribute to the formation of radicals in significant quantities.
A model developed to describe the co-oxidation of anisole, with an updated phenol and phenoxy submechanism and an ortho-cresol submechanism initially developed from toluene studies is able to predict the measured IDTs fairly without further refinement. It Figure 6. Rate of production analysis for the submechanism of ortho-cresol in the case of a stoichiometric 80/20 iso-octane/ortho-cresol blend at T C = 686 K. P C = 20 bar.Italic numbers correspond to 5% ortho-cresol conversion (i.e., first-stage ignition), bold numbers to 75% ortho-cresol conversion (i.e., final ignition).

Discussion
Insight on the inhibiting effect of ortho-cresol can be gained from the kinetic analysis. The inhibiting effect originates from a radical scavenging effect of ortho-cresol, which is mainly due to the implication of OH radicals into hydrogen abstraction reactions. OH radicals are typically formed in this case by the active low-temperature branching originating from iso-octane which induces first-stage ignition. This scavenging effect results in reduced heat release in the first-stage ignition, thereby limiting the temperature increase and increasing the duration of the second stage of ignition. This results in comparable first-stage ignition delays for blends with up to 30% ortho-cresol, as demonstrated in Figure 2, while the total ignition delay increases monotonically as the ortho-cresol content increases. When R2 is formed in the first stage of ignition, it does not lead to radical chain-branching since its most important conversion channels require higher temperatures. R1 can undergo further reactivity, but then yields phenol, cyclohexadienones, or ortho-hydroxybenzaldehyde, i.e., species which do not contribute to the formation of radicals in significant quantities.
A model developed to describe the co-oxidation of anisole, with an updated phenol and phenoxy submechanism and an ortho-cresol submechanism initially developed from toluene studies is able to predict the measured IDTs fairly without further refinement. It therefore constitutes a useful first step for future research on ortho-cresol. Despite the fair performance of the model, suggestions on further improvement of the model can be made: First, the model is lumped for all three isomers, but the bond dissociation energies are different for the three cresol isomers, as exemplified for the O-H bonds in m-, p-, o-cresol [22]: ortho-cresol = 86.1 < para-cresol = 87.2 < phenol = 88.0, meta-cresol = 89.4 in kcal/mol. This results in global rate constants for the thermal hydrocracking of cresols decreasing in the order: k(ortho) > k(para) > k(meta) [23]. The relative position of methyl and hydroxyl substituents, therefore, has an effect on the reaction rates, especially in the case of H-atom abstraction reactions that lead to the reactivity of ortho-cresol, as suggested by the present work in Figure 6. Second, as described above, most reaction rates specific to ortho-cresol in the model come from analogies with toluene and phenol and can certainly be refined with help from ab initio calculations. Finally, additional pathways could be added into the model, such as additional reaction pathways for R1 and R2 or H-atom abstractions from the aromatic ring. It is often expected that the H-atom abstraction reactions from an aromatic ring are not competitive with those happening on the substituents, especially in low-temperature conditions. Nevertheless, this class of reactions was found to be sensitive in the case of anisole [11], especially because of hydrogen bonding facilitating the formation of a pre-complex. Moreover, this work on the co-oxidation of anisole showed possible reaction pathways for the relevant methoxyphenyl radical with molecular oxygen, and their significance for the accurate modeling of the mole fraction profiles of reaction intermediates. These additional improvements of the kinetic model are likely to improve the description of the kinetics of co-oxidation of ortho-cresol.

Conclusions
In this study, the ignition delay times of iso-octane/ortho-cresol/air mixtures were measured in the ULille RCM, thereby constituting one of the first experimental studies of the combustion kinetics of ortho-cresol. The experimental data are used to check the performances of a model from the literature which has never been validated against the oxidation of ortho-cresol. The model is able to predict the FSIDT and IDT of the blends with an acceptable agreement.
ortho-cresol demonstrated an inhibiting effect on the low-temperature combustion oxidation mechanism of iso-octane, as evidenced by the experimental results. In the lowtemperature range, cresol drastically decreases the reactivity of iso-octane/air mixtures even at a low blending ratio. However, in the intermediate temperature range, the IDTs of the blends are mostly similar to those of pure iso-octane. This can result in an increase in octane number sensitivity when ortho-cresol is blended with iso-octane. The conversion of cresol is led by the H-atom abstraction by hydroxyl radicals, acting as a radical scavenger. The resulting fuel radicals can react with HO 2 reforming cresol, therefore, inhibiting the overall reactivity. In conclusion, ortho-cresol has a strong anti-knock effect as it strongly increases the ignition delay times when blended with iso-octane, however low-temperature reactivity is not fully suppressed. It should also be noted that the fact that the phenoxy radical is easily formed can result in the formation of large quantities of polycyclic aromatic hydrocarbons, as well as their oxygenated counterparts. More fundamental studies, including ab initio calculations of the relevant rate constants, would help gain further insight into the oxidation behavior of cresols.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/en14217105/s1. Table S1: List of the experimental conditions and IDT values, as well as volume profiles of use for the simulations. Funding: This work is a contribution to the CPER research project CLIMIBIO. The authors thank the ADEME, the French Ministère de l'Enseignement Supérieur et de la Recherche, the Hauts-de-France Region and the European Funds for Regional Economical Development for their financial support to this project.

Data Availability Statement:
The data presented in this study are available in the Supplementary Material.