导弹水下发射燃气泡计算

采用一维非定常气流场模型和轴对称理想水流场模型，对水下发射导弹的尾部非定常燃气泡内外流场进行了耦合数值求解。考虑了高温燃气与水介质之间的传热、汽化等对泡内气体流动的影响，数值模拟了燃气泡的生长和脱落过程，揭示了燃气泡中气体流动的基本规律，计算了燃气泡发展对导弹所受水动力的影响。     

基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。     

不同发射深度下导弹水下点火气水流体动力计算

从流体动力角度研究了不同发射深度下，导弹水下点火这一非定常非线性过程。整个系统分为外部水流场、喷管流场和燃气泡流场三个区域加以考虑。水流场采用不可压势流模型，用边界元方法求解；喷管内流场采用非定常一元流动模型，用特征线差分法求解，并设置了激波检测功能；燃气泡采用基于质量和能量守恒的零维计算模型。在时间域中用步进方法实现了三个流场的耦合求解。给出了四种发射深度下的数值计算结果，展示了导弹水下点火的一     

CE/SE方法在内弹道瞬态两相流中的应用

根据CE/SE方法思想,构造在变换坐标系下内弹道两相流动数值计算的 格式. 采用此格式研究了膛内瞬态二维两相流动点火过程. 数值计算结果反映实际 规律并和试验结果符合. 结果表明,对中心点火方式,在点传火过程初期,膛内压力在 径向和轴向上存在较大的压力梯度,当火药全部着火后,药室内压力在径向上逐渐趋于 一致. 这表明CE/SE方法能够完全模拟瞬态二维两相流动的情况.     

膛口反应流并行数值模拟

采用轴对称多组分N-S方程对含有高速运动弹丸的膛口反应流进行了数值模拟.控制方程采用时间分裂方法并在大型计算机上采用MPI方法进行多核并行求解,其中对流项采用二阶AUSM+格式和MUSCL插值方法进行处理,燃气采用氢气-空气混合气,反应机理为9组分19步基元反应.对于弹丸引起的网格运动,采用嵌套网格法处理.并行验证算例与串行计算结果一致,采用20个CPU计算时效率为64％.根据数值结果详细讨论了发射过程中的气体动力学和化学动力学过程,并且通过对两种条件下的计算结果比较分析了化学反应对膛口流场发展的影响.结果表明,上述算法能够较为正确地模拟弹丸和化学反应对膛口流场的影响,并大大提高了计算速度.     

燃爆压裂中压挡液柱运动规律的动力学模型

针对燃爆压裂过程中压挡液柱受冲击运动机理的复杂性,假设火药燃气与压挡液柱存在完全气液 接触界面,采用拉格朗日的微元分析方法,建立了由连续性方程、动量守恒方程、能量守恒方程组成的压挡液 柱运动规律动力学模型,并给出了该模型与火药燃爆模型的耦合数值解法。经程序编制和实例计算表明,在 综合考虑火药燃气对液柱的宏观推动作用、冲击压缩作用、液柱自身的动能分布及管壁对其摩擦阻力的影响 后,火药燃烧过程中气液界面上升高度有限(实例计算不足0.1m),可起到很好的持压作用;但全过程中最高 液柱位移较大(18.9m),水力振荡增效作用明显。研究成果对提高燃爆压裂的数值模拟精度具有一定促进 作用。     

考虑两相流内弹道的自行火炮发射动力学计算

为了更加准确地计算火炮弹丸的起始扰动,将两相流内弹道应用于发射动力学研究。建立了完整的自行火炮系统发射动力学方程组,包括火炮系统的体动力学方程组、弹丸在膛内运动的动力学方程以及两相流内弹道方程。编制了计算程序,实现了对某自行火炮发射过程的数值模拟,在准确计算内弹道过程的同时,获得了火炮的动力响应、弹丸膛内运动和起始扰动。部分模拟结果与实验实测结果吻合较好。计算表明,采用两相流内弹道模型将提高发射动力学计算精度。     

?????????????????????????????????

以某跨声速轴流压气机单转子为研究对象,应用数值模拟技术,采用全通道 计算方案,利用Jameson有限体积中心差分格式并结合Spalart-Allmaras 湍流模型获得进口畸变条件下该轴流压气机转子性能和内部流动细节,详细分析了进气周向 总压畸变对压气机转子内流场流动结构的影响. 并将计算结果与实验结果进行了比较, 结果表明,数值模拟的结果与实验结果符合较好,计算方案切实可行.     

静不稳定空空导弹导轨发射机弹分离数值模拟

将基于N-S方程求解的流动模拟与求解导弹六自由度运动方程和控制律相结合,进行静不稳定空空导弹导轨发射分离轨迹的数值模拟.结果表明,对静不稳定导弹,如果不施加控制,载机的横侧向气动干扰会使采用导轨发射的导弹在离轨后滚转姿态很快发散,俯仰和偏航姿态快速改变,施加姿态控制后,导弹离开载机的过程中姿态比较稳定,改善了发射安全性.     

双燃式超燃发动机冷态内流场的数值研究

研究了双燃式一体化通道(包含进气道、双燃式燃烧室和尾喷管)的冷态内流场特性.首次在激波风洞中对内流场进行纹影照相,用TVD格式求解三维全N-S方程对喷管和一体化通道进行分区数值模拟,并考察了几何参数对内流场的影响.结果表明对典型工况(h     

二级轻气炮内弹道的数值模拟与性能分析

为了优化二级轻气炮的设计参数和工作参数,必须进行数值模拟。本文介绍了准一维不定常流的Lagrangian算法,讨论了数值模拟中的若干问题,分析了气动过程。根据我们的数值模拟,归纳了影响气炮性能的一些基本规律,可供选择工作参数时综合分析。     

高压气体驱动二级轻气炮发射过程的数值模拟方法及应用

二级轻气炮是一种常见的超高速发射装置,多年来其数值研究大多采用简化一维模型,鲜有三维有限元模型。以14 mm口径高压气体驱动二级轻气炮为研究对象,采用耦合欧拉-拉格朗日（coupled Eulerian-Lagrangian, CEL）算法,根据膜片破裂与否,将二级轻气炮模型解耦为2个分级三维数值模型。为确定实验难以测得的参数（材料摩擦因数和膜片破膜压力）,设计正交试验,拟合确定活塞与泵管间摩擦因数为0.82,弹丸与发射管摩擦因数为0.30和膜片破膜压力为11.73 MPa。正交结果表明,摩擦因数对计算结果影响较大,在高压气体驱动二级轻气炮的计算中不应忽略。通过上述方法建立数字化高压气体驱动二级轻气炮,完整复现气炮发射过程,计算的弹丸终速与实验结果吻合度高。选取验证工况详细分析了气炮发射过程内流场变化,并呈现关键时刻的压力云图。该气炮简化方法、分级思想和关键参数确认方法可推广应用于固体发射药驱动、爆轰驱动等其他驱动形式的二级/多级轻气炮。     

A ballistic compressor-based experiment for the visualization of liquid propellant jet combustion above 100 MPa

 This paper describes the components and operation of an experimental setup for the visualization of liquid propellant (LP) jet combustion at pressures above 100 MPa. The apparatus consists of an in-line ballistic compressor and LP injector. The ballistic compressor, based on a modified 76 mm gun, provides high-pressure (ca. 55 MPa) clear hot gas for the jet ignition. A piston (projectile) is fired toward a test chamber beyond the barrel’s end, and its rebound is arrested in a transition section that seals the test chamber to the barrel. The LP jet is injected once the piston is restrained, and combustion of the jet further elevates the pressure. At a preset pressure, a disc in the piston ruptures and the combustion gas vents sonically into the barrel. If a monopropellant is used, the jet injection-combustion process then resembles liquid rocket combustion but at very high pressures (ca. 140 MPa). This paper discusses the ballistics of the compression and compares experimental results to those predicted by a numerical model of the apparatus. Experimentally, a pressure of 70 MPa was achieved upon a 12.5 volumetric compression factor by firing a 10 kg piston into 1.04 MPa argon using a charge of 75 g of small-grain M1 propellant. Received: 16 December 1996/Accepted: 15 July 1997     

Numerical Modeling of Interior Ballistics Processes in Light Gas Guns

An improved procedure for modeling the operation of a light-gas gun is proposed. The motion of working bodies in both the firing chamber and the light-gas chamber is studied within the framework of the mechanics of heterogeneous media. The problem of barrel heating taking into account its melting and removal of thermal ablation products into the medium inside the bore is solved in a coupled formulation. Heat and mass transfer and friction on the barrel surface are calculated using empirical dependences. The deformable piston is considered compressible and elastoviscoplastic. Allowance is made for the presence of a clearance between the lateral surface of the piston and the barrel bore walls and the associated gas flow between the firing and the light-gas chamber. Calculation results are given.     

Two-phase flows in gun barrel: Theoretical and experimental studies

In the paper theoretical and numerical model of two-phase flow of solid granular propellant and its products of combustion in the gun barrel during interior ballistic cycle is given. Two cases are considered: base ignition of propellant charge and ignition by igniter. The theoretical model includes the balance equations of mass, momentum and energy for both phases, as well as necessary constitutive laws. The igniter efflux in the propellant chamber is obtained by incorporation in the model the two-phase flow model of igniter function. The convergent, unconditionally stable, numerical procedure is formed to solve the system of equations of the theoretical model. An original procedure of numerical grid adaptation to the flow field increase, caused by the projectile motion down the gun bore, is developed. The TWOPIB code for the computation of whole interior ballistic cycle of ammunition is developed. Four kinds of experimental investigations were carried out:igniter function in open air, flamespreading through propellant charge in the fibreglass tube during base ignition or during ignition by igniter, and firing of 100 mm APFSDS projectile. Verification of the theoretical–numerical approach by the comparison with experimental data is carried out. The great number of computational results is presented for the parameters that can not be measured, but which are necessary for more complete understanding of examined processes. The presented theoretical–numerical access enables, not only the complete optimisation of propellant charges, but more successful solutions of many interior ballistic problems.     

Projectile acceleration in a single-stage gun at breech pressures below 50 MPa

Experimental studies were carried out to investigate projectile acceleration in a single-stage gun at breech pressures below 50 MPa. The gun was driven by firing either liquid or solid propellant. In-bore projectile velocity was continuously recorded using the well-known, precise VISAR interferometer technique so that accurate projectile acceleration data could be deduced. Both the attained projectile acceleration and muzzle exit velocity depend upon the charge-to-mass ratio and the pressure at which the blow-out disk ruptures. The results obtained from these experiments render information on the interplay between propellant combustion and projectile acceleration for low in-bore pressure regimes, and they provide the input data required for adequate numerical simulation.     

高速斜入水和水平入水气炮水箱实验系统

为开展模型高速斜入水和水中高速航行的水流场实验研究,研制了立式和卧式气炮与水箱组合的实验系统。通过快速阀和活塞阀控制气炮激发和驱动状态,一级气炮采用高压空气直接驱动弹托和模型,二级气炮采用高压空气驱动重活塞压缩使集气腔中产生高压气体,再驱动弹托和模型达到预定速度。通过调节水箱和发射管角度,使高速模型斜入水或水平入水。其中,立式可变发射角二级气炮可发射质量1～1000 g的模型至2500 m/s最大速度,卧式一级气炮可发射质量1～100 kg的模型至300 m/s最大速度。和小气室、高燃气压力火药驱动方式相比,新型气炮采用大体积、中低驱动压力气室,高压气体更接近等熵膨胀做功,调节高压气体压力,能较好地满足模型质量和速度的宽范围要求。结合光反射通断法测速、高速摄影和阴影流场显示等测量技术,得到立式气炮压缩管重活塞运动速度、压缩管末端压力时间曲线和模型倾斜与水平入水的流场阴影图像。结果表明:重活塞速度在膜片破裂前和理论计算值符合较好,但破膜后差异较大。立式气炮流场阴影图像反映了模型斜入水产生的空中和水中激波以及在气水界面的反射激波、空泡形成和侧向气水界面的破碎与飞溅等现象。从卧式气炮的模型水平入水阴影图像提取气泡轮廓,清楚地看出尾部气泡气水界面的波动和失稳。和商业计算软件Fluent计算结果相比,空泡上游区域基本重合,但尾流区域强湍流导致两者存在明显差异。和水洞实验相比,气炮水箱实验系统近真实地再现高速入水过程伴随的冲击和动态空化等物理现象和模型尺度效应。     

磁流体特性对磁流变火炮后坐阻尼器性能的影响

针对火炮后坐磁流变阻尼器的特点,以某单管25 mm火炮实验用磁流变阻尼器为研究对象,基于Herschel-Bulkley本构模型,建立了该炮磁流变后坐阻尼器的轴对称一维层流模型,获得了不同磁场作用下阻尼力随活塞速度的变化规律。运用火炮的后坐运动方程,计算了不同磁流体特性指数下,火炮的后坐位移和后坐速度。计算结果表明,磁流体特性指数的变化对磁流变后坐阻尼器的性能影响显著。     

The oblique-plate impact experiment

An experimental method is presented for the study of one-dimensional plane waves corresponding to combined pressure and shear. The experiment involves the impact of two skewed flat plates. A projectile plate is accelerated using a gas gun and made to impact a target plate in a vacuum chamber. The projectile and target plates are parallel, but inclined relative to the axis of the gun so that the particle velocity in the target has components both normal and parallel to the plane of impact. The particle velocity at the target rear (free) surface is recorded as a function of time. The normal velocity component is monitored using a laser velocity interferometer; the transverse motion is monitored using a shadow technique. The measured wave profiles can be compared to theoretical predictions based on one-dimensional-wave theory.