基于湍动能谱的自适应大涡模拟方法

推进技术 ›› 2013, Vol. 34 ›› Issue (9): 1214-1221.

基于湍动能谱的自适应大涡模拟方法

单繁立,朴英

清华大学航天航空学院,北京100084;清华大学航天航空学院,北京100084

发布日期:2021-08-15
作者简介:单繁立(1983—),男,博士生,研究领域为航空宇航推进系统、湍流燃烧。E-mail:cfl02@mails.tsinghua.edu.cn

Adaptive Large Eddy Simulation Based on Energy Spectrum

School of Aerospace,Tsinghua University,Beijing 100084,China;School of Aerospace,Tsinghua University,Beijing 100084,China

Published:2021-08-15

摘要/Abstract

摘要： 为了在计算量有限的前提下实现对不同尺度上湍流燃烧结构的直接求解,提出了新的基于湍动能谱的自适应大涡模拟的概念。在这一概念下,利用湍动能谱和自适应网格技术建立了跨尺度算法,从而实现DNS和LES两种方法的耦合。在湍流非预混氢气射流火焰的模拟中,该方法既能捕捉大尺度湍流结构,又能分辨小尺度燃烧结构,计算出的物理量的瞬态和时均分布和DNS结果吻合得较好。计算偏差主要体现在DNS和LES区域的过渡处,这与亚格子湍动能在过渡处的变化不够平滑有关,表明跨尺度算法需要优化,以完善对不同尺度间湍动能传输的模拟。 

关键词: 大涡模拟;自适应;湍动能谱;跨尺度算法;湍流燃烧 

Abstract: In order to directly resolve turbulent combustion at different scales with acceptable computational expenditure, a new method called Energy Spectrum Based Adaptive Large Eddy Simulation was proposed. A trans-scale algorithm based on energy spectrum and dynamic adaptive mesh technology was developed for this method to couple DNS and LES. It is shown in the simulation of a turbulent nonpremixed hydrogen jet flame that this method can resolve the turbulent coherent structures at large scales as well as the combustion structures at small scales. It is validated with good agreement on transient and temporal-average distribution of physical quantities between this method and DNS. The computational deviations of this method in the transition region between DNS and LES regimes are greater than those in other places. This is because the variation of sub-grid turbulent kinetic energy in the transition region is not sufficient smooth, indicating that the trans-scale algorithm has to be optimized for a better simulation of turbulent kinetic energy cascade between different scales. 

Key words: Large eddy simulation; Adaptive; Energy spectrum; Trans-scale algorithm; Turbulent combustion

单繁立,朴英. 基于湍动能谱的自适应大涡模拟方法[J]. 推进技术, 2013, 34(9): 1214-1221.

SHAN Fan-li,PIAO Ying. Adaptive Large Eddy Simulation Based on Energy Spectrum[J]. Journal of Propulsion Technology, 2013, 34(9): 1214-1221.

参考文献

[1] 曹红军, 张会强, 王兵,等.基于有限反应速率的扩散燃烧大涡模拟 [J].工程热物理学报, 2010,1(5) :871-874.
[2] Yoo C S, Lu T, Chen J H, et al.Direct Numerical Simulations of Ignition of a Lean n-Heptane/Air Mixture with Temperature Inhomogeneities at Constant Volume:Parametric Study [J].Combustion and Flame, 2011,8:1727-1741.
[3] Chen J H.Petascale Direct Numerical Simulation of Turbulent Combustion—Fundamental Insights Towards Predictive Models [J].Proceedings of the Combustion Institute, 2011,3:99-123.
[4] Chakravarthy V, Menon S.Subgrid Modeling of Premixed Flames in the Flamelet Regime [J].Flow, Turbulence and Combustion, 2000,5:23-45.
[5] El-Asrag H, Menon S.Large Eddy Simulation of Bluff-Body Stabilized Swirling Non-Premixed Flames [J].Proceedings of the Combustion Institute, 2007,1:1747-1754.
[6] Duan L.DNS of Hypersonic Turbulent Boundary Layers [D].Princeton:Princeton University, 2011.
[7] Duan L, Feldick A M, Martin M P, et al.Study of Turbulence-Radiation Interaction in Hypersonic Turbulent Boundary Layers [R].AIAA 2011-749.
[8] 单繁立, 朴英.湍流燃烧的动态自适应直接数值模拟研究 [J].推进技术, 2013,4(6):786-794.(SHAN Fan-li,PIAO Ying.Direct Numerical Simulation of Turbulent Combustion with Dynamic Adaptive Mesh Technology[J].Journal of Propulsion Technology,2013,4(6):786-794.)
[9] Kim W, Menon S.An Unsteady Incompressible Navier-Stokes Solver for Large Eddy Simulation of Turbulent Flows [J].International Journal for Numerical Methods in Fluids, 1999,1:983-1017.
[10] El-Asrag H.Large Eddy Simulation Subgrid Model for Soot Prediction [D].Atlanta:Georgia Institute of Technology, 2007.
[11] Liu S, Meneveau C, Katz J.On the Properties of Similarity Subgrid-Scale Models as Deduced from Measurements in a Turbulent Jet [J].Journal of Fluid Mechanics, 1994,5:83-119.
[12] Lilly D K.A Proposed Modification of the Germano Subgrid-Scale Closure Method [J].Physics of Fluids, 1992,4(3):633-635.
[13] Toro E F.Riemann Solvers and Numerical Methods for Fluid Dynamics [M].Berlin:Springer, 1999.
[14] Barth T J, Jespersen D.The Design and Application of Upwind Schemes on Unstructured Meshes [R].AIAA 89-0366.
[15] Williamson J H.Low-Storage Runge-Kutta Schemes [J].Journal of Computational Physics, 1980,5:48-56.
[16] Chen Z.Studies on the Initiation, Propagation, and Extinction of Premixed Flames [D].Princeton:Princeton University, 2009.
[17] Brown P N, Byrne G D, Hindmarsh A C.VODE:A Variable-Coefficient ODE Solver [J].SIAM Journal on Scientific and Statistical Computing, 1989,0(5):1038-1051.
[18] Sun M.Numerical and Experimental Studies of Shock Wave Interaction with Bodies [D].Sendai:Tohoku University, 1998.
[19] Pope S B.Turbulent Flows [M].New York:Cambridge University Press, 2000.
[20] James S, Jaberi F A.Large Scale Simulations of Two-Dimensional Nonpremixed Methane Jet Flames [J].Combustion and Flame, 2000,3:465-487.
[21] Connaire M O, Curran H J, Simmie J M, et al.A Comprehensive Modeling Study of Hydrogen Oxidation [J].International Journal of Chemical Kinetics, 2004,6:603-622.
[22] 杨武兵.湍流非预混燃烧简化的PDF模型的大涡模拟研究 [D].北京:清华大学, 2006.
[23] Poinsot T J, Lele S K.Boundary Conditions for Direct Simulations of Compressible Viscous Flows [J].Journal of Computational Physics, 1992,1:104-129.