Isothermal study of effusion cooling flows using a large eddy simulation approach

An isothermal numerical study of effusion cooling flow is conducted using a large eddy simulation(LES) approach.Two main types of cooling are considered,namely tangential film cooling and oblique patch effusion cooling.To represent tangential film cooling,a simplified model of a plane turbulent wall jet along a flat plate in quiescent surrounding fluid is considered.In contrast to a classic turbulent boundary layer flow,the plane turbulent wall jet possesses an outer free shear flow region,an inner near wall region and an interaction region,characterised by substantial levels of turbulent shear stress transport.These shear stress characteristics hold significant implications for RANS modelling,implications that also apply to more complex tangential film cooling flows with non-zero free stream velocities.The LES technique used in the current study provides a satisfactory overall prediction of the plane turbulent wall jet flow,including the initial transition region,and the characteristic separation of the zero turbulent shear stress and zero shear strain locations.Oblique effusion patch cooling is modelled using a staggered array of 12 rows of effusion holes,drilled at 30° to the flat plate surface.The effusion holes connect two channels separated by the flat plate.Specifically,these comprise of a channel representing the combustion chamber flow and a cooling air supply channel.A difference in pressure between the two channels forces air from the cooling supply side,through the effusion holes,and into the combustion chamber side.Air from successive effusion rows coalesces to form an aerodynamic film between the combustion chamber main flow and the flat plate.In practical applications,this film is used to separate the hot combustion gases from the combustion chamber liner.The numerical model is shown to be capable of accurately predicting the injection,penetration,downstream decay,and coalescence of the effusion jets.In addition,the numerical model captures entrainment of the combustion chamber mainstream flow towards the wall by the presence of the effusion jets.Two contra-rotating vortices,with axes of rotation along the stream-wise direction,are predicted as a result of this entrainment.The presence and characteristics of these vortices are in good agreement with previous published research.      

基于湍流边界层时均速度分布的脊状表面减阻规律研究

通过对试验测得的脊状表面湍流边界层时均速度分布的分析,获得了脊状表面的减阻范围和基本规律.脊状表面理论零点及壁面摩擦速度的计算采用基于湍流边界层Spalding壁面律公式的最小二乘拟合方法;模型板摩擦阻力计算采用基于边界层动量积分方程的方法.研究表明:脊状结构使得湍流边界层粘性底层增厚,近壁面法向速度梯度降低,过渡层与对数律区上移;V型脊状结构的无量纲宽度s+介于6至18时具有减阻效果;当s+≈12时,减阻量最大,最大值约7.7%.      

非结构网格下含冷却孔的涡轮转子三维流场计算

在转动坐标系里,讨论非结构网格下,含多排冷却孔的涡轮转子三维流场雷诺平均Navier-Stokes方程组的求解。计算表明:含冷却孔的转子流场计算要比静子流场复杂的多。为了得到较为满意的流场数值解,本文对叶片压力面与吸力面采用了不同的湍流模式、并且还对流场迭代中初始场选取方式以及转子叶尖间隙计算站的设置等问题进行了一系列的改进。几个典型流场算例表明:初始场选取、湍流模式的合理选用对转动坐标系下获得一个收敛的三维流场数值解至关重要;另外,转子叶尖间隙站的设置直接影响着间隙附近的流场,尤其对现代高负荷涡轮转子,间隙流动的研究就更为重要。      

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

三角翼跨声速动态失速与涡破裂特性研究

通过数值方法研究了高速气流中细长三角翼作匀速上仰机动飞行时的背风面分离流与涡破裂特性和气动力特性。结果表明,在三角翼匀速上仰非定常绕流中,涡破裂起始攻角与上仰速度关系比较复杂,在本计算速度范围内,涡破裂起始攻角变化很小;动态时涡破裂点发展速度随攻角变化规律与静态结果差别甚大;三角翼大攻角非定常绕流前缘涡涡轴附近的流动以不稳定涡和涡破裂为主要特征;在一定的上仰速度以后即使在中等攻角前非定常升力与定常升力仍有较大差别,与低速情况有所不同;上仰运动使涡破裂点发展速度被延缓,从而提高了失速攻角和最大升力,并使失速攻角与涡破裂起始攻角的间隔进一步拉大,同时不改变升阻比     

绕三角翼旋涡流动的数值模拟

本文通过高清晰度的Ｎ－Ｓ方程数值模拟对绕大后掠三角机翼的旋涡流动结构及旋涡破裂进行了研究,给出了机翼背风面上复杂的旋涡流动结构,绘制了流动横截面上的截面流线。算例的计算采用了隐式近似分解格式。     

使用Gao-Yong湍流方程计算翼体角偶流动

Gao-Yong湍流方程基于侧偏统计平均方法保留了湍流脉动量的一阶统计信息,并引入加权漂移速度对称性及正交各向异性,导出漂移流连续方程、动量方程和机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组.大量算例验证了此方程对广泛范围的复杂湍流问题的适定性.为了进一步的工程应用,将其加入开源计算流体力学软件OpenFOAM(Open Field Operation and Manipulation)中,并就翼体角偶流动进行了数值模拟.翼体角偶流动具有三维分离、马蹄涡等十分复杂的流动特征,通过与实验测量结果的对比研究表明,Gao-Yong湍流方程不论在定性还是定量上,都成功捕捉到这类复杂流动现象的主要特征,得到了令人满意的结果.     

坝陵河大桥桥位深切峡谷风剖面实测研究

用相控阵声雷达风廓线仪对坝陵河大桥桥址处深切峡谷中风剖面进行了实地观测,结果表明:在地形复杂的山区深切峡谷中,平均风剖面受峡谷地形影响较大,形态有时显得较为复杂而不规则,呈现锯齿形、正切变形或逆切变形等多样化.在峡谷的上半部,平均风速较大时其剖面相对较为规则,但不符合规范中描述平坦地貌平均风速剖面的幂函数形式,而具有e指数函数的变化规律.此外,统计显示,风向角在低空范围内受峡谷地形的影响要比在高空范围内严重,其在低空范围的变化幅度也要比在高空范围的变化幅度大.风迎角的散布范围和绝对值大小均随高度的增加呈减小的趋势.观测得到的峡谷湍流度大于平坦地貌湍流度的规范推荐值,同时峡谷湍流度具有一定的随机性,并且这种随机性随高度的增加而变大.     

进气畸变对燃烧室性能的影响

基于详细化学反应动力学和火焰面模型,以某型发动机燃烧室为计算模型,在压气机气流畸变来流条件下对燃烧室内流场进行了数值模拟,将计算结果与其实验结果进行了对比,并对比研究了均匀与非均匀入口气流条件下燃烧室燃烧性能的差异。结果表明:非均匀入口效应会对燃烧室内部流场产生很大的影响,流场的对称性受到破坏,燃料的掺混燃烧效果也会受到影响,进而导致燃烧室燃烧效率降低,影响整个燃烧室的工作性能。增大主燃孔面积,提高头部进气量可以很好地改善此型燃烧室的燃烧性能。      

螺旋管内超临界CO2流动方向对换热的影响

应用重正化群（RNG）k-ε湍流模型对超临界CO₂流体在内径为9 mm,有效受热长度为5.5 m,节距为32 mm,绕径为283 mm的竖直螺旋管中的加热过程开展了数值模拟。研究了质量流速、进口压力、热通量以及不同流向对超临界CO₂换热和压降的影响,并进一步分析变物性、浮升力和离心力在不同流动方向上对螺旋管中换热的耦合作用。结果表明：浮升力对超临界流体在螺旋管向上流动与向下流动的影响差别不大,而对水平流动的影响较大尤其是当流体在螺旋管的截面到入口截面的距离与管径之比为150～350之间（即临界温度附近）时,由于变物性、浮升力和离心力的耦合作用,导致水平流动方向上换热系数的震荡。