CN107247823B - 基于准一维化学动力学过程和自模性的喷焰流场预估方法 - Google Patents
基于准一维化学动力学过程和自模性的喷焰流场预估方法 Download PDFInfo
- Publication number
- CN107247823B CN107247823B CN201710342079.XA CN201710342079A CN107247823B CN 107247823 B CN107247823 B CN 107247823B CN 201710342079 A CN201710342079 A CN 201710342079A CN 107247823 B CN107247823 B CN 107247823B
- Authority
- CN
- China
- Prior art keywords
- flame
- flow field
- jet
- low
- template
- 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
Images
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)
Abstract
本发明公开了基于准一维化学动力学过程和自模性的喷焰流场预估方法,本发明针对低空喷焰的流场特点,可以划分为初始区、过渡区和核心区三个部分,在不同的喷焰流场分区,采用不同的计算方法实现对低空喷焰的快速仿真。本发明能够实现低空喷焰初始段、过渡段、核心区高精度的带化学动力反应过程的流场实时/准实时预估,为低空喷焰流场特性的快速准确计算提供一种方法,解决以往基于CFD方法喷焰流场计算效率低的问题。
Description
技术领域
本发明涉及基于准一维化学动力学过程和自模性的低空喷焰流场预估方法,特别是用于低空喷焰红外辐射特性仿真的流场快速计算方法。
背景技术
目前CFD(计算流体力学)方法是实现低空喷焰流场仿真的主要方法。对于动目标的低空喷焰流场仿真而言,CFD方法是针对剖分的网格求解考虑化学反应动力过程、输运方程、能量守恒的NS方程,为了保证收敛性和求解精度,一般网格数量至少在50-200万的量级,仿真获得某一典型状态下的流场计算时间成本较高,效率较低,远不能满足实际工程的需求。
发明内容
针对现有技术中存在的不足,本发明提供了基于准一维化学动力学过程和自模性的喷焰流场预估方法。
本发明是通过如下技术方案实现的:基于准一维化学动力学过程和自模性的喷焰流场预估方法,包括如下步骤:步骤1)低空喷焰流场初始段预估:
基于膨胀比、伴随流速度、喷焰出口的比热比三个核心参数建立大量的初始段模板,并以比热比、膨胀比、伴随流速度依次从模板中选取最接近的模板作为基础,构建任意状态下的低空喷焰初始段流场;
步骤2)低空喷焰流场核心区预估:
首先利用准一维气体动力学守恒定律,采用Arrhenius化学反应模型,计算轴线方向上化学动力过程的喷焰流场特性;然后利用自模性完成喷焰流场特性的二维扩展,从而构建低空喷焰的核心区流场。
上述技术方案中,所述低空喷焰流场初始段预估包括以下步骤:
第一步:建立典型状态的喷焰初始段模板,建立随膨胀比NPR、喷管出口的比热比γg、伴随流速度U∞变化的初始段流场的模板;
基于CFD++建立化学冻结过程,包含喷焰气体和环境气体的喷焰初始段模板,
第二步:喷焰初始段模板选择
对于基于准一维化学动力学过程和自模性的低空喷焰流场快速预估方法,喷焰初始段模板基于以下顺序选择相应的模型,首先匹配最接近的γg,接下去是NPR,最后是U∞,基于选择的模板,计算流场特性;
第三步:控制喷焰初始段流场尺度
得到高速伴随流对喷焰尺度Lt的影响函数g(U∞)。
上述技术方案中,初始段的长度Lt是从喷管出口到环境气体的质量分数减小到0.1 时的距离,喷焰的最大等效半径Rt是喷焰边界的温度降低到环境来流的1.05倍。
上述技术方案中,按下述公式计算流场特性:
上述技术方案中,根据给定的以下三式:
式1:式1:高速伴随流影响Lt的函数f(U∞)及高速伴随流影响Rt的函数g(U∞)”,
g(U∞)=1/U∞;
式2至式3:初始段喷焰流场的尺度控制如下所示:
其中,kL和kg是拟合常数,得到高速伴随流下喷焰的尺度L和R。
上述技术方案中,低空喷焰流场核心区预估包括:
第一步:准一维模型求解喷焰流场核心区:
采用准一维模型求解核心区轴线上的喷焰流场平均特性,通过求解下述方程,其中起始的长度单位以初始段的半径b0,
基于上述的假设,下面给出了从喷管出口到喷焰核心区的守恒关系式如下:
组分:Y0i=Yei
其中ε是喷管的扩张比,CFmax和CD是最大的推力系数和阻力系数;
第二步:自模性喷焰流场核心区扩展
通过自模性方法计算流场径向方向的流体参数,
对于流体特性如轴向速度、压强、温度和质量浓度在轴向任意位置上的径向分布近似满足高斯分布,如下式:
其中rφ是流场特性的半宽尺寸,其表征的中心线上值1/2处的流场宽度,σφ表征的是状态背离参数,
通过上述第一步的准一维方程的求解,获得任意切面处的流动属性参数和射流的半径b,假定自模性的流动参数趋近于伴随流的混合边界,即在r=b/2处存在φ=φ∞,喷焰气体密度在任意切面处是恒定的,考虑流动参数的守恒性,修正后的φm和如下式所示流动属性参数可认为rφ=b/6和σφ为0.5时,采用下述公式可求解全部的流场属性:
通过上述步骤,可实现低空喷焰流场的快速预估,给出喷焰流场的压强、温度、密度、组分浓度的分布特性,支撑喷焰红外辐射特性的快速仿真。
上述技术方案中,卷吸系数α是随马赫数变化的参数,其变化符合如下参数;
本发明具有如下有益效果:本发明能够实现低空喷焰初始段、过渡段、核心区高精度的带化学动力反应过程的流场实时/准实时预估,为低空喷焰流场特性的快速准确计算提供一种方法,解决以往基于CFD方法喷焰流场计算效率低的问题。
附图说明
图1是低空喷焰流场的结构分布图
图2是低空喷焰流场预估方法的流程图
图3是低空喷焰流场初始段模板图
图4是当γg=1.2时,Lt随NPR和U∞拟合曲线
图5是本方法与CFD方法仿真的轴线上温度与组分分布图
图6是本方法与CFD方法仿真的温度与组分分布云图
图7是不同高度下本方法与CFD方法仿真的轴线上温度与组分分布图
图8是不同高度下本方法与CFD方法仿真的温度与组分分布云图
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:参见图1至图4,基于准一维化学动力学过程和自模性的喷焰流场预估方法,包括如下步骤。
首先,获取低空喷焰的流体参数,包括发动机出口气体和环境气体的参数。
步骤1)进行低空喷焰流场初始段预估,具体如下。
第一步:建立典型状态的喷焰初始段模板
喷焰流场是由喷管出口流动参数和伴随流(即环境大气来流)决定的,可归结为膨胀比NPR(喷管出口压强与环境压强比)、喷管出口的比热比γg、伴随流速度U∞。由于喷管出口速度很高(一般大于1500m/s),外部环境中气体无法进入喷焰内部,所以初始段可以忽略化学反应的作用。建立随膨胀比NPR、喷管出口的比热比γg、伴随流速度U∞变化的初始段流场的模板。
基于CFD++建立化学冻结过程,包含喷焰气体和环境气体的喷焰初始段模板,初始段的长度Lt是从喷管出口到环境气体的质量分数减小到0.1时的距离,喷焰的最大等效半径Rt是喷焰边界的温度降低到环境来流的1.05倍。
第二步:喷焰初始段模板选择
对于基于准一维化学动力学过程和自模性的低空喷焰流场快速预估方法,喷焰初始段模板基于以下顺序选择相应的模型,首先匹配最接近的γg,接下去是NPR,最后是U∞,基于选择的模板,按下述公式计算流场特性:
第三步:控制喷焰初始段流场尺度
因为建立的模板无法覆盖所有的流动状态,基于典型模板统计分析,对于给定的γg和 U∞<U*/2(其中),喷焰尺度Lt和Rt随Ln(NPR)线性变化,随着伴随流速度U∞的增加,对应的Lt随之增加,当伴随流速度趋近于喷焰极限膨胀速度时,Lt达到其最大值,后随着伴随流增加,Lt减小。见图4。
下式给出了高速伴随流影响Lt的函数f(U∞)及高速伴随流影响Rt的函数g(U∞):
g(U∞)=1/U∞;
初始段喷焰流场的尺度控制如下所示:
其中,kL和kg是拟合常数,得到高速伴随流下喷焰的尺度L和R。
步骤2)进行低空喷焰流场核心区预估,具体如下。
第一步:准一维模型求解喷焰流场核心区
采用准一维模型求解核心区轴线上的喷焰流场平均特性。通过求解下述方程,其中起始的长度单位以初始段的半径b0。
基于上述的假设,下面给出了从喷管出口到喷焰核心区的守恒关系式如下:
组分:Y0i=Yei
其中ε是喷管的扩张比,CFmax和CD是最大的推力系数和阻力系数。
第二步:自模性喷焰流场核心区扩展
通过自模性方法计算流场径向方向的流体参数。
对于流体特性如轴向速度、压强、温度和质量浓度在轴向任意位置上的径向分布近似满足高斯分布,如下式:
其中rφ是流场特性的半宽尺寸,其表征的中心线上值1/2处的流场宽度,σφ表征的是状态背离参数。
通过上述第一步的准一维方程的求解,获得任意切面处的流动属性参数和射流的半径b。假定自模性的流动参数趋近于伴随流的混合边界,即在r=b/2处存在φ=φ∞,喷焰气体密度在任意切面处是恒定的,考虑流动参数的守恒性,修正后的φm和如下式所示
流动属性参数可认为rφ=b/6和σφ为0.5时,采用下述公式可求解全部的流场属性:
通过上述步骤,可实现低空喷焰流场的快速预估,给出喷焰流场的压强、温度、密度、组分浓度的分布特性,支撑喷焰红外辐射特性的快速仿真。
以下对本发明的实施例中的各个参数符号进行说明。
符号说明:
b为喷焰区域的半径,单位:m
cp定压比热,单位:J/(kg·K)
D喷管出口的直径,单位:m
h焓,单位:J
L初始段和过渡段的轴线长度,单位:m
M喷焰流场的马赫数
NPR为膨胀比,表征为喷管出口压强与环境压强比
P流体的压强,单位:Pa
R初始段和过渡段径向长度,单位:m;
气体常数r径向长度,单位:m
x轴向长度,单位:m
T喷焰流场温度,单位:K
U为伴随流速度,单位:m/s
u为喷焰流场的轴向速度,单位:m/s
v为喷焰流场的径向速度,单位:m/s
W为喷焰流场的摩尔质量,单位:kg/(kmol)
Y为质量分数
ρ为喷焰流场的密度,单位:kg/m3
说明书中各公式所用到的下标含义:
0喷焰初始参数 a环境参数
g喷焰气体参数 c燃烧室参数
e喷管出口参数 i喷焰组分的计数
j喷焰化学反应的计数 m轴线上喷焰参数
t模板参数 ∞环境来流条件
示例说明如下:
U∞为伴随流速度。
Lt为初始段的长度,表征为从喷管出口到环境气体的质量分数减小到0.1时的距离。
Rt为喷焰的最大等效半径,表征为喷焰边界的温度降低到环境来流的1.05倍时的长度。
Yig=YigeYgt其中Yig为第i组分的气体质量分数,Yige为模板中第i组分喷管出口处气体质量分数,Ygt为模板中第i组分的气体质量分数。
Yia=YiaYat其中Yia为环境来流第i组分气体质量分数,Yia为模板中环境来流第i组分气体质量分数,Yat为模板中环境来流第i组分的气体质量分数。
P=ρRT其中P为压强;ρ为密度;R为气体常数;T为温度。
本发明针对低空喷焰的流场特点,可以划分为初始区、过渡区和核心区三个部分,在不同的喷焰流场分区,采用不同的计算方法实现对低空喷焰的快速仿真。其中初始区的流场结构主要由发动机出口参数和环境压力所决定的,但是不同的导弹目标随高度喷焰流场初始段是结构相似的,采用模板匹配方法实现初始段和过渡段的快速预估,在初始段和过渡段由于喷焰速度较快,外部的空气难以进入喷焰内部,此段不存在复燃,因此分三步实现:①模板选择,根据发动机出口参数选择相对应的初始段流场结构;②设置流场参数性质;③流场尺度的标度律修正;核心区由于喷焰速度降低,环境中空气被卷吸进入喷焰中,与未燃的组分(H2、CO等)发现化学反应,进而产生高温复燃区,采用准一维化学反应模型和自模性扩展来估算喷焰的核心区。基于上述的方法,从而实现二维喷焰流场的快速计算,在快速计算过程中仅需求解一维的流体动力学过程,这相对于传统的求解方法,计算效率会有极大地提升,同时本方法也考虑的复燃效应,与传统仿真方法不会存在太大的误差。
本发明的主要创新点是针对喷焰流场的结构,构建了不同仿真方法,并基于准一维化学反应模型和自模性扩展的方法实现了核心区高效、快速仿真,极大地提升了低空喷焰流场的仿真效率。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。
Claims (3)
1.一种基于准一维化学动力学过程和自模性的喷焰流场预估方法,包括如下步骤:
进行低空喷焰流场初始段预估,具体包括:
基于膨胀比、伴随流速度、喷焰出口的比热比三个核心参数建立大量的初始段模板,并以比热比、膨胀比、伴随流速度依次从模板中选取最接近的模板作为基础,构建各个状态下的低空喷焰初始段流场;
进行低空喷焰流场核心区预估,包括:
利用准一维气体动力学守恒定律,采用Arrhenius化学反应模型,计算轴线方向上化学动力过程的喷焰流场特性;
利用自模性完成喷焰流场特性的二维扩展,构建低空喷焰的核心区流场。
2.如权利要求1所述的喷焰流场预估方法,其特征在于:所述低空喷焰流场初始段预估具体为以下步骤:
第一步:建立典型状态的喷焰初始段模板,建立随膨胀比NPR、喷管出口的比热比γg、伴随流速度U∞变化的初始段流场的模板,基于CFD++建立化学冻结过程,包含喷焰气体和环境气体的喷焰初始段模板;
第二步:喷焰初始段模板选择:
对于基于准一维化学动力学过程和自模性的低空喷焰流场快速预估方法,喷焰初始段模板基于以下顺序选择相应的模型,首先匹配最接近的γg,接下去是NPR,最后是U∞,基于选择的模板,计算流场特性;其中,流程特性包括流场温度、压强、密度、速度、组分质量和浓度分布;
第三步:控制喷焰初始段流场尺度:
得到高速伴随流对喷焰尺度Lt的影响函数g(U∞)。
3.如权利要求2所述的喷焰流场预估方法,其特征在于:初始段的长度Lt是从喷管出口到环境气体的质量分数减小到0.1时的距离,喷焰的最大等效半径Rt是喷焰边界的温度降低到环境来流的1.05倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710342079.XA CN107247823B (zh) | 2017-05-12 | 2017-05-12 | 基于准一维化学动力学过程和自模性的喷焰流场预估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710342079.XA CN107247823B (zh) | 2017-05-12 | 2017-05-12 | 基于准一维化学动力学过程和自模性的喷焰流场预估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107247823A CN107247823A (zh) | 2017-10-13 |
CN107247823B true CN107247823B (zh) | 2020-08-18 |
Family
ID=60016925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710342079.XA Active CN107247823B (zh) | 2017-05-12 | 2017-05-12 | 基于准一维化学动力学过程和自模性的喷焰流场预估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107247823B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108647419B (zh) * | 2018-05-03 | 2022-01-28 | 北京环境特性研究所 | 一种随高度变化低空喷焰红外辐射特性预估方法及装置 |
CN108595867B (zh) * | 2018-05-03 | 2022-03-08 | 北京环境特性研究所 | 一种伴随流影响下低空喷焰红外辐射特性预估方法及装置 |
CN110083861A (zh) * | 2019-03-18 | 2019-08-02 | 南京航空航天大学 | 一种构建碳氢燃料射流火焰拓扑结构的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005249418A (ja) * | 2004-03-01 | 2005-09-15 | Mazda Motor Corp | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
CN102880734A (zh) * | 2012-06-21 | 2013-01-16 | 中国人民解放军电子工程学院 | 基于cfd的飞机尾喷流大气扩散建模方法 |
CN103279616A (zh) * | 2013-06-03 | 2013-09-04 | 西安近代化学研究所 | 一种固体推进剂烟雾影响可见光制导信号虚拟试验方法 |
CN103569373A (zh) * | 2013-11-13 | 2014-02-12 | 中国航空工业集团公司西安飞机设计研究所 | 涡桨飞机在静止状态下尾喷流速度场的计算方法 |
CN104123414A (zh) * | 2014-07-19 | 2014-10-29 | 安徽农业大学 | 一种气动发动机进排气通道流场和热传递优化的参数化方法 |
-
2017
- 2017-05-12 CN CN201710342079.XA patent/CN107247823B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005249418A (ja) * | 2004-03-01 | 2005-09-15 | Mazda Motor Corp | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
CN102880734A (zh) * | 2012-06-21 | 2013-01-16 | 中国人民解放军电子工程学院 | 基于cfd的飞机尾喷流大气扩散建模方法 |
CN103279616A (zh) * | 2013-06-03 | 2013-09-04 | 西安近代化学研究所 | 一种固体推进剂烟雾影响可见光制导信号虚拟试验方法 |
CN103569373A (zh) * | 2013-11-13 | 2014-02-12 | 中国航空工业集团公司西安飞机设计研究所 | 涡桨飞机在静止状态下尾喷流速度场的计算方法 |
CN104123414A (zh) * | 2014-07-19 | 2014-10-29 | 安徽农业大学 | 一种气动发动机进排气通道流场和热传递优化的参数化方法 |
Non-Patent Citations (2)
Title |
---|
The Similarity Parameters of Subscale Test for Simulating Dynamics of Launch Combustion Gas Jet;Chen,J.等;《Acta Aerodynamic Sinica》;20050930;第23卷;第307-311页 * |
导弹喷焰流场与凝结尾迹的工程计算方法;郭琼等;《弹箭与制导学报》;20160430;第36卷;第83-86页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107247823A (zh) | 2017-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107247823B (zh) | 基于准一维化学动力学过程和自模性的喷焰流场预估方法 | |
Staats et al. | Active control of the corner separation on a highly loaded compressor cascade with periodic nonsteady boundary conditions by means of fluidic actuators | |
CN109460626B (zh) | 冲压发动机性能参数计算方法 | |
CN108182319B (zh) | 一种超声速一体化喷管设计方法 | |
CN102880734A (zh) | 基于cfd的飞机尾喷流大气扩散建模方法 | |
Natta et al. | Flow analysis of rocket nozzle using computational fluid dynamics (CFD) | |
CN107672821A (zh) | 基于吻切锥理论的变马赫数“并联”宽速域乘波飞行器设计方法 | |
CN105138787A (zh) | 基于特征线追踪的超声速流场设计方法 | |
Patel | Flow analysis and optimization of supersonic rocket engine nozzle at various divergent angle using Computational Fluid Dynamics (CFD) | |
Bons et al. | The effect of real turbine roughness with pressure gradient on heat transfer | |
CN111079325B (zh) | 基于代理模型的涡扇发动机喷流噪声实时计算及预测方法 | |
CN115329489A (zh) | 一种曲率连续的拉瓦尔喷管设计方法 | |
Llopis-Pascual | Supercritical Coanda jets for flight control effectors | |
Choi et al. | Efficient diffuser design for plasma wind tunnels with a large blockage model | |
CN115879216A (zh) | 一种内流道强波系干扰控制下的流场重构设计方法 | |
CN110162901A (zh) | 高超声速飞行器轴对称构型前体的优化设计方法及系统 | |
Babu et al. | Parametric study of the performance of two-dimensional scramjet intake | |
Cao et al. | Experimental study of the formation and diffusion of soap bubbles for wrapping and transporting indoor airborne pollution | |
Vyas et al. | Effect of transonic inlet design on the performance of a mico-turbojet engine | |
Chernyshev et al. | Research of the flow-over of the mechanical core of four pipes by a wind current, using Control Volume Method | |
Mouton | Simulation of atmospheric boundary layer for air inlet testing in crossflow | |
Tripathi et al. | Drag mitigation of trilobed airship hull through aerodynamic comparison with conventional single-lobed hull | |
Aziz et al. | High fidelity design optimization of a three dimension supersonic intake | |
Xun et al. | A new design and simulation of gas turbine exhaust ejector | |
Mern et al. | A study of numerical simulation of supersonic conical nozzle exhaust |
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 |