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

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和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流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

可重复使用热防护材料应用与研究进展

可重复使用热防护系统是为高速重复使用飞行器而发展的关键性技术,涵盖了地球大气环境及非地球大气环境下的弹道式再入、高马赫数巡航等应用场景。根据现有高马赫数飞行器热防护现状,对高马赫数飞行器的主要热防护系统类型、特点和使用场景进行了简要介绍。在此基础上,结合国外里程碑式可重复使用飞行器（X-15、SR-71、航天飞机、X-33、X-37B、Spaceliner等）,梳理了可重复使用热防护材料的应用与研究进展,论述了代表性可重复使用热防护材料的发展、性能、研制进度、特点及应用前景。对国外在可重复使用热防护材料研制中的设计及发展思路,以及所存在的主要问题进行了总结归纳,为可重复使用热防护材料未来的发展提供了思路。     

基于PIV技术的单圆孔脉冲射流流场特征

对稳态射流及脉冲射流冲击靶板时的流场特性结构进行了探索和分析。采用高频粒子图像测速技术,在射流管口到冲击靶板间距为6倍管径的条件下,对稳态射流进口雷诺数为6 000的稳态射流及脉冲频率为20 Hz的脉冲射流进行了实验测量,得到了射流核心区、壁面射流区及滞止区内的速度分布。研究发现:①由于射流剪切作用的影响,脉冲射流核心区的最大轴向脉动速度为稳态射流的3倍。②滞止区内,由于射流的剪切作用和壁面的滞止作用,导致了脉冲射流轴向速度梯度最大为稳态射流的2倍,同时,滞止区内的最大脉动速度是稳态射流脉动速度的3倍。③脉冲射流对壁面的卷吸以及旋涡的产生和传播过程,破坏了壁面射流区稳定的速度边界层。相比稳态射流,脉冲射流的流场增加了湍流相干结构的含能并产生周期性的大尺度卷吸涡。     

Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian (ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces. In addition, the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method. Based on this method, fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained. Numerical results such as maximum Root Mean Square (RMS) drag, general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results. This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.     

Sensitivity analysis of flowfield modeling parameters upon the flow structure and aerodynamics of an opposing jet over a hypersonic blunt body

Implementation of an opposing jet in design of a hypersonic blunt body significantly modifies the external flowfield and yields a considerable reduction in the aerodynamic drag. This study aims to investigate the effects of flowfield modeling parameters of injection and freestream on the flow structure and aerodynamics of a blunt body with an opposing jet in hypersonic flow. Reynolds-Averaged Navier-Stokes (RANS) equations with a Shear Stress Transport (SST) turbulence model are employed to simulate the intricate jet flow interaction. Through utilizing a Non-Intrusive Polynomial Chaos (NIPC) method to construct surrogates, a functional relation is established between input modeling parameters and output flowfield and aerodynamic quantities in concern. Sobol indices in sensitivity analysis are introduced to represent the relative contribution of each parameter. It is found that variations in modeling parameters produce large variations in the flow structure and aerodynamics. The jet-to-freestream total-pressure ratio, jet Mach number, and freestream Mach number are the major contributors to variation in surface pressure, demonstrating an evident location-dependent behavior. The penetration length of injection, reattachment angle of the shear layer, and aerodynamic drag are also most sensitive to the three crucial parameters above. In comparison, the contributions of freestream temperature, freestream density, and jet total temperature are nearly negligible.     

整体弹性结构法及在叶片流固耦合分析中的应用

针对三维叶片时域流固耦合振动响应计算普遍耗时的问题,采用一种假设整体结构的模态位移,求解固体时域响应和实现高效网格变形,发展了一种3维时域流固耦合分析的整体弹性结构方法,并应用于压气机叶片0°和180°相位角的气动稳定性分析中。结果表明,所发展方法的计算结果与传统双向时域算法和文献的结果较为吻合,而计算效率相比于传统算法显著提升;在所分析的两个相位角下,叶片振动的气动阻尼均随流量减小先增大后降低,相比于0°相位角,180°下叶片的气动稳定性大幅提高,表明该方法能有效应用于叶片的工程流固耦合研究。     

基于拓扑优化和形状优化的桨叶结构设计

直升机旋翼桨叶剖面结构设计是研究旋翼动力学设计的基础。提出一种联合拓扑优化和形状优化的 两级优化设计方法,可用于桨叶剖面结构的设计。第一级优化以变密度法(SIMP)作为拓扑优化的材料插值方 法,以桨叶体积最小化为设计目标,约束桨叶节点位移和应力,建立桨叶的拓扑优化求解模型;第二级优化以重 构的桨叶模型为基础开展形状优化,降低局部应力集中以及找到合理的边界节点位置。对优化后的模型进行 有限元分析,结果表明:通过拓扑优化和形状优化的两级优化,能够得到满足强度和稳定性要求的结构布局,为 桨叶结构设计提供指导方案。