CN112685893B - 激波诱燃冲压发动机燃烧室设计方法 - Google Patents

激波诱燃冲压发动机燃烧室设计方法 Download PDF

Info

Publication number
CN112685893B
CN112685893B CN202011594471.1A CN202011594471A CN112685893B CN 112685893 B CN112685893 B CN 112685893B CN 202011594471 A CN202011594471 A CN 202011594471A CN 112685893 B CN112685893 B CN 112685893B
Authority
CN
China
Prior art keywords
angle
wedge
oblique
combustion chamber
curve
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011594471.1A
Other languages
English (en)
Other versions
CN112685893A (zh
Inventor
覃建秀
许亮
杨武兵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN202011594471.1A priority Critical patent/CN112685893B/zh
Publication of CN112685893A publication Critical patent/CN112685893A/zh
Application granted granted Critical
Publication of CN112685893B publication Critical patent/CN112685893B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Combustion Methods Of Internal-Combustion Engines (AREA)

Abstract

本发明公开了一种激波诱燃冲压发动机燃烧室设计方法,方法包括如下步骤:根据燃烧室入口来流参数获取斜激波的激波角βOSW‑斜劈角度θ曲线、斜爆震波的爆震波角βODW‑斜劈角度θ曲线及斜劈角度θ‑斜激波后温度T曲线、斜劈角度θ‑爆震波后温度T曲线;根据爆震波角βODW‑斜劈角度θ曲线和斜劈角度θ‑斜激波后温度T曲线确定斜劈角度范围;在斜劈角度范围内选取多个斜劈角度,对多个斜劈角度分别行燃烧室流场计算;根据流场计算结果,结合燃烧室尺寸约束确定斜劈角度θ0;根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw。本发明提供了燃烧室斜劈角度和长度设计方法,确保燃烧室中的斜爆震波能够实现驻定和起爆。

Description

激波诱燃冲压发动机燃烧室设计方法
技术领域
本发明属于发动机技术领域,更具体地,涉及一种激波诱燃冲压发动机燃烧室设计方法。
背景技术
激波诱燃冲压发动机(shcramjet)是未来吸气式高超声速飞行器最理想的动力系统之一,该发动机在进气道进行喷注和混合,以混合气形式进入燃烧室,通过斜激波诱发燃烧,能有效弥补了传统超燃冲压发动机机身一体化设计带来的缺点,缩短燃烧室长度,降低机身结构重量,扩宽飞行马赫数范围。激波诱燃冲压发动机燃烧室如图1所示,该发动机在燃烧室中采用斜劈来诱发斜激波,从而诱导斜爆震波。斜爆震波起爆和驻定是其关键技术之一,而斜劈的角度和长度是斜爆震波起爆和驻定的关键影响因素,斜劈角度过大会导致斜爆震波无法驻定,而过小的斜劈角度或过短的斜劈长度无法起爆斜爆震波,过长的斜劈长度会导致燃烧室结构重量增加,降低发动机的推重比,所以,在燃烧室设计中,关键是斜劈角度和长度的选取。
因此,如何根据来流条件设计斜劈角度和长度是本领域亟待解决的问题。
发明内容
本发明的目的是提供一种激波诱燃冲压发动机燃烧室设计方法,以根据来流条件来确定激波诱燃冲压发动机燃烧室的斜劈角度和长度。
为了实现上述目的,本发明提供一种激波诱燃冲压发动机燃烧室设计方法,所述方法包括如下步骤:
根据燃烧室入口来流参数获取斜激波的激波角βOSW-斜劈角度θ曲线、斜爆震波的爆震波角βODW-斜劈角度θ曲线及斜劈角度θ-斜激波后温度T曲线、斜劈角度θ-爆震波后温度T曲线;
根据爆震波角βODW-斜劈角度θ曲线和斜劈角度θ-斜激波后温度T曲线确定斜劈角度范围;
在斜劈角度范围内选取多个斜劈角度,对多个斜劈角度分别行燃烧室流场计算;
根据流场计算结果,结合燃烧室尺寸约束确定斜劈角度θ0
根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw
可选地,所述燃烧室入口来流参数包括来流马赫数Ma0、来流温度T0和来流压力P0
可选地,所述斜劈角度范围为(θCJ,θT=1600),其中,θCJ为C-J爆震对应的斜劈角度,θT=1600为T=1600K对应的斜劈角度。
可选地,在斜劈角度范围内至少选取三个斜劈角度,并分别行燃烧室流场计算。
可选地,所述临界斜劈长度Lw通过如下公式获得:
其中,Ma1为斜劈尾部附近的马赫数,L1为斜激波-斜爆震波转变点在水平方面投影距离,β为斜激波角或斜爆震波角;
可选地,所述激波角βOSW(爆震波角βODW)-斜劈角度θ曲线和斜劈角度θ-斜激波(爆震波)后温度T曲线由斜激波(斜爆震波)关系式获得,所述关系式为:
其中,Ma0n=Ma0sinβ,Ma0为燃烧室入口来流马赫数,T0为来流温度,γ为来流混合气的比热比,Cp为来流混合气的定压比热,Q为混合气燃烧放热量,β为斜激波角或斜爆震波角。
可选地,所述根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw包括:根据斜劈角度θ0的流场计算结果,考虑膨胀波的影响,建立斜爆震波起爆准则。
可选地,所述燃烧室的斜劈长度L>Lw
本发明的有益效果在于:
1)本发明给出了激波诱燃冲压发动机燃烧室的斜劈角度和长度设计方法,确保燃烧室的斜爆震波能够驻定和起爆,并最大程度减小燃烧室的结构重量,提高发动机的推重比。
2)本发明在驻定斜劈角度范围内缩小了斜劈角度的选取范围。
3)本发明能够根据来流条件和斜劈角度,给出了临界斜劈长度。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1为激波诱燃冲压发动机燃烧室的示意图。
图2为临界斜劈长度的示意图。
图3为激波角βOSW(爆震波角βODW)-斜劈角度θ曲线图。
图4为斜劈角度θ-波后温度T曲线图。
图5为不同斜劈角度下的温度流场图(单位m)。
图6为不同斜劈长度下的温度和H2O流场图(黑色等值线代表H2O)。
图7为激波诱燃冲压发动机燃烧室的设计流程图。
附图标记说明
1、燃烧室;2、斜劈。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
本发明公开了一种激波诱燃冲压发动机燃烧室设计方法,方法包括如下步骤:
根据燃烧室入口来流参数获取斜激波的激波角βOSW-斜劈角度θ曲线、斜爆震波的爆震波角βODW-斜劈角度θ曲线及斜劈角度θ-斜激波后温度T曲线、斜劈角度θ-爆震波后温度T曲线;
根据爆震波角βODW-斜劈角度θ曲线和斜劈角度θ-斜激波后温度T曲线确定斜劈角度范围;
在斜劈角度范围内选取多个斜劈角度,对多个斜劈角度分别行燃烧室流场计算;
根据流场计算结果,结合燃烧室尺寸约束确定斜劈角度θ0
根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw
具体地,本发明给出了激波诱燃冲压发动机燃烧室的斜劈角度和长度设计方法,确保燃烧室的斜爆震波能够驻定和起爆,并最大程度减小燃烧室的结构重量,提高发动机的推重比。本发明在驻定斜劈角度范围内缩小了斜劈角度的选取范围。本发明能够根据来流条件和斜劈角度,给出了临界斜劈长度。
作为可选方案,燃烧室入口来流参数包括来流马赫数Ma0、来流温度T0和来流压力P0
作为可选方案,斜劈角度范围为(θCJ,θT=1600),其中,θCJ为C-J爆震对应的斜劈角度,θT=1600为T=1600K对应的斜劈角度。
作为可选方案,在斜劈角度范围内至少选取三个斜劈角度,并分别行燃烧室流场计算。
作为可选方案,临界斜劈长度Lw通过如下公式获得:
其中,Ma1为斜劈尾部附近的马赫数,L1为斜激波-斜爆震波转变点在水平方面投影距离,β为斜激波角或斜爆震波角;
作为可选方案,激波角βOSW(爆震波角βODW)-斜劈角度θ曲线和斜劈角度θ-斜激波(爆震波)后温度T曲线由斜激波(斜爆震波)关系式获得,关系式为:
其中,Ma0n=Ma0sinβ,Ma0为燃烧室入口来流马赫数,T0为来流温度,γ为来流混合气的比热比,Cp为来流混合气的定压比热,Q为混合气燃烧放热量,β为斜激波角或斜爆震波角。
作为可选方案,根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw包括:根据斜劈角度θ0的流场计算结果,考虑膨胀波的影响,建立斜爆震波起爆准则。
作为可选方案,燃烧室的斜劈长度L>Lw
具体地,当斜劈长度L>Lw,斜爆震波能成功起爆。
实施例
激波诱燃冲压发动机燃烧室如图1所示,包括燃烧室1和斜劈2。
本实施例中,燃烧室入口来流参数包括P0=16.261kPa,T0=721.8K,Ma0=6.0。
(1)假设斜激波(斜爆震波)前后混合气的比热和比热比不变,根据斜激波(斜爆震波)关系式获得斜激波的激波角βOSW-斜劈角度θ曲线、斜爆震波的爆震波角βODW-斜劈角度θ曲线,如图3所示;
(2)根据斜激波(斜爆震波)关系式获得斜劈角度θ-斜激波后温度T曲线、斜劈角度θ-爆震波后温度T曲线,如图4所示;
(3)根据爆震波角βODW-斜劈角度θ曲线和斜劈角度θ-斜激波后温度T曲线确定斜劈角度范围(θCJ,θT=1600);
(4)在斜劈角度范围内选取三个斜劈角度θ=14°、16°、18°,如图4所示;
(5)假设斜劈为无限长,采用Fluent软件针对三个斜劈角度分别进行流场计算,结果如图5所示,根据燃烧室尺寸约束(燃烧室长度小于12cm),选取θ0=18°;
(6)根据图5中θ0=18°的流场结果,Ma1=3.4,L1=0.079m,βOSW=26°,θ=18°,考虑膨胀波的影响,建立斜爆震波起爆准则,可得Lw=0.047,因此燃烧室斜劈长度L大于临界斜劈长度Lw=0.047m时,斜爆震波能成功起爆,在临界斜劈长度时,波系结构如图2所示;;
(7)针对0.05m、0.047m、0.044m三个不同斜劈长度,计算温度和H2O流场图(黑色等值线代表H2O),计算结果如图6所示,证实了斜劈长度大于0.047m时,斜爆震波才能成功起爆。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (7)

1.一种激波诱燃冲压发动机燃烧室设计方法,其特征在于,所述方法包括如下步骤:
根据燃烧室入口来流参数获取斜激波的激波角βOSW-斜劈角度θ曲线、斜爆震波的爆震波角βODW-斜劈角度θ曲线及斜劈角度θ-斜激波后温度T曲线、斜劈角度θ-爆震波后温度T曲线;
根据爆震波角βODW-斜劈角度θ曲线和斜劈角度θ-斜激波后温度T曲线确定斜劈角度范围;
在斜劈角度范围内选取多个斜劈角度,对多个斜劈角度分别行燃烧室流场计算;
根据流场计算结果,结合燃烧室尺寸约束确定斜劈角度θ0
根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw
所述斜劈角度范围为(θCJ,θT=1600),其中,θCJ为C-J爆震对应的斜劈角度,θT=1600为T=1600K对应的斜劈角度。
2.根据权利要求1所述的激波诱燃冲压发动机燃烧室设计方法,其特征在于,所述燃烧室入口来流参数包括来流马赫数Ma0、来流温度T0和来流压力P0
3.根据权利要求1所述的激波诱燃冲压发动机燃烧室设计方法,其特征在于,在斜劈角度范围内至少选取三个斜劈角度,并分别行燃烧室流场计算。
4.根据权利要求1所述的激波诱燃冲压发动机燃烧室设计方法,其特征在于,所述临界斜劈长度Lw通过如下公式获得:
其中,Ma1为斜劈尾部附近的马赫数,L1为斜激波-斜爆震波转变点在水平方面投影距离,β为斜激波角或斜爆震波角;
5.根据权利要求1所述的激波诱燃冲压发动机燃烧室设计方法,其特征在于,所述激波角βOSW-斜劈角度θ曲线和斜劈角度θ-斜激波后温度T曲线由斜激波关系式获得,所述爆震波角βODW-斜劈角度θ曲线和斜劈角度θ-爆震波后温度T曲线由斜爆震波关系式获得,所述关系式为:
其中,Ma0n=Ma0sinβ,Ma0为燃烧室入口来流马赫数,T0为来流温度,γ为来流混合气的比热比,Cp为来流混合气的定压比热,Q为混合气燃烧放热量,β为斜激波角或斜爆震波角。
6.根据权利要求1所述的激波诱燃冲压发动机燃烧室设计方法,其特征在于,所述根据斜劈角度θ0的流场计算结果,确定临界斜劈长度Lw包括:根据斜劈角度θ0的流场计算结果,考虑膨胀波的影响,建立斜爆震波起爆准则。
7.根据权利要求1所述的激波诱燃冲压发动机燃烧室设计方法,其特征在于,所述燃烧室的斜劈长度L>Lw
CN202011594471.1A 2020-12-29 2020-12-29 激波诱燃冲压发动机燃烧室设计方法 Active CN112685893B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011594471.1A CN112685893B (zh) 2020-12-29 2020-12-29 激波诱燃冲压发动机燃烧室设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011594471.1A CN112685893B (zh) 2020-12-29 2020-12-29 激波诱燃冲压发动机燃烧室设计方法

Publications (2)

Publication Number Publication Date
CN112685893A CN112685893A (zh) 2021-04-20
CN112685893B true CN112685893B (zh) 2023-08-01

Family

ID=75455167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011594471.1A Active CN112685893B (zh) 2020-12-29 2020-12-29 激波诱燃冲压发动机燃烧室设计方法

Country Status (1)

Country Link
CN (1) CN112685893B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117688697B (zh) * 2024-02-02 2024-04-26 中国人民解放军空军工程大学 旋转爆震发动机进气道设计方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110516310A (zh) * 2019-07-31 2019-11-29 中国空气动力研究与发展中心 旋转爆震反压的非定常数值模拟方法
CN111207009A (zh) * 2019-12-26 2020-05-29 中国空气动力研究与发展中心 利用外加瞬时能量源在超声速气流中起爆斜爆震波的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110516310A (zh) * 2019-07-31 2019-11-29 中国空气动力研究与发展中心 旋转爆震反压的非定常数值模拟方法
CN111207009A (zh) * 2019-12-26 2020-05-29 中国空气动力研究与发展中心 利用外加瞬时能量源在超声速气流中起爆斜爆震波的方法

Also Published As

Publication number Publication date
CN112685893A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
US6347509B1 (en) Pulsed detonation engine with ejector bypass
CN101881238B (zh) 一种吸气式脉冲爆震发动机及其起爆方法
CN112685893B (zh) 激波诱燃冲压发动机燃烧室设计方法
CN101975122B (zh) 带有磁流体能量旁路系统的驻定爆震发动机
CN108488004B (zh) 一种基于可变斜楔角的驻定爆震发动机
CN106837603B (zh) 一种超声速爆震发动机及其推进系统
CN109322761B (zh) 高马赫数飞行中发动机环形燃烧室及螺旋斜爆震燃烧方法
CN108869095B (zh) 一种超声速爆震稳定自持的边界抽吸控制方法
AU2018279792B2 (en) Flight vehicle air breathing propulsion system with isolator having obstruction
US11473500B2 (en) Method of reducing low energy flow in an isolator of a flight vehicle air breathing engine
CN104500272A (zh) 一种低流阻近壁小空间环形激波聚焦直接起爆装置
CN106930864A (zh) 一种超声速爆震发动机及其推进系统
CN110516410A (zh) 一种飞行器总体设计参数的确定方法及装置
CN115828412A (zh) 考虑外弹道约束的单室双推固体火箭发动机配比参数设计方法
CN108412618A (zh) 一种高超/超声速轴对称进气道唇口及其设计方法
CN112432204B (zh) 一种能够内流减阻的凹腔结构及超燃冲压发动机
CN109209681A (zh) 一种回流式的脉冲爆震发动机爆震管结构
CN113153569B (zh) 一种平稳排气的多管脉冲爆震发动机
CN112577070A (zh) 一种低阻高效的超燃冲压发动机推力室一体化设计方法
Strahle Theoretical consideration of combustion effects on base pressure in supersonic flight
KR101020596B1 (ko) 연소기 및 이를 포함하는 엔진
CN114673606B (zh) 采用塞式喷管的斜爆震发动机
JP2022514200A (ja) スパイラル状冷却チャネル付き燃焼室ライナー
CN109165411A (zh) 采用偏置且斜切结构喷管的固体发动机内弹道计算方法
Weidner A Pre-Mixed Shock-Induced-Combustion Approach to Inlet and Combustor Design for Hypersonic Applications

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant