旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

基于局部无网格的混合算法研究

研究了一种求解任意外形绕流的局部无网格混合算法.该算法计算区域整体采用了直角网格,只在物面附近嵌入局部无网格处理.整体直角网格的规则性与正交性,不仅使得网格生成快,而且计算格式实施简单,避免了与非正交性相关的计算项,从而引入误差相对较小.局部无网格处理只要求物面附近布点离散,具有灵活性,使得发展的算法适合处理任意外形.先以理想不可压流动问题为例,给出了其控制方程Laplace方程的具体求解实施过程,并给出了直槽道内柱体绕流的计算结果;接着把方法推广用于求解可压缩流动的Euler方程,对喷管等典型流动问题进行了数值模拟,并与传统有限体积法结果进行了比较,在捕捉的波系结构、表面压力系数分布等方面具有良好的一致性,表明发展的新方法是可行的.     

二维间断有限元混合网格混合算法研究

针对二维流场进行结构／非结构混合网格划分,在靠近物面处生成结构化网格,在结构化网格以外和远场之间以非结构化网格填充。使用间断有限元方法（DGM）离散Euler方程。在混合网格混合算法中,结构化网格区域使用结构DGM,非结构化网格区域使用非结构DGM,保证结构／非结构对接面处的通量守恒,并加入moment限制器消除激波前后伪震荡。对NACA0012和RAE2822翼型跨音速无粘流动的数值模拟结果说明了混合算法的有效性以及在捕捉间断方面的能力。     

基于“各向异性”四面体网格聚合的复杂外形混合网格生成方法

复杂外形粘性流动的计算网格生成一直是阻碍CFD实际工程应用的瓶颈问题之一,虽然三棱柱/四面体混合网格的出现极大地缓解了这一问题,但是在飞行器几何外形非常复杂的情况下,在粘性边界层内生成高质量的贴体三棱柱网格仍然是非常困难的事情。本文从"各向异性"四面体网格的生成算法中得到启示,发展了一种基于"各向异性"四面体网格聚合的混合网格生成方法。通过对DLR-F6翼身组合体和F16战斗机外形的混合网格生成及数值模拟表明,该方法生成的混合网格具有网格质量好、自动化程度高等特点,对提高复杂外形飞行器关键气动参数的模拟精度具有重要意义。     

预处理方法在混合笛卡尔网格中的应用研究

使用基于混合笛卡尔网格的预处理方法,对低马赫数下的定常/非定常流动问题进行了数值模拟研究。网格生成中,在物面附近生成贴体结构网格,其余计算区域使用笛卡尔网格进行填充,两套网格之间通过查找"贡献单元"实现流场数据间的传递。同时,使用基于密度的预处理方法,发展了一套可求解从低速(极低马赫数)到常规马赫数的N-S方程求解器,时间离散使用隐式双时间步LU-SGS格式,空间离散使用具有二阶精度的格心格式有限体积方法。非定常流动问题的数值模拟过程中,使用逆距离插值来进行新现网格单元上流场参数确定。通过对NACA0012翼型在低马赫数下的定常绕流和动态失速问题的数值模拟研究表明:使用预处理方法能够加快低马赫数下数值计算的收敛速度和提高计算的准确性,基于混合笛卡尔网格的N-S方程求解器能够模拟包含运动边界的低速不可压非定常流动问题。耦合预处理技术的混合笛卡尔网格方法给包含低速和常规马赫数的复杂运动边界问题的求解提供了新的思路。     

基于滑动网格的螺旋桨复杂绕流数值模拟研究

利用SST湍流模型,求解雷诺平均Navier-Stokes方程数值模拟螺旋桨粘性绕流。为了实现桨叶的相对运动,采用了滑动网格技术,即将计算区域分为静止区域和随桨叶一起运动的旋转区域。以课题组设计的1.2 m螺旋桨为例,分别对不同风速下的桨叶粘性流动进行了非定常数值模拟,计算得到的螺旋桨气动特性与实验值吻合良好,证明了方法的正确性。     

基于混合CFD方法的旋翼桨-涡干扰噪声计算

建立了一个基于Euler方程和自由尾迹分析的新的、高效的混合计算方法,用于旋翼桨-涡干扰(BVI)气动和噪声特性的研究.在该方法中,仅针对单片桨叶进行计算,通过求解Euler方程来捕捉流场中旋翼近尾迹的影响以及桨叶表面的非定常载荷,流场空间中旋翼远尾迹的位置由自由尾迹模型给出;提出了广义网格速度方法以将旋翼远尾迹的影响...     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

非定常旋涡-翼型相互作用问题数值模拟研究

针对直升机旋涡-桨叶相互作用问题,以其二维简化模型旋涡-翼型相互作用问题为研究对象,发展了基于雷诺平均Navier-Stokes方程和自适应混合笛卡尔网格的流场数值模拟方法,通过给定合适的自适应判据,动态网格自适应技术能够准确捕捉旋涡位置和保持旋涡强度。通过两类旋涡-翼型相互作用问题的数值模拟研究,揭示了旋涡与翼型相互作用过程中产生的复杂流动现象和非定常气动力变化。数值模拟结果与相关文献中的试验数据和数值计算结果吻合,且对三维桨叶的气动载荷研究和振动分析具有指导意义。     

DDES方法在复杂旋翼流场计算中的应用

应用延迟脱体涡模拟（DDES）方法,进行了复杂分离流动旋翼流场的数值模拟研究。为控制计算规模,更好地模拟旋翼桨叶近壁面附近区域流动分离及远场尾迹发展,采用了贴体与背景自适应直角网格嵌套的网格系统;同时为提高计算速度,采用了基于分布式存储的数据并行模式,实现了贡献单元搜索、洞切割和流场计算的并行化。计算工作包含一个共轴双旋翼悬停和单旋翼下降状态,目的是分别讨论DDES方法在旋翼表面附近空间区域是否存在大范围分离时计算结果与雷诺平均Navier-Stokes结果之间的差异。计算表明,对于下降状态,DDES方法预测出了更强的尾迹区分离流动;对于桨叶表面基本为附着流动的双旋翼中等桨距角情况,两种方法预测的流场结构相似,仅在下旋翼桨叶气动力和旋翼桨根下方空间区域的诱导速度计算结果上存在微细差异。     

空心气冷低压涡轮动叶气热耦合数值模拟

以空心气冷低压涡轮动叶为研究对象,采用高质量的流体域和固体域网格控制技术,带转捩模型的双方程SST湍流模型,开展了基于CFD方法的叶片气热耦合问题研究。获得了不同冷气流量比(分别为1.0%,1.38%,1.8%和2.2%)、温比(分别为2.1,2.25,2.3,2.4和2.5)和压比(分别为1.4,1.6和1.8)对叶片换热特性的影响规律,设计状态中截面按弧长平均的叶片壁面金属温度计算值较试验值偏小0.3%,气热耦合计算的叶片壁面温度分布与试验结果吻合良好,验证了气热强耦合计算方法的精度,为涡轮叶片温度场分析提供了一种有效的方法。     

周向弯曲叶片叶顶泄漏流随流量变化特性数值研究

对带有周向前弯和周向后弯叶片的低压轴流风机,叶顶泄漏流动随流量变化特性进行了数值研究和实验校核,探讨了随流量减小,周向弯曲对叶顶泄漏涡的起源和发展规律的影响.性能计算结果与试验测量较为吻合,计算结果表明:随着流量减小,泄漏涡起源向前缘方向迁移,周向弯曲使泄漏流与主流掺混加剧;周向前弯叶片泄漏涡的起源距叶片前缘最远,涡强度最弱,移动最缓慢;周向前弯稳定叶顶流动,降低端壁损失,扩大稳定工作范围.      

直升机的气动弹性问题

直升机的气动弹性问题与固定翼飞机不同,不仅要考虑单片桨叶,更要将旋翼视为一个整体,考虑其动态入流、尾迹影响以及旋翼与机身之间的相互耦合等。就单片桨叶而言,在结构动力学上,需要考虑离心力场、几何非线性以及桨叶的非线性挥舞-摆振-扭转耦合;在气动力上,需要考虑动态入流以及桨尖处可能的失速效应,本质上属于非线性气动弹性力学范畴。由于旋翼气动力通常是以周期形式通过旋翼轴传给机身,并引起机身振动,而机身运动又通过改变桨叶根部形态反过来影响旋翼的气动弹性特性,这种旋翼/机身耦合问题,也是近年来直升机气动弹性问题研究中的重要方向和热点之一。此外,随着旋翼流场数值分析方法的日趋成熟,采用动态重叠网格或滑移网格方法来实现桨叶运动,并通过动网格技术来实现桨叶的弹性变形,从而实现弹性旋翼流场的数值模拟,目前正呈现出勃勃生机,成为直升机气动弹性研究的又一重要方向和热点。随着各种新构型直升机的相继出现,如倾转旋翼机、前行桨叶概念旋翼(ABC)直升机和复合式直升机等,也带来了新的气动弹性问题。不断发现问题、解决问题,推动本学科持续发展,永远是气动弹性工作者终身奋斗的目标。     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

Aeroelastic stability analysis of a bird-damaged aeroengine fan assembly

Bird strike is a major consideration when designing fan blades for large-diameter aeroengines. Current methods rely on impact tests and structural optimisation but it is highly desirable to have predictive numerical models to assess the aerodynamic and aeroelastic stability of bird-damaged fan assemblies. The aim of this paper is to present such a methodology and to study a representative case. The particular fan assembly under investigation contained two consecutive blades with unequal impact damage, the so-called heavy-damage and medium-damage blades. A detailed finite element analysis of the dynamic behaviour revealed that the vibration modes were significantly different from those of the tuned assembly. The twin modes were found to be split into single modes, some with highly distorted modeshapes, the so-called rogue modes. A nonlinear viscous flow analysis revealed truly unsteady effects and time-accurate aeroelasticity analyses with vibratory blade motion were undertaken to investigate the flutter stability. The computational domain included both a whole-annulus fan assembly and an intake duct and the resulting mesh contained approximately 2,200,000 grid points. The investigation was conducted for two points on the compressor characteristic, the first one corresponding to higher mass flow/lower pressure ratio and the second one to lower mass flow/higher pressure ratio. At the higher mass flow point, the flow separation was restricted to the immediate surrounding passages and the forcing onto the downstream blades was relatively small. However, a rotating stall event was observed for the lower mass flow point and the subsequent unsteady aerodynamic forces on the blade were high. At both mass flow settings, the flutter stability of the damaged fan assembly was predicted to be worse than that of the undamaged reference assembly.     

低雷诺数下横流-射流中剪切涡的试验研究

为深入分析横流-射流(JICF)的流动特性及其中的复杂涡系结构,从流动机理上研究燃机叶片气膜冷却,揭示高温燃气流与冷却流的掺混机理,本文对横向流中单孔射流所形成的剪切涡开展了试验研究。主要研究了速度比、雷诺数及射流角对JICF所形成剪切涡的影响。结果表明:速度比、雷诺数以及射流角会改变主流与射流之间的掺混程度,从而改变射流轨迹的曲率、高度及垂向渗透能力,最终改变剪切涡的特性;迎风涡与背风涡分别是由射流边界层涡与主流边界层涡形成的,当主流边界层涡强度大于射流边界层涡时,背风涡是流场的主导结构,反之,迎风涡将成为流场的主要涡系结构。     

含裂纹叶片的轴流式压气机整体叶盘振动特性分析

某型轴流式压气机性能试验过程中出现转子叶片裂纹故障.结合该型压气机结构特点,分析了振动响应的时域信号、频域信号以及高频成分出现的频率波动现象,总结了含裂纹叶片的轴流式压气机整体叶盘振动特性.分析结果显示:含裂纹叶片的整体叶盘使得轴流压气机振动响应呈现非线性特性;振动频谱中出现共振峰分离现象,表明叶片裂纹引起轴流压气机1级转子叶盘出现结构失谐;随转速升高,两个共振峰对应频率均出现波动现象(波动周期分别约为8s和9.2s),叶片裂纹不断扩展且数量增加,整体叶盘失谐程度不断加深,整体叶盘系统出现低阶模态局部化现象.      

具有分流叶片跨声速离心叶轮激波及能量损失分析

对具有分流叶片跨声速离心叶轮,采用激波基本关系式对数值信息进行分析和判定,理论证实了叶轮流道中激波的存在,并确定流道内高速气流易在主叶片前缘、分流叶片前缘以及主叶片的吸力面与分流叶片的压力面通道位置形成激波.进一步对比分析了叶顶间隙和运行工况对激波位置分布特征的影响及引起的能量损失.结果表明:固定叶顶间隙下,随着流量减小,激波位置向流道上游移动;叶顶间隙导致激波线位置向流道下游偏移;分流叶片前缘处有最强激波并引起较明显能量损失.      

减压器关闭过程内部流场的动态仿真和特性分析

针对自主研制的大流量减压器,采用动网格技术和流固耦合仿真技术,将气体的流动和活动组件的运动耦合起来,实现了关闭过程中内部流场的动态仿真,描述了关闭过程减压器内压力和速度变化过程,掌握了流场中波系、回流区的演化情况。同时,通过对流场中不同位置参数进行监测,分析了监测点各参数随时间的变化情况,以及流量和开度的实时变化情况,研究了基于动态流场的关闭过程动态特性,揭示了影响动态特性的内在因素。     

叶尖泄漏掺混损失影响因素分析

借鉴Denton不可压叶尖泄漏损失的理论分析,推导出可压缩流动中叶尖泄漏掺混损失基本关系式,其集中反映了尖区叶表压强、叶尖间隙、泄漏流流量系数及总损失系数影响.在仅考虑泄漏掺混损失情况下,数值分析了叶尖基元负荷、间隙、泄漏流量系数整体水平及其弦向分布等因素对叶尖泄漏掺混损失的影响.结果展示了各因素对叶尖泄漏掺混损失的作用及机理.最后提出一个重要的组合参数,反映了叶片尖部载荷和叶尖间隙的流向分布对叶尖泄漏掺混损失的综合影响.