旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

受限空间内射流影响下主流内气固两相流动实验研究

为了研究壁面射流在受限空间内对不同粒径颗粒群分布的影响及其作用机理,搭建了横向矩形通道内的射流气固两相流动实验台。通过调整均匀流动段内的隔栅网数量和位置,在无射流条件下使不同粒径颗粒在高度方向上均匀充满整个横向通道。在主流雷诺数7.9×104,射流雷诺数1.8×104的条件下,采用PIV系统对流动通道中的气流场和颗粒分布进行测量。由瞬态涡结构图片发现,带主流受限空间内的射流涡结构与自由射流涡结构具有明显差异,存在不对称的连续涡结构发生、脱落与破裂过程。通过观察射流对不同尺度颗粒的影响,详细分析了射流对颗粒分布发生影响的两大机理,并由此预见了壁面射流控制手段在工程上的应用性。     

基于全隐式紧耦合算法的气动弹性数值仿真

为进行气动弹性问题的计算,提出了一种全隐式紧耦合算法,在子迭代过程中分别采用LU\|SGS隐格式和隐式线性多步法交替求解气动和结构方程,以获得物理时间域的高精度解。一种径向基函数和超限插值结合的方法被用来进行气动网格的快速变形。运用该算法,进行了Isogai wing和AGARD 445.6 wing的颤振分析,颤振边界的计算结果与文献值和实验值较符合,表明该全隐式紧耦合算法能够有效地计算气动弹性问题。
     

二维叶栅气固耦合颤振分析

引入叶片结构动力学方程,通过数值求解气动/结构耦合方程体系,研究了气固耦合条件下叶栅内的二维非定常粘性流动,并初步进行了叶栅的颤振分析。针对两种叶型分别计算了不同折合频率下的流场,通过比较振动叶栅相关参数随时间变化的曲线,表明折合频率是决定叶片颤振与否的重要因素。      

控制体有限元方法中的一种高阶插值格式

针对控制体有限元方法对流动方向不平行于坐标轴方向的回流问题数值模拟采用各种迎风格式都容易产生假扩散的问题,本文给出了一种由插值单元及其迎风相邻单元所确定的高阶插值格式。通过对Re=1000时三种不同形状的二维四边形空穴流模拟,表明该格式与局部斜迎风格式的适当组合可有效地减小控制体有限元方法对回流问题数值模拟的假扩散。      

区域分裂法在湍流流场数值解中的研究

应用区域分裂算法(简称DDM)对湍流流场进行数值研究。对湍流流动控制方程建立一种通用的DDM-1计算格式S_{I,ω,给出重叠区域内变量处理方法。DDM-1的SI,ω计算格式的适用性由突扩圆管分离流流场数值计算所验证,数值计算结果与实验数据比较表明:采用SI,ω计算格式所预测的再附点位置及平均轴向速度、湍动能强度比在整个计算区域内采用单一模型的方法能够更好地与实验值吻合。}