直接供气预旋转静系流动和换热数值模拟

采用CFD商业软件FLUENT分别从维度模型和湍流模型两个角度对接近真实的高压涡轮旋转盘腔进行数值模拟,研究了一种直接供气预旋转静系旋转盘腔模型内的流动和换热特性.通过计算研究发现:低位腔内的流动仅由中心入流控制,二维计算不能反映高位腔内中心入流与预旋入流的掺混;二维和三维计算所得的腔内压力分布相似,实验值介于两者之间;主盘面的换热主要受中心入流控制,二维和三维计算差别不大,Realizable k-ε湍流模型在定性和定量上均能更好地反映换热实验结果.      

湍流模型对三维翼梢涡流场数值模拟的影响

 梢涡流场数值模拟对减少梢涡的不良影响意义深远。应用雷诺平均Navier-Stokes(RANS)方法建立了三维翼梢涡流场数值模型,选用的湍流模型包括两种经典两方程模型和一种代数雷诺应力模型。计算内容包括翼梢涡流场速度分布、压力分布及涡核位置,为了准确模拟梢涡涡核内系统旋转和流线曲率影响,将旋转和曲率修正方法应用于选用的湍流模型中。计算时对翼梢涡涡核处网格进行了加密处理。将计算结果与实验结果对比后发现,修正后的湍流模型较标准湍流模型明显提高了预报精度,与实验结果吻合较好,具有工程应用价值。     

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

应用二阶雷诺应力封闭的湍流旋流计算

本文叙述了三种不同的湍流旋流流动有限差分的计算结果,即:涡流管中的流动、强自由旋转射流、轴对称,旋转突然扩张的湍流流动。湍流的模化是求解微分形式的雷诺应力分量输运方程的二阶封闭。计算的结果均与试验数据作了比较。比较的结果表明雷诺应力模型明显地优于常用的K—ε模型。而预示与试验数据的吻合程度也是令人满意的。     

轴流压气机子午流面计算中三维剪切作用的湍流模型

本文发展了一个轴流压气机子午流面二维计算中模拟三维剪切结构对子午流场作用的湍流模型。试验与这一模型进行了初步的比较。模型表明,在多级轴流压气机子午流场计算中,必须计入实际流场三维剪切结构对子午流场的湍流扩散作用。     

飞机增升装置的数值模拟研究

增升装置的设计是飞机设计的一项极其重要的内容.本文通过多块网格技术生成计算网格,运用全湍流SST(Shear stress transport)模型和带有转捩修正的SST湍流模型对飞机三维增升装置的气动力特性和流场特征进行了数值模拟研究,通过对数值计算结果与实验结果的对比分析,研究了两种湍流模型的特点及其在飞机增升装置...     

跨音速单转子压气机三维湍流流场的数值计算

用有限体积表示的三维雷诺平均N-S方程显式格式的时间推进法, 对一个单转子压气机的三维湍流流场进行了计算。方法中应用了Baldwin-Lomax湍流模型。这个单转子压气机是第39届国际燃气轮机会议作为考核三维流场数值方法的盲题。计算结果与实验实测的结果进行了比较, 证明本文提出的方法可以用来验证叶轮机械气动设计。      

双方程k-ωSST湍流模型的显式耦合求解及其在叶轮机械中的应用

双方程k-ω剪切应力输运(SST)湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax(BL)模型、Spalart-Allmaras(SA)模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

Flow and Heat Transfer Characteristics in Rotating Two-pass Channels Cooled by Superheated Steam

In a modern gas turbine,using superheated steam to cool the vane and blade for internal convection cooling is a promising alternative to traditional compressor air.However,further investigations of steam cooling need to be performed.In this paper,the three-dimensional flow and heat transfer characteristics of steam are numerically investigated in two-pass square channels with 45° ribbed walls under stationary and rotating conditions.The investigated rotation numbers are 0 and 0.24.The simulation is carried out by solving the Reynolds averaged Navier-Stokes equations employing the Reynolds stress turbulence model,especially considering two additional terms for Coriolis and rotational buoyancy forces caused by the rotating effect.For comparison,calculations for the air-cooled channels are done first at a Reynolds number of 25 000 and inlet coolant-to-wall density ratio of 0.13.The results are compared with the experiment data.Then the flow and heat transfer in steam-cooled channels are analyzed under the same operating conditions.The results indicate that the superheated steam has better heat transfer performance than air.Due to the combined effect of rotation,skewed ribs and 180° sharp turn,the secondary flow pattern in steam-cooled rotating two-pass channels is quite complex.This complex secondary flow pattern leads to strong anisotropic turbulence and high level of anisotropy of Reynolds stresses,which have a significant impact on the local heat transfer coefficient distributions.     

基于非接触式测量的旋转叶片动应变重构方法

基于叶端定时非接触式测量和振动响应传递比的概念,开展高速旋转叶片动应变重构方法的研究。在频域内推导了叶片任意测点位移与任意测点动应变的传递比,给出了单模态共振下响应传递比关于位移和应变模态振型的解析表达式;建立旋转叶片的三维(3D)有限元模型,开展考虑旋转预应力效应的叶片模态分析,提取位移和应变模态振型,获得任意转速下叶端位移与叶根关键点动应变的传递比。开展高速旋转叶片叶端定时非接触式测量实验,采用周向傅里叶算法对叶端定时信号进行处理,获得叶片在不同转速单模态共振下的叶端位移,结合响应传递比,重构5个旋转叶片的关键点动应变。结果表明:旋转叶片在9000r/min和13000r/min转速下发生1阶共振时,与应变片实测结果相比,叶根处应力最大点、次大点和边缘点3个关键点的动应变平均重构误差均小于15%,验证了旋转叶片动应变重构方法的有效性。     

叶轮机械三维粘性动静叶干涉的数值模拟

本文发展了 1套模拟叶轮机械三维粘性动静叶干涉的数值方法。在成熟的空间离散方法基础上 ,提出了 1种高效的双时间步隐式时间推进方法 ,并构造了一种保证守恒的动静叶交界面处的数值通量计算方法。本文对某三维亚音速涡轮的动静叶干涉进行了数值模拟 ,结果与实验符合较好。     

旋转叶片频率转向与振型转换特性

通过对典型旋转叶片的有限元分析, 讨论了旋转叶片动频特性的“频率转向”和“振型转换”问题, 尤其是关于多模态耦合的复杂“频率转向”和“振型转换”的现象和特点.首先讨论了旋转叶片动频特性中的两阶模态间“频率转向”和“振型转换”问题, 说明了这种耦合振动问题的机理, 然后研究了多阶模态相互耦合所导致的复杂“频率转向”和“振型转换”现象, 对其基本特点进行了详细的分析和说明.      

复合材料风扇叶片-机匣碰摩振动的数值模拟

基于全三维模型的复合材料风扇叶片-机匣碰摩动力学特性的数值研究。建立了复合材料三维风扇叶片-机匣实体有限元模型,考虑叶片离心刚度的影响,机匣由贝塞尔曲面拟合,风扇叶片采用3种不同铺层形式。基于该模型,在机匣二节径变形的情况下,计算了不同铺层叶片在不同转速下与机匣碰摩后的动力学响应。计算结果表明:铺层形式对叶片中、低转速下的碰摩振动形式影响较大,带有90°铺层的叶片的最高振幅及不稳定振动所在区域的转速相对较低,改变铺层形式可以对复合材料风扇叶片-机匣的碰摩动力学特性加以控制。当转速靠近由3倍频与叶片1阶模态造成的共振点附近,或由6倍频与叶片2阶模态造成的共振转速附近时,叶片与机匣的碰摩会导致非稳定振动的产生。该方法与结果对复合材料风扇叶片的碰摩动力学特性研究具有一定指导意义。      

SST湍流模型在高超声速绕流中的改进

为模拟高超声速湍流问题,对剪切应力输运(SST)湍流模型系数进行了修正。数值格式采用改进的总变差递减(TVD)格式,并对湍流模型的负值强制项进行了隐式处理。在此基础上计算了绕平板以及具有分离、再附、激波/边界层干扰等复杂流动结构的压缩拐角的高超声速流动。计算结果与试验数据及半经验公式的对比表明:SST湍流模型引入的雷诺剪切应力与湍动能之比为常数(Bradshaw数)在高超声速绕流中并不成立。Bradshaw数修正后的SST湍流模型与原模型相比,所计算的壁面压力、摩擦阻力和壁面热流分布更接近试验结果。     

风力机叶片表面压力的计算与外场测试分析

研究在外场工况下,对风力机叶片表面压力的测试方法,并将测试结果与CFD计算结果进行比较。为了获取叶片在外场非稳态工况下的压力信息,沿叶片展向选取7个典型段面布置带式压力传感器。在数值计算中,通过数码扫描得到试验翼型的几何形状并建立计算模型,用非压缩的N—S方程和SSTk-w。湍流模型耦合,分别对7个翼型的气动性能进行计算。通过对试验和计算结果的对比发现,因为三维旋转效应的存在,基于动量叶素理论的二维翼型计算常常低估了实际风轮动力的产生,旋转叶轮表面压力分布和二维翼型计算结果明显不同。     

湍流方程扩散项离散对某径向叶轮流场计算的影响

在相对旋转坐标系下采用Harten,Lax and van Leer contact （HLLC） 格式离散对流项,自行开发了基于多块结构化网格的有限体积程序,实现了对叶轮机械内部流场的数值求解.分别对半开式径向叶轮、闭式后弯叶轮展开数值模拟,程序和商业软件计算得到的不同叶高处表面压力数据,其相对差异不超过1%,验证了算法的正确性.针对湍流方程的扩散项,分别使用完全离散和略去交叉导数项离散,通过对湍流黏度等值线、气动轴向力和力矩的比较表明:在网格正交性较好的情况下,略去交叉导数项的离散对计算结果的影响小于1%,显著地减小湍流方程离散的计算量.      

A three-equation transition model with mechanical considerations

A three-equation transition model based on the transition V-model is proposed for subsonic flows in this study. Considering the mechanical approximation of the generation process of the pre-transitional vorticities, the value of laminar Reynolds shear stress related to the mean shear deformation was calculated in the original transition V-model. Then a new transition model, named V-SA model, was proposed, which considered the phenomenological process of transition and presented great results for flows with and without pressure gradient. It is well-known that the baseline Shear Stress Transport (SST) turbulence model shows excellent performance of accuracy and robustness in plentiful flow cases, but it is important to predict boundary layer transition. The current model (V-SST) successfully couples the V-model to the SST turbulence model by introducing the effective turbulent viscosity and additional correction terms into the transport equations. A thorough evaluation of its ability to predict transition features is performed versus the well-documented flat plate of ERCOFTAC, including T3A and T3B without pressure gradient, T3L2 and T3L3 with semi-circular leading edge, the three-dimensional 6:1 prolate-spheroid under two angles of attack, and the NLR-7301 airfoil under different Mach numbers. Numerical results show that the current model has an attractive and superior performance in the simulation of boundary layer transition processes     

基于响应传递比的转子叶片动应变重构

提出基于响应传递比的转子叶片动应变重构方法,从系统应变频响特性出发,在频域获得应变-应变响应传递比关于模态振型的解析表达式,建立叶片已知测点应变与不可测位置应变的映射关系;开展高速旋转叶片动应变测量实验,建立叶片有限元模型,并采用转子叶片前4阶固有频率的实验结果修正有限元模型;开展考虑转速影响的转子叶片模态分析,提取对应的应变模态振型,计算得到单模态振动下应变-应变响应传递比;综合应变响应传递比和实测动应变重构叶片不可测位置的动应变。实验结果表明:与应变片实测结果相比,基于响应传递比的转子叶片动应变重构结果相对误差小于10%。