跨音速压气机叶栅流动有粘-无粘相互作用的迭代计算

 本文利用较为先进的附面层内法向速度分布等半经验相关式,发展了一种适用于大逆压梯度的附面层积分方法。将激波近似地用大逆压梯度代替,对湍流附面层进行了快速计算,再将求得的附面层参数与无粘主流计算通过建立适当的数学模型结合起来,反复迭代后将计算结果与实验数据及无粘计算的结果比较,可以看出有粘-无粘迭代对无粘计算结果有明显改进。     

用隐式多重网格法计算实用外形无粘绕流

网格展弦比的好坏严重地影响多重网格的效率。文中细致地处理了多重网格算法的细节,了网格展弦比的影响,发展了一种计算实用外形绕流的隐式多重网格解法,计算效率提高了５倍。     

叶轮机械叶栅流道几何处理的计算方法

叙述了叶轮机械叶栅流道网格生成前的几何处理的计算方法,通过将处理后的数据保存为PLOT3D文件格式,可以在其他CFD软件中直接使用,或在商业网格生成软件中快速获得叶栅流道的实体造型,与传统的在CAD类软件中进行交互式几何处理相比,可以大大缩短CFD计算周期。     

隐式多重网格法求解叶轮机械三维跨声速湍流流场

本文提出了一种高效率求解叶轮机械三维跨声速湍流流场的隐工多重网格算法。它以矢通量分裂格式为基础，采用ＦＡＳ型多重网格加速技术。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

跨声速叶栅粘流计算的多级Taylor－Galerkin有限元法

本文对Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限无法在收敛速度和稳定性两方面作了改进，与多级有限元法相结合，构成了多级Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法。基于Ｎａｖｉｅｒ－Ｓｔｏ－ｋｅｓ方程，研究了平面跨声速叶栅流动，并与实验结果作了比较。计算结果表明新方法具有较好的稳定性及较快的收敛速度。     

两种多重网格法求解径流叶轮机械可压流动

为了加快数值求解叶轮机械内部流场的收敛速度,并改善计算结果,在数值求解径流叶轮机械可压流动的计算中,采用了由微元、二重和三重网格逐级转入和同时转入的两种多重网格法。对德国宇航局Ｋｒａｉｎ Ｈ设计的离心压气机内部流场的数值计算表明,两种多重网格法的计算结果差别不大,尽管在收敛速度方面存在差异,但均与测量数据吻合比较好。     

用N-S方程计算三维叶栅流

发展了基于微元中心有限体积方法,采用雷诺平均N-S方程求解三维叶栅流动的程序,特别注意了格式中人工粘性的加法;程序采用Baldwin-Lomax模型。格式经过平板层流及湍流流动的校验,层流与Blasius解、湍流与实验经验公式及Wilcox的计算值符合很好,最后以Goldman涡轮导向器为例进行了叶栅流动计算校核,结果符合于实验以及Chima等人的计算。      

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

涡轮叶栅流动和传热耦合计算

通过求解三维雷诺平均N-S(Navier-Stokes)方程和带转捩模型的二方程SST(Shear stresstransport)湍流模型,完成了带简单冷却的MARKⅡ涡轮导向叶栅流动和传热耦合计算.计算结果与试验数据的对比表明,发展的带转捩模型的数值方法明显地提高了叶栅温度和外换热系数的计算精度.同时研究发现,湍流模型对叶栅表面压力分布影响很小,而对叶栅表面温度和外换热系数影响很大;在其它边界条件相同的情况下进口湍流度对壁面压力几乎无影响,但对叶栅温度和外换热系数影响较大.      

用多重网格法模拟三维流动

采用多重网格法模拟三维流动。控制方程在三维直角或圆柱坐标系中通过有限体积方法进行离散 ,SIMPLE算法对速度与压力进行耦合 ,亚松弛策略提高离散方程的收敛速度。在压力方程的求解过程中 ,界面速度应用由Rhie、Chow和Majumdar等人提出的动量方程插值法来剔除不合理的压力波动。三维流动算例表明采用多重网格法模拟三维流动可以得到良好的计算精度和收敛性。     

平面大头叶栅跨声速流动的一种有效解法

准确数值预估叶片头部表面的压力分布,一直是叶片气膜冷却设计中十分重要而又尚未圆满解决的问题。本文针对这个问题提出了一个高准确度、高效率的计算方法。 在文献[1]的基础上,本文从任意正交曲线坐标系中的基本方程出发,引入Von Misses变换,导出了相应的流线控制方程,并具体地在拟边界层坐标系中计算了RA叶栅流场,结果与实验符合甚好,据作者所知,本文方法是已见到的计算同类叶栅方法中准确度最高、花费机时最少的。     

二维可压粘性叶栅流场计算的预测修正迭代多重网格法

 本文中提出了一种解可压N-S方程的快速数值法。它由预测步,修正步,迭代步以及多重网格法构成。叶栅出口采用了一种使数值解收敛性进一步加快的技术,我们称之为定常无反射边界条件技术。本文应用了k-ε双方程紊流模型。 对若干叶栅在-13°～8°的攻角范围内进行了数值试验,并与对应的物理实验结果进行了比较。结果是满意的。     

平面叶栅跨音流动的数值计算及变工况性能的研究

提出一种隐式矢通量分裂差分格式并用来直接求解Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组。该方法避开了近似因子分解及矩阵运算，具有精度高、稳定性好、计算量少等优点。在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

涡轮叶栅端壁区流动的实验研究

本文在大尺寸低速开式叶栅传热风洞中对一种高压涡轮导向叶栅中的流场进行了实验研究。采用五孔针对5个雷诺数下的叶栅端壁区三维流场进行了测量,并用线簇和小球浮动法对5个工况的流动进行了流场显示。实验结果表明:马蹄涡压力面分支在向吸力面运动的过程中,破坏了来流附面层的结构,在马蹄涡压力面分支之后,叶栅通道中产生了一个新的从压力面到吸力面的新附面层,新附面层的厚度小于来流附面层厚度;三维流动区约占叶栅通道的40%;雷诺数的增大将增强端壁区的三维流动。从流场显示图片可以观测叶片吸力面靠近端壁的角涡形成与发展,以及吸力面上的三角形区域;流场显示的通道涡大小与流场测量结果吻合。本文的实验结果有助于分析端壁表面和叶片表面换热特性的形成机理。      

跨声速叶栅粘流计算的多级Taylor—Galerkin有限元法

本文对Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法在收敛速度和稳定性两方面作了改进，与多级有限元法相结合，构成了多级Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法。基于Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，研究了平面跨声速叶栅流动，并与实验结果作了比较。计算结果表明新方法具有较好的稳定性及较快的收敛速度。     

带有机匣处理结构的叶栅流场数值计算

完成了一种数值计算带有机匣处理装置叶栅流动的方法。计算采用三维、非定常、可压缩Ｎ－Ｓ方程组及ｋ－ε紊流模型方程，选择３种典型的处理结构，周向槽式机匣、轴向缝式机匣、轴向倾斜缝式机匣为算例，并引入参数周向平均的概念，以与稳态实验测量的结果做定量对比，结果发现符合较好。本文为分析机匣处理的扩稳机理提供了理论基础，也为工程设计机匣处理装置提供了新的数值工具。     

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

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