Numerical study of transient interchamber processes when reaching the SPRE operational conditions

A comprehensive numerical study of transient interchamber processes occurring when reaching the SPRE operational conditions is carried out. A conjugate problem under consideration includes nonstationary operation of an igniter heating, ignition and subsequent nonstationary and turbulent burning of a solid-propellant charge; nonstationary three-phase homo-heterogeneous flow of combustion products in the combustion chamber, nozzle and downstream of the nozzle block of a rocket engine; motion of a nozzle block cover. The calculation results are presented.     

离心叶轮内稀释固液两相流动的数值模拟

 对固液两相流动进行了数值模拟。液相和固相分别采用 k- ε湍流模型和代数颗粒流( AP)湍流模型并考虑了液固两相间的相互作用。用此模型结合 SIMPLE- C算法对离心叶轮内的稀释固液两相流动进行了计算并研究了颗粒对流场的影响。计算结果同试验结果基本一致     

数值模拟纵向隔热屏加力室热态流场

数值研究纵向隔热屏加力燃烧室的在三维贴体坐标系下燃烧整体流场,所用的数学模型有:k-ε模型、EBU-Arrehenius和二阶矩紊流燃烧模型及六通量模型等.数值分析两种不同几何形状、进口气流参数分布及紊流燃烧模型等对加力室内各气流参数、隔热屏和加力室筒体壁面温度分布的影响,计算结果与试验数据比较表明:不同几何形状对加力室紊流燃烧流场的影响要比进口气流参数分布大些,二阶矩模型更适用模拟紊流燃烧流动.      

Numerical simulation of a boiling up fluid flow in the variable cross-section channels

A physical and mathematical model of a boiling-up fluid flow in the variable crosssection channels is proposed. The numerical model is based on the solution of the Navier-Stokes nonstationary equations for compressible fluid by the implicit finite-difference method. We present the results of the numerical modeling for two-phase flow characteristics at the initial stage of the evaporation process and the steady-state quasistationary flow of a two-phase medium.     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

中心进气的旋转盘平均换热特性研究

实际发动机的气冷涡轮盘结构被简化成具有中心进气的旋转盘模型, 以实验方法研究了旋转雷诺数Re、进气流量系数C_w、间隙比G和出气间隙比G_s对旋转盘表面平均换热特性的影响。在本实验范围内, Re较大时, 它对盘面的Nu有较为明显的影响;C_w的增大会提高转盘表面的换热能力;G和G_s对盘面的N_u影响较小。根据实验结果给出了一个计算旋转盘表面平均对流换热系数的经验关系式。      

紊流风场下的无人机飞行状态卡尔曼滤波估计

为了模拟紊流风场下的无人机飞行状态响应,减小紊流风场对无人机飞行速度的影响,使用数值法建立了紊流风场有色噪声模型,生成了符合大气紊流相关特性的紊流风场并对其进行数值仿真。在此基础上,通过状态扩增处理将系统噪声白化,设计出一种基于有色噪声的无人机飞行速度卡尔曼滤波器,该滤波器解决了基本卡尔曼滤波仅适用于白噪声的情况。仿真结果表明,使用此滤波器可有效减小紊流风场对无人机飞行速度的影响,进而满足飞行速度控制输入的精度要求。     

平板射流中大尺度涡的理论模型

根据流动稳定性理论,提出了一种解释平板射流中大尺度涡产生机理的二维不稳定波理论模型,用这种模型计算所得的流动结构和相干扰动量与实验结果吻合较好.说明不稳定波模型可以较好地描述湍流平板射流中的大尺度涡结构的主要特征.      

转捩对平流层飞艇外形优化的影响

采用梯度算法优化了平流层飞艇外形以减少气动阻力,目标函数飞艇总阻力系数通过求解Navier-Stokes方程获得.优化过程中采用了SSTκ-ω转捩湍流模型和Realizable κ-ε全湍流模型求解绕流流场,采用大涡模拟方法详细分析了获得的两种飞艇外形的气动性能和特点,并且比较了两种湍流模拟方法的流场求解对梯度法优化的目标函数以及飞艇形状特点的影响.结果表明:是否考虑转捩对优化的结果有明显的影响,考虑转捩优化方案的气动性能明显优于全湍流流动求解优化方案.采用转捩模型的优化由于层流区域的存在改变了飞艇尾部的形状,从而获得了尾部分离区较小的飞艇外形轮廓.      

GAO-YONG湍流模型对湍流传热的数值模拟

将GAO-YONG湍流模型应用于湍流传热的研究,分别计算了平板剪切湍流和二维平面冲击射流的湍流传热问题.边界层剪切湍流流动与换热的计算表明:与传统的湍流模型不同,GAO-YONG湍流模型不需要对近壁区域做任何特殊处理(比如壁面函数、低Reynolds数修正等)即可模拟出从壁面到主流区的全部流动与传热情况;另外,对于冲击射流Nusselt数的模拟也得到了与实验符合较好的计算结果,准确地捕捉到了2种冲击高度下流场换热的不同特征,表明了GAO-YONG湍流模型能够较高精度地计算湍流换热.      

k-ω SST两方程湍流模型中参数影响的初步分析

Menter的k-ωSST两方程湍流模型在流体力学计算中有良好表现。利用均匀试验设计方法对k-ωSST模型中八个参数在小攻角亚声速、跨声速和大攻角三种典型的翼型绕流流场计算中的影响规律进行了分析。分析结果表明在附着流中参数a1是对流场的影响最大的参数;而在分离流中,参数a1对流场的影响与湍流产生与耗散等项的作用相当,模型中的四个参数:σω2、β2、β*和a1的取值都会对计算结果有较大的影响。     

二维超音速湍流混合特性研究

着重讨论顺流喷射情况下的超音速混合的流场特征及影响因素。简要介绍了顺流喷射混合的研究现状。同时，探讨了不同的来流参数对流场特征的影响及混合机理。计算采用了一种二阶ＴＶＤ格式有限体积法及Ｂ－Ｌ湍流模型求解二维湍流Ｎ－Ｓ方程。结果表明：所采用的计算方法能成功地抑制非物理振荡，适用于二维超音速内流的湍流混合特性研究。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

贫燃料湍流预混燃烧火焰稳定性的数值研究

对突扩燃烧室内甲烷-空气预混燃烧进行了数值模拟,时均控制方程组的封闭采用RNG 湍流输运模型和premixed combustion模型.通过模拟研究了不同当量比及进口速度对甲烷-空气预混燃烧效果的影响.CFD数值模拟结果给出了燃烧室内湍流预混反应流的速度场及温度场分布.模拟结果表明当量比为0.6、进口速度为30 m/s时燃烧室内产生稳定预混火焰.该结果对于航空燃气轮机低NOx排放燃烧技术的研究提供了一定的参考条件.     

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.      

不可压缩湍流的多尺度模型

基于可变时间间隔平均方法,提出了一种不可压缩湍流多尺度模型及平均流动方程.与传统的单尺度湍流模型不同,该模型在建立的过程中,保留了湍流的多尺度特性,结合平均流动方程,可以更好地预测湍流流场特征.通过模拟平面后台阶流动和不对称平板扩压器流动,并将预测结果与标准k-ε模型的预测结果对比,初步验证了其可信性及优越性.结果表明:计算所得的平面后台阶流动的流向再附长度与台阶边压力系数比比标准k-ε模型提高精度约20%;平板扩压器流动的回流区位置误差约为7%、倾斜壁面摩擦因数误差约为5%,而标准k-ε模型未能预测出分离现象.可以看出该模型适用于典型的分离流动,在湍流流场的预测中表现优异,具有一定的工程应用价值      

非线性 k-ε 紊流模型在平面叶栅紊流计算中的应用

 采用Speziale的非线性k-ε模型计算了平面叶栅内的紊流流动,并且与标准的k-ε模型的计算结果作了对比,表明非线性模型同样可以应用于具有周期性边界条件的复杂流动,在紊流量的计算结果上表现出了明显的各向异性。     

槽流的亚临界转捩与局地湍流

在中低雷诺数时,处于线性稳定的槽道流若受到强扰动可发生亚临界转捩变为湍流。近年人们对该类转捩机理的研究取得重要进展。对于平板泊肃叶流,亚临界转捩之初是稀疏湍流态,其特征结构是远间隔的包含小尺度涡和高低速条带的大尺度的湍流带,可倾斜伸长。该阶段的湍流占比有上限但并非雷诺数的单值函数。随着雷诺数的增加转捩进入平衡局地湍流态,即存在统计定常态,其湍流占比是雷诺数单值函数,可由定向逾渗模型描述。进一步增加雷诺数,湍流带的分裂愈发频繁,最终流场会布满湍流带,在更高的雷诺数时变为均匀湍流。论文概述了为比拟亚临界转捩过程所提出的动力学模型,以及为定量表征管流、平板库艾特流和平板泊肃叶流转捩过程的相似性所提出的局地稳定性参数,并在最后对槽道流亚临界转捩研究的发展做了简要展望。     

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。     

来流条件对SFRJ燃速及自持燃烧性能的影响

为了研究来流条件对固体燃料燃速以及自持燃烧性能的影响规律,采用3阶单调迎风格式(MUSCL)重构方法,AUSMPW⁺通量分裂格式,k-w剪切应力输运(SST)湍流模型,7组分3反应有限速率化学反应模型以及2阶矩湍流燃烧模型,编制了二维轴对称湍流燃烧仿真程序。仿真结果显示:来流质量流率和总温的增加会提高燃速,温度和质量通量的影响指数大约为0.95~1.00和0.67;来流条件影响燃速的主要机理是壁面附近温度和湍流黏性的变化,其中湍流黏性是主要影响因素;随着质量流率的增加,特征速度逐渐降低,当入口马赫数接近1时,特征速度大约为850 m/s,固体燃料将无法自持燃烧。