Characteristic Regimes of Transitional Separation Bubbles in Unsteady Flow

This experimental investigation deals with transition phenomena of a separated boundary layer under unsteady inlet flow conditions. The main purpose of this investigation is to understand the influence of the rotor-stator interaction in turbomachinery on the subsequent, highly loaded boundary layer. The research project is divided into two phases. In the first phase, which has been completed recently, only the variation of mean velocity caused by upstream blades was simulated in the experiments while the free-stream turbulence intensity was retained at a constant low level. The experiments are carried out in an Eifel-type wind tunnel to investigate the laminar separated boundary layer of a flat plate under oscillating inlet conditions. The adverse pressure gradient, similar to that of turbomachines, is generated by the contoured upper wall. The unsteadiness is produced by a rotating flap located downstream of the test section. The reduced frequency, the amplitude and the mean Reynolds number are varied to simulate the conditions prevailing in turbomachines. In addition to the Kelvin–Helmholtz instability of the separated shear layer, a lower frequency instability was observed. This is frequently referred to as `free shear layer flapping' and results in two distinctly different ways of re-attachment, depending primarily on the Reynolds number. For low momentum thickness Reynolds numbers at the separation point, large-scale vortices locked to the frequency of the unsteady main flow are identified. They originate nearly at the top of the separation bubble and are ejected downstream. A fully turbulent boundary layer develops after these vortices mix out. For higher Reynolds numbers, transition is completed within a short length of the free shear layer and there-attachment region. The characteristic momentum thickness Reynolds number separating these two regimes in unsteady flow is about 125. The Strouhal number (reduced frequency) does not appear to have any significant effect. Based on the experimental results, this behaviour is discussed in some detail. This revised version was published online in July 2006 with corrections to the Cover Date.     

Experimental study on the effect of unsteadiness on boundary layer development on a linear turbine cascade

 The results from an experimental investigation of unsteady boundary layer behavior on a linear turbine cascade are presented in this paper. To perform a detailed study on unsteady cascade aerodynamics and heat transfer, a new large-scale, high-subsonic research facility for simulating the periodic unsteady flow has been developed. It is capable of sequentially generating up to four different unsteady inlet flow conditions that lead to four different passing frequencies, wake structures, and freestream turbulence intensities. For a given Reynolds number, two different unsteady wake formations are utilized. Detailed unsteady boundary layer velocity. turbulence intensity, and pressure measurements are performed along the suction and pressure surfaces of one blade. The results display the transition and development of the boundary layer, ensemble-averaged velocity, and turbulence intensity. Received: 23 September 1996/Accepted: 19 February 1997     

The Effect of wake Turbulence Intensity on Transition in a Compressor Cascade

Direct numerical simulations of separating flow along a section at midspan of a low-pressure V103 compressor cascade with periodically incoming wakes were performed. By varying the strength of the wake, its influence on both boundary layer separation and bypass transition were examined. Due to the presence of small-scale three-dimensional fluctuations in the wakes, the flow along the pressure surface undergoes bypass transition. Only in the weak-wake case, the boundary layer reaches a nearly-separated state between impinging wakes. In all simulations, the flow along the suction surface was found to separate. In the simulation with the strong wakes, separation is intermittently suppressed as the periodically passing wakes managed to trigger turbulent spots upstream of the location of separation. As these turbulent spots convect downstream, they locally suppress separation.     

低压涡轮内部流动及其气动设计研究进展

随着高空无人飞行器研究的不断升温, 高空低雷诺数条件下动力装置的研究越来越受到人们的重视.结合近年来国内外相关领域的研究工作, 对低雷诺数低压涡轮内部复杂流动机理的研究进展进行了介绍, 包括低雷诺数情况下低压涡轮内部非定常流动的特点, 叶片边界层分离及转捩现象机理, 上游周期性尾迹与下游叶片边界层相互作用机理等. 在此基础上给出了适合低雷诺数条件的低压涡轮气动设计方法:尾迹通过与边界层的相互作用, 能够抑制分离, 进而减小叶型损失, 在气动设计中有效引入非定常效应可以大幅度提高低压涡轮的气动负荷或降低气动损失, 最终达到提高性能的目的;数值及实验结果验证了这种设计方法的有效性.      

Computational analysis of impinging shock-wave boundary layer interaction under conditions of incipient separation

The interaction of an oblique shock wave with a turbulent boundary layer under conditions of incipient separation is analyzed by means of large-eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) turbulence models, with the objective to explore their predictive capabilities, in particular with respect to the unsteady features of the interaction. Consistent with earlier direct numerical simulations, we have found that the flow dynamics in the interaction zone is characterized by strong intermittency associated with the formation of scattered spots of flow reversal near the nominal position of the reflected shock. Comparison with experimental results (at much larger Reynolds number) show that the qualitative features of the interaction are predicted reasonably well by both LES and RANS models. RANS models supplemented with a semi-empirical closure are also found to provide reasonable estimate of the fluctuating pressure loads at the wall.     

Investigation of steady laminar flows with an initial disturbance

Results of a parametric study of steady asymmetric flows are analyzed. Three-dimensional unsteady equations of hydromechanics for a compressible medium are solved by a time-dependent method. The range of the characteristic Reynolds number Re = 60–350 is considered. It is shown that a symmetric flow becomes asymmetric at Re = 90. This value can be considered as a threshold value for air. In the examples considered, the upper separation region is always smaller than the lower separation region owing to flow asymmetry in the vicinity of the left boundary of the domain of integration. The dependence of the separation region size on the Reynolds number is found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 11–19, May–June, 2008.     

Parametric study of unsteady laminar flows

Results of a parametric study of unsteady laminar flows are analyzed. Three-dimensional unsteady equations of hydromechanics for a compressible medium are solved. The range of the characteristic Reynolds number Re = 400–900 is considered. It is demonstrated that the laminar flow in a plane channel ceases to be steady at Re = 415. As the Reynolds number increases, the unsteady processes become more intense, disturbances penetrate inward the channel, and separation zones lose their stability. In the vicinity of the channel exit, however, the flow tends to stabilize, though it remains unsteady. No transition to a turbulent flow occurs in the examined range of Reynolds numbers.     

Direct Computations of Boundary Layers Distorted by Migrating Wakes in a Linear Compressor Cascade

A Direct Numerical Simulation (DNS) of flow in the V103 Low-Pressure (LP) compressor cascade with incoming wakes was performed. The computational geometry was chosen largely in accordance with the setup of the experiments performed by Hilgenfeld and Pfitzner (J Turbomach 126:493–500, 2004) at the University of the Armed Forces in Munich. The computations were carried out on the NEC-SX8 in Stuttgart using 64 processors and 85 million grid points. The incoming wakes stemmed from a separate DNS of incompressible flow around a circular cylinder with a Reynolds number of Re _d  = 3300 (based on mean inflow velocity and cylinder diameter). The boundary layer along the suction surface of the blade was found to separate and roll up due to a Kelvin–Helmholtz instability triggered by the periodically passing wakes. Inside the rolls further transition to turbulence was found to occur. The boundary-layer flow along the pressure surface did not separate, instead it underwent by-pass transition.     

Numerical study of the unsteady flow and heat transfer in channels with periodically mounted square bars

 Numerical investigations of unsteady laminar flow and heat transfer in a channel of height H with periodically mounted square bars of height d = 0.2H arranged side by side to the approaching flow have been conducted for different transverse separation distances of the bars. Five cases with transverse separation distance of 0, 0.5, 1, 1.5 and 2d for a Reynolds number of 300 in a channel with a periodicity length of 2H were studied. The unsteady Navier–Stokes equations and the energy equation have been solved by a finite volume code with staggered grids combined with the SIMPLEC algorithm and a fine grid resolution. Due to the arrangement of bars detached from the channel walls the flow is unsteady with vortex shedding from the bars. The amplitude and mean values of the drag coefficients, skin friction coefficients, friction factor and Nusselt numbers have a strong dependence of the transverse separation distance of the bars. Received on 28 January 2000     

Time-resolved stereo PIV measurements of shock–boundary layer interaction on a supercritical airfoil

Time-resolved stereo particle-image velocimetry (TR-SPIV) and unsteady pressure measurements are used to analyze the unsteady flow over a supercritical DRA-2303 airfoil in transonic flow. The dynamic shock wave–boundary layer interaction is one of the most essential features of this unsteady flow causing a distinct oscillation of the flow field. Results from wind-tunnel experiments with a variation of the freestream Mach number at Reynolds numbers ranging from 2.55 to 2.79 × 10⁶ are analyzed regarding the origin and nature of the unsteady shock–boundary layer interaction. Therefore, the TR-SPIV results are analyzed for three buffet flows. One flow exhibits a sinusoidal streamwise oscillation of the shock wave only due to an acoustic feedback loop formed by the shock wave and the trailing-edge noise. The other two buffet flows have been intentionally influenced by an artificial acoustic source installed downstream of the test section to investigate the behavior of the interaction to upstream-propagating disturbances generated by a defined source of noise. The results show that such upstream-propagating disturbances could be identified to be responsible for the upstream displacement of the shock wave and that the feedback loop is formed by a pulsating separation of the boundary layer dependent on the shock position and the sound pressure level at the shock position. Thereby, the pulsation of the separation could be determined to be a reaction to the shock motion and not vice versa.     

Prediction of unsteady,separated boundary layer over a blunt body for laminar,turbulent, and transitional flow

The focus of this paper is to study the ability of unsteady RANS‐based CFD to predict separation over a blunt body for a wide range of Reynolds numbers particularly the ability to capture laminar‐to‐turbulent transition. A perfect test case to demonstrate this point is the cylinder‐in‐crossflow for which a comparison between experimental results from the open literature and a series of unsteady simulations is made. Reynolds number based on cylinder diameter is varied from 10⁴ to 10⁷ (subcritical through supercritical flow). Two methods are used to account for the turbulence in the simulations: currently available eddy–viscosity models, including standard and realizable forms of the k–ε model; and a newly developed eddy–viscosity model capable of resolving boundary layer transition, which is absolutely necessary for the type and range of flow under consideration. The new model does not require user input or ‘empirical’ fixes to force transition. For the first time in the open literature, three distinct flow regimes and the drag crisis due to the downstream shift of the separation point are predicted using an eddy–viscosity based model with transition effects. Discrepancies between experimental and computational results are discussed, and difficulties for CFD prediction are highlighted. Copyright © 2004 John Wiley & Sons, Ltd.     

Computation of a drastic flow pattern change in an annular swirling jet caused by a small decrease in inlet swirl

This paper investigates the flow pattern change in an annular jet caused by a sudden change in the level of inlet swirl. The jet geometry consists of an annular channel followed by a specially designed stepped‐conical nozzle, which allows the existence of four different flow patterns as a function of the inlet swirl number. This paper reports on the transition between two of them, called the ‘open jet flow high swirl’ and the ‘Coanda jet flow.’ It is shown that a small sudden decrease of 4% in inlet swirl results in a drastic and irreversible change in flow pattern. The objective of this paper is to reveal the underlying physical mechanisms in this transition by means of numerical simulations. The flow is simulated using the unsteady Reynolds‐averaged Navier–Stokes (URANS) approach for incompressible flow with a Reynolds stress turbulence model. The analysis of the numerical results is based on a study of different forces on a control volume, which consists of the jet boundaries. The analysis of these forces shows that the flow pattern change consists of three different regimes: an immediate response regime, a quasi‐static regime and a Coanda regime. The simulation reveals that the pressure–tangential velocity coupling during the quasi‐static regime and the Coanda effect at the nozzle outlet during the Coanda regime are the driving forces behind the flow pattern change. These physical mechanisms are validated with time‐resolved stereo‐PIV measurements, which confirm the numerical simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

Nature of the surface heart transfer fluctuation in a hypersonic separated turbulent flow

This paper presents the results of an experimental study of the unsteady nature of a hypersonic separated turbulent flow. The nomimal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5×10⁷/m. The separated flow was generated using finite span forward facing steps. An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make multi-channel measurements of the fluctuating surface heat trtansfer within the separated flow. Conditional sampling analysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness. The compression waves converge into a single leading shock beyond the boundary layer. The shock structure is unsteady and undergoes large-scale motion in the streamwise direction. The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave. There exists a wide band of frequency of oscillations of the shock system. Most of the frequencies are in the range of 1–3 kHz. The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave. This intermittent phenomenon is considered as the consequence of the large-scale shock system oscillations. Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed. The project supported by National Natural Science Foundation of China     

Investigation of flow and heat transfer in corrugated-undulated plate heat exchangers

 An experimental and numerical investigation of heat transfer and fluid flow was conducted for corrugated-undulated plate heat exchanger configurations under transitional and weakly turbulent conditions. For a given geometry of the corrugated plates the geometrical characteristics of the undulated plates, the angle formed by the latter with the main flow direction, and the Reynolds number were made to vary. Distributions of the local heat transfer coefficient were obtained by using liquid-crystal thermography, and surface-averaged values were computed; friction coefficients were measured by wall pressure tappings. Overall heat transfer and pressure drop correlations were derived. Three-dimensional numerical simulations were conducted by a finite-volume method using a low-Reynolds number k–ɛ model under the assumption of fully developed flow. Computed flow fields provided otherwise inaccessible information on the flow patterns and the mechanisms of heat transfer enhancement. Received on 5 February 1999     

Numerical modeling of three-dimensional shear flow in a cavity with moving lids

The results of a numerical calculation of the unsteady shear flow of a viscous incompressible fluid in a cavity with moving lids are given for various geometries and values of the Reynolds number. The paper describes the numerical method developed for solving the three-dimensional unsteady Navier-Stokes equations in a regular domain and also the direct method of parallelling the solution for a multiprocessor computer with distributed memory. Test results are compared with the results of other authors. Moscow, Nizhnii Novgorod, Marseille. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 41–49, May–June, 1998.     

Numerical Simulation of the Flow Around an Infinitely Long Circular Cylinder in the Transition Regime

The development of three-dimensional structures and the succeeding transition to turbulence occurs in the wake of a circular cylinder at Reynolds numbers 190≤Re≤330. This regime is investigated numerically by means of a spectral element method. Earlier numerical works aimed mainly at reproducing characteristic wake patterns observed in experiments. Small sizes of computational domains and short integration times were chosen to save computational resources. Consequently, the quantitative results show a considerable scatter. Within this work, a step by step approach to highly accurate direct numerical simulations is described. Thorough studies of the effect of resolution and blockage are performed in the laminar, two-dimensional regime, resulting in Reynolds number relationships that exactly reproduce experimental data. Based on these results, a stability analysis is performed to obtain wavelengths that are unstable against spanwise perturbations and the critical Reynolds number for the onset of three-dimensionality. The most unstable wavelengths of the “mode A” and “mode B” instabilities and its multiples are used as periodicity length for direct numerical simulations. Effects of integration time, resolution in streamwise as well as spanwise directions, and periodicity length on the flow quantities are studied. Numerically obtained Reynolds number relationships of Strouhal number and base-pressure coefficients that fit accurately within measured results are given for the first time. Curves for drag and lift coefficients are provided and compared with previous numerical studies. Furthermore, physical interpretations of the wake transition are discussed. Since the separation of physical features and effects of experimental arrangements are frequently an open question, our numerical results are able to supply a contribution to the understanding of the physics of cylinder flow. Received 12 September 2000 and accepted 26 June 2001     

Ligament formation in sheared liquid–gas layers

We perform numerical simulations of two-phase liquid–gas sheared layers, with the objective of studying atomization. The Navier–Stokes equations for two-dimensional incompressible flow are solved in a periodic domain. A volume-of-fluid method is used to track the interface. The density ratio is kept around 10. The calculations show good agreement with a fully viscous Orr–Sommerfeld linear theory over several orders of magnitude of interface growth. The nonlinear development shows the growth of finger-like structures, or ligaments, and the detachment of droplets. The effect of the Weber and Reynolds numbers, the boundary layer width and the initial perturbation amplitude are discussed through a number of typical cases. Inversion of the liquid boundary layer is shown to yield more readily ligaments bending upwards and is thus more likely to produce droplets.     

Reynolds Stress Budgets in Couette and Boundary Layer Flows

Reynolds stress budgets for both Couette and boundary layer flows are evaluated and presented. Data are taken from direct numerical simulations of rotating and non-rotating plane turbulent Couette flow and turbulent boundary layer with and without adverse pressure gradient. Comparison of the total shear stress for the two types of flows suggests that the Couette case may be regarded as the high Reynolds number limit for the boundary layer flow close to the wall. The limit values of turbulence statistics close to the wall for the boundary layer for increasing Reynolds number approach the corresponding Couette flow values. The direction of rotation is chosen so that it has a stabilizing effect, whereas the adverse pressure gradient is destabilizing. The pressure-strain rate tensor in the Couette flow case is presented for a split into slow, rapid and Stokes terms. Most of the influence from rotation is located to the region close to the wall, and both the slow and rapid parts are affected. The anisotropy for the boundary layer decreases for higher Reynolds number, reflecting the larger separation of scales, and becomes close to that for Couette flow. The adverse pressure gradient has a strong weakening effect on the anisotropy. All of the data presented here are available on the web [36]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Computational aeroelastic simulations of self-sustained pitch oscillations of a NACA0012 at transitional Reynolds numbers

The phenomenon of low amplitude self-sustained pitch oscillations in the transitional Reynolds number regime is studied numerically through unsteady, two-dimensional aeroelastic simulations. Based on the experimental data, simulations have been limited in the Reynolds number range 5.0×10⁴<Re_c<1.5×10⁵. Both laminar and URANS calculations (using the SST k–ω model with a low-Reynolds-number correction) have been performed and found to produce reasonably accurate limit cycle pitching oscillations (LCO). This investigation confirms that the laminar separation of the boundary layer near the trailing edge plays a critical role in initiating and sustaining the pitching oscillations. For this reason, the phenomenon is being labelled as laminar separation flutter. As a corollary, it is also shown that turbulence tends to inhibit their existence. Furthermore, two regimes of LCO are observed, one where the flow is laminar and separated without re-attachment, and the second for which transition has occurred followed by turbulent re-attachment. Finally, it is established that the high-frequency, shear instabilities present in the flow which lead to von Kármán vortex shedding are not crucial, nor necessary, to the maintaining mechanism of the self-sustained oscillations.