阻力精确分解的方法研究与数值模拟

根据阻力产生的物理机理将阻力分解为近场阻力和中场阻力,建立实现了阻力精确分解方法:近场阻力法和中场阻力法,前者对飞行器表面压强和切应力积分将阻力分解为压差阻力和摩擦阻力,后者通过区域划分式积分将阻力分解为波阻、粘阻、诱导阻力和数值耗散阻力.选取RAE2822翼型、M6机翼和某宽体标模进行数值模拟,验证阻力分解方法的正确性,对比两种方法进行阻力分解辨识的能力与不足.研究结果表明,中场阻力法可以给出详细的阻力构成,更有利于进行阻力产生分析和减阻优化设计;近场阻力和中场阻力计算可以很好地达到阻力平衡,但由于消除了数值耗散阻力,中场阻力法的计算结果要比近场阻力法更接近实验值.     

基于尾迹积分的阻力计算方法研究

利用尾迹积分法按阻力产生的物理机制,将阻力分解为熵增阻力和诱导阻力进行数值计算.翼型、机翼和翼身组合体标模算例表明,尾迹积分法可以清楚给出由粘性和激波产生的熵增阻力及诱导阻力的展向分布.熵增阻力计算时,尾迹积分截面要尽量靠近后缘,若尾迹积分截面向下游移动,由于数值耗散使翼尖涡动能转化为内能从而将部分诱导阻力转化为熵增阻力,使诱导阻力值偏小,需要进行熵增修正,修正后的总阻力与实验值吻合更好.与壁面积分法的对比分析表明,壁面积分法计算阻力误差大部分来源于压差阻力.利用尾迹积分计算了升力,与壁面积分法相比,给出的升力线斜率与实验值更接近.     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

冷却系统阻力对汽车总阻力影响的数值预测模型研究

采用CFD方法数值模拟了小型汽车发动机冷却系统对全车阻力的影响特性.汽车的内部及外部(发动机舱及舱内)均设置了数值模拟网格.分别计算了散热器及风扇对阻力的影响,采用本校实验结果(散热器压力损失)及制造商提供的数据(散热器及风扇的质量流率)确定初始及边界条件,为基于RNS的CFD提供输入.对汽车轮胎运动模式及行驶速度对散热系统阻力及汽车总阻的影响特性也进行了计算分析.预测结果表明,冷却系统阻力占汽车总阻的6.9％,与实测值吻合.     

战术导弹超音速零升阻力计算

本文介绍了一种导弹类型超音速复杂组合体零升阻力计算方法。此方法先用面积律求取干扰因子和用二级理论计算主弹体的波阻,再用积分关系式解法求取湍流边界层位移厚度,并用它对波阻进行修正。摩阻、底阻和附加物的阻力计算则采用经验方法。计算结果表明,此方法有相当好的准确性,而且适用范围广。     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

基于 CFD 和推阻分解技术的全机溢流阻力预测与分析

飞机在飞行过程中,根据性能需求,需要不断调整发动机活门流量系数,因此发动机短舱唇口的压力分布形态会发生很大变化,前方外流作用于进气道内流管上的合力也将改变,从而引起溢流阻力变化。本文基于某型号飞机,结合 CFD 动力模拟和推阻分解方法,获得不同流量系数下的溢流阻力,并分析流量系数、马赫数、高度、迎角对溢流阻力的影响。溢流阻力预测方法和影响研究可为飞机/涡扇发动机一体化设计、大涵道比短舱设计和气动力预测提供参考。     

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

复杂外形气动力计算的阻力修正方法

基于气动力计算的位流线化理论,发展了一种亚音速附面层阻力修正方法。采用联合流场求解方法获得复杂外形的表面压力分布,在此基础上,将绕复杂外形物体沿流向面元划分为片条,在片条上,把三维附面层转换为二维附面层,对二维附面层进行数值求解,得到附面层粘性阻力修正。算例表明,该方法计算量小,与实验和N-S数值解结果吻合较好。为复杂气动外形飞行器初期设计、方案选型及评估提供了较为理想的工程设计手段。     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

DLR-F4翼身组合体的阻力计算

为了考察自行研发的CFD软件的计算能力和阻力计算精度，本文采用LU—SGS方法、MUSCL差分格式和Baldwin—Lomax代数湍流模型，数值模拟了AIAA阻力计算工作室提供的DLR—F4翼身组合体的绕流流场，综合分析了easel和case2的气动力的计算结果，并与NASA Christopher L．Rumsey采用CFL3D6．0和AFRL/VAAC Don W．Kinsey采用Cobalt60提供的两组计算结果以及AGARD提供的两种不同风洞的测力试验结果作了比较。计算结果表明，本文计算精度与国外CFD软件相当。为了提高激波，边界层干扰的模拟精度，今后要重点加强湍流模型的应用研究。     

DLR-F6翼身组合体数值计算

使用自行研制的混合网格亚跨超声速流场解算器程序MFlow计算了AIAA第三届阻力会议提供的DLR-F6及其带整流装置FX2B情况下的阻力。重点分析了两者的网格收敛特性、阻力极曲线以及压力分布等,并与阻力会议的各个软件的计算结果进行比较。详细分析了DLR-F6翼身结合处后缘附近网格精度对分离气泡计算的影响。计算结果表明,MFlow程序的计算精度与国外软件相近,能够为阻力计算提供比较精确的结果。     

基于阻力气动特性计算的飞艇艇身外形研究

采用SIMPLE算法、S-A湍流模型求解了三维不可压雷诺平均N-S方程,数值模拟了绕飞艇的低速不可压粘性流动,计算了不同艇身外形的阻力特性.本文方法计算结果与实验结果或其他文献的计算结果符合良好,可以用于飞艇的气动特性计算.通过对6个不同外形的相同体积飞艇阻力特性计算,得到了最佳艇身外形,表明本文方法可以用于艇身外形的选型设计.     

重心位置对飞机阻力及其飞行性能的影响

为了降低配平状态下飞机的阻力问题,推导了飞机升致阻力与重心位置间的关系式。以一架典型的后平尾飞机为例。计算了向后移动飞机重心的位置,放宽对飞机固有静稳定度的要求对飞机配平升致阻力和总阻力以及对部分性能的影响。结果表明,当向后移动重心位置时,其配平升致阻力下降,并在某一重心位置时达最小值。     

较低雷诺数零升阻力估算方法研究

根据层流和湍流的摩擦阻力估算公式,建立较低雷诺数条件下当量摩擦阻力系数的计算模型来估算零升阻力系数。通过对估算方法进行验证和分析,研究较低雷诺数时零升阻力系数随雷诺数的变化趋势。结果表明,建立的估算方法合理,并且准确地反映了零升阻力系数在较大雷诺数范围内的变化规律。     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

发动机短舱溢流阻力的数值模拟

 为分析发动机溢流阻力的产生原因并为进气道-发动机匹配研究提供一定的技术支持,以典型的发动机短舱为研究对象,分析了短舱上的进气道阻力,并详细介绍了相应溢流阻力的数值计算方法。在对发动机短舱外流模拟时,通过引入两个和发动机参数相关的数值边界条件,即一个是定流量的发动机进气道入流边界,另一个是定总温总压的发动机喷流边界,避免了对发动机复杂内流的模拟。利用NACA-1-81-100发动机进气道作为标准算例,验证了数值方法的可行性。通过模拟发动机短舱在不同工作流量下的流动,认为发动机短舱的溢流阻力产生归因为溢流造成的发动机外表面的阻力变化,以及捕获流管变细导致的附加阻力增加。     

减小翼型激波阻力的鼓包流动控制技术

针对2020年使用的N+2代民用飞机的翼身融合(BWB)布局发展需要,以减小激波阻力为目标,采用计算流体力学(CFD)方法,开展弱化激波、减小激波阻力的鼓包流动控制技术研究.提出了λ形激波结构“强干扰”和等熵压缩“弱干扰”两种鼓包激波减阻流动控制原理,给出了两种鼓包基本形状设计方法和工程应用的可行性分析,指出λ形激波结...