Transition region where the large-scale and very large scale motions coexist in atmospheric surface layer: wind tunnel investigation

The turbulent structures in atmospheric surface layer (ASL) are investigated in wind tunnel with hot-wire anemometers in present study. The results show that there exist two typical turbulent structures characterised by their streamwise length scales, i.e. large-scale motions (LSMs) and very large scale motions (VLSMs) as reported recently in pipe flow, channel flow, zero-pressure-gradient turbulent boundary layers and near-neutral ASLs. Moreover, a transition region containing both LSMs and VLSMs is first reported in present investigation, and this region depends on the Reynolds numbers. In the transition region, as the wall-normal distance increases, the turbulent energy ratio of LSMs is gradually weakened but VSLMs strengthened.     

压缩拐角激波/湍流边界层干扰特性分析

用大涡模拟方法对Mach数3.0下的压缩拐角激波/湍流边界层干扰问题进行数值研究.对拐角上游平板区域边界层转捩及湍流进行模拟,设定平板区域长度,使得转捩过程于平板区域发生并充分完成,从而在拐角处产生激波/湍流边界层相互干扰,研究激波/湍流边界层的作用机理.研究表明:流场能够在非定常扰动激励下迅速转捩,并于平板区发展为完全湍流;湍流边界层与激波相互作用过程中,拐角附近分离区较层流情况明显减小;展向不同区域分离区大小差异较大,局部区域分离现象消失.     

Controlled detonation initiation in hypersonic flow

Experimental evidence of controlled detonation initiation and propagation in a hypersonic flow of premixed hydrogen-air is presented. This controlled detonation initiation is created in a hypersonic facility capable of producing a Mach 5 flow of hydrogen-air. Flow diagnostics such as high-speed schlieren and OH* chemiluminescence results show that a flame deflagration-to-detonation transition occurs as a combined result of turbulent flame acceleration and shock-focusing. The experimental results define three new distinct regimes in a Mach 5 premixed flow: deflagration-to-detonation transition (DDT), unsteady compressible turbulent flames, and shock-induced combustion. A two-dimensional implicit-LES (ILES) simulation, which solves the compressible, reactive Navier-Stokes equations on an adapting grid is conducted to provide additional insight into the local physical mechanism of detonation transition and propagation.     

Large Eddy Simulation of a Vortex Ring Impacting a Bump

A vortex ring impacting a three-dimensional bump is studied using large eddy simulation for a Reynolds number Re=$4$x$10^4$ based on the initial diameter and translational speed of the vortex ring. The effects of bump height and vortex core thickness for thin and thick vortex rings on the vortical flow phenomena and the underlying physical mechanisms are investigated. Based on the analysis of the evolution of vortical structures, two typical kinds of vortical structures, i.e., the wrapping vortices and the hair-pin vortices, are identified and play an important role in the flow state evolution. The boundary vorticity flux is analyzed to reveal the mechanism of the vorticity generation on the bump surface. The circulation of the primary vortex ring reasonably elucidates some typical phases of flow evolution. Further, the analysis of turbulent kinetic energy reveals the transition from laminar to turbulent state. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the flow evolution and the flow transition to turbulent state.     

转捩与湍流对激波边界层干扰及底部流动结构的影响

为研究转捩与湍流对激波边界层干扰及底部流动结构的影响,文章选取了二维与三维高超声速双斜面进气道模型与大钝头着陆器模型,并使用γ-Re_θ转捩模型开展数值模拟研究.研究表明,对于二维进气道模型,随着前缘钝度的增加,激波边界层干扰位置前移,分离区变大,与层流流动情况相比,有转捩流动发生时,激波边界层干扰位置后移,同时分离流动强度变弱,分离区缩小;对于三维进气道模型,其拐角附近的分离流动呈现明显的三维特征,转捩流动也存在三维流动结构,与静风洞状态相比,噪音风洞状态下,有转捩流动发生,对壁面热流影响较大,对激波系影响很小.对于着陆器模型,底部流动发生转捩,使得底部流动由不稳定非定常的流动结构变为稳定定常的流动结构,这有益于姿态控制设计.      

不可压湍流的控制方程

使用侧偏统计平均方法理性推导了不可压湍流的封闭方程组．侧偏统计平均提供了湍流脉动量的一阶统计信息，在引入加权漂移速度对称性及正交各向异性后，使用动量传输链概念模化封闭了整个方程组．方程组不含任何经验系数，不使用壁面函数，保留了ＮＳ方程的均化的非线性特性．其级数形式的能量方程与非线性现象多尺度层次现象相对应，具备了描述湍流统计平均流动及拟序结构流动的双重功能．用平面射流，圆射流，层流－湍流转捩及后向台阶流等算例初步验证了方程对真实湍流的适定性．     

Investigation of Turbulent Transition in Plane Couette Flows Using Energy Gradient Method

The energy gradient method has been proposed with the aim of better understanding the mechanism of flow transition from laminar flow to turbulent flow. In this method, it is demonstrated that the transition to turbulence depends on the relative magnitudes of the transverse gradient of the total mechanical energy which amplifies the disturbance and the energy loss from viscous friction which damps the disturbance, for given imposed disturbance. For a given flow geometry and fluid properties, when the maximum of the function $K$ (a function standing for the ratio of the gradient of total mechanical energy in the transverse direction to the rate of energy loss due to viscous friction in the streamwise direction) in the flow field is larger than a certain critical value, it is expected that instability would occur for some initial disturbances. In this paper, using the energy gradient analysis, the equation for calculating the energy gradient function $K$ for plane Couette flow is derived. The result indicates that $K$ reaches the maximum at the moving walls. Thus, the fluid layer near the moving wall is the most dangerous position to generate initial oscillation at sufficient high $\operatorname{Re}$ for given same level of normalized perturbation in the domain. The critical value of $K$ at turbulent transition, which is observed from experiments, is about 370 for plane Couette flow when two walls move in opposite directions (anti-symmetry). This value is about the same as that for plane Poiseuille flow and pipe Poiseuille flow (385-389). Therefore, it is concluded that the critical value of $K$ at turbulent transition is about 370-389 for wall-bounded parallel shear flows which include both pressure (symmetrical case) and shear driven flows (anti-symmetrical case).     

二维槽道湍流拟序结构的大涡模拟

本文采用大涡模拟的方法,对二维槽道湍流流动进行了数值模拟。采用Chorin的分步投影法求解大尺度涡运动的Navier-Stokes方程,小尺度涡采用三种亚格子（SGS）模式分别模拟,给出了不同亚格子涡粘性模式下的模拟结果。对固壁面采用了壁函数。模拟结果再现了二维槽道流动拟序结构的发展演变过程。通过对不同入口速度下的瞬态流场的比较,揭示了入口速度分布对流场的影响。     

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

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

Experimental verification of the early crisis of drag using a single sphere

The recently discovered phenomenon of drag anomaly (early crisis) observed for droplets in a two-phase turbulent flow at transition Reynolds numbers (Re = 10–100) is experimentally supported in the case of a gas flow about a single hard sphere as well. Using a specially designed torsion balance, the drag force of a foam plastic sphere 3 mm in diameter placed in a blower-produced turbulent jet is measured. It is found that the turbulence of the flow about the sphere must be high for the drag anomaly to arise.     

3维Mortar谱元法模拟Kelvin-Helmhtz不稳定性混合层

采用Mortar谱元法和多处理器并行计算技术模拟了Kelvin-Helmhtz界面不稳定性湍流的混合发展过程,通过对混合层动量厚度、能谱和总动能的计算,评估了Kelvin-Helmhtz混合层的演化机理。计算结果表明:3维Mortar谱元法具有高计算精度和光滑区域的指数收敛特性,可以有效模拟混合层流动的湍流混合和演化,能够捕捉到涡的合并现象和大涡到小涡的级联过程;初期的混合层层流运动发展成具有连续谱结构的湍流运动过程,实现了Kelvin-Helmhtz界面不稳定性混合层流动从2维发展到3维的转捩特征,总湍流统计动能的变化反映了粘性耗散过程的作用。通过对Kelvin-Helmhtz 3维界面不稳定性混合层流动和3维层流向湍流转捩过程的数值模拟,程序的有效性得到了验算,表明谱元法应用于湍流混合模拟是可行的。     

Turbulence in boundary flow of superfluid 4He triggered by free vortex rings

The transition to turbulence in the boundary flow of superfluid 4He is investigated using a vortex-free vibrating wire. At high wire vibration velocities, we found that stable alternating flow around the wire enters a turbulent phase triggered by free vortex rings. Numerical simulations of vortex dynamics demonstrate that vortex rings can attach to the surface of an oscillating obstacle and expand unstably due to the boundary flow of the superfluid, forming turbulence. Experimental investigations indicate that the turbulent phase continues even after stopping the injection of vortex rings, which is also confirmed by the simulations.     

An unsteady pseudoshock model for barotropic gas flow

A mathematical unsteady pseudoshock model describing the continuous transition from supersonic to subsonic flow is constructed for a barotropic gas flow in a long flat channel or a nozzle. The model is based on a two-layer scheme of flow with mass transfer including a potential supersonic core and a turbulent boundary layer.     

Raman phase conjugate resonator for intracavity aero-optic turbulence aberration correction

A cw degenerate phase conjugate resonator (PCR) with intracavity aero-optic turbulent flow is demonstrated using sodium vapor pumped with a low intensity cw laser. A turbulent helium jet at a high flow velocity of 250 m/s and forcing frequency of 18 kHz generated the spatio-temporal aberrations. The intracavity normalized peak-to-peak intensity fluctuations at 18 kHz measured with a pinhole-detector are ten times smaller in the corrected beam than in the aberrated beam. The threshold pump intensity for the resonator is 4.5 W/cm² in the absence of the turbulent flow. The pump frequency is tuned well within the Doppler width of the Na D₁ transition. Evidence of a degenerate (laboratory frame) double-Λ mechanism involving Zeeman sublevels is also shown.     

Friction factor and wall heat transfer for laminar and turbulentflow in a cylindrical duct with a wall-stabilized arc

Friction factor and wall heat-transfer data for axially symmetric flow in a wall-stabilized arc analysis are given. Heat transfer results revealed three significant modes: laminar flow with an unstable arc, laminar flow with a wall-stabilized arc, and turbulent flow with a wall-stabilized arc. It was shown that the critical Reynolds number increases, in comparison with the case of the flow without Joule heating. A friction factor from Reynolds-number dependence peculiarities was not discovered for laminar to turbulent flow transition     

Role of large scale flow features on cycle-to-cycle variations of spark-ignited flame-initiation and its transition to turbulent combustion

This study investigates the influence of large-scale flow features, including flow structure and velocity magnitude, on the early-burn period variability in a homogenous-charge spark-ignited engine fueled with premixed propane-air mixture. Particle image velocimetry and in-cylinder pressure measurement data from a previous study - were processed to enable simultaneous flow characterization and flame-front tracking as well as apparent heat-release analysis. By combining probability analysis of flame development with conditional sampling of fast and slow early-burn cycles using 10% fuel mass fraction burned, it is shown that an undesirable flow structure produces an asymmetric flame development at the initial flame growth period. This asymmetric flame structure persists through the whole initial-to-turbulent transition period until the flame becomes fully turbulent. The undesirable flow condition is characterized by large-scale convective flows near spark plug rather than flows that lead to increased flame spread in multiple directions. The simultaneous flow and flame characterization enables the quantifications of flame-front propagation speed, unburned fuel-air mixture velocity ahead of flame front and local burning velocity at flame surface. Here the local burning velocity is referred to as laminar or turbulent flame speed. A simplified approach is introduced to derive integrated values for these quantities per crank-angle-degree, enabling the quantitative comparison of the trend-wise difference in these integrated metrics between fast and slow early-burn cycles. It is revealed that for the transition period, the CCV in the velocity magnitude of unburned fuel-air mixture ahead of the flame front accounts for nearly 50% to the variability of flame propagation speed. The burning velocity provides the remaining source of the flame propagation variability in this period. The flame propagation variations in the initial flame growth and fully turbulent periods are smaller than those in the transition period and are primarily dependent on the variability of large-scale flow features.     

L-mode to H-mode transition triggered by sawtooth-induced heat flux in EAST

H-modes induced by sawtooth events can be often observed in discharges with marginal auxiliary power injection in EAST. Poloidal flow shear at the very plasma edge, increasing ∼25% up to the threshold value, is observed just before the L-H transition by means of a fast reciprocating probe array in EAST. This suddenly risen poloidal flow shear, caused by the increased turbulent driven Reynolds force, is motived by the heat pulse originally released by a sawtooth crash at the plasma core. Associated with the critical poloidal flow shear, the local turbulent decorrelation rate increases significantly. The increased turbulent decorrelation rate compensated by nonlinear energy transfer rate from the turbulence to the low-frequency shear flows, exceeding the turbulence energy input rate, is sustained for several hundred microseconds till the turbulence quench happening.     

How does flow in a pipe become turbulent?

The transition to turbulence in pipe flow does not follow the scenario familiar from Rayleigh-Benard or Taylor-Couette flow since the laminar profile is stable against infinitesimal perturbations for all Reynolds numbers. Moreover, even when the flow speed is high enough and the perturbation sufficiently strong such that turbulent flow is established, it can return to the laminar state without any indication of the imminent decay. In this parameter range, the lifetimes of perturbations show a sensitive dependence on initial conditions and an exponential distribution. The turbulence seems to be supported by threedimensional travelling waves which appear transiently in the flow field. The boundary between laminar and turbulent dynamics is formed by the stable manifold of an invariant chaotic state. We also discuss the relation between observations in short, periodically continued domains, and the dynamics in fully extended puffs.     

Ginzburg–Landau description of laminar-turbulent oblique band formation in transitional plane Couette flow

Plane Couette flow, the flow between two parallel planes moving in opposite directions, is an example of wall-bounded flow experiencing a transition to turbulence with an ordered coexistence of turbulent and laminar domains in some range of Reynolds numbers [R _g, R _t] . When the aspect-ratio is sufficiently large, this coexistence occurs in the form of alternately turbulent and laminar oblique bands. As R goes up trough the upper threshold R _t, the bands disappear progressively to leave room to a uniform regime of featureless turbulence. This continuous transition is studied here by means of under-resolved numerical simulations understood as a modelling approach adapted to the long time, large aspect-ratio limit. The state of the system is quantitatively characterised using standard observables (turbulent fraction and turbulence intensity inside the bands). A pair of complex order parameters is defined for the pattern which is further analysed within a standard Ginzburg–Landau formalism. Coefficients of the model turn out to be comparable to those experimentally determined for cylindrical Couette flow.