S_1流面H型与C型双重网格计算

本文发展了一种S_1流面的H型和C型双重网格计算方法,即首先采用H型网格进行全流场计算,再在叶片前缘附近采用局部C型网格进行加密计算,在生成C型网格时采用了一种组合二次曲线拟合方法以便适用于各种压气机及透平叶型.文中给出的算例表明,本方法是提高S_1流面亚声速绕流计算精度和预测叶栅变工况性能的有效手段。     

后掠冀跨声流场S_1/S_2流面全三元迭代解

一、前言 叶轮机械三元流动通用理论提出的使用两类流面迭代求解叶轮机械三元流动正反问题的方法,已成为内流数值计算的一个基本方法。随着计算机速度的提高、内存的扩大,越来越多的三元直接解投入了使用。但是两类流面迭代求解仍然是叶轮机械内流计算的一个重要方法。     

帧转移型数字相机帧转移模糊的复原

帧转移型面阵CCD相机因其诸多优点而得到了广泛的使用 ,但它也存在帧转移模糊的问题。介绍了CCD相机的几种面阵类型及其各自的一些特点 ,重点介绍了帧转移型面阵器件的工作过程 ,分析了面阵CCD器件帧转移模糊形成的原因 ,并提出了解决帧转移模糊的方法。介绍了用该方法消除帧转移模糊的实例 ,验证了方法的有效性     

非正交曲线坐标S_1流面流函数反问题松弛计算

本文在文献[1—3]工作的基础上,从叶轮机械S_1流面反问题提法之一(给定叶栅吸力面速度分布及叶片厚度分布求解叶型坐标)出发,推导了流函数反问题主方程及有限差分方程.这方程是以计算网格坐标为主变量的二阶偏微分形式的动量方程,解决了文献[4—7]所未能解决的使用有旋的运动方程求解的问题.此方程与有旋的S_2流面流函数方程的一致性保证了叶轮机械三元求解的收敛性.进一步完善了叶轮机械使用两类流面的三元流设计方法.编制了计算机程序对典型的叶型作了计算例子,结果是理想的.     

离心压气机转子内部流场S_1／S_2全三元迭代解

为了比较准三元解和全三元解的差异，验证准三元解在计算离心压气机转子内部流场的准确程度，研究离心压气机转子内部流场全三元流动特性，本文对一有激光测量结果的高压比离心压气机［１］叶轮内部流场进行了全三元迭代计算，分析了Ｓ１／Ｓ２两类流面在叶轮通道内分布形态，比较了两类流面准三元解与全三元解的计算结果，讨论了无粘二次流的分布。并进一步和激光测量值及Ｎ－Ｓ三元直接解进行了详细的比较。     

Chauvenet准则和Laeda准则的分析

一、引言Chauvenet 准则和 Laeda 准则,都是判断实验数据的异常值是否包含粗差的常用统计方法。可是,对同一个样本,这两种准则的判断结果却往往相反,例如,用 Cha-uvenet 准则去判断文献[1]、[2]中那些正态样本的异常值,均得出样本的其余观测值是来自一个正态母体,而样本的异常值却来自另一个正态母体;而用 Laeda 准则去判     

使用非正交曲线坐标与速度分量S_2流面跨音速流函数解

基于叶轮机械两类流面迭代计算理论,在非正交曲线坐标上建立了S_2流面上弱守恒型流函数方程.使用人工密度修正方法求解S_2流面跨音流动正问题,用速度积分方法避免了密度双值问题,并编制了相应的计算机程序.     

使用非正交曲线坐标与非正交速度分量的跨声速S_2流面势函数方程人工可压缩性迭代解

本文推导建立了适于求解跨声速轴流式压气机转子中S_2流面正问题的全守恒型势函数方程,方程的求解采用人工密度的方法和一种新的Φ-ρ(Γ)迭代方法,能在S_2流面上自动捕获激波.用本方法编制的计算机程序对西德DFVLR单级跨声速轴流式压气机转子的一个最高效率点实验工况进行了验算,并将计算结果与实验结果作了比较。     

翘曲S_1流面跨声速流函数解

一、前言 在叶轮机械工程设计、计算中,往往使用只计算一个中心S_1流面和若干个迴转S_1流面的准三元迭代解。为了得到更准确的全三元解,文献[1]在全三元迭代计算中使用了翘曲的S_1流面计算机程序。文献[3]则发展了使用曲面拟合方法的翘曲S_1流面程序。在跨声流动存在强烈激波间断时,流面形状会在激波处发生折转,流片厚度也会突变。由于这种三元效应的存在,有必要发展任意翘曲S_1流面跨声程序,进行全三元跨声迭代解。本文在文献[5]的基础上发展了翘曲S_1流面跨声计算机程序。     

CAS压气机转子跨声速流场S_1/S_2流面全三元迭代解

作为跨声速单级压气机基础研究工作的一部分,本文在过去完成的跨声速流场准三无迭代解的基础上,发展了在任意翘曲的S_1流面和若干个S_2流面上跨声流函数解相互迭代的全三无计算程序体系。对CAS跨声单转子内部流场进行了验算。经过七轮迭代,第一次得到了三元跨声流场(不假定无旋)全三元迭代解,验证了跨声速全三元迭代计算的收敛性。     

透平S_2流面气动半反命题有限元解

一、变分原理 吴仲华教授的三元流动理论已在国内外透平机械计算中获得了广泛的应用。但基于吴氏理论的变分有限元计算实属少见。近几年来刘高联基于两个流面理论对透平机械中多种气动命题找到了泛函表达式,并编制了一些有限元程序。但S_2流面的计算尚少。 基于S_2流面半反命题的气动基本方程组,引入流函数后可得下列主方程:     

WP6涡轮部件设计状态及实验状态的数值模拟

本文对WP6涡轮部件按照折合转速及膨胀比相等的准则设计了两套实验方案,并对设计状态及设计的两套实验方案进行了S2正问题计算和S1流面计算。计算表明,所没计的实验方案与设计状态马赫数相似,达到了实验设计的目的。为校验所用S2正问题计算和S1流面计算程序,对设计状态进行了S2正问题计算、S1流面计算及FLUENT计算,总参数如流量、效率等吻合较好,叶表马赫数分布趋势也很相似。     

使用非正交曲线座标与速度分量S_2流面反问题流场线松弛解

针对叶轮机械S_2流面反问题的计算,介绍了使用任意非正交曲线座标和非正交速度分量的S_2流面反问题流场线松弛解法计算机程序.并对该程序与目前一般常用的速度推广法(流线曲率法)程序和矩阵直接解法程序的不同之处,作了简要的比较和评论.指出流场线松弛法的优点,特别是在采用叶片三维设计计算方法时,当叶片区沿流线方向必须设立更多的计算站时,流场线松弛方法是一个值得推荐的好方法.     

非正交曲线坐标叶轮机械S_2流面正问题计算机程序

对于非正交曲线坐标叶轮机械S_2流面正问题(亚声速情况)在计算机上实现了一种收敛算法.对一个高载荷轴流式两级风扇的验算说明计算结果与实验数据相比很一致. 主导方程采用新推得的非线性形式——方程(8),与传统沿用的线性形式——方程(13)相比,当流场中高亚声速区较大时,新解法改善了数值求解过程的收敛性.     

对亚音速静子使用一个中心S_2流面与多个S_1迴转面的三元迭代解

根据五十年代初开始建立起来的三元流动通用理论,利用文献[2—5]提供的计算机程序,用一个中心S_2流面和几个S_1迴转面交叉迭代的方法进行了高亚音静子的三元气动设计(反问题)计算。在交叉迭代过程中,静叶前S_1流片的形状和厚度变化及前、后远方的边界条件是一个不易确定的问题。本文通过新近提出的一个满足静叶前间隙中规定气体参量的相当流场的计算方法得到S_1流片的形状和厚度以及远前方周向均匀的进气条件。在迭代过程中熵的处理的一致性也是一个需要考虑的问题。计算结果表明S_1计算中必须适当考虑熵增的影响,以便与S_2流场中的熵值相一致,并能得到较高的迭代精度。这样做后,本文经过四轮迭代,就得到了S_1和中心S_2对应点上的各项参数相对差别都在2％以内的收敛解,并同时得到了满足流量和厚度分布的叶型及其叶面上的压力分布和速度分布,因之是一个可以实用的三元流叶片设计方法。     

多级跨音速压气机S_2流面正反问题混合计算

多级跨音速压气机Ｓ＿２流面正反问题混合计算苏杰先，冯国泰，魏沫，胡昆（哈尔滨工业大学）一、前言大型航空发动机跨音速压气机的发展有力地促进了小型弹用发动机的飞速进步。美国的＂战斧＂式巡航导弹用发动机Ｆ１０７－ＷＲ－１００以及法国的弹用发动机ＴＲＩ－６０...     

S_1流面跨声流场流函数矩阵解

跨声速叶栅流的计算,可采用时间相关法求解Euler方程,或用松弛方法求解势函数方程和流函数方程。一般说来,时间相关法耗费机时较多,势函数方法仅对无旋流适用。流函数方法适用于二元有旋流的计算,并且边界条件也较为简单,可方便地进行S_1和S_2两类流面迭代得到三元解。流函数方法的跨声计算最大的困难是密度双值问题     

交替方向积分法解S_1流面跨音流函数方程

在求解S_1流面跨音流函数方程中,提出并采用了交替方向积分法,一方面避开了ρ的双值性矛盾,同时也使基本方程的关系在整个求解域中有可能均衡地趋于满足.本方法利用了非正交曲线座标和人工可压缩性技术,原则上既可计算亚音来流的跨音问题,也可计算超音来流的跨音问题。     

CAS压气机转子中跨声速三元流动的S_1/S_2流面准三元迭代解

本文对我所跨声速单级轴流压气机0.8和0.9设计转速下转子通道内及其上、下游的流场,进行了S_1、S_2两类流面迭代计算,分析了三元跨声速流场,并与使用双焦点激光测速仪的测量结果进行了比较。在分析计算中,以激光测量的出口绝对气流角分布,常规测量的出口滞止压力分布及通过扭矩测量值校正得到的转子效率为输入参数,进行了6轮迭代计算,得到了跨声流场的收敛解。计算结果与测量值较接近。     

液态金属Ga急冷凝固中微观结构转变的模拟研究

采用分子动力学方法对液态金属Ga的快速凝固过程进行计算机模拟跟踪研究,运用HA键型指数法和原子团类型指数法分析了金属原子Ga的成键类型和形成的基本原子团结构.发现:与通常的液态金属(如Al)相反,随着温度的降低,二十面体及与二十面体相关的1551, 1541键型数目明显下降;与立方体(fcc,bcc),六角密集结构相关的1421,1422键型数目明显增加;而与菱面体相关的1321,1311键型的数目却显著增加,逐渐占据优势.最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc,bcc)、六角密集(hcp)等团簇结构所组成的非晶态结构.这正是非晶态金属Ga的g(r)曲线分裂的第二峰的顺序为前低后高,而与非晶态金属Al的g(r)曲线(其分裂的第二峰的顺序为前高后低)明显不同的微观结构上的物理根源.