Effect of span length and temperature on the 3-D confined flow around a vortex promoter

This article presents a numerical study on the influence of span length and wall temperature on the 3-D flow pattern around a square section vortex promoter located inside a micro-channel in the low Reynolds number regime. The first objective of the work is to quantify the critical Reynolds number that defines the onset of vortex shedding and to identify the different regimes that appear as a function of the channel aspect ratio (span to height ratio). We found that the critical Reynolds number for the onset of the Karman street regime increases as the aspect ratio decreases. In particular, for the aspect ratio of 1/2 the critical Reynolds number is nearly six times the critical Reynolds number of the 2-D problem. An intermediate oscillating regime between the steady and the Karman street solutions was also found to exist within a rather wide range of Reynolds numbers for small channel aspect ratios. The second objective was to investigate the influence of the vortex promoter wall temperature on both vortex shedding and flow pattern. This has practical engineering implications because the working fluid considered in the article is water that has a viscosity that depends significantly on temperature and promotes a strong coupling between the momentum and energy equations that influences the system behaviour. Results indicate that high surface temperature on the prism promotes the onset of the Karman street, suggesting design guidelines for micro-channel based heat sinks that make use of vortex promoters.     

Experimental study on the confined 3D laminar flow past a square prism with a high blockage ratio

An experimental PIV study is presented that addresses the confined 3D laminar flow behaviour past a square prism. The Reynolds number (Re), based on prism cross-section height varies between 100 and 256. The channel aspect ratio is 1/1 and the blockage ratio is 1/2.5. This geometry is representative of a passive method to enhance mixing in otherwise highly ordered laminar channel flow. It is found that the lateral walls exert a strong effect on the flow behaviour with two main consequences: (a) the onset of vortex shedding is delayed to a Re in the vicinity of 170, as opposed to the unconfined case where the critical Re is reported to be between 50 and 60 and (b) transition from the steady closed recirculation bubble regime to the vortex shedding regime is not abrupt. In particular, there is a range of Re for which the closed recirculation bubble pulsates with increasing amplitude prior to the onset of the Karman street regime. The experimental results are supported by numerical computations, and the relation of the results with the practical design of engineering systems is also discussed.     

Effects of aspect ratio and orientation on the wake characteristics of low Reynolds number flow over a triangular prism

The wake characteristics of unconfined flows over triangular prisms of different aspect ratios have been numerically analysed in the present work. For this purpose, a fixed Cartesian-grid based numerical technique that involves the porous medium approach to mimic the effect of solid blockage has been utilised. Correspondingly, laminar flow simulations ranging from the sub-critical regime (before the onset of vortex shedding) to the super-critical regime have been considered here within the limits of two-dimensionality. In the sub-critical regime, correlations relating the wake bubble length with Reynolds number (Re) have been proposed for various aspect ratios. Also, the effects of aspect ratio and Reynolds Number on the drag force coefficient (C_D) have been characterised for two different geometrical orientations of the prism (base or apex facing the flow). Subsequently, the critical Reynolds number at the onset of vortex shedding has been predicted for each of the aspect ratio considered, by an extrapolation procedure. The unsteady flow characteristics of the super-critical regime are finally highlighted for different aspect ratios and triangular orientations considered in the study.     

有限长正方形棱柱绕流的双稳态现象

运用大涡模拟(large eddy simulation, LES) 方法对有限长正方形棱柱绕流进行了研究. 棱柱高宽比为5,一端固定于平板上,另一端为自由端. 平板表面边界层厚度可忽略不计. 基于自由流速度和棱柱宽度的雷诺数为3 900. 由计算结果发现,有限长棱柱的尾流在自由端后下扫流的作用下具有显著的三维性. 尾流中会交替出现两种典型的流动状态:一种是柱体两侧所脱落的涡在尾流中形成类似卡门涡街的交错排列状态;另一种是柱体尾流中涡团对称分布且无明显的周期性. 两种典型流动状态的交替出现,对有限长棱柱气动力特性有直接的影响.     

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.     

Instantaneous structure and statistical feature of unsteady flow in a channel obstructed by a square rod

A two-dimensional numerical computation has been made for an unsteady flow in a channel obstructed by an inserted square rod. The results of the computation made for the flow with a parabolic inlet velocity profile at a specific value of channel Reynolds number are analyzed in detail. The obtained results reveal that momentum transfer is enhanced due to the apparent shear stress resulting from the nonzero value of cross-correlation between the streamwise and normal components of fluctuating velocity, , just as in turbulent shear flows, although the studied flow is quite different from turbulent flows in the sense that it is highly periodical and therefore free from randomness. This periodicity leads to a quick recovery of the velocity defect in some region of the wake of the rod. Special attention is paid to the time variation of flow structure. The crisscross motion of the Karman vortex previously found to occur is discussed again, and how it appears is explained in terms of the interaction between the Karman vortex and the disturbed wall shear layer. In the discussion of this relationship, wavering motion of the separation vorticity layers formed on both sides of the rod and the periodic formation of an isolated vortex island from the lifted tip of the wall vorticity layer are analyzed. The vortex island is found to play an important role not only for the occurrence of the crisscross motion of Karman vortex but also for the generation of the nonzero value of .     

圆柱尾迹涡的三维演化及结构特征

应用无粘涡丝运动学理论和局部诱导近似(LIA)方法,以Lagrange观点跟踪涡丝在背景流场中运动,用数值方法研究了中等Re数(≈10~3)下圆柱分离尾迹中Kármán涡和涡辫区涡丝的三维演化的机制和动力学过程,及其结构特征。背景流场考虑为尾迹时间平均速度流场和Kármán涡街流场。初始展向小扰动为指数形式和谐波形式。结果指出:Kármán涡和涡辫区中的涡丝具有展向不稳定性,形成流向涡量。在尾迹的初期输运过程中,表现出有序的大、小尺度涡结构。并进一步分析了其产生的机理。     

Side‐wall effects on the bifurcation of the flow through a sudden expansion

Three‐dimensional computations have been performed to study the flow through a symmetric sudden expansion with an expansion ratio of 3 at low Reynolds numbers. The aspect ratio of the flow channel is allowed to vary within a wide range to examine its influence on the flow which bifurcates from a symmetric state to an asymmetric state. The results reveal that the critical Reynolds number of the symmetry‐breaking bifurcation increases while the aspect ratio is reduced. The flow behaviour near the side walls is illustrated by using limiting streamlines. The origin of the singular points identifiable on the side wall can be traced back to the recirculating flows and the relevant reattachment/separation points in the core of the channel. It is seen that the determination of the exact critical Reynolds number is not trivial because it depends on how to define asymmetric flow. Computations have also been conducted to show that a slight asymmetry in the channel geometry causes a smooth transition from symmetric to non‐symmetric states. Copyright © 2007 John Wiley & Sons, Ltd.     

On the interplay between inertial and viscoelastic effects for the flow in weakly modulated channels

The flow inside a spatially modulated channel is examined for viscoelastic fluids of the Oldroyd‐B type. The lower wall is flat and the upper wall is sinusoidally modulated. The modulation amplitude is assumed to be small. Thus, a regular perturbation expansion of the flow field coupled to a variable‐step finite‐difference scheme is used to solve the problem. Convergence and accuracy assessment against earlier experimental results indicate that there is a significant range of validity of the perturbation approach. The influences of wall geometry, inertia and viscoelasticity on the flow kinematics and stresses are investigated systematically. In particular, the interplay between the flow and fluid parameters effects on the conditions for the onset of backflow, number of vortices, their size and location is revealed. The distance between the flow separation and reattachment locations identifies the vortex size. Non‐monotonic dependence of the vortex size on elasticity is reported. The critical conditions for the onset of negative elasticity effects on vortex size are identified. The critical Reynolds number for the onset of backflow initially decreases then levels off or even increases as elasticity increases. For highly elastic fluid and large enough Reynolds number, more than one vortex appear near the lower wall. Copyright © 2005 John Wiley & Sons, Ltd.     

The ‘void’ structure in the wake of a self-oscillating flexible circular cylinder

We discuss the experimental vortex wake of a flexible circular cylinder undergoing vortex-induced vibration at low Reynolds number and a large cylinder aspect ratio. Hydrogen bubbles formed on the cylinder track the von Karman vortex cores. They show a characteristic ‘void’ structure. We propose a vortex skeleton model that includes a pinch-off of opposite-signed cores. Voids occurred at a node in streamwise vibration when close to an antinode in transverse cylinder vibration. A vibration model predicts the ratio of shedding frequency to natural cylinder vibration frequency necessary for void formation at specific spanwise locations.     

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

Period-doubling route to chaos for a swirling flow in an open cylindrical container with a rotating disk

Conclusion The period-doubling route to chaos for a swirling vortex flow in an open cylindrical container, at an aspect ratio of 2, driven by a rotating bottom disk was recognized by using laser-Doppler velocimetry. The onset of periodic motion for the flow was found when Re was in the range between 1850 and 1900. The flow was subharmonically bifurcated into a double-period motion when Re was about 2150. When the Reynolds number was in the range from 2300 to 2400, the flow bifurcated again through the period-doubling mode. When the Reynolds number was further increased, the flow eventually showed chaotic motion. The existence of a free surface promotes the onset of periodicity, and the difference of the critical Reynolds number with and without a free surface was estimated to be about 600.This work was supported by the National Science Council of the Republic of China under grant No. NSC-82-0410-E-002-191     

Effect of homogeneous stratification of the fluid on the dynamics of a kármán vortex street behind a circular cylinder

The experimental results of studying the effect of homogeneous stratification of the fluid on the conditions of generation of a Kárman vortex street [1] developing in the wake of a cylinder in steady horizontal motion are described. In a homogeneous medium at Reynolds numbers Re >5 two symmetrical regions of vorticity of opposite sign are formed behind the cylinder and move together with the latter. As the speed of the cylinder increases, the link between the vortices and the cylinder grows weaker, the vortices are stretched out along the flow and at Re > 40 begin to separate alternately, forming a vortex street in the wake. At first, the frequency of vortex separation increases sharply with increase in Re, but then levels off. It is found that in a uniformly stratified fluid the onset of vortex separation from the moving cylinder is delayed. The dependence of the critical Reynolds number (onset of vortex separation) on the internal Froude number is obtained. The effect of stratification of the fluid on the frequency of separation of the vortices in the Kármán street is investigated. The effect of the Froude number on the dependence of the Strouhal number on the Reynolds number is established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 83–86, January–February, 1986.In conclusion the authors wish to thank A. T. Onufriev for his interest in their work and useful discussions of the results.     

Numerical simulation of confined laminar jet flows

Two‐dimensional incompressible jet development inside a duct has been studied in the laminar flow regime, for cases with and without entrainment of ambient fluid. Results have been obtained for the flow structure and critical Reynolds number values for steady asymmetric jet development and for the onset of temporal oscillations, at various values of the duct‐to‐jet width ratio (aspect ratio). It is found that at low aspect ratios and Reynolds numbers, jet development inside the duct is symmetric. For larger aspect ratios and Reynolds numbers, the jet flow at steady state becomes asymmetric with respect to the midplane, and for still higher values, it becomes oscillatory with respect to time. When entrainment is present, the instabilities of asymmetric development and temporal oscillations occur at a much higher critical Reynolds number for a given aspect ratio, indicating that the stability of the jet flow is higher with entrainment. Copyright © 2000 John Wiley & Sons, Ltd.     

Formation of vortex breakdown in conical–cylindrical cavities

Numerical simulations in confined rotating flows were performed in this work, in order to verify and characterize the formation of the vortex breakdown phenomenon. Cylindrical and conical–cylindrical geometries, both closed, were used in the simulations. The rotating flow is induced by the bottom wall, which rotates at constant angular velocity. Firstly the numerical results were compared to experimental results available in references, with the purpose to verify the capacity of the computational code to predict the vortex breakdown phenomenon. Further, several simulations varying the parameters which govern the characteristics of the flows analyzed in this work, i.e., the Reynolds number and the aspect ratio, were performed. In these simulations, the limits for the transitional regime and the vortex breakdown formation were verified. Steady and transient cases, with and without turbulence modeling, were simulated. In general, some aspects of the process of vortex breakdown in conical–cylindrical geometries were observed to be different from that in cylinders.     

Three-dimensional numerical simulation of the wake flow of an afterbody at subsonic speeds

We numerically investigate the wake flow of an afterbody at low Reynolds number in the incompressible and compressible regimes. We found that, with increasing Reynolds number, the initially stable and axisymmetric base flow undergoes a first stationary bifurcation which breaks the axisymmetry and develops two parallel steady counter-rotating vortices. The critical Reynolds number (Re _cs) for the loss of the flow axisymmetry reported here is in excellent agreement with previous axisymmetric BiGlobal linear stability (BiGLS) results. As the Reynolds number increases above a second threshold, Re _co, we report a second instability defined as a three-dimensional peristaltic oscillation which modulates the vortices, similar to the sphere wake, sharing many points in common with long-wavelength symmetric Crow instability. Both the critical Reynolds number for the onset of oscillation, Re _co, and the Strouhal number of the time-periodic limit cycle, St_sat, are substantially shifted with respect to previous axisymmetric BiGLS predictions neglecting the first bifurcation. For slightly larger Reynolds numbers, the wake oscillations are stronger and vortices are shed close to the afterbody base. In the compressible regime, no fundamental changes are observed in the bifurcation process. It is shown that the steady state planar-symmetric solution is almost equal to the incompressible case and that the break of planar symmetry in the vortex shedding regime is retarded due to compressibility effects. Finally, we report the developments of a low frequency which depends on the afterbody aspect ratio, as well as on the Reynolds and on the Mach number, prior to the loss of the planar symmetry of the wake.     

Numerical analysis on transitions and symmetry-breaking in the wake of a flapping foil

Flying and marine animals often use flapping wings or tails to generate thrust. In this paper, we will use the simplest flapping model with a sinusoidal pitching motion over a range of frequency and amplitude to investigate the mechanism of thrust generation. Previous work focuses on the Karman vortex street and the reversed Karman vortex street but the transition between two states remains unknown. The present numerical simulation provides a complete scenario of flow patterns from the Karman vortex street to reversed Karman vortex street via aligned vortices and the ultimate state is the deflected Karman vortex street, as the parameters of flapping motions change. The results are in agreement with the previous experiment. We make further discussion on the relationship of the observed states with drag and thrust coefficients and explore the mechanism of enhanced thrust generation using flapping motions.     

Numerical analysis of pitching-rotor aerodynamics

The influence of periodic blade pitching on rotor aerodynamics is numerically investigated at a Reynolds number typical of micro-air vehicles. Blade pitching motion is parameterized using three variables, giving rise to a large parameter space that is explored through 74 test cases. Results show that a relevant tuning of pitching variables can lead to an increase in rotational efficiency and thrust, which is found to be primarily related to the occurrence of reversed von Karman street, leading edge vortex (LEV) formation and dynamic stall phenomenon. In addition, for cases where reversed von Karman street occurs, the flow is found to be quasi-two-dimensional, suggesting that quasi-two-dimensional approaches can provide relevant approximations of the global aerodynamics. Overall, the analysis demonstrates that blade pitching can be beneficial to the aerodynamic performance of micro-air vehicles and helps draw guidelines for further improvements of flapping-rotor concepts.     

Self-sustained oscillations of a confined impinging jet

The present paper investigates the dynamics of a laminar plane jet impinging on a flat plate in a channel. An experimental parametric study is carried out to determine the flow regimes at different levels of confinement and Reynolds numbers. For very confined jets, the flow is steady whatever the Reynolds number. The overall structure of the flow is symmetric with respect to the jet axis and is characterized by the presence of recirculation zones at the channel walls. The dynamics is radically different for less confined jets. Above a critical Reynolds number, the flow bifurcates in the form of an oscillating flapping mode of the impinging jet. Analyses of the experimental results provide with a quantitative characterization of this regime in terms of amplitude, wavelength and frequency. This self-oscillating bifurcated flow induces strong sweepings of the target plate by the jet and intense vortex dipole ejections from the impacted wall. Such a regime is expected to be particularly useful in the enhancement of the local heat transfer at relatively low cost in terms of flow rate.