CN111173644B - 一种冲压发动机性能确定方法 - Google Patents

一种冲压发动机性能确定方法 Download PDF

Info

Publication number
CN111173644B
CN111173644B CN202010009715.9A CN202010009715A CN111173644B CN 111173644 B CN111173644 B CN 111173644B CN 202010009715 A CN202010009715 A CN 202010009715A CN 111173644 B CN111173644 B CN 111173644B
Authority
CN
China
Prior art keywords
section
outlet
airflow
flow
shock
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
CN202010009715.9A
Other languages
English (en)
Other versions
CN111173644A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010009715.9A priority Critical patent/CN111173644B/zh
Publication of CN111173644A publication Critical patent/CN111173644A/zh
Application granted granted Critical
Publication of CN111173644B publication Critical patent/CN111173644B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02KJET-PROPULSION PLANTS
    • F02K7/00Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof
    • F02K7/10Plants in which the working fluid is used in a jet only, i.e. the plants not having a turbine or other engine driving a compressor or a ducted fan; Control thereof characterised by having ram-action compression, i.e. aero-thermo-dynamic-ducts or ram-jet engines

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

本发明公开一种冲压发动机性能确定方法,包括如下步骤:获取来流截面的气流参数;根据来流截面的气流参数得到发动机上各激波对应的激波截面的气流参数;根据激波截面的气流参数得到隔离段出口截面的气流参数;根据隔离段出口截面的气流参数通过热力循环分析结合化学平衡的方法得到燃烧室出口截面的气流参数;根据燃烧室出口截面的气流参数得到尾喷管出口截面的气流参数;根据尾喷管出口截面的气流参数与来流截面的气流参数得到冲压发动机的性能参数。在确定燃烧室出口截面气流参数的过程中,采用等压燃烧假设结合化学平衡的方法,提高了最终得到的发动机理论性能的准确性。

Description

一种冲压发动机性能确定方法
技术领域
本发明涉及冲压发动机技术领域,具体是一种冲压发动机性能确定方法。
背景技术
冲压发动机是超声速和高超声速飞行器的动力装置,其利用发动机的前向运动来压缩空气,根据进入燃烧室的气流为超声速还是亚声速将其分为亚燃冲压发动机和超燃冲压发动机。现有的亚燃冲压发动机的飞行马赫数一般在3左右,超燃冲压发动机目前大都在研制中,其飞行马赫数一般在6以上。飞行马赫数的不同决定了该两种类型冲压发动机进气道构型的差异,同时由于燃料燃烧平衡温度的限制,进而影响到燃烧室的燃烧组织形式。
目前亚燃及超燃冲压发动机理论性能的确定方法是主要通过一些模型简化和基本假设,计算几个典型截面的参数,进而得到飞行器整体的性能,例如在处理进气道时采用激波系假设,处理燃烧室时采用等压或等面积假设,处理尾喷管时采用等熵膨胀假设。其中,现有技术中一般将燃料的热值设为定值,该热值一般都是在常温常压下测定的,但是通常燃烧室中的温度和压力并不是常温常压,因此该方法存在一定的问题,尤其是当飞行器飞行马赫数较高,燃烧室入口静温超过2000K时,恒定热值会导致燃烧室温度超过平衡温度,使得计算结果没有意义,在工程技术上达不到计算得到的参数值,进而影响最终得到的发动机理论性能的准确性,不能为发动机前期的性能设计及实验后的性能评估提供有效的指导。
发明内容
针对上述现有技术中的不足,本发明提供一种冲压发动机性能确定方法。
为实现上述目的,本发明提供一种冲压发动机性能确定方法,包括如下步骤:
步骤1,获取来流截面的气流参数;
步骤2,根据来流截面的气流参数得到发动机上各激波对应的激波截面的气流参数;
步骤3,根据激波截面的气流参数得到隔离段出口截面的气流参数;
步骤4,根据隔离段出口截面的气流参数通过采用等压燃烧假设结合化学平衡的方法得到燃烧室出口截面的气流参数;
步骤5,根据燃烧室出口截面的气流参数得到尾喷管出口截面的气流参数;
步骤6,根据尾喷管出口截面的气流参数与来流截面的气流参数得到冲压发动机的性能参数。
作为上述技术方案的进一步改进,所述来流截面的气流参数、激波截面的气流参数、隔离段出口截面的气流参数、燃烧室出口截面的气流参数、尾喷管出口截面的气流参数均包括气流的静压、静温、总温、总压、流速、密度中的至少一种。
作为上述技术方案的进一步改进,步骤2中,发动机上各激波对应的激波截面的气流参数为:
Figure BDA0002356680830000021
Figure BDA0002356680830000022
Figure BDA0002356680830000023
Figure BDA0002356680830000024
Figure BDA0002356680830000025
Figure BDA0002356680830000026
Figure BDA0002356680830000027
Figure BDA0002356680830000028
Figure BDA0002356680830000029
Figure BDA0002356680830000031
Figure BDA0002356680830000032
Figure BDA0002356680830000033
式中,a1.1表示第一个激波截面的进气道偏折角,s1.1表示第一个激波截面的斜激波角,γ为比热比,R为空气气体常数;M0为来流截面的气流马赫数,P0为来流截面的气流静压,T0为来流截面的气流静温;M1.1为第一个激波截面的气流马赫数,P1.1为第一个激波截面的气流静压,T1.1为第一个激波截面的气流静温,Tt1.1为第一个激波截面的气流总温,ρ1.1为第一个激波截面的气流密度;M1.i-1为第i-1个激波截面的气流马赫数,P1.i-1为第i-1个激波截面的气流静压,T1.i-1为第i-1个激波截面的气流静温;a1.i表示第i个激波截面的进气道偏折角,s1.i表示第i个激波截面的斜激波角,M1.i为第i个激波截面的气流马赫数,P1.i为第i个激波截面的气流静压,T1.i为第i个激波截面的气流静温,Tt1.i为第i个激波截面的气流总温,ρ1.i为第i个激波截面的气流密度;其中,i=2,3,```,n;n表示激波截面的总数或发动机上的激波数。
作为上述技术方案的进一步改进,步骤3中,若冲压发动机为亚燃冲压发动机,发动机上隔离段出口截面的气流参数为:
Figure BDA0002356680830000034
Figure BDA0002356680830000035
Figure BDA0002356680830000036
式中,M2为隔离段出口截面的气流马赫数,P2为隔离段出口截面的气流静压,T2为隔离段出口截面的气流静温;M1.n为第n个激波截面的气流马赫数,P1.n为第n个激波截面的气流静压,T1.n为第n个激波截面的气流静温。
作为上述技术方案的进一步改进,步骤3中,若冲压发动机为超燃冲压发动机,发动机上隔离段出口截面的气流参数为:
Figure BDA0002356680830000037
Figure BDA0002356680830000038
Figure BDA0002356680830000041
Tt2=Tt1,n
Figure BDA0002356680830000042
Reθ=U1,nθ/v
Rex=U1,nx/v
Figure BDA0002356680830000043
式中,Tt2为隔离段出口截面的气流总温,Tt1,n为第n个激波截面的气流总温,D表示特征直径,θ表示动量边界层厚度,x表示隔离段长度,Reθ表示动量边界层雷诺数,Rex表示气流雷诺数,v表示气流运动粘度,μ为气流动力粘度,μ0为273.11K时空气的动力粘度,ρ1.n为第n个激波截面的气流密度。
作为上述技术方案的进一步改进,步骤4中,所述采用等压燃烧假设结合化学平衡的方法得到燃烧室出口截面的气流参数,具体为:
基于等压燃烧假设得到发动机上燃烧室出口截面的气流静压与气流流速:
P3=P2
U3=U2
Figure BDA0002356680830000044
式中,U2为隔离段出口截面的气流流速;U3为燃烧室出口截面的气流流速,P3为燃烧室出口截面的气流静压;
基于燃烧室中的燃烧化学反应的元素守恒与热力学第一、第二定律得到发动机上燃烧室出口截面的气流静温T3及燃烧产物的化学组成;
基于燃烧室出口截面的气流静压、气流流速、气流静温得到燃烧室出口截面的气流总温、气流总温总压与气流总温马赫数:
Figure BDA0002356680830000045
Figure BDA0002356680830000046
Figure BDA0002356680830000047
式中,Tt3为燃烧室出口截面的气流总温,Pt3为燃烧室出口截面的气流总压,M3为燃烧室出口截面的气流马赫数。
作为上述技术方案的进一步改进,步骤5中,通过理想膨胀假设方法,得到发动机尾喷管出口截面的气流参数为:
P4=P0
Pt4=Pt3
Tt4=Tt3
Figure BDA0002356680830000051
Figure BDA0002356680830000052
Figure BDA0002356680830000053
式中,P4为尾喷管出口截面的气流静压,Tt4为尾喷管出口截面的气流总温,Pt4为尾喷管出口截面的气流总压,M4为尾喷管出口截面的气流马赫数.U4为尾喷管出口截面的气流流速。
作为上述技术方案的进一步改进,步骤6中,所述冲压发动机的性能参数为比推力与比冲,具体为:
Figure BDA0002356680830000054
Isp=[(1+f)U4-U0]|/(fg)
式中,Fsp为比推力,Isp为比冲,f为燃空比,g为重力加速度,F为推力,m0为来流空气质量流量。
本发明提供的一种冲压发动机性能确定方法,通过获得的来流截面的气流参数进而依次确定发动机上激波截面、隔离段出口截面、燃烧室出口截面、尾喷管出口截面,最终得到飞行器整体的性能,其中,在确定燃烧室出口截面气流参数的过程中,采用等压燃烧假设结合化学平衡的方法,进而得到燃烧室出口截面的平衡温度,提高了最终得到的发动机理论性能的准确性,为发动机前期的性能设计和实验后的性能评估提供了有效的方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
图1为本发明实施例中亚燃冲压发动机构型示意图;
图2为本发明实施例中超燃冲压发动机构型示意图;
图3为本发明实施例中冲压发动机性能确定方法的流程示意图;
图4-7为采用本实施中冲压发动机性能确定方法得到的发动机在不同因素影响下的性能参数变化示意图;
图8为燃烧室产物随发动机工作当量比的分布示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明,本发明实施例中所有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
另外,在本发明中如涉及“第一”、“第二”等的描述仅用于描述目的,而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“连接”、“固定”等应做广义理解,例如,“固定”可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接,还可以是物理连接或无线通信连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
另外,本发明各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
本实施例公开了一种冲压发动机性能确定方法,利用该方法可以得到燃烧产物及相关气动参数,从而实现对飞行器性能的评估。参考图1-2,为本实施例中给出的发动机构型,其中,图1为一种亚燃冲压发动机构型示意图,图2为一种超燃冲压发动机构型示意图。图1与图2中均划分出多个截面,其中,截面1为来流截面、截面1.1为第一个激波截面、截面1.2为第二个激波截面、截面2为隔离段出口截面、截面3为燃烧室出口截面、截面4为尾喷管出口截面。本实施例中冲压发动机中进气道的构型只具有两道斜激波,因此划分出两个激波截面,即n=2;当进气道的构型发生变化导致斜激波的数量增加或减少时,激波截面的数量也需要同步增加或减少。、
本实施例中,来流截面的气流参数、激波截面的气流参数、隔离段出口截面的气流参数、燃烧室出口截面的气流参数、尾喷管出口截面的气流参数均包括气流的静压、静温、总温、总压、流速、密度中的至少一种。
参考图3,本实施例中冲压发动机性能确定方法具体包括如下步骤:
步骤1,获取来流截面的气流参数,具体包括来流的气流静压、气流静温与气流流速,其中,来流的气流流速即来流的马赫数。
步骤2,基于本实施例中冲压发动机进气道的构型,建立激波系结构,根据来流截面的气流参数采用斜激波公式得到第一激波截面与第二激波截面的气流参数,具体为:
发动机上各激波对应的激波截面的气流参数为:
Figure BDA0002356680830000071
Figure BDA0002356680830000072
Figure BDA0002356680830000073
Figure BDA0002356680830000081
Figure BDA0002356680830000082
Figure BDA0002356680830000083
Figure BDA0002356680830000084
Figure BDA0002356680830000085
Figure BDA0002356680830000086
Figure BDA0002356680830000087
Figure BDA0002356680830000088
Figure BDA0002356680830000089
Figure BDA00023566808300000810
Figure BDA00023566808300000811
式中,a1.1表示第一个激波截面的进气道偏折角,s1.1表示第一个激波截面的斜激波角,γ为比热比;M0为来流截面的气流马赫数,P0为来流截面的气流静压,T0为来流截面的气流静温;M1.1为第一个激波截面的气流马赫数,P1.1为第一个激波截面的气流静压,T1.1为第一个激波截面的气流静温,Pt1.1为第一个激波截面的气流总压,Tt1.1为第一个激波截面的气流总温,ρ1.1为第一个激波截面的气流密度;a1.2表示第二个激波截面的进气道偏折角,s1.2表示第二个激波截面的斜激波角,M1.2为第二个激波截面的气流马赫数,P1.2为第二个激波截面的气流静压,T1.2为第二个激波截面的气流静温,Pt1.2为第二个激波截面的气流总压,Tt1.2为第二个激波截面的气流总温,ρ1.2为第二个激波截面的气流密度,R为空气气体常数,取值为R=287J/kg·K;
经过式(1)-式(14)即能得到第二个激波截面的气流马赫数M1.2、气流静压P1.2、气流静温T1.2、气流密度ρ1.2与气流总温Tt1.2
步骤3,在得到第二个激波截面的气流马赫数、气流静压、气流静温与气流密度后,即能进一步得到隔离段出口截面的气流参数:
若冲压发动机为亚燃冲压发动机,即图1中所示的情况,通过采用正激波公式得到发动机上隔离段出口截面的气流参数,具体为:
Figure BDA0002356680830000091
Figure BDA0002356680830000092
Figure BDA0002356680830000093
式中,M2为隔离段出口截面的气流马赫数,P2为隔离段出口截面的气流静压,T2为隔离段出口截面的气流静温。
若冲压发动机为超燃冲压发动机,即图2中所示的情况,通过采用隔离段经验公式得到发动机上隔离段出口截面的气流马赫数,采通过采用等熵关系式和理想气体方程得到隔离段出口截面的气流静压与气流静温,具体为:
Figure BDA0002356680830000094
Figure BDA0002356680830000095
Figure BDA0002356680830000096
Tt2=Tt1,2 (21)
Figure BDA0002356680830000097
Reθ=U1,nθ/v (23)
Rex=U1,nx/v (24)
Figure BDA0002356680830000098
式中,Tt2为隔离段出口截面的气流总温,D表示特征直径,θ表示动量边界层厚度,x表示隔离段长度,Reθ表示动量边界层雷诺数,Rex表示气流雷诺数,v表示气流运动粘度,μ为气流动力粘度,μ0为273.11K时空气的动力粘度,取值为μ0=1.7894×10-5N·s/m2
经过式(18)-式(25)即能得到隔离段出口截面的气流马赫数M2、气流静压P2与气流静温T2
步骤4,根据隔离段出口截面的气流参数通过采用等压燃烧假设结合化学平衡的方法得到燃烧室出口截面的气流参数:
首先,无论是亚燃冲压发动机还是超燃冲压发动机均为等压燃烧模型,即隔离段出口截面和燃烧室出口截面处的静压相等,根据动量守恒也能得到隔离段出口截面和燃烧室出口截面处气流流速也相等,即满足式(26)-(28):
P3=P2 (26)
U3=U2 (27)
Figure BDA0002356680830000101
式中,U2为隔离段出口截面的气流流速,c为声速;U3为燃烧室出口截面的气流流速,P3为燃烧室出口截面的气流静;
燃烧室中的燃烧化学反应的元素守恒与热力学第一、第二定律得到发动机上燃烧室出口截面的气流静温T3
例如,发动机的燃烧物为化合物CxHyOzNuSv,其中x、y、z、u、v分别表示C、H、O、N、S的含量,发动机的燃烧产物含有CO、CO2、H2O、N2、O2、H2、SO2
燃烧室中的燃烧化学反应为:
CxHyOzNuSv+m(O2+3.76N2)→n1CO+n2CO2+n3H2O+n4N2+n5O2+n6H2+n7SO2
其中,(O2+3.76N2)为助燃剂,m为可根据当量比确定的助燃剂含量,n1、n2、n3、n4、n5、n6、n7为燃烧产物的化学计量系数;
通过元素守恒与CO2、H2O的逆反应以及热力学第一定律和第二定律得到n1、n2、n3、n4、n5、n6、n7:和燃烧室出口截面的气流静温T3
n1+n2=x (29)
n1+2n2+n3+2n5+2n7=z+2m (30)
2n3+2n6=y (31)
2n4=3.76m (32)
n7=v (33)
Figure BDA0002356680830000102
Figure BDA0002356680830000103
Figure BDA0002356680830000111
Figure BDA0002356680830000112
式中,
Figure BDA0002356680830000113
Figure BDA0002356680830000114
为逆反应的化学平衡常数,可通过吉布斯自由能ΔG确定,
Figure BDA0002356680830000115
Figure BDA0002356680830000116
可以根据物质热力学库查询,Pstandard为热力学中ΔG对应的测量压力,通常为1atm,Rconst为通用气体常数,取Rconst=8.314J/(mol·K)。
经过对式(29)-式(37)的求解即能得到燃烧室出口截面的气流静温T3。其中,求解的具有过程为首先将隔离段截面出口的气流静温T2给定为T3的初值,即T3=T2;随后逐步增加T3的取值,例如每次增加0.1K,此时每个T3值均会存在相应含量的燃烧产物,直至求出的n1、n2、n3、n4、n5、n6、n7既满足热力学第一定律也满足热力学第二定律,此时求出的T3即为平衡温度,即是燃烧室出口截面的气流静温。
再例如,发动机的燃烧物为氢气H2,,发动机的燃烧产物含有H2O、N2、O2、H2、H、O;
燃烧室中的燃烧化学反应为:
H2+m(O2+3.76N2)→n1H2O+n2N2+n3O2+n4H2+n5H+n6O+n7OH+n8NO
其中,(O2+3.76N2)为助燃剂,m为可根据当量比确定的助燃剂含量,n1、n2、n3、n4、n5、n6、n7为燃烧产物的化学计量系数;
通过元素守恒与H2O、OH、NO、O2、H2的逆反应以及热力学第一定律和第二定律得到n1、n2、n3、n4、n5、n6、n7、n8:和燃烧室出口截面的气流静温T3
n1+n4+0.5n5+0.5n7=1 (38)
0.5n1+n3+0.5n6+0.5n7+0.5n8=m (39)
n2+0.5n8=3.76m (40)
Figure BDA0002356680830000117
Figure BDA0002356680830000118
Figure BDA0002356680830000119
Figure BDA00023566808300001110
Figure BDA00023566808300001111
式中,
Figure BDA0002356680830000121
KOH、KNO、KO、KH为逆反应的化学平衡常数,可通过吉布斯自由能ΔG确定。
经过对式(38)-式(45)的求解即能得到燃烧室出口截面的气流静温T3
在得到燃烧室出口截面的气流静温T3后,即能通过常规流体公式得到燃烧室出口截面气流的总温、总压与马赫数,具体为:
Figure BDA0002356680830000122
Figure BDA0002356680830000123
Figure BDA0002356680830000124
式中,Tt3为燃烧室出口截面的气流总温,Pt3为燃烧室出口截面的气流总压,M3为燃烧室出口截面的气流马赫数;
经过式(26)-式(48)即能得到燃烧室出口截面的气流静温T3、气流静压P3、气流总温Tt3、气流总压Pt3
步骤5,采用理想膨胀假设,根据燃烧室出口截面的气流参数得到尾喷管出口截面的气流参数,具体为:
P4=P0 (49)
Pt4=Pt3 (50)
Tt4=Tt3 (51)
Figure BDA0002356680830000125
Figure BDA0002356680830000126
Figure BDA0002356680830000127
式中,P4为尾喷管出口截面的气流静压,Tt4为尾喷管出口截面的气流总温,Pt4为尾喷管出口截面的气流总压,M4为尾喷管出口截面的气流马赫数.U4为尾喷管出口截面的气流流速。
经过式(49)-式(54)即能得到尾喷管出口截面的气流流速U4
步骤6,根据尾喷管出口截面的气流参数与来流截面的气流参数得到冲压发动机的性能参数,其中,冲压发动机的性能参数为比推力与比冲,具体为:
Figure BDA0002356680830000128
Isp=[(1+f)U4-U0]/(f·g)
式中,Fsp为比推力,Isp为比冲,f为燃空比,g为重力加速度,F为推力,m0为来流空气质量流量。
优选的,在步骤4中采用了理想膨胀假设,在实施过程中亦可采用给定喷管膨胀比的方式进行计算,此时需根据下求得燃烧室出口截面的马赫数M3,再结合式(51)求得燃烧室出口截面处的其他气流参数。
为了验证本实施例中冲压发动机性能确定方法的可靠性,采用了公开文献中HyShotII的飞行实验和地面实验的数据,将其与本实施例中冲压发动机性能确定方法的计算结果对比。本实施例中冲压发动机性能确定方法所计算的发动机在不同因素影响下的性能参数变化如图4-7所示,对比结果如表1所示。从表1可知,采用本实施例中方法得到的计算结果与公开文献中的实验数据误差较小,基本都保持在5%以内,充分说明了本实施例中冲压发动机性能确定方法的准确性。
参考图4-7为采用本实施中冲压发动机性能确定方法得到的发动机在不同因素影响下的性能参数变化示意图,图8为燃烧室产物随发动机工作当量比的分布示意图,通过图8有助于分析不同因素影响下的燃烧室燃烧状况,进而可帮助理解发动机性能参数为什么这样变化。
表1数据对比结果
Figure BDA0002356680830000131
以上所述仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是在本发明的发明构思下,利用本发明说明书及附图内容所作的等效结构变换,或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。

Claims (7)

1.一种冲压发动机性能确定方法,其特征在于,包括如下步骤:
步骤1,获取来流截面的气流参数;
步骤2,根据来流截面的气流参数得到发动机上各激波对应的激波截面的气流参数;
步骤3,根据激波截面的气流参数得到隔离段出口截面的气流参数;
步骤4,根据隔离段出口截面的气流参数通过采用等压燃烧假设结合化学平衡的方法得到燃烧室出口截面的气流参数,具体为:
基于等压燃烧假设得到发动机上燃烧室出口截面的气流静压与气流流速:
P3=P2
U3=U2
Figure FDA0003031199980000011
式中,U2为隔离段出口截面的气流流速;U3为燃烧室出口截面的气流流速,P3为燃烧室出口截面的气流静压,P2为隔离段出口截面的气流静压,γ为比热比,R为空气气体常数,T2为隔离段出口截面的气流静温,M2为隔离段出口截面的气流马赫数;
基于燃烧室中的燃烧化学反应的元素守恒与热力学第一、第二定律得到发动机上燃烧室出口截面的气流静温T3及燃烧产物的化学组成,具体为:
基于燃烧室中燃烧物的化学反应的燃烧物与燃烧产物之间的逆反应以及热力学第一定律和第二定律建立方程组,对方程组求解即得到发动机上燃烧室出口截面的气流静温T3,其求解过程为:首先将隔离段出口截面的气流静温T2给定为T3的初值,即T3=T2;随后逐步增加T3的取值,此时每个T3值均会存在相应含量的燃烧产物,直至求出的n1~X既满足热力学第一定律也满足热力学第二定律,此时的T3即为平衡温度,即是燃烧室出口截面的气流静温,其中,n1~X表示各燃烧产物的化学计量系数;
基于燃烧室出口截面的气流静压、气流流速、气流静温得到燃烧室出口截面的气流总温、气流总压与气流马赫数:
Figure FDA0003031199980000012
Figure FDA0003031199980000013
Figure FDA0003031199980000021
式中,Tt3为燃烧室出口截面的气流总温,Pt3为燃烧室出口截面的气流总压,M3为燃烧室出口截面的气流马赫数;
步骤5,根据燃烧室出口截面的气流参数得到尾喷管出口截面的气流参数;
步骤6,根据尾喷管出口截面的气流参数与来流截面的气流参数得到冲压发动机的性能参数。
2.根据权利要求1所述冲压发动机性能确定方法,其特征在于,所述来流截面的气流参数、激波截面的气流参数、隔离段出口截面的气流参数、燃烧室出口截面的气流参数、尾喷管出口截面的气流参数均包括气流的静压、静温、总温、总压、流速、密度中的至少一种。
3.根据权利要求2所述冲压发动机性能确定方法,其特征在于,步骤2中,发动机上各激波对应的激波截面的气流参数为:
Figure FDA0003031199980000022
Figure FDA0003031199980000023
Figure FDA0003031199980000024
Figure FDA0003031199980000025
Figure FDA0003031199980000026
Figure FDA0003031199980000027
Figure FDA0003031199980000028
Figure FDA0003031199980000029
Figure FDA00030311999800000210
Figure FDA00030311999800000211
Figure FDA00030311999800000212
Figure FDA0003031199980000031
式中,a1.1表示第一个激波截面的进气道偏折角,s1.1表示第一个激波截面的斜激波角,γ为比热比,R为空气气体常数;M0为来流截面的气流马赫数,P0为来流截面的气流静压,T0为来流截面的气流静温;M1.1为第一个激波截面的气流马赫数,P1.1为第一个激波截面的气流静压,T1.1为第一个激波截面的气流静温,Tt1.1为第一个激波截面的气流总温,ρ1.1为第一个激波截面的气流密度;M1.i-1为第i-1个激波截面的气流马赫数,P1.i-1为第i-1个激波截面的气流静压,T1.i-1为第i-1个激波截面的气流静温;a1.i表示第i个激波截面的进气道偏折角,s1.i表示第i个激波截面的斜激波角,M1.i为第i个激波截面的气流马赫数,P1.i为第i个激波截面的气流静压,T1.i为第i个激波截面的气流静温,Tt1.i为第i个激波截面的气流总温,ρ1.i为第i个激波截面的气流密度;其中,i=2,3,……,n;n表示激波截面的总数或发动机上的激波数。
4.根据权利要求3所述冲压发动机性能确定方法,其特征在于,步骤3中,若冲压发动机为亚燃冲压发动机,发动机上隔离段出口截面的气流参数为:
Figure FDA0003031199980000032
Figure FDA0003031199980000033
Figure FDA0003031199980000034
式中,M2为隔离段出口截面的气流马赫数,P2为隔离段出口截面的气流静压,T2为隔离段出口截面的气流静温;M1.n为第n个激波截面的气流马赫数,P1.n为第n个激波截面的气流静压,T1.n为第n个激波截面的气流静温。
5.根据权利要求3所述冲压发动机性能确定方法,其特征在于,步骤3中,若冲压发动机为超燃冲压发动机,发动机上隔离段出口截面的气流参数为:
Figure FDA0003031199980000035
Figure FDA0003031199980000036
Figure FDA0003031199980000041
Tt2=Tt1,n
Figure FDA0003031199980000042
Reθ=U1,nθ/v
Rex=U1,nx/v
Figure FDA0003031199980000043
式中,Tt2为隔离段出口截面的气流总温,Tt1,n为第n个激波截面的气流总温,D表示特征直径,θ表示动量边界层厚度,x表示隔离段长度,Reθ表示动量边界层雷诺数,Rex表示气流雷诺数,v表示气流运动粘度,μ为气流动力粘度,μ0为273.11K时空气的动力粘度,ρ1.n为第n个激波截面的气流密度。
6.根据权利要求1所述冲压发动机性能确定方法,其特征在于,步骤5中,通过理想膨胀假设方法,得到发动机尾喷管出口截面的气流参数为:
P4=P0
Pt4=Pt3
Tt4=Tt3
Figure FDA0003031199980000044
Figure FDA0003031199980000045
Figure FDA0003031199980000046
式中,P4为尾喷管出口截面的气流静压,Tt4为尾喷管出口截面的气流总温,Pt4为尾喷管出口截面的气流总压,M4为尾喷管出口截面的气流马赫数,U4为尾喷管出口截面的气流流速,P0为来流截面的气流静压。
7.根据权利要求6所述冲压发动机性能确定方法,其特征在于,步骤6中,所述冲压发动机的性能参数为比推力与比冲,具体为:
Figure FDA0003031199980000047
Isp=[(1+f)U4-U0]/(fg)
式中,Fsp为比推力,Isp为比冲,f为燃空比,g为重力加速度,F为推力,m0为来流空气质量流量。
CN202010009715.9A 2020-01-06 2020-01-06 一种冲压发动机性能确定方法 Active CN111173644B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010009715.9A CN111173644B (zh) 2020-01-06 2020-01-06 一种冲压发动机性能确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010009715.9A CN111173644B (zh) 2020-01-06 2020-01-06 一种冲压发动机性能确定方法

Publications (2)

Publication Number Publication Date
CN111173644A CN111173644A (zh) 2020-05-19
CN111173644B true CN111173644B (zh) 2021-06-04

Family

ID=70652479

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010009715.9A Active CN111173644B (zh) 2020-01-06 2020-01-06 一种冲压发动机性能确定方法

Country Status (1)

Country Link
CN (1) CN111173644B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114722743B (zh) * 2022-05-24 2022-11-01 中国人民解放军国防科技大学 基于燃烧室化学平衡的超燃冲压发动机一维性能估计方法
CN115618171B (zh) * 2022-06-06 2023-10-24 北京理工大学 一种基于同伦算法的推进剂燃烧平衡产物求解方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103870683A (zh) * 2014-03-03 2014-06-18 北京动力机械研究所 一种超燃冲压发动机燃烧室性能预估方法
CN103970957A (zh) * 2014-05-20 2014-08-06 北京航空航天大学 一种弹性乘波体高超声速飞行器仿真方法
CN105184015A (zh) * 2015-09-29 2015-12-23 北京动力机械研究所 基于功能样机的冲压发动机验证方法及验证系统
CN106407571A (zh) * 2016-09-22 2017-02-15 北京机械设备研究所 一种高超声速吸气式冲压发动机气动推力的分析方法
CN109063407A (zh) * 2018-10-29 2018-12-21 南京航空航天大学 一种超燃冲压发动机稳态模型的建模方法
CN110516410A (zh) * 2019-10-21 2019-11-29 湖南德雅坤创科技有限公司 一种飞行器总体设计参数的确定方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103870683A (zh) * 2014-03-03 2014-06-18 北京动力机械研究所 一种超燃冲压发动机燃烧室性能预估方法
CN103970957A (zh) * 2014-05-20 2014-08-06 北京航空航天大学 一种弹性乘波体高超声速飞行器仿真方法
CN105184015A (zh) * 2015-09-29 2015-12-23 北京动力机械研究所 基于功能样机的冲压发动机验证方法及验证系统
CN106407571A (zh) * 2016-09-22 2017-02-15 北京机械设备研究所 一种高超声速吸气式冲压发动机气动推力的分析方法
CN109063407A (zh) * 2018-10-29 2018-12-21 南京航空航天大学 一种超燃冲压发动机稳态模型的建模方法
CN110516410A (zh) * 2019-10-21 2019-11-29 湖南德雅坤创科技有限公司 一种飞行器总体设计参数的确定方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吸气式高超声速飞行器推进系统耦合建模与分析;唐硕等;《飞行力学》;20130630;第31卷(第3期);第244-249页 *
唐硕等.吸气式高超声速飞行器推进系统耦合建模与分析.《飞行力学》.2013,第31卷(第3期), *

Also Published As

Publication number Publication date
CN111173644A (zh) 2020-05-19

Similar Documents

Publication Publication Date Title
Nicholls et al. Intermittent detonation as a thrust-producing mechanism
Herrmann et al. Experimental study of boundary-layer bleed impact on ramjet inlet performance
CN111173644B (zh) 一种冲压发动机性能确定方法
Canning et al. Ballistic-range technology
CN110362783A (zh) 一种火箭基组合循环发动机推力计算方法
CN105446167B (zh) 高超声速超燃冲压发动机实时模型、仿真方法
CN111157248B (zh) 基于地面直连试验的冲压发动机及其燃烧室性能评估方法
CN115618765A (zh) 一种基于asum格式的控制用超燃冲压发动机建模方法
CN115060504A (zh) 冲压发动机燃烧模态及隔离段气流参数实时确定方法
Faget et al. Free-jet Tests of a 6.5-inch-diameter Ram-jet Engine at Mach Numbers 1.81 and 2.00
Hillyer et al. Extension of test time in Ludwieg Tunnels
CN115982943A (zh) 一种超燃冲压发动机容积动力学一维建模方法
CN112985681B (zh) 一种rbcc引射模态内推力实验测量评估方法
Sterling et al. Numerical investigations of pulse detonation wave engines
Lee et al. Starting characteristics of the hypersonic wind tunnel with the mach number variation
Masumoto et al. Experimental study on effect of boundary layer on combustion modes in a supersonic combustor
Cooper et al. Detonation tube impulse in subatmospheric environments
CN116562194B (zh) 冲压旋转爆震发动机推力评估方法及系统
Dyer et al. Reaction ratio and nozzle expansion effects on the PDE
Khashwain et al. SCRAMJET: a supersonic transport
Tsuboi et al. Numerical study and performance evaluation for pulse detonation engine with an aerospike nozzle
Ismail et al. Thrust Chamber Dynamic and Propulsive Performance of Biofuels Under Detonation Combustion
Cox Experimental facilities for hypersonic research
Vargas et al. Propagation Of A Detonation Wave In Supersonic Crossflows Using A Pre-Detonator
Nekkanti Analysis of thrust development in a pulse detonaton engine

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