DPWⅢ机翼和翼身组合体构型数值模拟

采用亚跨超CFD软件平台(TRIP)数值模拟了DPW Ⅲ提供的DPW_W1/W2两种机翼构型和DPW-F6/F6_FX2B两种翼身组合体构型,主要目的是通过两种机翼构型和两种翼身组合体构型的数值模拟,研究网格密度对运输机构型气动特性计算结果的影响。数值模拟采用的多块对接网格来自AIAA CFD Drag Prediction Workshop Ⅲ(DPW Ⅲ),采用National Transonic Facility(NTF)的试验结果和CFL3D的计算结果作对比。详细研究了网格密度对两种机翼构型和翼身组合体的总体气动特性和压力分布的影响。采用SST两方程模型计算两种构型均得到了网格收敛结果,网格密度主要影响压差阻力而对摩擦阻力影响较小,计算结果较好地预测了机翼和翼身组合体外形优化前后总体气动特性的变化量。     

TRIP2.0软件的确认：DPWⅡ复杂组合体的数值模拟

 采用“亚跨超声速计算流体力学软件平台”（TRIP2.0）数值模拟了阻力预测小组（AIAA CFD Drag Prediction Workshop Ⅱ,DPWⅡ）翼/身/架/舱复杂组合体运输机构型,数值模拟采用的多块对接网格、测压和测力的试验结果均来自DPWⅡ,对比计算采用了CFL3D的结果。重点针对DLR-F6翼/身/架/舱复杂组合体构型,详细研究了网格密度和湍流模型对总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致。采用SA一方程和SST两方程模型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布有一定的影响。     

基于广义Richardson外插方法的颤振模拟耦合时间精度研究

为更准确地评估颤振问题时域模拟在实际计算中表现出的时间精度,采用不同精度的耦合方法对Isogai Wing和AGARD 445.6 Wing的颤振问题进行时域模拟。参照网格收敛性分析方法,提出了基于广义Richardson外插方法的颤振问题时域模拟耦合时间精度分析方法,分析时域计算结果的时间步长收敛性,计算数值模拟结果时间精度并获得时间步长无关解。分析表明,对于Isogai Wing和AGARD 445.6 Wing颤振问题,计算结果具有良好的时间收敛性,采用广义Richardson外插方法对各耦合方法分析所得精度与理论分析结果基本一致。在合适的时间步长区间内,可忽略具体时间步长选取对广义Richardson外插方法分析结果的影响,验证了提出的气动弹性模拟耦合时间精度分析方法。     

采用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方法分析了雷诺数效应对某运输机全机构型阻力的影响。基于对某运输机的翼/身/架/舱组合体绕流进行的计算分析,归纳总结出了一种适应于复杂构型的低雷诺数到高雷诺数的阻力修正方法。通过验证可知,该修正方法修正效果良好,为工程中的阻力预测提供了一种简洁有效的方法。     

包含支撑装置和机翼变形的CRM-WB构型气动特性数值模拟

CRM(Common Research Model)是第四届~第六届AIAA阻力预测研讨会(DPW)选择的参考外形。第四届和第五届DPW会议气动特性数值模拟结果的统计分析表明试验与计算之间存在明显差异。采用CFD方法和流固耦合(FSC)方法数值模拟了带支撑装置的CRM翼身组合体(CRM-WBS)构型的气动特性,以评估模型支撑装置和静气动弹性变形对CRM翼身组合体(CRM-WB)构型气动特性数值模拟结果的影响。通过与不包含支撑装置的CRM-WB构型的数值模拟结果和NASA Langley NTF风洞试验结果的对比分析,CRM-WBS构型的CFD计算结果表明,支撑装置导致机翼上翼面激波位置前移,升力系数、阻力系数、低头力矩系数下降。CRM-WBS构型的流固耦合计算结果表明,静气动弹性变形主要影响机翼上表面激波位置并显著降低外侧机翼激波位置前的负压值,进一步导致升力系数、阻力系数、低头力矩系数下降。包含支撑装置和静气动弹性变形的CRM-WB构型气动特性数值模拟结果更加接近试验结果。     

NASA共同研究模型翼身组合体数值模拟研究

采用结构网格求解器WiseManPlus软件,对DPW5提供的NASA共同研究模型CRM翼身组合体构型进行数值模拟研究.采用雷诺平均N-S方程,选择SA一方程湍流模型和SST两方程湍流模型,开展了网格收敛性分析及抖振分析计算研究,评估了该软件对此类民机构型的阻力预测能力.按照DPW5的要求,采用会议统一提供的基准结构网格进行计算分析,研究结果表明,计算的巡航设计点的阻力及大攻角下的气动力系数与试验值吻合较好,达到会议统计分析平均水平.计算结果表明,所采用的结构网格求解器WiseManPlus软件计算精度较高,适用于大型民机高速巡航及抖振现象的数值模拟研究.     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

三维侧压式进气道的减阻研究

为实现三维侧压式进气道的减阻,采用了数值模拟和来流马赫数Ma₀=5.3风洞试验相结合的方法研究了两种三维侧压式进气道,其中之一为基准进气道,其二是在基准进气道的基础上采用某些减阻措施的低阻进气道.结果表明:所采用的减阻措施是有效的,在Ma₀=5.3的条件下,进气道的压差阻力可以减少4.7%左右,摩擦阻力可以减少约4.9%,附加阻力可以转化为附加推力,总阻力可以减少4.7%左右.与此同时,进气道的总体性能也得到明显改善.      

内外环对高速球轴承中球体拖曳阻力系数的影响

在液体润滑的高速球轴承中,球体公转运动会受到润滑流体的拖曳阻力,对球轴承的动力及产热等性能造成显著影响。针对球轴承中球形滚动体公转的流体拖曳阻力问题,建立平板夹层单球体绕流模型,通过数值模拟和实验研究的方法,研究内外环对球体绕流的流动特征及压力分布的影响,分析单球体绕流的拖曳阻力系数变化规律,并与开放空间球绕流模型进行对比。研究表明,建立的数值计算方法和实验结果吻合良好,误差不超过10%。在相同雷诺数下,由于内外壁面的限流作用,球体与壁面接触区域压差阻力明显增大,导致夹层空间球绕流中球体拖曳阻力系数明显高于开放空间球绕流,在103Re105内约为开放空间球体阻力系数的两倍。     

二元高超声速进气道内部阻力特性分析与研究

以二元高超声速进气道为研究对象,以数值模拟为手段,详细分析了该进气道内壁面上的压差阻力、摩擦阻力、总阻力随飞行马赫数的变化。接着用流线跟踪得到该二元高超声速进气道的分流面,研究了作用在该分流面上的附加阻力随飞行马赫数的变化。为了更准确地评价该进气道的阻力特性,采用了进气道的推力效率和推阻比两个概念,并且探讨了它们随飞行马赫数的变化。最后,分析讨论了作用在该进气道捕获流管上的实际总阻力随飞行马赫数的变化。     

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

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

非线性壁板颤振分析

利用一种改进的计算流体力学与计算结构动力学(CFD/CSD)耦合方法研究了由气动和结构几何非线性引起的壁板颤振问题。在非定常气动力计算中,考虑了通量分裂格式、隐式时间推进方法和几何守恒律;二维和三维壁板的结构几何非线性建模则采用了有限元的协同旋转理论,并利用一种近似能量守恒算法求解结构的非线性响应。流场和结构求解器采用二阶松耦合方法联立求解,并将其应用于壁板在超声速、跨声速和亚声速的颤振计算中。当考虑结构几何非线性和气动非线性时,出现了典型的极限环振荡现象,并对颤振边界和极限环振荡幅度进行了对比分析。     

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

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

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

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

Preliminary study on heat flux measurement data of TT-0 flight test

With the explosive development of aerospace science, the design of the new generation airliner at higher speeds is attracting more attentions. To achieve this goal, it is necessary to achieve accurate prediction of the aerodynamic heating / force loads and successful reduction of drag and heat flux. As a remedy for the existing studies which are based upon the CFD and wind tunnel tests, this study presents a flight test for the drag and heat reduction spike technology. The principal goals of this flight test were to provide reference for verifying the accuracy of the prediction technology on ground and promote the development of the drag and heat reduction technology. By adopting the OS-X rocket, the TT-0 test vehicle designed by Shenyang Aircraft Design & Research Institute reached a maximum Mach number of 5.8 and a maximum altitude of 38 km. Hypersonic and supersonic pressure data by pressure scanning valves and heat fluxes by gauges at different locations were obtained successfully. Also, heat fluxes obtained by in-house CFD code are illustrated in comparison with the flight data. The results indicate that the numerical errors are large in most cases. More technologies, such as more CFD codes and more numerical procedures, should be adopted to conduct studies on this issue in the future.     

运输机大上翘机身后体参数影响研究

采用CFD方法研究了外形差异较大的三种大上翘后体机身的阻力与流动特性,分析了不同后体形状及几何参数的影响。压差阻力虽仅占机身总阻力的15~20%,但却决定总阻力的大小。在上翘角和长细比接近的情况下,扁平度是影响流动分离和压差阻力的主要因素。给出了后体设计的参数选择及应采取措施。     

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

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

二维气动问题中Kriging代理模型精度影响因素

以二维跨声速临界翼型的阻力特性为对象,探讨样本点数目、Kriging代理模型参数及其类型等对模型精度的影响.阻力系数采用计算流体力学（CFD）方法得到.模型精度的验证采用交叉验证方法,采用平均误差、最大误差和标准交叉验证残差来衡量Kriging代理模型的精度.研究结果表明:①Kriging代理模型预测气动阻力效果较好.②模型精度随样本点的增多而提高,剔除与样本点响应趋势不相符的“奇异点”后,模型精度显著提高,平均误差减小5%~38%,最大误差减小13%~77%.③核函数类型对模型精度的影响最大,相关参数次之,回归模型的影响最小.采用高斯相关函数、2阶多项式回归模型,以及合适的相关参数值时,Kriging代理模型的精度最高.