CN116562193B - 旋转爆震发动机燃烧效率分析方法及系统 - Google Patents
旋转爆震发动机燃烧效率分析方法及系统 Download PDFInfo
- Publication number
- CN116562193B CN116562193B CN202310835904.5A CN202310835904A CN116562193B CN 116562193 B CN116562193 B CN 116562193B CN 202310835904 A CN202310835904 A CN 202310835904A CN 116562193 B CN116562193 B CN 116562193B
- Authority
- CN
- China
- Prior art keywords
- outlet
- combustion chamber
- specific heat
- ave
- combustion
- 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
Links
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 401
- 238000005474 detonation Methods 0.000 title claims abstract description 36
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 74
- 230000003068 static effect Effects 0.000 claims abstract description 48
- 238000012360 testing method Methods 0.000 claims abstract description 19
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 16
- 239000007789 gas Substances 0.000 claims description 114
- 239000007921 spray Substances 0.000 claims description 64
- 239000000446 fuel Substances 0.000 claims description 30
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 18
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 9
- 239000001569 carbon dioxide Substances 0.000 claims description 9
- 239000000295 fuel oil Substances 0.000 claims description 8
- 229910052757 nitrogen Inorganic materials 0.000 claims description 8
- 239000000203 mixture Substances 0.000 claims description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 239000001301 oxygen Substances 0.000 claims description 6
- 229910052760 oxygen Inorganic materials 0.000 claims description 6
- 239000003921 oil Substances 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- UBAZGMLMVVQSCD-UHFFFAOYSA-N carbon dioxide;molecular oxygen Chemical compound O=O.O=C=O UBAZGMLMVVQSCD-UHFFFAOYSA-N 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 239000003350 kerosene Substances 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 2
- 230000002452 interceptive effect Effects 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 6
- 238000001816 cooling Methods 0.000 abstract description 4
- 230000004044 response Effects 0.000 abstract description 4
- 239000002737 fuel gas Substances 0.000 description 13
- 230000008859 change Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000009530 blood pressure measurement Methods 0.000 description 2
- 230000000977 initiatory effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000009529 body temperature measurement Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 230000001590 oxidative effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/04—Testing internal-combustion engines
- G01M15/05—Testing internal-combustion engines by combined monitoring of two or more different engine parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明提供了一种旋转爆震发动机燃烧效率分析方法及系统,涉及旋转爆震发动机技术领域,采用沿燃烧室出口周向均匀间隔分布的多个测试位置进行压力检测,并由此计算燃烧室出口静压,相较于温度检测响应速度更快,并且可以解决旋转爆震发动机因内部流场时间非定常性和空间非均匀性导致的燃烧效率难以准确分析的技术问题,分析燃烧效率仅通过压力检测,根据气体动力学函数计算得出燃烧室出口总温,并根据温升法计算燃烧室燃烧效率,无需布置热电偶,进而无需设置水冷装置,燃烧效率分析更加便捷、准确。
Description
技术领域
本发明涉及旋转爆震发动机技术领域,尤其是涉及一种旋转爆震发动机燃烧效率分析方法及系统。
背景技术
旋转爆震发动机运行时,燃料与氧化剂的混合气在燃烧室内经旋转爆震波燃烧后,燃料中的化学能转变为热能加热工质,形成高压高温燃气,燃气经喷管流出发动机形成推力。燃料燃烧时,绝大部分燃料释放热量,但有一部燃料来不及燃烧直接随气流流出喷管。将实际加热工质的热量与该燃料在绝热条件下完全燃烧时所释出的热量之比称为燃烧效率,通过燃烧效率可评价旋转爆震发动机的燃料有效利用率,并以此作为旋转爆震发动机的重要性能参数。
现有的燃烧效率分析方法,主要是通过测量温度和燃气成分分析,但旋转爆震发动机流场复杂,流场存在空间上非均匀性,具有周向和轴向速度分量,难以准确测温,进而导致原有燃烧效率分析方法难以适用。此外,温度传感器响应速度慢,在时间上存在非定常性,不能反映瞬态的燃烧效率。
发明内容
本发明的目的在于提供一种旋转爆震发动机燃烧效率分析方法及系统,以缓解旋转爆震发动机燃烧效率分析结果存在时间和空间非定常性的技术问题。
第一方面,本发明提供的旋转爆震发动机燃烧效率分析方法包括以下步骤:
获取沿燃烧室出口周向均匀间隔分布的多个测试位置在爆震波稳定运动时的压力pi,t,其中,i为测试位置序号,t为压力采集时间;
根据公式,计算燃烧室出口静压ps4,其中,x为压力测试位置总数;
根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4;
根据温升法计算燃烧室燃烧效率算值ηc。
结合第一方面,本发明提供了第一方面的第一种可能的实施方式,其中,根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4的步骤包括:
计算燃烧室内气体质量平均气体常数Rave;
获取燃烧室出口比热比初值γ4、喷管出口比热比初值γ8和喷管出口马赫数Ma8;
获取燃烧室出口面积A4和排气喷管出口面积A8;
根据公式,计算燃烧室出口马赫数Ma4,其中,Ma8=1;
根据公式,计算燃烧室出口总压Pt4;
根据公式,计算燃烧室出口总温Tt4;
根据公式,计算燃烧室出口静温Ts4。
结合第一方面的第一种可能的实施方式,本发明提供了第一方面的第二种可能的实施方式,其中,计算燃烧室内气体质量平均气体常数Rave的步骤包括:
获取燃烧效率初值η0;
根据公式,计算参与燃烧的燃油流量mf,c;
根据公式,计算参与燃烧的空气流量ma,c,L0为完全燃烧单位燃油需要的理论空气量;
根据公式,计算未参与燃烧的燃油流量mf,nc;
根据公式,计算未参与燃烧的空气流量ma,nc;
设定燃油平均化学式为CmHn,且燃烧产物为二氧化碳和水蒸气;
根据公式,计算燃烧产物中二氧化碳的流量/>,其中,为二氧化碳分子质量,/>为燃油分子质量;
根据公式,计算燃烧产物中的水蒸气流量/>,其中,为水蒸气分子质量;
根据公式,计算未参与燃烧空气中的氮气流量/>,其中,/>为氮气的质量分数;
根据公式,计算未参与燃烧空气中的氧气流量/>,其中,为氧气的质量分数;
根据公式,计算燃烧室内各组分气体的质量占比/>,其中,j组分包含氮气、氧气、二氧化碳、水蒸气和煤油蒸汽,n0=5;
根据公式,计算燃烧室内气体质量平均气体常数Rave,其中,Rj为j组分的气体常数。
结合第一方面,本发明提供了第一方面的第三种可能的实施方式,其中,根据温升法计算燃烧室燃烧效率算值ηc的步骤包括:
计算燃烧室出口气体质量平均定压比热C4p,ave;
根据公式,计算燃烧室出口比热比算值γ4,ave,其中,C4p,ave为燃烧室出口气体质量平均定压比热,C4v,ave为燃烧室出口气体质量平均定容比热,且C4v,ave=C4p,ave-Rave;
计算喷管出口总压Pt8、喷管出口总温Tt8和喷管出口静温Ts8;
计算喷管出口气体质量平均定压比热C8p,ave;
根据公式,计算喷管出口比热比算值γ8ave,其中,C8p,ave为喷管出口气体质量平均定压比热,C8v,ave为喷管出口气体质量平均定容比热,且C8v,ave=C8p,ave-Rave;
根据公式,计算燃烧室燃烧效率算值ηc,其中,/>为空气的定压比热,/>为燃油的定压比热,Tt2为喷注段入口空气总温,ma为进入燃烧室的空气流量,mf为进入燃烧室的燃油流量,mg为燃烧室内的油气混合物流量,mg=ma+mf,LHV为燃油的低热值。
结合第一方面的第三种可能的实施方式,本发明提供了第一方面的第四种可能的实施方式,其中,对燃烧室出口比热比、喷管出口比热比和燃烧室燃烧效率三者中至少其一进行迭代;
对燃烧室出口比热比进行迭代的步骤包括:比较燃烧室出口比热比初值γ4与燃烧室出口比热比算值γ4,ave的差值是否小于第一预设差值,若否则将燃烧室出口比热比算值γ4,ave赋值于燃烧室出口比热比初值γ4,重新计算燃烧室出口气体质量平均定压比热C4p,ave和燃烧室出口比热比算值γ4,ave,直至燃烧室出口比热比初值γ4与燃烧室出口比热比算值γ4ave的差值小于第一预设差值;
对喷管出口比热比进行迭代的步骤包括:比较喷管出口比热比初值γ8与喷管出口比热比算值γ8,ave的差值是否小于第二预设差值,若否则将喷管出口比热比算值γ8,ave赋值于喷管出口比热比初值γ8,重新计算喷管出口气体质量平均定压比热C8p,ave和喷管出口比热比算值γ8,ave,直至喷管出口比热比初值γ8与喷管出口比热比算值γ8,ave的差值小于第二预设差值;
对燃烧室燃烧效率进行迭代的步骤包括:比较燃烧室燃烧效率算值ηc与燃烧效率初值η0的差值是否小于第三预设差值,若否则将燃烧室燃烧效率算值ηc赋值于燃烧效率初值η0,根据燃烧效率初值η0计算燃烧室内各种物质组分,重新计算燃烧室内气体质量平均气体常数Rave、燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4、燃烧室出口静温Ts4、燃烧室出口气体质量平均定压比热C4p,ave、燃烧室出口比热比算值γ4,ave、喷管出口气体质量平均定压比热C8p,ave、喷管出口比热比算值γ8ave和燃烧室燃烧效率算值ηc,直至燃烧室燃烧效率算值ηc与燃烧效率初值η0的差值小于第三预设差值。
结合第一方面的第三种可能的实施方式,本发明提供了第一方面的第五种可能的实施方式,其中,计算燃烧室出口气体质量平均定压比热C4p,ave的步骤包括:
根据公式,计算燃烧室出口气体质量平均定压比热C4p,ave,其中,/>为燃烧室内i组分气体的质量占比,/>为燃烧室出口处i组分的定压比热,,a4、b4、d4、e4和f4分别为多项式计算参数。
结合第一方面的第三种可能的实施方式,本发明提供了第一方面的第六种可能的实施方式,其中,计算喷管出口总压Pt8、喷管出口总温Tt8和喷管出口静温Ts8的步骤包括:
根据公式,计算排气喷管出口总压Pt8,其中,Ma8=1;
根据公式,计算燃烧室出口总温Tt8,其中,A8为排气喷管出口面积;
根据公式,计算燃烧室出口静温Ts8。
结合第一方面的第三种可能的实施方式,本发明提供了第一方面的第七种可能的实施方式,其中,计算喷管出口气体质量平均定压比热C8p,ave的步骤包括:
根据公式,计算喷管出口气体质量平均定压比热C8p,ave,其中,/>为燃烧室内i组分气体的质量占比,/>为喷管出口处i组分的定压比热,,a8、b8、d8、e8和f8分别为多项式计算参数。
第二方面,本发明提供的旋转爆震发动机燃烧效率分析系统,包括:
多个压力传感器,多个压力传感器分别连接于火焰筒,且多个压力传感器沿燃烧室出口周向均匀间隔分布;
交互器,用于输入初始条件,并输出分析结果,初始条件包括:燃烧室出口比热比初值γ4、喷管出口比热比初值γ8、喷管出口马赫数Ma8、燃烧效率初值η0、获取燃烧室出口面积A4和排气喷管出口面积A8,分析结果包括燃烧室燃烧效率算值ηc;
处理器,交互器和多个压力传感器分别与处理器连接,处理器用于计算获得燃烧室出口静压ps4、燃烧室内各种物质组分、燃烧室内气体质量平均气体常数Rave、燃烧室出口气体质量平均定压比热C4p,ave、燃烧室出口气体质量平均定容比热C4v,ave、燃烧室出口比热比算值γ4,ave、喷管出口气体质量平均定压比热C8p,ave、喷管出口气体质量平均定容比热C8v,ave、喷管出口比热比算值γ8,ave和燃烧室燃烧效率算值ηc。
结合第二方面,本发明提供了第二方面的第一种可能的实施方式,其中,旋转爆震发动机燃烧效率分析系统还包括:流量计和温度传感器,流量计和温度传感器分别与处理器连接;
流量计安装于供油管路,温度传感器安装于旋转爆震发动机的喷注段。
本发明实施例带来了以下有益效果:采用沿燃烧室出口周向均匀间隔分布的多个测试位置进行压力检测,并由此计算燃烧室出口静压,相较于温度检测响应速度更快,并且可以解决旋转爆震发动机因内部流场时间非定常性和空间非均匀性导致的燃烧效率难以准确分析的技术问题,分析燃烧效率仅通过压力检测,根据气体动力学函数计算得出燃烧室出口总温,并根据温升法计算燃烧室燃烧效率,无需布置热电偶,进而无需设置水冷装置,燃烧效率分析更加便捷、准确。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或相关技术中的技术方案,下面将对具体实施方式或相关技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的旋转爆震发动机燃烧效率分析方法的流程示意图;
图2为旋转爆震发动机的剖视图;
图3为旋转爆震发动机的示意图。
图标:1-喷注段;2-扩压器;3-火焰筒;4-排气喷管;5-喷嘴;6-起爆装置;7-压力传感器。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。公式中的物理量,如无单独标注,应理解为国际单位制基本单位的基本量,或者,由基本量通过乘、除、微分或积分等数学运算导出的导出量。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
本发明实施例提供的旋转爆震发动机燃烧效率分析方法,包括:获取沿燃烧室出口周向均匀间隔分布的多个测试位置在爆震波稳定运动时的压力pi,t,其中,i为测试位置序号,t为压力采集时间;根据公式,计算燃烧室出口静压ps4,其中,x为压力测试位置总数;根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4;根据温升法计算燃烧室燃烧效率算值ηc。
采用上述旋转爆震发动机燃烧效率分析方法,即使爆震波在燃烧室内沿周向高频传播,导致燃烧室在空间上存在高度的非均匀性,在时间上存在强非定常性,通过沿燃烧室出口周向均匀间隔分布的多个测试位置分别进行压力采集,并计算燃烧室出口静压ps4,不仅可以解决静压测量的时间非定常性问题,而且相较于以往通过测温分析燃烧效率的方式压力响应速度更快,通过计算燃烧室出口静压,根据气体动力学函数,得到燃烧室出口总温Tt4和燃烧室出口静温Ts4,根据温升法计算燃烧室燃烧效率算值ηc,无需布置热电偶,进而无需设置水冷装置,燃烧效率分析更加便捷、准确。
需要说明的是,可在本方法的初始阶段获取输入燃烧室出口面积A4、排气喷管出口面积A8、空气质量流量ma、空气定压比热Cpa、燃油质量流mf、燃油定压比热Cpf和燃烧室入口总温Tt2和燃烧效率初值η0。其中,燃烧室出口面积A4和排气喷管出口面积A8可提前测算获得,燃油质量流mf通过安装于供油管路的流量计测得,压缩空气进入电加热器前采用流量函数法,测量总温总压,结合孔板吼道面积,间接测出进入发动机的空气质量流量ma,在喷注段入口通过温度传感器测量获得燃烧室入口总温Tt2,燃烧效率初值η0可根据相似型号的旋转爆震发动机预估获得,亦可由其他燃烧效率分析方法获得,并代入本实施方式提供的分析方法,从而获得更佳准确的燃烧室燃烧效率算值ηc。此外,在根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4的步骤中,通过迭代计算燃烧室内气体质量平均气体常数Rave和比热比γave,可以克服燃气的比热和比热比等物性参数,随燃气温度和燃气组分变化而产生较大变化的问题,能够实现以燃烧室出口静压和燃烧室流量为约束,迭代求解燃气的比热比等物性参数和燃气的温度等气动热力学参数,进而精准求解燃烧室燃烧效率。上述燃气不仅指燃油,进入燃烧室的空气和燃料的混合气统称为燃烧室内的燃气。
此外,在计算燃烧室出口静压ps4的步骤中,∆t大于等于10ms,每个测试位置压力pi,t的采样频率大于等于10kHz。通过求解每个压力传感器7∆t时间内的算数平均压力,并将沿燃烧室出口周向均匀间隔分布的多个压力值求平均,从而得到准确的燃烧室出口静压ps4,可以避免分析结果存在时间和空间非定常性的技术问题。
在本发明实施例中,根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4的步骤包括:计算燃烧室内气体质量平均气体常数Rave;获取燃烧室出口比热比初值γ4、喷管出口比热比初值γ8和喷管出口马赫数Ma8;获取燃烧室出口面积A4和排气喷管出口面积A8;在发动机排气喷管达到临界状态时进行测试,根据公式,计算燃烧室出口马赫数Ma4,其中,Ma8=1;根据公式/>,计算燃烧室出口总压Pt4;根据公式,计算燃烧室出口总温Tt4;根据公式/>,计算燃烧室出口静温Ts4。
进一步的,计算燃烧室内气体质量平均气体常数Rave的步骤包括:获取燃烧效率初值η0;根据公式,计算参与燃烧的燃油流量mf,c;根据公式/>,计算参与燃烧的空气流量ma,c,L0为完全燃烧单位燃油需要的理论空气量;根据公式,计算未参与燃烧的燃油流量mf,nc;根据公式/>,计算未参与燃烧的空气流量ma,nc;设定燃油平均化学式为CmHn,且燃烧产物为二氧化碳和水蒸气;根据公式/>,计算燃烧产物中二氧化碳的流量/>,其中,/>为二氧化碳分子质量,/>为燃油分子质量;根据公式/>,计算燃烧产物中的水蒸气流量/>,其中,/>为水蒸气分子质量;根据公式,计算未参与燃烧空气中的氮气流量/>,其中,/>为氮气的质量分数;根据公式/>,计算未参与燃烧空气中的氧气流量/>,其中,/>为氧气的质量分数;根据公式/>,计算燃烧室内各组分气体的质量占比/>,其中,j组分包含氮气、氧气、二氧化碳、水蒸气和煤油蒸汽,n0=5;根据公式/>,计算燃烧室内气体质量平均气体常数Rave,其中,/>为j组分的气体常数。
进一步的,根据温升法计算燃烧室燃烧效率算值ηc的步骤包括:计算燃烧室出口气体质量平均定压比热C4p,ave;根据公式,计算燃烧室出口比热比算值γ4,ave,其中,C4p,ave为燃烧室出口气体质量平均定压比热,C4v,ave为燃烧室出口气体质量平均定容比热,且C4v,ave=C4p,ave-Rave;计算喷管出口总压Pt8、喷管出口总温Tt8和喷管出口静温Ts8;计算喷管出口气体质量平均定压比热C8p,ave;根据公式/>,计算喷管出口比热比算值γ8ave,其中,C8p,ave为喷管出口气体质量平均定压比热,C8v,ave为喷管出口气体质量平均定容比热,且C8v,ave=C8p,ave-Rave;根据公式/>,计算燃烧室燃烧效率算值ηc,其中,/>为空气的定压比热,/>为燃油的定压比热,Tt2为喷注段入口空气总温,ma为进入燃烧室的空气流量,mf为进入燃烧室的燃油流量,mg为燃烧室内的油气混合物流量,mg=ma+mf,LHV为燃油的低热值。
进一步的,对燃烧室出口比热比、喷管出口比热比和燃烧室燃烧效率三者中至少其一进行迭代;对燃烧室出口比热比进行迭代的步骤包括:比较燃烧室出口比热比初值γ4与燃烧室出口比热比算值γ4,ave的差值是否小于第一预设差值,若否则将燃烧室出口比热比算值γ4,ave赋值于燃烧室出口比热比初值γ4,重新计算燃烧室出口气体质量平均定压比热C4p,ave和燃烧室出口比热比算值γ4,ave,直至燃烧室出口比热比初值γ4与燃烧室出口比热比算值γ4ave的差值小于第一预设差值;对喷管出口比热比进行迭代的步骤包括:比较喷管出口比热比初值γ8与喷管出口比热比算值γ8,ave的差值是否小于第二预设差值,若否则将喷管出口比热比算值γ8,ave赋值于喷管出口比热比初值γ8,重新计算喷管出口气体质量平均定压比热C8p,ave和喷管出口比热比算值γ8,ave,直至喷管出口比热比初值γ8与喷管出口比热比算值γ8,ave的差值小于第二预设差值;对燃烧室燃烧效率进行迭代的步骤包括:比较燃烧室燃烧效率算值ηc与燃烧效率初值η0的差值是否小于第三预设差值,若否则将燃烧室燃烧效率算值ηc赋值于燃烧效率初值η0,根据燃烧效率初值η0计算燃烧室内各种物质组分,重新计算燃烧室内气体质量平均气体常数Rave、燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4、燃烧室出口静温Ts4、燃烧室出口气体质量平均定压比热C4p,ave、燃烧室出口比热比算值γ4,ave、喷管出口气体质量平均定压比热C8p,ave、喷管出口比热比算值γ8ave和燃烧室燃烧效率算值ηc,直至燃烧室燃烧效率算值ηc与燃烧效率初值η0的差值小于第三预设差值。
参见图1,第一预设差值、第二预设差值和第三预设差值皆可取值为e,e≤0.001。当满足收敛条件时结束相应计算步骤,并进行下一步骤,在燃烧室出口比热比、喷管出口比热比和燃烧室燃烧效率三者皆进行迭代计算时,燃烧室燃烧效率算值ηc的精度达到最高。
进一步的,计算燃烧室出口气体质量平均定压比热C4p,ave的步骤包括:根据公式,计算燃烧室出口气体质量平均定压比热C4p,ave,其中,/>为燃烧室内i组分气体的质量占比,/>为燃烧室出口处i组分的定压比热,,a4、b4、d4、e4和f4分别为多项式计算参数,可通过查询工程热力学文献获得。
进一步的,计算喷管出口总压Pt8、喷管出口总温Tt8和喷管出口静温Ts8的步骤包括:根据公式,计算排气喷管出口总压Pt8,其中,Ma8=1;根据公式/>,计算燃烧室出口总温Tt8,其中,A8为排气喷管出口面积;根据公式/>,计算燃烧室出口静温Ts8。
进一步的,计算喷管出口气体质量平均定压比热C8p,ave的步骤包括:根据公式,计算喷管出口气体质量平均定压比热C8p,ave,其中,/>为燃烧室内i组分气体的质量占比,/>为喷管出口处i组分的定压比热,,a8、b8、d8、e8和f8分别为多项式计算参数,可通过查询工程热力学文献获得。
如图1所示,流程示意图将上述分析方法的步骤整合、简化,可将其转换为相应的计算机语言植入分析系统的处理器。
如图2和图3所示,旋转爆震发动机,包括:喷注段1、扩压器2、火焰筒3、排气喷管4、喷嘴5和起爆装置6,喷注段1、扩压器2、火焰筒3和排气喷管4同轴并沿轴向依次设置,多个喷嘴5设置于喷注段1,且多个喷嘴5沿喷注段1的周向间隔设置。火焰筒3接近扩压器2的位置安装有起爆装置6,通过起爆装置6点火起爆,爆震波沿火焰筒3周向传递,以火焰筒3内环形腔作为燃烧室,并向排气喷管4传递,以使尾气经排气喷管4喷出。本发明实施例提供的旋转爆震发动机燃烧效率分析系统,包括:多个压力传感器7,多个压力传感器7分别连接于火焰筒3,且多个压力传感器7沿燃烧室出口周向均匀间隔分布;交互器,用于输入初始条件,并输出分析结果,初始条件包括:燃烧室出口比热比初值γ4、喷管出口比热比初值γ8、喷管出口马赫数Ma8、燃烧效率初值η0、获取燃烧室出口面积A4和排气喷管出口面积A8,分析结果包括燃烧室燃烧效率算值ηc;处理器,交互器和多个压力传感器7分别与处理器连接,处理器可植入上述任一实施方式记载的分析方法,并用于计算获得燃烧室出口静压ps4、燃烧室内各种物质组分、燃烧室内气体质量平均气体常数Rave、燃烧室出口气体质量平均定压比热C4p,ave、燃烧室出口气体质量平均定容比热C4v,ave、燃烧室出口比热比算值γ4,ave、喷管出口气体质量平均定压比热C8p,ave、喷管出口气体质量平均定容比热C8v,ave、喷管出口比热比算值γ8,ave和燃烧室燃烧效率算值ηc。
需要说明的是,压力传感器7的数量理论上越多,从而克服流场周向上非均匀性的能力更强,考虑到运算量和测试精度需要,优选将压力传感器7的数量设置为8个。
在本发明实施例中,旋转爆震发动机燃烧效率分析系统还包括:流量计和温度传感器,流量计和温度传感器分别与处理器连接;流量计安装于供油管路,温度传感器安装于旋转爆震发动机的喷注段。需要说明的是,流量计和温度传感器可采用旋转爆震发动机自带的器件,通过相应的接口引出其信号,并连接处理器。
采用上述旋转爆震发动机燃烧效率分析系统,开展测试时,确保发动机排气喷管达到临界状态,利用多个压力传感器7高频检测,可以实现时间和空间上的平均,进而获得平均静压。由可压缩流动的气体动力学函数,得到排气喷管出口的总压,进而计算获得排气喷管总温。根据燃烧室进口总温和排气喷管喉道总温,得到工质经爆震燃烧后的焓变,进而得到旋转爆震发动机的燃烧效率。不仅可以解决旋转爆震发动机工况时间和空间非定常性、难以准确分析燃烧效率的技术问题,而且可以提高燃烧效率分析的便捷性和准确性。
如图1、图2和图3所示,采用上述旋转爆震发动机燃烧效率分析方法及系统,具有以下有益效果:
(1)由于爆震波在燃烧室内沿周向高频传播,导致燃烧室在空间上存在高度的非均匀性,在时间上存在强非定常性,上述旋转爆震发动机燃烧效率分析方法及系统在燃烧室出口布置多个压力传感器7采集燃烧室出口静压,解决了燃烧室出口静压的空间非均匀性问题。对多个压力传感器7采集的多个周期静压进行时间平均,得到燃烧室出口静压,解决了静压测量的时间非定常性问题。
(2)通过测量静压,并根据气动动力学函数,得到燃烧室出口总温,根据温升法得到燃烧室的燃烧效率。测量仪器只用到了压力传感器7,不需要布置热电偶,省去了水冷装置,测试装置简单,燃烧效率分析更准确。
(3)克服了燃气的比热和比热比等物性参数,随燃气温度和燃气组分的变化而产生较大变化的技术问题。上述实施方式提供的旋转爆震发动机燃烧效率分析方法及系统以燃烧室出口静压和燃烧室流量为约束,迭代求解燃气的比热比等物性参数和燃气的温度等气动热力学参数,进而求解燃烧室燃烧效率。考虑温度对物性参数比热比和气体常数的影响,燃烧室出口总压、总温计算更加准确,燃烧室燃烧效率计算准确度更高。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (9)
1.一种旋转爆震发动机燃烧效率分析方法,其特征在于,包括以下步骤:
获取沿燃烧室出口周向均匀间隔分布的多个测试位置在爆震波稳定运动时的压力pi,t,其中,i为测试位置序号,t为压力采集时间;
根据公式,计算燃烧室出口静压ps4,其中,x为压力测试位置总数;
根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4;
根据温升法计算燃烧室燃烧效率算值ηc;
所述根据气体动力学函数计算燃烧室出口马赫数Ma4、燃烧室出口总压Pt4、燃烧室出口总温Tt4和燃烧室出口静温Ts4的步骤包括:
计算燃烧室内气体质量平均气体常数Rave;
获取燃烧室出口比热比初值γ4、喷管出口比热比初值γ8和喷管出口马赫数Ma8;
获取燃烧室出口面积A4和排气喷管出口面积A8;
根据公式,计算所述燃烧室出口马赫数Ma4,其中,Ma8=1;
根据公式,计算所述燃烧室出口总压Pt4;
根据公式,计算所述燃烧室出口总温Tt4,其中,mg为燃烧室内的油气混合物流量;
根据公式,计算所述燃烧室出口静温Ts4。
2.根据权利要求1所述的旋转爆震发动机燃烧效率分析方法,其特征在于,所述计算燃烧室内气体质量平均气体常数Rave的步骤包括:
获取燃烧效率初值η0;
根据公式,计算参与燃烧的燃油流量mf,c;
根据公式,计算参与燃烧的空气流量ma,c,L0为完全燃烧单位燃油需要的理论空气量;
根据公式,计算未参与燃烧的燃油流量mf,nc;
根据公式,计算未参与燃烧的空气流量ma,nc;其中,ma为进入燃烧室的空气流量,mf为进入燃烧室的燃油流量;
设定燃油平均化学式为CmHn,且燃烧产物为二氧化碳和水蒸气;
根据公式,计算燃烧产物中二氧化碳的流量/>,其中,为二氧化碳分子质量,/>为燃油分子质量;
根据公式,计算燃烧产物中的水蒸气流量/>,其中,为水蒸气分子质量;
根据公式,计算未参与燃烧空气中的氮气流量/>,其中,/>为氮气的质量分数;
根据公式,计算未参与燃烧空气中的氧气流量/>,其中,/>为氧气的质量分数;
根据公式,计算燃烧室内各组分气体的质量占比/>,其中,j组分包含氮气、氧气、二氧化碳、水蒸气和煤油蒸汽,n0=5;
根据公式,计算所述燃烧室内气体质量平均气体常数Rave,其中,Rj为j组分的气体常数。
3.根据权利要求1所述的旋转爆震发动机燃烧效率分析方法,其特征在于,所述根据温升法计算燃烧室燃烧效率算值ηc的步骤包括:
计算燃烧室出口气体质量平均定压比热C4p,ave;
根据公式,计算燃烧室出口比热比算值γ4,ave,其中,C4p,ave为燃烧室出口气体质量平均定压比热,C4v,ave为燃烧室出口气体质量平均定容比热,且C4v,ave=C4p,ave-Rave,其中,Rave为燃烧室内气体质量平均气体常数;
计算喷管出口总压Pt8、喷管出口总温Tt8和喷管出口静温Ts8;
计算喷管出口气体质量平均定压比热C8p,ave;
根据公式,计算喷管出口比热比算值γ8ave,其中,C8p,ave为喷管出口气体质量平均定压比热,C8v,ave为喷管出口气体质量平均定容比热,且C8v,ave=C8p,ave-Rave;
根据公式,计算燃烧室燃烧效率算值ηc,其中,为空气的定压比热,/>为燃油的定压比热,Tt2为喷注段入口空气总温,ma为进入燃烧室的空气流量,mf为进入燃烧室的燃油流量,mg为燃烧室内的油气混合物流量,mg=ma+mf,LHV为燃油的低热值,Tt4为燃烧室出口总温。
4.根据权利要求3所述的旋转爆震发动机燃烧效率分析方法,其特征在于,对燃烧室出口比热比、喷管出口比热比和燃烧室燃烧效率三者中至少其一进行迭代;
所述对燃烧室出口比热比进行迭代的步骤包括:比较所述燃烧室出口比热比初值γ4与所述燃烧室出口比热比算值γ4,ave的差值是否小于第一预设差值,若否则将所述燃烧室出口比热比算值γ4,ave赋值于所述燃烧室出口比热比初值γ4,重新计算所述燃烧室出口气体质量平均定压比热C4p,ave和所述燃烧室出口比热比算值γ4,ave,直至所述燃烧室出口比热比初值γ4与所述燃烧室出口比热比算值γ4ave的差值小于所述第一预设差值;
所述对喷管出口比热比进行迭代的步骤包括:比较所述喷管出口比热比初值γ8与所述喷管出口比热比算值γ8,ave的差值是否小于第二预设差值,若否则将所述喷管出口比热比算值γ8,ave赋值于所述喷管出口比热比初值γ8,重新计算所述喷管出口气体质量平均定压比热C8p,ave和所述喷管出口比热比算值γ8,ave,直至所述喷管出口比热比初值γ8与所述喷管出口比热比算值γ8,ave的差值小于所述第二预设差值;
所述对燃烧室燃烧效率进行迭代的步骤包括:比较所述燃烧室燃烧效率算值ηc与燃烧效率初值η0的差值是否小于第三预设差值,若否则将所述燃烧室燃烧效率算值ηc赋值于所述燃烧效率初值η0,根据所述燃烧效率初值η0计算燃烧室内各种物质组分,重新计算所述燃烧室内气体质量平均气体常数Rave、所述燃烧室出口马赫数Ma4、所述燃烧室出口总压Pt4、所述燃烧室出口总温Tt4、所述燃烧室出口静温Ts4、所述燃烧室出口气体质量平均定压比热C4p,ave、所述燃烧室出口比热比算值γ4,ave、所述喷管出口气体质量平均定压比热C8p,ave、所述喷管出口比热比算值γ8ave和所述燃烧室燃烧效率算值ηc,直至所述燃烧室燃烧效率算值ηc与所述燃烧效率初值η0的差值小于第三预设差值。
5.根据权利要求3所述的旋转爆震发动机燃烧效率分析方法,其特征在于,所述计算燃烧室出口气体质量平均定压比热C4p,ave的步骤包括:
根据公式,计算燃烧室出口气体质量平均定压比热C4p,ave,其中,Yi为燃烧室内i组分气体的质量占比,/>为燃烧室出口处i组分的定压比热,,a4、b4、d4、e4和f4分别为多项式计算参数,Rj为j组分的气体常数。
6.根据权利要求3所述的旋转爆震发动机燃烧效率分析方法,其特征在于,所述计算喷管出口总压Pt8、喷管出口总温Tt8和喷管出口静温Ts8的步骤包括:
根据公式,计算所述排气喷管出口总压Pt8,其中,Ma8=1,γ8为喷管出口比热比初值;
根据公式,计算所述喷管出口总温Tt8,其中,A8为排气喷管出口面积;
根据公式,计算所述喷管出口静温Ts8。
7.根据权利要求3所述的旋转爆震发动机燃烧效率分析方法,其特征在于,所述计算喷管出口气体质量平均定压比热C8p,ave的步骤包括:
根据公式,计算喷管出口气体质量平均定压比热C8p,ave,其中,为燃烧室内i组分气体的质量占比,/>为喷管出口处i组分的定压比热,,a8、b8、d8、e8和f8分别为多项式计算参数,Rj为j组分的气体常数。
8.一种旋转爆震发动机燃烧效率分析系统,其特征在于,包括:
多个压力传感器(7),多个所述压力传感器(7)分别连接于火焰筒(3),且多个所述压力传感器(7)沿燃烧室出口周向均匀间隔分布;
交互器,用于输入初始条件,并输出分析结果,所述初始条件包括:燃烧室出口比热比初值γ4、喷管出口比热比初值γ8、喷管出口马赫数Ma8、燃烧效率初值η0、获取燃烧室出口面积A4和排气喷管出口面积A8,所述分析结果包括燃烧室燃烧效率算值ηc;
处理器,所述交互器和多个所述压力传感器(7)分别与所述处理器连接,所述处理器用于计算获得燃烧室出口静压ps4、燃烧室内各种物质组分、燃烧室内气体质量平均气体常数Rave、燃烧室出口气体质量平均定压比热C4p,ave、燃烧室出口气体质量平均定容比热C4v,ave、燃烧室出口比热比算值γ4,ave、喷管出口气体质量平均定压比热C8p,ave、喷管出口气体质量平均定容比热C8v,ave、喷管出口比热比算值γ8,ave和燃烧室燃烧效率算值ηc;
其中,计算获得燃烧室出口静压ps4和燃烧室燃烧效率算值ηc的步骤包括:
获取多个所述压力传感器(7)在爆震波稳定运动时的压力pi,t,其中,i为测试位置序号,t为压力采集时间;
根据公式,计算燃烧室出口静压ps4,其中,x为压力测试位置总数;
计算燃烧室内气体质量平均气体常数Rave;
获取燃烧室出口比热比初值γ4、喷管出口比热比初值γ8和喷管出口马赫数Ma8;
获取燃烧室出口面积A4和排气喷管出口面积A8;
根据公式,计算所述燃烧室出口马赫数Ma4,其中,Ma8=1;
根据公式,计算所述燃烧室出口总压Pt4;
根据公式,计算所述燃烧室出口总温Tt4,其中,mg为燃烧室内的油气混合物流量;
根据公式,计算所述燃烧室出口静温Ts4;
根据温升法计算燃烧室燃烧效率算值ηc。
9.根据权利要求8所述的旋转爆震发动机燃烧效率分析系统,其特征在于,还包括:流量计和温度传感器,所述流量计和所述温度传感器分别与所述处理器连接;
所述流量计安装于供油管路,所述温度传感器安装于旋转爆震发动机的喷注段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310835904.5A CN116562193B (zh) | 2023-07-10 | 2023-07-10 | 旋转爆震发动机燃烧效率分析方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310835904.5A CN116562193B (zh) | 2023-07-10 | 2023-07-10 | 旋转爆震发动机燃烧效率分析方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116562193A CN116562193A (zh) | 2023-08-08 |
CN116562193B true CN116562193B (zh) | 2023-10-10 |
Family
ID=87488372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310835904.5A Active CN116562193B (zh) | 2023-07-10 | 2023-07-10 | 旋转爆震发动机燃烧效率分析方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116562193B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539480A (zh) * | 2009-04-30 | 2009-09-23 | 哈尔滨工业大学 | 超燃冲压发动机的燃烧效率的一维评价方法 |
CN109404163A (zh) * | 2018-12-13 | 2019-03-01 | 西安航天动力研究所 | 测量冲压发动机燃烧室气膜冷却套内流速的装置及方法 |
CN115901268A (zh) * | 2022-11-08 | 2023-04-04 | 中国航发沈阳发动机研究所 | 一种在发动机上准确获取燃烧室总压损失系数的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2447707B (en) * | 2007-03-23 | 2009-11-11 | Thermo Fisher Scientific Inc | Combustion tube and method for combusting a sample for combustion analysis |
US11725815B2 (en) * | 2019-12-06 | 2023-08-15 | GTL Company | Apparatuses, systems, and methods for optimizing acoustic wave confinement to increase combustion efficiency |
-
2023
- 2023-07-10 CN CN202310835904.5A patent/CN116562193B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101539480A (zh) * | 2009-04-30 | 2009-09-23 | 哈尔滨工业大学 | 超燃冲压发动机的燃烧效率的一维评价方法 |
CN109404163A (zh) * | 2018-12-13 | 2019-03-01 | 西安航天动力研究所 | 测量冲压发动机燃烧室气膜冷却套内流速的装置及方法 |
CN115901268A (zh) * | 2022-11-08 | 2023-04-04 | 中国航发沈阳发动机研究所 | 一种在发动机上准确获取燃烧室总压损失系数的方法 |
Non-Patent Citations (3)
Title |
---|
Experimental investigation on combustion efficiency of a partially premixed kerosene-air rotating detonation combustor;Shida Xu 等;《Fuel》;全文 * |
不同当量比下喷管对旋转爆震特性的影响研究;王顺利 等;《爆炸与冲击》;全文 * |
吸气式旋转爆震发动机热力循环过程分析与性能计算;李冬 等;《推进技术》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116562193A (zh) | 2023-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bach et al. | Performance analysis of a rotating detonation combustor based on stagnation pressure measurements | |
Zhang et al. | Effects of inner cylinder length on H2/air rotating detonation | |
AU777846B2 (en) | System and method for direct non-intrusive measurement of corrected airflow | |
CN108562440B (zh) | 空气流量分区控制的燃烧室试验装置及其试验方法 | |
JP2008107346A (ja) | 燃焼挙動の時間領域での解析法及びシステム | |
Bach et al. | Impact of outlet restriction on RDC performance and stagnation pressure rise | |
CN111157248A (zh) | 基于地面直连试验的冲压发动机及其燃烧室性能评估方法 | |
Hu et al. | Investigations on quantitative measurement of heat release rate using chemiluminescence in premixed methane-air flames | |
Fievisohn et al. | Experimental measurements of equivalent available pressure-lessons learned | |
Nguyen et al. | Study of mixing in a swirling jet | |
Zhang et al. | Experimental and numerical investigation of ultra-wet methane combustion technique for power generation | |
Dovizio et al. | RANS simulations of a series of turbulent V-shaped flames using conditional source-term estimation | |
CN112943484B (zh) | 一种研究粗糙度对喷管壁面传热影响的实验装置 | |
CN116562193B (zh) | 旋转爆震发动机燃烧效率分析方法及系统 | |
Hield et al. | Comparison of open and choked premixed combustor exits during thermoacoustic limit cycle | |
CN116127815B (zh) | 一种带引射喷管涡扇发动机的建模方法 | |
RU2609819C1 (ru) | Устройство для проведения высокотемпературных газодинамических испытаний проточных элементов турбомашин | |
CN111024359B (zh) | 一种短时气体喷注流量测量方法 | |
CN116562194B (zh) | 冲压旋转爆震发动机推力评估方法及系统 | |
Weber et al. | Choked Flow in a Converging/Diverging Rotating Detonation Engine Exit | |
Beard et al. | Mass flow rate measurement in a transonic turbine test facility with temperature distortion and swirl | |
CN116593168B (zh) | 冲压旋转爆震发动机耗油率评估方法及系统 | |
CN113432149B (zh) | 一种燃油喷嘴改型的试验验证方法和系统 | |
Cadiou | Experimental method for the combustion efficiency calculation in a reheat duct | |
Gallis et al. | Development and Validation of a 0-D/1-D Model to Evaluate Pulsating Conditions from a Constant Volume Combustor |
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 |