膜片非理想打开行为对自由活塞激波风洞运行的影响

2020-03-29张冰冰余亦甫

实验流体力学 2020年1期

朱浩, 张冰冰, 余亦甫

(中国航天空气动力技术研究院, 北京 100074)

0 引言

自由活塞激波风洞是获得高焓值高密度试验气流的重要地面设备。随着高超声速技术的发展，此类风洞日益受到重视[1-3]。重活塞压缩器，作为自由活塞激波风洞的核心部件，通过重活塞对驱动气体实施压缩，使之达到足够的温度和压力，从而获得强激波。在自由活塞激波风洞中，膜片打开是形成激波的先决条件。理论上一般假定膜片打开是一个瞬态过程，这使很多问题得以简化，但同时也掩盖了激波形成的复杂过程。基于试验观测，Heshall[4]对膜片打开后激波的形成进行了更为细致的刻画；Hickman[5]通过可视化技术观测到膜片打开时存在拟定常的自由射流；Outa等[6]发现了斜激波和多尺度膨胀波之间的相互作用。一些研究者还通过数值模拟勾勒出膜片打开过程中诸多难以观测到的流动细节[7-8]。日本NAL-CTA自由活塞发射机构采用CCD高速摄像机测量活塞速度，讨论了无量纲速度和压力恢复系数之间的关系，其中涉及膜片开口口径的影响[9]。

当前的研究重点关注主膜片“非理想”打开行为对活塞运动和激波形成的影响。“非理想”一词包含两层意思：一是膜片打开是一个非瞬态过程；二是膜片打开的最后几何形态也存在差异。本文基于膜片的剪切-应变模型和活塞动力学模型，得到了膜片打开-活塞运动的耦合方程组，分析了膜片非理想打开行为对活塞末端速度的影响规律，同时也研究了膜片非理想打开行为对激波形成和行进的影响。

1 膜片打开模型

在激波管实验中，膜片打开的具体情况在很大程度上依赖于膜片材质本身的延展性：膜片延展性越好，膜片打开得越规则[10]。在高压力气体作用下，初始缝隙出现于膜片痕道，并在痕道交点(膜片中心)形成穿孔；穿孔逐步扩大，形成若干膜瓣，最终膜片打开，形状接近球面(图1)。膜片初始打开过程比较缓慢，约占整个打开时间的20%[10]。厚度0.157～0.447 cm钢质膜片的打开时间约为800 μs左右[11]，该膜片厚度范围也为自由活塞激波风洞所采用。基于此，自由活塞激波风洞常压力驱动时间约为膜片打开时间的数倍，两者属于相同量级。

图1 膜片打开示意图

膜片打开是一个动态剪切-应变过程，且膜瓣周围的绕流流场十分复杂。假定膜瓣上的应力分布一致，则膜瓣打开角度θ可以表示为[11]：

(1)

式中：h为膜瓣宽度，ρd、w分别为膜片密度和厚度，pHA,r为膜片打开时驱动气体(氦/氩混合气体)的压力，σ为弯曲应力。在本文中，下标“0”和“r”分别表示活塞运动初始时刻、膜片打开时刻，下标“HA”和“A”分别表示驱动气体、空气，下标“p”表示活塞。

(2)

式中：时间t的单位为μs；p4表示驱动压力(平均意义下约等于pHA,r，可近似以pHA,r代替)。利用式(2)可知，当t=0时，膜片满足θ=0、dθ/dt=0，则膜瓣打开至θ=π/3位置约需400 μs，较为接近文献[11]提供的实际测量值。

2 膜片动态打开对活塞运动的影响

在压缩管中，重活塞在巨大的压力差下加速向下游运动，直至达到峰值速度Umax(在某个典型状态下，活塞峰值速度Umax接近300 m/s)，此时活塞前后压力相等；此后，被压缩气体压力进一步增高，活塞速度快速下降。活塞运动的相关研究参见文献[12-13]。以下讨论膜片打开这一动态过程对活塞压缩带来的影响，藉此完善对活塞压缩过程的刻画。

不失一般性，认为膜片面积等于激波管截面积，即S0=πh2/4；膜片形成孔的截面积S仅仅是膜瓣打开角度θ的函数，即S=S(θ)，可以得到：

S(θ)=S0(1-cosθ)2

(3)

驱动气体下泄方程和活塞运动方程可表示为[1,12]：

(4)

(5)

式中，mHA、ρHA、pHA和γHA分别为驱动气体的质量、密度、压力和比热比，Mp为活塞质量，x为活塞位移，pA为活塞后脸所受压力，D为压缩管直径，f为摩擦力。联立式(1)、(3)、(4)和(5)，就完成了对膜片打开过程对活塞运动影响的刻画。

为便于研究，引入无量纲时间τ=(t×aA,0)/D、无量纲位移ξ=(L-x)/D、活塞马赫数φ=up/aA,0和驱动气体下泄质量比μ=mHA/mHA,r，其中，up表示活塞沿x方向的运动速度，aA,0为活塞运动初始时刻空气的声速，L为压缩管长度，对式(1)、(3)、(4)和(5)进行无量纲化，可得到如下简化方程：

(6)

在上述方程组中：b3=(b1×pA,r)/pA,0，b4=b1b2，

若膜片为瞬态打开，在方程组(6)中令θ=0即可，这可视为方程组(6)的简化。但膜片打开动态过程对活塞末端运动的作用直接影响着自由活塞风洞的运行实践，而简化方程组无法对此进行刻画，方程组(6)的应用价值更加凸显。

3 数值结果

方程组(6)具有很强的非线性，不易获得解析解，需要通过数值方法加以研究。算例中，压缩管长度L为33.372 m，直径D为0.55 m，激波管直径为0.15 m，活塞质量Mp为200 kg；膜片材质为不锈钢，膜片厚度w为4.0 cm；驱动气体为氦/氩混合气体，体积比为7∶3，初始压力为19.0 kPa；高压空气压力0.65 MPa。

图2显示了膜片动态打开和瞬态打开情况下活塞速度和驱动气体压力2个重要参数的差异。从图中可以看到：在两种情况下，活塞均可以充分减速；但膜片动态打开过程在一定程度上减小了驱动气体下泄速率，对活塞减速有利，其代价是驱动气体压力(无量纲)上升更快，常压力驱动时间(驱动气体压力峰值下浮10%所需的时间)将减小。在该算例中，膜瓣从0°打开至75°用时接近1.2 ms。若膜片打开速度增加一倍，活塞减速将变缓(但并不明显)，而驱动气体无量纲压力峰值将明显降低，具体细节见图3。

实际运行中，由于膜片材料和痕道刻划的差异，即使在相同运行状态下，同一批次的膜片也未必能够精确地在某个固定压力下打开，而且膜片最终张开的几何形状也存在差别。膜片打开时间的提前或者滞后，对驱动气体状态和活塞运动都将产生影响。

基于方程组(6)，可以从理论上分析膜片提前(或者滞后)打开所带来的影响。图4比较了压缩比λ为30和40的情况下活塞速度和驱动气体压力的变化。相对于λ=40的情况而言，膜片在λ=30时打开可视为“提前”。从图4可以看到，当膜片提前打开时，活塞仍然能够充分减速，但减速距离稍长。膜片提前打开使得驱动气体压力增长相当明显而且陡峭(这意味着常压力驱动时间的减少)。理论计算进一步显示：在λ为30和40情况下，当膜片刚出现缝隙时，驱动气体压力分别为5.51和8.90 MPa；但随着活塞进一步压缩，两种情况下的驱动气体峰值压力相当接近，分别为11.55和11.82 MPa。结合图3和4，可以看出膜片动态打开过程将在一定程度上弱化膜片提前打开所带来的不利影响。

图2 膜片动态打开/瞬态打开情况下的活塞速度和驱动气体压力(λ=40)

Fig.2 Piston speed and pressure in two cases (λ=40)

图3 膜片不同打开速度(或时间)下的活塞速度和驱动气体压力(λ=40)

Fig.3 Piston speed and pressure under different rupture of diaphragm (λ=40)

图4 膜片正常/提前打开情况下的活塞速度和驱动气体压力(压缩比40、30)

Fig.4 Piston speed and driver gas pressure with different compression ratios

4 试验结果

在自由活塞激波风洞FD-21中开展相关试验，进一步研究膜片打开过程中非理想行为对激波形成和行进的影响规律。压缩管中布置的光电传感器可以实时准确地捕捉活塞轨迹和速度信息[14]。图5(a)对比了活塞速度理论预测值和风洞4个车次的试验测量值(4个车次的初始条件见表1)，两者之间相差小于5%，且活塞末端速度理论预测值更加贴合实际测量值。在4个车次中，活塞均可实现“软着陆”，这进一步证实了理论分析的结果，验证了方程组(6)的合理性。图5(b)展示了4个车次膜片打开后的情况，这也意味着膜片打开过程中非理想行为具有一定的频度。

表1 风洞试验初始条件及相关参数Table 1 Initial conditions and related parameters of wind tunnel tests

图5 FD-21风洞重活塞速度历程及膜片打开情况

由于FD-21风洞膜片夹持机构的结构限制，无法采用常规流场显示技术获取激波的形成和演化过程。在试验中，依靠激波管管壁布置的压力传感器数据计算激波马赫数，以此间接获取激波形成和行进的信息。换言之，激波马赫数差异是膜片打开非理想行为带来的直接结果。图6展示了4个车次中激波马赫数Ma随位置的变化。可以看到：膜片打开行为的差异对早期激波马赫数影响较大，但这一影响在激波行进过程中逐渐式微，在激波管末端，4个车次的激波马赫数差异不足5%。因此，在相同的运行参数下，膜片打开过程的差异不会对风洞贮室条件带来明显影响。早期激波马赫数的差异暗示了膜片打开行为的差异对激波形成的复杂影响。保持适当的激波管长度，可以有效抑制膜片非理想打开行为的不利影响。

图6显示激波马赫数峰值所在位置非常一致。由于压力传感器位置固定，无法准确区分不同车次的激波形成距离(激波马赫数达到峰值时所对应的位移)。一个初步的结论是：这4个车次的激波形成距离偏差不会超过2.5 m。