基于晶体塑性理论的疲劳裂纹起始数值模拟

在疲劳荷载作用早期,材料在晶体尺度出现裂纹的萌生,为了研究疲劳短裂纹的发展,利用Monte Carlo法建立了多晶体晶粒集合的Voronoi有限元模型,并在Ansys Usermat子程序接口下编写了晶体塑性本构方程子程序,修正了以拉伸硬化为主的疲劳裂纹起裂的计算机模拟方法,并结合TANAKA和MURA的位错偶极子模型,模拟了晶粒集合在疲劳荷载作用下的裂纹萌生,最后与文献实验数据对比,这一修正的方法更加符合实验结果和宏观现象.     

基于OpenFOAM的真实流体模型库及均相求解器开发

为了正确认识高压环境下的流体非理想性特征,及其对于低温射流喷雾过程的影响作用,将多种不同真实流体模型写入到开源CFD程序OpenFOAM仿真平台中,编译形成具有广泛通用性的真实流体热物理模型库;并将高压下真实流体特征的影响引入到压力修正算法中;开发出适用于描述超临界环境下低温流体流动过程的均相求解器。对高压下流体热物理性质的非理想性,以及超临界环境下的低温液氮射流喷注过程进行研究。结果表明:高压下的真实流体效应主要体现在低温区域,本文所建立的模型库可以在任意热力学状态范围内正确预测流体热物性;基于此模型库所开发的均相求解器可正确预测超临界环境下低温射流的详细流动特征和结构,而理想气体模型则会大幅低估射流核心区密度分布,中心线上误差高达68.9%;综合考虑计算精度和效率,PR EoS及相关模型更适合于这类超临界流动问题的求解;低温射流跨临界喷注中局部区域所发生的虚假沸腾现象会导致射流温升速率降低、体积膨胀速率加快;射流核心区压缩因子明显1,在数值仿真时必须考虑高压环境下的真实流体效应。     

预燃式等离子体射流点火器工作特性实验研究

为解决发动机点火包线小于飞行包线的实际问题,提供等离子体点火技术在航空发动机中的工程应用新思路,论文设计了一种预燃式等离子体射流点火器,实验研究了放电特性和射流特性。结果表明,预燃式等离子体射流点火器与空气等离子体射流点火器相比,在提升射流能量降低电源功率方面有着较大的优势,电流相同时通入甲烷在较大流量时可减小驱动电源功率,总流量为44L/min时,减幅可达14.99%;同时预燃式等离子体射流较空气等离子体射流稳定,且射流长度增加,扩大了点火面积,有利于点火。     

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

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

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

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

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.     

高超声速湍流高效模拟算法

采用了直接过滤的Navier-Stokes(N-S)方程组对高速可压缩湍流进行研究。针对高超声速湍流的非线性流动特性,对N-S方程直接过滤推导了大尺度湍流流场的控制方程,更精确地反映高速湍流的可压缩性,建立了可压缩湍流的大涡模拟TDM模型。使用传统的Smagorinsky模型的非线性推广,采用基于非Favre过滤的超声速可压缩湍流的大涡模拟模型,应用Caylay-Hamilton定理,建立可压缩湍流大涡模拟的非线性亚格子模型,并发展为动力学模式,模型中的两个常数通过当地流场动态的确定,消除了可调经验常数的影响。针对构造的高超声速湍流大涡模型开发相应的高效并行算法。