单级跨音压气机整圈三维动静叶干涉的数值模拟

本文使用商业软件FLUENT,对某单级跨音压气机的整圈三维非定常动静叶干涉流场进行了数值模拟,计算并比较了拉开动静叶轴向间隙和正常轴向间隙情况下的不同流动状况。非定常整圈三维计算的结果模拟了清晰的动静叶干涉图像,并显示了静叶三个连续通道出口流量的周期性变化和滞后现象,表明可以利用成熟的CFD商业软件研究叶轮机械内部的复杂物理现象。     

动静圆柱坐标系下N-S方程的分区算法及子域干涉处理研究

耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程,统一了动子域和静子域的控制方程;在子域计算中,采用高分辨率的三阶ENN格式及LU-SGS隐式解法;在子域干涉处理中,研究并对比了动静子域干涉面的重合处理与重叠处理;最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气机叶栅内三维紊流流场的数值模拟,验证了该分区算法的有效性,计算结果与实验数据吻合较好。     

低速轴流压气机中前后静叶对动叶顶部区域流动的影响

数值研究了某二级低速轴流压气机在高负荷工况下,动叶前后动静干涉对动叶顶部区域非定常流动的影响。结果显示一次泄漏涡和二次泄漏涡是动叶顶部的主要流动特征。下游静叶干涉增加了动叶顶部间隙流的非定常波动强度和泄漏损失。上游静叶干涉减少了动叶顶部二次泄漏流强度。同无动静干涉相比,同时存在上下游动静干涉使动叶顶部区域总压损失减少9.1%,其收益主要来自上游静叶的非定常干涉。     

叶轮机内部转/静干涉作用的数值模拟

本文对叶轮机械内部转/静干涉问题的非定常数值模拟进行了研究.采用双重时间步方法和Dornain-Scaling方法求解二维N-S方程组,数值模拟了某涡轮叶栅内部的非定常流场.通过计算展示了在转/静干涉作用下的叶型尾迹,转子尾迹的传输和尾迹的切割等非定常现象.     

跨声压气机动静干涉效应的数值研究

采用数值方法对一轴流跨声压气机在设计点的非定常流场进行了模拟, 对级内动静干涉进行了深入分析, 研究了尾迹和位势作用等对动静叶表面气动负荷的影响.计算结果表明:在压气机中, 上游动叶的尾迹等对静叶通道内部流动, 叶片表面静压的波动, 以及边界层流动损失的发生、发展和输运产生明显的影响, 需要在设计中加以考虑.      

某小型风扇内进口总压畸变引起流动损失分析

采用数值模拟方法对方波形式的进口总压畸变对风扇流场的影响进行了全三维非定常数值模拟.详细分析了畸变造成的动静叶栅流道内部流动损失的形式和原因.分析表明, 动静叶流道内的损失在畸变和非畸变区域向下游疏运的方式不同.在动叶根部和中径区域流动损失的主要形式是疏运损失和掺混损失.畸变会使多种损失相叠加.在静叶中, 畸变会在中径以下区域造成大范围的漩涡, 造成回流, 从而使损失增加.      

阵风与叶栅干涉噪声的数值模拟

采用计算气动声学方法求解二维线化欧拉方程对阵风与叶栅的干涉噪声进行了数值模拟。为减少频散和耗散误差,采用七点四阶频散相关保持有限差分格式进行空间离散;时间推进采用优化的4/6低频散低耗散龙格库塔法;在各边界处均选择或建立了适当的边界条件。首先通过模拟阵风在自由空间中的传播验证了无反射边界条件的可靠性;然后对阵风与平板叶栅的干涉问题进行了模拟,数值结果与精确解符合得很好;最后给出了阵风与NACA-0012翼型叶栅干涉的数值模拟结果,显示了数值求解方法的可行性。      

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

前置圆柱列涡轮静叶栅流场特性的数值研究

数值模拟了前置圆柱列涡轮静叶栅的三维流动,详细分析了圆柱列动静状态下叶栅流向的气动特性变化,发现不同时刻上游尾迹扫过下游叶片排不同位置,是造成多级涡轮非定常特性的主要原因。通过圆柱列动静两种状态下叶片表面压力脉动、流场变化特性的比较,发现静叶表面的压力分布存在差异,静叶前缘受上游圆柱列影响较大;随着流动的继续至叶栅流道下游,差异逐渐变小,直至消散;而叶栅前缘的差异主要通过叶片自身具有的大转折角逐步被抹平。     

涡轮叶片凹槽状叶顶非定常流动传热特性的数值研究

为了研究涡轮叶片凹槽状叶顶的非定常流动传热特性,以GE-E3第一级静叶和动叶为研究对象,采用ANSYS-CFX数值求解三维Reynolds-Averaged Navier-Stokes(RANS)和标准k-ω紊流模型。数值预测的叶顶换热系数分布与实验数据吻合良好,从而验证了数值方法的可靠性。数值计算结果表明:静叶尾迹对动叶顶部的流动和换热特性影响显著。压力面侧前缘区域和吸力面中间位置的流场受动静干涉影响显著。叶顶表面的换热系数脉动主要出现在靠近前缘的凹槽底部表面和再附着线附近及吸力面侧肩壁。靠近动叶尾缘区域的换热系数脉动同时受动静干涉作用和下游流场的影响。定常计算得到的换热系数在前缘冲击区和分离线附近高于非定常时均值,在压力面侧肩壁附近小于非定常时均值。定常计算得到的平均换热系数要高出非定常计算结果3.5%。     

一种改进的流场计算加速收敛方法

探讨了二维非结构网格无粘流场计算的加速方法。在对残差平均和多重网格这两种主要加速方法进行分析后,提出了一种改进的兼具多重网格思想的残差平均方法。经对绕NACA0012翼型无粘流动计算的验证,收效明显。在几乎不增加计算量的情况下,比原有的残差平均方法可减少迭代步25%～30%.该思想可推广用于三维、有粘流动计算和结构化网格中。     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

跨音速压气机在非设计工况下动静叶相互作用的数值模拟

本文计算了某跨音速压气机的一级在非设计工况下动静叶干扰产生的非定常流动。计算以准三维的Navier-Stokes方程为模型,湍流模式采用Baldwin-Lomax代数模式。主控方程用高分辨率的MUSCL型Roe格式离散,离散方程用高效的LU-SSOR方法求解。本文提出了一种完全守恒的动、静叶交界面处理方法。周期性边界处理采用Rai的方案。根据计算结果,讨论了非设计工况下动静叶干扰的机理。      

可压混合层涡并放热效应的涡动力分析

 基于双元单步不可逆燃烧动力学简化模型 ,采用紧致差分算法直接求解二维可压 Navier-Stokes方程 ,研究了时间发展非预混扩散火焰平面反应混合层流动。拟序结构由数值线性稳定性分析方法设定的粘性可压最不稳定小扰动初场促生。随放热强度增大 ,得到了基频涡卷撕裂现象。涡动力分析表明 ,放热主要通过膨胀和斜压这两种强非线性过程抑制大尺度剪切涡的演化     

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

轴流压气机转静二维非定常薄层N-S数值模拟

采用非定常薄层Navier-Stokes紊流运动方程来描述一跨声轴流压气机中径处的二维转-静非定常流场，速度脉动的相关分析表明转子尾迹对静子的热流流场几乎没有什么影响，而粘性的影响取决于叶片排之间的轴向距离。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

球柱三维分离及涡流的数值模拟

对三维可压非定常Ｎ－Ｓ方程进行有限体积中心离散,并用Ｒｕｎｇｅ－Ｋｕｔａ显式格式进行数值求解。数值模拟了球柱分离及涡流的形态结构,给出了具有激波、边界层、大分离区、涡流及其间相互干扰的复杂流型     

直升机旋翼悬停流场的欧拉方程计算

描述了三维Ｅｕｌｅｒ方程数值模拟旋翼悬停流场的过程,欧拉方程的求解采用了风格中心有限体积法、五步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进格式。