数值模拟入射斜激波／平板湍流边界层干扰流动

应用GAO-YONG可压缩湍流模型数值计算了入射斜激波／平板湍流边界层相互干扰现象。计算程序中的对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程。计算中包含了无分离流动、初始分离流动以及较大分离流动等多个情况,比较了平板壁面压力、法向平均速度剖面、壁面摩阻系数Cf以及壁面斯坦顿数St等的计算结果与实验值。结果发现：GAO-YONG可压缩湍流模型能够很好地预测入射斜激波／平板湍流边界层相互干扰下的无分离以及小分离流动,对高马赫数下的大分离流动也能得到较合理的结果,但再附点之后的壁面摩阻系数以及斯坦顿数的计算值不够理想。     

旋转矩形通道内湍流雷诺应力分析

采用直接数值模拟方法对旋转矩形通道内的湍流流动和换热进行了研究.非稳态N-S方程的空间离散采用二阶中心差分法,时间推进采用二阶显式Adams-Bashforth格式.通过分析近壁湍流雷诺应力输运方程各项的分布以及比较壁面附近湍流生成和耗散强度,讨论了科氏力对湍流脉动的影响.结果表明:系统旋转效应增强了压力面附近的湍流生成、黏性耗散、黏性扩散、湍流扩散及压力应变相关项和压力速度扩散项,而在吸力面附近,雷诺应力输运方程中各项因旋转而受到强烈抑制.     

K-∈湍流模型下门桥水阻力的数值计算

采用RANS方程和RNGκ-ε二方程湍流模型及VOF算法对绕门桥的自由表面流动进行了数值计算。数值计算时，使用有限体积法离散控制方程，对流项采用二阶迎风差分格式，扩散项采用中心差分格式。对门桥模型进行了一系列速度下的数值模拟；计算结果与试验测量数据相比较，该计算方法具有较好的计算精度，可用于门桥性能分析和优化。     

轴流泵水力模型内部流场数值模拟

基于Reynolds时均N-S方程,采用标准k-ε湍流模型,压力、速度耦合使用SIMPLEC算法,离散采用具有二阶精度的隐式格式差分,对轴流泵过流部件内部流场进行三维定常湍流数值模拟,得到泵内流动的速度和压力矢量分布图,以及其他一些流动的信息.数值模拟结果表明,设计工况下的流场分布和流态总体较好,但叶轮叶片背面进口靠近轮缘处出现局部低压,导叶外壁区域负荷大,说明该轴流泵水力模型还有进一步改进和对其性能进一步提高的必要.     

轴流式叶栅内部流动变工况性能的计算分析

应用有限元方法求解不可压黏性流动Ｎ－Ｓ耦合方程,对ＮＡＣＡ叶型的叶栅在不同冲角的条件下,叶栅内部的流动作了较为完整的计算分析。结果包括叶片表面的压力分布、叶片表面法向的速度分布以及表面摩擦系数。在计算中,采用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型和壁面函数以提高计算的速度。结果表明：冲角不仅对叶片表面的压力分布,而且对叶片附近的流动状态,特别是在叶栅流道的后半部,有很大的影响。     

燃气轮机涡轮叶片气场耦合仿真研究

采用有限差分方法开发了气冷涡轮气弹耦合求解器。该求解器分别在流动区域求解时间平均N-S方程,在固体区域求解振动方程,在流固边界则施加了气弹耦合边界条件。在流动问题求解中,采用了B-L代数湍流模型封闭时均N-S方程,采用三阶AUSMPW 差分格式离散对流项,并采用隐式LU-SGS格式求解离散后的代数方程;振动问题则采用显式格式求解;采用代数网格法进行固体区域动态网格生成。采用该耦合求解器对某涡轮级动叶的振动响应问题进行了数值仿真,研究表明：所开发的耦合求解器能够用于振动响应问题的分析,同时本算例中涡轮动叶的振动不存在发散现象,并且叶片振动对流动的影响很小     

燃气轮机涡轮级气弹耦合仿真研究(英文)

采用有限差分方法开发了气冷涡轮气弹耦合求解器。该求解器分别在流动区域求解时间平均N-S方程。在固体区域求解振动方程,在流固边界则施加了气弹耦合边界条件。在流动问题求解中,采用了B-L代数湍流模型封闭时均N-S方程,采用三阶AUSMPW+差分格式离散对流项,并采用隐式LU-SGS格式求解离散后的代数方程;振动问题则采用显式格式求解;采用代数网格法进行固体区域动态网格生成。采用该耦合求解器对某涡轮级动叶的振动响应问题进行了数值仿真,研究表明:所开发的耦合求解器能够用于振动响应问题的分析,同时本算例中涡轮动叶的振动不存在发散现象,并且叶片振动对流动的影响很小。     

旋转对平板通道内湍流流动和换热影响的直接模拟

直接模拟了展向旋转数等于0~15时平板通道内的湍流流动和换热,分析了科氏力对湍流平均特性、湍流统计特性、雷诺应力输运和湍流结构的影响.采用有限差分法离散非稳态N S方程,湍流雷诺数和普朗特数分别为150和0 71,网格数为64×128×64.计算结果表明,随着旋转数的增加,压力面附近湍流强度增大,换热增强,而吸力面附近换热减小.雷诺应力输运方程中的旋转项在各应力中的能量分配作用远大于压力分配项,且旋转项使得雷诺应力输运方程的各项分布和相对关系发生了较大的变化.由于科氏力的作用,压力面附近的条带间距减小,条带变细、变短,而吸力面流动脉动衰减,有层流化的趋势.     

球形汇聚-反射激波作用下湍动能方程的模化

球形激波与湍流间的相互作用在惯性约束聚变(ICF)、超新星爆发、冲击波碎石等许多工程问题中有重要意义.本文通过直接数值模拟,研究了球形激波的汇聚与反射诱导的湍流及其相互作用问题,分析了这一流动的特点和湍流的产生机制.重点讨论了该流动的湍动能系综平均方程中各项的作用和起主导作用的物理机制,发现传统的剪切湍流的产生机制和黏性项量级普遍偏小,激波才是诱发湍流的主要因素,因而在湍动能方程中起主导作用.由于激波的主导作用,压力脉动对湍动能方程有重要影响,与压力脉动有关的各个项在湍动能方程中均起重要作用.一些在一般湍流中可以忽略的项,例如质量流量项、压力体胀项等,在球形汇聚-反射激波与湍流相互作用问题中起重要作用.本文还探讨了涡黏性模式的有效性,发现满足可实现性原理的涡黏性模式可以较好地模拟雷诺应力.同时,本文揭示了现有湍动能方程封闭模式在本问题的应用中存在很大误差,并尝试对这些项提出新的封闭模式.     

极小展弦比后加载叶栅气动特性三维数值分析和实验

以数值手段分析了极小展弦比下后加载透平叶栅内的定常3D粘性流动及损失，并对该透平叶栅在不同出口Ma和进气角下进行平面叶栅吹风实验．数值方法采用3D可压缩粘性流动压力修正算法求解N-S方程及Baldwing-Lomax代数湍流模型．通过数值计算分析和平面叶栅吹风实验研究表明，后加载技术应用于极小展弦比叶栅中，可以起到控制3D叶栅损失与二次流生成和发展的作用．     

单孔堵塞对均速管流量计测量精度的影响

对均速管流量计单孔堵塞后的内部流动进行数值模拟,得到压差,然后算出流量系数并与正常结果相比较,最后分析孔堵塞时对测量精度的影响。计算中引用k-ε湍流模型,应用有限体积法对控制方程离散和求解,在前处理软件GAMBIT中将计算区域划分为236多万个计算单元。     

大曲率结构化表面湍流流场特性对比研究

湍流模型选择对于小尺寸大曲率圆管内软性磨粒流湍流流场模拟计算具有重要影响.针对此问题,提出使用标准k-ε模型与考虑两相磨粒流平均流动中的旋转及旋流流动情况的重整化群(RNG)k-ε模型进行仿真对比研究的方法.基于颗粒动力学理论的欧拉模型,对软性磨粒流的控制方程进行求解.通过使用两种湍流模型,对小尺寸大曲率圆管内湍流流动进行仿真计算,得到了两湍流模型下软性磨粒流流场的速度、压力、湍动能等数值模拟结果.仿真结果对比分析表明:标准k-ε模型模拟得到的速度值与实验值存在一定误差,而RNG k-ε模型仿真的速度值误差相对较小;标准k-ε模型对湍动能的估值计算比RNG k-ε湍流模型高,难以计算有滞止点的流动过程;RNG k-ε湍流模型能更有效地模拟有大曲率带分离的湍流流动,从而证明该模型更适合小尺寸大曲率圆管内软性磨粒流湍流流动的研究.     

机舱通风口冲击射流湍流特征分析及验证

以机舱个性化通风口为背景,对自由射流加壁面冲击射流,采用几种湍流模型数值分析了10倍冲击高度、雷诺数为9,000的气流遇到冲击板后湍流脉动、流动转折等流动特征.探讨了湍流模型对冲击射流关键位置平均速度、脉动速度以及雷诺应力等流动参数特征的模拟能力,通过与实验测量结果的对比,确定了湍流模型在冲击射流中的适用性和局限性.综合来看,Relizable k-ε模型可以得到与实验值较好的符合度,但在近壁面区域会略有偏差.分析结果认为将RANS和LES组合将是开展实际机舱内通风口射流研究的最佳策略.     

紊流脉动对生物反应器运行影响的数值模拟

为研究紊流脉动对生物反应器运行的影响,针对时均速度为零的振动格栅反应器,构建气-液两相流模型并对其进行验证。借助Fluent软件对反应器中气和水的流动进行数值模拟,得出10种工况下紊流脉动强度与气含率、氧转移系数和雷诺应力的关系。结合实验测得的生物膜理化数据,分析紊流脉动强度对生物膜生长的影响。结果表明,增大紊流脉动可有效提高反应器内的气含率和氧转移系数,促进气液两相进行高效传质。对于本研究所用振动格栅生物反应器,当紊流强度q<2.7 cm/s时,生物膜厚度和密度随着紊流强度的增加而增加。当q>2.7 cm/s时,雷诺应力的作用显著,生物膜逐渐变薄变密。紊流强度q=3.4~4.8 cm/s时,生物膜生长和污染物去除效果最佳。     

绕平板气液两相流数值计算方法

为探究现行多相流模型和湍流模型对通气两相流动的适应性及所体现的数理本质,选取典型的绕平板气液两相流为研究对象,通过对比不同多相流模型、湍流模型对绕平板气液两相流的数值计算结果与实验结果,建立了适用于绕平板气液两相流的数值模型。 结果表明,基于单流体理论的Mixture多相流模型忽略了流场内客观存在的气-液相界面,计算结果不能真实反映实际流动现象,与实验结果不符;Level Set方法由于考虑了气液界面的表面张力的影响,能合理反映气液界面波动的特性,可以获得和实验较为一致的计算结果;滤波器湍流模型在标准k-ε模型的基础上,通过对湍流黏性项的修正,降低了计算流场中的湍流黏度,可以更精细地描述与时间相关的气液多相复杂流动现象,进而有效提高对非定常流动计算过程的预测精度。      

絮凝搅拌器内部流场特性数值模拟

在气液交界面形态对比和LDV试验验证基础上,基于CFX 侵入式实体模型,应用标准k-ε湍流模型对容积为1 L的絮凝搅拌器内部流场特性进行了数值研究.结果表明:下循环区的流体运动速度和速度梯度较大;搅拌强度增加,速度梯度增大;转速对叶轮附近切向速度梯度影响最大,该区域切向速度梯度值与其他区域最大相差两个数量级;下循环区湍流耗散率、湍流动能和涡流黏度较高;在研究转速范围内,转子转速对湍流动能影响相对较大,对涡流黏度影响相对较小.     

栅中水翼的空化流体动力特性

采用数值模拟的方法对空化条件下Clark-y水翼和栅中水翼的水动力特性进行了研究.计算中,空化模型采用基于Rayleigh-Plesset方程的Kubota模型,湍流模型则采用基于滤波器的混合湍流模型.通过实验结果对计算结果的可靠性进行了验证.结果表明,受上下翼型的影响,栅中水翼吸力面的低压区和压力面的高压区相对水翼情况有所收缩,造成平均升力系数变小;与水翼相比,栅中水翼的升阻力变化周期更长,波动范围更大.     

应用三维代数应力通量模型模拟分层流

根据分层流紊动扩散能力在水平与垂直方向存在明显差异的特征,建立了三维代数应力通量模型。模型以有限体积法离散,SIMPLEC方法求解速度、压力和温度耦合方程组。对某分层流试验数值模拟计算表明,计算结果与实测值吻合良好;对于各向异性的分层流,代数应力通量模型与双方程k-ε模型相比,更好地保留与模拟了其基本物理特征,所得结果更趋合理。     

Mixed Model for Silt-Laden Solid-Liquid Two-Phase Flows

The kinetic theory of molecular gases was used to derive the governing equations for dense solid-liquid two-phase flows from a microscopic flow characteristics viewpoint by multiplying the Boltzmann equation for each phase by property parameters and integrating over the velocity space. The particle collision term was derived from microscopic terms by comparison with dilute two-phase flow but with consideration of the collisions between particles for dense two-phase flow conditions and by assuming that the particle-phase velocity distribution obeys the Maxwell equations. Appropriate terms from the dilute two-phase governing equations were combined with the dense particle collision term to develop the governing equations for dense solid-liquid turbulent flows. The SIMPLEC algorithm and a staggered grid system were used to solve the discretized two-phase governing equations with a Reynolds averaged turbulence model. Dense solid-liquid turbulent two-phase flows were simulated for flow in a duct. The simulation results agree well with experimental data.     

关于湍流多尺度统计理论方法的论证

用黑河化音测站的实测风、温脉动资料对一种湍流多尺度统计理论方法进行了分析验证,并应用此方法通过对相关系数、湍流时间尺度、湍谱和扩散参数等湍流特征量的研究,分析了黑河戈壁地区的湍流微结构和扩散特征.结果表明:多尺度分析方法正适合研究具有多尺度特征的大气湍流现象,可以比传统方法提供更多的有关湍流属性的信息,能够更直观地表现出不同尺度湍涡对扩散作用的大小,是发展湍流统计理论的一种有效工具