摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

对转螺旋桨流场气动干扰数值模拟

基于结构网格动态面搭接技术,通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程,对单独螺旋桨(SRP)滑流流场进行了数值模拟,结果与风洞试验数据吻合良好,在此基础上深入研究了对转螺旋桨（CRP）的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征,适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因,后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下,对转螺旋桨气动力会产生周期性波动,波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速,因此在相同工况下,对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍,效率提升约1.5%。      

离心式压气机模型级内非定常流动的数值试验

应用同位网格下的压力校正算法和滑移网格技术 ,在贴体曲线坐标系下 ,对离心式压气机模型级内部由于叶轮与叶片扩压器相互作用引起的三维非定常湍流流场进行了数值模拟。计算结果表明 ,在设计工况下 ,由于叶轮与叶片扩压器的相互作用引起的非定常干涉流动 ,对整级内部流动有非定常影响。速度场的非定常特性主要表现在叶轮与叶片扩压器间的径向间隙范围内和扩压器的入口处 ,在叶轮和扩压器叶道内部衰减很快。而压力脉动同时向流道的上下游传播 ,在叶轮和扩压器内部流动仍具有显著的压力非定常脉动。计算结果还模拟了扩压器入口处运动涡系与扩压器叶片的相互作用及其周期性衍灭过程     

对转桨扇三维流场特性数值研究

对转桨扇（Contra-rotating propfan,CRP）是下一代民用航空推进备选方案开式转子发动机最重要的气动部件,其气动性能对整机性能影响显著。本文对不同进距比下的对转桨扇三维流场进行数值模拟,结合压气机及螺旋桨相关理论分析了对转桨扇内部流场及其滑流区涡结构和滑流特征。结果表明,对转桨扇后排流动特征及性能参数变化幅度均超过前排。对转桨扇实际进口气流角受到轮毂附面层、诱导速度、抽吸效应的共同影响,可根据不同叶高位置轴向速度的分布规律判断三种影响分别起主导作用的位置。在对转桨扇滑流区中,桨尖涡是导致损失的主要原因,径向涡量衰减比周向和轴向涡量衰减更快。对转桨扇滑流在径向上影响至3.5倍叶高位置。气流出后排桨扇后会持续加速直至静压达到环境压力,加速区域长度约为5倍桨扇半径。     

直升机前飞旋翼桨尖涡定常质量射流控制计算

针对直升机前飞状态开展了旋翼非定常桨尖涡的模拟,采用定常桨尖空气质量射流(TAMI)控制方式对旋翼桨尖涡进行了控制。建立了一个适用于前飞旋翼桨尖涡高精度捕捉和质量射流控制模拟的数值方法,在该方法中,采用有限体积法进行空间离散,应用5阶Roe-WENO(weighted essentially non-oscillatory)格式进行流场重构及控制面对流通量计算;采用双时间方法进行时间推进,伪时间步上应用隐式LU-SGS(lower upper symmetric Gauss-Seidal)格式;引入射流边界条件对质量射流进行模拟;采用运动嵌套网格方法对前飞旋翼桨尖的挥舞运动进行模拟,并对桨叶网格和背景网格进行针对性加密。基于所建立的方法对前飞状态旋翼非定常桨尖涡及其质量射流控制进行了模拟,计算结果表明:前飞状态下旋翼桨尖涡存在较大的前后差异,桨盘前侧的桨尖涡涡核强度远小于桨盘后侧;桨盘前侧旋翼桨尖涡的涡核强度很难由定常质量射流控制来减弱,而桨盘后侧的旋翼桨尖涡则可以通过定常质量射流得到有效控制。     

二维扩压叶栅非定常涡面分析

轴流压气机内流动是典型的非定常旋涡流动,而现有压气机分析设计体系在定常假设下以流场空间结构为对象.研究压气机内真实流动的途径之一是分析流场时空结构.把转/静叶片交错排列简化为单叶片排,再简化成2-D扩压叶栅.从叶栅非定常涡面的角度来描述非定常旋涡流场时空结构.作为举例,利用合成射流控制涡面演化,改善流场时空结构.此涡面简化模型可为进一步研究叶轮机叶排间干扰复杂流场的时空结构提供分析工具.      

螺旋桨黏性流场非定常数值模拟

基于多块结构滑移网格,采用多参考系模型（MFR）,对某六叶螺旋桨黏性流场进行了非定常数值模拟.分别建立包围螺旋桨的旋转区域和螺旋桨以外的静止区域,很好地解决了螺旋桨的相对旋转问题.计算结果与实验结果吻合良好,验证了计算方法的正确性.对前进比为0.90和0.75情况下的桨后复杂流场进行分析,得到了螺旋桨尾涡的一些静态和动态特征,螺旋桨尾涡中包含正负两层,与实验观察到的现象一致.获得了螺旋桨对桨后轴向速度分布的影响规律:螺旋桨对桨叶45%和70%半径处轴向速度增加远大于90%半径处轴向速度;随前进比减小,该区域内轴向速度继续增加,90%半径处轴向速度几乎不变.      

桨扇三维流场的椭圆型微分方程组计算网格生成法

为了对桨扇周围的跨音流场有深入的理解,用计算流体力学方法在巨型计算机上进行流场数值模拟是一条重要途径。因此引起学术界的重视。桨扇叶片具有又弯又扭的三维形状,就决定了桨扇周围流场的三维性质,给桨扇跨音流场的数值模拟带来很大的困难。为了数值求解过程中改善解的稳定性和收敛性,以及为了提高数值解的精度,就必须妥善解决三维网格系统生成问题。     

短舱对螺旋桨滑流影响的IDDES数值模拟

基于非结构重叠网格技术,对短舱与螺旋桨滑流间的相互作用进行了非定常数值模拟研究。为了更好地捕捉螺旋桨尾涡的细节信息,计算采用基于Spalart-Allmaras模型的改进延迟脱体涡模拟(IDDES)方法,并在非定常计算过程中运用网格自适应技术以提高流场特征的空间分辨率。研究结果表明:IDDES方法获得的拉力系数计算值与实验值吻合良好,短舱的存在会增大螺旋桨的拉力系数;短舱对螺旋桨桨毂涡的结构影响较大,但对桨尖涡的螺旋结构影响较小;对单独螺旋桨算例来说,桨尖涡与桨毂涡的失稳发展过程都具有周期性,且在有/无短舱情况下桨尖涡的失稳位置相同,失稳后桨尖涡之间配对融合过程一致,从而说明桨毂涡对桨尖涡的失稳没有影响。     

基于DDES方法的叶栅分离旋涡的非定常流动数值研究

为了研究叶栅分离旋涡的非定常流动特性,选取亚声速叶栅为研究对象,进行了基于延迟分离涡模拟方法 (DDES)的非定常数值模拟。针对进口Ma0.3的情况,对比分析了进口攻角、叶栅稠度、叶片弯角、叶片中弧线弯度分布和叶片厚度等因素对于叶栅分离旋涡运动的影响。结果表明,当攻角增大时,尾迹附近流场会逐渐由定常模式转化为非定常对涡脱落模式,同时对涡脱落频率不断减小;随着攻角的进一步增大,叶背分离强度不断增加,叶背分离涡开始单独脱落,并主导流场旋涡运动。强烈的叶背分离涡脱落引起了叶片气动力的剧烈脉动,其幅值是小攻角状态下由尾迹对涡脱落引起叶片气动力脉动的30倍以上。叶栅几何参数的改变对于流场旋涡运动同样有很大影响,流场同样呈现出"定常尾迹"、"对涡脱落"、"叶背分离涡脱落"这三种主要的旋涡运动形式之一。叶背分离点的相对位置则是影响旋涡非定常运动形式和旋涡脱落频率的主要内在因素。     

非结构CFD软件MPI+OpenMP混合并行及超大规模非定常并行计算的应用

常规工程应用中，非定常数值模拟（如多体分离）的计算量十分巨大，如果为了达到更高的计算精度，加密网格或者采用高精度方法将会使得计算量进一步增大，导致非定常数值模拟在CFD工程应用中成为十分耗时和昂贵的工作，因此，提高非定常数值模拟的可扩展性和计算效率十分必要。为充分发挥既有分布内存又有共享内存的多核处理器的性能和效率优势，对作者团队开发的非结构网格二阶精度有限体积CFD软件（HyperFLOW）进行了混合并行改造，在计算节点间采用MPI消息传递机制，在节点内采用OpenMP共享内存的MPI+OpenMP混合并行策略。首先分别实现了两种粒度（粗粒度和细粒度）的混合并行，并基于国产in-house集群采用CRM标模（约4 000万网格单元）定常湍流算例对两种混合并行模式进行了测试和比较。结果表明，粗粒度在进程数和分区数较少的小规模并行时具有效率优势，16线程时效率较高；而细粒度混合并行在大规模并行计算时具有优势，8线程时效率较高。其次，验证了混合并行在非定常计算情况下的可扩展性，采用机翼外挂物投放标模算例，分别生成3.6亿和28.8亿非结构重叠网格，采用对等的（P2P）网格读入模式和优化的重叠网格隐式装配策略，网格读入和重叠网格装配耗时仅需数十秒；采用3.6亿网格，完成了非定常状态效率测试及非定常分离过程的湍流流场计算，在in-house集群上12 288核并行效率达到90%（以768核为基准），在天河2号上12 288核并行效率达到70%（以384核为基准），数值模拟结果与试验结果符合良好。最后，在in-house集群上采用28.8亿非结构重叠网格进行了4.9万核的并行效率测试，结果显示，4.9万核并行效率达到55.3%（以4 096核为基准）。     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

倾转旋翼典型飞行状态气动特性的CFD分析

针对倾转旋翼特有的悬停、巡航、过渡3种飞行模式,基于CFD方法分别对倾转旋翼的上述3种典型飞行状态进行了数值模拟.首先,基于运动嵌套网格思想建立了适合倾转旋翼运动模拟的多层运动嵌套网格系统.然后,基于RANS(Reynolds averaged Navier-Stokes)方程建立了适用于非定常流场分析的CFD模拟方法.为了提高整个流场的计算效率,时间推进采用高效的隐式LU-SGS(lower upper symmetric Gauss Seidel)方法,并运用SPMD(single program multiple data)并行计算模式.最后,采用所建立的网格生成方法和流场求解器,进行了不同飞行状态下的倾转旋翼气动特性数值分析,验证了建立CFD方法的有效性并得出了多种飞行状态下的气动特性.结果表明:悬停状态下倾转旋翼相对于直升机旋翼有更加明显的下洗流影响区域.巡航状态两旋翼中心间距为2.5R时,倾转旋翼之间存在明显的气动干扰.过渡状态下,倾转旋翼提供的拉力在倾转角为45°前下降较快,且在倾转角较小时,倾转旋翼桨叶表面拉力分布仍能表现出直升机旋翼的特性.      

渐变反推力平面叶栅端壁流动仿真与试验

针对某型反推力装置中的渐变反推力平面叶栅,开展了性能预测和流动试验的研究与验证。结合数值仿真和性能试验,运用表面油流等测量手段,探索了反推力叶栅流线分布和流动的基本结构,揭示由不同弯角组成的反推力叶栅旋涡结构和流动特征。试验与仿真结果表明:反推力叶栅端壁呈现出3维非定常有旋的流动特征,其发展和演变的过程与涡轮叶栅中的流动过程相似,依据仿真结果和流道收敛的特征,确立了反推力平面叶栅的基本测量方法,并通过了试验结果的验证,对今后的反推力叶栅数值仿真和试验工作有很大的借鉴意义。     

叶栅安装角异常的非定常流场数值模拟

采用商用软件NUMECA对叶栅安装角异常进行二维流场定常和非定常数值模拟,并进行对比分析,得到了因叶栅安装角异常造成通道堵塞对该叶片以及相邻叶栅通道流场结构和非定常气动力的影响.结果表明:安装角异常角度大于0°时,主要对沿吸力面方向叶栅通道流场结构影响较大,随着安装角异常角度的增大,气流不断恶化,出现大面积的附面层分离,尾涡脱落主频率随着安装角异常角度的增大而减小,次频率随着安装角异常角度的增大而增大;并导致叶片非定常气动力相对脉动量迅速增大,这可能是导致叶片疲劳破坏的原因之一.     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

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

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