The improvement of turbulence modeling for the aerothermal computation of hypersonic turbulent boundary layers

The engineering computation of turbulent flows is mainly based on turbulence modeling,however,accurate aerothermal computation of hypersonic turbulent boundary layers is still a not well-solved problem. Aerothermal computation for turbulent boundary layers on a supersonic or hypersonic blunt cone with small bluntness is done firstly by using both direct numerical simulation and BL model,and seven different cases are investigated. Then the results obtained by the two methods are compared,and the reason causi...     

Coupled radiation and turbulent multiphase flow in an aluminised solid propellant rocket engine

We present in this paper numerical simulations of coupled radiative transfer and turbulent flows at high temperature and pressure, typical of multiphase flows encountered in aluminised solid propellant rocket engines. The radiating medium is constituted of gases and of liquid or solid particles of oxidised aluminum. The turbulent flow of the gaseous phase is treated by using a four equation, low Reynolds number, boundary-layer-type turbulence model. The distributions of concentrations, temperatures, and temperature fluctuation variances of particles are calculated from a Lagrangian approach and a turbulence dispersion model. Thermal and mechanical non-equilibrium between the gas and different classes of particles is allowed. A locally one dimensional, iteratively based, radiative transfer solver is developed to compute wall fluxes and radiative source terms. It is shown that the thermal boundary layer attenuates significantly the radiative fluxes coming from the outer regions. Particle radiation is found to be much more important than gas radiation. Turbulent dispersion of particles in the boundary layer induces a decrease of particle concentration in the region of maximum turbulent kinetic energy, and then, decreases the attenuation effect of wall fluxes due to the boundary layer. The effects of turbulent temperature fluctuations are found to be small in the problem under consideration.     

近壁区理论耗散率k-ε模型的研究

提出在湍流边界层近壁区采用共振三波的理论模型描述湍流相干结构,根据理论模型计算了ε的分布。并且在传统k-ε模式基础上依照理论ε值计算了平均速度分布。在粘性作用层理论值与直接数值模拟符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。     

二维高超声速后台阶表面传热特性实验研究

高超声速后台阶流动是大气层内高速飞行器发动机设计、表面热防护以及高超声速拦截器红外成像窗口气动光学效应校正等诸多先进高超声速技术研发过程中所涉及的一类基础流动问题. 研究高超声速后台阶流动特性对有效提升飞行器综合性能, 进一步掌握高超声速流动机理具有重大基础 意义. 本文以二维高超声速后台阶流动为研究对象, 在KD-01高超声速激波风洞中测量了二维后台阶模型表面传热系数和表面静压, 并将实测台阶下游表面传热系数分布同采用高超声速边界层理论所得估计值进行了比较. 为进一步验证实验结果, 使用NPLS技术测量了其中一种实验状态下台阶周围流动结构. 研究发现, 对于二维高超声速后台阶流动, 台阶下游表面传热分布受台阶处边界层外缘流动特性的直接影响; 在台阶下游分离区和再附区内, 气体黏性占主导作用; 在台阶下游远场区域, 边界层流动特性趋同于平板边界层; 下游边界层基本结构取决于台阶处边界层相对厚度. 对高超声速后台阶流动, 若使用数值模拟方法研究气动热问题, 应当使用湍流模型.     

An explicit algebraic model of the planetary boundary layer turbulence: test computation of the neutrally stratified atmospheric boundary layer

An explicit algebraic model of Reynolds stresses and the turbulent heat flux vector for the planetary boundary layer in a neutrally stratified boundary layer of the atmosphere above a homogeneous rough surface is tested. The version of the algebraic model under consideration is constructed on the physical principles of the RANS (Reynolds-averaged Navier?Stokes) approximation for describing stratified turbulence, it employs three forecasting equations, and a correct reproduction of the main characteristics of a neutral atmospheric boundary layer — the components of the mean wind velocity, the wind turn angle, and the turbulent statistics is shown. Test computations show that the proposed model may be used for goal-oriented investigations of the atmospheric boundary layer.     

火星进入器小攻角飞行的静不稳定性

针对火星探路者号进入段飞行过程,求解三维流体力学Navier-Stokes方程,采用高温真实气体模型和等效比热比完全气体模型,对气动力特性进行预测和分析,对发现的小攻角静不稳定现象开展机理分析.比较海盗号数值计算结果、文献数据和飞行试验数据,结果符合很好,验证了物理化学模型及数值方法;探路者号的升阻特性和配平特性预测数据表明,真实气体模型和LAURA数据符合很好,等效比热比模型可以较准确地计算升力和阻力系数;2°攻角飞行时,进入器沿轨道出现静不稳定,完全气体模型无法模拟,分析认为迎风区声速线在肩部的移动和背风区亚声速区的分布变化以及亚声速区泡的出现,导致激波层压力下降和肩部膨胀区的压力下降向上游的传递过程发生变化,造成表面压力分布改变,最终诱发静不稳定.     

吸气式高超声速飞行器气动力气动热的数值模拟方法及应用

对吸气式高超声速飞行器而言,物面热流和摩阻的准确预测对飞行器设计及安全十分关键.介绍采用CFD准确预测气动力和气动热的方法,包括流动的控制方程、湍流模型及湍流的先进壁面函数边界条件,介绍流动的数值求解方法.对典型超声速层流和湍流流动的摩擦阻力和热流进行详细的验证与确认,考察CFD工具在使用先进壁面函数边界条件后,湍流计算的法向网格无关性能力.对设计的一种吸气式高超声速飞行器的气动力和气动热进行数值模拟,为飞行器的气动设计及热防护提供了可靠的数据.     

近壁区湍流脉动压力相关项的理论模型

提出在湍流边界层近壁区采用三维波的理论模型描述湍流相干结构,根据理论模型计算了Ｒｅｙｎｏｌｄｓ应力输运方程中的脉动速度与脉动压力梯度的相关项,理论计算结果与直接数值模拟（ＤＮＳ）符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。     

<Emphasis Type="Italic">β</Emphasis>-distribution for Reynolds stress and turbulent heat flux in relaxation turbulent boundary layer of compression ramp

A challenge in the study of turbulent boundary layers (TBLs) is to understand the non-equilibrium relaxation process after sep-aration and reattachment due to shock-wave/boundary-layer interaction. The classical boundary layer theory cannot deal with the strong adverse pressure gradient, and hence, the computational modeling of this process remains inaccurate. Here, we report the direct numerical simulation results of the relaxation TBL behind a compression ramp, which reveal the presence of intense large-scale eddies, with significantly enhanced Reynolds stress and turbulent heat flux. A crucial finding is that the wall-normal profiles of the excess Reynolds stress and turbulent heat flux obey a β-distribution, which is a product of two power laws with respect to the wall-normal distances from the wall and from the boundary layer edge. In addition, the streamwise decays of the excess Reynolds stress and turbulent heat flux also exhibit power laws with respect to the streamwise distance from the corner of the compression ramp. These results suggest that the relaxation TBL obeys the dilation symmetry, which is a specific form of self-organization in this complex non-equilibrium flow. The β-distribution yields important hints for the development of a turbulence model.     

Radiative-Convective Heat Transfer in Discharge Chamber at Turbulent Gas Flow

In the paper radiative-convective heat transfer in electric arc flowed through by Ar in the discharge chamber is studied. The equation system for heat, mass and momentum transfer in boundary layer approximation is considered, and the turbulence model with radiation pulsation and electrodynamic values is suggested. By solution of the equation system with the help of finite differences method distributions of temperature, heat flows and turbulent transfer coefficients over chamber length were obtained. The radiation and turbulence influence on the electric arc characteristics is analysed.     

基于温敏漆技术的圆锥高超声速大攻角绕流背风面流动结构实验研究

高超声速流动中, 大攻角下圆锥背风面边界层会存在流动分离与再附、边界层转捩等多种流动现象, 进而对圆锥表面温度分布产生显著的影响。为了对这一复杂流动规律及其对表面温升分布的影响进行讨论, 研究基于温敏漆技术, 得到了在Mach数为6的低湍流度来流条件下, 攻角为10°的圆锥背风面温升分布结果。通过对不同位置、不同方位角处温升分布曲线的分析, 对大攻角下圆锥背风面边界层流动发展过程及不同发展阶段的流动特征进行了讨论。同时, 通过对来流总压的调节, 得到了不同Reynolds数下的圆锥背风面温升分布结果, 总结了Reynolds数对流动的影响规律。研究发现, 高超声速大攻角圆锥背风面边界层流动发展过程中会依次出现层流分离、定常横流涡影响、转捩以及湍流分离与再附等流动特征, 而在不同的Reynolds数下, 各个流动特征产生影响的范围不同, 随着Reynolds数的降低, 层流范围和定常横流涡影响范围均有所增加, 而从观察到横流影响到转捩开始发生的范围基本相同。      

Interaction between coherent structures and surface temperature and its effect on ground heat flux in an unstably stratified boundary layer

Surface layer plumes, thermals, downdrafts and roll vortices are the most prominent coherent structures in an unstably stratified boundary layer. They contribute most of the temperature and vertical velocity variance, and their time scales increase with height. The effects of these multi-scale structures (surface layer plumes scale with surface layer depth, thermals scale with boundary layer height and the resulting roll vortices scale with convective time scale) on the surface temperature and ground heat flux were studied using turbulence measurements throughout the atmospheric boundary layer and the surface temperature measurements from an infrared camera. Plumes and thermals imprint on the surface temperature as warm structures and downdrafts imprint as cold structures. The air temperature trace shows a ramp-like pattern, with small ramps overlaid on a large ramp very close to the surface; on the other hand, surface temperature gradually increases and decreases. Turbulent heat flux and ground heat flux show similar patterns, with the former lagging the latter. The maximum values of turbulent heat flux and ground heat flux are 4 and 1.2 times the respective mean values during the ejection event. Surface temperature fluctuations follow a similar power-law exponent relationship with stability as suggested by surface layer similarity theory.     

Direct numerical simulation methods of hypersonic flat-plate boundary layer in thermally perfect gas

High-temperature effects alter the physical and transport properties of air such as vibrational excitation in a thermally perfect gas,and this factor should be considered in order to compute the flow field correctly.Herein,for the thermally perfect gas,a simple method of direct numerical simulation on flat-plat boundary layer is put forward,using the equivalent specific heat ratio instead of constant specific heat ratio in the N-S equations and flux splitting form of a calorically perfect gas.The results calculated by the new method are consistent with that by solving the N-S equations of a thermally perfect gas directly.The mean flow has the similarity,and consistent to the corresponding Blasius solution,which confirms that satisfactory results can be obtained basing on the Blasius solution as the mean flow directly in stability analysis.The amplitude growth curve of small disturbance is introduced at the inlet by using direct numerical simulation,which is consistent with that obtained by linear stability theory.It verified that the equation established and the simulation method is correct.     

γ-Re_θ转捩模型在内部冷却叶片中的应用

燃气轮机内部通流冷却叶片的流动传热特性对于燃机的整体效率有着重要的影响,而叶片表面的转捩流场一直是阻碍对叶片冷却性能进行精确分析的主要因素.本文采用SST湍流模型及γ-Re_θ转捩模型对MARK Ⅱ叶型进行数值模拟,通过分析叶栅通道内传热特性以考察该转捩模型的适应性,肯定了该转捩模型对于湍流边界层发生位置判断的正确性,以及对转捩发生点上游叶片表面温度变化趋势的预测的正确性.同时本文指出γ-Re_θ转捩模型与湍流模型结合时可能存在对湍动能过渡抑制的问题,有待进一步研究.     

基于Ludwieg管的高超声速边界层转捩实验

高超声速边界层层/湍流转捩是高超声速飞行器气动力和气动热设计中的难点和热点问题.为了降低开展高超声速边界层不稳定性与转捩实验研究的门槛,研究基于Ludwieg管原理设计并建造了一座Mach 6高超声速管风洞,重点对Ludwieg管风洞的启动和运行过程开展了数值模拟,分析了储气段弯管布局对试验段流场的影响;之后,对该高超...     

Velocity,temperature, and Reynolds-stress scaling in the wall region of turbulent boundary layer on a permeable surface

A finite total number of flow parameters in the wall region of a turbulent boundary layer points to universal behavior of turbulent shear stress as a function of mean-velocity gradient and turbulent heat flux as a function of both mean-velocity and mean-temperature gradients. Combined with dimensional arguments, this fact is used to reduce the momentum and heat equations to first-order ordinary differential equations for temperature and velocity profiles amenable to general analysis. Scaling laws for velocity and temperature in boundary layer flows with transpiration are obtained as generalizations of well-known logarithmic laws. Scaling relations are also established for shear stress and rms transverse velocity fluctuation. The proposed method has substantial advantages as compared to the classical approach (which does not rely on fluid-dynamics equations [1–3]). It can be applied to establish scaling laws for a broader class of near-wall turbulence problems without invoking closure hypotheses.     

Reply to the note on the flow behaviour of combustion gas in the boundary layer of an MHD generator

The temperature profile in Fig. 3 of Ref. 1 is better approximated by a relation for a turbulent boundary layer than a laminar boundary layer. The critical Reynolds number is expected to be lowered in the flow train of MHD combustion gas because of large turbulence in the combustor and by roughness along the duct.     

再入尾迹湍流对雷达散射截面影响分析

通过探讨高超声速再入湍流尾迹等离子体场中电磁波的散射机制,推导出在工程上描述湍流亚密等离子体雷达散射的一阶畸变波Born近似模型,分析了该模型在充分发展湍流尾迹等离子体场中的适用性,完成了适用于三维尾迹等离子体场的程序设计.以已有的湍流尾迹等离子体流场数据为基础,分析了再入尾迹湍流等离子体流动对雷达散射截面的影响.选取考察的几个有代表性的因素为:湍流模型、转捩过程、湍流尺度、电子组分脉动初始条件等.由结果可以看到,湍流转捩过程和湍流尺度对雷达散射截面值影响不大;电子组份脉动强度初始值影响较明显;湍流模型在特定条件下影响亦不大.     

LES of spatially developing turbulent boundary layer over a concave surface

We revisit the problem of a spatially developing turbulent boundary layer over a concave surface. Unlike previous investigations, we simulate the combined effects of streamline curvature as well as curvature-induced pressure gradients on the turbulence. Our focus is on investigating the response of the turbulent boundary layer to the sudden onset of curvature and the destabilising influence of concave surface in the presence of pressure gradients. This is of interest for evaluating the turbulence closure models. At the beginning of the curve, the momentum thickness Reynolds number is 1520 and the ratio of the boundary layer thickness to the radius of curvature is δ₀/R = 0.055. The radial profiles of the mean velocity and turbulence statistics at different locations along the concave surface are presented. Our recently proposed curvature-corrected Reynolds Averaged Navier-Stokes (RANS) model is assessed in an a posteriori sense and the improvements obtained over the base model are reported. From the large Eddy simulation (LES) results, it was found that the maximum influence of concave curvature is on the wall-normal component of the Reynolds stress. The budgets of wall-normal Reynolds stress also confirmed this observation. At the onset of curvature, the effect of adverse pressure gradient is found to be predominant. This decreases the skin friction levels below that in the flat section.