平面COUETTE流动的湍流斑研究

湍流斑是近壁剪切流动中的重要现象,它们的生成和发展与流动转捩和湍流的形成密切相关.采用Navier-Stokes方程直接数值模拟和壁面脉冲模型研究平面Cou-ette流动中湍流斑的生成和演化.高精度高分辨率的三维耦合紧致差分格式被用于数值计算.探讨了平面Couette流动中湍流斑的重要特征,包括高频脉动的形成、雷诺应力的产生、扰动的急剧增长和湍流斑形状的变化,特别是流向涡的复杂演化过程.结果表明层流流动中的湍流斑具有湍流流动中一些基本特征.     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

三维管内激波/湍流附面层干扰流场的数值模拟

给出了可压缩性及曲率修正两方程湍流模型，用它对跨音速三维喷管管内激波／湍流附面层干扰流场进行了数值模拟；给出了不同纵向截面上的激波结构、等马赫线图和四周固壁上的摩擦力线谱。将计算所得的摩擦力线谱、激波结构与实验结果进行了比较。分析了横截面上的流动结构，计算与实验吻合较好，进一步证实可压缩性及曲率修正两方程湍流模型能够较好地模拟激波／湍流附面层干扰流场，满足工程设计和分析要求。     

从壁面脉冲到湍流相干结构演化的数值模拟

发展了高效求解基于Navier-Stokes方程的扰动方程高阶紧致差分方法，研究了壁面强脉冲到湍流相干结构的复杂演化过程。该方法采用高精度、高分辨率的三维耦合中心差分格式，以及和高阶时间分裂法相结合因而更为有效，比一般谱方法更适用于较为广泛的流动区域和边界条件的湍流相干结构研究。中模拟了槽道流动中壁面强脉冲诱导的湍流相干结构的产生、发展和演化过程，显示了湍流相干结构的基本特征和发展规律．计算结果表明，所用壁面强脉冲模型比共振三波模型更有效．     

非正交曲线坐标湍流旋流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴体坐标系统，用坐标变换的方法处理锥形燃烧室内线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。紊流模型采用工程上常用的ｋ－ε模型。还运用贴坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提出了理论依据。     

后台阶流动的大涡模拟

利用大涡流模拟,对后台阶中的流动进行了数值计算,并和实验数据进行了比较,吻合较好,表明了所采取的亚网格湍流模型的合理性,但是需要改进。     

微小矩形通道内流动与换热特性

采用试验和数值模拟方法研究了当量直径为0.633～1.079 mm的微小通道内的流动与换热特性.工质为 液体水,雷诺数Re为100~3 000.试验结果表明,该尺寸范围的通道内层流向紊流转变的临界雷诺数约为1700,层流表观阻力系数和努塞尔数的试验值与数值计算结果吻合良好.最后通过数值模拟对恒定泵功及恒定流 量条件下不同结构参数通道的换热进行了计算,给出了最佳换热结构.     

竖直楔形流道的单气泡上升过程的实验研究

通过实验对竖直楔形流道内单气泡上升过程的行为特性进行了研究.利用高速摄影设备和计算机图像处理技术准确地追踪了气泡上升过程中运动轨迹和形态的变化.通过对气泡上升过程中旋转角、瞬时速度、瞬时加速度等参数的分析,找到了楔形流道内气泡上升过程的运动规律,并且通过推导获得了气泡上升过程中拽力系数的计算公式.     

跨声速风洞扩散段数值模拟

运用时间相关的全隐ＡＦ方法,计算二维及轴对称亚声速高速扩散段内流场。流动控制方程使用全Ｎ－Ｓ方程。湍流模型选择了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,和两层κ－ε二方程模型。经过同国外实验结果的比较,证明本文计算结果合理。在此基础上,计算了中国气动力研究与发展中心２．４ｍ×２．４ｍ引射式跨声速风洞第二扩散段在不考虑喷流干扰及常规运行状态下的流场。     

二元内流通道内湍流附面层分离的数值与实验研究

本文应用激光多普勒测速仪测量了二维扩压器中不可压湍流附面层分离流动,得到了时均速度和表面压强分布,表面摩擦系数和时均速度分布经验关系式做了比较。预测中采用准同时粘性/无粘流匹配模型,并比较了两种形状因子H_1经验关系式所得预测结果。在最大反流以前,所预测的附面层位移厚度及形状因子与实验吻合。     

钝体燃烧器后能量逆转和紊流参数的数学模拟

本文用非定常数值法研究了钝体燃烧器后的紊流参数分布和能量逆转,并对Boussinesq的雷诺应力模式和双方程k-ε紊流模型提出了改进方案。这样不仅计算出了能量逆转区,而且紊流参数分布也与实验相符。     

飞机座舱内空气速度和温度分布的数值模拟

飞机座舱内空气的流场和温度场直接影响到乘员的舒适性。本文应用雷诺平均Navier-Stokes方程和Launder-Sharma低雷诺数模型计算了飞机座舱内的流场和温度场。文中对对流－扩散方程中的对流项分别采用低阶和高阶精度格式进行了离散，计算结果表明高阶精度格式更适合飞机座舱内空气分布的数值模拟。     

非正交曲线坐标湍流族流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴你坐标系统，用坐标变换的方法处理锥形燃烧室曲线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。亲流模型采用工程上常用的ｋ－ε模型。还运用贴体坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提供了理论依据。     

来流马赫数对座舱气动加热影响的数值模拟

采用数值手段研究了来流马赫数对飞机座舱附近区域气动加热的影响,计算的来流马赫数为0.8~2.0,基于翼展的雷诺数为107。通过分析飞机蒙皮表面的温度分布、表面摩擦因数、壁面法向速度梯度和热流量等物理量,给出了飞机头部座舱附近区域的气动热分布情况。计算结果表明,随着来流马赫数的增加,壁面法向速度梯度增大,飞机对称面上的平均温度、表面摩擦因数和热流量随之增大。此外,针对某型号飞机的计算结果表明,座舱附近区域上表面中心线上的平均温度与来流马赫数、自由来流温度之间呈一定的函数关系。     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型,采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式,对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟,探讨了内外喷流及喷管内的流动特征,并与实验数据进行了比较,计算结果与实验结果吻合较好,验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性,从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

建筑物风环境三维定常数值仿真

针对建筑物的风环境问题，采用计算流体动力学CFD软件对两幢建筑物进行了数值模拟，得出了建筑物的风场分布，以及风载的分布规律，结果表明，气流在来流方向第一幢建筑物的顶部、两侧的前沿产生分离，建筑物之间产生强烈的漩涡，顶部的局部负压最大。     

轴流压气机过失速行为的数值模拟

基于轴流压气机级特性,建立了带燃烧室的多级轴流压气机过失速的一维逐级分析模型,发展了分析轴流压气机过失速响应的动态滞后方法,确定了压气机过失速响应的纺一时间常数。并将该模拟技术应用到某压气机中,通过调节敛喷嘴面积和燃烧室燃油流量等措施来诱导压气机失速。压气机旋转失速后,压气机流量下降,增压能力降低,燃烧室温度升高,可以通过调节燃油流量等措施来使压缩系统恢复正常工作。     

槽道湍流相干结构的研究

采用高阶差分方法求解非定常不可压Navier-Stokes方程，进行槽道流相干结构的数值研究，与常用的谱方法相比，不需限制边界条件，可用于更一般的粘性流体的流动问题，文中所用时间分裂法和三维耦合高阶差分格式适用于包括邻近边界点在内的整个计算区域，克服了三维各自用高阶中心差分格式不适用邻近边界点的困难。通过算例模拟了槽道流动中单个相干结构的演化发展规律，结果与实验一致。     

旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。