重活塞驱动立式气炮流场的参数计算.docVIP

重活塞驱动立式气炮流场的参数计算* 张伟伟 谢武俊 徐胜利 黄生洪 （中国科学技术大学力学与机械工程系 合肥 230026） 摘 要 本文目的是认识重活塞驱动气炮流场特征和影响性能主要因素。理论估算基于等熵和非定常气体动力学关系式，弹丸速度估算采用等底压假设。采用商业软件FLUENT，计算重活塞驱动立式气炮三维流场，得到活塞速度随压缩管长、弹丸出口速度随时间变化以及压缩终态的集气腔流场参数分布。结果表明：活塞最大速度估算值和计算值偏差小于20％。活塞加、减速近似为线性规律。对应不同破膜压力，活塞最大速度点位于管长3/4处。提高破膜压力，弹丸速度提高幅值有限。弹丸初期加速规律近似满足等底压假设，但后期偏差较大。被驱动气体压缩终态近似为准静态，集气腔流场三维效应不明显。本文条件下，弹丸出口速度大于1.5km/s，满足设计要求。 关键词：自由活塞，气炮，发射管，数值模拟 1. 引 言 以重活塞为驱动器的气炮在高速碰撞实验研究中有广泛应用，其核心为自由活塞驱动技术，简称重活塞驱动器（piston driver）。目前，国内外气炮或弹道靶都采用水平卧式布置。为研究高超声速弹丸斜入水现象，需建造立式布置气炮。液面受重力作用始终处于水平位置，卧式气炮不适合研究弹丸和液面斜撞击问题。Ray Stalker教授1966年提出自由活塞驱动技术，已成为高超声速空气动力学和高速碰撞研究的常用驱动器。澳大利亚国立大学应用物理系（T3）和昆士兰大学高超声速气体动力学中心（T4，X1-X3）[1]、日本NAL Kakuda研究中心（HIEST）[2]、德国哥廷根超高速空气动力学研究所（HEG）、美国加州理工学院航天工程系（T5）[3]、俄罗斯科学院西伯利亚分院理论和应用力学研究所（ITAM, RAS）均建造活塞驱动器的激波风洞和膨胀管。研究表明：该驱动器可提供高质量、物理压缩的高温高压气体。国内还没有研制活塞驱动器的高焓和高速碰撞设备。 气炮性能参数计算包括压缩管和膨胀管，分别依据平衡热力学和非定常空气动力学理论。对水平卧式的轴对称气炮结构可得到近似解。本文气炮布置为“L”型结构，与轴对称结构相比，流场存在较强三维效应。首先，本文用等熵和非定常气体动力学方法估算活塞和弹丸速度。然后利用FLUENT商业软件求解非定常Navier-Stokes方程，得到气炮流场数值解。对比两者结果，认识求解精度的差别，为设备调试和状态参数确定奠定基础。 2 设备结构简图和主要参数 图1给出了设备简图。压缩管轴线离地面5m，设备呈L型结构。压缩管右端采用配重活塞，为减小重活塞运动导致管体后退（recoil），集气腔和压缩管采用滑动配合，以消除发射管和压缩管相互影响。利用压缩管和发射管面积变化，可获得较大驱动增益。主要设计参数：压缩管长6m、内径 125mm。发射管长为4－6m、内径40mm。弹丸出口速度大于1.5km/s。 重活塞质量40kg。弹丸和弹托总质量为40g（弹丸质量为20g）。破膜压力大于15.0Mpa。驱动气体（He）终态压力、温度大于15.0Mpa和2000K。 图1 立式气炮示意图 1配重活塞 2进气腔 3重活塞 4集气腔 5缓冲活塞 6夹膜机构 7弹托捕捉器 8 实验舱 9观察窗 10 弹丸 11高压气罐 12压缩管 13发射管 14 水面 3 气炮参数理论计算 图1给出了气炮简图。气炮工作过程为：高压气体r驱动活塞，管内预先充入初态i的气体被压缩至终态m。此时，膜片瞬时破裂，终态m气体膨胀并推动弹丸运动，弹丸底压状态为ms。首先考虑压缩管气体计算。即找出驱动活塞高压气体和被驱动气体初、终态之间的关系。假设活塞速度低于气体声速、压缩终态温度较低（3000K），忽略活塞、壁面在压缩过程中的热传导和摩擦效应。认为被驱动气体由I到m为等熵压缩，并处于热力学平衡，因此，可利用平衡热力学关系式，得到如下关系式 ， （1） 其中，为比热比，P、V分别对应压力和体积。脚标m,i,r对应被驱动气体初、终态和高压气体状态。 图2 气炮参数理论估算示意图 （1）发射管计算：确定破膜后弹丸底压状态m。分析模型为一维变截面管道（见图2），当选定破膜压力，考虑到变截驱动增益，即变截面处的定常、非定常流膨胀，利用非定常波系的等熵关系，可确定弹丸底部状态m和驱动气体状态i关系。 图3 压缩管中Pm/Pi和Tm/Ti和Pr/Pi关系曲线 图4 重活塞速度和运动时间随压缩管长度变化 （Xi=6m，Pr/Pi=9，Pi=0.132MPa，氦气驱动） 图3给出了驱动气体终态压力、温度随高压气体压力的变化。当PR相同，不同比热比对应状态m参数差别较大。选择较大（1.67）He作为驱动气体，较小PR得