螺旋桨非定常滑流的高效数值模拟研究

基于自主研发的大型"亚跨超CFD软件平台"(TRIP3.0),开发了针对螺旋桨飞机的旋转机械非定常求解模块。控制方程采用曲线坐标系下的非定常Euler方程,综合运用动态拼接结构网格技术和双时间步推进方法,重点模拟了双发螺旋桨飞机在螺旋桨转动时的非定常滑流效应。为了保证计算模拟的高效性,引入多重网格技术和大规模并行计算技术。通过模拟某单独螺旋桨验证算例,验证了该方法的可行性。模拟了不同迎角和螺旋桨转速下全机的气动特性,分析了滑流对流场的影响情况,结果表明:螺旋桨滑流对全机流场和气动特性影响明显,并使得升力系数增大。     

风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0)，采用刚性运动网格技术和动态拼接网格技术，开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究，影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程，采用有限体积法离散控制方程，离散方程的时间方向采用“双时间步”方法求解，空间方向无粘项离散采用 MUSCL-Roe 格式，湍流模型采用 SA 和 SST 湍流模型，并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构，进而影响吸力面的压力分布，而对压力面的流动结构和压力分布基本没有影响；采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

基于动态混合网格的多体分离数值模拟方法

在以往发展的动态混合网格技术的基础上,耦合六自由度运动方程计算和非定常流场计算,建立了一种多体分离问题的非定常数值模拟方法。首先采用刚体运动的六自由度运动方程,建立刚体运动轨迹计算模块,然后以该模块为纽带将以往建立的非定常计算方法与改进的基于Delaunay背景网格插值方法和局部网格重构方法相结合的动态混合网格生成方法耦合起来,形成了精度较高、效率较快的多体分离数值模拟方法。利用该方法对典型多体分离问题进行了数值模拟,计算结果与风洞实验结果吻合较好。     

微型扑翼飞行器非定常运动对平尾的影响

以西北工业大学自行研制的微型扑翼飞行器ASN211为研究对象,利用其简化的二维扑翼及平尾串列翼模型进行了非定常数值模拟,分析了扑翼俯仰运动及沉浮运动对平尾气动性能的影响。在数值模拟模块中,模型的俯仰运动及沉浮运动由动网格技术实现。通过计算流体力学(CFD)软件Fluent对此非定常流场进行数值计算,重点研究了扑翼非定常运动尾流对平尾气动效率的影响。定常状态与非定常时均条件下平尾升力曲线的对比分析表明,扑翼的非定常运动能够增大平尾的失速迎角及最大升力系数,因而使平尾的失速特性得到改善。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

跨声速机翼操纵面非定常气动力Euler方程计算

采用代数法与椭圆型方程优化相结合的结构生成方法为带操纵面F5机翼生成面对接块结构网格;弹性变形技术生成动态网格;有限体积法、双时间推进法求解非定常Euler方程.给出一种适于面对接块结构网格计算的守恒量插值方法.数值计算了F5机翼操纵面作简谐运动时的非定常气动力,计算值与试验值吻合较好.     

基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

CFD在螺旋桨飞机滑流影响研究中的应用

通过求解绝对坐标系下的非定常雷诺平均Navier-Stokes方程，计算了非定常滑流对某飞机巡航构型低速状态的气动特性影响。为描述螺旋桨桨叶的相对运动，采用了动态重叠网格技术，并在并行环境下采用全隐式双时间步方法和多重网格技术来保证时间精度和提高计算效率。将计算结果与试验结果相比较，在线性段，计算得到的升力系数与试验值总体偏差不超过2%，阻力系数相对试验值而言整体偏大，基本控制在8%左右，俯仰力矩系数的计算结果与试验值趋势一致。总结了滑流对螺旋桨飞机宏观气动特性的影响规律，由于滑流在空间发展过程中不断与周围空气相混合，使得滑流边界逐渐"模糊"，因此滑流流管在空间发展过程中的巨大畸变对滑流影响的分析工作带来了难度。针对这种现象，本文将有/无滑流的流场进行对比，通过当地动压增量来定义滑流的加速效应边界，以及通过当地气流角增量来定义滑流的洗流效应边界。该方法能较好地捕捉和解释由于滑流对飞机部件干扰而使得飞机方向安定性呈现的非线性现象，初步揭示了螺旋桨滑流复杂尾迹流动的特点，在螺旋桨飞机非定常滑流影响机理研究方面具有一定的参考价值。     

基于动网格的非定常粘性空气动力学计算

提出了一种将阵面推进法和半扭转弹簧法相结合的动网格生成方法。运用阵面推进法决定附面层动网格的运动规律,半扭转弹簧法生成附面层外的非结构动网格。使得算法能够同时适用于非定常粘性与非粘性空气动力学计算。文中对于非定常可压缩N-S方程,在非结构动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代,计算了二维NACA0012翼型俯仰简谐运动的低雷诺数非定常空气动力,计算结果与参考文献提供的结果吻合良好。此外,本文还对三维1/4椭圆薄板的低雷诺数非定常空气动力进行了计算。     

基于预估校正和嵌套网格的虚拟飞行数值模拟

 针对导弹虚拟飞行数值模拟问题,发展了空气动力学/飞行力学数值计算方法和软件。控制方程为非定常雷诺时均Navier-Stoker(RANS)方程和刚体六自由度运动方程;流场求解器为有限体积法结构网格求解器,时间推进采用双时间步法,湍流模型为Spalart-Allmaras一方程模型;采用Adams预估校正法实现飞行力学方程与流场控制方程的耦合计算;使用嵌套网格方法模拟多体运动。首先模拟了美国国家航空航天局(NASA)窄条翼导弹模型纵向虚拟飞行,研究耦合方式和时间步长的影响。仿真结果表明,双时间步三阶Adams耦合方法,同等精度下可以显著增大时间步长,缩短仿真时间。最后,采用该方法模拟了导弹自由摇滚特性和纵向虚拟飞行,模拟结果与试验值吻合较好。     

绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

Improvements on computations of high speed propeller unsteady aerodynamics

This paper presents the application of the CANARI flow solver to the computation of unsteady effects in the aerodynamic interaction of a high speed propeller with the aircraft. The method is first validated on the APIAN isolated propeller test case by comparison with experiment at M=0.7. The method is then applied to the time accurate 3D Euler computation of a generic transport aircraft at M=0.68. Analysis of the results shows significant unsteady effects both on the propeller forces and on the wing aerodynamic flows, by comparison with steady computations.     

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

后台阶三维缝隙气膜冷却的非稳态数值模拟

为了进一步深入了解涡轮叶片尾缘冷却结构的气体流动情况及冷却特性,在原有稳态计算模型的基础上建立了非稳态的计算模型,研究了不同吹风比下(0.5,2.0)的出口壁面冷却效率的分布情况。计算结果表明:(1)非稳态效应使得出口下游的湍流度增大,非稳态时均冷却效率的计算结果比稳态的要低一些。(2)吹风比为2.0时,二次流对出口附近流动起决定作用,并且冷气的横向掺混充分,主流二次流的上下掺混缓慢;吹风比为0.5时,主流与二次流的上下掺混剧烈,非稳态的计算结果在出口肋后附近的冷却效率比稳态结果有所提高。(3)非稳态的计算结果比稳态的计算结果更接近实验结果。     

2级涡轮内部流动定常与非定常计算差异研究

为获取进而认识涡轮内部流动状态,以某2级约化形式的动力涡轮为研究对象,分别对其进行定常和非定常数值计算和分析.研究表明:定常计算与非定常计算对涡轮内部流动的模拟结果,如叶片表面的压力分布、叶排进出口的气流角、叶片通道中的二次流流向涡、展向涡、叶片通道中的损失等,均存在差异;流动的非定常性越强,定常与非定常计算结果的差异越大,且该差异大小对于静叶与动叶呈相反的展向分布规律.     

考虑压气机动静干涉效应的非设计工况性能数值模拟

1引言非定常流动是叶轮机械内部固有的本质属性,而现有的设计体系中不论是在设计或分析时基本上采用定常流动的假设,从而忽略了非定常流动对叶轮机械在气动、结构强度等方面的影响。目前,国内外对非定常流动的研究比较多,但其中采用非定常流动数值模拟方法预估压气机非设计工况     

叶轮机转子内流动的三维跨声速Navier-Stokes方程组数值计算及试验验证

应用有限体积法及Baldwin-Lomax代数湍流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛.采用隐式算法对对流通量进行求解,而耗散通量采用显式算法.采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算耗散通量.关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解.应用该算法模拟了转子Rotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等.在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验有较大的差别.