基于欧拉方法的2 维翼型冰晶结冰数值计算

针对现有航空发动机冰晶结冰试验方法存在数据采集困难,且难以揭示冰晶结冰机理问题,对某楔形翼型进行冰晶结 冰数值计算。采用S-A 湍流模型获得流场结果,应用欧拉方法获得冰晶和液滴轨迹结果,基于Messinger 模型获得冰形,并以 NASA-NRC 的第139 号试验结果验证了数值方法的可行性。结果表明：在其他工况相同时,压力越低（飞行高度越高）,冰晶/ 液滴蒸 发越强烈,结冰情况越严重;分析不同液态水含量与总水含量比例（CLW/CTW）对结冰的影响可知,液态水过少或过多都不利于冰层形成。     

3 m×2 m结冰风洞热流场品质提高及评估

结冰风洞热流场品质符合性是大型结冰风洞适航应用的基础.为明晰制冷系统性能升级优化对3 m×2 m结冰风洞热流场品质的影响,开展了热流场符合性验证试验,评估了热交换器出口和试验段两位置处热流场品质,给出了气流总温修正关系,形成了热流场控制包线.结果表明:热交换器出口和试验段模型区内热流场品质在主要试验工况下均优于SAE ...     

结冰实验段空气-水混合流动过程的数值模拟

针对结冰实验段内低温空气与高温雾化水滴的混合流动过程进行了数值模拟,采用拉格朗日方法建立了空气-游离水滴/冰粒多相混合流动与传热传质的一维分析模型,考虑了水滴粒径分布、相变及空气-水滴/冰粒传热传质等因素,求解了流道内气流温度/湿度/速度,水滴温度/速度/半径等参数沿流动方向的变化,分析了水滴粒径分布、气流温度、液态水含量、相对湿度等参数对流动与传热过程的影响,与已有数据的对比验证了模型及结论的正确性.      

结冰翼型表面明冰的压力分布和应力计算

在不同来流速度和攻角下,比较分析了冰型外部压力分布的变化规律,计算分析了外部压力分布对表面冰型的力学特性影响规律,研究了冰型和蒙皮界面间最大剪切应力的变化规律,并结合剪切强度随界面温度变化的实验曲线,分析了最大剪切应力对冰层脱落的贡献.计算结果表明:最大剪切应力是剪切强度的5.48%,因此气动力对冰层的脱落贡献不大;当开启热除冰装置后,由于蒙皮表面温度的不断上升,冰型和蒙皮界面间剪切强度的下降,当表面温度为-2℃时,最大剪切应力超过剪切强度的10%,此位置可能率先产生破裂,从而加速冰层的脱落.      

基于设计变量相关性的翼型设计方法

针对传统设计方法缺乏物理背景等不足,利用流动数值仿真及相关性分析的手段,对几何外形与流动变化的相关性规律进行经验总结,并将这种设计经验进行数值化描述后,引入到优化设计中,建立新的气动外形优化设计模型,使设计模型具有一定的物理意义,从而提高设计效率,改善设计结果。通过高速自然层流翼型HSNLF0213的优化设计,利用相关性的分析结果,建立具有物理意义的优化模型,开展了基于单点优化设计且同时满足多点设计要求的设计方法的应用研究,对这种设计方法及思路进行了初步的探索。设计结果证明了该方法的可行性。     

Gurney襟翼改善翼型动态失速特性研究

采用CFD数值方法,研究了NACA0012翼型、加装传统Gurney襟翼及改进不对称Gurney襟翼后翼型的动态失速特性。给出了传统襟翼对翼型动态失速特性的影响,并带来较大动态低头力矩的不足,基于传统襟翼的不足,提出了改进的不对称Gurney襟翼方案。研究表明,不对称Gurney襟翼可较好改善翼型的动态失速特性,在增加动态升力的同时,俯仰低头力矩明显减小,可能是直升机旋翼的较理想翼型。     

民用飞机二维增升装置设计

为满足民机市场对飞机高效、经济的需求,飞机在巡航飞行阶段必须具有飞行速度高、阻力小等特性,具体到气动设计方面,就要求飞机的机翼具有较大的翼载荷,即减小机翼面积以减小阻力,这给飞机起飞着陆阶段使用的增升装置的设计提出更高的要求。运用计算流体力学的方法进行了民用飞机二维增升装置,即多段翼型的设计,得到了良好的设计结果,并验证了在二维增升装置设计过程中采用前缘缝翼外形设计——带前缘缝翼进行后缘襟翼外形设计——前后缘增升装置缝道参数优化——起飞、着陆构型协调的设计路线图的可行性与有效性。     

二维翼型结冰数值计算

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

30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。