低速压气机静叶非定常Clocking效应数值模拟

采用三维粘性非定常数值模拟方法研究了重复级低速压气机内部非定常流动,其中包括静叶Clocking(时序)效应以及尾迹-边界层、尾迹-尾迹和尾迹-泄漏流的相互作用.对静叶时序效应的结果分析显示,当两级静叶相对周向位置改变时最大效率变化只有0.1%,表明低速压气机中静叶时序效应不明显.另外,非定常结果流动图画显示,在非定常环境下上游周期性尾迹的通过对下游叶片边界层、尾迹以及叶尖泄漏流的发展都有一定影响.     

旋转弹气动力建模与飞行轨迹仿真

运用一种自回归滑动平均(ARMA)的时域气动力建模方法,以计算流体力学与刚体动力学(CFD/RBD)耦合仿真的输出结果为样本,对旋转弹的非定常气动力进行建模。利用建立的气动力模型与刚体动力学方程求解模块耦合,实现了旋转弹轨迹的快速仿真,并讨论了不同的建模方式对仿真精度的影响。算例结果表明:采用气动力模型与刚体动力学方程耦合仿真技术可以在不同初始发射条件下进行旋转弹飞行姿态与运动轨迹预测,且与CFD/RBD仿真结果吻合较好,证明ARMA气动力建模方法可以在保证旋转弹轨迹预测精度的同时大幅缩短仿真时间,节省计算资源。     

锯齿尾缘叶片气动特性和绕流流场的数值研究

以基于NACA 0018翼型的锯齿尾缘仿生叶片为研究对象,采用大涡模拟的方法研究锯齿相对齿宽与相对齿高对锯齿尾缘叶片的气动特性和非定常绕流流场的影响规律和机制.研究表明,尾缘锯齿参数对叶片气动性能的影响是复杂的非线性过程,在一定来流攻角范围内能提高升阻比,但失速提前.如在9.4°~14.8°来流攻角范围内,不同相对齿宽系列叶片的升阻比高于原始叶片,升阻比与锯齿相对齿宽之间没有线性关系.研究还表明,锯齿尾缘能延迟边界层分离,加速尾迹的流动掺混和能量扩散,改变非定常涡结构和涡脱落频率.相对齿高的变化对非定常流动特性的影响更为显著.尾缘锯齿诱导的二次湍流射流和吸力面侧反向涡对改变了原始叶片的绕翼环量,进而影响锯齿尾缘叶片的气动特性和绕流流场特性.      

壁温对隔离段激波串振荡的影响

为了探索存在换热条件下超燃冲压发动机隔离段内部的复杂流动机理,采用2阶精度的非定常数值模拟程序研究了壁温对隔离段内部激波串振荡特性的影响。结果表明:压力振荡曲线与激波串前缘位置之间存在紧密关系。此外,随着壁温升高,隔离段壁面压力振荡的频率减小、振幅增加。研究还发现,隔离段压力振荡标准差最大值位于激波串前缘激波及其与上下壁面边界层相互作的位置,而湍动能最大值则分布在分离区后侧剪切层区域内。      

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

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

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

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

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法,建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题,旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法,通过添加冗余约束,避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块,对基于运动嵌套网格 的空间流场进行求解,湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块,基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解,验证了网格变形模块和CFD模块的有效性,然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上,计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷,并与飞行测试数据进行了对比。计算结果表明,文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析,提高了旋翼气弹载荷的预测精度。     

空射导弹发射初始弹道数值仿真

导弹从载机上发射时,载机对导弹的流场存在干扰,从而影响导弹在发射初始阶段的姿态和运动轨迹。为了确定导弹的初始弹道参数,本文利用CFD商业软件Fluent6.1计算导弹发射时的非定常流场,并通过该软件的接口嵌入导弹的六自由度弹道方程,同时计算了导弹在载机流场干扰下的非定常气动力和初始弹道。计算结果清晰地表现出导弹与载机的分离过程和载机流场的干扰形式,并可估算出载机的流场干扰范围。     

地面效应对尾涡消散的影响研究

采用大涡模拟方法可以分析和计算近地阶段尾涡消散和运动规律,但模型繁琐、计算复杂。本研究基于镜像涡方法建立了形式相对简单、计算迅速的近地尾涡运动模型和消散模型,计算结果与激光雷达测量数据和大涡模拟数据的偏差都不超过5％,满足动态尾涡间隔研究的需要。     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。