基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

超声速条件下亚毫米液滴的变形破碎模态

以超声速气流中液滴变形破碎行为为研究内容,对水平激波管中承受激波冲击的亚毫米水液滴(0.44~1.09 mm)变形破碎过程进行了观测,实验激波马赫数范围为1.07~2.11。利用纹影法,结合高分辨率高速相机对不同破碎模态下液滴的变形破碎特征进行了记录,得到了袋状、多模态、剪切和爆炸式等破碎模式下的液滴纹影图像,分析了液滴运动参数的时空关系。得出了液滴变形阶段,液滴无量纲横向变形宽度以及液滴无量纲迎风面位移随无量纲时间的变化发展规律,并且得出在液滴初始直径相同时,不同液滴破碎模式的无量纲最大横向变形宽度的变化,其中袋状、多模态、剪切破碎模式的无量纲横向最大变形宽度均在1.15~1.61范围内变化,爆炸式破碎模式的无量纲横向最大变形宽度均在0.21~0.68范围内变化。     

一种正向乘波前体设计方法初步研究

以定楔角乘波体设计方法为基础,研究了影响高超/超声速乘波体"乘波"的主要因素,给出了前体前缘实际气流压缩角的确定方法及影响因素,可知在相同的来流马赫数和压缩角δ下,随着前缘角θ和气流与前缘夹角α的增加,实际气流偏转角γ减小。据此,基于幂函数进气道前体构形,给出了前缘激波不脱体的限制条件及具体的判定方法,分析了乘波体典型几何特征参数对前缘激波不脱体的影响规律,结果显示在相同的来流马赫数和压缩角度下,增大前缘形状因子n,减小前体的长宽比L/W及增大前缘角均有利于激波不脱体。根据给出的前体几何参数对前缘激波脱体的影响规律曲线,对一种"前体几何外形构造+前缘激波附体条件限制"的正向前体乘波器工程设计方法进行了研究,给出了具体设计流程,并进行了初步的数值仿真验证,表明通过该方法设计的乘波前体流动特征与预期的结果吻合,说明文中所给出的激波附体条件及影响规律是可信的,乘波前体设计方法是可行的。     

Shock protection of penetrator-based instrumentation via a sublimation approach

It is often necessary for space-borne instrumentation to cope with substantial levels of shock acceleration both in the initial launch phase, as well as during entry, descent and landing in the case of planetary exploration. Current plans for a new generation of penetrator-based space missions will subject the associated on-board instrumentation to far greater levels of shock, and ways must therefore be found to either ruggedize or else protect any sensitive components during the impact phase. In this paper, we present an innovative method of shock protection that is suited for use in a number of planetary environments, based upon the temporary encapsulation of said components within a waxy solid which may then be sublimated to return the instrument back to its normal operation. We have tested this method experimentally using micromachined silicon suspensions under applied shock loads of up to 15,000g, and found that these were able to survive without incurring damage. Furthermore, quality factor measurements undertaken on these suspensions indicate that their mechanical performance remains unaffected by the encapsulation and subsequent sublimation process.     

Structure and kinetic properties of slow-mode shocks associated with magnetic reconnection in the near-Earth magnetotail

We study the structure and kinetic properties of slow-mode shocks near the plasma sheet boundary layer (PSBL) associated with magnetic reconnection by Cluster observation. The presence of slow-mode shocks is confirmed by traditional Rankine–Hugoniot (RH) analysis and Monte-Carlo shock fitting method. The Walén analysis, applied to the tailward flow associated with slow-mode shocks, also supports that plasma was accelerated across a Petschek-type slow-mode shock connected to the diffusion region. Back-streaming ions were observed on the shock layer, and cold ions were accelerated and heated by slow-mode shocks. In addition, whistler and electrostatic solitary waves were observed around the slow-mode shocks. These waves might be excited by the observed field-aligned electron beams near the shocks.     

基于剪裁技术考虑的某冲击试验分析

利用环境工程剪裁技术对某冲击试验的结果进行分析,结果表明,按照引用标准评定为"超差"的试验,当考虑到试验件结构的特性并进行必要的数据转换后,试验结果是有效。首先分析造成超差的原因为天线组件的剩余响应,然后通过频域数字滤波进行数据转换。     

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。