带制退器的膛口射流噪声数值模拟与实验研究

为了研究膛口装置对膛口噪声气动特性的影响，对带膛口制退器的某小口径武器的膛口射流噪声进行了数值模拟和实验研究。采用计算流体力学CFD (computational fluid dynamics)-计算气动声学CAA (computational aeroacoustics)耦合算法对膛口噪声进行数值模拟，即对膛口流场进行瞬态CFD模拟，获取流场数据，然后利用所得到的结果采用声学方程模拟声源信息求解声场。基于数值模拟结果，分析了膛口流场变化及噪声的指向性分布，并与实验结果进行了对比。研究表明:膛口制退器的安装改变了膛口流场结构，影响了膛口射流噪声的指向性分布。计算结果与实验结果的误差小于9%，验证了该计算方法的可行性。研究结果可为膛口射流噪声的预测及膛口制退器的设计提供一定的参考。     

膛口流场的数值模拟

本文使用严格二维的FCT算法数值模拟美国M16 5.36毫米自动步枪的非定常的轴对称膛口流场。用此方法不仅可得到膛口冲击波的发展过程,而且可得到不少数值方法不易得到的马赫盘的发展过程。本文计算结果与Schmidt的实验结果作了比较,二者吻合较好。这一结果为进一步开展更复杂(考虑弹丸、膛口制退器等)的膛口流场的数值模拟奠定了基础。     

基于ALE方程的动网格膛口流场数值研究

基于ALE(Arbitrary Lagrangian-Eulerian Equation)方程的有限体积法,采用高精度Roe方法及结构化动网格,利用嵌入网格技术及动边界条件,对弹丸由膛内高压气体推动射出到完全飞高初始流场的整个过程进行了数值模拟.根据数值结果绘制了膛口流场密度和压力分布的时序图,形象地再现了膛口流场中初始激波、弓形激波及膛口冲击波,以及接触闻断、漩涡和剪切层等的动力学发展过程.结果基本反映了膛口流场的典型变化特征,为进一步研究真实膛口流场(如带化学反应、湍流等)奠定了基础.     

含化学反应膛口流场的无网格数值模拟

基于无网格方法,对包含大位移运动边界和非平衡化学反应的膛口流场进行了数值模拟。所发展算法是基于线性基函数最小二乘显式无网格方法,忽略黏性及湍流的影响,对流场采用ALE(arbitrary Lagrangian-Eulerian)形式的Euler方程描述,对流通量和化学反应源项采用多组分HLLC(Harten-Lax-van Leer-Contact)格式和有限速率反应模型计算,对于运动边界造成的点云畸形采用局部点云重构方法处理,重构过程中采用虚拟边阵面推进。对圆柱绕流和激波诱导燃烧流场进行了数值模拟,验证了重构方法和化学反应计算的有效性。最后对12.7 mm口径机枪膛口流场进行了模拟,结果同实验照片、非结构网格方法结果吻合较好,数值结果清晰地再现了膛口初始冲击波、膛口冲击波、欠膨胀射流波系结构的动力学发展过程,以及膛口焰的时间、空间分布特征。     

弹道枪不同水深下全淹没式发射膛口流场的数值分析

为了解弹道枪水下全淹没发射时，水深对膛口流场演化特性的影响，建立了二维轴对称非稳态膛口流场模型。采用流体体积函数多相流模型、标准k-ε湍流模型和Schnerr-Sauer空化模型，结合动网格及用户自定义函数技术，对水下全淹没发射膛口流场演变全过程进行了数值模拟。搭建了弹道枪水下可视化射击实验平台，对12.7 mm口径弹道枪在水中全淹没式发射时膛口流场演化过程进行了观测，并验证了数值模型的合理性。在此基础上，对比了不同水深下（h=1～100 m）膛口流场的演化特性。通过对比发现:在不同水深条件下，在膛口流场影响范围内，弹丸膛外行程随时间的变化均满足指数函数规律；水越深，膛口流场典型波系结构形成所需时间越长，且燃气在膛口轴向马赫盘处的温度和压力峰值越低，压力振荡幅度也越小，更快趋于平稳，但在径向上，水越深，压力振荡持续时间越长。     

膛口反应流并行数值模拟

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

内含运动物体的非定常流场计算

本文改造了高阶精度的MUSCL(MonotonicUpstream-centered Scheme for Conservation Laws)差分格式,对含有膛口装置和运动弹丸的火炮膛口流场进行了数值模拟,得到了由空气冲击波和射流近场结构所组成的全流场流动图象。     

火炮在不同介质中发射的膛口流场特性分析

为研究水下炮密封式发射膛口流场及在不同介质中的膛口流场分布特性,建立了水下密封式发射二维轴对称膛口多相流数值模型.采用VOF(volume of fluid)模型、标准k-ε湍流模型,结合用户自定义函数及动网格技术,分别对水下密封式发射与空气中发射膛口流场演化过程进行了数值模拟与对比.计算结果表明,火炮在水下发射时的膛口流场与空气中发射时有明显差异.水下密封式发射时的最大膛压与空气中基本相同,弹丸初速较空气中发射降低了32 m/s,而膛口压力与温度有明显的升高;水下密封式发射时大约在140μs初步形成马赫盘,而空气中发射时马赫盘形成较晚,约在320 μs;与空气中发射相比,水下发射时的激波核心区面积更小,且弹丸头部不存在冠状冲击波.水下密封式发射时,马赫盘距离膛口轴向位移随时间变化呈指数增长;空气中发射时,马赫盘距离膛口轴向位移随时间变化呈线性增长.     

含内插管扩张型抑制装置的膛口流场数值研究

采用ＴＶＤ格式求解Ｎ－Ｓ方程，数值模拟了含内插管扩张型抑制装置的膛口流场，展示了由膛口冲击波包围着的特征多变的射流结构的演变过程，模拟结果与实验照片相符。分析了膛口冲击波与射流结构的相互作用，揭示了锥状冲击波和锥状燃气团的形成机理。     

小口径膛口射流噪声的数值模拟

为了研究小口径武器的膛口气动噪声特性,采用CFD-CAA耦合算法对7.62 mm枪的射流噪声场进行了数值模拟。由于膛口流场结构复杂,在目前的计算发展水平下还不足以采用CAA直接法,因而本文中采用混合方法,即首先采用CFD方法计算7.62 mm枪的膛口流场,然后利用所得结果,采用声学方程计算射流噪声,具体为膛口近场采用LES进行计算,远场声场采用FW-H声拟法计算。通过对比验证实验,验证了该计算方法的可行性。然后,对7.62 mm枪射流噪声进行了数值模拟,分析了噪声指向性,绘制了声压级云图。研究表明:在本文的计算条件下,射流噪声强度主要集中在近膛口区域;且射流最大噪声主要分布在与轴线方向成30°~60°范围内。     

Application of an explicit TVD scheme for unsteady,axisymmetric, muzzle brake flow

Operator splitting in the presence of source terms is necessary in order to apply Harten's second-order accurate, total-variation-diminishing (TVD), shock-capturing scheme to higher-dimensional problems. By employing Godunov boundary treatment at the muzzle brake, such splitting is applied to the problem of muzzle brake flow simulation for the case of blasts which are impeded, by vertical and slanted baffles, respectively. Results from numerical studies of various types of wall boundary condition treatment which are consistent with Harten's TVD scheme are indicated.     

阻旋栅对梳齿密封动静特性影响研究

阻旋栅可改变密封进口流体周向流动与进口预旋，是提高系统稳定性的主要方法之一. 本文作者应用计算流体力学方法研究了阻旋栅几何参数对梳齿密封动静特性的影响，计算分析了阻旋栅在不同长度、间隙、周向个数及不同进口预旋比下密封流场分布与动力特性系数，并与无阻旋栅梳齿密封进行对比. 研究表明:阻旋栅能够有效抑制密封进口周向流动、降低密封腔室周向压力；随着阻旋栅周向个数与阻旋栅间隙的减小，其抑制效果增强，阻旋栅长度的增加对周向速度影响则越来越小；提高预旋比将使密封内流体周向速度增加. 与传统梳齿密封相比，具有阻旋栅的梳齿密封直接阻尼增加，交叉刚度降低，进而有效阻尼提高. 阻旋栅间隙s=0.20 mm、长度l=3.25 mm、数量n=90时密封有效阻尼较大，系统稳定性最好.      

Minimal models for squealing of railway block brakes

Block brakes have been used to brake railway vehicles for approximately 200 years. During the last 40 years, disk brakes have replaced the block brakes at passenger cars but block brakes are still used for all freight wagons. One problem of block brakes is that they show an enormous tendency to squeal. Although the block brake is a very common and old technical component, there exists almost no scientific work on its noise behavior regarding squeal. On the other hand, a lot of work has been done on the problem of disk brake squeal especially concerning the modeling of the excitation mechanism. The goal of this paper is to investigate whether and how models from disk brake squeal can be modified to model block brake squeal. The starting point of these investigations is a minimal model for an automotive disk brake introduced by von Wagner et al. (J Sound Vibration, 51(1–2):223–237, 2007) that is adapted to the block brake problem. It can be shown that such a simple converted model does not show any instability at all. A deeper analysis suggests that the reasons for squeal in block brakes could originate from in-plane vibrations of the brake disk or specific geometrical properties of the railway wheel. The self-excited vibrations explaining the squeal occur at relatively low rotational speeds far below the first critical rotor speed which has rarely been observed in rotor dynamics.     

NUMERICAL ANALYSIS OF ARC MOTION IN RAILGUN MUZZLE UNDER THE EFFECT OF BACKFLOW

在电磁轨道炮发射过程中,膛口外的高温气团以远高于声速的运动速度向内膛运动,这会使得膛口的电弧形态发生改变以及电弧位移。为了实现膛口电弧运动过程的数值计算,本文基于磁流体理论,引入了一种和电弧等离子体自身电导率密切相关的电位边界条件处理方法,建立了膛口位置电弧仿真模型。对回流效应下的电弧运动规律进行了分析,并且研究了运动过程中电弧对内膛流场的温度分布产生的影响。结果表明电弧的运动加剧了膛内温升。     

膛口流场一维简化非傅里叶传热模型

针对膛口流场复杂的传热现象提出一个简化的工程计算模型。根据该模型得出:在一定条件下，流场内存在非线性波传热模式，并把运动的膛口焰火球近似看成一个非线性波，所得计算结果与实验数据相符。     

化爆冲击波和炮口冲击波对生物致伤效应的对比研究

通过TNT化爆和火炮现场试验,探讨了化爆冲击波和炮口冲击波的异同点,结果发现两者的频谱大致相似。但化爆冲击波只有一次激波,而炮口冲击波有两次激波,有时甚至可见到三次激波。生物效应致伤的靶器官也有所不同,化爆冲击波以肺损伤较为多见,而炮口冲击波似乎以上呼吸道更为敏感,结果提示化爆可以模拟炮口冲击波,但又不完全等同于炮口冲击波,必要的火炮现场试验是必不可少的。