有限长度的过渡领域微槽道流动（一分析预测和试验的比较）

在连续介质和滑流区有限长度槽道内的稳定气流，通过整体质量守恒规律得到了沿槽道轴向的压力分布和质量流率的分析表达式。在过渡领域区，也提供了槽道流动压力分布的严格的动理学理论解。在Poiseuille流动的质量流率动理学理论解给定的情况下，压力分布的控制方程被证实是被前一作者从轴承润滑问题移植到解微槽道问题的退化Reynolds方程。同时得到了质量流率分析表达式，它们和现有试验数据的比较相符很好。     

混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

疏水微槽道内层流减阻的实验研究

研究了光滑和带有横向凹槽结构的疏水微槽道内层流的流动特性和表面滑移效应。在硅片上加工了矩形截面微槽道,利用十八烷基三氯硅烷（octadecyltrichlomsilane OTS）在槽道内壁形成疏水薄膜。实验结果表明在光滑疏水微槽道内的水流压降比无滑移条件下的理论值减少8％。对于侧壁带有凹槽结构的疏水微槽道,流动阻力可以降低10％～30％。笔者采用micro—PIV测量得到的壁面表观滑移速度约为槽道中心速度的8％,滑移长度约为2μm。实验结果与滑移壁面条件下三维槽道内层流的解析解吻合,同时得到了带有凹槽结构的疏水微槽道内的流速分布。     

固体发动机喷管内跨音速旋转流动研究

对有旋转情况下的喷管内跨音速流动进行了理论分析, 建立了流动模型, 导出了壅塞截面上的控制方程, 并用数值方法对方程求解, 较细致地刻画了有旋转情况下壅塞截面上参数分布。初步研究表明旋转能够有效地减小喷管质量流率。      

奇异积分方程方法在接触压力分析中的应用

总结回顾了奇异积分方程的数值解法,对于第二类奇异积分方程使用分段连续函数法进行了求解;对于两个弹性体接触问题,通过接触体之间的滑移函数和间隙函数,建立求解接触压力的奇异积分方程.分别针对圆柱体与弹性半空间体、抛物线型压头与弹性半空间体、榫头与榫槽3类接触问题,确立奇异积分方程的具体表达式,而后使用分段连续函数方法进行求解,获得接触面上的接触压力.最后将计算所得的接触压力分别与理论解和有限元解进行了对比.对于圆柱体与弹性半空间体接触问题,奇异积分方程法的最大接触压力与理论解和有限元解的相对误差分别为0.3%和0.5%;对于榫头与榫槽接触问题,奇异积分方程方法计算所得的最大接触压力与有限元解的相对误差为1.8%,验证了奇异积分方程方法的有效性.      

不可压湍流反应流直接模拟和RANS反应模型检验

用谱方法对三维不可压槽道湍流反应流动进行了直接数值模拟,得到了温度和质量分数耦合的瞬态数据库.结果显示温度与质量分数脉动在近壁区都有条带结构.推导了雷诺平均方法中湍流二阶矩反应模型中所求关联量的精确输运方程,发现耗散项是化学反应率系数-质量分数关联量封闭的关键,耗散项需要考虑化学反应影响.基于数据库的统计结果,对关联量模型方程中各项进行了先验研究,发现产生项和耗散项的贡献最大,扩散项和反应项的贡献较小,化学反应对各项大小和分布形状有明显影响.在算例中,直接模拟统计得到的化学反应率系数-质量分数关联和用代数二阶矩模型的模拟值很接近,说明ASOM模型具有一定的合理性.      

利用上游槽道对角区马蹄涡的控制

在平板上放置圆柱形成角区流动,利用布置在圆柱上游平板上的二维和三维槽道来控制或削弱角区马蹄涡,采用风洞试验和数值模拟开展研究。结果表明,二维和三维槽道均能推迟边界层的分离,使圆柱根部马蹄涡的强度减弱、尺度减小;同时槽道上游压力和逆压梯度均有所下降,槽道下游压力显著升高而逆压梯度总体降低。二维槽道对马蹄涡强度的削弱为61.15%~66.51%,而三维槽道对其削弱为66.65%~80.93%。讨论了三维槽道参数(包括槽道宽度、深度以及其中心线与圆柱中心距离)对控制效果的影响。槽道与圆柱的距离在对马蹄涡的控制中起主导作用。槽道控制的机理是,由于槽道的抽吸效应使得其上游靠近壁面的边界层中涡量较高的流体被卷吸入槽道形成槽道涡,槽道涡由三维槽道输运到下游。同时,随着槽道与圆柱的距离减小,更多的边界层流体流入槽道内。正是上述"槽道效应"使得槽道下游的逆压梯度降低,马蹄涡强度减弱,分离区范围减小。     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

槽道流动中湍斑的数值模拟

采用Navier-Stokes方程直接数值模拟槽道流动中湍斑的形成和发展.建立高精度、高分辨率的三维耦合中心差分格式,并和时间分裂法、高阶迎风差分格式一起组成高效的高阶紧致差分方法;克服一般高阶中心差分格式不适用于边界邻域的困难,可用于包括边界邻域在内的整个计算区域;在边界上满足Navier-Stokes方程的无反射和压力条件,保证了格式的精度和一致性.通过模拟槽道流动中从壁面脉冲到湍斑的演化过程,显示出湍斑内雷诺应力和速度、压力脉动等湍流的基本特征,表明了湍斑的生成与壁面脉冲的强度和Re数大小直接相关.     

考虑分子转动自由度的离散统一气体动理学格式

离散统一气体动理学格式（DUGKS）是一种适用于全流域模拟的有限体积方法。之前的研究考虑了分子平动自由度，验证了DUGKS在多尺度问题中的准确性及稳定性。本文基于Rykov模型方程构造了离散统一气体动理学格式，并采用Landau-Teller-Jeans转动能量松弛模型，可用于双原子气体从连续流动到稀薄流动的多尺度问题计算。测试了激波结构、超声速平板绕流以及超声速圆柱绕流等非平衡流动问题，计算结果显示出双原子气体分子中存在平动自由度与转动自由度对应的能量交换过程，并与统一气体动理学格式（UGKS）、直接蒙特卡罗（DSMC）方法的解以及实验值吻合较好。     

长微管道气体流动的近似解析解

对长微管道中压力驱动的气体流动，考虑到稀薄气体效应，给出了可压缩气体流动简化控制议程的近似解析解，得到了流动的速度分布、密度分布以及近似质量流量公式，分析了影响流量的因素。     

稀薄流到连续流的气体运动论统一数值算法初步研究

从非线性模型Ｂｏｌｔｚｍａｎｎ方程出发,引入简化速度分布函数、使用离散速度坐标法对速度空间进行离散、降维,去掉分布函数对速度分量的连续依赖性;采用时间分裂法,将简化速度分布函数松驰变化方程分解为源项碰撞变化方程、对流运动方程,进行耦合计算,应用ＮＮＤ耗散差分方法直接模拟气体分子速度分布函数;发展离散速度数值积分法,通过宏观取矩获取物理空间各点的流动参数,从而建立一套能有效模拟各流域气动问题的简化的     

直接供气预旋转静系流动特性的实验

用实验方法研究直接供气转静系预旋腔内的流动特性,得到了静盘表面的静压分布、预旋孔的排气系数以及盘缘出口流量的分配规律.结果表明:流量改变对静盘压力分布形状的影响较大,而转速改变的影响则较小;在几何条件不变以及预旋进气和中心进气流量相等的前提下,预旋孔排气系数主要由预旋进口流量系数决定;对于该实验模型和研究工况来说,盘缘出口流量系数比在0.7以上.      

基于基尔霍夫面的四维声学频域公式

流动诱发噪声问题是实际工程领域极为普遍的问题之一,经典的声比拟模型仅以声学压力为参考来评估声场特征分布还远远不够。从声压和声学速度矢量为变量的四维线性波动方程出发,选择包围非线性声源的基尔霍夫面为积分面,并结合对流格林函数,给出均匀运动介质的四维声学频域积分公式;针对静止、旋转单极子源和偶极子源开展数值预测研究。结果表明：本文获得的声压和声学速度分布与解析解吻合,均匀来流情形下静止点源的声场分布表现出典型的对流效应;受均匀来流、点源的自激频率、谐波阶次和旋转频率的共同影响,旋转点源的声场分布则表现出明显的多普勒效应和对流效应。本文对均匀流矢量气动噪声开展的精细化研究能够为声能量的评估和传输路径预测提供技术支撑,进一步为降噪研究提供理论依据。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

涡轮叶片尾缘扰流柱通道流动换热计算

为了深入研究涡轮叶片尾缘扰流柱区域流动的工程计算方法,用开发的基于流体网络计算的一维工程算法,对工作叶片尾缘扰流柱通道进行了计算,并与实验结果进行了对比。计算结果表明:工程计算得出的弦向出流量分布、压力分布从规律和数值上与实验结果都吻合较好,工程算法能够比较准确地计算出扰流柱通道的流动;扰流柱区的气体流动方向基本为弦向,但也存在着一定的径向流动,所以有必要在流体网络计算中设置表示径向流动的流路。     

迷宫密封的理论与实验研究

 根据相似原理。通过实验回归分析法,得出压力损失系数ζ与齿腔几何特性和流动特性的关系式。然后应用能量定律、连续方程、状态方程和Bernoulli方程与其联立,求解后可获得迷宫密封的全部特性。并可归纳成流量系数α、穿越系数β和迷宫函数ψ对全压比λ的特性曲线。应用这类曲线,可方便地计算出迷宫密封的泄漏率。     

基于优化算法的压气机叶片气动设计

将数值最优化方法与叶轮机械正问题流场计算相结合,实现对风扇/压气机叶片积叠线、子午面流道和沿径向弦长分布组合的优化设计,数值最优化采用基于局域网的并行遗传算法.运用该软件对NASArotor67风扇转子进行重新设计,设计目标为设计点压比和流量不变、效率最高以及堵点流量不变.与原风扇转子比较,稳定工作流量范围略增大并且在整个工作流量范围内优化转子压比近于不变、效率提高约0.5个百分点.      

径向冷却通道截面特征长度对气冷稳定器性能影响

通过研究径向冷却通道截面特征长度的改变对气冷稳定器的冷态流场、冷却性能及流动损失的影响，为其进一步的优化设计提供参考依据。采用经过实验验证的CFD方法，对某气冷稳定器进行了数值模拟研究，在所研究的特征长度范围内，结果表明:特征长度的增加，减少了径向冷却通道下部的进气阻力，使得风兜引气率增加了667%；下部出气孔出流由于受到回流涡的影响，出流量明显低于上部出气孔，特征长度的改变对出流量分布的影响不大；特征长度的增加会加强径向稳定器近环形稳定器处的外壁面冷却效果，但是会减弱径向稳定器近中心锥处的外壁面冷却效果；特征长度的增加，使得出口总压恢复系数下降了0.086%，沿程总压恢复系数曲线向下移动。      

塞式喷管底部压强模型

为了工程计算的需要，首先分析了已有的底部压强计算模型，找出了影响底部压强的主要因素，即塞锥出口截面的压强、马赫数，塞锥出口倾角和底部二次流流率；然后根据数值计算得到的底部压强和冷流实验测出的底部压强，分别拟合底部压强与塞锥出口压强和马赫数的关系式，底部压强与塞锥倾角的关系式，底部压强与底部二次流流率的关系式，最后把这些关系式几何平均得到了所有这些因素与底部压强的关系式。与实验数据比较，该拟合公式不但很好地反映了底部压强的变化规律，而且能够满足工程估算的需要。