CN113790899A - 液体火箭发动机喷注器燃烧状态评估装置及方法 - Google Patents
液体火箭发动机喷注器燃烧状态评估装置及方法 Download PDFInfo
- Publication number
- CN113790899A CN113790899A CN202111061147.8A CN202111061147A CN113790899A CN 113790899 A CN113790899 A CN 113790899A CN 202111061147 A CN202111061147 A CN 202111061147A CN 113790899 A CN113790899 A CN 113790899A
- Authority
- CN
- China
- Prior art keywords
- thermocouple
- heat flow
- hole
- injector
- temperature
- 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.)
- Granted
Links
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000007788 liquid Substances 0.000 title claims abstract description 25
- 238000003825 pressing Methods 0.000 claims abstract description 18
- 239000007921 spray Substances 0.000 claims abstract description 15
- 238000007789 sealing Methods 0.000 claims abstract description 6
- 238000011156 evaluation Methods 0.000 claims description 17
- 239000000463 material Substances 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 230000001052 transient effect Effects 0.000 claims description 9
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 4
- 238000009529 body temperature measurement Methods 0.000 claims description 4
- 229910052802 copper Inorganic materials 0.000 claims description 4
- 239000010949 copper Substances 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 229910001220 stainless steel Inorganic materials 0.000 claims description 3
- 239000010935 stainless steel Substances 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000009530 blood pressure measurement Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 abstract description 3
- 239000003380 propellant Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011900 installation process Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000005304 optical glass Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Combustion & Propulsion (AREA)
- Chemical & Material Sciences (AREA)
- Computing Systems (AREA)
- Testing Of Engines (AREA)
Abstract
本发明提供一种液体火箭发动机喷注器燃烧状态评估装置及方法,解决现有采用高速相机通过观察窗拍摄火焰,存在成本高、密封难度大、准备周期长的问题。装置包括连接组件、N个热电偶及依次同轴设的前法兰、身部、喷管和后法兰,N为≥2整数;身部上设有燃烧室,喷管上设有喉道;连接组件用于紧固前法兰和后法兰,实现身部和喷管的夹紧;身部沿径向开设有轴向均布的N个热电偶放置孔,热电偶放置孔包括外孔和直径小于外孔的内孔,内孔底面与燃烧室内壁面之间存在距离h,N个热电偶分别设在N个热电偶放置孔内,每个热电偶上设有与外孔配合的限位部以及套装有位于限位部外侧的弹簧;身部外壁设有压板,用于通过弹簧将热电偶压紧在热电偶放置孔中。
Description
技术领域
本发明涉及液体火箭发动机喷注器燃烧状态评估技术,具体涉及一种液体火箭发动机喷注器燃烧状态评估装置及方法。
背景技术
在液体火箭发动机中,喷注器是用来组织推进剂混合燃烧的,对发动机的整体性能有着至关重要的影响。喷注器的燃烧状态通常是用燃烧室压力来进行评估,但是燃烧室压力实际上是燃烧状态的一个综合参数,并不能反映燃烧状态的细节。因此,有学者提出以燃烧火焰长度来评估燃烧状态,火焰长度反映了推进剂从混合至完全燃烧所需要的轴向空间,对于燃烧室长度设计具有很重要的参考价值。但是火箭发动机燃烧室高温高压的恶劣环境对火焰长度的观测提出了巨大挑战。目前,通常采用在燃烧室侧壁开观察窗,用高速相机直接拍摄火焰。该方式虽可以直观的观察火焰长度,但是也存在观察窗所用光学玻璃价格昂贵、密封难度大、准备周期长等缺点。
发明内容
为了解决现有采用高速相机通过观察窗拍摄火焰,存在成本高、密封难度大、准备周期长的技术问题,本发明提供了一种液体火箭发动机喷注器燃烧状态评估装置及方法。
为实现上述目的,本发明提供的技术方案是:
一种液体火箭发动机喷注器燃烧状态评估装置,其特殊之处在于:包括连接组件、N个热电偶以及依次同轴设置的前法兰、身部、喷管和后法兰,N为大于等于2的整数;身部上设置有燃烧室,喷管上设置有与燃烧室同轴且连通的喉道;
所述连接组件用于紧固前法兰和后法兰,实现身部和喷管的夹紧;
所述身部外壁开设有与燃烧室连通的燃烧室测压孔,以及沿径向开设有轴向均布的N个热电偶放置孔,热电偶放置孔包括外孔和直径小于外孔的内孔,内孔底面与燃烧室内壁面之间存在距离h,h满足:
其中,p为燃烧室工作过程中的最高压力,d为内孔的直径,n为安全系数,可取2~5,σ为燃烧室材料的抗剪切强度,以上参数单位均为标准单位制。h在满足上述条件的前提下应尽量取小值;
相邻2个热电偶放置孔的轴线间隔为a,a满足:a≥2d′;
其中,d′为外孔的直径,在不影响安装热电偶的情况下,a的取值应尽量小。
所述N个热电偶分别设置在N个热电偶放置孔内,每个热电偶上设有与外孔配合的限位部,以及套装有位于限位部外侧的弹簧;
所述身部外壁设有热电偶放置孔的位置设有压板,用于通过弹簧将热电偶压紧在热电偶放置孔中。
进一步地,所述N个热电偶的外端部均伸出压板;
所述限位部为设置在热电偶外圆面的环形限位凸起,环形限位凸起设置在外孔上,弹簧的两端分别通过环形限位凸起外端面和压板内端面限位;
所述弹簧上套装有限位套筒,限位套筒的一端止靠于身部外壁,其另一端端与压板内端面之间存在距离。
进一步地,所述连接组件包括沿周向布置的M个锁紧件,M为大于等于2的整数,每个锁紧件包括长螺杆和2个第一螺母,长螺杆设置在前法兰和后法兰上,且两端分别伸出前法兰和后法兰后,分别与第一螺母连接。
进一步地,所述前法兰和后法兰均采用不锈钢材料,身部和喷管采用紫铜材料;
所述后法兰上安装有可拆卸的堵盖,用于喷注器点火试验前的密封性检测。
进一步地,所述前法兰上设有火花塞放置孔;
所述前法兰与身部之间、身部和喷管之间、喷管和后法兰之间均采用定位槽和定位块配合方式。
同时,本发明提供了一种液体火箭发动机喷注器燃烧状态评估方法,其特殊之处在于,包括以下步骤:
1)将待评估喷注器与评估装置的前法兰对接,热电偶安装到位,进行设定时间段的点火试验,测得所有热电偶在各采样时间点的瞬态壁温;
2)对瞬态壁温进行热流反演计算,获得燃烧状态稳定后的燃烧室热流分布特征曲线;
3)根据燃烧室热流分布特征曲线,确定热流分布特征曲线的平稳段起点,该平稳段起点为火焰结束位置,根据火焰结束位置获得火焰长度,进而评估喷注器燃烧状态。
进一步地,步骤2)具体如下:
2.1)假设身部的内壁面热流在圆周方向均匀分布,将用于热流反演的计算域简化为二维,时间步长取热电偶的采样时间间隔;
2.2)第一个时间步的起始时刻记为点火时刻,其初始温度场已知;
第一个时间步每个热电偶对应位置处的热流分别给定一个假设值,然后对计算域应用二维非稳态导热方程进行温度场计算,得到每个热电偶的计算温度,
根据二维非稳态导热方程计算得到计算温度,进而得到计算温度与热电偶实测瞬态壁温的温度残差,比较温度残差与设定阈值,若温度残差小于设定阈值,获得第一时间步的热流分布;若否,修改每个热电偶对应位置处热流的假设值进行迭代,直至温度残差小于设定阈值;
2.3)以第一时间步结束时刻的温度场为第二时间步的初始温度场,采用步骤2.2)的方法获得第二时间步的热流分布;
2.4)利用步骤2.3)的方法,获得剩余每个时间步的热流分布;
2.5)根据所有时间步的热流分布,获得燃烧状态稳定后的燃烧室热流分布特征曲线。
进一步地,步骤2.2)中;
计算温度残差具体为:设N个热电偶对应于N个位置的热流分布q=(q1,q2,...qN),N个热电偶的实测温度为Tm=(Tm1,Tm2,...TmN),在给定热流分布q条件下的计算温度为Tc(q)=(Tc1(q),Tc2(q),...TcN(q));温度残差按照下式计算:
J=[Tm-Tc(q)][Tm-Tc(q)]T;
设定阈值按照下式计算:
ε=N*(ΔT)2*Δt;
其中,ΔT为热电偶的测温精度,Δt为热电偶采样时间间隔。
修改每个热电偶对应位置处热流的假设值进行迭代具体为:根据下式计算修改后的热流分布qk+1;
S·[qk+1]T=[Tm-Tc(qk)]T+S·[qk]T
式中,qk为修改前的热流分布;S为温度对热流的响应矩阵,计算公式如下:
进一步地,步骤2.1)中,所述设定时间段为2s~5s。
进一步地,步骤3)中,根据火焰结束位置获得火焰长度,进而评估喷注器燃烧状态具体为:火焰结束位置越靠近上游,火焰长度越短,则喷注器燃烧越高效。
与现有技术相比,本发明的优点是:
1、本发明评估装置及方法仅依靠布置于燃烧室身部的热电偶,获得燃烧室内壁面热流分布,来估测火焰长度,即可实现对喷注器燃烧状态的评估,具有简单易行,试验成本低,周期短的特点。
2、本发明热电偶设计限位部,由安装于限位部上的弹簧对热电偶施加压紧力,弹簧上端的压板可用于调节压紧力,确保热电偶与壁面接触良好,同时也避免压紧力过大损坏测温端头。
3、本发明热电偶端头采用小尺寸露端设计,可减小热电偶温度延迟,提高相应速度。
4、本发明评估装置采用模块化设计,便于更换模块,适应不同的喷注器。
5、本发明评估方法采用热流迭代方式获得热流分布,相比于传统方法,具有收敛速度快,收敛稳定性好的特点,当数据量较大时,在时间效益方面具有明显优势。
附图说明
图1为本发明液体火箭发动机喷注器燃烧状态评估装置的结构示意图(未示出热电偶);
图2为本发明液体火箭发动机喷注器燃烧状态评估装置中热电偶放置孔处的结构示意图;
图3为本发明液体火箭发动机喷注器燃烧状态评估装置中热电偶处的结构示意图;
图4为本发明液体火箭发动机喷注器燃烧状态评估方法中热流反演用计算域示意图;
图5为本发明液体火箭发动机喷注器燃烧状态评估方法中热流反演流程图;
其中,附图标记如下:
1-喷注器,2-前法兰,3-身部,31-燃烧室,4-喷管,41-喉道,5-后法兰,6-堵盖,7-长螺杆,8-热电偶放置孔,81-外孔,82-内孔,9-燃烧室测压孔,10-火花塞放置孔,11-热电偶,12-第一螺母,13-压板,14-螺柱,15-第二螺母,16-弹簧,17-限位套筒,18-限位部。
具体实施方式
以下结合附图和具体实施例对本发明的内容作进一步详细描述。
如图1所示,本发明一种液体火箭发动机喷注器燃烧状态评估装置,包括N个热电偶11以及同轴设置的前法兰2、身部3、喷管4和后法兰5,N为大于等于2的整数;身部3和喷管4被前法兰2和后法兰5夹紧在中间,前法兰2和后法兰5通过连接组件实现紧固。本实施例为了实现前法兰2、身部3、喷管4、后法兰5相互之间定位的准确性,前法兰2与身部3之间、身部3和喷管4之间、喷管4和后法兰5之间均采用定位槽和定位块配合方式。
前法兰2和后法兰5均采用不锈钢材料;身部3和喷管4采用导热系数较高的紫铜材料,均为圆筒形,则身部3具有燃烧室31,喷管4具有与燃烧室31同轴且连通的喉道41。
本实施例连接组件包括沿周向布置的多个锁紧件,本实施例锁紧件采用4个,每个锁紧件包括1个长螺杆7和2个第一螺母12,长螺杆7设置在前法兰2和后法兰5上,且两端分别伸出前法兰2和后法兰5后,与第一螺母12连接,实现前法兰2和后法兰5的连接。
前法兰2用于与待评估喷注器1连接,喷注器1可替换,喷注器1通过螺钉与评估装置的前法兰2连接,并用铜垫片进行密封;后法兰5上安装有可拆卸的堵盖6,用于喷注器1点火试验前的密封检漏,在进行点火试验前拆除堵盖6。
前法兰2上设置有火花塞放置孔10,如果是自燃推进剂,则不需要点火装置,可以封堵该火花塞放置孔10。
身部3外壁开设有与燃烧室31连通的燃烧室测压孔9,以及沿径向开设有轴向均布的N个热电偶放置孔8。
如图2所示,热电偶放置孔8为盲孔,包括外孔81和内孔82两部分,外孔81的孔径较大,用于热电偶11的定位,内孔82较小,用于放置热电偶11的感温端头,测温位置处于内孔82的底部,内孔82底面与燃烧室31内壁面之间存在距离h,h满足:
其中,p为燃烧室工作过程中的最高压力,d为内孔的孔径,n为安全系数,可取2~5,σ为燃烧室材料的抗剪切强度,以上参数单位均为标准单位制。h在满足上述条件的前提下应尽量取小值。
相邻2个热电偶放置孔8的轴线间隔为a,a取值应大于等于外孔81直径的2倍,以保证热电偶11有足够的安装空间;在不影响安装热电偶11的情况下,a的取值应尽量小。本实施例h为1mm,a为10mm。
如图3所示,N个热电偶11分别设置在N个热电偶放置孔8内,每个热电偶11上设有与外孔81配合的限位部18,本实施例限位部18为设置在热电偶11外圆面的环形限位凸起,环形限位凸起的一端设置在外孔81上,其另一端伸出外孔81;为了提高热电偶11的响应速度,热电偶11采用露端形式,且采用较小的端头尺寸,端头直径约1mm。
身部3外壁设有热电偶放置孔8的位置设有压板13,热电偶11的外端部伸出压板13,热电偶11上套装有位于限位部18和压板之间的弹簧16,压板13将弹簧16压紧,实现热电偶11用弹簧16顶压在热电偶放置孔8中,保证测温过程中感温端头与内孔82底面接触良好。为了防止弹簧16受压后出现横向变形,本实施例在弹簧16外周设有限位套筒17。热电偶11的具体安装过程为:将两根螺柱14安装在燃烧室上,将带弹簧16的热电偶放入限位套筒17,一起放入热电偶放置孔8中,将压板13套在螺柱14上,压在弹簧16外端,调整合适的弹簧16压紧力,安装第二螺母15固定压板13,此时,热电偶11的环形限位凸起一端位于外孔81上,另一端伸出外孔81,限位套筒17的一端止靠于身部3外壁,另一端与压板13端面之间存在距离,实现热电偶11的压紧安装。
上述评估装置对喷注器1燃烧状态的评估方法,包括以下步骤:
1)将待评估喷注器1与评估装置的前法兰2对接,热电偶11安装到位,进行设定时间段的点火试验,测得所有热电偶11在各采样时间点的瞬态壁温;
其中,设定时间段不宜过短或过长,过短有可能燃烧状态尚未稳定,过长则燃烧室31无法承受,因此,本实施例设定时间段(试验时间)为2s~5s;
2)对通过热电偶11测得的瞬态壁温进行热流反演计算,获得燃烧状态稳定后的燃烧室热流分布特征曲线;
3)根据燃烧室热流分布特征曲线,确定热流分布特征曲线的平稳段起点,该平稳段起点为火焰结束位置,根据火焰结束位置获得火焰长度,进而评估喷注器1燃烧状态;其中,燃烧状态稳定时,燃烧室31热流沿轴向通常呈现先增加然后趋向平稳的分布特征,认为热流平稳段的起点就是火焰结束的位置,该位置越靠近上游,则说明推进剂混合至完全燃烧所需的轴向距离(火焰长度)越短,喷注器1燃烧组织越高效。
步骤2)中,瞬态壁温进行热流反演计算方法如下:
热流反演流程如图5所示。假设身部3的内壁面热流在圆周方向均匀分布,则用于热流反演的计算域可简化成二维,如图4所示。时间步长取热电偶11的采样时间间隔,从第一个时间步开始进行热流反演;
第一个时间步的起始时刻(点火时刻)的温度场是已知的,先假设每个热电偶11对应位置处的热流,然后对计算域应用二维非稳态导热方程进行温度场计算,若计算得到的计算温度与热电偶11实测温度有残差,则修改之前假设的热流,直至计算温度与实测温度的温度残差小于设定阈值,即符合收敛标准,这样就获得了第一个时间步的热流分布。
进行下一时间步的热流反演时,以上一时间步结束时刻的温度场作为初始温度场,重复利用上述第一个时间步的热流分布的获取步骤,依次类推即可获得每个时间步的热流分布。
根据所有时间步的热流分布,获得燃烧状态稳定后的燃烧室热流分布特征曲线。
其中,在每一个时间步,不满足收敛标准时,修改假设的热流分布进行迭代的方法如下:
a)设N个热电偶11,对应于N个位置的热流分布q=(q1,q2,...qN),N个热电偶11的实测温度为Tm=(Tm1,Tm2,...TmN),在给定热流分布q条件下的计算温度为Tc(q)=(Tc1(q),Tc2(q),...TcN(q));
温度残差按照下式计算:
J=[Tm-Tc(q)][Tm-Tc(q)]T;
设定阈值按照下式计算:ε=N*(ΔT)2*Δt,其中,ΔT为热电偶的测温精度,Δt为热电偶采样时间间隔。
b)计算温度对热流的响应矩阵S:
c)根据下式计算修改后的热流分布qk+1;
修改前的热流分布qk和修改后的热流分布qk+1满足以下方程:
S·[qk+1]T=[Tm-Tc(qk)]T+S·[qk]T。
以上仅是对本发明的优选实施方式进行了描述,并不将本发明的技术方案限制于此,本领域技术人员在本发明主要技术构思的基础上所作的任何变形都属于本发明所要保护的技术范畴。
Claims (10)
1.一种液体火箭发动机喷注器燃烧状态评估装置,其特征在于:包括连接组件、N个热电偶(11)以及依次同轴设置的前法兰(2)、身部(3)、喷管(4)和后法兰(5),N为大于等于2的整数;身部(3)上设置有燃烧室(31),喷管(4)上设置有与燃烧室(31)同轴且连通的喉道(41);
所述连接组件用于紧固前法兰(2)和后法兰(5),实现身部(3)和喷管(4)的夹紧;
所述身部(3)外壁开设有与燃烧室(31)连通的燃烧室测压孔(9),以及沿径向开设有轴向均布的N个热电偶放置孔(8),热电偶放置孔(8)包括外孔(81)和直径小于外孔(81)的内孔(82),内孔(82)底面与燃烧室(31)内壁面之间存在距离h,h满足:
其中,p为燃烧室(31)工作过程中的最高压力,d为内孔(82)的直径,n为安全系数,σ为燃烧室材料的抗剪切强度;
相邻2个热电偶放置孔(8)的轴线间隔为a,a满足:a≥2d′;
其中,d′为外孔(81)的直径;
所述N个热电偶(11)分别设置在N个热电偶放置孔(8)内,每个热电偶(11)上设有与外孔(81)配合的限位部(18),以及套装有位于限位部(18)外侧的弹簧(16);
所述身部(3)外壁设有热电偶放置孔(8)的位置设置有压板(13),用于通过弹簧(16)将热电偶(11)压紧在热电偶放置孔(8)中。
2.根据权利要求1所述液体火箭发动机喷注器燃烧状态评估装置,其特征在于:所述N个热电偶(11)的外端均伸出压板(13);
所述限位部(18)为设置在热电偶(11)外圆面的环形限位凸起,环形限位凸起设置在外孔(81)上,弹簧(16)的两端分别通过环形限位凸起外端面和压板(13)内端面限位;
所述弹簧(16)上套装有限位套筒(17),限位套筒(17)的一端止靠于身部(3)外壁,另一端与压板(13)内端面之间存在距离。
3.根据权利要求2所述液体火箭发动机喷注器燃烧状态评估装置,其特征在于:所述连接组件包括沿周向布置的M个锁紧件,M为大于等于2的整数,每个锁紧件包括长螺杆(7)和2个第一螺母(12),长螺杆(7)设置在前法兰(2)和后法兰(5)上,且两端分别伸出前法兰(2)和后法兰(5),并与第一螺母(12)连接。
4.根据权利要求3所述液体火箭发动机喷注器燃烧状态评估装置,其特征在于:所述前法兰(2)和后法兰(5)均采用不锈钢材料,身部(3)和喷管(4)采用紫铜材料;
所述后法兰(5)上安装有可拆卸的堵盖(6),用于喷注器(1)点火试验前的密封性检测。
5.根据权利要求1至4任一所述液体火箭发动机喷注器燃烧状态评估装置,其特征在于:所述前法兰(2)上设有火花塞放置孔(10);
所述前法兰(2)与身部(3)之间、身部(3)和喷管(4)之间、喷管(4)和后法兰(5)之间均采用定位槽和定位块配合方式。
6.一种液体火箭发动机喷注器燃烧状态评估方法,其特征在于,包括以下步骤:
1)将待评估喷注器(1)与评估装置的前法兰(2)对接,热电偶(11)安装到位,进行设定时间段的点火试验,测得所有热电偶(11)在各采样时间点的瞬态壁温;
2)对瞬态壁温进行热流反演计算,获得燃烧状态稳定后的燃烧室热流分布特征曲线;
3)根据燃烧室热流分布特征曲线,确定热流分布特征曲线的平稳段起点,该平稳段起点为火焰结束位置,根据火焰结束位置获得火焰长度,进而评估喷注器(1)燃烧状态。
7.根据权利要求6所述液体火箭发动机喷注器燃烧状态评估方法,其特征在于,步骤2)具体如下:
2.1)假设身部(3)的内壁面热流在圆周方向均匀分布,将用于热流反演的计算域简化为二维,时间步长取热电偶(11)的采样时间间隔;
2.2)第一个时间步的起始时刻记为点火时刻,其初始温度场已知;第一个时间步每个热电偶(11)对应位置处的热流分别给定一个假设值,然后对计算域应用二维非稳态导热方程进行温度场计算,得到每个热电偶(11)的计算温度,根据二维非稳态导热方程计算得到计算温度,进而得到计算温度与热电偶(11)实测瞬态壁温的温度残差,比较温度残差与设定阈值,若温度残差小于设定阈值,获得第一时间步的热流分布;若否,修改每个热电偶(11)对应位置处热流的假设值进行迭代,直至温度残差小于设定阈值;
2.3)以第一时间步结束时刻的温度场为第二时间步的初始温度场,采用步骤2.2)的方法获得第二时间步的热流分布;
2.4)利用步骤2.3)的方法,获得剩余每个时间步的热流分布;
2.5)根据所有时间步的热流分布,获得燃烧状态稳定后的燃烧室热流分布特征曲线。
8.根据权利要求7所述液体火箭发动机喷注器燃烧状态评估方法,其特征在于,步骤2.2)中;
计算温度残差具体为:设N个热电偶(11)对应于N个位置的热流分布q=(q1,q2,...qN),N个热电偶(11)的实测温度为Tm=(Tm1,Tm2,...TmN),在给定热流分布q条件下的计算温度为Tc(q)=(Tc1(q),Tc2(q),...TcN(q));
温度残差按照下式计算:
J=[Tm-Tc(q)][Tm-Tc(q)]T;
设定阈值按照下式计算:
ε=N*(ΔT)2*Δt;
其中,ΔT为热电偶的测温精度,Δt为热电偶采样时间间隔;
修改每个热电偶(11)对应位置处热流的假设值进行迭代具体为:根据下式计算修改后的热流分布qk+1;
S·[qk+1]T=[Tm-Tc(qk)]T+S·[qk]T;
式中,qk为修改前的热流分布;S为温度对热流的响应矩阵。
9.根据权利要求8所述液体火箭发动机喷注器燃烧状态评估方法,其特征在于:步骤2.1)中,所述设定时间段为2s~5s。
10.根据权利要求9所述液体火箭发动机喷注器燃烧状态评估方法,其特征在于:步骤3)中,根据火焰结束位置获得火焰长度,进而评估喷注器(1)燃烧状态具体为:火焰结束位置越靠近上游,火焰长度越短,则喷注器(1)燃烧越高效。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111061147.8A CN113790899B (zh) | 2021-09-10 | 2021-09-10 | 液体火箭发动机喷注器燃烧状态评估装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111061147.8A CN113790899B (zh) | 2021-09-10 | 2021-09-10 | 液体火箭发动机喷注器燃烧状态评估装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113790899A true CN113790899A (zh) | 2021-12-14 |
CN113790899B CN113790899B (zh) | 2024-02-09 |
Family
ID=79183000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111061147.8A Active CN113790899B (zh) | 2021-09-10 | 2021-09-10 | 液体火箭发动机喷注器燃烧状态评估装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113790899B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114876674A (zh) * | 2022-07-08 | 2022-08-09 | 太原理工大学 | 一种实验液体火箭发动机用水冷喷管 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030018224A (ko) * | 2001-08-27 | 2003-03-06 | 현대모비스 주식회사 | 액체로켓 인젝터시험장치 |
CN101608953A (zh) * | 2008-06-19 | 2009-12-23 | 北京航空航天大学 | 一种燃烧室内壁面温度和热流分布的测量方法及装置 |
CN102022225A (zh) * | 2010-12-30 | 2011-04-20 | 北京航空航天大学 | 富氢/富氧燃气多喷嘴气气喷注器试验装置 |
CN106370433A (zh) * | 2016-08-12 | 2017-02-01 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种评估发动机主燃烧室火焰筒出口温度场品质的方法 |
CN111879443A (zh) * | 2020-07-14 | 2020-11-03 | 西安近代化学研究所 | 火箭发动机内气-液两相热流密度测量工装 |
-
2021
- 2021-09-10 CN CN202111061147.8A patent/CN113790899B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030018224A (ko) * | 2001-08-27 | 2003-03-06 | 현대모비스 주식회사 | 액체로켓 인젝터시험장치 |
CN101608953A (zh) * | 2008-06-19 | 2009-12-23 | 北京航空航天大学 | 一种燃烧室内壁面温度和热流分布的测量方法及装置 |
CN102022225A (zh) * | 2010-12-30 | 2011-04-20 | 北京航空航天大学 | 富氢/富氧燃气多喷嘴气气喷注器试验装置 |
CN106370433A (zh) * | 2016-08-12 | 2017-02-01 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种评估发动机主燃烧室火焰筒出口温度场品质的方法 |
CN111879443A (zh) * | 2020-07-14 | 2020-11-03 | 西安近代化学研究所 | 火箭发动机内气-液两相热流密度测量工装 |
Non-Patent Citations (2)
Title |
---|
刘占一;刘计武;汪广旭;石晓波;HAIDN O J;: "气氧甲烷单喷嘴燃烧室壁面热流的测量和数值模拟", 推进技术, no. 05 * |
孙冰;刘迪;王太平;: "液体火箭发动机燃烧室壁面热流测量方法研究", 推进技术, no. 09 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114876674A (zh) * | 2022-07-08 | 2022-08-09 | 太原理工大学 | 一种实验液体火箭发动机用水冷喷管 |
CN114876674B (zh) * | 2022-07-08 | 2022-11-15 | 太原理工大学 | 一种实验液体火箭发动机用水冷喷管 |
Also Published As
Publication number | Publication date |
---|---|
CN113790899B (zh) | 2024-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101608953B (zh) | 一种燃烧室内壁面温度和热流分布的测量方法及装置 | |
CN105547578A (zh) | 一种带冷却水套的高温燃气动态压力引压装置 | |
CN109900478B (zh) | 一种中小型无人机火箭推力测试装置及测试方法 | |
US7819578B2 (en) | Fluid temperature measurement device | |
CN113790899A (zh) | 液体火箭发动机喷注器燃烧状态评估装置及方法 | |
CN110186583B (zh) | 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法 | |
US4843887A (en) | Method and device for measuring pressures in tubular bodies | |
CN105675643A (zh) | 复合固体推进剂烤燃实验装置 | |
CN104457457A (zh) | 雷管延期体延时测试新方法 | |
US7962307B2 (en) | Integrated apparatus for measuring static pressure | |
CN113125025B (zh) | 一种真空热场测温铠装热电偶拆卸装置及测温方法 | |
CN111795762B (zh) | 一种热流计动态特性校准装置及方法 | |
CN110146550B (zh) | 基于电阻抗成像的复合材料高温部件氧化程度监测方法 | |
CN112943484A (zh) | 一种研究粗糙度对喷管壁面传热影响的实验装置 | |
CN111289561A (zh) | 一种气凝胶隔热工装结构及隔热板材耐热性考核方法 | |
CN116519170A (zh) | 一种适用于爆炸环境下的热电偶动态校准装置 | |
CN104977158A (zh) | 一种管内燃烧超高温在线测压液动高温三通球阀试验装置 | |
Venkataraman et al. | Thin-wire thermocouple design for exhaust gas temperature pulse measurements in internal combustion engines | |
CN2932341Y (zh) | 裂解炉cot热电偶抗冲刷套管 | |
CN114837854B (zh) | 一种模拟发动机慢速烤燃的低成本试验装置及方法 | |
CN114776483B (zh) | 三维瞬态燃速测量装置及燃速测量方法 | |
US4595824A (en) | Heater block assembly for use in thermal oxidation testing of jet fuel | |
CN110208386A (zh) | 高温内耗值测试方法及装置 | |
CN109443681A (zh) | 一种适用于高真空和高温高压条件下的测试接头 | |
CN116562193B (zh) | 旋转爆震发动机燃烧效率分析方法及系统 |
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 |