基于粒子轨迹的结冰风洞收缩段优化设计数值模拟

为了研究收缩段型面和液滴轨迹特性对结冰风洞流场品质与结冰试验效果的重要影响,获得收缩段的优化设计参数,采用数值模拟方法,使用二维简化模型对3m×2m结冰风洞收缩段的流场特性进行了研究,对几种典型收缩曲线的流场进行了计算,对比了收缩曲线对沿轴向压力分布、出口截面速度分布、不同粒径云雾粒子的极限释放距离及试验段入口云雾粒子分布的影响,结果表明:三次曲线与x_m为0.4的五次曲线不能满足流场合格指标,其他曲线流场特性差别不大,维氏曲线具有最小极限释放距离,在收缩段出口,三次曲线和五次曲线粒子包络面积最大,综合比较,x_m为0.5的五次曲线能有效兼顾流场均匀性、极限释放距离和收缩段出口粒子包络面积比.      

相对流速对非球形铝滴旋转受力的影响

针对固体火箭发动机中非球形铝滴的旋转受力问题进行数值模拟,详细分析流速变化对颗粒的受力影响。计算结果表明,非球颗粒在流场中的升力系数与无量纲转速之间已经不再满足经典的理论关系式。雷诺数Re1时,周向升力系数存在着交替出现的2组极大值与极小值;升力系数均为正值,且周向变化较为平缓。随来流速度增加,阻力系数变得均衡,系数值向0值逼近;升力系数也急速降低。当相对流速由0.5 m/s提高到5 m/s时,颗粒所受阻力锐降89%。颗粒的阻力系数与雷诺数之间满足新的关系式。     

多段翼构型的积冰数值模拟

基于非结构网格提出一种针对多段翼构型的积冰数值模拟方法.通过求解气体欧拉方程得到多段翼的绕流流场,基于空气绕流流场的气流速度分布,采用迭代法对水滴控制方程进行数值求解,得到多段翼各翼面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对某典型三段翼霜冰和明冰的结冰情况进行了预测研究,根据多段翼的流场特点提出一种明冰冻结模型的求解思路,并分析了积冰对多段翼性能的影响.      

双路离心喷嘴雾化特性的PDPA实验研究

采用相位多普勒粒子分析仪（PDPA）对某型航空发动机双路离心喷嘴的雾化特性进行了实验研究。PDPA可直接测得测点处的喷雾液滴的尺寸分布和速度大小,并据此求出了测点处的索特尔平均直径SMD和液滴的平均速度。在喷雾锥三个横截面上进行了测量,得到了SMD的空间分布,据此得到了喷雾锥的锥角,并与光学照相和计算机图像处理测得喷雾锥角进行了对比。实验结果表明：液滴尺寸随着供油压力的增大而减小,当压力增大到一定程度时,液滴尺寸趋于不变;当主、副油路分别单独工作时,随测量横截面与喷口之间距离Z的增加,SMD减小;在供油压力不变时,同一个测量横截面内,随着径向距离X的增加SMD值变化不大;喷雾锥角基本不随供油压力改变而变化。     

二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

考虑喷注流强分布的纵向稳定性建模与分析

喷注器流强分布是除声学阻尼装置以外使液体火箭发动机稳定工作的重要措施,通过控制喷注流强分布使推进剂的燃烧尽量远离主要振型的波腹区,减少热声耦合振荡的能量源,从而达到抑制高频振荡的效果。因此,建立喷注流强分布与稳定性之间的综合分析模型,研究流强分布的不稳定性抑制特性具有重要实际意义。针对采用自击式喷嘴器、液滴蒸发作为燃烧速率控制过程的某自燃推进剂缩尺燃烧室高频纵向燃烧不稳定问题,近似采用蒸发速率峰值区代替集中燃烧释热区,引入燃烧室三维声学控制方程以考虑多喷嘴条件下燃烧响应空间分布,建立了针对喷注器流强分布条件下的高频纵向燃烧稳定性分析模型,并对喷注流强的稳定性抑制特性进行了分析,给出了不同分布流强下燃烧室一阶纵向信号的增长率变化规律。研究表明,喷注流强分布有利于燃烧室稳定,"驼峰区"喷注孔径的增大对改善高频纵向不稳定性更为显著,"驼峰区"流强增加30%,相应的增长率降低15%。     

水滴撞击特性的粒子统计法研究

基于水滴撞击特性研究的拉格朗日方法,提出了适合于考察非定常情况的粒子统计法.通过统计非定常流场中不同时刻撞击在防冰表面上不同区域的粒子数,获得防冰表面上粒子浓度分布随时间的变化情况,与基于染色法水滴撞击试验获得的色度分布具有直接对应关系;采用基于染色法的动态图像测量方法对某型发动机帽罩缩比模型开展了水滴撞击特性试验,对所提粒子统计法进行了验证.基于拉格朗日方法的粒子统计法获得的计算结果与试验结果符合较好,与采用拉格朗日轨迹法获得的结果具有较好的一致性.粒子统计法所模拟的物理过程与试验过程一致,在处理湍流随机性和来流液态水含量分布不均等方面具有优势.     

风力机叶片翼型的结冰数值模拟研究

研究了风力机叶片翼型的结冰数值模拟方法,及翼型结冰后对其气动性能的影响。求解雷诺平均N-S方程,引入k-ωSST湍流模型封闭方程,获得风力机叶片流场;采用拉格朗日法计算风力机叶片翼型周围的水滴运动轨迹,同时考虑了多尺寸水滴分布的影响,获得翼型表面的局部水收集系数分布;根据质量守恒和能量守恒原理,计算翼型表面的各项热流,获得翼型表面的结冰速率和结冰冰形;考虑到翼型结冰是时间的动态函数,采用多时间步长法完成结冰数值模拟。计算了风力机叶片翼型在不同环境条件和气象条件下的结冰冰形,同时模拟了风力机翼型结冰后周围流场的变化,并与干净翼型的气动特性进行了对比。结果表明,环境温度较高时形成的明冰对翼型气动性能的影响较大,结冰导致翼型升力下降,升阻比减小,最大减小幅度达到61%,同时结冰后的翼型会提前进入失速区,导致桨叶气动性能恶化。     

翼型表面结冰准定常数值模拟

对翼型表面的结冰过程进行了准定常数值模拟。在每个时间步长内完成网格随着壁面边界的移动而更新、周围流场和水滴撞击特性重新计算、冰形计算及壁面边界的重构工作,如此循环直至所需的结冰计算时间。采用拉格朗日轨迹追踪法获得水滴撞击特性,应用考虑壁面粗糙度影响的边界层积分法计算壁面的对流换热系数,在此基础上求解Messinger结冰热力学模型,冰层始终沿壁面外法线方向生长。对NACA0012翼型在不同环境下生成的三种典型冰形进行了预测,并与实验结果进行了比较,表明本文所述方法的有效性。     

镁基水冲压发动机燃烧室长度及进水距离设计方法

针对两次进水的镁基水冲压发动机,对发动机总体的燃烧室长度和一次、二次进水距离设计方法进行了研究.根据理论分析和试验测温结果,提出金属镁以熔融液滴形式进入燃烧室,镁滴燃烧为蒸气相扩散燃烧过程.采用简化的D3/2定律计算镁滴燃烧时间,并以镁滴燃烧和雾化水滴蒸发所需时间为依据,建立了燃烧室长度和一次、二次进水距离设计方法.参...