结构/非结构混合网格数值模拟栅格翼

为了详细了解栅格翼内部的流动特性,设计满足飞行器性能要求的栅格翼,本文开展了多种不同单栅格外形的气动特性数值计算,对单个栅格的外形进行了初步选型研究;然后采用结构网格、结构与非结构混合网格求解NS方程,对带有两片栅格翼的导弹进行了数值模拟,结构网格计算结果与结构非结构混合网格计算结果一致,并与国外文献的计算结果和风洞试验结果进行了对比,三者的结果吻合较好.采用混合网格比结构网格的网格规模小许多,可以节省大量计算机内存和计算时间,是一种比较好的栅格翼数值模拟方法.     

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

混合网格方法研究栅格翼亚跨超声速气动特性

采用基于结构/非结构混合网格的CFD方法对栅格翼的亚跨声速气动特性进行了研究,计算方法经过试验数据的验证,可达到工程精度。对单独栅格翼的研究表明,跨声速壅塞时,流动通过栅格前方亚声速气流减速从栅格外侧溢流而实现流量调节;超声速壅塞时,流动通过在栅格翼上产生网状脱体激波使来流减速至亚声速来进行流量调节,且计算结果与理论估算上下临界马赫数范围一致。对某栅格翼气动布局导弹的计算结果表明,在跨声速壅塞区间内,壅塞有自我调节作用,通过溢流使栅格通道内维持相似流动,从而使栅格气动性能保持基本平稳;同时由于溢流使流量下降,栅格翼升力效率下降,导弹静稳定性减弱。     

栅格冀导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术，对栅格翼导弹黏性流场进行数值模拟，预测其气动特性。计算的马赫数为0．7～2．5。用有限体积法和LU-SGS算法求解N-S方程。计算与实验数据以及栅格翼附近流场结构进行了比较。     

栅格翼导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术,对栅格翼导弹黏性流场进行数值模拟,预测其气动特性.计算的马赫数为0.7～2.5.用有限体积法和LU-SGS算法求解N-S方程.计算与实验数据以及栅格翼附近流场结构进行了比较.     

栅格绕流数值模拟研究

栅格长宽比是影响栅格翼气动特性的重要参数,为研究它对栅格翼气动特性的影响规律,本文通过求解NS方程数值模拟研究单栅格长宽比对栅格气动特性的影响,并且对栅格边框厚度和形状对阻力的影响,在亚、跨、超声速多个马赫数下进行了数值模拟.结果显示,在所选择的长宽比范围内,长宽比对阻力的影响不大,对法向力有明显的影响,选择适当的栅格边框外形和厚度,可以大大减小栅格的阻力.     

栅格翼外形特征对减阻影响的研究

栅格翼相对平板翼有其独特的优越性,对其气动特性进行优化很有必要。采用风洞实验和数值计算的方法,分别对不同翼弦格宽比的栅格翼及不同后掠方式的栅格翼进行了研究,风洞实验结果显示,栅格翼的翼弦格宽比存在一个气动性能的最佳值,使得升阻比最大;数值计算结果证明栅格翼前缘局部后掠能有效减小波阻,是一种新的减小栅格翼阻力的方式。     

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

带栅格翼导弹超声速阶段滚转阻尼导数的数值研究

滚转阻尼导数是判断导弹动稳定性的重要气动参数。采用求解定常流场的方法,对超声速阶段平板翼翼身组合体和栅格翼翼身组合体的滚转特性进行数值模拟,并与平板翼翼身组合体的实验数据相比,符合较好。计算结果表明,栅格翼的滚转阻尼导数变化复杂,由于弹体背风面分离流和栅格翼翼面失速的影响,滚转阻尼导数随着攻角增加有两次明显转折;平板翼的滚转阻尼导数随着马赫数增加逐渐减小,而栅格翼的滚转阻尼导数随马赫数同样呈现两次转折,在Ma=3.5达到最大值。     

结冰对翼型及翼身组合体气动特性影响研究

飞机结冰是影响飞机飞行安全的重要因素之一,研究结冰对飞机气动性能的影响具有重要意义。应用多块结构化网格生成方法,通过求解雷诺平均N-S方程分别对翼型及DLR-F4翼身组合体干净外形和三种不同带冰外形的流场进行计算,分析不同冰型对翼型及翼身组合体绕流流场及气动特性的影响。     

喷流流动控制增升效应数值模拟

飞机的增升系统是影响其起降性的重要因素,为了研究飞机下表面喷流流动控制增升效应,采用数值模拟方法与结构网格求解不可压Navier-Stokes方程,对NACA0012翼型和机翼进行研究。通过改变喷气口的参数来研究其对翼型和机翼气动特性的影响。对于翼型,主要研究喷气口速度和位置对翼型气动特性的影响;对于机翼,主要研究喷气口区域沿展向分布和喷气速度对机翼气动特性和能量利用率的影响。计算结果表明:采用喷气控制可以获得较好的气动特性,增升效果也比较明显,控制效果和上述几个参数均有一定关系。本次研究为机翼喷气口设计以及机翼局部喷气控制展向分布的选取提供了依据。     

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

基于五阶空间离散精度的WCNS格式,开展了CRM翼/身/平尾组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测方面的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细3套网格,对应的y⁺分别为1、2/3和4/9。通过对CRM模型的计算和分析,研究了网格规模对气动特性、压力分布和翼根后缘局部分离区的影响。通过与DPW4统计结果和部分实验结果的对比,高阶精度数值模拟结果表明,本文的气动特性计算结果与统计平均结果吻合较好;网格密度对机翼内侧的展向站位压力分布影响较小,对机翼外侧展向压力分布影响较大;网格密度对翼根后缘局部分离区略有影响。     

旋转机翼飞机旋翼飞行验证平台设计

为了更好地研究旋转机翼飞机旋翼模式下的飞行性能,针对旋转机翼飞机的特点,开展旋转机翼飞机旋翼飞行验证平台的设计,主要包括旋转机翼的结构设计、动力系统设计.通过数值计算的方法,对不同相对厚度椭圆翼型气动特性进行分析研究.建立旋转机翼悬停气动特性的估算方法,地面试验结果和计算结果的对比表明所建立的估算方法是有效的,该方法对今后开展旋转机翼飞机研究具有一定的参考价值和意义.     

半展长襟翼梯形翼构型数值模拟技术研究

基于雷诺平均的 Navier-Stokes 方程和拼接结构网格技术,采用 MUSCL 格式和 SST 湍流模型,研究了网格密度对半展长襟翼梯形翼高升力构型的数值模拟结果的影响。相应的风洞试验是1998年在 NASA Ames 12英尺增压风洞(PWT)中完成的,试验结果包括了总体气动特性、压力分布。研究内容主要包括网格密度对收敛历程、气动力特性、压力分布和表面流线的影响,以及气动力特性随迎角的变化。研究表明,Ma=0.15,α=16.7°时,网格密度对收敛历程、典型站位压力分布和表面流态基本没有影响,气动力特性随网格密度单调变化;采用不同密度的网格,典型剖面的压力分布与试验结果吻合良好;与修正后的试验数据相比较,数值模拟得到的失速迎角前的气动力系数与试验结果吻合良好。     

悬停状态倾转旋翼/机翼干扰流场及气动力的CFD计算

基于一套高效通用的多层运动嵌套网格技术,建立了适合悬停状态倾转旋翼机旋翼/机翼气动干扰特性分析的高效混合计算流体力学（CFD）方法。在倾转旋翼/机翼贴体网格区采用可压雷诺平均Navier-Stokes（RANS）方程作为主控方程,过渡/背景网格区选用Euler方程,湍流模型选用Spalart-Allmaras模型。时间推进上采用双时间推进格式进行非定常求解,并在方法中运用了SPMD（Single Program Multiple Data）模式的并行加速技术。在此基础上,首先分别采用UH-60A直升机旋翼及XV-15倾转旋翼机旋翼作为数值算例,验证了CFD方法的有效性。然后着重对倾转旋翼/机翼的非定常干扰流场及气动力分布特征进行了数值研究,模拟得到与实际情况相符的“喷泉效应”干扰现象。计算结果表明,干扰作用使得倾转旋翼相对于孤立旋翼拉力数值减小了3%,但总的拉力系数损失达到了17%,证明悬停状态下气动干扰对飞行器气动性能有重要影响。     

螺旋桨滑流对菱形翼布局无人机气动的影响

采用基于混合网格技术及k-k_L-ω转捩模型求解雷诺平均Navior-Stokes （RANS）方程的多重参考系（MRF）方法对带桨状态的低雷诺数菱形翼布局无人机（UAV）的气动特性进行了准定常数值模拟,通过对无人机带桨状态和干净构型的气动力系数及流场结构特性进行对比分析,研究了螺旋桨在不同安装位置时其滑流对菱形翼布局无人机气动特性的影响。研究结果表明,螺旋桨滑流并不总能提高无人机的升力特性;螺旋桨安装在机头及前翼时后翼受到螺旋桨滑流形成的组合涡系的影响,其气动性能有较大的变化;螺旋桨滑流对机翼的增升作用还受到机翼掠角和螺旋桨旋转方向的影响;受布局特性影响,当安装位置远离焦点时,螺旋桨滑流对无人机的俯仰力矩特性影响较大。     

基于伴随方法的机翼多设计点气动反设计方法

针对三维机翼多点多约束气动反设计问题,基于伴随理论方法和粘性流雷诺平均N-S方程,通过粘性流数值模拟、伴随方程与梯度精确数值求解、计算网格高效算法及梯度类优化算法等有效结合,并采用考虑多设计点梯度权重系数的并行计算近似模式,开展了一种三维机翼多设计点多约束气动反设计方法研究,进行了典型算例验证。研究表明:所发展的机翼多设计点气动反设计方法具有较好的鲁棒性及优化效率。     

多面体网格在静气动弹性计算中的应用

介绍了多面体网格技术在静气动弹性问题中的应用.基于M6算例和径向基函数插值法,探究了多面体网格的CFD计算效率和变形能力,使用AGARD445.6机翼进行了静气动弹性数值计算.结果表明,多面体网格具有计算效率高、变形能力强、计算结果准确等优点,可适用于静气动弹性问题的数值计算.     

Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode

A rotor CFD solver is developed for simulating the aerodynamic interaction phenomenon among rotor, wing and fuselage of a tilt rotor aircraft in its helicopter mode. The unsteady Navier–Stokes equations are discretized in inertial frame and embedded grid system is adopted for describing the relative motion among blades and nacelle/wing/fuselage. A combination of multi-layer embedded grid and ‘‘extended hole fringe" technique is complemented in original grid system to tackle grid assembly difficulties arising from the narrow space among different aerodynamic components, and to improve the interpolation precision by decreasing the cell volume discrepancy among different grid blocks. An overall donor cell searching and automatic hole cutting technique is used for grid assembly, and the solution processes are speeded up by introduction of Open MP parallel method. Based on this solver, flow fields and aerodynamics of a tilt rotor aircraft in hover are simulated with several rotor collective angles, and the corresponding states of an isolated rotor and rotor/wing/fuselage model are also computed to obtain reference solution.Aerodynamic interference influences among the rotor and wing/fuselage/nacelle are analyzed,and some meaningful conclusions are drawn.