CN112539941B - 考虑真实气体效应的液体火箭发动机热试验参数设置方法 - Google Patents
考虑真实气体效应的液体火箭发动机热试验参数设置方法 Download PDFInfo
- Publication number
- CN112539941B CN112539941B CN202011405279.3A CN202011405279A CN112539941B CN 112539941 B CN112539941 B CN 112539941B CN 202011405279 A CN202011405279 A CN 202011405279A CN 112539941 B CN112539941 B CN 112539941B
- Authority
- CN
- China
- Prior art keywords
- gas
- turbine
- compression factor
- turbine power
- thermal test
- 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
Classifications
-
- 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/14—Testing gas-turbine engines or jet-propulsion engines
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供了一种考虑真实气体效应的液体火箭发动机热试验参数设置方法,旨在避免出现涡轮功率计算的偏差。该方法以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,通过气体压缩因子的线性加权对理想气体的涡轮功率进行修正;其中,气体压缩因子基于RKS状态方程进行求解;基于计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。本发明充分考虑了真实气体效应对涡轮功率的影响,设置的发动机联试和热试工况参数更为准确,从而能够降低试验风险。
Description
技术领域
本发明涉及一种液体火箭发动机热试验参数设置方法,特别涉及高温高压下涡轮功率有效的工程计算。
背景技术
大推力液氧煤油发动机是我国未来深空探测、载人登月和开展大型空间活动合理的动力选择之一。在燃气发生器/涡轮泵联式之前,需要确定试车工况参数,即进行液体火箭发动机热试验参数设置,包括涡轮功率、涡轮泵转速、燃气发生器温度等。其中,首先需要准确设置涡轮功率,进而再根据流量、压力、功率平衡计算其他热试验参数。
随着液氧煤油发动机推力和性能的提升,燃气发生器的压力和温度相应增加,导致燃气的性质偏离理想气体,真实气体效应突出。现有的液体火箭发动机热试验参数设置方法没有充分考虑到真实气体效应对液体火箭发动机热试验参数带来的影响。
发明内容
本发明的目的是为了准确设置发动机联试和热试工况参数,降低试验风险。为此,本发明充分考虑了真实气体效应对涡轮功率的影响,探寻高温高压下涡轮功率有效的工程计算方法。
真实气体与理想气体的偏差可以用压缩因子来描述。在高温高压或低温条件下,与气体性质和工作条件相关,气体的压缩因子向更大或者更小的方向偏离。当压力和温度很高时,不考虑真实气体效应将导致涡轮功率的偏差,进而导致涡轮设计和涡轮泵功率平衡的偏差,使得涡轮泵转速、燃气发生器温度等参数偏差过大,进而影响发动机联试和热试工况参数设置以及试验的安全防护***设计。
本发明的方案如下:
方案一:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,利用涡轮膨胀过程的压缩因子的线性加权对理想气体的涡轮功率进行修正;其中,气体压缩因子基于RKS状态方程进行求解;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
上述步骤1)的涡轮功率计算方式如下:
Pt,r=(αZ1+(1-α)Z2)Pt,i (1)
式(1)中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;
压缩因子Z1、Z2均基于RKS状态方程来计算,具体如下:
Tc为临界温度,K;pc为临界压力,kPam=0.480+1.574ω-0.176ω2,ω为偏心因子;Tr=T/Tc;
所述压缩因子Z1、Z2统一用Z表示,Z=pv/(RT);式(2)转换为以下式(3),分别相应计算得出涡轮入口气体的压缩因子Z1和涡轮出口气体的压缩因子Z2;
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式(3)中,A=ap/(RT)2,B=bp/(RT)。
进一步地,所述权系数α根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化确定。
进一步地,对于临界范围工作的涡轮,权系数α范围为0.4~0.6。
可选地,工作气体为氧气,则偏心因子ω取0.0222。
方案二:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功和质量流量整体的影响,利用涡轮膨胀过程的压缩因子对理想气体的涡轮功率进行修正;其中,对于真实气体效应对涡轮绝热功的影响,是通过压缩因子线性加权进行修正;对于真实气体效应对涡轮质量流量的影响,具体是考虑真实气体效应对涡轮临界参数(即马赫数为1时的气体状态参数)的影响;压缩因子基于RKS状态方程进行求解;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
上述步骤1)的涡轮功率计算方式如下:
式中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;Zcr为临界状态的压缩因子;
该式中,权系数α的优化方法以及压缩因子Z1、Z2、Zcr的计算方法与上述方案一相同。
对于临界状态的压缩因子Zcr,参照以上式(2)、(3)计算,其中的p、T、v相应取工作气体的临界参数pc、Tc、vc。若工作气体为氧气,则偏心因子ω取0.0222,pc取5.043MPa,Tc取154.58K,vc=1/ρc,ρc=436.1kg/m3。
本发明具有以下有益效果:
本发明充分考虑了真实气体效应对涡轮功率的影响,设置的发动机联试和热试工况参数更为准确,从而能够降低试验风险。
具体实施方式
以下通过实施例进一步详述本发明。
对于高压涡轮而言,基于理想气体的涡轮理论绝热功偏差达到10%-20%。这会导致涡轮功率的偏差,进而导致液体火箭发动机静态参数平衡偏差,使得涡轮泵转速以及涡轮出口温度偏差过大,从而导致试验数据与理论设定工况参数偏差过大。
为此,以下实施例充分考虑了真实气体效应对涡轮功率的影响,基于理想气体的涡轮功率的计算进行了修正。
涡轮功率计算方法一:
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,从而进行修正。修正公式为气体压缩因子的线性加权,权系数处于0~1之间。最优系数通过模拟退火算法进行确定。压缩因子由RKS状态方程进行求解。
Pt,r=(αZ1+(1-α)Z2)Pt,i (1)
式(1)中,Pt,r为考虑真实气体效应的涡轮功率,Pt,i为基于理想气体的涡轮功率,α为权系数,Z1为涡轮入口气体的压缩因子,Z2为涡轮出口气体的压缩因子;
权系数α范围为0~1。根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化求解。对于临界范围工作的涡轮,优化后,α范围为0.4~0.6。
对于高压富氧燃气为工质的涡轮,压缩因子Z1、Z2利用RKS状态方程(真实气体状态方程)来描述精度更好,计算量小适合工程计算;具体如下:
m=0.480+1.574ω-0.176ω2,Tr=T/Tc,b=0.08664RTc/pc
式中:p-压力,kPa;R-通用气体常数(8.314kJ/(kmolK);T-温度,K;v-比容(m3/kmol);pc-临界压力,kPa;Tc-临界温度;K;ω-偏心因子,为常数(若工作气体为氧气,则取0.0222);
对于涡轮入口气体的压缩因子Z1,这里的参数p、T、v即涡轮入口气体的压力、温度、比容;同理,对于涡轮出口气体的压缩因子Z2,这里的参数p、T、v即涡轮出口气体的压力、温度、比容。
用压缩因子Z统一表示,式(2)可改写为
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式中:A=ap/(RT)2,B=bp/(RT),Z=pv/(RT)。
涡轮功率计算方法二:
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功和质量流量整体的影响,从而进行修正。首先计算理想气体的涡轮功率,然后利用涡轮膨胀过程的压缩因子进行修正。修正公式分为两部分,一部分考虑真实气体效应对涡轮绝热功的影响进行修正,第二部分为考虑真实气体效应对涡轮质量流量影响的修正。第一部分与上述计算方法一相同,第二部分考虑真实气体效应对涡轮临界参数的影响,从而修正涡轮质量流量。
所以,以上公式(1)改为:
式(4)中,Zcr为临界状态的压缩因子;其他参数符号含义、权系数α的优化方法以及压缩因子的计算方法与上述计算方法一相同。
对于临界状态的压缩因子Zcr,参照以上式(2)、(3)计算,其中的p、T、v相应取工作气体的临界参数pc、Tc、vc。若工作气体为氧气,则偏心因子ω取0.0222,pc取5.043MPa,Tc取154.58K,vc=1/ρc,ρc=436.1kg/m3。
基于以上计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。这些液体火箭发动机静态参数平衡计算的具体方法与常规方法相同,故不再赘述。
本实施例根据液氧煤油发动机联试测量的涡轮流量、入口温度、入口压力、出口压力、效率等参数计算考虑真实气体效应的涡轮功率,并与涡轮功率试验数据对比,验证了计算方法的正确性。
Claims (8)
1.考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,包括:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功的影响,利用涡轮膨胀过程的压缩因子的线性加权对理想气体的涡轮功率进行修正;其中,气体压缩因子基于RKS状态方程进行求解;
涡轮功率计算方式如下:
Pt,r=(αZ1+(1-α)Z2)Pt,i (1)
式(1)中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;
压缩因子Z1、Z2均基于RKS状态方程来计算,具体如下:
Tc为临界温度,K;pc为临界压力,kPa;m=0.480+1.574ω-0.176ω2,ω为偏心因子;Tr=T/Tc;
所述压缩因子Z1、Z2统一用Z表示,Z=pv/(RT);式(2)转换为以下式(3),分别相应计算得出涡轮入口气体的压缩因子Z1和涡轮出口气体的压缩因子Z2;
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式(3)中,A=ap/(RT)2,B=bp/(RT);
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
2.根据权利要求1所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,所述权系数α根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化确定。
3.根据权利要求2所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,对于临界范围工作的涡轮,权系数范围为0.4~0.6。
4.根据权利要求1所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,工作气体为氧气,则偏心因子ω取0.0222。
5.考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,包括:
1)计算涡轮功率
以理想气体的涡轮功率为基础,考虑真实气体效应对涡轮绝热功和质量流量整体的影响,利用涡轮膨胀过程的压缩因子对理想气体的涡轮功率进行修正;其中,对于真实气体效应对涡轮绝热功的影响,是通过压缩因子线性加权进行修正;对于真实气体效应对涡轮质量流量的影响,具体是考虑真实气体效应对涡轮临界参数的影响;压缩因子基于RKS状态方程进行求解;
涡轮功率计算方式如下:
式(1)中,Pt,r为考虑真实气体效应的涡轮功率;Pt,i为基于理想气体的涡轮功率;α为权系数,范围为0~1;Z1为涡轮入口气体的压缩因子;Z2为涡轮出口气体的压缩因子;Zcr为临界状态的压缩因子;
压缩因子Z1、Z2均基于RKS状态方程来计算,具体如下:
Tc为临界温度,K;pc为临界压力,kPa;m=0.480+1.574ω-0.176ω2,ω为偏心因子;Tr=T/Tc;
所述压缩因子Z1、Z2统一用Z表示,Z=pv/(RT);式(2)转换为以下式(3),分别相应计算得出涡轮入口气体的压缩因子Z1和涡轮出口气体的压缩因子Z2;
Z3-Z2+(A-B-B2)Z-AB=0 (3)
式(3)中,A=ap/(RT)2,B=bp/(RT);
对于临界状态的压缩因子Zcr,参照以上式(2)、(3)计算,其中的p、T、v相应取工作气体的临界参数pc、Tc、vc;
2)计算其他参数
基于步骤1)计算得到的涡轮功率,再根据流量、压力、功率平衡计算其他热试验参数。
6.根据权利要求5所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,所述权系数α根据液体火箭发动机涡轮入口压力和温度的范围,利用模拟退火算法进行优化确定。
7.根据权利要求6所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,对于临界范围工作的涡轮,权系数范围为0.4~0.6。
8.根据权利要求5所述的考虑真实气体效应的液体火箭发动机热试验参数设置方法,其特征在于,工作气体为氧气,则偏心因子ω取0.0222,pc取5.043MPa,Tc取154.58K,vc=1/ρc,ρc=436.1kg/m3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011405279.3A CN112539941B (zh) | 2020-12-02 | 2020-12-02 | 考虑真实气体效应的液体火箭发动机热试验参数设置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011405279.3A CN112539941B (zh) | 2020-12-02 | 2020-12-02 | 考虑真实气体效应的液体火箭发动机热试验参数设置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112539941A CN112539941A (zh) | 2021-03-23 |
CN112539941B true CN112539941B (zh) | 2023-01-20 |
Family
ID=75015855
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011405279.3A Active CN112539941B (zh) | 2020-12-02 | 2020-12-02 | 考虑真实气体效应的液体火箭发动机热试验参数设置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112539941B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114970394B (zh) * | 2022-06-02 | 2023-03-21 | 西安航天动力研究所 | 一种高压补燃发动机混合燃气涡轮绝热功计算方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4845941A (en) * | 1986-11-07 | 1989-07-11 | Paul Marius A | Gas turbine engine operating process |
CN101949324A (zh) * | 2010-08-17 | 2011-01-19 | 清华大学 | 一种内燃机涡轮增压通流匹配方法 |
CN103116705A (zh) * | 2013-02-06 | 2013-05-22 | 中国航天科技集团公司第六研究院第十一研究所 | 一种补燃循环火箭发动机故障仿真分析方法 |
CN104272207A (zh) * | 2012-03-01 | 2015-01-07 | 诺沃皮尼奥内股份有限公司 | 用于实时燃气涡轮性能报告的方法和*** |
CN105074413A (zh) * | 2013-02-05 | 2015-11-18 | 西门子股份公司 | 用于燃气轮机的自动测试*** |
CN107702854A (zh) * | 2017-07-10 | 2018-02-16 | 北京航天动力研究所 | 一种满足火箭发动机涡轮泵装配及使用的高速动平衡方法 |
KR20180119417A (ko) * | 2017-04-25 | 2018-11-02 | 한국전력공사 | 터빈 시험장치 |
CN109376445A (zh) * | 2018-11-07 | 2019-02-22 | 北京动力机械研究所 | 燃气涡轮发动机起动建模方法 |
CN111473941A (zh) * | 2020-04-29 | 2020-07-31 | 北京动力机械研究所 | 一种闭式循环涡轮特性试验方法 |
CN111652418A (zh) * | 2020-05-25 | 2020-09-11 | 西安西热节能技术有限公司 | 一种火电机组动态精细化复合参数滑压曲线生成方法 |
-
2020
- 2020-12-02 CN CN202011405279.3A patent/CN112539941B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4845941A (en) * | 1986-11-07 | 1989-07-11 | Paul Marius A | Gas turbine engine operating process |
CN101949324A (zh) * | 2010-08-17 | 2011-01-19 | 清华大学 | 一种内燃机涡轮增压通流匹配方法 |
CN104272207A (zh) * | 2012-03-01 | 2015-01-07 | 诺沃皮尼奥内股份有限公司 | 用于实时燃气涡轮性能报告的方法和*** |
CN105074413A (zh) * | 2013-02-05 | 2015-11-18 | 西门子股份公司 | 用于燃气轮机的自动测试*** |
CN103116705A (zh) * | 2013-02-06 | 2013-05-22 | 中国航天科技集团公司第六研究院第十一研究所 | 一种补燃循环火箭发动机故障仿真分析方法 |
KR20180119417A (ko) * | 2017-04-25 | 2018-11-02 | 한국전력공사 | 터빈 시험장치 |
CN107702854A (zh) * | 2017-07-10 | 2018-02-16 | 北京航天动力研究所 | 一种满足火箭发动机涡轮泵装配及使用的高速动平衡方法 |
CN109376445A (zh) * | 2018-11-07 | 2019-02-22 | 北京动力机械研究所 | 燃气涡轮发动机起动建模方法 |
CN111473941A (zh) * | 2020-04-29 | 2020-07-31 | 北京动力机械研究所 | 一种闭式循环涡轮特性试验方法 |
CN111652418A (zh) * | 2020-05-25 | 2020-09-11 | 西安西热节能技术有限公司 | 一种火电机组动态精细化复合参数滑压曲线生成方法 |
Non-Patent Citations (1)
Title |
---|
真实气体效应对高压涡轮理论绝热功的影响;秦艳平 等;《中国航天第三专业信息网第四十届技术交流会暨第四届空天动力联合会议论文集》;20190831;第527-533页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112539941A (zh) | 2021-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104481704B (zh) | 一种实现燃气轮机起动过程中燃料实时控制方法及装置 | |
CN103069130B (zh) | 优化化学计量燃烧的***和方法 | |
CN112539941B (zh) | 考虑真实气体效应的液体火箭发动机热试验参数设置方法 | |
CN109902426A (zh) | 空冷式高参数汽轮机高中压缸性能与结构设计及监控方法 | |
Trasino et al. | Modeling and performance analysis of the Rolls-Royce fuel cell systems limited: 1 MW plant | |
CN104160131A (zh) | 构造为调整电力输出的燃气涡轮机发动机 | |
CN113536192B (zh) | 一种基于加力燃烧室油气比计算加力总燃油流量的方法 | |
Zhou et al. | Performance analysis and optimized control strategy for a three-shaft, recuperated gas turbine with power turbine variable area nozzle | |
CN114060143B (zh) | 基于旁通阀流通特性的变海拔增压压力稳定性控制方法 | |
Dellenback | Improved gas turbine efficiency through alternative regenerator configuration | |
CA3000664C (en) | Caes combustor test facility | |
CN112115576A (zh) | 一种轴流压缩机缩尺过程多变效率修正方法和*** | |
CN109858129B (zh) | 一种关于三联供***的燃气轮机动态仿真方法 | |
WO2024120432A1 (zh) | 一种舰用燃气轮机起动性能模型建模方法 | |
De Paepe et al. | A study on the performance of steam injection in a typical micro Gas Turbine | |
Wang et al. | Control strategy for actual constraints during the start–stop process of a supercritical CO2 Brayton cycle | |
Larosa et al. | SOFC/MGT Coupling: Different Options With Standard Boosters | |
Ying et al. | Study on flow parameters optimisation for marine gas turbine intercooler system based on simulation experiment | |
CN115169048A (zh) | 一种基于多领域组件建模的重型燃气轮机建模方法 | |
Ishiyama et al. | Operating scenario of 3GWth class FFHR power plant with bypass controlled supercritical CO2 gas turbine power generation system | |
CN113806680B (zh) | 一种计算汽轮机进汽压损对汽轮机内效率修正量的方法 | |
Raggio et al. | Transient Analysis of a Micro Gas Turbine With Fuel Composition Change | |
CN108644020A (zh) | 一种改进的燃气轮机排气温度控制方法 | |
CN204283626U (zh) | 一种燃气轮机起动过程中燃料闭环控制装置 | |
Ali et al. | A review on application of hydrogen in gas turbines with intercooler adjustments |
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 |