CN111368430A - 一种垂直发射航行体水中点火燃气后效估计方法 - Google Patents
一种垂直发射航行体水中点火燃气后效估计方法 Download PDFInfo
- Publication number
- CN111368430A CN111368430A CN202010142835.6A CN202010142835A CN111368430A CN 111368430 A CN111368430 A CN 111368430A CN 202010142835 A CN202010142835 A CN 202010142835A CN 111368430 A CN111368430 A CN 111368430A
- Authority
- CN
- China
- Prior art keywords
- navigation body
- equation
- model
- phase
- pressure
- 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
- 238000000034 method Methods 0.000 title claims abstract description 51
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims description 28
- 230000000694 effects Effects 0.000 title description 3
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims description 27
- 239000000203 mixture Substances 0.000 claims description 24
- 230000033001 locomotion Effects 0.000 claims description 17
- 238000013508 migration Methods 0.000 claims description 8
- 230000005012 migration Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 8
- 239000012530 fluid Substances 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 239000013535 sea water Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 4
- 230000016507 interphase Effects 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 4
- 210000000476 body water Anatomy 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 abstract description 7
- 239000007789 gas Substances 0.000 description 23
- 239000002737 fuel gas Substances 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 6
- 230000005484 gravity Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 239000007921 spray Substances 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Images
Landscapes
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明首先构建垂直发射航行体水中点火数值仿真模型,然后利用航行体水下点火压力和姿态数据对数值仿真模型校验和完善,在此基础上,研究航行体水中点火过程中尾部流场变化情况,分析不同点火深度点火发动机射流流场对发射平台壁面、发射筒的影响,同时对不同点火对航行体水下发射能力的影响进行估计。
Description
技术领域
本发明属于垂直发射航行体水中点火燃气后效的预测方法,尤其涉及燃气后效中的压力、温度对水中航行体、发射筒及发射平台影响的估计方法。
背景技术
对于采用水中点火的垂直发射航行体而言,通过及时对航行体的水中弹道和姿态进行控制,可提高航行体在水中运动的稳定性,并有助于获得有利的出水姿态,从而增强航行体出水后的姿态控制能力。通常,水下发射航行体在其尾部离开发射筒口一定距离处点火进行发射。如果航行体尾部与发射筒口距离较近处点火,则发动机可以利用从发射筒内溢出并附着在航行体尾部的燃气泡,作为发动机喷流建立初期燃气的受纳空间,减轻发动机直接在水中点火所造成的冲击,提高发动机工作安全性,但是由于水下点火处距离发射平台较近,发动机高速喷流会对筒盖及周围发射平台结构产生破坏,影响后续发射。如果航行体远离发射筒口处(简称筒口)点火,附着在航行体尾部的燃气泡体积变小,同时燃气泡内渗人水,发动机若直接在此环境中点火,由于水的巨大惯性,则发动机喷流会受到水的阻碍,导致喷管内压强过高,从而威胁发动机工作安全性。显然,在决定采用水中点火方案时,必须明确水中点火对航行体载荷的影响,并充分估计水中点火对发射筒及发射平台所造成的威胁。本发明中,燃气后效由垂直发射航行体水下点火后产生的压力和温度参数来表征。流场参数是指航行体周边空间中不同位置点的温度和压力。航行体运动参数是指航行体的速度和位置。
发明内容
本发明首先构建垂直发射航行体水中点火数值仿真模型,然后利用航行体水下点火压力和姿态数据对数值仿真模型校验和完善,在此基础上,研究航行体水中点火过程中尾部流场变化情况,分析不同点火深度下发动机射流流场对发射平台壁面和发射筒的影响,同时对水中点火燃气后效进行估计。
本发明的具体步骤为:
(一)航行体动力方程及仿真模型构建。建立水下发射过程中的航行体动力模型,通过实时求解航行体受力,通过动力学方程求得加速度和角加速度,然后进行数值积分,获得航行体当前速度、位移及姿态等。航行体受力如图1所示,图1中oxyz为发射坐标系,相对于地球静止,ox1y1z1为航行体坐标系,固联于航行体;v为航行体运动速度;F为航行体所受的浮力;G为自身重力;f为航行体轴向阻力;ox1为航行体轴向,ox1与v的夹角为α;oy1与v的夹角为T为发动机产生的总推力;φ为发动机喷管为纠正航行体出水姿态而偏转的摆角。
航行体在水下发射过程中,首先利用高压燃气将航行体弹射出筒,航行体尾部离开发射筒后,筒内燃气进入外界环境,与周围的水汽混合并附着在尾部形成燃气泡。发动机接收到点火指令后,航行体运动的发射动力主要由两部分组成:弹射力和发动机提供的轴向推力。
航行体在水下发射过程中,其航行体动力模型为:
式中:t为航行体运动时间,F′为发射筒高压燃气弹射力,T′为发动机点火后航行体受到的轴向推力,m为航行体质量,s为航行体的迎流面积,C为阻力系数,ρ为海水密度,λ为附加质量;α和值利用航行体加速度计测量的各个方向的加速度值进行解算求取。式(1)中的航行体运动参数值为后续计算燃气后效参数提供动态输入。
本发明利用计算流体动力学(Computational Fluid Dynamics,CFD)模型,在CFD模型中引入经验参数,基于已有的实验数据对CFD模型参数进行合理调整,形成合理的计算模型。
(二)对计算域进行网格划分并进行CFD求解。将式(1)和式(2)所示的航行体动力模型输入到CFD求解器中,利用CFD求解器对航行体水下运动参数进行求解。在得到航行体运动参数后对网格节点参数进行更新,进而对航行体能量方程及控制方程进行数值离散求解,最终得出航行体近体流场压力及温度分布参数。由于模型的对称性,计算域模型采取二维轴对称模型。在建模时不考虑航行体尾部产生的空化现象,且由于研究的是尾部流场变化,头部对仿真流场的影响较小,因此将计算模型简化为尾部只有喷管扩张段的无限长圆柱体,综上所述,计算模型采用尾部只有喷管扩张段的无限长二维轴对称圆柱体。CFD计算模型具体采用二维轴对称的无限长圆柱体模型。CFD求解器计算域如图2所示。网格划分的原则是,全部流场区域采用结构化网格。该方式能够消除计算结果对网格的依赖性。在划分喷管时,在轴向方向、喷管及尾部附近处的网格大小接近一致。外场划分梯度网格,靠近航行体位置网格细密,远场较为稀疏。按照以上要求划分网格,得到航行体尾部局部网格如图3所示。
(三)利用CFD方法求解燃气后效参数。本发明采用Mixture多相流模型对流场进行求解,其控制方程基本形式如下:
式中,ρm为混合物密度:
式中,αk为第k相的体积分数。利用Mixture多相流模型进行求解时,分别设置了主相(Primary phase)和副相(Secondary phase),所有相体积分数之和为1,副相可为多个,对于副相p的体积分数求解方法为:
动量方程是牛顿第二定律在流体动力学中的体现。对于Mixture多相流模型的动量方程,可通过将所有相各自的动量方程相加的方法获得,其表示式为:
能量方程是能量守恒定律在流体动力学中的体现,Mixture多相流模型中混合物的能量方程为:
式中,keff为有效传导系数,表达式为:
其中,kt为湍动热传导系数,kk为第k相的湍动能,其取值取决于所选择的湍流模型。式(13)右端项中第一项代表了传导引起的能量迁移,SE为其余的体积热源项。
对于可压缩流动,Ek的表达式为:
对于不可压缩流动,表达式为:
Ek=hk (16)
其中,hk为第k相的焓。
上述偏微分方程组,包括式(3)至式(16),通过网格节点在计算域流场内进行离散化求解,具体步骤如下:
(1)在计算过程的每个时间步,利用CFD方法联合求解Mixture模型控制方程、动量方程、能量方程和航行体动力模型,获得流场和航行体运动参数,通过更新网格将二者进行耦合,计算完毕后转入下一时间步;
(2)当第一步得到的速度在局部不满足连续方程时,从连续性方程和线化动量方程推导出压力校正的泊松方程,然后解出压力校正方程,获取压力和速度场;
(3)用上一个时间步长更新的、除了压力速度和温度外的其它变量值,解出湍流、能量和辐射等标量;
(4)当存在相间耦合时,用离散相轨迹计算来更新连续相的源项;
(5)根据计算的残差曲线检查设定方程的收敛性,当所有变量的残差值都降到10-3时,就认为计算收敛,即完成了对压力和温度的离散求解。
表1给出了一组计算初始参数设置。
表1主要参数设置
在水域边界设置中,初始压力输入数值通过重力梯度法获得。重力梯度法就是根据远场边界网格节点坐标和航行体发射水深来计算原场边界压力输入的方法。设标准大气压为P0=101325Pa,重力加速度g=9.81m/s2,海水密度ρ=1.02Kg/m3,则初始远场压力总压,简称为总压,为P=P0+ρgh,其中,h为压力边界网格节点高度与航行体发射深度的差值。总温及回流总温初始值由发射筒内温度传感器读取。
在CFD模型设置中,湍流模型采用RNGk-ε,可以达到较好的收敛性;由于流场变量在壁面附近存在很大梯度的流动,因此壁面函数采用非平衡壁面函数法;流场初始化过程中,先对整流场按表1进行参数初始化,而后标记出全部为水的区域并对此区域按表1中初始项给出的参数修正,再迭代计算;在计算步长选择上,首先保证计算精度,其次要保证单时间步迭代计算的收敛性,最后要考虑计算的经济性。
在对航行体近体流场进行离散求解后,即可得出流场属性的关键参数(如压力、密度和温度等)。该模型具体求解流程如图4所示。
综上所示,本发明的具体步骤为:
一种垂直发射航行体水中点火燃气后效估计方法,其步骤包括:
首先建立航行体在水下发射过程中的航行体动力模型,
其中,式中,t为航行体运动时间,F′为发射筒高压燃气弹射力,T′为发动机点火后航行体受到的轴向推力,m为航行体质量,s为航行体的迎流面积,C为阻力系数,ρ为海水密度,λ为附加质量;α和值利用航行体加速度计测量的各个方向的加速度值进行解算求取;
将航行体动力模型输入到CFD求解器中,利用CFD求解器对航行体水下运动参数进行求解;在得到航行体运动参数后对网格节点参数进行更新,进而对航行体能量方程及边界条件控制方程进行数值离散求解,最终得出航行体近体流场压力及温度分布参数;
采用Mixture模多相流型对流场进行求解,其控制方程基本形式为:
其中,为按质量平均的速度:式中,ρm为混合物密度:式中,αk为第k相的体积分数;利用Mixture多相流模型进行求解时,分别设置了主相和副相,所有相体积分数之和为1,副相为多个,对于副相p的体积分数求解方法为:
式中,keff为有效传导系数,表达式为:其中,kt为湍动热传导系数,kk为第k相的湍动能,其取值取决于所选择的湍流模型;代表了传导引起的能量迁移,SE为其余的体积热源项;对于可压缩流动,Ek的表达式为:对于不可压缩流动,表达式为:Ek=hk,其中,hk为第k相的焓。上述偏微分方程组通过网格节点在计算域流场内进行离散化求解,具体步骤为:在计算过程的每个时间步,利用CFD方法联合求解Mixture模型、动量方程、能量方程和航行体动力模型,获得流场和航行体运动参数,通过更新网格将二者进行耦合,计算完毕后转入下一时间步;当第一步得到的速度在局部不满足连续方程时,从连续性方程和线化动量方程推导出压力校正的泊松方程,然后解出压力校正方程,获取压力和速度场;用上一个时间步长更新的、除了压力速度和温度外的其它变量值,解出湍流、能量和辐射等标量;当存在相间耦合时,用离散相轨迹计算来更新连续相的源项;根据计算的残差曲线检查设定的方程的收敛性,当所有变量的残差值都降到10-3时,就认为计算收敛,即完成了对压力和温度的离散求解。
本发明的有益效果主要集中在以下三个方面:
①本发明对航行体水中点火产生的燃气后效现象进行了分析,通过建立航行体水下运动方程组及轴向阻力计算模型,基于Fluent软件,采用动网格技术对航行体离开发射筒及发动机水中点火燃气后效过程进行了分析。
②本发明通过分析不同点火距离燃气后效产生的压力、温度对发射筒及发射平台的影响,给出了具体的估计结果,为发射筒的设计及发射平台安全性评估提供了依据。
③本发明分析了不同点火距离发动机燃气射流流量和速度的变化规律,综合对比评估发动机不同时刻喷管轴线压力特征,为水下航行体选择合适的水中点火距离提供科学的数据支撑和理论依据。
附图说明
图1为发射坐标系和航行体坐标系定义。
图2为CFD模型的求解区域。
图3为航行体尾部局部网格。
图4为数值仿真模型求解流程。
图5为发射平台壁面监测点压力变化规律。
图6为发射平台壁面监测点温度变化规律。
具体实施方式
这里结合实施例对本发明过程进行具体说明。
一种垂直发射航行体水中点火燃气后效估计方法,其步骤包括:
首先建立航行体在水下发射过程中的航行体动力模型,
其中,式中,t为航行体运动时间,F′为发射筒高压燃气弹射力,T′为发动机点火后航行体受到的轴向推力,m为航行体质量,s为航行体的迎流面积,C为阻力系数,ρ为海水密度,λ为附加质量;α和值利用航行体加速度计测量的各个方向的加速度值进行解算求取;
将航行体动力模型输入到CFD求解器中,利用CFD求解器对航行体水下运动参数进行求解;在得到航行体运动参数后对网格节点参数进行更新,进而对航行体能量方程及边界条件控制方程进行数值离散求解,最终得出航行体近体流场压力及温度分布参数;
采用Mixture模多相流型对流场进行求解,其控制方程基本形式为:
其中,为按质量平均的速度:式中,ρm为混合物密度:式中,αk为第k相的体积分数;利用Mixture多相流模型进行求解时,分别设置了主相和副相,所有相体积分数之和为1,副相为多个,对于副相p的体积分数求解方法为:
式中,keff为有效传导系数,表达式为:其中,kt为湍动热传导系数,kk为第k相的湍动能,其取值取决于所选择的湍流模型;代表了传导引起的能量迁移,SE为其余的体积热源项;对于可压缩流动,Ek的表达式为:对于不可压缩流动,表达式为:Ek=hk,其中,hk为第k相的焓。上述偏微分方程组通过网格节点在计算域流场内进行离散化求解,具体步骤为:在计算过程的每个时间步,利用CFD方法联合求解Mixture模型、动量方程、能量方程和航行体动力模型,获得流场和航行体运动参数,通过更新网格将二者进行耦合,计算完毕后转入下一时间步;当第一步得到的速度在局部不满足连续方程时,从连续性方程和线化动量方程推导出压力校正的泊松方程,然后解出压力校正方程,获取压力和速度场;用上一个时间步长更新的、除了压力速度和温度外的其它变量值,解出湍流、能量和辐射等标量;当存在相间耦合时,用离散相轨迹计算来更新连续相的源项;根据计算的残差曲线检查设定的方程的收敛性,当所有变量的残差值都降到10-3时,就认为计算收敛,即完成了对压力和温度的离散求解。
不同位置点火时,发动机喷管所黏附的燃气泡体积不同,如果燃气泡内点火,则可以起到缓冲作用,从而减弱对航行体、发射平台的冲击载荷。本发明中,通过设置不同的点火距离,分析燃气后效对点火、发射平台壁面和发射筒的影响。
采用如表1所示的参数来设置CFD计算模型。假设发动机距筒口点火的距离分别为2.2、4.3、6.5(数值进行了无量纲化处理),对应时间为0.4、0.6、0.8(数值进行了无量纲化处理)。在点火时刻,Fluent设置中将喉部的边界wall改为压力入口(PressureInlet),通过压力入口参数设置选项设置航行体发动机喉部入口压力和温度参数如表2所示,通过Fluent自定义函数(UDF)定义航行体质量和转动惯量等参数,通过网格节点离散求解航行体水动力及压力等参数,湍流模型选用增强RNGk-ε模型,计算迭代步长为1E-5。
表2喉部入口参数(数值进行了无量纲化处理)
项目 | 量值 |
总压(Gauge Total Pressure) | 6 |
静压(Initial Gauge Pressure) | 3 |
总温(Total Temperature) | 3300 |
这里对近筒口距离2.2处点火燃气后效进行了估计。
通过自定义函数UDF设置压力监测点位置,利用本发明方法得到的发射平台壁面的压力与温度的变化规律曲线分别如图5和图6所示。
从图5中可以看出,发射平台壁面各个测点的压力脉动变化规律基本一致,且离筒口中心越远,脉动幅值越小,当t>0.7后,压力值基本维持在当地水深压力的0.55倍;在点火后,测点Pb1的压力值出现3次压力峰值,这是由于Pb1是发射筒与发射平台壁面的交界点,发动机燃气射流与发射筒内高压燃气同时对Pb1产生作用。从图6中可以看出,发射平台壁面各测点温度沿壁面的传递作用与压力传递不同,热传递的作用距离很近,除了测点Pb1,其余测点温度基本不受影响,均保持在300左右;对于测点Pb1的温度变化,主要是由于喷管喷出的燃气其高温主要集中分布在轴线处,温度传播通过燃气之间的传递,在t=0.2时刻以后,发射筒的Tb1点温度才开始升高,此后出现几个温度上升的过程;在t=0.641时刻,燃气已经收缩至发射筒内,此后测点温度回归300。
Claims (6)
1.一种垂直发射航行体水中点火燃气后效估计方法,其步骤包括:
首先建立航行体在水下发射过程中的航行体动力模型,
式中,t为航行体运动时间,F′为发射筒高压燃气弹射力,T′为发动机点火后航行体受到的轴向推力,m为航行体质量,s为航行体的迎流面积,C为阻力系数,ρ为海水密度,λ为附加质量;α和θ值利用航行体加速度计测量的各个方向的加速度值进行解算求取;
将航行体动力模型输入到CFD求解器中,利用CFD求解器对航行体水下运动参数进行求解;在得到航行体运动参数后对网格节点参数进行更新,进而对航行体能量方程及控制方程进行数值离散求解,最终得出航行体近体流场压力及温度分布参数;
采用Mixture模多相流型对流场进行求解,其控制方程基本形式为:
其中,为按质量平均的速度:式中,ρm为混合物密度:式中,αk为第k相的体积分数;利用Mixture多相流模型进行求解时,分别设置了主相和副相,所有相体积分数之和为1,副相为多个,对于副相p的体积分数求解方法为:
式中,keff为有效传导系数,表达式为:其中,kt为湍动热传导系数,kk为第k相的湍动能,其取值取决于所选择的湍流模型;代表了传导引起的能量迁移,SE为其余的体积热源项;对于可压缩流动,Ek的表达式为:对于不可压缩流动,表达式为:Ek=hk,其中,hk为第k相的焓;
上述偏微分方程组通过网格节点在计算域流场内进行离散化求解,具体步骤为:在计算过程的每个时间步,利用CFD方法联合求解Mixture模型控制方程、动量方程、能量方程和航行体动力模型,获得流场和航行体运动参数,通过更新网格将二者进行耦合,计算完毕后转入下一时间步;当第一步得到的速度在局部不满足连续方程时,从连续性方程和线化动量方程推导出压力校正的泊松方程,然后解出压力校正方程,获取压力和速度场;用上一个时间步长更新的、除了压力速度和温度外的其它变量值,解出湍流、能量和辐射等标量;当存在相间耦合时,用离散相轨迹计算来更新连续相的源项;根据计算的残差曲线检查设定方程的收敛性,当所有变量的残差值都降到10-3时,就认为计算收敛,即完成了对压力和温度的离散求解。
2.一种如权利要求1所述的垂直发射航行体水中点火燃气后效估计方法,其特征在于,CFD计算模型采用二维轴对称的无限长圆柱体模型,全部流场区域采用结构化网格。
3.一种如权利要求2所述的垂直发射航行体水中点火燃气后效估计方法,其特征在于,在划分喷管时,在轴向方向、喷管及尾部附近处的网格大小接近一致。
4.一种如权利要求2所述的垂直发射航行体水中点火燃气后效估计方法,其特征在于,外场划分梯度网格,靠近航行体位置网格细密,远场较为稀疏。
5.一种如权利要求1所述的垂直发射航行体水中点火燃气后效估计方法,其特征在于,在CFD模型设置中,湍流模型采用RNGk-ε。
6.一种如权利要求1所述的垂直发射航行体水中点火燃气后效估计方法,其特征在于,在CFD模型设置中,壁面函数采用非平衡壁面函数法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010142835.6A CN111368430B (zh) | 2020-03-04 | 2020-03-04 | 一种垂直发射航行体水中点火燃气后效估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010142835.6A CN111368430B (zh) | 2020-03-04 | 2020-03-04 | 一种垂直发射航行体水中点火燃气后效估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111368430A true CN111368430A (zh) | 2020-07-03 |
CN111368430B CN111368430B (zh) | 2022-09-30 |
Family
ID=71208475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010142835.6A Active CN111368430B (zh) | 2020-03-04 | 2020-03-04 | 一种垂直发射航行体水中点火燃气后效估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111368430B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113947039A (zh) * | 2021-09-24 | 2022-01-18 | 哈尔滨工程大学 | 一种对航行体出水的运动与其尾空泡的发展进行预报的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO20140079A1 (no) * | 2014-01-24 | 2015-07-27 | Ledaflow Technologies Da | Fremgangsmåte for bestemmelse av strømningsparametere for en flerfasestrømning i en rørledningsseksjon |
-
2020
- 2020-03-04 CN CN202010142835.6A patent/CN111368430B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO20140079A1 (no) * | 2014-01-24 | 2015-07-27 | Ledaflow Technologies Da | Fremgangsmåte for bestemmelse av strømningsparametere for en flerfasestrømning i en rørledningsseksjon |
Non-Patent Citations (3)
Title |
---|
刘可等: "发射速度对航行体出水姿态影响数值仿真研究", 《舰船电子工程》 * |
燕国军等: "水环境下固体火箭发动机喷流流场数值研究", 《固体火箭技术》 * |
王鹏等: "航行体水下点火三维流场数值模拟", 《西北工业大学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113947039A (zh) * | 2021-09-24 | 2022-01-18 | 哈尔滨工程大学 | 一种对航行体出水的运动与其尾空泡的发展进行预报的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111368430B (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103970957B (zh) | 一种弹性乘波体高超声速飞行器仿真方法 | |
Qu et al. | Numerical investigation of the aerodynamics of a delta wing in ground effect | |
CN111709196A (zh) | 一种水下高速航行体抗空化空蚀能力评估方法 | |
CN112651083A (zh) | 一种基于fluent仿真的爆震激光推力器分析方法 | |
CN111368430B (zh) | 一种垂直发射航行体水中点火燃气后效估计方法 | |
Ocokoljić et al. | Aerodynamic shape optimization of guided missile based on wind tunnel testing and computational fluid dynamics simulation | |
CN116702653A (zh) | 火箭燃气射流冲击垂直发射装置流场仿真方法 | |
CN113515900A (zh) | 模拟连续旋转爆震发动机尾流的数值仿真方法 | |
Volkov et al. | Transverse injection of a jet from the surface of a flat plate into the supersonic flow over it | |
Cheng et al. | Modeling tail bubble dynamics during the launch of an underwater vehicle using the boundary element method | |
Yang et al. | Water-exit process modeling and added-mass calculation of the submarine-launched missile | |
Cummings et al. | Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies | |
Zheng et al. | Flutter analyses of complete aircraft based on hybrid grids and parallel computing | |
CN114077770A (zh) | 一种用于飞行器水上迫降sph模拟的直接造波方法 | |
Li et al. | 3-D Simulation of external cooling of aero-optical side window | |
Madhulaalasa et al. | Effect of the Leading-Edge Sweep on Wrap-Around Fins | |
Kwak et al. | Convergence study of inlet buzz frequency with computational parameters | |
Bibin et al. | Flow analysis of secondary air injection in conical rocket nozzle | |
Zhang et al. | 2D axisymmetric CFD simulation of underwater torpedo launch tube flow | |
Liu et al. | Unsteady numerical simulation of jet interaction by coupling aerodynamics and flight dynamics | |
Duan et al. | On the transient numerical simulation of solid rocket motor by coupling quasi one-dimension internal flow with three-dimension propellant grain combustion | |
Gu et al. | Self-Turning Process and Aerodynamic Characteristics of Slender Bodies with Different Center-of-Gravity Positions | |
Zhang et al. | Research on underwater vehicle based on multiphase flow control | |
CN117313587B (zh) | 一种内孤立波与背景剪切流相互作用模拟方法及其系统 | |
Liu et al. | Underwater ejection multifield coupling model and response characteristics |
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 |