基于浸入式边界法的叶栅颤振数值模拟

基于浸入式边界法建立了求解振荡叶栅非定常流动的快速计算模型,并在模型中添加湍流模型,使之能够考虑更加接近实际流动的情况。具体的方法是求解目前为止最常用的雷诺时均Navier-Stokes方程并添加低雷诺数湍流模型(Lam-Bremhorst模型,k-ε模型的一种)进行计算,为了验证该方法的正确性,分别对层流边界层和湍流边界层进行了数值模拟,计算结果与布拉修斯解和壁面率吻合得很好,证明了湍流模型的可靠性。在此基础上,对高雷诺数条件下的振荡叶栅进行了数值模拟。结果表明,折合速度是影响叶栅振荡的重要因素,这与层流计算的结论类似。值得注意的是耦合过程没有生成贴体网格,减小了计算时间,可以准确快速地模拟真实的叶栅流动情况。     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

压气机转子叶片颤振边界的预测方法

通过发展的基于能量法的颤振数值预测方法得到了压气机转子叶片的颤振边界.将计算结构动力学(CSD)分析得到的叶片表面节点位移插值到耦合面的流体网格点上,并将设计的多层动网格技术应用于计算流体动力学(CFD)方法,实现叶栅振荡作用下的非定常分析,得到叶片表面的非定常气动功以及模态气动阻尼比.以某第一级压气机转子叶片为例,对...     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

叶间相位角对叶片颤振的影响

发展了一种基于能量法并计及叶间相位角(IBPA)影响的数值计算方法.建立了全环振荡叶栅模型,将每一个扇区分为叶片表面附近的可动域及其外围的固定域,采用有限元形函数将计算结构动力学(CSD)网格点位移的传递到计算流体动力学(CFD)的耦合面网格点上,求解了全环振荡叶栅由k-ε湍流模型封闭的Reynolds平均Navier-Stokes(RANS)方程.以NASA 67转子叶片为例,在设计转速下,基于能量法分析了叶间相位角对叶片颤振的影响.计算结果表明:叶间相位角对叶片颤振有显著影响.      

振荡叶栅三维非定常流动数值分析

以全三维、非线性、非定常欧拉方程计算作用于振荡叶栅的非定常气动载荷。欧拉方程的数值离散采用了具有时、空二阶精度的MacCormack显式分裂的两步格式,用时间推进方法得到其数值解。应用动网格技术捕获叶片的运动,保证叶片振动过程中网格边界与叶片表面保持一致。以相移周期条件处理具有不同相位的叶片振动问题,同时发展了伪边界形状修正方法以简化周期边界条件处理。     

航空发动机叶片颤振机理探究

航空燃气涡轮发动机中,由颤振引发的叶片疲劳失效是导致发动机机械故障的主要原因之一。为了研究叶片颤振这一非定常气动现象,在能量法的基础上采用弱耦合的计算方法,模拟了发动机风扇叶片在不同转速和不同出口反压条件下的振动,得到了叶片的颤振边界。其中流场与叶片建模用到三维线性插值法转化振动位移,流场动网格建模用到了多层动网格技术。计算表明,叶片表面气动功的分布与叶栅内的激波分布有很大的关系,激波在叶栅通道中的振荡是引起叶片颤振的原因之一。激波在叶栅中前后振荡,使得叶片表面上单位面积做正功的区域发生变化。做正功的区域扩大时,叶片表面气动功绝对值减小,从而气动阻尼比减小,更靠近颤振边界。     

浸入式贴体网格边界方法的对转桨扇非定常流动数值模拟

应用所提出的浸入式贴体网格边界方法和开发的CFD程序,进行了对转桨扇跨声速全周叶片流场的非定常数值模拟.介绍了程序流程、网格生成策略和求解器原理.分析了前、后桨叶尖涡、尾缘涡的非定常干扰对叶片表面流场的影响.统计了对转桨扇的推力、功率和推进效率的非定常波动,并同NUMECA模拟的定常结果进行了对比,结果发现两种模拟得到...     

错频对叶片的气动弹性稳定性影响

采用求解动网格下的非定常可压缩Navier-Stokes方程模拟了振荡叶片下的气动弹性问题,研究了叶片的频率错频对叶片气动弹性稳定性的影响.通过数值模拟平面叶栅的气动弹性第10标准算例,验证了气动弹性的数值模拟方法,计算了不同叶片间相位角和折合频率下的气动阻尼系数,研究了叶片振动频率改变对叶片气动弹性稳定性的影响.计算结果说明:频率错频是提高气动弹性稳定性的有效方法,其主要作用是减小叶片间振动的耦合效应和叶片间相位角的影响,并且随着错频量增大叶片稳定性增强;通过模拟三维环形叶栅的气动弹性第4标准算例,计算了气动阻尼系数随错频量变化的规律,验证了错频量和气动弹性稳定性增强的规律.      

风扇/压气机叶型厚度对颤振特性的影响

采用计算流体力学与结构动力学相结合的方法,数值模拟了三维振荡叶栅非定常黏性流场;通过对叶片表面非定常气动力及其所做非定常气动功的计算分析,采用能量法对叶片颤振与否进行预估判断。针对三维直叶片和风扇转子叶片,通过调整叶片的相对厚度,研究了叶片厚度变化对风扇/压气机颤振特性的影响。此外,通过对振荡叶栅非定常流场结构的研究,发现了叶片吸力面气流分离与叶片振动之间的耦合关系。本文的研究在风扇/压气机设计中,可用于评估最大相对厚度等叶型结构设计参数对气弹稳定性的影响,对叶轮机颤振机理研究具有一定的参考意义。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动.该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递.分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法.流场...     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。     

吸气／喷气叶栅的三维贴体网格生成技术

对某吸气/喷气叶栅的复杂的多联通域边界,采用孔斯双参数样条曲面模型,通过一维问题边界条件控制疏密度,选择合理的参数来强调正交性;分别在双参数域和绝对坐标系下求解微分方程,得到二维、三维网格。应用上述方法,为某叶栅叶片前缘和端壁开多排孔方案生成了三维贴体网格,获得了理想的效果。     

基于高阶间断伽辽金数值方法的浸没边界法

为了降低复杂边界条件下贴体网格生成难度，近年来基于笛卡尔网格的浸没边界法逐渐成为研究此类问题的主要数值方法之一，然而计算精度和计算效率仍然是此类方法目前面临的挑战。与传统的空间2阶精度的有限体积格式相比，空间精度为3阶或3阶以上的高精度方法具有空间精度高、数值分辨率高、数值耗散小的优点，而作为高精度数值方法之一的间断伽辽金方法在浸没边界方面的应用仍较少。本文将高阶间断伽辽金方法的高精度优势与浸没边界法无需贴体网格的优势结合起来，提出了适用于可压缩流动的高精度浸没边界法。其中边界条件采用体积惩罚方式实现，同时采用牛顿法迭代以及MPI并行提高计算效率，物面的数据重建采用插值点处的逆距离权重（IDW-IP）方法代替高阶格式下多项式插值方法。本文基于笛卡尔网格测试验证了二维定常和非定常情况下的数值模拟效果，并与传统贴体网格的计算结果进行了对比。     

三维振荡叶栅内部非定常流动数值模拟研究

进一步发展了非定常雷诺平均N-S方程求解程序,并基于影响系数法对三种不同振型下的三维涡轮振荡叶片绕流问题进行了数值模拟研究.数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;振型对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,其中叶片运动导致叶栅各通道中流量的再分配以及出口气流角的变化等对叶片的负荷分布会产生较为显著的影响.     

一种新的适用于叶栅绕流计算的网格生成法及其应用

本文对平面叶栅气动计算网格生成方法作了简要评述。在此基础上,根据叶栅槽道流动的特点,提出一种以等比级数为基础,用几何方法构作H网格的数值方法。用这种网格进行流场计算的结果表明,所生成的网格在叶片前后缘区有很高的分辨能力,能够算出速度分布的细微结构。本文的计算结果与亚音速及跨音速叶栅速度图解符合得很好。     

叶栅稠密度及进出口气流角对反推性能的影响

导流叶栅几何及气动参数的选取直接关系到叶栅式反推力装置的性能,本文选取冲击式对称叶栅进行数值模拟,采用混合网格生成技术,研究叶栅稠密度、叶片进口气流角和出口气流角对轴向反推力系数和流量系数的影响。计算结果表明,叶栅稠密度在一定范围内变大时,流量系数减小,轴向反推力系数增加;在反推气流不被进气道吸入的情况下,当气流角增大,流量系数变大,轴向反推力系数降低。双曲非等厚叶片的气动性能要好于单曲等厚叶片。     

非失速二维振荡叶栅非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序,借助快速网格生成方法,对周期性俯仰运动下NACA0012叶栅的非失速绕流问题进行了数值模拟研究,并分析了振荡叶栅的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法是可靠的;结果还表明,随着振幅,折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.