数值模拟方法对NASA CRM模型阻力预测的影响

在复杂工程外形的数值模拟中，网格类型、规模和分布、湍流模型、数值格式等都会不同程度影响模拟结果。如何评估这些模拟方法对模拟结果的影响，并识别对模拟结果影响较显著的因素，对于CFD的发展方向有积极的借鉴指导意义。为了综合研究不同因素对商用运输机外形阻力预测的影响，以AIAA第六届阻力预测会议外形NASA CRM为研究对象，同时考虑了网格、湍流模型、无黏通量格式和体心梯度求解方法等因素对阻力预测的影响。分析中采用枚举法和正交试验设计两种方法，并使用了基于聚类分析的定性敏感性分析方法和基于Mckay主影响分析的定量方法，识别出对阻力预测影响较大的因素，这为数值模拟方法的发展指明了方向。     

桨间距对某型共轴对转螺旋桨气动性能的影响

为掌握共轴对转螺旋桨的桨间气动干扰规律，降低桨间气动干扰强度，提升共轴对转螺旋桨的气动性能，基于非定常雷诺平均Navier-Stocks方程耦合湍流模型的计算方法，并使用了滑移网格技术，研究了4种不同桨间距的6×6构型对转螺旋桨的桨间气动干扰对其气动性能的影响，桨间距选择了4种不同方案。研究结果表明：在4种不同的轴向桨间距中，当桨间距为0.25倍螺旋桨直径时，共轴对转螺旋桨的平均推进效率最高，并且气动干扰导致的效率脉动幅度较小；随着桨间距的增大，前、后排桨受到的气动干扰强度都会减小，相比于后桨，前桨因气动干扰造成的脉动对桨间距更加敏感。可见共轴对转螺旋桨的桨间距会对螺旋桨的气动干扰，及气动性能产生较为明显的影响，在设计共轴对转螺旋桨时选择合适的桨间距，有利于提高螺旋桨气动性能。     

气液针栓式喷注器液膜破碎过程和结构的数值研究

为了更好地认识针栓式喷注器雾化场的结构，基于网格自适应加密技术以及VOF (volume of fraction)方法追踪气液的分界面，采用realizable k-ε湍流模型模拟整个流动过程，还原了不同时刻气/液撞击的初次破碎过程，数值模拟结果与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。进一步对针栓式喷注器气/液撞击的初次破碎过程、内部流场涡结构、速度场进行分析，研究了初次破碎雾化的动力学过程和机理。研究结果表明：液桥的形成主要是由液洞的扩展和拉伸、合并而形成，而液滴主要是由中心液膜拉伸、液丝断裂以及液桥断裂而形成，液膜破碎阶段形成的涡结构是造成液膜断裂的主要原因。     

计算流体力学非结构混合网格高精度格式研究现状

介绍了航空CFD技术中的非结构混合网格高精度格式的基本理论和方法,国内外研究现状,几种非结构混合网格高精度格式。     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

基于Fibonacci搜索方法的航空发动机性能寻优

介绍了Fibonacci搜索方法应用于发动机性能优化中的原理,建立带两个出口面积、气流马赫数可变的引射器的涡喷发动机的总体性能计算模型,应用Fibonacci搜索方法在优化参数区间内对该发动机进行了最大推力、最低油耗、最低涡轮前温度模式的寻优.举例说明了Fibonacci搜索方法应用的具体过程.给出了不同高度、马赫数下的寻优结果.结果表明:最优目标性能满足的约束边界与飞行高度、马赫数密切相关.不同的飞行高度、马赫数对应了不同的最优控制策略.但不论飞行任务如何,根据该方法进行的优化控制所获得的目标性能是最佳的.该优化方法不仅适用于该类涡喷发动机,也为其他类型航空发动机的优化控制打下一定的基础.      

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

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

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

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

外挂物投放问题的二维非结构动态网格技术研究

以动态网格守恒型方程为控制方程,在二维非结构网格基础上发展了弹簧式光顺法和局部网格重构法两种动态网格变形技术。结合运动轨迹计算方程,将动态网格法应用于飞机外挂物投放问题中,成功模拟了导弹自由投放的整个过程,获得了详细的包括弹体运动速度、运动轨迹和受力情况等在内的丰富的流场信息。成功的应用不仅验证了动态网格法的可行性;而且也表明,动态网格法的进一步完善与发展将成为解决外挂物投放问题的有力工具。     

关于变形网格“几何守恒律”概念的讨论

在国内外文献中应用包含有变形的动网格进行流动数值模拟时,要满足所谓的几何守恒律条件;通过对这一概念最早提出推导过程和近年来应用发展情况分析以后发现,逻辑上存在不适当之处.从流体力学基本理论出发,对包含有变形的动网格技术的理论基础进行讨论,证实几何守恒律是流体力学控制方程的伴随方程或退化方程;通过简单模型分析有限体积方法离散过程,发现计算过程中不满足几何守恒律所引起的非物理现象本质上是目前有限体积离散处理中微元体积计算方法不符合物理定律所致;因此在变形动网格计算方法中不存在必须要满足的几何守恒律.最后根据理论提出物理背景清晰的微元体积计算方法,数值验证可行.