一般曲线坐标系下概率密度函数方法及其在超声速燃烧中的应用

doi:10.7527/S1000-6893.2022.26677

流体力学与飞行力学

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

一般曲线坐标系下概率密度函数方法及其在超声速燃烧中的应用

关清帝, 梁剑寒(), 张林, 陈文武, 陈玉俏

国防科技大学空天科学学院临近空间技术研究所，长沙 410073

收稿日期:2021-11-19 修回日期:2021-12-10 接受日期:2022-01-19 出版日期:2023-02-25 发布日期:2022-01-26
通讯作者: 梁剑寒 E-mail:jhleon@vip.sina.com
基金资助:
国家自然科学基金(91641201)

Probability density function method in general curvilinear coordinate system and its application in supersonic combustion

Qingdi GUAN, Jianhan LIANG(), Lin ZHANG, Wenwu CHEN, Yuqiao CHEN

Science and Technology on Scramjet Laboratory，College of Aerospace Science and Engineering，National University of Defense Technology，Changsha 410073，China

Received:2021-11-19 Revised:2021-12-10 Accepted:2022-01-19 Online:2023-02-25 Published:2022-01-26
Contact: Jianhan LIANG E-mail:jhleon@vip.sina.com
Supported by:
National Natural Science Foundation of China(91641201)

摘要/Abstract

摘要：

输运型概率密度函数（PDF）方法能够精确求解湍流/化学反应相互作用，有望为超燃冲压发动机中复杂的超声速湍流燃烧模拟提供可靠的计算模型。目前超声速流中的概率密度函数方法研究大多局限在简单构型下，严重限制了其在实际超声速燃烧问题中的应用。本文旨在发展基于多块结构网格的曲线坐标概率密度函数方法，使其能适用于具有复杂构型的实际超声速燃烧问题模拟。推导了曲线坐标下的概率密度输运方程和随机微分方程，建立了一般曲线坐标系下的概率密度函数方法及其数值求解框架。提出了多块结构网格下的曲线坐标粒子追踪方法，通过多块网格间的拓扑信息进行粒子局部坐标转换与通信，解决了曲线坐标下粒子局部坐标与全局坐标不一致的问题，成功实现了粒子在全局网格的连续追踪。在一维激波管和二维时间混合层中，通过与精确解和直接数值模拟结果的对比，初步验证了该方法和模型的准确性。进一步在二维圆柱绕流和三维球头绕流算例中，测试了其处理复杂构型的能力。最后开展了典型的超声速同轴射流火焰的数值模拟，与实验结果进行了定量对比。结果显示预测的反应组分与实验符合较好，全面验证了一般曲线坐标系下的概率密度函数方法在实际超声速燃烧问题中的适用性和准确性。

关键词: 概率密度函数方法, 大涡模拟, 粒子追踪, 曲线坐标, 超声速燃烧

Abstract:

The transported Probability Density Function （PDF） method can accurately solve turbulence/chemical reaction interactions， and can be used as a reliable computational model for complex supersonic turbulent combustion simulations in scramjet engines. At present， research on the PDF method for the supersonic flow is mostly confined to simple configuration， which seriously limits its application in practical supersonic combustion. The objective of this paper is to develop a curvilinear PDF method based on multi-block structured grids， which is applicable to the simulation of real supersonic combustion problems with complex configurations. The PDF equations and stochastic differential equations in curvilinear coordinates are derived， and the PDF method in a general curvilinear coordinate system and its numerical solution framework are established. A method for particle tracking in curvilinear coordinates on the multi-block structural grid is proposed. The topological information between multiple grid blocks is used for particle local coordinate transformation and communication， which solves the problem of inconsistency between particle local coordinates and global coordinates in curvilinear coordinates. Continuous tracking of particles on the global grid is successfully achieved. The accuracy of the curvilinear coordinate method and model is verified by comparing the results with the exact solution and direct numerical simulations in the one-dimensional shock tube and two-dimensional temporal mixing layer. The ability of the method to handle complex configurations is further tested in the simulations of the flow around a two-dimensional cylinder and a three-dimensional spherical blunt. Finally， numerical simulation of a typical supersonic coaxial jet flame is carried out， and the numerical results are quantitatively compared with the experimental results. The comparison results show that the predicted reaction components are in good agreement with those obtained in the experiment， which comprehensively verifies the applicability and accuracy of the PDF method in the general curvilinear coordinate system in actual supersonic combustion problems.

Key words: probability density function method, large eddy simulation, particle tracking, curvilinear coordinates, supersonic combustion

中图分类号:

V231.21

关清帝, 梁剑寒, 张林, 陈文武, 陈玉俏. 一般曲线坐标系下概率密度函数方法及其在超声速燃烧中的应用[J]. 航空学报, 2023, 44(4): 126677-126677.

Qingdi GUAN, Jianhan LIANG, Lin ZHANG, Wenwu CHEN, Yuqiao CHEN. Probability density function method in general curvilinear coordinate system and its application in supersonic combustion[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2023, 44(4): 126677-126677.

图/表 18

图 1

图 2

图 3

图 4

图 5

图 6

图 7

图 8

图 9

图 10

图 11

图 12

图 13

图 14

图 15

表 1

图 16

图 17

参考文献 31

1	杨越，游加平，孙明波. 超声速燃烧数值模拟中的湍流与化学反应相互作用模型［J］. 航空学报， 2015， 36（1）： 261-273.
	YANG Y， YOU J P， SUN M B. Modeling of turbulence-chemistry interactions in numerical simulations of supersonic combustion［J］. Acta Aeronautica et Astronautica Sinica， 2015， 36（1）： 261-273 （in Chinese）.
2	许爱国，单奕铭，陈锋，等. 燃烧多相流的介尺度动理学建模研究进展［J］. 航空学报， 2021， 42（12）： 625842.
	XU A G， SHAN Y M， CHEN F， et al. Progress of mesoscale modeling and investigation of combustion multiphase flow［J］. Acta Aeronautica et Astronautica Sinica， 2021， 42（12）： 625842 （in Chinese）.
3	PIERCE C D， MOIN P. Progress-variable approach for large-eddy simulation of non-premixed turbulent combustion［J］. Journal of Fluid Mechanics， 2004， 504： 73-97.
4	BERGLUND M， FEDINA E， FUREBY C， et al. Finite rate chemistry large-eddy simulation of self-ignition in supersonic combustion ramjet［J］. AIAA Journal， 2010， 48（3）： 540-550.
5	KERSTEIN A R. A linear-eddy model of turbulent scalar transport and mixing［J］. Combustion Science and Technology， 1988， 60（4-6）： 391-421.
6	陈崇沛，梁剑寒，关清帝，等. 守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用［J］. 航空学报， 2021， 42（S1）： 726364.
	CHEN C P， LIANG J H， GUAN Q D， et al. Conservative compressible one-dimensional turbulence method and its application in supersonic scalar mixing layer［J］. Acta Aeronautica et Astronautica Sinica， 2021， 42（S1）： 726364 （in Chinese）.
7	POPE S B. PDF methods for turbulent reactive flows［J］. Progress in Energy and Combustion Science， 1985， 11（2）： 119-192.
8	PANT T， JAIN U， WANG H F. Transported PDF modeling of compressible turbulent reactive flows by using the Eulerian Monte Carlo fields method［J］. Journal of Computational Physics， 2021， 425： 109899.
9	GIVI P. Model-free simulations of turbulent reactive flows［J］. Progress in Energy and Combustion Science， 1989， 15（1）： 1-107.
10	POPE S B. Computations of turbulent combustion： progress and challenges［J］. Symposium （International） on Combustion， 1991， 23（1）： 591-612.
11	GAO F， O'BRIEN E E. A large-eddy simulation scheme for turbulent reacting flows［J］. Physics of Fluids A： Fluid Dynamics， 1993， 5（6）： 1282-1284.
12	JAISHREE J. Lagrangian and Eulerian probability density function methods for turbulent reacting flows ［D］. State College： Pennsylvania State University， 2011.
13	YANG Y， WANG H F， POPE S B， et al. Large-eddy simulation/probability density function modeling of a non-premixed CO/H₂ temporally evolving jet flame［J］. Proceedings of the Combustion Institute， 2013， 34（1）： 1241-1249.
14	JABERI F A， COLUCCI P J， JAMES S， et al. Filtered mass density function for large-eddy simulation of turbulent reacting flows ［J］. Journal of Fluid Mechanics， 1999， 401： 85-121.
15	BANAEIZADEH A， LI Z R， JABERI F A. Compressible scalar filtered mass density function model for high-speed turbulent flows［J］. AIAA Journal， 2011， 49（10）： 2130-2143.
16	KOMPERDA J， GHIASI Z， LI D R， et al. Simulation of the cold flow in a ramp-cavity combustor using a DSEM-LES/FMDF hybrid scheme［C］∥54th AIAA Aerospace Sciences Meeting. Reston： AIAA， 2016.
17	ZHANG L， LIANG J H， SUN M B， et al. An energy-consistency-preserving large eddy simulation-scalar filtered mass density function （LES-SFMDF） method for high-speed flows［J］. Combustion Theory and Modelling， 2018， 22（1）： 1-37.
18	ZHANG L， LIANG J H， SUN M B， et al. A conservative and consistent scalar filtered mass density function method for supersonic flows［J］. Physics of Fluids， 2021， 33（2）： 026101.
19	VALIDI A， SCHOCK H， JABERI F. Turbulent jet ignition assisted combustion in a rapid compression machine［J］. Combustion and Flame， 2017， 186： 65-82.
20	RANADIVE H D. High-order compressible formulation for LES/PDF simulations of turbulent reacting flows［D］. Sydney： University of New South Wales， 2019.
21	CHENG T S， WEHRMEYER J A， PITZ R W， et al. Raman measurement of mixing and finite-rate chemistry in a supersonic hydrogen-air diffusion flame［J］. Combustion and Flame， 1994， 99（1）： 157-173.
22	ABDULRAHMAN H， VALIDI A， JABERI F. Large-eddy simulation/filtered mass density function of non-premixed and premixed colorless distributed combustion［J］. Physics of Fluids， 2021， 33（5）： 055118.
23	YOSHIZAWA A， HORIUTI K. A statistically-derived subgrid-scale kinetic energy model for the large-eddy simulation of turbulent flows［J］. Journal of the Physical Society of Japan， 1985， 54（8）： 2834-2839.
24	LI P B， LI C Y， WANG H B， et al. Distribution characteristics and mixing mechanism of a liquid jet injected into a cavity-based supersonic combustor［J］. Aerospace Science and Technology， 2019， 94： 105401.
25	LIU C Y， SUN M B， WANG H B， et al. Ignition and flame stabilization characteristics in an ethylene-fueled scramjet combustor［J］. Aerospace Science and Technology， 2020， 106： 106186.
26	SHU C W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws： NASA/CR-97-206253［R］. Washington， D.C.： NASA， 1997.
27	GIKHMAN I I， SKOROKHOD A V. The theory of stochastic processes III［M］. Berlin， Heidelberg： Springer， 2007.
28	WANG H F， POPOV P P， POPE S B. Weak second-order splitting schemes for Lagrangian Monte Carlo particle methods for the composition PDF/FDF transport equations［J］. Journal of Computational Physics， 2010， 229（5）： 1852-1878.
29	EVANS J， SCHEXNAYDER C J， BEACH H. Application of a two-dimensional parabolic computer program to prediction of turbulent reacting flows： NASA-TP-1169［R］. Washington， D.C.： NASA， 1978.
30	LIU C Y， WANG N， YANG K， et al. Large eddy simulation of a supersonic lifted jet flame in the high-enthalpy coflows［J］. Acta Astronautica， 2021， 183： 233-243.
31	CONAIRE M Ó， CURRAN H J， SIMMIE J M， et al. A comprehensive modeling study of hydrogen oxidation［J］. International Journal of Chemical Kinetics， 2004， 36（11）： 603-622.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

实验参数	数值
实验参数	内射流	外射流
马赫数	2.0	1.9
温度/K	251	1 495
速度/（m∙s^-1）	2 432	1 510
压力/kPa	100	100
H₂质量分数	1.0	0
O₂质量分数	0	0.261
H₂O质量分数	0	0.271
N₂质量分数	0	0.468

[1]	乔竑玮, 梁剑寒, 张林, 孙明波, 陈玉俏. 超声速燃烧中的概率密度函数方法研究进展[J]. 航空学报, 2024, 45(8): 28802-028802.
[2]	刘小勇, 王明福, 刘建文, 任鑫, 张轩. 超燃冲压发动机研究回顾与展望[J]. 航空学报, 2024, 45(5): 529878-529878.
[3]	纪鉴恒, 蔡尊, 王泰宇, 孙明波, 王振国. 宽速域超燃冲压发动机流动燃烧过程研究进展[J]. 航空学报, 2024, 45(3): 28696-028696.
[4]	张育鸣, 戴玉婷, 黄广靖, 杨超, 蒋树杰. 柔性变形后缘翼型阵风减缓及气动噪声分析[J]. 航空学报, 2024, 45(10): 129219-129219.
[5]	张子佩, 赵钟, 陈坚强, 刘健, 邓小兵. 风雷软件LES开发设计与验证[J]. 航空学报, 2023, 44(6): 127171-127171.
[6]	施方成, 高振勋, 田雨岩, 蒋崇文, 王田天, 李椿萱. 超声速理想膨胀喷流噪声的大涡模拟[J]. 航空学报, 2023, 44(2): 626266-626266.
[7]	孙明波, 朱家健, 罗天罡, 李沁远, 田轶夫, 万明罡, 孙永超. 非稳态超声速燃烧电激励调控技术研究进展[J]. 航空学报, 2023, 44(15): 528787-528787.
[8]	赵宾宾, 张恒, 李杰. 翼型结冰状态复杂分离流动数值模拟综述[J]. 航空学报, 2023, 44(1): 627211-627211.
[9]	张钧铎, 左青海, 林孟达, 黄伟希, 潘卫军, 崔桂香. ARJ21飞机尾涡在侧风条件下的近地演化数值模拟[J]. 航空学报, 2022, 43(5): 125043-125043.
[10]	谢晨月, 王建春, 万敏平, 陈十一. 基于人工神经网络的可压缩湍流大涡模拟模型[J]. 航空学报, 2021, 42(9): 625723-625723.
[11]	杨昌和, 李彦达, 张江, 王昉, 袁晓如. 高性能流场并行粒子追踪数据管理系统[J]. 航空学报, 2021, 42(9): 625748-625748.
[12]	王萍, 郑晓静. 风沙两相流数值模拟研究进展[J]. 航空学报, 2021, 42(9): 625767-625767.
[13]	王德鑫, 褚佑彪, 刘难生, 李祝飞, 杨基明. 高背压进气道中内外流耦合作用的大涡模拟[J]. 航空学报, 2021, 42(9): 625754-625754.
[14]	刘远强, 王艳冰, 项松, 王梦琪. 面向电动飞机应用的不同叶尖螺旋桨噪声特性[J]. 航空学报, 2021, 42(3): 624567-624567.
[15]	王方, 王煜栋, 姜胜利, 陈军, 唐军, 徐华胜, 李象远, 邢竞文, 高东硕, 金捷. AECSC-JASMIN湍流燃烧仿真软件研发和检验[J]. 航空学报, 2021, 42(12): 625003-625003.

一般曲线坐标系下概率密度函数方法及其在超声速燃烧中的应用

Probability density function method in general curvilinear coordinate system and its application in supersonic combustion

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 18

参考文献 31

相关文章 15

编辑推荐

Metrics

本文评价