扰流板作用下湍流/非湍流界面特性

在零压梯度平板湍流边界层流动中安装垂直流向高度为h的扰流板,诱导流场产生横向涡,研究横向涡影响下湍流/非湍流界面特性沿流向的发展。结果表明,在本实验条件下,整个流动经历了从湍流边界层到流动分离和再附,再向湍流边界层恢复的过程。在扰流板下游约18h距离后,扰流板尾迹的影响逐渐衰退,壁面剪切对湍流强度的贡献开始逐渐恢复,在扰流板下游约55h距离后,湍流边界层再次充分发展。与此同时,由于扰流板后流场流动结构拟序性的增强,湍流/非湍流界面的分形维度受扰流板影响而减小,表明脱落涡有使界面多尺度特性、三维性减弱的趋势。此外,界面高度的概率密度分布受扰流板影响呈现显著的右偏,主要与扰流板增强喷射运动强度,使得界面更容易抬升相关。流动结构及界面特性受扰流板影响后的流向演化有同步变化的模式,扰流板对界面特性影响主要集中于（-5~18）h的流向范围。      

各向异性涡黏模型在通道分离流动中的应用

本文针对燃机带肋冷却通道发展了一种考虑近壁特性和远壁流动各向异性特点的代数各向异性涡黏模型。本文针对肋片处分离流动的各向异性特点,结合隐式代数应力模型和各向同性涡黏模型,推导各向异性因子修正三个正向涡黏系数,反映正雷诺应力的各向异性特点。同时,通过近壁分析完善k和ε封闭形式。将各向异性涡黏模型应用到,1)后台阶流动预测中,湍流流动预测精度提高约10%;2)带90°肋通道流动传热预测中,平均流动预测精度得到改善,传热预测精度提高6%~10%。结果表明,本文针对内冷通道的分离流动作出的各向异性修正是有效的。     

不同亚格子模式在后台阶湍流流动大涡模拟中的应用

本文用大涡模拟方法研究了湍流后台阶流场中的大涡演变过程,并在此基础上研究了目前大涡模拟中比较常用的六种亚格子模式。在相同的流动几何参数以及计算条件下,给出了不同亚格子模式下湍流流动瞬时压力场以及流场中瞬时粘性大小的分布,尤其给出了流场瞬时演变的大尺度涡结构。将六种亚格子模式从回流区长度、计算时间、计算结果的准确性和流场稳定性等不同角度进行了比较。为选取合适的亚格子模式深入研究湍流耗散机理奠定了基础。     

带喷流超声速后台阶流场精细结构及其运动特性研究

采用基于纳米示踪的平面激光散射技术(NPLS)对带超声速喷流的后台阶流动精细结构进行了研究. 来流马赫数为3.4, 喷流实测马赫数为2.45, 而名义马赫数为2.5. 结果清晰地揭示了激波、剪切层、混合层、Kelvin-Helmholtz涡、羊角涡及湍流大尺度结构等大量典型流场结构. 基于大量流场精细结构图像, 对典型位置处的结构进行了空间两点相关性分析, 在喷流混合层前端涡结构小于湍流充分发展的尾端, 结构角相对也小. 喷流工作时, 模型台阶下游表面由一薄层气膜覆盖. 获得了模型流向和不同高度展向平面内的流场结构, 对照纹影试验结果, 分析了流动特点及时间演化规律. 采用微型压力扫描系统测试了模型表面的压力系数分布, 靠近喷流下游处压力系数区域0.0146. 针对NPLS图像做了流动的分形维数的分析, 发现在流动初始阶段分形维数接近于1, 越靠下游分形维数越高.     

湍流混合层流场的PIV测量

本文使用PIV对在坚直通道内放置一个特殊设计的隔板所形成的湍流混合层流动进行测量,高低侧速比为4:1,基于两股流体速度差和管道水力半径的Re数范围4400～158400.发现混合层中大涡拟序结构的尺度随雷诺数的增加而增大,而后又随雷诺数的继续增大而减小,气泡的加入会延缓或阻碍大涡拟序结构的发展.对雷诺应力、湍流强度、涡量、旋涡强度在混合层流场内随雷诺数的变化和分布规律进行分析,发现混合层内雷诺应力、湍流强度、涡量及旋涡强度均集中分布在隔板下游一个较窄的锥形区域内,雷诺应力和湍流强度随雷诺数的增大先增人后减小,随离开隔板距离的增大而减小.涡量及旋涡强度随雷诺数的增大而增加,随离开隔板距离的增大而减小.     

燃气透平叶栅端壁传热特性的数值研究

采用三维数值求解方法,对透平叶栅端壁区域的流动和传热特性进行了研究.利用试验数据考核了相应的数值方法,分析了网格数目和湍流模型对叶栅端壁附近流动传热特性计算的影响,比较了不同进口雷诺数和湍流度条件下端壁传热特性的变化。结果表明;马蹄涡和通道涡等二次流动直接影响端壁区传热,传热强度分布规律基本与当地流动的湍动能保持一致。湍流模型对端壁压力场的计算影响较小,但对端壁传热特性的求解的精度影响较大。采用v~3—f湍流模型能较好地预测端壁传热分布。来流雷诺数和湍流度的变化改变了端壁边界层厚度和涡系结构,使得端壁传热强度和梯度分布发生变化。     

带控制的超声速二维后台阶流动显示

在Mach数3.4的来流条件下,对二维后台阶流动精细结构开展了实验研究.实验分为后台阶上游无控制加粗糙带扰动及微涡流发生器（micro-vortex generator,MVG）扰动3种状态,采用基于纳米示踪的平面激光散射（nano-tracer based planar laser scattering,NPLS）方法获得了流向和展向切面内的高时空分辨率流动显示图像,并测量了模型表面静压分布.对大量NPLS图像取平均,研究了流场结构的时间平均规律,对比不同时刻的瞬态流场精细结构图像,发现不同状态下的湍流大尺度结构的特征时间.有粗糙带状态相对无粗糙带台阶下游回流区压力更低,而下游压力较高,台阶上游区别不大;受MVG控制后台阶下游附近区域压力突增;MVG对流动的控制改变能力较强,粗糙带能调整台阶上下游附近流动平稳过渡,流场壁面压力没有突变.      

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

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

低速二维后向台阶流动中的壁面压强脉动

为了能排除三维结构对涡脱和剪切层旋涡的影响，使用实验方法研究了一个小展高比（AR=0.125）的后向台阶流动.该实验台类似纯二维的Hele-Shaw Cell.流动被局限在两个平行且距离为5 mm的有机玻璃板之间.台阶高度H为40 mm，扩张比2:3.在台阶下游中央沿流动方向安装16个麦克风组成的传感器阵列采集壁面脉动压强.来流速度U₀在9~26 m/s之间连续可调.通过计算脉动压强分布、频谱，不同位置的相关性和相干性系数，发现并分析流动存在一个临界Reynolds数.流场在临界Reynolds数前后存在明显不同的流动特征.实验结果表明在低Reynolds数下依然存在剪切层的低频摆动；当Reynolds数大于临界Reynolds数时，分离后流动由涡结构传播的特性主导.      

DES模型在压气机亚音转子中的应用探讨

本文运用分离涡模拟(DES)方法研究了不同工况下压气机亚音转子的流动情况,分析了其时均与瞬时流场中顶部间隙泄漏流动和根部角区的流动分离.通过与S-A模型计算结果的对比表明,DES模型在模拟顶部泄漏流动及二次泄漏、泄漏流在转子下游与尾迹的干涉时能够捕捉到更强的旋涡结构,在模拟转子根部角区的分离时也能获得更为丰富的流动现象.对不同工况的DES计算表明负荷的上升会使泄漏涡形成的位置向上游移动,从而导致并加剧二次流动,并对叶栅下游泄漏涡与尾迹的干涉产生影响.对设计工况下瞬时流场的分析表明,泄漏涡在叶栅下游体现出周期性的强弱变化,近叶根分离区也体现出明显的非定常性.     

An algebraic variational multiscale–multigrid method for large-eddy simulation of turbulent variable-density flow at low Mach number

An algebraic variational multiscale–multigrid method is proposed for large-eddy simulation of turbulent variable-density flow at low Mach number. Scale-separating operators generated by level-transfer operators from plain aggregation algebraic multigrid methods enable the application of modeling terms to selected scale groups (here, the smaller of the resolved scales) in a purely algebraic way. Thus, for scale separation, no additional discretization besides the basic one is required, in contrast to earlier approaches based on geometric multigrid methods. The proposed method is thoroughly validated via three numerical test cases of increasing complexity: a Rayleigh–Taylor instability, turbulent channel flow with a heated and a cooled wall, and turbulent flow past a backward-facing step with heating. Results obtained with the algebraic variational multiscale–multigrid method are compared to results obtained with residual-based variational multiscale methods as well as reference results from direct numerical simulation, experiments and LES published elsewhere. Particularly, mean and various second-order velocity and temperature results obtained for turbulent channel flow with a heated and a cooled wall indicate the higher prediction quality achievable when adding a small-scale subgrid-viscosity term within the algebraic multigrid framework instead of residual-based terms accounting for the subgrid-scale part of the non-linear convective term.     

Investigations on the Effect of Backward-Facing and Forward-Facing Steps on Turbulent Mixed-Convection Flow Over a Flat Plate

The effects of backward-facing and forward-facing steps on a turbulent buoyancy-dominated mixed-convection flow over a flat plate are examined experimentally. Air velocity and temperature distributions and their turbulent fluctuations are measured simultaneously by using a two-component laser-Doppler velocimeter and a cold wire anemometer, respectively. The experiment was carried out for a step (backward-facing/forward-facing) height of 22 mm, a temperature difference, ΔT, of 30°C between the heated walls and the free-stream air (corresponding to a local Grashof number Gr _xi = 6.45 × 10¹⁰), and a free-stream velocity of 0.48 m/s. It was found that the introduction of backward- and forward-facing steps increases the turbulence intensity of the velocity and temperature fluctuations downstream of the step. The present results also reveal that the maximum local Nusselt number occurs in the vicinity of the reattachment zone, and it is approximately twice for the case of the backward-facing step and three times for the case of the forward-facing step than that of the flat plate value at similar flow and thermal conditions.     

On flow unsteadiness induced by structural vibration

In this paper, the effects of structural vibration on flow unsteadiness are investigated numerically. A fully coupled model, that solves the unsteady flow equations as well as the dynamic equations of the structure, is used. Numerical experiments are carried-out for flow over a backward-facing step, where a large number of numerical and experimental data exist for comparisons. The flexible structure is upstream of the step and is excited by a plane acoustic wave from the side opposite to the flow. Three Reynolds number cases are studied: 300 for a laminar flow, 3000 for a transitional flow and 15 000 for a turbulent flow. The results obtained are in good agreement with experimental observations and show the strong coupling between structural vibration and the resulting flow unsteadiness.     

TURBULENT NATURAL CONVECTION FLOW OVER A BACKWARD-FACING STEP

Measurements of turbulent natural convection boundary-layer air flow over a two-dimensional, vertical, backward-facing step are reported. The step geometry consists of an adiabatic backward-facing step, an upstream wall, and a downstream wall. Both the upstream and downstream walls are heated to a uniform and constant temperature. Laser-Doppler velocimeter and cold-wire anemometer were used, respectively, to measure simultaneously the time-mean velocity and temperature distributions and their turbulent fluctuations. The experiment was carried out for step heights of 0, 11, and 22 mm and a temperature difference, Delta T, of 30 degrees C between the heated walls and the ambient air. The present results reveal that the turbulence intensity of both the streamwise and transverse velocity fluctuations, and the intensity of temperature fluctuations downstream of the step, increase as the step height increases. Also, it was found that both the reattachment length and the heat transfer rate from the downstream heated wall increase with increasing step height.     

Flame-acoustic coupling of combustion instability in a non-premixed backward-facing step combustor: the role of acoustic-Reynolds stress

Combustion instability in a laboratory scale backward-facing step combustor is numerically investigated by carrying out an acoustically coupled incompressible large eddy simulation of turbulent reacting flow for various Reynolds numbers with fuel injection at the step. The problem is mathematically formulated as a decomposition of the full compressible Navier–Stokes equations using multi-scale analysis by recognising the small length scale and large time scale of the flow field relative to a longitudinal mode acoustic field for low mean Mach numbers. The equations are decomposed into those for an incompressible flow with temperature-dependent density to zeroth order and linearised Euler equations for acoustics as a first order compressibility correction. Explicit coupling terms between the two equation sets are identified to be the flow dilatation as a source of acoustic energy and the acoustic Reynolds stress (ARS) as a source of flow momentum. The numerical simulations are able to capture the experimentally observed flow–acoustic lock-on that signifies the onset of combustion instability, marked by a shift in the dominant frequency from an acoustic to a hydrodynamic mode and accompanied by a nonlinear variation of pressure amplitude. Attention is devoted to flow conditions at two Reynolds numbers before and after lock-on to show that, after lock-on, the ARS causes large-scale vortical rollup resulting in the evolution of a compact flame. As compared to acoustically uncoupled simulations at these Reynolds numbers that show an elongated flame with no significant roll up and disturbance in the upstream flow field, the ARS is seen to alter the shear layer dynamics by affecting the flow field upstream of the step as well, when acoustically coupled.     

A one-equation turbulence model for recirculating flows

A one-equation turbulence model which relies on the turbulent kinetic energy transport equation has been developed to predict the flow properties of the recirculating flows. The turbulent eddy-viscosity coefficient is computed from a recalibrated Bradshaw’s assumption that the constant a₁ = 0.31 is recalibrated to a function based on a set of direct numerical simulation (DNS) data. The values of dissipation of turbulent kinetic energy consist of the near-wall part and isotropic part, and the isotropic part involves the von Karman length scale as the turbulent length scale. The performance of the new model is evaluated by the results from DNS for fully developed turbulence channel flow with a wide range of Reynolds numbers. However, the computed result of the recirculating flow at the separated bubble of NACA4412 demonstrates that an increase is needed on the turbulent dissipation, and this leads to an advanced tuning on the self-adjusted function. The improved model predicts better results in both the non-equilibrium and equilibrium flows, e.g. channel flows, backward-facing step flow and hump in a channel.     

两相湍流统一色噪声法概率密度函数模型

统一色噪声近似方法对简单一维色噪声问题研究较为充分, 本文将统一色噪声法应用到高度复杂的多维气固两相湍流系统之中.首先从颗粒运动Langevin方程出发, 利用统一色噪声法获得两相湍流Fokker-Planck方程, 然后以此为基础建立颗粒轨道两阶矩模型.文中建立的新模型成功应用于后台阶两相湍流流场的数值模拟, 预报合理正确.研究表明, 对于多维两相湍流系统, 统一色噪声法仍然行之有效.     

Investigation on the flameholding mechanisms in supersonic flows: backward-facing step and cavity flameholder

Abstract  

As effective devices to extend the fuel residence time in supersonic flow and prolong the duration time for hypersonic vehicles cruising in the near-space with power, the backward-facing step and the cavity are widely employed in hypersonic airbreathing propulsive systems as flameholders. The two-dimensional coupled implicit RANS equations, the standard k-ε turbulence model, and the finite-rate/eddy-dissipation reaction model have been used to generate the flow field structures in the scramjet combustors with the backward-facing step and the cavity flameholders. The flameholding mechanism in the combustor has been investigated by comparing the flow field in the corner region of the backward-facing step with that around the cavity flameholder. The obtained results show that the numerical simulation results are in good agreement with the experimental data, and the different grid scales make only a slight difference to the numerical results. The vortices formed in the corner region of the backward-facing step, in the cavity and upstream of the fuel injector make a large difference to the enhancement of the mixing between the fuel and the free airstream, and they can prolong the residence time of the mixture and improve the combustion efficiency in the supersonic flow. The size of the recirculation zone in the scramjet combustor partially depends on the distance between the injection and the leading edge of the cavity. Further, the shock waves in the scramjet combustor with the cavity flameholder are much stronger than those that occur in the scramjet combustor with the backward-facing step, and this causes a large increase in the static pressure along the walls of the combustor.     

Assessment of four inflow conditions on large-eddy simulation of a gently curved backward-facing step

In this study, we investigated the impact of four different inflow generation methods used in large-eddy simulation on the spatially developing boundary-layer before a recirculation bubble is formed over a gently curved backward-facing step at a Reynolds number 13,700 based on the step height. The configuration under study is very challenging because the separation is caused by a weak adverse pressure gradient, thus making a very sensitive and reliable assessment for evaluating the different inflow conditions used. The first method is a precursor-based simulation in which the velocity data on a certain plane is recorded in a library and then used as the inlet condition for the primary simulation. The second method used is the random noise generation method. The third method is based on generation of turbulent spots which incorporates the distribution of Reynolds stress tensor, and the last one is the so-called rescaling/recycling method proposed by Lund and colleagues. All these methods are compared together in terms of separation and reattachment locations of the recirculation bubble. The flow structures are represented by qualitative criteria, and also streamwise Reynolds stresses and production of turbulent kinetic energy of the flow are assessed and compared together at different stations before and after separation to illuminate how the developing structures within the boundary-layer can affect the locations of separation and reattachment. Distribution of pressure coefficient for different methods showed that there is a relation between production of turbulent kinetic energy and favourable pressure gradient of each method before the separation occurs. Finally, spectra of pressure fluctuation revealed how each inflow condition influences the shedding-instability frequencies.