突扩后突缩的二维流场计算

本文模拟固体火箭发动机药柱后端流场,建立了简化的计算模型,采用压力耦合方程半隐式解法,对突扩后突缩的二维粘性流场进行了计算分析.在计算过程中,对网格划分及控制体的选择、复杂边界条件的处理、差分格式的选取、方程的离散化及方程的解法等做了分析研究.计算结果表明:在突扩台阶后得到的二次涡与Finash,F.和Ghaddar,N.K.等人的结论吻合,进一步证实了Pan,F.和Acriuos,A.的嵌套涡族理论.     

突扩间隙比对单通道扩压器内流动力特性影响的试验研究

为了研究突扩间隙比与扩压器内流动力特性的影响关系,采用水流模拟流动显示及PIV(Particle Image Velocimetry)流场测速试验方法,对不同突扩间隙比的内流瞬态流场、时均流场、雷诺切应力等关键流场信息所表征的扩压器内流动力特性进行了研究。结果表明:突扩间隙比为1.6～1.7时,静压恢复系数与总压损失系数比值达到最大,此时扩压性能最优;扩压器突扩间隙比为1.64时,主流从前置扩压器流出时由于强剪切作用卷起扩散涡,为维持突扩区域空间内流动稳定性,进而诱导出消耗涡;总压损失变化是突扩区消耗涡、回流扰动及火焰筒前缘背压反流共同作用的结果,雷诺切应力与总压损失系数成正比;扩压器出口速度分布对于突扩间隙比变化的不敏感性体现出较好的内流流动稳定性。     

跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

平面突扩管道中的湍流流动及预混燃烧计算

应用随机涡法(Random Vortex Method)对平面后台阶突扩管道中的湍流流动及预混燃烧过程进行了数值计算。与平均流场所显示的气流结构不同,瞬时流场表现出大尺度的多涡结构。由于这些大尺度旋涡之间的相互作用,回流区长度在平均值附近发生拟周期性变化。预混燃烧计算结果表明,大尺度旋涡对火焰面在流场中的分布有重要影响。     

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

侧面突扩燃烧室冷态流场可视化研究

燃烧室流场中旋涡的不是稳定是造成整体式压式发动机侧面扩燃烧室燃烧振荡的重要原因。在透明矩形侧面突扩燃烧室模型上进行了水流模拟显示实验。显示出了燃烧室流场中的振荡涡系和稳定涡系。实验发现：燃烧室头部旋涡非常稳定；进口射流剪切层存在“马蹄涡”的周期性脱落；射流在燃烧室通道内卷绕形成涡强较大的二次流螺旋柱状涡对，涡对相撞又使螺旋涡失稳、振荡和破碎。旋涡不稳定性是侧面突扩燃烧室燃烧振荡的流体力学原因，分流     

轴对称突扩直管流场计算

本文是为喷管喉部沉积问题的研究而进行的喷管突扩流场研究工作的一部分,重点研究了产生突扩流效应的轴对称直管流场.在等粘度和k-ε紊流模型下,分别计算了三种速度分布下的流场,和国内外相关实验结果及Spalding,D.B.的计算结果进行比较,取得了满意的结果.同时考察了动边界下流场的变化,实现了突扩流效应的理论预测.本文考察了SIMPLE方法.在源项处理、方程的非线性、方程的解法等方面做了分析研究和计算工作,解决了一些具体应用问题.     

旋翼非定常平行桨-涡干扰流场的数值模拟

建立了一个基于Euler方程的旋翼平行桨-涡干扰流场的数值模拟方法。该方法在时间离散上使用隐式的双时间法,以避免因物理时间步长过小而引起的数值不稳定问题;为减小由于空间离散格式精度和网格密度引起的涡数值耗散,采用了预定涡方法。以Kitaplioglu的平行桨-涡干扰试验模型为算例,计算了表面瞬时压强分布,并与可得到的试验结果相比较,验证了方法的有效性,且描述了干扰流场的气动力变化特性以及桨叶表面扰动压力波的产生和传播过程,然后分析了马赫数、涡强度、干扰距离等参数对桨-涡干扰流场特性的影响。在此基础上,得出了一些有意义的结论。     

旋流-突扩燃烧室冷态流场研究

本文对带有收扩段的轴对称旋流-突扩燃烧室模型的冷态流场进行了试验和数值计算.当旋流角为45°时,涡破碎是一个较弱的转变过程,流场受尾喷口收缩的影响较小.所以,带有收扩段和45°旋流角的旋流燃烧室比起其他结构形式的同轴旋流-突扩燃烧室具有类似超临界流动的特点.试验表明:当旋流角大于45°时,中心回流区反而减小.因此,对于旋流燃烧室来说,旋流强度增大,中心回流区并不一定增大.本文采用扩编的TEACH程序和SIMPLE方法对旋流流场进行了数值模拟计算.数值计算结果比较接近试验结果.     

用Tam-Thies湍流模型研究带小突片的喷管流场特性

 首先,采用标准k-ε模型、可实现（realizable）k-ε模型、重正化群（RNG）k-ε模型和Tam-Thies模型4种湍流模型模拟了不同工况下的两种轴对称喷管和带有4片小突片的喷管流场,并将计算结果与实验值进行对比。结果发现：Tam-Thies模型所得出的模拟结果与实验值符合最好。然后,用Tam-Thies模型模拟了小突片后倾角不同的喷管的流场,和没有小突片的喷管的计算结果对比,发现：带有小突片的喷管的尾喷流的核心区长度变短,喷管出口下游气体的混合加剧;而当小突片后倾角增大时,核心区长度先减小后增大;在每个小突片下游产生一对方向相反、强度相同的流向涡,当小突片后倾角增大时,流向涡强度呈现出先增大后减小的趋势,而喷管的引射系数增益和推力损失系数都逐渐减小。     

降落伞的流场特性研究

降落伞的流场特性直接关系到降落伞的气动力性能,一直以来是降落伞研究中的核心问题。本文采用计算流体动力学(CFD)方法对降落伞的流场进行数值模拟。通过采用不同的网格生成方法对降落伞的绕流流场进行研究,表明生成网格质量对流场计算精度有很大的影响。在一系列数值实验的基础上,证明采用多块贴体坐标的网格处理方法能够获得良好的网格质量,同时还具有易于寻找伞衣内外面压力,计算稳定性、可靠性高等优点。在此基础上,采用了七孔探针对降落伞绕流流场进行了测量,实验结果和计算结果吻合较好。     

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

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

共用尾喷管多管脉冲爆震发动机数值模拟研究

采用基元反应模型和显式迎风TVD差分格式数值模拟方法,研究共用尾喷管多管脉冲爆震发动机(Pulse Detonation Engine,简称PDE)各爆震管之间相互影响.采用点隐算法解决化学反应引起的刚性问题,模拟了共用尾喷管内的流场,以及爆震波退化为激波遇到尾喷管收敛斜面反射及反射波向爆震管上游传播的过程.研究结果为多管PDE的共用尾喷管设计提供了理论基础.     

螺旋式发射管内火箭燃气冲击流场数值模拟

用改进的MUSCL格式解三维、可压缩平均雷诺纳维尔-斯托克斯方程组,湍流模型为Spalart-Allmaras代数模型,计算和分析了火箭发射管内燃气冲击流场.在建立差分算法时,将有限体积离散和黎曼解算器相结合,简化了计算工作.首先,以小喉截面超声速喷管流动问题开展数值实验,数值结果与实验数据吻合很好;其次,计算了火箭发射管内燃气流场,获得发射管内气流压强和温度等物理参数分布.     

二维收扩喷管流场计算

 <正> 研究尾喷流热辐射对隐身技术发展有着重要意义。实验表明,二维喷管红外抑制效能比传统的轴对称喷管好。二维喷管流场计算是红外辐射场计算的一部分。     

再入体实验模型热化学非平衡绕流流场的数值分析

利用数值求解一个和多个振动温度热化学非平衡Navier-Stokes方程的CFD计算程序,对爆轰风洞球锥试验模型的热化学非平衡绕流流场进行了数值模拟,分析了分子组分振动温度在全流场(头身部流场和底部流场)中的分布规律.研究结果表明:(1)计算的压力、电子数密度以及流场的光辐射数据与试验数据符合较好;(2)在再入体身部,大多数分子的振动温度峰值大大高于平动温度的峰值;在再入体近尾出现振动温度冻结并大大高于平动温度的现象;CO2的两个振动态的振动温度分布非常接近平动温度分布.     

多段翼型粘性绕流计算研究

针对多段翼型外形结构复杂,结构网格生成难度大的缺点,运用割补法对重叠网格技术进行了研究,生成贴体及与边界正交的高质量的结构化网格,在识别贡献单元方面,介绍了一种基于kd树的找重方法.流场求解应用雷诺平均Navier-Stokes方程结合Spalart-Allmaras一方程湍流模型在所生成的重叠网格上对多段翼型的粘性绕流进行了模拟,为了兼顾流场计算效率问题,在重叠网格的各个子块上采用多重网格方法.计算结果表明,用该方法可以很好地预计压力分布,并且计算效率得到很大的提高.     

旋翼地面效应的气动建模与特性

地面对旋翼气动特性影响明显，且导致旋翼流场更加复杂。为分析地面效应下的旋翼桨尖涡和流场变化特性，基于涡面和无滑移边界条件，求解第2类Fredholm方程获得地面涡面矢量分布，且将涡面矢量按涡扩散方程扩散到流体中，建立考虑黏性效应的地面气动模型，并耦合非定常面元/黏性涡粒子混合法以体现旋翼桨叶气动特性和旋翼尾迹的非定常效应，构建旋翼地面效应气动分析方法。通过计算Lynx尾桨地面效应下的性能和桨尖涡轨迹，并计算Maryland大学模型旋翼和NASA缩比旋翼地面效应下的垂向、径向速度分布，且与试验和CFD计算结果对比，验证了本文方法能较好捕捉地面效应下的旋翼尾迹变化特性和复杂速度场特性，且结果表明本文方法能较好模拟地面效应下旋翼桨尖涡的收缩、扩散、井喷、地面射流等物理现象。     

三角翼受迫俯仰滚转耦合运动的气动特性研究

基于刚性动网格的技术,选用B-L湍流模型,利用有限控制体积法对N-S方程进行数值离散,对76°大后掠三角翼的受迫俯仰滚转耦合运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析。计算结果表明:俯仰滚转耦合运动时,三角翼上表面的涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大。     

雷诺数效应对小流量多级轴流压气机的性能影响

通过试验和数值模拟,对一台设计压比为6,效率0.84的小流量轴流压气机进行了雷诺数效应的分析和研究。介绍了这台压气机在设计转速下的性能曲线、不同雷诺数下折合转速n=0.85的试验结果,以及二维平面叶栅和三维数值模拟分析结果。试验和数值模拟的分析表明,随着雷诺数降低,会增加附面层厚度,改变压气机的三维流场结构,降低压气机的性能。根据折合转速n=0.85的试验和数值模拟分析结果,对这台压气机在折合转速n=1.0性能进行了雷诺数修正。