串列钝体绕流的数值计算

用二阶投影法求解二维不可压粘性流体的N -S方程 ,计算了高雷诺数Re=1× 10 5下串列圆柱的非定常绕流 ,得到的阻力系数及斯特劳哈尔数与试验结果吻合良好 ,结合工程实际 ,用同一程序计算预测了钢管混凝土拱桥哑铃型拱肋的阻力系数和涡脱频率     

跨声速旋成体弹丸气动力的数值计算

用三阶MUSCL TVD格式解三维可压缩流动N-S方程组,湍流模式为Sparlart-Allmaras模型,数值模拟了来流马赫数Ma∞=0.96、攻角α=0°和Ma∞=0.91、攻角α=2°条件下的旋成体弹丸跨声速流场,获得的弹体表面压力系数分布数值结果与实验数据吻合良好。同时还给出了流场的马赫数等值线、弹丸底部流线分布图谱,并对计算的弹底涡流结构作了扼要的分析。     

三角形沟槽旋成体表面减阻性能的数值模拟

基于仿生微小非光滑表面具有减黏降阻特性的基本思想,在高速转动旋成体表面布置不同深度和间距的三角形沟槽.采用RNGκ-ε模型对其三维流场进行模拟,分别计算表面光滑旋成体与表面具有三角形沟槽的旋成体壁面阻力系数,对比两者壁面剪应力大小可知,将三角形沟槽布置于高速旋转的旋成体表面,可降低旋成体在高速转动时壁面的空气阻力,从而降低动力消耗,并且沟槽深度和间距均对旋成体壁面阻力产生不同影响.与光滑旋成体相比,三角形沟槽旋成体最大减阻率为12.060%.     

二维机翼振动非定常流场的数值模拟

准确数值模拟机翼动态失速时非定常流场主要有湍流模型和数值算法两方面的难点。在研究湍流模型的基础上 ,采用了 q-ω低 Reynolds数双方程湍流模型 ,并采用了一种高精度、高分辨率的数值计算方法 (L U- SGS- GE隐式算法和四阶 MUSCL TVD格式 )数值模拟了振动机翼周围的附着流动、轻度失速和深度失速非定常流场。结果表明 :采用本方法计算附着流动和轻微失速状态流场时 ,其升力系数迟滞曲线和阻力系数迟滞曲线与实验结果吻合得令人满意。深度失速的计算精度比前人工作有明显改善     

弹丸跨声速粘性绕流数值计算

用数值模拟方法研究弹丸跨声速绕流的流动特性,数值模拟的出发方程为平均雷诺数的Ｎ－Ｓ方程,湍流模型有用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型,数值格式为四阶的ＭＵＳＣＬ ＴＶＤ格式。     

高雷诺数下双圆柱绕流诱发振动的数值模拟

基于表面涡方法研究了不同间隙比及归一化阻尼参数下,串列和并列双圆柱的流体诱导振动问题,并计算了圆柱的振动响应、流体力、振动频率等.结果表明:圆柱为刚性时的计算模拟结果与相关文献的实验结果吻合良好;在文中的研究条件下,串列双弹性圆柱具有屏蔽流的特征,涡脱落频率比刚性圆柱的大,其值主要取决于归一化阻尼而非间隙比;流体从并列双弹性圆柱分离的涡是对称的,与并列双刚性圆柱相比其平均力系数略小,但由于流体诱导振动造成的脉动力系数较大,故流体诱导振动对并列圆柱的动态响应有重要的影响.研究结果为亚临界雷诺数范围内热交换器管束的复杂流体诱导振动等问题提供了一种可能的计算方法.     

齿轮传动非稳态热弹流润滑数值计算模型的研究

根据弹流润滑理论的基本原理，建立了齿轮传动非稳态热弹流润滑数值计算模型。该模型不仅考虑了齿轮传动中系统振动引起的动载荷、卷吸速度、曲率半径随时间和坐标的变化，以及流体的可压缩性等非稳态效应．还考虑了齿轮传动中的热效应和流体的非牛顿性。综合考虑这些因素可准确模拟齿轮传动的实际工况，更有效地指导实际生产。     

旋成体滚转过程中Bistable流场的数值模拟与实验

通过风洞测压与数值模拟研究了旋成体滚转过程。结果表明,侧向力的滚转角效应与模型加工精度有关;位于模型头部尖端一定量级的微扰动不仅能触发非对称,而且促使流场呈现Bistable状态。然而法向力无明显的滚转角效应。研究流场表明Bistable的左右两类涡系在结构上无明显差异,也并非完全的镜像状态。通过控制微扰动的周向位置可以获得期望的非对称多重涡流场。涡系诱导作用增强,对侧向力和法向力都有增值。     

跨音速压气机间隙流与处理机匣相互作用分析

基于对压气机转子顶部间隙泄漏流的深刻认识，针对某跨音速轴流压气机转子，设计了一种新型的处理机匣结构，并对带处理机匣的压气机转子内部流动进行了全三维非定常数值模拟，数值计算所获得的总性能（实壁）与试验结果符合较好．该种新型处理机匣结构的引入能在不降低压气机设计点效率的前提下有效地提高压气机的失速裕度．对处理机匣与顶部间隙泄漏流之间的相互作用机制进行了详细分析．结果表明：处理机匣结构能抑制间隙泄漏涡破裂现象的发生，并将间隙泄漏涡破裂后导致的阻塞区抽吸进入处理槽，从而有效地提高了跨音速压气机的失速裕度．     

高超声速双锥体绕流的数值计算与流场分析

 双锥几何外型的高超声速绕流问题涉及激波-激波相互干扰以及激波-边界层相互作用等复杂流动现象,是计算流体力学领域的热点问题,常被作为此类复杂流动的典型算例。本文从守恒型Navier-Stokes方程组出发,考虑了非平衡态气体的振动激发过程,采用高分辨率TVD有限差分格式研究了5种来流马赫数（15.56,10.34,9.59,8.20和8.06）,2种气体环境（纯N2和纯O2）,以及2种双锥外型（尖头体和钝头体）的气动热力学特性,较详细地给出了上述各种情况下壁面热流密度分布以及气动力系数等重要参数,所得结果与国外相关实验数据吻合较好,初步显示了本文所编程序的有效性。对差分格式精度和典型限制器（minmod、van Leer限制器）进行比较,发现数值计算所得分离区大小受格式耗散性影响很大,其中使用van Leer限制器所得到的分离区大小与实验结果吻合得较好。     

摇荡物体周围粘性流的数值模拟

提出了一种新的求解非定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值方法，在每个时间步长上采用人工可压性技术，通过对动量方程取散度，从而获得离散化的Ｐｏｉｓｓｏｎ方程，当关于虚拟时间的迭代过程达到收敛时，便得到了无散速率场，采用该方法模拟了一个纵摇椭圆柱体周围的非定常不可压性流动。     

任意旋成面跨音速叶栅流动的一种高效率解法

提出了任意旋成面叶栅跨音速绕流的一种高效率解法。导出了以Ｖｏｎ Ｍｉｓｓｅｓ 坐标为自变量的任意旋成面叶栅中流动的流线控制方程；简单而有效地处理了速度的双 值问题。用数值方法研究了流面半径和流片厚度随轴向变化对于流动的影响。     

二维管道内交变流动的格子-Boltzmann方法模拟研究

采用格子-Boltzmann方法（LBM）对二维管道内不可压缩交变流动——Womersley流的流动与换热进行了数值模拟研究．结果表明：管道内压力梯度呈周期性变化，在周期较小时，管道的绝大部分区域速度分布相当平坦，在靠近壁面附近出现了速度的峰值，且靠近壁面处速度变化剧烈；在周期较大时，速度分布近似于抛物线，管道中心速度总大于壁面附近速度，且不同周期下的模拟结果与解析解都吻合得很好．另外，当速度交变时，温度表现为波动变化．速度与温度的相位差、温度的分布特性、温度的波动特性等都与压力梯度交变周期及幅值有关．研究结果表明，LBM可以用于交变流动的模拟．鉴于LBM的众多优点，该方法有望成为求解复杂流动与换热的一种有效的数值模拟手段．     

跨音速下旋成体干扰流场计算

应用粘性-无粘性干扰方法计算了亚跨音速下旋成体的轴对称绕流流场,并着重研究了船尾处的后体流场.将流场划分为无粘的位流区和粘性影响集中的边界层区,利用有限差分法和积分关系式法分别求解之后,以边界层位移厚度为粘性匹配条件进行耦合计算,求出弹身表面的压力分布.该方法可直接用于计算导弹燃气喷流产生的“阻塞”效应对船尾-后体流场的影响。     

单方柱绕流的直接数值模拟研究

为了揭示柱体绕流的湍流流动机理,采用直接数值模拟方法对雷诺数为1200的单方柱绕流工况进行研究。首先通过与文献中斯特劳哈尔数、平均流速和表面压强系数等统计量进行对比,验证了数值方法的可靠性。其次采用本征正交分解方法系统地提取流场中的相干结构,结果表明：对于速度场的模态分解,第一阶模态代表平均速度场的特征,第二、三阶模态提取的是流场中的低频大尺度旋涡特征,第四、第五阶模态提取的是流场中的高频小尺度旋涡特征。最后基于联合概率密度函数分析速度梯度张量第二、第三不变量,发现方柱下游大致可分为两个流动阶段：发展阶段,流场以涡流层结构和耗散作用为主,涡流管结构逐渐生成;成熟阶段,流场中湍流结构伴随着高涡量拟能和高能量耗散率。     

Numerical Simulation of Rotor-stator Interactions in a Transonic Compressor Stage

ＩｎｔｒｏｄｕｃｔｉｏｎＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｒｏｔｏｒｓｔａｔｏｒｉｎｔｅｒａｃｔｉｏｎｓｈａｓｒｅｃｅｉｖｅｄｗｉｄｅｓｐｒｅａｄａｔｔｅｎｔｉｏｎｉｎｔｈｅｌａｓｔｆｅｗｙｅａｒｓｂｅｃａｕｓｅｏｆｔｈｅｉｎｃｒｅａｓｉｎｇｉｎｔｅｒｅｓｔｉｎｕｎｄｅｒｓｔａｎｄｉｎｇｔｈｅｕｎｓｔｅａｄｙｎａｔｕｒｅｏｆｆｌｏｗｓｉｎｔｕｒｂｏｍａｃｈｉｎｅｒｙ.Ｔｈｅｒｅａｌｆｌｏｗｓｉｎｓｉｄｅａｔｕｒｂｉｎｅｏｒａｃｏｍｐｒｅｓｓｏｒａｒｅｕｎｓｔｅａｄｙａｎｄａｒｅｃｈａｒａ…     

部分流泵非定常流动分析

在设计工况下采用滑移网格技术对部分流泵进行了整机非定常流动数值计算,分析了叶轮、蜗壳内典型的非定常流动规律.计算表明:部分流泵内流的非定常特性明显,蜗壳喉部尤为剧烈;瞬时扬程的变化是由于叶轮与蜗壳喉部强烈的非定常流动特性所决定的,只有采用非定常计算方法才能反映其内流的实质.该计算为进一步研究部分流泵的内流现象、提高效率、减少水力损失提供了一定的理论依据.     

不同头型回转体低速倾斜入水过程流场特性数值模拟

基于有限体积法和VOF多相流模型求解气、水两相流动的RANS方程,结合动网格技术,对回转体低速倾斜入水过程进行数值模拟研究.通过将数值计算结果与试验数据进行比较,验证了数值计算方法的有效性.基于该方法对不同头型回转体低速倾斜入水过程进行分析,得到不同头型条件下回转体入水空泡形态发展规律、回转体运动特性及流体动力特性变化规律.研究结果表明:同一入水深度,入水空泡直径和流场中最小压力值随着回转体头部锥角的增加而增大;不同头型回转体,锥头母线同一位置处压力随着头部锥角的增加而增大;回转体阻力系数与其头部锥角大小直接相关,锥角较大时,阻力系数也较大,速度衰减也较快.