湍流附面层方程的有限差分解

本文用三点隐式有限差分法数值求解了轴对称流动的层流、转捩和湍流可压缩附面层方程,湍流附面层用两层概念处理,而每一层用适当的涡粘性模型代替Reynolds应力项,用常值湍流Prandtl数将湍流热流项与Reynolds应力关联起来。支配方程用迭代法求解。计算结果表明,采用本文改进的涡粘性模型和涡粘性分布的平滑办法,对文中所计算的几种情况,都给出了满意的结果。     

受限空间强旋射流扩散火焰的数值模拟

通过求解包括辐射换热项的二维轴对称法富尔平均的NS方程,采用涡耗散的燃烧模型,对公开实验数据的受限空间内强旋射流的无预混火焰进行了数值模拟,得到了流场的平均速度、主要标量的空间分布;通过和实验数据对比,分析了目前常见的湍流模型及数值求解技术在模拟复杂流场中的应用。     

采用新型非轴对称端壁优化设计方法提高涡轮性能的数值研究

为提高非轴对称端壁控制端壁处二次流的潜能,进一步提高涡轮性能,发展了一种新型的非轴对称端壁优化设计方法,并以高压涡轮导叶为研究对象,采用端壁参数化、三维N-S方程流场求解与基于人工神经网络的遗传算法相结合的方法进行非轴对称端壁优化,分析了优化后的非轴对称端壁造型对涡轮导叶流场的影响。结果表明:优化后的非轴对称端壁改善了涡轮导叶的流场,延迟了通道涡的生成和发展,削弱了角涡的强度,降低了导叶通道内的流动损失,涡轮导叶出口处的总压损失系数降低了3.724%。此外,非轴对称下端壁造型对高压涡轮导叶上半叶高流场的影响不大。     

超声速大迎角分离流中三种湍流模型的比较研究

通过对雷诺平均N-S方程的数值求解,模拟了细长旋成体在超声速大迎角下的流场,比较了S-A,SST和EASM三种湍流模型。不同轴向位置的物面周向压力分布和实验结果的对比结果表明:湍流模型对物面压力分布影响较大。横截面流态分析表明,空间主体涡结构受湍流模型影响不大,靠近壁面的二次分离涡的模拟精度和湍流模型密切相关。弹体积分气动力受湍流模型影响不大,和实验结果相当吻合。     

层板冷却结构流阻特性的实验与数值模拟

在吸气式风洞上对两种不同内部流动结构的单层多孔层板进行了流阻实验，采用沿程阻力关系式得出了其流阻特性。应用软件FLUENT，采用隐式有限体积法求解Navier—Stokes雷诺时均湍流方程，对层板冷却结构的内部流动进行了数值模拟，湍流模型采用Realizable κ-ε双方程模型，近壁处湍流利用壁面函数法处理，采用SIM—PLE算法求解速度与压力的耦合。计算表明，层板内部流场结构十分复杂，射流冲击后在扰流柱前反卷形成驻涡。层板结构内部存在流速很低的区域。计算得出的流阻特性与实验结果进行了比较，两者符合较好。     

机身前体横截面形状对三角翼涡破裂位置的影响

通过求解三维定常雷诺时均方程,采用剪切应力输运(SST)湍流模型,在亚声速范围内,分别对融合体前体-三角翼组合体和旋成体前体-三角翼组合体流场中翼涡破裂现象进行了数值模拟.模拟结果显示:与旋成体前体相比,在中、大迎角时,融合体前体的分离涡,涡量集中、强度高,进入机翼上方流场后,能够与翼涡密切耦合,彻底地改变翼涡强度的分...     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

平板湍流转捩过程中Λ涡和环状涡的周围流场研究

采用高精度湍流直接数值模拟技术,对平板流动转捩过程中的涡系结构,尤其是Λ涡和环状涡附近的流场进行了详细的研究.高分辨率的数值模拟给出了实验手段无法得到的全流场信息,以及雷诺平均纳维斯托克斯方程方法(RANS)无法给出的精细涡结构信息.数值模拟在精确捕捉Λ涡和环状涡等转捩过程中典型涡系结构、上喷和下扫等典型流动现象的同时,发现了一些与传统认识所不同的现象:1)Λ涡并不是当地涡量集中的地方;2)环状涡的附近存在一个速度显著大于远场来流速度的区域.这些现象对传统湍流转捩理论提出了挑战,并引发人们对湍流生成与维持机理的进一步研究.     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

侧斜射流与横向流相互干涉的试验研究

本文对侧斜射流与横向流的相互干涉流场作了详细测量.分析了流场参数的空间分布特征和沿程变化规律,对涡特征参数的变化规律及其影响因素也作了分析与讨论.根据测量结果,拟合归纳出适用于一般侧斜射流状态、速度比在1—5范围内的涡参数变化的经验公式.     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

环形燃烧室冷态流场数值模拟中的数学方法

针对带双级径向旋流器环形燃烧室火焰筒,研究不同紊流模型(standardk-εmodel,RNG(renormalization group)k-εmodel和代数应力模型(algebraic stress model))及算法(SI MPLE(semi-i mplicitmethod for pressure-linked equations),PISO(pressure i mplicit split-operator))对复杂流场求解准确性及效率的影响.火焰筒三维贴体网格采用型线定点与偏微分方程相结合的方式生成,用壁面函数法处理近壁区域流动.计算结果表明,PISO算法提高求解效率;不同紊流模型预测结果均与试验数据符合较好,其中代数应力模型最为准确,更适用于复杂流场的数值预测.      

解决湍流圆/平射流异常现象的数值研究

圆射流/平面射流异常现象是湍流研究领域的一个著名的难题,是验证湍流微分方程模型的试金石。本文用GAO-YONG不可压湍流控制方程对圆射流进行了数值验证,成功地解决了该难题,进一步证明了该方程对真实湍流的适定性和通用性。      

涵道尾桨的CFD模拟与验证

采用CFD(计算流体力学)方法分析在悬停和侧飞时直升机涵道尾桨的流场与性能。在轴对称的圆柱坐标系中,用有限容积法和SIMPLE(压力耦合方程的半隐式法)算法求解定常不可压的湍流N-S方程。在分析中,旋转的螺桨被描绘成沿螺桨桨叶展向分布的、与本地流动参数相关的、时间平均的动量源项,通过在N-S方程中加入此动量源项来替代螺桨对流体的作用。计算方法还包括涵道壁的阶梯型近似,原始变量的交错网格,k-ε湍流模型和涵道壁上的壁面函数法等措施。涵道尾桨的流场分析和性能预测同实验测量数据的良好的一致性表明,这种CFD方法可以有效的分析涵道尾桨的具体设计问题。      

区域分裂法在湍流流场数值解中的研究

应用区域分裂算法(简称DDM)对湍流流场进行数值研究。对湍流流动控制方程建立一种通用的DDM-1计算格式S_{I,ω,给出重叠区域内变量处理方法。DDM-1的SI,ω计算格式的适用性由突扩圆管分离流流场数值计算所验证,数值计算结果与实验数据比较表明:采用SI,ω计算格式所预测的再附点位置及平均轴向速度、湍动能强度比在整个计算区域内采用单一模型的方法能够更好地与实验值吻合。}      

脱体涡模拟方法研究及其应用

基于现有的非定常雷诺平均N-S方程的求解方法,研究与发展了脱体涡模拟方法（DES）。在紊流附面层内用RANS方法,湍流模型采用S-A模型,在其他区域结合S-A模型运用Smagorinski的LES模型,并将S-A模型中的物面耗散项进行适当的改进,使得改进的模型既充当了RANS中的S-A紊流模型,又充当了Smagorinski大涡模拟模型。为验证所发展方法的有效性,将研究的DES方法应用到一些典型流场计算中,并与已有结果或实验结果进行比较。     

应用S-A模型的自由射流和冲击射流数值模拟

为探索轴对称射流准确的数值模拟方法，根据可压缩的轴对称N-S方程，利用二阶精度的有限体积法对轴对称的自由射流和垂直冲击射流进行数值模拟。结果表明采用Spalart-Allmaras湍流模型，可很好描述射流流动。收缩喷嘴射流计算应从出口上游的内流开始，而不能在出口截面上给定声速条件。计算结果与实验吻合较好，所获得的喷嘴内外部流场的流动状况以及冲击射流参数分布，可用于优化喷嘴内部结构设计。     

受限共轴射流的湍流混合过程数值模拟

本文以下Ｓｔｉｒｌｉｎｇ发动机引射管设计为研究对象，对亚声速受限共轴射流中的两种气体混合过程，用数值模拟方法求解了轴对称雷诺平均Ｎ－Ｓ方程组。提出的一种改进的显式有限体积法，包括基于网格角点的具有高分辨的中心／迎风ＱＵＩＣＫ算法，四阶人工耗散、局部时间步长和隐式余数光滑技术的应用。为加速收敛，采用一种类似于总焓阻尼的方法，引入质量流量阻尼和气流中化学组份质量分数阻尼，以快速获得足够近似的解，适合工