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

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

一种叶轮机三维叶型黏性反问题方法

提出了一种应用于叶轮机三维叶片气动设计的黏性反问题方法,该方法将叶片表面边界视为是运动的壁面,通过建立并求解壁面法向动量方程,从而可根据给定的压力分布或压力面与吸力面之间的压力差分布直接求解叶型几何.结果表明:该方法可与定常雷诺平均Navier-Stokes(RANS)模拟配合进行求解,鲁棒性强、计算耗时短,且符合真实黏性流动环境,无需对叶型进行黏性修正.算例验证结果显示,该方法在单排和多排叶片环境下都能获得反问题解,计算稳定性和收敛性良好.      

基于边界涡量流的二维叶型优化建模与求解

针对二维叶片整体造型优化问题,提出一种基于边界涡量流的二维叶型优化模型,并对该模型进行求解。建立以含参数的Bezier曲线定义的二维叶型型线为边界、以边界涡量流正峰值的最小化为目标函数的优化模型,通过分析边界涡量流的生成机制,把叶型型线边界信息代入边界涡量流中,将边界涡量流用叶型型线参数和型线上的气动参数表示,进而获得新的目标函数表达式。该目标函数表达式可直接对叶型型线参数求导,仅使用简单的定步长梯度算法即可对优化模型求解。结果表明:优化后的叶型型线上的边界涡量流的正峰值显著降低。     

PSE方法在后部加载叶栅绕流稳定性分析及转捩预测中的应用

采用数值求解与实验测量相结合的方法, 分析后加载叶栅绕流的稳定性并预测转捩.推导出正交曲线坐标系下三维扰动波的抛物化稳定性方程, 在风洞实验中测量叶栅叶片表面的静压分布, 并结合测量结果对设计叶栅的流动稳定性进行特征值分析.分析结果表明:(1)后加载叶型能有效推迟转捩, 提高叶栅的变冲角特性;(2)当冲角变化时叶顶流动稳定性的变化最显著, 负冲角流动稳定, 正冲角流动最不稳定.      

一种设计亚音速叶型的解析方法

本文在近似于势、流函数的正交面上采用解析方法进行亚音速叶型设计。将改写后的无旋条件与连续方程联立求解,可在计算面上将叶型设计问题用标准拉普拉斯方程来表示。这一方法的优点是在流场计算中用解析解取代数值迭代计算过程,求解简便,节省计算机时,并为从理论上分析各种不同叶型的设计特性提供了便利条件。文中提供的算例表明,方法是有效和可行的。      

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

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

叶型反设计计算研究

以某型号发动机压气机叶型生产坐标为设计输入,寻找一种切实可行的叶型反设计方法,以便对该压气机进行反设计计算研究,评估其真实性能,从而为改进设计提供依据.具体反设计过程为:根据等高面叶型生产坐标、设计经验和分析判断给出压气机各级气动参数分布及其子午面布局,得出其S2流面计算结果,然后冉根据S2流面计算结果将生产坐标反设计到流线面上,通过流线上的叶型反设计得出一维、二维及三维计算分析程序所需的计算输入并进行计算分析,对比计算和试验结果,判断反设计方法的有效性.     

轴流压气机叶片准三维自动优化设计

为提高轴流压气机叶片气动设计效果,将叶片准三维设计方法与叶型自动优化技术相结合,构建轴流压气机叶片准三维自动优化设计平台。在叶片三维数值计算结果中提取总压损失系数径向分布,反馈到S_2流面计算程序,替代原有的基于试验数据的损失模型,重新进行通流设计,提高S_2流面流场计算精度;S_1流面叶型设计采用并行遗传算法与叶轮机CFD方法相结合的叶型自动优化方法,以数值寻优代替叶型气动外形的人工调整,减小对设计经验的依赖。运用设计平台,对一设计流量6.3kg/s,设计压比2.07的轴流压气机第一级转子进行气动设计,得到径向三维叶片,设计点流量、总压比均接近设计目标,等熵效率接近90%,失速裕度达15.5%,结果令人满意。采用三维优化方法进一步进行叶片积叠线弯、掠设计,在失速点到设计点的稳定工作范围内效率均得到了1%的提升。     

复杂外形飞行器多块结构贴体网格生成

通过求解椭圆型方程生成三维贴体坐标网格,采用分块技术生成多块结构网格体系。每个子区域转换成计算空间内的立体域,具有自身的曲线坐标系,从而所形成的整个网格系统不会过分扭曲和不规则,这将有利于对于复杂边界的真实飞行器绕流流场的求解。文中讨论了控制函数对网格生成的影响。     

三维环境下离心/斜流压气机二维叶型优化设计

考虑离心/斜流压气机转子叶片通道内流动的强三维性,提出在三维环境下进行二维叶型优化设计.通过对能量方程中黏性耗散项改进,解决了Denton黏性体积力方法模拟离心/斜流叶轮三维流场效率偏高的不足.将改进流场计算模块与并行遗传算法寻优模块相结合,构成离心/斜流压气机二维叶型优化设计软件.采用所研制的软件,分别对离心叶轮和斜流叶轮叶尖处叶型进行优化设计.优化叶片基本达到目标流量和压比,在整个工作范围内效率都提高明显:在设计点离心叶轮效率由0.938提高到0.947,斜流叶轮效率由0.899提高到0.918.      

周向平均方法在某风扇/增压级分析中的应用简

 通过对三维（3D）Navier-Stokes方程进行周向平均,得到了通流模型的控制方程,对其采用时间推进有限体积方法进行数值求解。为实现风扇/增压级在设计初期的快速性能评估,考察了周向平均方法在风扇/增压级分析中的准确性。分别利用NUMECA三维数值模拟软件和周向平均通流模型（CAM）对某高通流风扇/增压级进行了性能分析,从对比结果来看,周向平均通流模型在近设计点给出了与三维数值模拟十分接近的特性参数,最大误差不超过2.0%。在风扇转子中,由于周向平均通流模型能捕获通道激波,其物理本质与三维平均结果有所区别,因此径向参数分布与三维有所差异。而在亚声速流动下的增压级及外涵道各叶片排出口参数的径向分布与三维数值模拟结果都能很好地吻合。     

周向平均方法在某风扇/增压级分析中的应用

通过对三维（3D）Navier-Stokes方程进行周向平均,得到了通流模型的控制方程,对其采用时间推进有限体积方法进行数值求解。为实现风扇/增压级在设计初期的快速性能评估,考察了周向平均方法在风扇/增压级分析中的准确性。分别利用NUMECA三维数值模拟软件和周向平均通流模型（CAM）对某高通流风扇/增压级进行了性能分析,从对比结果来看,周向平均通流模型在近设计点给出了与三维数值模拟十分接近的特性参数,最大误差不超过2.0%。在风扇转子中,由于周向平均通流模型能捕获通道激波,其物理本质与三维平均结果有所区别,因此径向参数分布与三维有所差异。而在亚声速流动下的增压级及外涵道各叶片排出口参数的径向分布与三维数值模拟结果都能很好地吻合。     

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

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

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

CAD技术在磨叶片叶尖工具中的应用

以生产准备中磨叶片工具的计算为例，使用三维CAD系统的自由曲面建模技术，精确计算出叶片工具设计中所需叶片通道数据。通过与传统计算方法的对比，可以看出：在工艺准备中使用三维CAD技术，不但能准确快速的解决生产准备过程中的复杂计算，而且能有效减少设计的错误，缩短工艺准备周期。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。