Formulation of a Jet Fuel Surrogate and Its Kinetic Chemical Mechanism by Emulating Physical and Chemical Properties of Real Jet Fuel

: The application of jet fuel in gas turbines and diesel engines adheres to the Army’s single-fuel forward policy, streamlining supply chains. To ensure precise engine combustion numerical studies, surrogate fuels and mechanisms should faithfully replicate real fuel properties and combustion traits. In this work, a new four-component jet fuel surrogate containing 39.05% n-dodecane/21.79% isocetane/11.49% decalin/27.67% toluene by mole fraction is formulated based on a property optimizer. The new-formulated fuel surrogate can satisfactorily emulate the chemical and physical properties of real jet fuel, including cetane number (CN), threshold sooting index (TSI), molecular weight (MW), lower heating value (LHV), the ratio of hydrogen and carbon (H/C), liquid density, viscosity, and surface tension. Furthermore, a reduced and robust kinetic chemical mechanism (containing 124 species and 590 reactions) that could be directly employed in practical engine combustion simulations has also been developed for the proposed surrogate jet fuel. The mechanism is validated through comprehensive experimental data, including ignition delay time (IDT) determined in shock tubes and rapid compression machines (RCMs), species mole fractions measured in premixed flames and jet stirred reactors (JSRs), and laminar flame speeds. Generally, the property deviations of the jet fuel surrogate are less than 2% except for MW (10.73%), viscosity (5.88%), and surface tension (8.71%). The comparison results between the predictions and measurements are in good agreement, indicating that the current kinetic mechanism is capable of reflecting the oxidation process of real jet fuel. The current mechanism can accurately capture variations in the ignition delay time in the negative temperature coefficient (NTC) region as well. In the future, the proposed surrogate jet fuel could be applied in practical engine computational fluid dynamic (CFD) simulations.


Introduction
Jet fuel has been widely adopted as a general fuel for aircraft and ground vehicles in battlefields under the auspices of the U.S. army single-fuel policy [1], which aims to significantly reduce the costs of fuel supply.Furthermore, jet fuel is expected to reduce NOx and particle matter emissions in diesel engines, and thus it would be beneficial to investigate the effects of using jet fuel on the combustion and emission characteristics of diesel engines.However, due to the limitation of current computational capacity, the real jet fuel, which comprises hundreds of compounds, was unable to be directly employed in a practical engine simulation [2,3].Therefore, surrogate jet fuels, which normally consist of several well-characterized compounds, were proposed [4][5][6].The research on surrogate jet fuels is expected to achieve sustainability in the aviation fuel industry.
Previous research studies have been conducted on the development of jet fuel surrogate and their chemical mechanisms.In order to reproduce real fuels as well as possible, surrogates are required to closely emulate the thermo-physical properties and combustion characteristics of target jet fuels.Violi et al. [7] developed a JP-8 surrogate that contained six hydrocarbons and could emulate the distillation curve and thermo-physical properties of the target fuel.A detailed chemical mechanism of the surrogate was also presented in their work and validated by experimental data.Vasu et al. [8] evaluated the aforementioned chemical mechanism by comparing the measured and simulated IDT in a shock tube, and they reported that the mechanism could provide accurate IDT predictions above a temperature of 1000 K while predicting IDT poorly at low temperatures.Dagaut et al. [9] compared the kinetic modelling results of four jet fuel surrogates (one-and three-component mixtures of n-decane, n-propylbenzene, and n-propylcyclohexane) with JSR experimental data.The comparison results showed that the three-component surrogate was the most reliable substitute for the target fuel among all the test fuel surrogates.Gokulakrishnan et al. [10] developed a detailed kinetic model for the four-component kerosene fuel surrogate (n-decane/n-propylcyclohexane/n-propylbenzene/decene), which exhibited good performance in predicting the measured species concentrations.However, the IDT validations were not comprehensively conducted by them.Eddings et al. [11] proposed two six-component jet fuel surrogates named Hex-11 and Hex-12 to emulate a Jet-8 pool fire.The comparison results indicated that the surrogates could accurately capture the burning rate, emissive power, flame height, and puffing frequency of the steady-state pool fire.Dooley et al. [12] formulated a three-component jet fuel surrogate, 'MURI1 (n-decane/isooctane/toluene), based on chemical group theory.The CN and H/C of the target fuel were regarded as the target properties of MURI1.Several devices, including a flow reactor and shock tube, were employed to validate the combustion properties of the proposed surrogate.They further developed an improved jet fuel surrogate named MURI2 comprising n-dodecane, iso-octane, n-propylbenzen, and 1,3,5-trimethylbenzene [13].Apart from CN and H/C, TSI and MW were also selected as the emulating metrics.According to the ignition characteristic study conducted by Malewicki et al. [14], MURI2 exhibited a satisfactory performance on predicting the species profiles of small molecules including oxygen, carbon monoxide, carbon dioxide, and C1-C3 hydrocarbons.Kim et al. [15] proposed two four-component aviation fuel surrogates named UM1 and UM2 based on a model-based optimizer.The physical and chemical properties of the target fuel were reproduced by the two surrogates to some extent.
Despite the extensive research efforts directed towards the creation of various jet fuel surrogates and their accompanying kinetic mechanisms, notable limitations remain apparent in this scientific domain: (1) A recurring challenge is the inability of these surrogates to comprehensively replicate the intricate amalgamation of the physical and chemical properties characteristic of the target fuel.Achieving a harmonious equilibrium between these multifaceted attributes remains an elusive goal.(2) Furthermore, the chemical kinetic mechanisms developed for these surrogates have often exhibited tendencies towards impracticality.They tend to be excessively voluminous or overly intricate, thus impeding their direct application within the confines of practical engine combustion simulations.(3) In addition, the validation of these proposed surrogate mechanisms has not been exhaustive, especially when it comes to crucial parameters such as ignition delay time (IDT), species mole fractions, and flame speeds.The inadequacy of rigorous validation procedures raises pertinent concerns about the reliability and robustness of these surrogates in real-world applications.Addressing these limitations represents an imperative for the advancement of this field.
With the aim of addressing the mentioned issues, a four-component fuel surrogate incorporating the physical and chemical properties of real jet fuel was developed using a property optimizer.Eight properties, including CN, TSI, MW, LHV, H/C, liquid density, viscosity, and surface tension, were selected as the target physical and chemical properties.A reduced kinetic mechanism, comprising 124 species and 590 reactions, for the surrogate was created which can prove advantageous in practical engine combustion simulations.Massive amounts of experimental data, including IDT, species concentration, and laminar flame speed, were collected to validate the mechanism from two aspects: the surrogate mixture and its components.

Target Fuel and Its Properties
The fuel combustion process inside an engine is very complex and always involves intricate physical and chemical processes.As shown in Figure 1, the fuel jet is injected into the cylinder via an injector, breaks down to small liquid droplets under the action of air shear force, and then evaporates and mixes with air to form a combustible fuel-air mixture.The processes are heavily affected by fuel physical properties such as viscosity, liquid density, and surface tension [16,17].After that, the spontaneous ignition and combustion processes of the mixed fuel/air occur, where the chemical properties of the fuel have a significant influence on the events.Herein, CN dominates fuel ignitability, and MW plays a critical role in the diffusive transport process between liquid and the vaporized fuel.H/C influences the local ratio of fuel/air and adiabatic flame temperature, while LHV represents the energy generated from the fuel oxidation.Considering the above processes, eight properties, including CN, TSI, MW, LHV, H/C, liquid density, viscosity, and surface tension, were chosen as target properties for the surrogate in order to comprehensively reproduce the chemical and physical characteristics of real fuel in engine.
rogate was created which can prove advantageous in practical engine comb lations.Massive amounts of experimental data, including IDT, species conce laminar flame speed, were collected to validate the mechanism from two asp rogate mixture and its components.

Target Fuel and Its Properties
The fuel combustion process inside an engine is very complex and alw intricate physical and chemical processes.As shown in Figure 1, the fuel jet i the cylinder via an injector, breaks down to small liquid droplets under the shear force, and then evaporates and mixes with air to form a combustible ture.The processes are heavily affected by fuel physical properties such as uid density, and surface tension [16,17].After that, the spontaneous ignition tion processes of the mixed fuel/air occur, where the chemical properties of a significant influence on the events.Herein, CN dominates fuel ignitabil plays a critical role in the diffusive transport process between liquid and t fuel.H/C influences the local ratio of fuel/air and adiabatic flame temperatur represents the energy generated from the fuel oxidation.Considering the abo eight properties, including CN, TSI, MW, LHV, H/C, liquid density, viscosity tension, were chosen as target properties for the surrogate in order to com reproduce the chemical and physical characteristics of real fuel in engine.In this work, Jet-A POSF-4658 (supplied by Edwards, AFRL-WP, Wrig was selected as the target real fuel for two reasons: First, it is a representative prising different jet-A batches, and its components and various properties h tensively studied.Second, the experimental data on the IDT, laminar flam species concentration profiles of the fuel are available in the literature, so com be more readily conducted.However, since the data for LHV and surface ten In this work, Jet-A POSF-4658 (supplied by Edwards, AFRL-WP, Wright, OH, USA) was selected as the target real fuel for two reasons: First, it is a representative jet fuel comprising different jet-A batches, and its components and various properties have been extensively studied.Second, the experimental data on the IDT, laminar flame speed, and species concentration profiles of the fuel are available in the literature, so comparisons can be more readily conducted.However, since the data for LHV and surface tension of POSF-4658 are not available, the LHV of JP-8 fuel and surface tension of Jet-A were used instead as references in this study, similar to the previous studies [15,18].The relevant properties of POSF-4658 were listed in Table 1.
Table 1.The properties of the target jet fuel [18].

Target Properties
Jet Table 1.The properties of the target jet fuel [18].

Surrogate Fuel Components
Although real jet fuels contain a large number of compounds, according to Ref. [19] conventional jet fuels normally comprise 60% chain paraffins, 20% cycloalkanes, and 20%

Surrogate Fuel Components
Although real jet fuels contain a large number of compounds, according to Ref. [19], conventional jet fuels normally comprise 60% chain paraffins, 20% cycloalkanes, and 20% aromatics.Hence, it is reasonable to adopt the well-characterized compounds recognized in real jet fuels as the main surrogate components.Moreover, the components of jet fuel mainly belonged to C9-C16 hydrocarbons, which indicates that HC class and molecule size are also important criteria for surrogate components.The components whose chemical kinetic mechanisms have been well established are more preferred in our study.Based on these considerations, four components, including n-dodecane, isocetane, decalin, and toluene, were finally selected.The properties (CN, TSI, MW, LHV, and H/C) and formulae of the four components are listed in Table 2.The temperature-dependent properties, density, viscosity, and surface tension, are plotted in Figure 2. It can be easily seen that a single component is unable to represent real jet fuel because of the significant property differences between the single component and real fuel.

Formation of Jet Fuel Surrogate
In this work, we improved the optimization algorithm from our previous study [18] to obtain the optimum component proportions of a jet fuel surrogate.The calculations for the target properties are showed in Table 3.
As reported in [20], the CN of the mixture can be calculated using the volume fraction average of the component CNs.The TSI of a mixture is obtained by calculating the mole fraction average of the pure component's TSI [23], while the LHV is estimated using the mass fraction average of the LHV.H/C and MW were directly determined by the component formula.The mixture density is also computed by the volume fraction average of component density.The viscosity is obtained using the Grunberg-Nissan equation [24,25], whereas surface tension is derived from the parachor correlation [26].Finally, the equations proposed by Kim et al. [15] were employed as the merit functions: where i and j represent the numbers of properties and experimental data points, respectively.Num data,i refers to the total number of measurements of the ith property, and Num target represents the total number of target properties.Q i,j,cal and Q i,j,exp are the computed and measured properties, respectively.
Table 3.The estimation methods [18] of the properties of fuel surrogate.

MW Average of mole fraction:
C i H i is the number of hydrogen atoms of component i;

Viscosity
Grunberg-Nissan equation [24]: G ij is the binary interaction parameter

Surface tension
Parachor correlation: P L,mix = 0.5 ∑ ∑ X i X j P i + P j σ(T) mix is liquid surface tension, P is parachor, ρ(T) L,mix,molar is liquid mixture molar density A flowchart of the whole optimization process is given in Figure 3.It should be noted that the component will be removed when its mole fraction is less than 0.3% during the optimization process, the same as in our previous study [18].Based on the above calculations, the optimum jet fuel surrogate was achieved, and it comprised 39.05% ndodecane, 21.79% isocetane, 11.49% decalin, and 27.67% toluene by mole fraction and was abbreviated as JFS.
Table 4 shows the comparison results among the JFS surrogate and other jet fuel surrogates, including UM1, UM2, MURI2, S5 [27], and HEX12 [11].In general, the property deviations of the surrogate JFS are less than 2% except for MW (10.73%), viscosity (5.88%), and surface tension (8.71%).The high deviations of MW and surface tension can be regarded as acceptable compared to other surrogates.As shown in Table 4, the deviation of the MW of S5 is 12.11%, and the averaged deviations of surface tension of UM2, Hex12, and S5 are larger than 15%, which are all larger than the corresponding property deviations of the surrogate JFS.It is also observed that the CN, MW, H/C and, LHV of the surrogate UM1 show good agreement with the target fuel; however, the deviations of TSI (−21.52%) and viscosity (−21.2%) are too large.The surface tension of UM2 deviates from the target property by 15.8%, even though CN, LHV, and density are well captured by the surrogate UM2.The LHV, MW, CN, and H/C are well emulated by MURI2; however, the deviation of viscosity (34.608%) is too large.As for the properties of HEX12 and S5, the majority of them greatly deviate from the target properties.Overall, the surrogate JFS performs better than other surrogates at reproducing the chemical and physical properties of real jet fuel; hence, it is more likely to be adopted for practical engine combustion simulation.
that the component will be removed when its mole fraction is less than optimization process, the same as in our previous study [18].Based on t tions, the optimum jet fuel surrogate was achieved, and it comprised 39. 21.79% isocetane, 11.49% decalin, and 27.67% toluene by mole fraction ated as JFS.Table 4 shows the comparison results among the JFS surrogate and rogates, including UM1, UM2, MURI2, S5 [27], and HEX12 [11].In gene deviations of the surrogate JFS are less than 2% except for MW (10.73%), and surface tension (8.71%).The high deviations of MW and surface t garded as acceptable compared to other surrogates.As shown in Table 4 the MW of S5 is 12.11%, and the averaged deviations of surface tension and S5 are larger than 15%, which are all larger than the corresponding tions of the surrogate JFS.It is also observed that the CN, MW, H/C and rogate UM1 show good agreement with the target fuel; however, the (−21.52%)and viscosity (−21.2%) are too large.The surface tension of UM

Methodology
In this study, to facilitate the simulation of a practical engine combustion, a so-called decoupling methodology [5,18] was adopted for developing the JFS chemical mechanisms.The chemical kinetic mechanism was established through a systematic formation process, starting from C16 reactions and progressing towards H 2 /O 2 /C1 reactions.The subsequent steps involved incorporating a NOx sub-mechanism and polycyclic aromatic hydrocarbon (PAH).The reduction and optimization procedures employed in this process are briefly outlined as follows: (1) The initial stage of the reduction and optimization process involves conducting a reaction pathway analysis to identify the key reactions.Subsequently, unimportant species and reactions are eliminated from the initial kinetic model.Simultaneously, the rate of production (ROP) and sensitivity analyses are performed to evaluate the remaining species and reactions.This allows for a clear understanding of the impact of each reaction on the oxidation process.( 2) Subsequently, the reaction rate constants were optimized to improve the agreement between the simulated and experimental data of fuel ignition delay time (IDT).
(3) Afterward, the concentrations of species and laminar flame speeds predicted by the reduced mechanism were compared to the corresponding measurements.This allowed for references to fine-tune the reaction rate constants further.(4) Steps 1 to 3 were iteratively repeated until the desired size and accuracy of the mechanism were attained.

Toluene Sub-Mechanism
The current chemical kinetic model mainly comprises four sub-mechanisms, including the toluene sub-mechanism, decalin sub-mechanism, n-dodecane sub-mechanism, and isocetane sub-mechanism (Figure 4).

Toluene Sub-Mechanism
The current chemical kinetic model mainly comprises four sub-mechanisms, including the toluene sub-mechanism, decalin sub-mechanism, n-dodecane sub-mechanism, and isocetane sub-mechanism (Figure 4).The toluene sub-mechanism was tuned from the reduced TRF-PAH mechanism [28], and it was employed as the initial mechanism based on two reasons: (1) The size of the TRF-PAH mechanism is compact, which is advantageous for three-dimensional engine combustion simulation; and (2) The mechanism has good performance in predicting spe- The toluene sub-mechanism was tuned from the reduced TRF-PAH mechanism [28], and it was employed as the initial mechanism based on two reasons: (1) The size of the TRF-PAH mechanism is compact, which is advantageous for three-dimensional engine combustion simulation; and (2) The mechanism has good performance in predicting species concentrations and flame speeds.
Figure 4 illustrates the principal reactions of the reduced toluene mechanism.Toluene is predominantly consumed through three pathways: the H-atom abstraction reactions leading to the production of C 6 H 4 CH 3 (R6), C 6 H 5 CH 2 (R1), and the reaction R5.Among these, the formation of C 6 H 5 CH 2 is the primary pathway for toluene consumption, while C6H5OH ranks as the second most abundant product.Through some intermediate species, such as C 6 H 5 CH 2 O, C 6 H 5 , and C 6 H 5 O, the toluene molecule ultimately converts into the small molecules.Some species, such as C 6 H 5 CH 2 OO, C 6 H 5 CO, and OC 6 H 4 CH 3 , which are absent in Figure 4, were also incorporated into the toluene sub-mechanism.
The NOx and PAH sub-mechanisms of the original TRF-PAH mechanism were retained in the current mechanism.The NOx sub-mechanism includes four species and 12 reactions, and it contains thermal reactions and N 2 O-intermediate reactions.The PAH mechanism is able to reproduce PAH formation up to a four-ring PAH.Herein, H-atom abstraction reactions, methyl substitution reactions, and hydrogen abstraction acetylene addition (HACA) make great contributions to PAH formation, and hence the reactions of MAHs (e.g., A 1 C 2 H and A 1 C 2 H 3 ) were considered.As reported in Ref. [28], the predictions obtained from the PAH mechanism demonstrated outstanding agreement with the measured data collected from eight different flames, which indicated that the present PAH mechanism is satisfactory to some extent.

Decalin Sub-Mechanism
The decalin mechanism was reduced from the detailed decalin oxidation mechanism proposed by Dagaut et al. [29].As shown in Figure 4, the current decalin mechanism contains both low-and high-temperature mechanisms, and so the mechanism exhibits the capability to accurately describe decalin oxidation across a broad temperature range.The reactions framed by the dashed box (Figure 4) belong to the low-temperature mechanism.
Decalin radicals can be diverse considering the symmetrical structure of a decalin molecule.Three C 10 H 17 radicals lumped as 'RDECALIN' were produced via decalin consumption reactions with small molecules including O 2 , H, and HO 2 , as shown in Figure 4 (R12-R15).In accordance with Yu et al. [30], the low-temperature and hightemperature mechanisms exhibit significant differences in their core reactions.In the low-temperature regime, the primary consumption pathway for RDECALIN involves the O 2 addition reaction (R16), while other consumption reactions are negligible.The formation of alkylperoxy radicals (RDECOO) via R16 is followed by their isomerization to hydroperoxy alkyl (QDECOOH).QDECOOH, in turn, undergoes a conversion to ZDECA (lumped C 10 H 17 O 4 isomers) through the reaction R18.ZDECA is then consumed, leading to the formation of C 10 H 16 O 3 isomers (KHDECA), and the decalin chain-branching sequence is completed by the reaction R20, which represents the KHDECA decomposition at low temperatures.Conversely, in the high-temperature regime, the decomposition of RDECALIN primarily occurs through reactions that tend to generate C 5 H 8 species and cyclopentylcyclopentene (DCYC5).Then C 5 H 8 decomposes to cyclopentadiene (CYC 5 H 6 ) via the H-abstraction reaction (R23).Eventually, the big hydrocarbon molecules are converted into small C2-C4 molecules.

N-Dodecane Sub-Mechanism
The n-dodecane mechanism was reduced via our previous work [31] and was briefly introduced here.As shown in Figure 4, n-dodecane is consumed by three reactions (R29-R31) to produce dodecyl radicals (C 12 H 25 ) at low temperatures.C 12 H 25 then transforms to peroxy dodecyl radicals (C 12 H 25 OO) via the pathway reaction R32, which is significant for the formation of the negative temperature coefficient (NTC) region.Hydroperoxy dodecyl radicals (C 12 H 24 OOH) are produced by a isomerization reaction (R33).C 12 H 24 OOH transforms into O 2 C 12 H 24 OOH by reacting with oxygen (R34).C 12 KET radicals are produced from the reaction R35, and they further decompose into a few species, including C 5 H 11 CO, with the reaction R36.The small species will be generated through the decomposition reaction R37.However, at high temperatures, dodecyl radicals first convert to C 12 H 24 by a reaction with oxygen (R38) and then decompose into C2-C4 molecules via the reaction R39 or directly produce the small hydrocarbon molecules by R40.

Isocetane Sub-Mechanism
The isocetane sub-mechanism was also taken from our previous work [32].A short depiction of the isocetane sub-mechanism is given here.The consumption of isocetane is completed by the reactions R41-R43 at low temperatures (Figure 4).Ketohydroperoxide is produced via the reactions R44-R47, while the decomposition of ketohydroperoxide is completed by the reactions R48 and R49.At high temperatures, all reactions are represented by the three reactions R50, R51, and R52.

Results and Discussion
The experimental data, including IDTs, species concentrations, and laminar flame speeds determined by various devices, were used for validating the current kinetic mechanism.The mechanism was first validated for each component and then for the surrogate fuel mixture.The validation of the individual components will give a more comprehensive assessment on the accuracy and reliability of the JFS mechanism.The simulations were performed on CHEMKIN-PRO Version-19.0 software [33].

IDT
It is necessary to validate the current mechanism on IDT because the fuel autoignition behaviors characterized by IDT are important for engine combustion and emission performances [34].Figure 5 exhibits a comparison of the measured and simulated IDTs for toluene in a shock tube under varying conditions of temperature, pressure, and equivalence ratios.The measurements were determined by Shen et al. [35], and the predicted IDT with the current mechanism and the based mechanism [28] were exhibited.Both the reduced mechanism and the based mechanism demonstrated strong agreement with the measurements.The current mechanism showed better performance under an equivalence ratio of 1.0.Compared to the based mechanism, the IDTs predicted by the current mechanism are closer to the measurements at low pressure.In addition, the predicted IDTs of toluene are compared with the determined IDTs in an RCM [36].The corresponding results are presented in Figure 6.The toluene IDTs were slightly underestimated, but the discrepancy is acceptable.

Species Concentration
As emphasized by Pitz et al. [37], it is equally crucial to validate the profiles of major species concentrations alongside the validation of fuel IDT.Hence, the current mechanism was validated for the major species profiles during the oxidation process.In Figure 7, the concentration profiles of critical species in the premixed toluene flame are displayed.These species include reactants (toluene and oxygen), an inactive gas (argon), and products (such as carbon monoxide and hydrogen).Li et al. [38] experimentally investigated premixed toluene flames at three equivalence ratios (Φ = 0.75, 1.0, and 1.5) under a pressure of 10 atm.A satisfactory agreement was observed between the predicted and observed species profiles for all flames, with a notable consistency observed in the case of the stoichiometric toluene flame.The argon species concentration initially decreased and then stabilized further downstream of the burner.The model accurately predicted this trend, demonstrating its capability to replicate the mole expansion effects observed in the toluene flame, as reported in reference [38].Additionally, the model successfully captured the overall profiles of the final products, specifically water and carbon dioxide, and the maximum deviations between the measurements and predictions for final products are less than 10%.Carbon monoxide displayed a profile resembling that of hydrogen, gradually surpassing carbon dioxide as the predominant carbon product.Despite the discrepancies in the predicted species concentrations of toluene and oxygen compared to the corresponding measurements at equivalence ratios of 0.75 and 1.5, the model successfully represented the depletion positions of these species with accuracy.

IDT
It is necessary to validate the current mechanism on IDT because the fuel auto-ignition behaviors characterized by IDT are important for engine combustion and emission performances [34].Figure 5 exhibits a comparison of the measured and simulated IDTs for toluene in a shock tube under varying conditions of temperature, pressure, and equivalence ratios.The measurements were determined by Shen et al. [35], and the predicted IDT with the current mechanism and the based mechanism [28] were exhibited.Both the reduced mechanism and the based mechanism demonstrated strong agreement with the measurements.The current mechanism showed better performance under an equivalence ratio of 1.0.Compared to the based mechanism, the IDTs predicted by the current mechanism are closer to the measurements at low pressure.In addition, the predicted IDTs of toluene are compared with the determined IDTs in an RCM [36].The corresponding results are presented in Figure 6.The toluene IDTs were slightly underestimated, but the discrepancy is acceptable.[35]) and predicted (solid lines: current mechanism, dash lines: the based mechanism [28]) toluene IDTs determined in shock tube.

Species Concentration
As emphasized by Pitz et al. [37], it is equally crucial to validate the profiles of major species concentrations alongside the validation of fuel IDT.Hence, the current mechanism was validated for the major species profiles during the oxidation process.In Figure 7, the   Figure 8 shows the toluene oxidation process in a JSR [39] at a residence time of 0.6 s and a pressure of 10 atm.The concentration evolution of both the major species (CO, CO 2 , and C 7 H 8 ) and the minor species (C 2 H 2 and C 4 H 4 ) are presented in Figure 8.The current mechanism rendered good performance in predicting the concentration profiles of C 7 H 8 , CO, and CO 2 .C 2 H 2 is esteemed as a significant precursor of PAH and soot formation [40].Moreover, the proposed model exhibits its capability to predict the toluene concentration evolution.Furthermore, the mechanism successfully predicted the species profiles of C 4 H 4 under different conditions, except for the scenario with an equivalence ratio of 1.5.

Laminar Flame Speed
The laminar flame speed (LFS) is a valuable parameter that provides insights into the reactivity, heat release rate, and propagation characteristics of actual fuels.Figure 9 shows the comparison between the simulated and measured LFSs of counter-flow toluene flames.The measurements were conducted under two conditions: p = 1 atm, Tu = 298 K [41], and p = 3 atm, Tu = 450 K [42], where Tu denotes the temperature of the unburned mixture.It can be observed that the mechanism performed well in predicting the flame behavior under atmospheric conditions.However, the discrepancy of the predictions and measurements were larger under fuel-rich conditions compared to oxygen-rich conditions.Therefore, we can conclude that the current mechanism is more suitable for modeling lean-burn engine combustion.

Verifications of Decalin 4.2.1. IDT
Figure 10 displays the comparisons of the simulated and measured ignition delay times for decalin in a shock tube.The experimental data were obtained from previous studies [43,44], while the simulation results were obtained with both the present mechanism and the detailed mechanism [29].Except for the case of Φ = 0.5 and p = 12 atm, the current mechanism consistently outperformed the detailed mechanism across all other conditions, as demonstrated in Figure 10.However, considering the results achieved by the detailed mechanism, the level of discrepancy observed under the mentioned condition is still within an acceptable range.These results confirm that the current mechanism effectively captures the dependence of IDT on pressure and equivalence ratio, particularly at two specific conditions: (1) p = 20 atm, Φ = 0.5; (2) p = 40 atm, Φ = 1.0.
flames.The measurements were conducted under two conditions: p = 1 atm, Tu = 298 K [41], and p = 3 atm, Tu = 450 K [42], where Tu denotes the temperature of the unburned mixture.It can be observed that the mechanism performed well in predicting the flame behavior under atmospheric conditions.However, the discrepancy of the predictions and measurements were larger under fuel-rich conditions compared to oxygen-rich conditions.Therefore, we can conclude that the current mechanism is more suitable for modeling lean-burn engine combustion.

IDT
Figure 10 displays the comparisons of the simulated and measured ignition delay times for decalin in a shock tube.The experimental data were obtained from previous studies [43,44], while the simulation results were obtained with both the present mechanism and the detailed mechanism [29].Except for the case of Φ = 0.5 and p = 12 atm, the current mechanism consistently outperformed the detailed mechanism across all other conditions, as demonstrated in Figure 10.However, considering the results achieved by the detailed mechanism, the level of discrepancy observed under the mentioned condition is still within an acceptable range.These results confirm that the current mechanism effectively captures the dependence of IDT on pressure and equivalence ratio, particularly at two specific conditions: (1) p = 20 atm, Φ = 0.5; (2) p = 40 atm, Φ = 1.0.
Figure 11 depicts the ignition delay time (IDT) of decalin over a temperature range of 770 K to 1250 K. Non-thermal combustion (NTC) behaviors were observed within this temperature range.At low-pressure conditions (p = 20 atm), the NTC region started at around 910 K, while at high temperatures, the NTC region initiated at a higher temperature of approximately 940 K.The shift in the NTC region suggests that pressure significantly influences the NTC behavior, which aligns with similar findings reported in Ref. [45].  Figure 11 depicts the ignition delay time (IDT) of decalin over a temperature range of 770 K to 1250 K. Non-thermal combustion (NTC) behaviors were observed within this temperature range.At low-pressure conditions (p = 20 atm), the NTC region started at around 910 K, while at high temperatures, the NTC region initiated at a higher temperature of approximately 940 K.The shift in the NTC region suggests that pressure significantly influences the NTC behavior, which aligns with similar findings reported in Ref. [45].
In Figure 12, the decalin IDT is presented as determined in a rapid compression machine over a broad temperature range.The measurements were conducted at various equivalence ratios and a 15 atm compression pressure [30].The results clearly demonstrate that IDTs decrease with increasing equivalence ratios.The variations in IDTs among different conditions were relatively small within the temperature range of 631 to 750 K, indicating a weak influence of equivalence ratio on IDTs at low temperatures.As the temperature increased, the differences in the ignition delay times (IDTs) became more pronounced.This can be attributed to the intensified chain-branch reactions occurring between the decalin molecules and hydroperoxyl radicals in both the negative-temperature-coefficient and high-temperature zones.The current mechanism successfully captures these observed behaviors.In Figure 12, the decalin IDT is presented as determined in a rapid compres chine over a broad temperature range.The measurements were conducted at equivalence ratios and a 15 atm compression pressure [30].The results clearly strate that IDTs decrease with increasing equivalence ratios.The variations in IDT different conditions were relatively small within the temperature range of 631 t indicating a weak influence of equivalence ratio on IDTs at low temperatures.As perature increased, the differences in the ignition delay times (IDTs) became m nounced.This can be attributed to the intensified chain-branch reactions occur tween the decalin molecules and hydroperoxyl radicals in both the negative-temp coefficient and high-temperature zones.The current mechanism successfully these observed behaviors.

Species Concentration
Figure 13 exhibits the species profiles of the premixed decalin flam dictions and measurements showed excellent agreement across variou tios.According to the ROP analysis, CO is predominantly generated t involving HCO and HCCO, while its consumption occurs through its droxyl radicals, leading to the formation of CO2.The formation of H2O enced by H-atom abstraction reactions as well as the reaction between produced through a reaction between decalin and formaldehyde (CH2 of depletion for decalin species varied with different equivalence ratios, the stoichiometric decalin flame, the depletion position was observed t

Species Concentration
Figure 13 exhibits the species profiles of the premixed decalin flames [46].The predictions and measurements showed excellent agreement across various equivalence ratios.According to the ROP analysis, CO is predominantly generated through reactions involving HCO and HCCO, while its consumption occurs through its reaction with hydroxyl radicals, leading to the formation of CO 2 .The formation of H 2 O is notably influenced by H-atom abstraction reactions as well as the reaction between OH and H 2 .H 2 is produced through a reaction between decalin and formaldehyde (CH 2 O).The positions of depletion for decalin species varied with different equivalence ratios, and in the case of the stoichiometric decalin flame, the depletion position was observed to be closer to the burner surface.Similar observations were reported in the referenced study [38].
produced through a reaction between decalin and formaldehyde (CH2O).The positio of depletion for decalin species varied with different equivalence ratios, and in the case the stoichiometric decalin flame, the depletion position was observed to be closer to t burner surface.Similar observations were reported in the referenced study [38].The species concentration profiles of decalin were also examined in the JSR [29].C 2 H 2 and CH 4 were selected for validating the accuracy of the current mechanism since they are the significant intermediate hydrocarbons produced during decalin oxidation.The current mechanism accurately predicted the concentration variations of C 2 H 2 and CH 4 , although it slightly overestimated the decalin concentration.It was observed that the concentration of CO initially increased and then reached a steady level above 900 K.The concentration of CO 2 exhibited a continuous increase with rising temperature.Overall, the current mechanism demonstrated good performance in reproducing these trends.Furthermore, the concentration profiles of CH 2 O, a major aldehyde produced during the oxidation process of decalin, were also accurately captured by the current mechanism (Figure 14).current mechanism demonstrated good performance in reproducing these trends.Furthermore, the concentration profiles of CH2O, a major aldehyde produced during the oxidation process of decalin, were also accurately captured by the current mechanism (Figure 14).[29]) and predicted (lines) concentration profiles of species resulting from the oxidation of decalin in the JSR at an equivalence ratio of 1.0 and a compression pressure of 10 atm.

Laminar Flame Speed
Figure 15 illustrates the determined and predicted LFSs for the counter-flow decalin flame at atmospheric pressure [47].The initial mixture temperature for the flame is 443 K.It is evident that the current mechanism accurately predicts the laminar flame speed of decalin, although the deviation between measured and predicted LFSs becomes larger under high equivalence ratios.[29]) and predicted (lines) concentration profiles of species resulting from the oxidation of decalin in the JSR at an equivalence ratio of 1.0 and a compression pressure of 10 atm.

Laminar Flame Speed
Figure 15 illustrates the determined and predicted LFSs for the counter-flow decalin flame at atmospheric pressure [47].The initial mixture temperature for the flame is 443 K.It is evident that the current mechanism accurately predicts the laminar flame speed of decalin, although the deviation between measured and predicted LFSs becomes larger under high equivalence ratios.
thermore, the concentration profiles of CH2O, a major aldehyde produced idation process of decalin, were also accurately captured by the current me ure 14).[29]) and predicted (lines) concentration profil sulting from the oxidation of decalin in the JSR at an equivalence ratio of 1.0 and pressure of 10 atm.

Laminar Flame Speed
Figure 15 illustrates the determined and predicted LFSs for the counte flame at atmospheric pressure [47].The initial mixture temperature for the It is evident that the current mechanism accurately predicts the laminar fl decalin, although the deviation between measured and predicted LFSs b under high equivalence ratios.Figure 16 shows the predicted and measured n-dodecane IDTs at two equivalence ratios (0.5 and 1.0) under three pressures (14 atm, 20 atm, and 40 atm).The shock tube measurements were taken from the Refs.[48][49][50].The results confirmed that the current mechanism was capable of capturing the variation trends of the measured IDTs, except in the case of p = 40 atm.The NTC behaviors observed in the temperature range from 750 K to 900 K were well reflected by the mechanism.In addition, the experimental phenomenon that IDT decreases with an increase in pressure was also captured by the mechanism.
Figure 16 shows the predicted and measured n-dodecane IDTs at two eq ratios (0.5 and 1.0) under three pressures (14 atm, 20 atm, and 40 atm).The s measurements were taken from the Refs.[48][49][50].The results confirmed that th mechanism was capable of capturing the variation trends of the measured IDTs the case of p = 40 atm.The NTC behaviors observed in the temperature range fr to 900 K were well reflected by the mechanism.In addition, the experimental phe that IDT decreases with an increase in pressure was also captured by the mecha

Species Concentration
Malewicki et al. [14] conducted JSR experiments to investigate the n-dodecane process.The experiments were carried out at a residence time of 1 s, a pressure of 10 three equivalence ratios.The comparison between the computed and measured shown in Figure 17.The prediction for the variation trend of O2 and n-dodecane (C temperature was consistent with the measurements under all conditions.H2O w underestimated, but the deviation is still at an acceptable level.There were some between the simulated and determined mole fractions of C2H2 and CO2; nevertheless all simulated profiles of C2H2 and CO2 were satisfactory.

Species Concentration
Malewicki et al. [14] conducted JSR experiments to investigate the n-dodecane oxidation process.The experiments were carried out at a residence time of 1 s, a pressure of 10 atm, and three equivalence ratios.The comparison between the computed and measured results is shown in Figure 17.The prediction for the variation trend of O 2 and n-dodecane (C 12 H 26 ) with temperature was consistent with the measurements under all conditions.H 2 O was slightly underestimated, but the deviation is still at an acceptable level.There were some deviations between the simulated and determined mole fractions of C 2 H 2 and CO 2 ; nevertheless, the overall simulated profiles of C 2 H 2 and CO 2 were satisfactory.

Laminar Flame Speed
The measured data of the n-dodecane laminar flame speed were adopted here from the experiment conducted by Kumar et al. [51].The measurement was performed at atmospheric pressure and unburned mixture temperatures of 400 K and 470 K.It can be easily seen from Figure 18 that the predictions showed fairly good agreement with the experimental data under all conditions.

Verifications of Isocetane 4.4.1. IDT
The IDTs of isocetane were measured in a shock tube under two pressures (10 atm and 40 atm) and three equivalence ratios (0.5, 1.0, and 1.5) by Oehlschaeger et al. [52].Figure 19 presents the predictions and measurements of isocetane IDT, and the maximum deviations between the measured and predicted IDT appeared under the condition of p = 10 atm and Φ = 1.0, confirming the derived mechanism can accurately predict IDT.The decreasing trend of IDT with increasing pressure was perfectly captured by the mechanism.Due to the absence of experimental data, the NTC behaviors were not obviously observed at low pressure and low equivalence ratios.Therefore, it is necessary to do further IDT validations for isocetane in the future.
temperature was consistent with the measurements under all conditions.H2O was s underestimated, but the deviation is still at an acceptable level.There were some dev between the simulated and determined mole fractions of C2H2 and CO2; nevertheless, th all simulated profiles of C2H2 and CO2 were satisfactory.[14]) and predicted (lines) mole fractions of the n-dodecane oxidation in JSR experiments at different equivalence ratios.

Laminar Flame Speed
The measured data of the n-dodecane laminar flame speed were adop the experiment conducted by Kumar et al. [51].The measurement was pe mospheric pressure and unburned mixture temperatures of 400 K and 47 easily seen from Figure 18 that the predictions showed fairly good agree experimental data under all conditions.

Species Concentration
Figure 20 shows the concentrations profiles of the major species produced during isocetane oxidation in a JSR [53].The measurements were performed at a residence time of 1 s, under a pressure of 10 atm, and under three equivalence ratios (0.5, 1.0, and 2.0).It can be seen that the mechanism nicely predicts the consumption trend of isocetane under all conditions.The species mole fractions of H2O and CO were also well predicted.The concentration of isocetane was overestimated in the high temperature range when Φ = 0.5.For intermediate hydrocarbon products, the species profiles of CH4 were not excellently reproduced by the mechanism, especially in the case of Φ = 1.0.

Species Concentration
Figure 20 shows the concentrations profiles of the major species produced during isocetane oxidation in a JSR [53].The measurements were performed at a residence time of 1 s, under a pressure of 10 atm, and under three equivalence ratios (0.5, 1.0, and 2.0).It can be seen that the mechanism nicely predicts the consumption trend of isocetane under all conditions.The species mole fractions of H 2 O and CO were also well predicted.The concentration of isocetane was overestimated in the high temperature range when Φ = 0.5.For intermediate hydrocarbon products, the species profiles of CH 4 were not excellently reproduced by the mechanism, especially in the case of Φ = 1.0.

Verifications of Jet Fuel
In this section, comparisons of ignition delay times, species mole fractions, and laminar flame speeds were further carried out for the newly developed JFS surrogate and the target jet fuel POSF-4658.

IDT
The shock tube experiments of jet fuel (POSF-4658) were conducted by S. Vasu et al. [8] and Wang et al. [54], as shown in Figure 21.[53])and predicted (lines) major species profiles of isocetan JSR experiments at different equivalence ratios.

Verifications of Jet Fuel
In this section, comparisons of ignition delay times, species mole fractions, and la inar flame speeds were further carried out for the newly developed JFS surrogate and target jet fuel POSF-4658.The measured (symbols [53])and predicted (lines) major species profiles of isocetane in JSR experiments at different equivalence ratios.

Verifications of Jet Fuel
In this section, comparisons of ignition delay times, species mole fractions, and lam inar flame speeds were further carried out for the newly developed JFS surrogate and the target jet fuel POSF-4658.

IDT
The shock tube experiments of jet fuel (POSF-4658) were conducted by S. Vasu et al [8] and Wang et al. [54], as shown in Figure 21.The comparisons between the determined and simulated IDTs of jet fuel at three pres sures are exhibited in Figure 21a, while Figure 21b shows the results at three equivalence The comparisons between the determined and simulated IDTs of jet fuel at three pressures are exhibited in Figure 21a, while Figure 21b shows the results at three equivalence ratios: Φ = 0.5, 1.0, and 1.5.The overall trends of the jet fuel IDTs were well reflected by the current mechanism, especially at low pressure and low equivalence ratios, which can be attributed to the reason that the H/C, CN, and LHV of the proposed JFS surrogate are similar to those of the target fuel.Although the predicted NTC region underwent a slight shift to the left under the conditions of Φ = 1.5 and 20 atm pressure, the NTC behavior was satisfactorily captured via the developed mechanism.Furthermore, the maximum deviation between the predicted and measured IDTs is less than 0.001 s, and the small deviations between the predicted and measured IDTs can be attributed to two reasons: (1) Some reactions and species are removed from the reduced mechanism, which inevitably makes the mechanism unable to reproduce the measurements as well as the detailed mechanism.(2) There are still some differences between the properties of the JFS surrogate and the target fuel, which in turn influence the performance of the associated mechanism.
Figure 22 shows the comparison results between the measured IDTs and the IDTs predicted from the current JFS mechanism as well as the UM1 and UM2 mechanisms at Φ = 1.0 and p = 20 atm.Both the surrogates UM1 and UM2 contain four components, and the components of UM2 and the current JFS surrogate are identical.The detailed mechanisms of UM1 and UM2 were developed by Kim et al. [15].In general, all of the mechanisms give good predictions of the measured IDTs of jet fuels.The predictions still suffer some deviations in IDT against the experimental data within the intermediate temperature range (from 855 to 1030 K), even with the detailed mechanisms, and the maximum deviation between the predicted and measured IDTs is about 0.0015 s.It is worth mentioning that the reduced JFS mechanism has a comparable emulation capability to the UM1 and UM2 mechanisms, although the UM1 and UM2 mechanisms are much more detailed compared to the current mechanism.
Sustainability 2023, 15, 13792 ratios: Φ = 0.5, 1.0, and 1.5.The overall trends of the jet fuel IDTs were well refle the current mechanism, especially at low pressure and low equivalence ratios, wh be attributed to the reason that the H/C, CN, and LHV of the proposed JFS surro similar to those of the target fuel.Although the predicted NTC region underwent shift to the left under the conditions of Φ = 1.5 and 20 atm pressure, the NTC behav satisfactorily captured via the developed mechanism.Furthermore, the maximum tion between the predicted and measured IDTs is less than 0.001 s, and the sma tions between the predicted and measured IDTs can be attributed to two reasons: ( reactions and species are removed from the reduced mechanism, which inevitabl the mechanism unable to reproduce the measurements as well as the detailed mec (2) There are still some differences between the properties of the JFS surrogate target fuel, which in turn influence the performance of the associated mechanism Figure 22 shows the comparison results between the measured IDTs and t predicted from the current JFS mechanism as well as the UM1 and UM2 mechanis = 1.0 and p = 20 atm.Both the surrogates UM1 and UM2 contain four components, components of UM2 and the current JFS surrogate are identical.The detailed mech of UM1 and UM2 were developed by Kim et al. [15].In general, all of the mechanis good predictions of the measured IDTs of jet fuels.The predictions still suffer som ations in IDT against the experimental data within the intermediate temperatur (from 855 to 1030 K), even with the detailed mechanisms, and the maximum d between the predicted and measured IDTs is about 0.0015 s.It is worth mention the reduced JFS mechanism has a comparable emulation capability to the UM1 an mechanisms, although the UM1 and UM2 mechanisms are much more detailed co to the current mechanism.

Species Concentration
Figure 23 shows the species profiles of the JFS surrogate in comparison with perimental data determined in a shock tube at the temperature range of 890 K to by Malewicki et al. [14].

Species Concentration
Figure 23 shows the species profiles of the JFS surrogate in comparison with the experimental data determined in a shock tube at the temperature range of 890 K to 1680 K by Malewicki et al. [14].The concentration profiles of CH4 and C2H2 were well captured by the current mechanism, except for the slight overprediction for CH4 at the lean condition (Φ = 0.46), which implies that the mechanism is reliable in predicting the soot and PAH formulation in engine combustion.The decaying trend of O2 and the increasing trend of CO2 were both greatly reflected by the mechanism.As shown in Figure 23, the current mechanism also showed good performance in predicting the CO concentration except for a slight underestimation under the rich-and high-temperature condition (Figure 22b).Overall, the current mechanism exhibited a satisfactory prediction of the major species concentrations of the JFS surrogate, although the maximum deviation between the measurements and the predictions of CH4 was larger than the measured value.In fact, similar deviations were also observed even using a detailed mechanism, as reported by Malewicki et al. [14].

Laminar Flame Speed
The comparison results between the simulated and measured laminar flame speeds is shown in Figure 24.The laminar flame speeds of jet fuel were measured by Dooley et al. [13] at atmospheric pressure and the unburned mixture temperatures of 400 K and 470 K.The predicted values showed a surprising consistency with the measurements, and the maximum discrepancy between them was less than 9.4%.According to Li et al. [55], H2, CO, and small hydrocarbons have decisive impacts on flame propagation; hence, great consistency can be attained through the adoption of a detailed C1/H2/CO sub-mechanism in the current JFS mechanism.The concentration profiles of CH 4 and C 2 H 2 were well captured by the current mechanism, except for the slight overprediction for CH 4 at the lean condition (Φ = 0.46), which implies that the mechanism is reliable in predicting the soot and PAH formulation in engine combustion.The decaying trend of O 2 and the increasing trend of CO 2 were both greatly reflected by the mechanism.As shown in Figure 23, the current mechanism also showed good performance in predicting the CO concentration except for a slight underestimation under the rich-and high-temperature condition (Figure 22b).Overall, the current mechanism exhibited a satisfactory prediction of the major species concentrations of the JFS surrogate, although the maximum deviation between the measurements and the predictions of CH 4 was larger than the measured value.In fact, similar deviations were also observed even using a detailed mechanism, as reported by Malewicki et al. [14].

Laminar Flame Speed
The comparison results between the simulated and measured laminar flame speeds is shown in Figure 24.The laminar flame speeds of jet fuel were measured by Dooley et al. [13] at atmospheric pressure and the unburned mixture temperatures of 400 K and 470 K.The predicted values showed a surprising consistency with the measurements, and the maximum discrepancy between them was less than 9.4%.According to Li et al. [55], H 2 , CO, and small hydrocarbons have decisive impacts on flame propagation; hence, great consistency can be attained through the adoption of a detailed C1/H 2 /CO sub-mechanism in the current JFS mechanism.The concentration profiles of CH4 and C2H2 were well captured by the current mechanism, except for the slight overprediction for CH4 at the lean condition (Φ = 0.46), which implies that the mechanism is reliable in predicting the soot and PAH formulation in engine combustion.The decaying trend of O2 and the increasing trend of CO2 were both greatly reflected by the mechanism.As shown in Figure 23, the current mechanism also showed good performance in predicting the CO concentration except for a slight underestimation under the rich-and high-temperature condition (Figure 22b).Overall, the current mechanism exhibited a satisfactory prediction of the major species concentrations of the JFS surrogate, although the maximum deviation between the measurements and the predictions of CH4 was larger than the measured value.In fact, similar deviations were also observed even using a detailed mechanism, as reported by Malewicki et al. [14].

Laminar Flame Speed
The comparison results between the simulated and measured laminar flame speeds is shown in Figure 24.The laminar flame speeds of jet fuel were measured by Dooley et al. [13] at atmospheric pressure and the unburned mixture temperatures of 400 K and 470 K.The predicted values showed a surprising consistency with the measurements, and the maximum discrepancy between them was less than 9.4%.According to Li et al. [55], H2, CO, and small hydrocarbons have decisive impacts on flame propagation; hence, great consistency can be attained through the adoption of a detailed C1/H2/CO sub-mechanism in the current JFS mechanism.

Conclusions
In this study, a new four-component jet fuel surrogate (n-dodecane/iso-octane/isocetane/decalin/toluene) was formulated with a property optimizer.The chemical and physical properties of the target jet fuel were considered during the formulation of the fuel surrogate.Compared to the previous-existing jet surrogates, the current jet fuel surrogate showed better performance in emulating the real jet fuel properties.The deviations in CN, TSI, LHV, and density are −0.35,−1.29, 0.87, and 1.84, which are extremely in a low level.The property deviations of MW (10.73%), viscosity (5.88%), and surface tension (8.71%) are slightly larger, but these remain within acceptable ranges.
A compact reduced kinetic mechanism consisting of 124 species and 590 reactions was proposed for the newly formulated surrogate.This mechanism is highly suitable for practical engine combustion simulations due to its reduced size.The mechanism is composed of four sub-mechanisms covering toluene, decalin, n-dodecane, and isocetane.It was systematically developed, starting from the C5-C16 reactions and transitioning to the C2-C4 reactions, followed by the H 2 /O 2 /C1 reactions, PAH reactions, and NOx reactions.To validate the proposed mechanism, experimental data including ignition delay time (IDT), species concentrations, and laminar flame speeds were utilized, encompassing both the surrogate mixture and its individual components.The results demonstrated a favorable concordance between the proposed mechanism and the experimental data.The current mechanism could accurately capture the variations of the ignition delay time in the negative temperature coefficient (NTC) region.Consequently, the proposed mechanism emerges as a viable choice for effective utilization in practical engine simulations.

Figure 1 .
Figure 1.The physical and chemical processes in compression ignition engine.

Figure 1 .
Figure 1.The physical and chemical processes in compression ignition engine.

Figure 3 .
Figure 3.The flow chart of optimization process for the jet fuel surrogate.

Figure 3 .
Figure 3.The flow chart of optimization process for the jet fuel surrogate.

Figure 4 .
Figure 4.The major reaction pathways for the JFS kinetic mechanism.

Figure 4 .
Figure 4.The major reaction pathways for the JFS kinetic mechanism.

Figure 7 .
Figure 7.The major species profiles of premixed toluene flames under three equivalence ratios (0.75, 1.0, and 1.5) and a pressure of 10 atm.The determined results [38] are represented by symbols, while the simulated values from the current mechanism are depicted by lines.

Figure 7 .
Figure 7.The major species profiles of premixed toluene flames under three equivalence ratios (0.75, 1.0, and 1.5) and a pressure of 10 atm.The determined results [38] are represented by symbols, while the simulated values from the current mechanism are depicted by lines.

Figure 7 .
Figure 7.The major species profiles of premixed toluene flames under three equivalence ratios (0.75, 1.0, and 1.5) and a pressure of 10 atm.The determined results [38] are represented by symbols, while the simulated values from the current mechanism are depicted by lines.

Figure 8 .Figure 8 .
Figure 8.The concentrations of species resulting from the oxidation of toluene in the JSR are presented for various equivalence ratios at 0.6 s residence time and a pressure of 10 atm.The measured Figure 8.The concentrations of species resulting from the oxidation of toluene in the JSR are presented for various equivalence ratios at 0.6 s residence time and a pressure of 10 atm.The measured values [39] are denoted by symbols, while the simulation results obtained by the current mechanism are represented by lines.

Figure 9 .
Figure 9.The simulated and measured LFSs in counter-flow toluene flames.

Figure 10 .
Figure 10.The measured (symbols) and predicted IDTs (solid lines: current mechanism; dash lines: the detailed mechanism) for decalin in a shock tube.

Figure 10 .
Figure 10.The measured (symbols) and predicted IDTs (solid lines: current mechanism; d the detailed mechanism) for decalin in a shock tube.

Figure 12 .
Figure 12.The measured (symbols) and predicted (lines) decalin IDTs in a RCM at three equivalence ratios under 15 atm compression pressure.

Figure 13 .
Figure 13.The premixed decalin flames species profiles at various equivalence ratios.The measur values [46] are denoted by symbols, while the simulated values from the current mechanism represented by lines.

Figure 13 .
Figure 13.The premixed decalin flames species profiles at various equivalence ratios.The measured values [46] are denoted by symbols, while the simulated values from the current mechanism are represented by lines.

Figure 14 .
Figure 14.The measured (symbols[29]) and predicted (lines) concentration profiles of species resulting from the oxidation of decalin in the JSR at an equivalence ratio of 1.0 and a compression pressure of 10 atm.

Figure 15 .
Figure 15.The laminar flame speeds of the counter-flow decalin flame.The line stands for the simulations, and dots represent the experimental data [47].

Figure 14 .
Figure 14.The measured (symbols[29]) and predicted (lines) concentration profiles of species resulting from the oxidation of decalin in the JSR at an equivalence ratio of 1.0 and a compression pressure of 10 atm.

Figure 14 .
Figure14.The measured (symbols[29]) and predicted (lines) concentration profil sulting from the oxidation of decalin in the JSR at an equivalence ratio of 1.0 and pressure of 10 atm.

Figure 15 .
Figure 15.The laminar flame speeds of the counter-flow decalin flame.The line sta ulations, and dots represent the experimental data [47].

Figure 15 .
Figure 15.The laminar flame speeds of the counter-flow decalin flame.The line stands for the simulations, and dots represent the experimental data [47].

Figure 16 .
Figure 16.The measured (symbols) and predicted (lines) IDTs of n-dodecane in shock tu

Figure 17 .Figure 17 .
Figure 17.The measured (symbols [14]) and predicted (lines) mole fractions of the species during n-dodecane oxidation in JSR experiments at different equivalence ratios.

Figure 18 .Figure 18 .
Figure 18.The laminar flame speeds of the n-dodecane flame.The line stands for and dots represent the experimental data [51].

Figure 20 .
Figure 20.The measured (symbols[53])and predicted (lines) major species profiles of isocetan JSR experiments at different equivalence ratios.

Figure 21 .
Figure 21.The comparisons between the determined and simulated IDTs of jet fuel under differen (a) pressures and (b) equivalence ratios in a shock tube.

Figure 21 .
Figure 21.The comparisons between the determined and simulated IDTs of jet fuel under different (a) pressures and (b) equivalence ratios in a shock tube.

Figure 22 .
Figure22.The comparison between the determined IDTs of jet fuel[12] and the simulated the current mechanism, as well as the UM1 and UM2 mechanisms[15].

Figure 22 .
Figure22.The comparison between the determined IDTs of jet fuel[12] and the simulated IDTs by the current mechanism, as well as the UM1 and UM2 mechanisms[15].

Figure 23 .
Figure 23.The measured and predicted species profiles of jet fuel.Symbols are the experimental data [14]; lines are the predictions made by the current mechanism.

Figure 24 .
Figure24.The comparison between the measured and predicted laminar flame speeds.Symbols represent the experimental data[13], lines are the values predicted by the current mechanism.

Figure 23 .
Figure 23.The measured and predicted species profiles of jet fuel.Symbols are the experimental data [14]; lines are the predictions made by the current mechanism.

Sustainability 2023, 15 , 13792 22 of 25 Figure 23 .
Figure 23.The measured and predicted species profiles of jet fuel.Symbols are the experimental data [14]; lines are the predictions made by the current mechanism.

Figure 24 .
Figure24.The comparison between the measured and predicted laminar flame speeds.Symbols represent the experimental data[13], lines are the values predicted by the current mechanism.

Figure 24 .
Figure24.The comparison between the measured and predicted laminar flame speeds.Symbols represent the experimental data[13], lines are the values predicted by the current mechanism.

Table 2 .
The properties of the surrogate components.

Table 4 .
The comparison results among JFS and other jet fuel surrogates.