非定常分离流动的Largrange方法数值模拟

本文对Ｌａｇｒａｎｇｅ坐标系下的边界层方程进行了数值模拟，研究了静止流场中的线涡诱导的非定常边窠层流动，数值结果表明，随着边界层接近分离，确实出现了有限时间奇性。流动在驻定点附近很窄的区域内堆积，并在垂直壁面方向迅速增长，形成一脉冲似的突跃，基于此现象讨论了非定常分离的性质和三种分离模式，采用Ｌａｇｒａｎｇｅ方法能很好的说明有限时间奇性，便于精确地给出分离的定义和判据，而且物理意义非常清楚。Ｌａｇ     

非定常动静叶干涉作用的二维数值模拟

采用涡方法研究叶轮机内动、静叶间的相互作用,它利用全三维的、通道时间平均的定常解为基础解,并由此给出非定常扰动场的初始解。为计算叶片对扰动场的响应过程,采用拉格朗日方法追踪扰动涡团的对流流动过程,用确定性涡方法来描述流体的粘性扩散过程。用该方法对某型压气机的一级（动、静叶栅距比为１．７４）的动、静叶相互作用,进行了二维有粘问题的数值模拟。数值结果表明,该方法具有较好的收敛性和收敛速度;转子出口的尾迹和各种旋涡等非定常因素在静叶通道中的切割、迁移过程,非定常性对静叶攻角的影响以及脉动压力的振幅衰减特性也在文中作了分析。     

二维凹槽超声速湍流流动数值模拟

本文采用时间、空间均为二阶精度的NND差分格式,以及Menterk ωSST二方程湍流模型,数值模拟了二维开式凹槽超声速粘性流动。用组分扩散方程计算了凹槽内粒子驻留时间。本文给出了M∞=3,长深比L/D=3,后壁倾斜角为90°和30°两种几何形状凹槽的计算结果,与相应的实验符合甚好。计算亦表明,该模型较高地估计了涡粘性,从而使流动更稳定。     

多级压气机非定常流动的数值模拟

进行了多级压气机动静叶片排干扰的三维非定常无粘流动数值模拟研究，采用高分辨率格式求解非定常Euler方程组。对某以级压气机的非定常无粘流场的计算结果表明，即使在不考虑粘性影响的情况下，压气机动静叶排干扰的非定常效应仍然是比较显著的，各叶片排所受到的非定常影响是不相同的，计算得到的叶片表面压力周期性波动的结果可为叶片强度计算提供较为准确的数据。     

管道系统激波运动的数值模拟

对激波在管道，活塞腔系统内运动的两种数值模型－－一维非定常气体运动模型的二维轴对称气体运动模型，进行了数值模拟，其计算曲线与相应的实验曲线符合较好，表明基于Ｇｏｄｕｎｏｖ格式的激波运动计算方法能成功地处理含有激波且边界条件复杂的管道系统内气体运动问题，因而能较精确地预测管道系统内详细的和变化过程，特别是预测激波在管道系统内的传播过程。     

二维粘性不可压缩流场的分块耦合求解方法

本文用分块耦合的方法，并采用S.Abdallah^[1,2]和M.L.Mansour^[3]的在非交错网络上隐式求解原始变量的Navier-Stokes方程的提法和离散格式，成功地解决了二维复杂流场的数值求解问题，获得了空腔流、T型叉管流、NACA0020翼型绕流和方块钝体绕流的稳态解，并进行了讨论。此方法还可推广到三维复杂流场的计算问题。     

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

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

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

低马赫数流动数值模拟方法的研究

在文献[1，2]提出的一种低马赫数流动数值模拟新方法的基础上，本文作者进一步优化了算法，确定了该方法适用的马赫数范围，讨论了参数σ对计算的影响。大量的定常／非定常、低雷诺数流动的计算结果表明这种方法具有计算精度高、收敛速度快以及计算花销小的特点。     

二维旋转湍流边界层内气固两相流动的数值模拟

本文用修正的ＤＯＲＯＤ－ＦＥＭ方法对不可压流体二维旋转湍流边界层进行了数值模拟，得到了与实验数据较符合的边界层形状因子和壁面摩擦系数的计算结果。采用单向耦合（Ｏｎｅ－ｗａｙ Ｃｏｕｐｌｉｎｇ）方法，计算了固粒在湍流边界层和无粘流中的运动轨迹，给出了壁面被固粒撞击后的冲蚀磨损分布，分析了近壁流场性态与壁面冲蚀磨损的关系。     

二维流场、热、结构一体化数值模拟

流场、热、结构问题,在传统的计算中,通常是先算流场、再算热响应和热结构问题,这显然是人为地分割了本身统一的、连续的过程,使得计算结果存在较多不足。本文给出一种把流场、热、结构耦合起来进行一体化数值模拟的方法。算例证明,本方法较好地模拟了动态的物理过程。     

进口热斑对涡轮级影响的非定常数值模拟

通过求解二维N-S方程对涡轮级进口有温度畸变时的流场和温度场进行了非定常数值模拟。由于叶排间的相互运动,使叶片表面压力和叶栅通道内的流场呈周期性地变化;当涡轮进口有热斑时,热斑中的气流总温要高于周围气流温度,在通过涡轮级时,热斑中的高温气流会向转子压力面上迁移,导致转子压力面上产生热点。热斑的存在对涡轮级中的压力分布影响很小,但是严重影响温度场的分布。同时进口热斑相对导叶的位置对叶片表面温度的分布也有很大影响。      

高压涡轮全环非定常流动数值模拟

为了进一步准确地掌握涡轮内部复杂流动特征,采用全环非定常数值模拟方法研究了某型高压涡轮内部非定常流动,重点分析了上游导叶尾迹和激波在动叶流动中的非定常作用机制,探讨了引起动叶通道内及叶片表面压力分布周期性变化的非定常扰动。数值结果表明:与尾迹作用相比,由导叶激波引起的非定常效应更为明显,尾迹引起的非定常扰动在动叶0.5倍轴向弦长位置下游较为突出。在动叶流场中及叶片表面上的静压变化显示,在动叶通道内还存在低频扰动,对局部流场的非定常扰动明显。最后指出,在非定常计算时采用区域缩放法对流场中非定常扰动的预测存在一定误差。     

模拟二维周期性非定常流动的谐波分析方法

介绍了一种从频域出发模拟二维周期性非定常翼型绕流的数值方法.考虑到很多工程中遇到的非定常流动问题都是周期性的,因此非定常流场中的守恒变量可以展开为傅立叶级数,其系数随空间位置的不同而改变.将这样的谐波假设引入欧拉方程,经过处理的方程不显式地含有时间导数项,实际求解的方程将具有定常方程的形式.在使用多重网格技术后,计算可以很快达到收敛.对NACA64A010翼型典型的周期性非定常绕流进行的数值模拟表明,这种方法在能够给出较为准确的计算结果的同时,与传统的双时间方法相比,可以大幅提高计算效率.     

绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

进口气流角对三维振荡涡轮叶片非定常流动影响的数值模拟研究

进一步发展了非定常雷诺平均N—S方程求解程序,并基于影响系数法,在三种不同振型下就进口气流角对三维涡轮振荡叶片绕流的影响进行了数值模拟研究。数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;进口气流角对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,且在不同模态下,进口气流角对振荡叶片非定常流动的影响规律表现不同。     

尾迹对压气机转子性能影响的非定常数值模拟

采用三维粘性非定常数值模拟方法,研究了上游尾迹对轴流压气机转子性能及其尖部非定常流动的影响.结果表明,在一定情况下,上游静子尾迹与转子内部流动的非定常相互作用,有可能改善近失速点的气动性能,如转子压比和效率升高,工作范围增大.其原因主要为:上游静子尾迹使转子尖区一次泄漏涡强度减弱,减少了二次泄漏涡强度或抑制了二次泄漏涡的产生,最终导致尖区损失减少;此外,尾迹使尖区激波位置后移,改变了尖部弦向的负荷分布,最终导致压气机稳定工作范围增大.     

用CSCM格式模拟二维收缩管道内的流动

 采用Lombard等人提出的CSCM(Conservative Supra-Characteristics Method)差分方法数值模拟二维收缩管道内波系干扰及激波边界层干扰等问题。所得结果同已往计算结果进行了比较。通过压力场和速度场的计算结果可看到Ma数对流场的影响。在Ma数较小时激波角较大,可看到激波经两次相交一次反射逐渐变弱的情况。Re数改变对流场也有明显影响,且给出由于激波干扰引起边界层分离的计算结果。最后还给出收缩管道内激波由正激波逐步形成斜激波的非定常过程。     

双级风扇三维粘性非定常流动的数值解

针对一台双级跨音速轴流风扇动静叶片排相互干扰的三维非定常流动,进行了N-S方程的数值求解。首先计算风扇定常流场,以定常流收敛结果为起始场计算非定常流动。计算结果表明,非定常流动对各叶片排进口气流角有较大的影响,其中静叶的变化大于动叶,与定常流结果相比,非定常的瞬态攻角增大的量大于减少的量,最多时增大近10°左右。静子叶片受到的非定常气动力变化幅度也大于转子叶片,叶片排所处的轴向位置不同,非定常流动的影响会有较大的差别。