The work performed in this project includes the following aspects:
(1) The NOx formation mechanism in a 2D thermodiffusively unstable premixed hydrogen flame was investigated through DNS, and the performance of the flamelet model in predicting the NOx species in the thermodiffusively unstable premixed hydrogen flame were assessed through an a priori analysis;
(2) A large-scale 3D DNS of thermodiffusively unstable premixed hydrogen flames in a sufficiently large computational domain was performed, requiring about 67 Million CPU hours;
(3) The characteristic patterns of thermodiffusively unstable premixed hydrogen flames were quantified, including the global burning velocity, flame surface area and stretch factor;
(4) The NOx formation mechanism of thermodiffusively unstable premixed hydrogen flames was investigated through a reaction pathway analysis;
(5) A flamelet tabulation method was proposed to predict NOx formation in thermodiffusively unstable premixed hydrogen flames, accounting for curvature effects.
Based on the large-scale DNS dataset, we found that the global burning velocity in the 3D computational domain is about 70% higher than the 2D simulation, which is mainly related to the increased flame surface area. In particular, the range of positive curvature in the 3D simulation is much wider than in the 2D simulation. The peak concentration and production rate of H radicals in the 3D simulation are around two and five times higher than in the 2D simulation, respectively. Over 90% NO is formed in the positively-curved regions, with the NNH reaction pathway being dominant in the 3D thermodiffusively unstable premixed hydrogen flame. The thermal-NO reaction pathway is overall negligible for both 2D and 3D simulations. Compared to the conventional flamelet model, the radicals (e.g. H, N, NH, NNH) that are sensitive to the local curvature value can be accurately predicted by the new flamelet model. However, the prediction accuracy of NOx species mass fractions and their production rates did not significantly improve compared to the conventional flamelet model, which is due to their slow chemistry and the non-uniqueness of the flamelet table.