A petascale direct numerical simulation study of the modelling of flame wrinkling for large-eddy simulations in intense turbulence

详细信息	查看全文 | 推荐本文 |

• 作者：Evatt R. Hawkes^a ; ^b ; ^{evatt.hawkes@unsw.edu.au} ; Obulesu Chatakonda^a ; Hemanth Kolla^c ; Alan R. Kerstein^c ; Jacqueline H. Chen^c
• 关键词：Direct numerical simulation ; Premixed combustion ; Large eddy simulation ; Flame speed ; Fractal
• 刊名：Combustion and Flame
• 出版年：2012
• 期刊代码：21_00102180
• 类别：ce
• 出版时间：August, 2012
• 卷：159
• 期：8
• 页码：2690-2703
• 文件大小：1165 K

摘要

A new set of petascale direct numerical simulations (DNS) modelling lean hydrogen combustion with detailed chemistry in a temporally evolving slot-jet configuration is presented as a database for the development and validation of models for premixed turbulent combustion. The jet Reynolds number is 10,000, requiring grid numbers up to nearly seven billion, which was achieved by computation on 120,000 processor cores. In contrast to many prior DNS studies, a mean shear exists that drives strong turbulent mixing within the flame structure. Three cases are simulated with different Damk枚hler numbers, while Reynolds number is held fixed. Basic statistics are presented showing that integrated burning rates up to approximately six times the laminar burning rate are obtained. It is shown that increased flame surface area accounts for most of the enhanced burning while increases in the burning rate per unit area also play an important contribution.

The database is then used to assess a new model of flame wrinkling intended for large-eddy simulations (LES). The approach draws on concepts from fractal geometry, requiring the modelling of an inner cut-off scale representing the smallest scale of flame wrinkling, and the fractal dimension controlling the resolution dependence of the unresolved flame surface area. In contrast to previous modelling, it is argued that the inner cut-off should be filter-size invariant in an inertial range. Then, dimensional and physical reasoning together with Damk枚hler鈥檚 limiting scaling laws for the turbulent flame speed are used to infer the cut-off and fractal dimension in limiting regimes. Two methods of determining the fractal dimension are proposed: a static, algebraic expression or a dynamic approach exploiting a Germano-type identity. Finally the model is compared against the DNS in a priori tests and is found to give excellent results, quantitatively capturing the trends with time, space, filter size and Damk枚hler number.