Effect of velocity ratio on the streamwise vortex structures in the wake of a stack

The time-averaged velocity and streamwise vorticity fields within the wake of a stack were investigated in a low-speed wind tunnel using a seven-hole pressure probe. The experiments were conducted at a Reynolds number, based on the stack external diameter, of Re_D=2.3×10⁴. The stack, of aspect ratio AR=9, was mounted normal to a ground plane and was partially immersed in a flat-plate turbulent boundary layer, where the ratio of the boundary layer thickness to the stack height was δ/H≈0.5. The jet-to-cross-flow velocity ratio was varied from R=0 to 3, which covered the downwash, crosswind-dominated and jet-dominated flow regimes. In the downwash and crosswind-dominated flow regimes, two pairs of counter-rotating streamwise vortex structures were identified within the stack wake. The tip vortex pair located close to the free end of the stack, and the base vortex pair located close to the ground plane within the flat-plate boundary layer, were similar to those found in the wake of a finite circular cylinder, and were associated with the upwash and downwash flow fields within the stack wake, respectively. In the jet-dominated flow regime, a third pair of streamwise vortex structures was observed, referred to as the jet-wake vortex pair, which occurred within the jet-wake region above the free end of the stack. The jet-wake vortex pair had the same orientation as the base vortex pair and was associated with the jet rise. The peak vorticity and strength of the streamwise vortex structures were functions of the jet-to-cross-flow velocity ratio. For the tip vortex structures, their peak vorticity and strength reduced as the jet-to-cross-flow velocity ratio increased.     

Exponential vs Algebraic Growth and Transition Prediction in Boundary Layer Flow

For applications regarding transition prediction, wing design andcontrol of boundary layers, the fundamental understanding of disturbancegrowth in the flat-plate boundary layer is an important issue. In thepresent work we investigate the energy growth of eigenmodes andnon-modal optimal disturbances. We present a set of linear governingequations for the parabolic evolution of wavelike disturbances validboth for the exponential and algebraic growth scenario. The base flow istaken as the Falkner–Skan similarity solution with favorable, adverseand zero pressure gradients. The optimization is carried out over theinitial streamwise position as well as the spanwise wave number andfrequency. The exponential growth is maximized in the sense that theenvelope of the most amplified eigenmode is calculated. In the case ofalgebraic growth, an adjoint-based optimization technique is used. Wefind that the optimal algebraic disturbance introduced at a certaindownstream position gives rise to a larger growth than for the optimaldisturbance introduced at the leading edge. The exponential andalgebraic growth is compared and a unified transition-predictionmethod based on available experimental data is suggested.     

Effect of finite cavity width on flow oscillation in a low-Mach-number cavity flow

The current study is focused on examining the effect of the cavity width and side walls on the self-sustained oscillation in a low Mach number cavity flow with a turbulent boundary layer at separation. An axisymmetric cavity geometry is employed in order to provide a reference condition that is free from any side-wall influence, which is not possible to obtain with a rectangular cavity. The cavity could then be partially filled to form finite-width geometry. The unsteady surface pressure is measured using microphone arrays that are deployed on the cavity floor along the streamwise direction and on the downstream wall along the azimuthal direction. In addition, velocity measurements using two-component Laser Doppler Anemometer are performed simultaneously with the array measurements in different azimuthal planes. The compiled data sets are used to investigate the evolution of the coherent structures generating the pressure oscillation in the cavity using linear stochastic estimation of the velocity field based on the wall-pressure signature on the cavity end wall. The results lead to the discovery of pronounced harmonic pressure oscillations near the cavity’s side walls. These oscillations, which are absent in the axisymmetric cavity, are linked to the establishment of a secondary mean streamwise circulating flow pattern near the side walls and the interaction of this secondary flow with the shear layer above the cavity.     

Analysis of Linear Viscous Flow with a Free Surface in a Circular Cylinder Subjected to Small-Amplitude Circular Oscillations

A fluid flow with a free surface inside a circular cylinder subjected to horizontal, circular oscillation was analyzed theoretically and numerically under the assumption of small-amplitude oscillation and high Reynolds number. It was shown that the nature of the oscillatory flow is of a standing-wave type when projected onto the axial plane and of a progressive-wave type when projected onto the azimuthal plane. The Stokes drift motion in the azimuthal direction and the steady streaming velocity at the edge of the bottom- and side-wall boundary layers are then used in the numerical computation for the steady axisymmetric recirculatory flow outside the boundary layers. We have found that the solutions can be well predicted by asymptotic analysis for the full Navier–Stokes equations in the low streaming-Reynolds-number limit. A simple experiment on flow visualization revealed a good agreement in the surface flow pattern on the bottom wall. It also provided steady recirculatory flows that were not much different from the numerical results on the whole. Received 16 August 1999 and accepted 14 January 2000     

Viscous and Inviscid Instabilities of Flow Along a Streamwise Corner

Here we consider the stability of flow along a streamwise corner formed by the intersection of two large flat plates held perpendicular to each other. Self-similar solutions for the steady laminar mean flow in the corner region have been obtained by solving the boundary layer equations for zero and nonzero streamwise pressure gradients. The stability of the mean flow is investigated using linear stability analysis. An eigensolver has been developed to solve the resulting linear eigenvalue problem either in a global mode to obtain an approximation to all the dominant eigenmodes or in a local mode to refine a particular eigenmode. The stability results indicate that the entire spectrum of two-dimensional and oblique viscous modes of a two-dimensional Blasius boundary layer is active in the case of a corner layer as well, but away from the cornerline. In a corner region of finite spanwise extent, the continuous spectrum of oblique modes degenerates to a discrete spectrum of modes of increasing spanwise wave number. The effect of the corner on the two-dimensional viscous instability is small and decreases the growth rate. The growth rate of outgoing oblique disturbances is observed to decrease, while the growth rate of incoming oblique disturbances is enhanced by the corner. This asymmetry between the outgoing and incoming viscous modes increases with increasing obliqueness of the disturbance. The instability of a zero pressure gradient corner layer is dominated by the viscous modes; however, an inviscid corner mode is also observed. The critical Reynolds number of the inviscid mode rapidly decreases with even a small adverse streamwise pressure gradient and the inviscid mode becomes the dominant one. Received 17 March 1998 and accepted 28 April 1999     

Interaction between longitudinal vortices and flat plate boundary layer

The interaction between longitudinal vortices and flat plate boundary layer has been studied numerically for both laminar and turbulent flow situations. The vortices are assumed to be placed in an otherwise two-dimensional boundary layer flow. The flow is assumed to be incompressible and steady. Considering the fact that the velocity, vorticity and temperature gradients in the transverse directions are much larger than the longitudinal (streamwise) gradients for these flows, the original Navier Stokes equations are parabolized in the streamwise direction. A simple model, based on Boussinesq hypothesis, is used for turbulent flow. The discretized equations are then solved step by step in the streamwise direction, using an iterative procedure at each station. Numerical solutions have been obtained for different parameters, such as the Reynolds number, the circulation and the initial position of the vortices. The computed flow patterns and the skin friction coefficient and Stanton number are found to be qualitatively consistent with available experimental results. It is shown that the interaction between the vortices and the boundary layer may severely disturb the boundary layer flow field and thus considerably increase the local skin friction and heat transfer rate on surface of an aircraft.     

Control of vortex formation from a vertical cylinder in shallow water: Effect of localized roughness elements

Vortex formation from a vertical cylinder in shallow water is controlled by placement of a narrow transverse strip of roughness elements on the bed (bottom surface). A technique of high-image-density particle image velocimetry is employed to obtain global, instantaneous representations of the flow patterns, which lead to phase- and time-averaged patterns of streamline topology and Reynolds stress on planes at and above the bed. Near the bed, the overall form of the streamline topology is maintained, even at larger heights of the roughness elements. With increasing height of the elements, the downstream saddle point is further displaced in the streamwise direction. Correspondingly, the streamwise extent of the negative pocket of the streamwise velocity component, i.e., the region of reverse flow along the bed surface, increases substantially in the streamwise direction. The Reynolds stress in the very near-wake, at locations upstream of the roughness elements, is significantly attenuated, even for small height of roughness. This attenuation occurs not only near the bed surface, but also at the midplane of the shallow water wake, and thereby indicates that the consequence of localized roughness is to exert a global influence. In fact, corresponding patterns of instantaneous velocity and vorticity indicate that consistent formation of large-scale vortices in the very near-wake region is attenuated with relatively small surface roughness on the bed. Downstream of the roughness elements, the patterns of Reynolds stress near the bed surface, as well as at the midplane of the water layer, are significantly altered relative to the case of no roughness. Near the bed, highly concentrated patterns of positive and negative Reynolds stress in the absence of roughness give way to lower-level regions of Reynolds stress in the form of alternating concentrations; the particular pattern depends on the height of the roughness elements. At the midplane of the water layer, the Reynolds stress patterns maintain their same overall form, but the extrema of the Reynolds stress concentrations are attenuated in magnitude and are shifted in the downstream direction, with increasing height of the roughness elements. These observations are complemented by patterns of instantaneous velocity and vorticity.     

Mode spectra of natural disturbances in a circular jet and the effect of acoustic forcing

A modal spectrum technique was used to study coherent instability modes (both axisymmetric and azimuthal) triggered by naturally occurring disturbances in a circular jet. This technique was applied to a high Reynolds number (400,000) jet for both untripped (transitional) and tripped (turbulent) nozzle exit boundary layers, with both cases having a core turbulence level of 0.15%. The region up to the end of the potential core was dominated by the axisymmetric mode, with the azimuthal modes dominating further downstream. The growth of the azimuthal modes was observed closer to the nozzle exit for the jet with a transitional boundary layer. Whether for locally parallel flow or slowly diverging flow, even at low levels of acoustic forcing, the inviscid linear theory is seen to be inadequate for predicting the amplitude of the forced mode. In contrast, the energy integral approach reasonably predicts the evolution of the forced mode.     

A bypass wake induced laminar/turbulent transition

The process of laminar to turbulent transition induced by a von Karman vortex street wake, was studied for the case of a flat plate boundary layer. The boundary layer developed under zero pressure gradient conditions. The vortex street was generated by a cylinder positioned in the free stream. An X-type hot-wire probe located in the boundary layer, measured the streamwise and normal to the wall velocity components. The measurements covered two areas; the region of transition onset and development and the region where the wake and the boundary layer merged producing a turbulent flow. The evolution of Reynolds stresses and rms-values of velocity fluctuations along the transition region are presented and discussed. From the profiles of the Reynolds stress and the mean velocity profile, a ‘negative' energy production region along the transition region, was identified. A quadrant splitting analysis was applied to the instantaneous Reynolds stress signals. The contributions of the elementary coherent structures to the total Reynolds stress were evaluated, for several x-positions of the near wall region. Distinct regions in the streamwise and normal to the wall directions were identified during the transition.     

Global instabilities in diverging channel flows

A global stability study of a divergent channel flow reveals features not obtained hitherto by making either the parallel or the weakly non-parallel (WNP) flow assumption. A divergent channel flow is chosen for this study since it is the simplest spatially developing flow: the Reynolds number is constant downstream, and for a theoretical Jeffery?CHamel flow, the velocity profile obeys similarity. Even in this simple flow, the global modes are shown to be qualitatively different from the parallel or WNP. In particular, the disturbance modes are often not wave-like, and the local scale, estimated from a wavelet analysis, can be a function of both streamwise and normal coordinates. The streamwise variation of the scales is often very different from the expected linear variation. Given recent global stability studies on boundary layers, such spatially extended modes which are not wave-like are unexpected. A scaling argument for why the critical Reynolds number is so sensitive to divergence is offered.     

Instability of cylinder wake under open-loop active control

Instability of a wake controlled by a streamwise Lorentz force is investigated through a Floquet stability analysis. The streamwise Lorentz force, which is a two-dimensional control input created by an electromagnetic actuator located on the cylinder surface,adjusts the base flow to affect the three-dimensional wake instability and achieve wake stabilization and transition delay. The instability mode at a Reynolds number Re = 300 can be transformed from B to A with N = 1.0, where N is an interaction number representing the strength of the Lorentz force relative to the inertial force in the fluid. The wake flow is Floquet stable when N increases to 1.3. The spanwise perturbation wavelengths are 3.926 D and 0.822 D in the modes A and B, respectively, where D is the cylinder diameter. In addition, the oscillating amplitudes of drag and lift are reduced with the increase in the interaction number. Particle tracing is used to explore the essential physical mechanism for mode transformation. The path lines show that suppression of flow separation hinders the fluid deformation and rotation, leading to the decrease in elliptic and hyperbolic instability regions, which is the material cause of mode transformation.All of the results indicate that wake stabilization and transition delay can be achieved under open-loop active control via the streamwise Lorentz force.     

Elementary solutions for streaky structures in boundary layers with and without suction

The temporal evolutions of small, streamwise elongated disturbances in the asymptotic suction boundary layer (ASBL) and the Blasius boundary layer (BBL) are compared. In particular, initial perturbations localized (δ-functions) in the wall-normal direction are studied, corresponding to an axi-symmetric jet coming out of a plane parallel to the flat plate. Analytical solutions are presented for the wall-normal and streamwise velocities in the ASBL case whereas both analytical and numerical methods are used for the BBL case. The initial position of the perturbation and its spanwise wave number are varied in a parameter study. We present results of maximum amplitudes obtained, the time to reach them, their position and optimal spanwise scales. Free-stream disturbances are shown to migrate towards the wall and reach their (negative) optimum inside the boundary layer. The migration is faster for the ASBL case and a larger amplitude is reached than for the BBL. For perturbations originating inside the boundary layer the amplitudes are overall larger and show the phenomenon of overshoot, i.e. positive amplitudes moving out of the boundary layer. The overall largest amplitudes are obtained for the BBL case, as in other studies, but it is shown that for free-stream disturbances initiated somewhere downstream the leading edge streak growth may be amplified due to suction since in the BBL the disturbance mainly advects above the boundary layer.     

Stability analysis of the onset of vortex shedding for wakes behind flat plates

Above a critical Reynolds number, wake flows behind flat plates become globally unstable, the leading modal instability in this case is known as Kelvin–Helmholtz mechanism. In this article, both local and BiGlobal linear instability analyses are performed numerically to study the onset of the shedding process. Flat plates with different base shapes are considered to assess geometry effects, and the relation between the critical shedding Reynolds number, \(Re_\mathrm{cr}\), and the boundary layer thickness is studied. Three types of base shapes are used: square, triangular and elliptic. It is found that the base shape has a great impact on the growth rate of least stable disturbance mode, thus would influence \(Re_\mathrm{cr}\) greatly, but it has little effect on the vortex shedding frequency. The shedding frequency is determined mainly by boundary layer thickness and has little dependence on the Reynolds number and base shape. We find that for a fixed Reynolds number, increasing boundary layer thickness acted in two ways to modify the global stability characteristics: It increases the length of the absolute unstable region and it makes the flow less locally absolutely unstable in the near-wake region, and these two effects work against each other to destabilize or stabilize the flow.     

Global pressure relaxation for laminar two-dimensional internal flow

This study extends the reduced Navier–Stokes (RNS) global pressure relaxation procedure developed by Rubin and co-workers for external flow to internal flow applications. The streamwise pressure gradient is split into a backward-differenced or initial value component, as in boundary layer marching, and a forward-differenced or boundary value component that represents the elliptic downstream effects. The streamwise convection terms are upwind-differenced and all other streamwise derivatives are backward-differenced. We thus obtain a standard boundary layer marching technique imbedded in a conventional line relaxation technique. For compressible flow the pressure iteration determines the interior flow interation as well as the inlet mass flux that is consistent with the outflow pressure boundary condition. Results have been computed for incompressible flow in both rectangular and curved channels, and for subsonic compressible flow in the simulation of an aerofoil in a wind tunnel. Converged solutions were obtained over a range of Reynolds numbers generating small to moderately large separation bubbles.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Global two-dimensional stability measures of the flat plate boundary-layer flow

The stability of the two-dimensional flat plate boundary-layer is studied by means of global eigenmodes. These eigenmodes depend both on the streamwise and wall-normal coordinate, hence there are no assumptions on the streamwise length scales of the disturbances. Expanding the perturbation velocity field in the basis of eigenmodes yields a reduced order model from which the stability characteristics of the flow, i.e. the initial condition and forcing function leading to the largest energy growth, are extracted by means of non-modal analysis. In this paper we show that, even when performing stability analysis using global eigenmodes, it is not sufficient to consider only a few of the least damped seemingly relevant eigenmodes. Instead it is the task of the optimization procedure, inherent in the non-modal analysis, to decide which eigenmodes are relevant. We show that both the optimal initial condition and the optimal forcing structure have the form of upstream tilted structures. Time integration reveals that these structures gain energy through the so called Orr mechanism, where the instabilities extract energy from the mean shear. This provides the optimal way of initiating Tollmien–Schlichting waves in the boundary layer. The optimal initial condition results in a localized Tollmien–Schlichting wavepacket that propagates downstream, whereas the optimal forcing results in a persistent Tollmien–Schlichting wave train.     

The linearized pressure Poisson equation for global instability analysis of incompressible flows

The linearized pressure Poisson equation (LPPE) is used in two and three spatial dimensions in the respective matrix-forming solution of the BiGlobal and TriGlobal eigenvalue problem in primitive variables on collocated grids. It provides a disturbance pressure boundary condition which is compatible with the recovery of perturbation velocity components that satisfy exactly the linearized continuity equation. The LPPE is employed to analyze instability in wall-bounded flows and in the prototype open Blasius boundary layer flow. In the closed flows, excellent agreement is shown between results of the LPPE and those of global linear instability analyses based on the time-stepping nektar++, Semtex and nek5000 codes, as well as with those obtained from the FreeFEM++ matrix-forming code. In the flat plate boundary layer, solutions extracted from the two-dimensional LPPE eigenvector at constant streamwise locations are found to be in very good agreement with profiles delivered by the NOLOT/PSE space marching code. Benchmark eigenvalue data are provided in all flows analyzed. The performance of the LPPE is seen to be superior to that of the commonly used pressure compatibility (PC) boundary condition: at any given resolution, the discrete part of the LPPE eigenspectrum contains converged and not converged, but physically correct, eigenvalues. By contrast, the PC boundary closure delivers some of the LPPE eigenvalues and, in addition, physically wrong eigenmodes. It is concluded that the LPPE should be used in place of the PC pressure boundary closure, when BiGlobal or TriGlobal eigenvalue problems are solved in primitive variables by the matrix-forming approach on collocated grids.     

Dynamics of vortical structures in cylinder/wall interaction with moderate gap ratio

The interaction between the wake of a transverse circular cylinder and the underlying flat-plate boundary layer with a moderate gap ratio G/D=1.0 is investigated using both hydrogen-bubble-based and PIV-based visualization techniques. The spanwise rollers in the cylinder wake are found to be capable of inducing secondary vortices in the near-wall region. The mutual induction from the counter-clockwise rollers, which are closer to the wall, plays a primary role, so that these secondary vortices present linear lift-up motion at first. Their subsequent evolution dominantly determines the characteristics of the wake/boundary-layer interaction. Two different vortex interaction scenarios are observed: the secondary vortices can be either entrained into the rollers or pushed down towards the wall. This leads to a rapid three-dimensional destabilization process, through which streamwise vortices are generated. And it is suggested that these streamwise vortices are the dominant structures to promote the following boundary layer transition.     

Study of the Interaction of a Streamwise-Oscillating Shock with the Boundary Layer on a Flat Plate

The distinctive features of the flow in the region of interaction between an oscillating shock and a flat-plate boundary layer are studied for the laminar, transitional, and turbulent flow regimes. The flow patterns and the pressure and heat transfer distributions in the interaction region are analyzed at different intensities, frequencies, and amplitudes of the oscillating shock. The results of the interaction of oscillating and steady shocks with the flat-plate boundary layer are compared.