航空发动机及燃气轮机等动力装备燃烧室广泛采用的贫燃燃烧方式经常遇到破坏性的非定常热声耦合振荡燃烧问题。非定常振荡燃烧数值预测是一个长期的研究热点和难题。发展了针对振荡燃烧的耦合直接求解数值模拟方法，包括优化的动态模型参数的高精度自适应湍流模型（SATES），耦合可压缩的详细化学反应建表FGM燃烧模型。选取的3种湍流燃烧模型包括有限速率模型（W1）及火焰面密度封闭方法中的Zimont（W2）和Fureby（W3）2种褶皱因子模型。针对经典的LIMOUSINE燃烧室多个部分预混振荡燃烧工况开展了数值研究，发现自适应湍流模拟框架下的3种燃烧模型均准确预测到了振荡燃烧的振荡频率，与试验相比，误差<6.4%；对于振荡燃烧压力脉动振幅的预测结果，有限速率模型（W1）和Zimont（W2）模型结果显著大于试验值，误差>380%；Fureby（W3）模型结果与试验值吻合较好，误差<17.9%。表明振荡燃烧的数值预测对不同的湍流及湍流燃烧模型具有较高的敏感性。不同的工况结果表明，振荡燃烧存在完全振荡模态和过渡模态，完全振荡模态中数值预测的特征主频在燃烧室上下游多个位置趋于一致；过渡振荡模态中上下游不同位置预测的特征主频存在双峰，即双频模式。本文发展的数值模拟方法对振荡燃烧的预测具有较高的计算精度和可靠性。

The lean-fuel combustion method is widely used in the combustion chamber of aircraft engines and gas turbines. However， this method often encounters destructive unsteady thermoacoustic combustion instability problems. Numerical prediction of unsteady combustion instability is a long-term research hotspot and challenging problem. In this study， a coupled direct numerical simulation method for combustion instability is developed， coupling a high-fidelity Self-Adaptive Turbulence Eddy Simulation （SATES） method with optimized dynamic model parameters and a compressible tabulated FGM combustion model with detailed chemical reaction. Three turbulent combustion models are selected including the finite rate model （W1） and two Flame Surface Density （FSD） models of Zimont （W2） and Fureby （W3）. Numerical studies were carried out for the classical partially-premixed LIMOUSINE combustion chamber with combustion instability. It is found that the three combustion models under the SATES framework accurately predict the instability frequency of oscillating combustion， and the difference from the experimental data is small， which less than 6.4%. As to the prediction results of the instability amplitude of oscillating pressure， the finite rate model （W1） and Zimont （W2） model results are significantly larger than the experimental data， with the difference larger than 380%， and the Fureby （W3） model results are in good agreement with the experiments with the difference less than 17.9%. It is shown that the numerical prediction of combustion instability has high sensitivity to different turbulence and turbulent combustion models. The results from different test cases show that there are complete instability modes and transition modes in the combustion instability. The dominant frequency predicted by the numerical simulations in the complete instability mode tends to be consistent at multiple positions upstream and downstream of the combustion chamber. There are bimodal peaks in the predicted peak frequencies at different positions upstream and downstream in the transition instability mode， that is， dual-frequency mode. These results demonstrate that the numerical simulation method developed in the present study has high computational accuracy and reliability for the prediction of combustion instability.