过冷大水滴动力学行为对翼型结冰的影响分析

过冷大水滴（SLD）是极端危险的飞行环境之一。因大粒径水滴独特的动力学行为变形破碎、飞溅反弹,传统结冰计算方法难以准确地反映SLD结冰情况。采用Navier-Stokes方法求解流场、Euler方法计算水滴撞击、Shallow Water模型模拟结冰,并与NASA实验结果进行了对比验证方法可信。结果表明：SLD动力学行为对结冰和冰形影响较大。其中,变形破碎改变了水滴运动轨迹和撞击范围,降低了水滴撞击极限,导致上下结冰极限减小2.83%、2.13%;飞溅降低了驻点附近水滴收集率,导致前缘积冰量减少8.09%;反弹显著降低了水滴撞击极限,导致上下结冰极限减小30.69%、20.01%;撞击后飞溅反弹二次水滴再入流场使得上下结冰极限增加6.14%、3.71%。同时,与干净翼型相比带冰翼型空气动力学性能严重退化,在相同迎角下,升力更小、阻力更大、气动效率更低。     

密集型空气雾化流场破碎特征数学模型研究

低速液体射流在高速湍流气体作用下的气液同轴射流雾化流场是瞬态密集型喷雾场,高速湍流中影响和控制雾化的因素很多。针对空气雾化流场液核分布与特征进行分析,开发了模拟液核的随机浸入体模型,结合大涡模拟方法对同轴射流空气雾化喷嘴下游流场进行数值模拟。模拟结果与实验结果对比表明：随机浸入体模型可快速捕捉液核的长度和位置信息。当气液动量比M为3~10000时,能够较准确地预测液核长度,当M>10时,其预测结果远优于唯象模型;该模型能够捕捉回流区、大尺度涡等流场结构。同时,可以准确地预测喷嘴附近的液滴粒径,特别是在气流速度较大（>60m/s）时,液滴平均直径预测误差<10%。     

液体圆柱射流在气流中的破碎特性实验研究

本文采用高速相机对低速横向气流作用下的圆柱射流表面波发展及液柱断裂和破碎进行观察研究。实验喷嘴为直射式，孔径为1mm，长径比为20。工质采用水和空气；工况为:温度293K，液体射流速度为2~20m/s，雷诺数为2400~22400，横向气流速度为10~40m/s，气流韦伯数为1.6~25.6，液气动量比为5~127。高速相机帧幅为2000，曝光时间为16s。通过实验观察到横向气流气体韦伯数的变化导致射流破碎形式呈现不同形式变化，液体射流的无量纲表面波波长与气流韦伯数的-0.31幂指数方成正比；主液柱断裂点沿横向气流方向的距离随着液气动量比的增大而减小，而沿初始射流方向的距离随液气动量比的增大而增大；断裂后产生的液滴在沿横向气流方向的速度分量为横向气流速度的0.1倍左右，而沿初始液体射流方向的速度分量先呈现出与液气动量比线性增长关系，直到其变为射流初始速度的0.8倍左右并保持在这一水平。在上述研究基础上，本文拟合了低速射流表面波的波长与气流韦伯数间关系式以及射流破碎位置、射流轨迹及液柱断裂产生液滴的速度与射流初始条件间的数学关系。     

剪切气流中无黏液体横向射流破碎机理

航空发动机燃烧室中,液体燃料的雾化过程、特别是初始雾化过程非常复杂,至今无法建立准确的初始雾化模型。忽略液体射流黏性,采用线性不稳定性分析法对无黏液体（工质为水）在二维剪切横向气流中的破碎机理进行研究。通过建立射流的色散方程,根据其表面波的增长率及波数的发展情况对射流破碎进行预测。当气体韦伯数或液体韦伯数增大时,射流表面波增长率显著增加,最佳波长明显减小。液气动量比大于临界值时,Kelvin-Helmholtz （K-H）不稳定性占主导作用;反之,Rayleigh-Taylor （R-T）不稳定性占主导作用。二维剪切气流在射流方向上具有速度梯度,只改变横向气动力对射流表面波的作用。气体韦伯数与液体韦伯数对射流破碎的作用与均匀气流相似;保持气流流量相同时,负梯度剪切气流可以加速射流的破碎。     

小型涡扇发动机燃油控制规律

以小型涡扇发动机DGEN380燃油控制规律为研究目标,使用解析法基于MATLAB平台建立其整机部件级稳态及动态数学模型并进行仿真计算。在旋转部件建模过程中使用工具获取特性数据,并绘制可以即时显示当前工作点/线的部件特性图;在稳态仿真时使用Newton-Raphson法和遗传算法2种方法求解发动机工作方程组,并对二者进行评估对比,针对遗传算法在发动机模型求解过程中的局限性对一般基础模型进行优化改进;在动态仿真时,采用欧拉法计算微分项。基于模型分析发动机燃油控制规律,并以发动机高度特性为例进行仿真计算。试验结果表明:Newton-Raphson法具有高精度和高速度,改进模型后遗传算法具有更好的适应性,模型及燃油控制规律的仿真结果较好地符合试验数据和理论趋势,对比误差小于3%。      

V型尾缘喷管流场特性数值模拟

模拟计算了一系列V型尾缘(Chevron)收敛喷管的详细涡结构(包括流向涡和展向涡)、引射系数、推力损失等性能,并与相关实验结果进行了对比验证。采用有限体积QUC IK格式和标准k-ε湍流模型,结合壁面函数求解雷诺平均的N-S方程。通过模拟计算与实验结果的对比,验证了计算方法的可行性。通过计算得出:V型尾缘喷管齿弯角越大,核心区长度越短,喷管出口下游的流向涡和展向涡的强度和尺度也越大。V型尾缘的齿弯角每增加1°,对应的喷管的堵塞比约增加0.98%,喷管的推力损失约增加1.71%,引射系数约增加3.45%。     

基于统计学特征的密集型气液雾化随机模型优化方法研究

在高气液动量比的空气雾化流场中，液滴和液丝从液核剥离的过程具有高自由度、分布密集的特点，传统理论模型难以对其准确预测。本文对结合大涡模拟方法的随机雾化模型进行优化，在初始雾化过程，提出使用液滴统计平均温度来表征液滴碰撞统计学特性的改进方法。液滴的统计平均温度分别采用气液相对动能模型的粒子追踪法和亚网格动能模型的粒子追踪法。研究表明，使用改进的气液雾化随机模型预测密集型空气雾化流场，大幅改善了传统雾化随机模型在初始雾化区域过预测的缺陷，平均动能的相对误差为15.5%，平均索特尔直径的相对误差为7.2%，与未改进前的模拟结果相比，误差降低了41.1%和15.0%。此外，本文还探究了喷雾张角模型对雾化流场预测结果的影响，分别将实验所得经验公式法、相界面气液动量平衡所得模拟法与亚网格动能模型的粒子追踪法结合。结果表明喷雾张角经验公式预测结果更为准确，在平均索特尔直径预测方面准确性提高了17.3%。     

黏性液体横向射流破碎机理

航空发动机中液体燃料雾化过程十分复杂，特别是雾化初始阶段，至今无法建立准确的雾化模型。本文利用线性不稳定性分析法研究均匀气场和二维剪切气场中，不同黏性液体横向射流破碎过程。在圆柱坐标系下，建立有黏液体横向射流色散方程，利用Muller法求解得到射流表面表面波的不稳定增长率随波数的变化情况。当来流为均匀气流时，考虑液体黏性影响，射流表面扰动波增长率的减小量与增长率的比值较不考虑黏性时的至少增大1 000倍。黏性对Kelvin-Helmholtz (K-H)和Rayleigh-Taylor (R-T)不稳定性均起到削弱作用，抑制了射流破碎，且表面波波数越多，黏性对破碎的抑制作用越强，但并不影响射流的不稳定波数范围。当黏性系数增大500倍时，表面波最大增长率降低80.37%,最佳波数减小40%,黏性力对表面波的抑制作用十分明显。横向来流为二维剪切气流时，横向气动力和剪切速度促进液体射流表面波生成进而产生破碎，而液体的表面张力和黏性力则会抑制表面波生成。液气动能比越大，气流剪切作用对射流不稳定性的影响越大，K-H不稳定性主导射流表面波的生成。进一步研究液体黏性对射流不稳定性的抑制作用，发现黏性...