叶轮机械颤振稳定性的工程预测方法

通过求解雷诺平均N-S方程并采用影响系数法,对三个基本正交模态的三维振荡涡轮叶片绕流问题进行了研究。结果表明数值模拟所得非定常气动力与实验吻合较好;基于刚体运动假设和模态叠加法,得到了可用于叶轮机械设计阶段颤振稳定性评估的稳定性参数图/气动阻尼图,即所有可能的刚体运动及其相应气动阻尼的关系图;通过对稳定性参数图/气动阻尼图的分析表明,在一定情况下振型对颤振稳定性有重要影响,应将其作为颤振稳定性设计的重要参数之一。     

单转子风扇的三维反问题气动设计

采用准三维流函数反问题设计方法与三维 N-S求解方法的相互迭代对单转子风扇进行气动设计。在叶片初步计算中得到叶栅进出口气流角沿径向分布 ,并将它作为本文的目标函数。采用准三维反问题求解方法 ,依次构造出各个 S1流面上的叶片几何形状和气流角分布。然后再采用 N-S方程的求解方法 ,对叶片进行全三维流场的数值计算。通过 N-S方程计算结果与目标函数的对比 ,重新修正叶片出口气流角分布 ,并作为下一次反问题设计的目标函数。经过反问题与 N-S方程求解的反复迭代 ,最终得到满足设计要求的叶型。     

压气机叶片流固耦合数值计算

采用流固耦合方法对压气机单排跨声叶片进行气动弹性数值模拟.首先对压气机流场进行三维定常求解,并与试验结果进行比对.在定常流场分析的基础上,进行流固耦合计算,采用非定常方法分别求解流体域Navier-Stokes(N-S)方程和固体域结构动力学方程.耦合过程中,流体域和固体域在物理边界上进行能量传递.计算结果表明:所采用的流固耦合数值方法用于判断叶轮机械叶片气弹稳定性是可行的,可以得到叶片的瞬态响应,从而判断叶片是否发生颤振.      

叶轮机械颤振稳定性工程预测方法在跨声风扇中的进一步探讨

通过求解雷诺平均Navier-Stokes(N-S)方程并采用影响系数法,对三个基本正交模态的三维振荡跨声风扇叶片绕流问题进行了研究,基于刚体运动假设和模态叠加法,得到了可用于叶轮机械设计阶段颤振稳定性评估的稳定性参数图.结果表明,跨声风扇内部,激波对非定常气动力的分布具有主导作用;在一定情况下振型对颤振稳定性有重要影响,应将其作为颤振稳定性设计的重要参数之一;跨声风扇和亚声低压涡轮的稳定性参数图对比分析表明,两者稳定性参数分布形式具有相似性.      

求解叶轮机械三维粘性可压流的初步研究

给出了用“参数多项式方法”求解叶轮机械三维粘性可压缩流的基本思路、方程和方法，采用圆柱坐标系下的三次参数多项式。中还应用此方法以一个透平和一个压气机叶片为例计算了三维粘性流场。     

T型尾翼民机深失速气动特性的计算研究

以非结构混合网格上的雷诺平均N-S方程(RANS)求解方法为基础,对采用T型尾翼的某民机外形在大迎角深失速情况下的粘性绕流进行了数值模拟,获得了该机在巡航和降落外形在深失速情况下的气动特性变化规律.通过对T型尾翼深失速情况下的流动特点进行分析,解释了T型尾翼布局飞机在深失速情况下俯仰力矩系数会随迎角出现不利于纵向稳定性的变化趋势的原因.     

钝体前激波的数值模拟

本文通过直接求解N-S方程,研究激波问题,对是否可以和如何用直接求解N-S方程来处理超跨声速流进行了探索,在对粘性项的处理建立可压缩流体的压力校正方程和计算方法等方面均提出并实施了新的方法,最终得到了基本合理的数值计算结果。     

基于CFD/CSD耦合的叶轮机叶片失速颤振计算

为了研究叶轮机叶片的失速颤振特性,发展了一种计算流体力学与计算结构力学(CFD/CSD)时域耦合方法。该方法通过每一物理时刻CFD和CSD的循环迭代实现了耦合计算。在CFD分析中,采用鲁棒性较好的空间离散格式AUSM+-UP,并基于延迟脱体涡模型(DDES)模拟了带分离流动。在结构分析中,通过模态法构建了旋转叶片动力学方程并运用杂交多步方法进行求解。以孤立转子Rotor37为例,计算了不同工况下流场总体与细节参数,与实验结果的对比验证了CFD算法的精度。对某转子叶片进行了颤振特性研究,计算所得的广义位移时间响应曲线表明该叶片在近失速工况下会发生失速颤振,其表现形式为一阶弯曲模态发散且各阶模态之间不耦合。分析表明,流场不稳定和非定常效应是引起失速颤振的关键因素,同时折合频率的降低也会导致原本气动弹性稳定的叶片发生失速颤振。     

跨声速风扇的流固耦合的颤振分析(英文)

发展了时间推进的叶片颤振的数值方法,采用了时间推进求解流体和固体相互作用的过程。气动模型是基于求解三维雷诺平均的N-S方程,采用了多块结构化网格的有限体积格式,对流和耗散通量的计算使用了二阶迎风格式和中心格式。在变形的动网格上流体运动的守恒型方程的求解采用双时间步,隐式格式和多重网格方法。叶片振动采用了振型叠加的线性气动弹性模型。该方法在气动弹性标准算例4进行了验证,并用于求解跨音速风扇的颤振问题。     

重型燃气轮机压气机叶片气动弹性稳定性初期预判

通过求解雷诺平均N-S方程并采用影响系数法,对某重型燃气轮机压气机叶片三个基本正交模态的叶片绕流问题进行了研究;同时结合有限元分析和稳定性参数分析,对压气机叶片气动弹性稳定性进行了初期预判。结果表明:以振型为判断标准的稳定性参数图,可作为压气机设计中气动弹性稳定性的预判工具;本文所述颤振稳定性初期预判方法,可为重型燃气轮机压气机气动弹性稳定性初期设计提供技术支持。     

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

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

分离流动的特性分析

本文分析了三维分离流的特性,且从基本物理现象出发,以数值求解NavierStokes方程所给出的三维粘性分离的结果和三维分离流的实验结果为依据,给出了三维粘性绕流产生分离的必要条件和物面上分离线的定义。     

联合应用现代优化方法与S2流面计算进行多级涡轮的气动设计

本文联合应用S2流面正问题计算和多级局部优化设计对某三级涡轮进行多级气动优化设计。优化联合采用人工神经网络和遗传算法,流场计算采用全三维粘性流N-S方程求解,计算网格采用H-O-H型网格,即入口段、出口段采用H型网格,叶片区域采用O型网格。通过优化,总效率提高1.1%,总体性能提高,达到设计要求。     

基于伴随方法的机翼多设计点气动反设计方法

针对三维机翼多点多约束气动反设计问题,基于伴随理论方法和粘性流雷诺平均N-S方程,通过粘性流数值模拟、伴随方程与梯度精确数值求解、计算网格高效算法及梯度类优化算法等有效结合,并采用考虑多设计点梯度权重系数的并行计算近似模式,开展了一种三维机翼多设计点多约束气动反设计方法研究,进行了典型算例验证。研究表明:所发展的机翼多设计点气动反设计方法具有较好的鲁棒性及优化效率。     

振荡翼型非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序.并借助运动网格生成方法,对周期性俯仰运动下NACA0012翼型的非失速绕流和动态失速问题进行了数值模拟研究,分析了振荡翼型的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法在模拟非失速绕流时是可靠的;结果还表明,随着振幅、折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.在模拟动态失速问题时.本文方法在深度失速下结果与实验符合较好,而在轻度失速下符合较差.经分析表明这是不同条件下非定常流动对湍流模型和转捩的敏感程度不同所致.     

失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

叶片失速颤振发作机理的探讨

 本文对于现有的叶片失速颤振计算模型进行了一些改进,以期使其更接近于真实物理图画。计算结果表明,就槽道激波以及物面气流分离对于叶片气动弹性稳定性所起的作用来看,在单自由度纯弯振荡和双自由度弯曲主导振荡两种情况之下是具有实质性差别的。在文中对于一个实际转子的失速颤振边界进行了数值预测,并与相应的实验结果进行了对比。     

二元喷管三维粘性流场的数值模拟

本文用数值方法对二元喷管的三维粘性流场进行模拟,将 Mac Cormack有限差分格式拓广到全三维 N-S方程的求解。在数值迭代中使用加快收敛速度的当地时间步长,并根据格式压缩计算机存储量和计算量。进行了两个算例的数值模拟,计算结果与实验相符。      

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。