CN109408915A - 固体火箭超燃冲压发动机燃烧流场仿真方法 - Google Patents
固体火箭超燃冲压发动机燃烧流场仿真方法 Download PDFInfo
- Publication number
- CN109408915A CN109408915A CN201811181262.7A CN201811181262A CN109408915A CN 109408915 A CN109408915 A CN 109408915A CN 201811181262 A CN201811181262 A CN 201811181262A CN 109408915 A CN109408915 A CN 109408915A
- Authority
- CN
- China
- Prior art keywords
- combustion
- gas
- flow field
- model
- fuel
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Engines (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明涉及一种固体火箭超燃冲压发动机燃烧流场仿真方法,属于发动机技术领域。本发明提供的一种固体火箭冲压发动机燃烧仿真方法,考虑二次燃烧化学平衡,通过二次燃烧热力学计算,得到生成物标准状态摩尔生成焓,确保燃料燃烧释放的能量与实际情况更为接近,能够更加准确的模拟燃烧温度,保证仿真温度的模拟更为准确。
Description
技术领域
本发明属于发动机技术领域,具体涉及一种固体火箭超燃冲压发动机燃烧流场仿真方法。
背景技术
数值仿真是发动机燃烧流动研究的重要手段,固体火箭超燃冲压发动机燃烧流动仿真十分复杂,燃气发生器内富燃料推进剂燃烧产生一次燃气,一次燃气喷入补燃室后,一次燃气中的可燃成分与来流空气发生燃烧化学反应。
对于固体火箭超燃冲压发动机燃烧流场仿真,难点主要在于模拟一次燃气燃烧的质量流率和成分,以及模拟在高速气流中的燃烧过程。
现有的固体燃料冲压发动机燃烧流场仿真方法见表1,其采用商用软件自带模型,其一次燃烧分解产物采用简化单一物质,二次燃烧化学反应物未考虑化学平衡,而是直接反应成最终产物,热值全部释放,而实际的燃烧过程存在一个化学平衡,热值不会全部释放,因此按照上述方法计算的燃烧温度将远远高于实际燃烧温度,燃烧温度对于超声速燃烧的模态判断及流场分析影响重大。
表1现有仿真方法简表
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何提供一种固体火箭冲压发动机燃烧仿真方法,提高固体火箭超燃冲压发动机燃烧仿真精度。
(二)技术方案
为了解决上述技术问题,本发明提供了一种固体火箭超燃冲压发动机燃烧流场仿真方法,该方法包括以下步骤:
步骤1、计算一次燃气成分,即富燃燃气成分
采用最小吉布斯自由能法,进行富燃料推进剂一次燃烧热力学计算,通过推进剂组分计算得到一次燃气总质量、各项成分及质量分数;提取质量分数大于1%的可燃物质,其余物质均视为非可燃物质;
步骤2、设置二次燃烧物性参数;
采用最小吉布斯自由能法,进行富燃料推进剂二次燃烧热力学计算,得到二次燃烧理论燃烧温度;根据步骤1得到的可燃物质,选择总包反应,确定燃烧反应生成物,根据能量方程计算生成物的标准状态摩尔生成焓,使燃料完全燃烧对应的温度等于理论燃烧温度,更接近实际流场,修改标准状态摩尔生成焓,并按照物质的原有属性参数填写其他二次燃烧物性参数;
步骤3、根据步骤1的计算结果,以及计算任务需要设置进口和出口计算的边界条件,并对流场的各参数进行初始化;
步骤4、建立湍流模型,在FLUENT软件中选取该湍流模型;
步骤5、建立气相燃烧模型,计算气相化学反应速率;
步骤6、通过步骤2设置的物质属性,以及步骤4、步骤5的模型,求解气相控制方程。
优选地,步骤4中所述湍流模型采用Realizablek-ε模型。
优选地,步骤5中所述气相燃烧模型为涡团耗散模型。
优选地,步骤6具体为:求解FLUENT软件自带的雷诺N-S方程组,对流体的速度、温度、密度进行求解;如果方程组收敛,则所有流场参数为最终结果,如果方程组不收敛,对所有初始参数加松弛因子后重新返回步骤1,直至方程收敛。
(三)有益效果
本发明提供的一种固体火箭冲压发动机燃烧仿真方法,考虑二次燃烧化学平衡,通过二次燃烧热力学计算,得到生成物标准状态摩尔生成焓,确保燃料燃烧释放的能量与实际情况更为接近,能够更加准确的模拟燃烧温度,保证仿真温度的模拟更为准确。
附图说明
图1是本发明的仿真方法流程图。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
本发明提供了一种固体火箭超燃冲压发动机燃烧流场仿真方法,如图1所示,该方法包括以下步骤:
步骤1、计算一次燃气成分(富燃燃气成分)
采用最小吉布斯自由能法,进行富燃料推进剂一次燃烧热力学计算,通过推进剂组分计算得到一次燃气总质量、各项成分及质量分数;提取质量分数大于1%的可燃物质,其余物质均视为非可燃物质,命名为G。
步骤2、设置二次燃烧物性参数;
采用最小吉布斯自由能法,进行富燃料推进剂二次燃烧热力学计算,得到二次燃烧理论燃烧温度;根据步骤1得到的可燃物质,选择总包反应,确定燃烧反应生成物,根据能量方程计算生成物的标准状态摩尔生成焓,使燃料完全反应(燃烧)对应的温度等于理论燃烧温度,更接近实际流场。除标准状态摩尔生成焓需要修改之外,其他二次燃烧物性参数均按照物质的原有属性参数填写。
步骤3、根据步骤1的计算结果,以及计算任务需要设置进口和出口计算的边界条件,并对流场的各参数进行初始化;
步骤4、建立湍流模型,包括:
采用Realizablek-ε模型,在FLUENT软件中选取该模型;
步骤5、建立气相燃烧模型,包括:
采用涡团耗散模型,计算气相化学反应速率;
步骤6、通过步骤2设置的物质属性,以及步骤4到5的模型,求解气相控制方程,包括:
求解FLUENT自带的雷诺N-S方程组,对流体的速度、温度、密度进行求解;如果方程组收敛,则所有流场参数为最终结果,如果方程组不收敛,对所有初始参数加松弛因子后重新进行第1步的计算,重复步骤全过程,循环往复上述步骤,直至方程收敛。
超音速燃烧对仿真温度的计算精确度要求很高,本发明提供的一种固体火箭冲压发动机燃烧仿真方法,考虑二次燃烧化学平衡,通过二次燃烧热力学计算,得到生成物标准状态摩尔生成焓,确保燃料燃烧释放的能量与实际情况更为接近,能够更加准确的模拟燃烧温度,保证仿真温度的模拟更为准确。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (4)
1.一种固体火箭超燃冲压发动机燃烧流场仿真方法,其特征在于,该方法包括以下步骤:
步骤1、计算一次燃气成分,即富燃燃气成分
采用最小吉布斯自由能法,进行富燃料推进剂一次燃烧热力学计算,通过推进剂组分计算得到一次燃气总质量、各项成分及质量分数;提取质量分数大于1%的可燃物质,其余物质均视为非可燃物质;
步骤2、设置二次燃烧物性参数;
采用最小吉布斯自由能法,进行富燃料推进剂二次燃烧热力学计算,得到二次燃烧理论燃烧温度;根据步骤1得到的可燃物质,选择总包反应,确定燃烧反应生成物,根据能量方程计算生成物的标准状态摩尔生成焓,使燃料完全燃烧对应的温度等于理论燃烧温度,更接近实际流场,修改标准状态摩尔生成焓,并按照物质的原有属性参数填写其他二次燃烧物性参数;
步骤3、根据步骤1的计算结果,以及计算任务需要设置进口和出口计算的边界条件,并对流场的各参数进行初始化;
步骤4、建立湍流模型,在FLUENT软件中选取该湍流模型;
步骤5、建立气相燃烧模型,计算气相化学反应速率;
步骤6、通过步骤2设置的物质属性,以及步骤4、步骤5的模型,求解气相控制方程。
2.如权利要求1所述的方法,其特征在于,步骤4中所述湍流模型采用Realizablek-ε模型。
3.如权利要求2所述的方法,其特征在于,步骤5中所述气相燃烧模型为涡团耗散模型。
4.如权利要求3所述的方法,其特征在于,步骤6具体为:求解FLUENT软件自带的雷诺N-S方程组,对流体的速度、温度、密度进行求解;如果方程组收敛,则所有流场参数为最终结果,如果方程组不收敛,对所有初始参数加松弛因子后重新返回步骤1,直至方程收敛。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811181262.7A CN109408915B (zh) | 2018-10-11 | 2018-10-11 | 固体火箭超燃冲压发动机燃烧流场仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811181262.7A CN109408915B (zh) | 2018-10-11 | 2018-10-11 | 固体火箭超燃冲压发动机燃烧流场仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109408915A true CN109408915A (zh) | 2019-03-01 |
CN109408915B CN109408915B (zh) | 2022-10-14 |
Family
ID=65467516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811181262.7A Active CN109408915B (zh) | 2018-10-11 | 2018-10-11 | 固体火箭超燃冲压发动机燃烧流场仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109408915B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110727906A (zh) * | 2019-09-02 | 2020-01-24 | 中国地质大学(武汉) | 基于多组分热解气体的固体可燃物燃烧特性求解方法 |
CN111079235A (zh) * | 2019-12-11 | 2020-04-28 | 内蒙动力机械研究所 | 一种固体火箭发动机内流场仿真快速收敛方法 |
CN112417775A (zh) * | 2020-10-30 | 2021-02-26 | 中国运载火箭技术研究院 | 固体火箭发动机推力终止泄压过程内弹道性能计算方法 |
CN115618171A (zh) * | 2022-06-06 | 2023-01-17 | 北京理工大学 | 一种基于同伦算法的推进剂燃烧平衡产物求解方法 |
CN113836647B (zh) * | 2021-07-16 | 2024-04-26 | 西安流固动力科技有限公司 | 弹射内弹道仿真分析系统、方法及电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6944583B1 (en) * | 2000-02-02 | 2005-09-13 | Centric Software, Inc. | Multi-threaded frame safe synchronization of a simulation |
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真系统 |
CN105354401A (zh) * | 2015-12-24 | 2016-02-24 | 中国人民解放军装备学院 | 一种多喷管火箭或导弹尾焰流场计算方法 |
CN106295082A (zh) * | 2016-09-27 | 2017-01-04 | 华中科技大学 | 一种平板式固体氧化物燃料电池的数值模拟方法 |
CN107346353A (zh) * | 2017-06-05 | 2017-11-14 | 民政部零研究所 | 一种固体燃烧物燃烧过程仿真方法及服务器 |
CN107391436A (zh) * | 2017-06-20 | 2017-11-24 | 北京航空航天大学 | 基于火焰面/进度变量模型的超声速湍流燃烧流动计算方法 |
-
2018
- 2018-10-11 CN CN201811181262.7A patent/CN109408915B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6944583B1 (en) * | 2000-02-02 | 2005-09-13 | Centric Software, Inc. | Multi-threaded frame safe synchronization of a simulation |
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真系统 |
CN105354401A (zh) * | 2015-12-24 | 2016-02-24 | 中国人民解放军装备学院 | 一种多喷管火箭或导弹尾焰流场计算方法 |
CN106295082A (zh) * | 2016-09-27 | 2017-01-04 | 华中科技大学 | 一种平板式固体氧化物燃料电池的数值模拟方法 |
CN107346353A (zh) * | 2017-06-05 | 2017-11-14 | 民政部零研究所 | 一种固体燃烧物燃烧过程仿真方法及服务器 |
CN107391436A (zh) * | 2017-06-20 | 2017-11-24 | 北京航空航天大学 | 基于火焰面/进度变量模型的超声速湍流燃烧流动计算方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110727906A (zh) * | 2019-09-02 | 2020-01-24 | 中国地质大学(武汉) | 基于多组分热解气体的固体可燃物燃烧特性求解方法 |
CN110727906B (zh) * | 2019-09-02 | 2023-08-29 | 中国地质大学(武汉) | 基于多组分热解气体的固体可燃物燃烧特性求解方法 |
CN111079235A (zh) * | 2019-12-11 | 2020-04-28 | 内蒙动力机械研究所 | 一种固体火箭发动机内流场仿真快速收敛方法 |
CN112417775A (zh) * | 2020-10-30 | 2021-02-26 | 中国运载火箭技术研究院 | 固体火箭发动机推力终止泄压过程内弹道性能计算方法 |
CN112417775B (zh) * | 2020-10-30 | 2024-03-26 | 中国运载火箭技术研究院 | 固体火箭发动机推力终止泄压过程内弹道性能计算方法 |
CN113836647B (zh) * | 2021-07-16 | 2024-04-26 | 西安流固动力科技有限公司 | 弹射内弹道仿真分析系统、方法及电子设备 |
CN115618171A (zh) * | 2022-06-06 | 2023-01-17 | 北京理工大学 | 一种基于同伦算法的推进剂燃烧平衡产物求解方法 |
CN115618171B (zh) * | 2022-06-06 | 2023-10-24 | 北京理工大学 | 一种基于同伦算法的推进剂燃烧平衡产物求解方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109408915B (zh) | 2022-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109408915A (zh) | 固体火箭超燃冲压发动机燃烧流场仿真方法 | |
Xisto et al. | The efficiency of a pulsed detonation combustor–axial turbine integration | |
CN105354401B (zh) | 一种多喷管火箭或导弹尾焰流场计算方法 | |
Badami et al. | Experimental and numerical analysis of a small-scale turbojet engine | |
Seldner | Generalized simulation technique for turbojet engine system analysis | |
Bottone et al. | The numerical simulation of diesel spray combustion with LES-CMC | |
CN109657401A (zh) | 一种固体燃料冲压发动机燃烧流场数值仿真方法 | |
Bedarev et al. | Simulating the regimes of oblique detonation waves arising at detonation initiation by a small-diameter projectile | |
CN110348078A (zh) | 一种涡轴发动机容积动力学结合热惯性效应的建模方法 | |
CN110442934A (zh) | 一种考虑固体发动机尾喷流辐射的高精度气动热计算方法 | |
CN112562793B (zh) | 一种针对燃料爆震燃烧的两步反应模型计算方法 | |
CN107247823A (zh) | 基于准一维化学动力学过程和自模性的喷焰流场预估方法 | |
Völzke et al. | The influence of the initial temperature on DDT characteristics in a valveless PDC | |
Zubrilin et al. | Simulation of CO emission in primary and secondary zone of a small gas turbine combustion chamber using CFD and reactors network | |
Cuciumita et al. | Evaluation of the Cycle averaged performances of a pulsed detonation engine based on thermodynamic cycle computations | |
Ponizy et al. | On modeling of pulse combustors | |
CN110530642A (zh) | 一种燃烧器试验台架及基于燃烧器的混合器稳态结晶试验方法 | |
Liu et al. | Numerical and experimental investigation on emission performance of a fuel staged combustor | |
CN115221638A (zh) | 亚燃冲压发动机非稳态过程的性能时间响应分析方法 | |
Irannezhad | A numerical study of reacting flows using finite rate chemistry | |
Kennedy et al. | Prediction of buoyancy controlled turbulent wall diffusion flames | |
Yentsch et al. | Exploratory simulations of the HIFiRE 2 scramjet flowpath | |
Boyde et al. | Validation of an ignition and flame propagation model for multiphase flows | |
Juniper et al. | The effect of damköhler number on the stand-off distance of cross-flow flames | |
Nusca et al. | Modeling Hypergolic Ignition in the Army's Impinging Stream Vortex Engine Including Injection Throttling |
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 |