平板附面层等离子体流动控制的数值模拟

通过求解电场中的拉普拉斯方程和赫姆霍兹方程的变形形式,成功地将等离子体激励对平板附面层流动的影响,以体积力向量的形式引入到NS方程之中.借助求解的电荷密度,通过将计算结果与实验数据的对比,可确定德拜长度、最大电荷密度、形状因子等可调参数的取值原则,建立起平板附面层等离子体流动控制的数值模拟方法,为将等离子体流动控制方法应用于外流及内流场中的强剪切流动控制,奠定关键性的技术支撑.      

驻室式处理机匣——内部流动及扩稳机理的数值研究

通过定常模拟的方法,研究了驻室式处理机匣对高负荷、高马赫数孤立转子性能的影响.结果表明:相同的处理机匣结构在不同的转速下取得的效果存在较大差异,在设计转速时,处理机匣并没有使风扇转子的失速裕度得到改善,并且效率下降约2.8%;在非设计转速时,处理机匣扩稳作用得到充分体现,转子稳定性得到较大程度的加强,并且没有导致转子性能的下降.通过对处理机匣内部流场的详细分析,揭示了驻室式处理机匣的扩稳机理.      

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。      

微通道气体单相流动特性的数值分析

微通道是各种微系统结构中一个基本的组成部分.考虑了可压缩与粘性加热的影响, 针对氮气在3种不同尺寸的微通道(水力直径范围:30～3 000μm)内的流动特性进行了数值研究.分别就不同尺寸、不同雷诺数(Re)下压力降、摩擦因子和摩擦常数的规律进行了研究, 着重考察了长径比、入口段出口段影响、可压缩性影响及粘性加热影响.计算分析揭示出转捩雷诺数提前, 在300μm管道转捩雷诺数提前到1 200附近.      

碳氢燃气超声速剪切流动燃烧数值模拟

数值模拟碳氢富燃燃气与空气二维超声速剪切流动的湍流燃烧.NND对流项差分格式高精度捕捉激波, B-L和普朗特混合长度模型模拟湍流粘性系数, Arrhenius反应动力学公式计算化学反应速率, 模拟非平衡态燃烧.计算的壁面压力值较好地符合了实验值.分析了剪切流驻焰稳定的波系结构, 从预测的流动参数和组分浓度分布判断燃烧状况的好坏, 从而为改善优化燃烧室结构提供工程参考.      

风机结构形式对其性能影响的数值与实验研究

通过数值模拟和实验研究,对不同结构形式的轴流风机进行了对比分析。找出了由于风机生产条件限制而造成的几种不当结构影响风机整体性能的原因,以及风机中主要损失产生的位置。结果表明,在风机原始设计过程中,应从整体优化的角度出发,预先考虑工厂中产品实际配装的几何结构形式,然后调整相关气动设计参数,从而要使气动设计与产品实际结构相匹配。      

超声速气流中横向射流雾化实验和数值模拟

为了对超声速气流中的横向射流雾化过程进行研究,在直连式实验台上进行了超声速气流中的横向射流雾化实验,并对实验工况进行了数值模拟.超声速来流的马赫数为2,实验喷雾的动压比范围为1.～11.7,实验工质为水.采用纹影法对射流雾化过程进行了拍摄,得到了有雾化的超声速流场结构和穿透深度拟合公式.同时采用欧拉-拉格朗日两相流计算方法对实验工况进行了数值模拟,雾化模型采用了一种新型K-H(Kelvin-Helmholtz)和R-T(Rayleigh-Taylor)混合模型,数值模拟与实验结果符合较好.      

有预旋进气转静盘腔中的流动和换热特性数值研究

运用RNG湍流模型对有静盘外缘预旋进气和轴向中心进气、转盘外缘轴向出气和外围屏径向出气的转-静盘腔中流动和传热过程进行了数值研究。在盘腔结构中,静止内隔片和旋转内隔片将盘腔分为内转-静盘腔室和外预旋腔室。由于盘腔流动结构具有周向周期性,取盘腔的1/30作为计算域。研究发现:旋转雷诺数Re_θ、预旋喷嘴进气无量纲流动速度C_w,_p、静盘中心轴向进气的无量纲流动速度C_w,_d和预旋比β_p都不同程度地影响盘腔内流动和换热。预旋和隔片的协同作用可以很好地起到阻止热燃气入侵转-静盘腔和提供低静温冷却气体流入涡轮叶片冷却通道。      

再生冷却推力室耦合传热数值模拟

对大推力液体火箭发动机再生冷却推力室内部燃气、室壁和再生冷却剂进行了耦合传热数值计算.采用二维轴对称N-S方程描述推力室内部燃气的湍流流动与传热,对冷却剂流动采用简化的一维模型,通过室壁的偶合传热采用一维传热模型.N-S方程的求解采用贴体坐标系下的有限容积法,速度和压强的耦合采用可压缩的SIMPLE算法.湍流的模拟采用可压缩的标准κ-ε模型,辐射传热采用热流模型计算.研究表明,本文方法可较好地模拟燃气的二维流动,同时能快速计算壁面热流密度、壁温和冷却套温升,计算结果对大推力发动机推力室的设计具有一定的指导意义.     

双级风扇三维粘性非定常流动的数值解

针对一台双级跨音速轴流风扇动静叶片排相互干扰的三维非定常流动,进行了N-S方程的数值求解。首先计算风扇定常流场,以定常流收敛结果为起始场计算非定常流动。计算结果表明,非定常流动对各叶片排进口气流角有较大的影响,其中静叶的变化大于动叶,与定常流结果相比,非定常的瞬态攻角增大的量大于减少的量,最多时增大近10°左右。静子叶片受到的非定常气动力变化幅度也大于转子叶片,叶片排所处的轴向位置不同,非定常流动的影响会有较大的差别。      

多股自耦合射流相干流动和换热特性数值研究

采用数值模拟的方法对多股自耦合射流在相干条件下的流动和换热特性进行了研究.研究结果表明:三股自耦合射流都在出口附近耦合成一股射流向下游运动,两侧的射流受到中间一股射流的卷吸融合在一起,融合后仍然是一股沿中心线对称的自耦合射流.在短轴中心截面自耦合射流型面较宽,射流喷出后在下游的扩展角度大,而长轴中心截面射流的型面较窄,扩展角度较小.在冲击靶板上,整个冲击区域呈现出中心温度最低的近似圆形的对称性分布,对流换热系数曲线呈现出中心对称的波峰分布形状,表明自耦合射流驻点区换热最强.对流换热系数沿靶板不同方向的分布趋势基本保持一致.      

竹节形扰流元对流动与换热特性影响的数值研究

采用数值计算方法, 研究了扰流柱形状对通道中流体流动和换热的影响, 模拟了带竹节形扰流柱通道的流场.结果表明:相对于圆柱形扰流柱, 装有圆柱竹节形和椭圆竹节形扰流柱阵列矩形通道的压力损失系数分别为62%～77%和25%～27%, 而恒热流壁面的平均奴塞尔数分别降低了9%22%和22%24%.随着扰流柱排列间距的减小, 压力损失增加, 而改变扰流柱排列横向间距引起的压力损失的变化要比改变流向间距显著.      

减阻表面换热特性的实验

采用实验方法对带"V"型沟槽表面通道在不同雷诺数下的流动与换热特性进行了研究.研究结果表明:宽与深均为0.5 mm的沟槽板通道在实验范围内具有明显的减阻效果, 减阻效果随着雷诺数的增大而减小;相比平板通道流动, 低雷诺数范围下(10000相似文献   

微阵列射流冲击肋化表面传热特性研究

本文以去离子水为工质,对微阵列射流冲击光滑和肋化表面的传热特性进行了实验研究,并以光滑表面传热性能为基准,讨论了肋化表面强化因子(ε)变化规律。实验结果表明：(1) 浸没和自由射流的传热系数均受无量纲射流距离(H/d)影响;提高射流Re数和减小无量纲孔间距(S/d)都能够增强换热。(2)浸没射流ε受H/d影响大,而自由射流ε基本不受H/d影响;两种射流方式的ε都随Re数增大,且渐趋于一个常数。(3)对流热阻在各种射流情况下均随流量不断下降,但下降趋势逐渐平缓,当流量大到一定时,热阻基本不再降低。     

三种不同结构超声速燃烧室流场的数值研究

采用有限体积法求解可压缩N-S方程及组分方程,对三种不同结构的超声速燃烧室进行了数值研究.研究结果表明,凹槽改变了燃烧室波系结构,促进氢气和空气的混合,利于燃烧,所形成的相对稳定的回流区可用于点火和稳定火焰;喷射二次氢气可消除仅一次喷氢所带来问题,进一步促进氢气和空气的混合及燃烧,但带来比较大的总压损失.      

带三角形突片气膜冷却结构换热特性的数值研究

运用RNG k-ε湍流模型对带平行于冷却壁面、且覆盖部分气膜孔的三角形突片气膜冷却结构,在不同吹风比和突片几何参数下的换热特性进行了研究。结果表明:在吹风比小于1时,突片对提高气膜冷却有效温比和增强换热的作用不大;当吹风比大于1时,突片对于提高气膜冷却的有效温比和换热增强比具有显著的作用,而且随着突片尺寸的加大,气膜冷却有效温比和换热增强比的提高幅度越大。      

扰流柱直径对倾斜出气孔层板内流动和换热影响的数值研究

用数值模拟的方法, 研究了倾斜出气孔层板结构最小单元体内扰流柱直径对于层板内流阻和换热的影响.采用了非结构化网格, 并对气流转角较大区域做了局部加密;选择k-ε湍流模型方程, 配以壁面函数率, 求解稳态无旋转的层板;采用了流体域与固体域耦合求解法计算耦合换热.给出了计算结果并分析了流动和换热情况.研究结果表明适当增加扰流柱直径有利于降低流阻、提高冷却效率.