Numerical simulation for rotating internal weakly viscoelastic flows in rectangular ducts

The present work develops a numerical method for the solution of rotating internal weakly viscoelastic flows in rectangular ducts for dimensionless parameters such as the Reynolds, Rossby and Weissenberg numbers, taken respectively in the intervals between 171 and 12000, 0.047 and 1/12 and up to 1/10000. It is shown that the usual counter‐rotating double‐vortex configuration of secondary flow breaks down with the increase of the Reynolds number (over the threshold of 171). For higher Reynolds numbers such as 7500 and 12000 the secondary flow diffuses to the interior of the duct where it assumes a fully developed configuration and the transition to the turbulence structure is observed. The Sobolev norms increase almost proportionally to the increase of the Reynolds number, and play an essential role for more complex problems involving transition to turbulence modelling. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Circular Couette flow with pressure-driven axial flow and a porous inner cylinder

In a rotating filter separator a suspension is introduced at one end of the annulus between a rotating porous inner cylinder and a fixed impermeable outer cylinder. The filtrate is removed through the inner cylinder and the concentrate is removed from the opposite end of the annulus from which the suspension entered. The flow in a rotating filter separator is circular Couette flow with a pressure-driven axial flow and a suction boundary condition at the inner cylinder. Flow visualization was used to determine the effect of the Taylor number, axial Reynolds number, and radial Reynolds number on the types of flows present in the annulus. A rich variety of secondary vortical flows appear, depending upon the flow parameters. The radial inflow at the inner cylinder delays the appearance of supercritical circular Couette flow and prevents the appearance of certain flow regimes that have a helical vortex structure. Nevertheless, the average azimuthal velocity measured using laser Doppler velocimetry indicates that the velocity profile is nearly the same for all supercritical flow regimes.This work was supported by a grant from The Whitaker Foundation     

单柱单锥型液—液旋流分离管内流场的LDV诊断

应用二维激光多普勒仪（ＬＤＶ）对一种单柱单锥型液－液旋流分离管内流场进行了测量,考察了流量、溢流比、压力比和气芯等参数对流场的影响。测量结果表明：切向速度分布呈典型的Ｒａｎｋｉｎｅ涡结构,沿轴向衰减很少,表明所用锥角是合适的;因该旋流管的水力直径较大,切向速度的总体水平较低,由于对了离特性带来了不利影响。此外,没有观察到切向速度分布的的双峰分布现象。轴向速度的总体水平较低,尤其是在锥形管的上游更为     

Asymmetric flows in planar symmetric channels with large expansion ratio

A continuation method has been used with a finite element grid and a geometric perturbation to compute two successive symmetry breaking flow transitions with increasing Reynolds number in flow of generalized Newtonian fluids through a sudden planar expansion. With an expansion ratio of 16, the onset Reynolds number is particularly sensitive to small geometric asymmetry and the critical Reynolds numbers for the two successive flow transitions are found to be very close. These transitions are delayed to higher onset Reynolds numbers by increasing the degree of pseudoplasticity. This trend is observed experimentally as well in this work and may be attributed to the competing effects of shear thinning and inertia on the size of the corner vortex before the symmetry breaking flow transition. After the second transition with an expansion ratio of 16, the two large staggered vortices on opposite walls occupy most of the transverse dimension so that the core flow between the vortices appears as a thin jet oscillating along the flow direction. This is more pronounced for the pseudoplastic liquid. After the second transition, the degree of flow asymmetry at a given location downstream of the expansion plane is larger for the pseudoplastic liquid than for the Newtonian liquid at comparable Reynolds numbers. The last feature is also evident in the experimentally observed velocity profiles. Copyright © 2002 John Wiley & Sons, Ltd.     

PIV MEASUREMENT OF CLOSE IMPINGING JET ON FLAT PLATE 1)

运用时间分辨粒子成像测速系统(time-resolved particle image velocimetry, TR-PIV)对近距离下射流冲击平板时的流场进行了直接测量, 通过对两个正交的平面流场开展测量, 揭示了冲击距离和雷诺数对射流间隙内三维流动特征及涡系结构演化规律的影响. 结果表明: 射流间隙存在三种典型的涡系结构, 分别为双涡环模式、单涡环模式和卷吸模式, 但在大流量湍流状态下, 射流可能会冲破涡环, 形成随机的高速出流, 各流动模式的出现主要与射流流态及壁面约束作用有关. 运用涡量分析对三种典型涡系结构的能量传递和损失特性进行了比较. 结果表明: 近距离冲击时, 射流的能量通过涡环模式向外传递. 在双涡环模式下, 两个涡环的旋向相反, 端面的约束作用使得两个涡环都被严格约束在射流棒端面之内, 且一次涡环强度显著大于二次涡环强度. 最后, 运用本征正交分解方法对射流间隙内的流动模态及其能量分布进行了分析. 单涡和双涡模式前十阶模态分析结果表明: 能量脉动在较低阶时即以配对的模式出现, 这表明一次涡环与二次涡环均具有良好的对称性, 同时在双涡模式中, 一次涡环是占主导作用的大尺度流动结构. 卷吸模式的前三阶模态分析表明: 射流的能量集中在射流上游, 能量随紊动扩散急剧衰减.     

近距离下射流冲击平板PIV实验研究

运用时间分辨粒子成像测速系统(time-resolved particle image velocimetry, TR-PIV)对近距离下射流冲击平板时的流场进行了直接测量, 通过对两个正交的平面流场开展测量, 揭示了冲击距离和雷诺数对射流间隙内三维流动特征及涡系结构演化规律的影响. 结果表明: 射流间隙存在三种典型的涡系结构, 分别为双涡环模式、单涡环模式和卷吸模式, 但在大流量湍流状态下, 射流可能会冲破涡环, 形成随机的高速出流, 各流动模式的出现主要与射流流态及壁面约束作用有关. 运用涡量分析对三种典型涡系结构的能量传递和损失特性进行了比较. 结果表明: 近距离冲击时, 射流的能量通过涡环模式向外传递. 在双涡环模式下, 两个涡环的旋向相反, 端面的约束作用使得两个涡环都被严格约束在射流棒端面之内, 且一次涡环强度显著大于二次涡环强度. 最后, 运用本征正交分解方法对射流间隙内的流动模态及其能量分布进行了分析. 单涡和双涡模式前十阶模态分析结果表明: 能量脉动在较低阶时即以配对的模式出现, 这表明一次涡环与二次涡环均具有良好的对称性, 同时在双涡模式中, 一次涡环是占主导作用的大尺度流动结构. 卷吸模式的前三阶模态分析表明: 射流的能量集中在射流上游, 能量随紊动扩散急剧衰减.      

Numerical study of adverse pressure gradient generation over a flat plate using a rotating cylinder

Generating an adverse pressure gradient (APG), using a rotating cylinder in the proximity of a plane wall under a laminar freestream flow, is studied numerically in this work. The magnitude of the generated APG is a function of the gap, G, between the cylinder and the wall, and the rotational speed of the cylinder, Ω. The flow in such a configuration is characterized by periodic transient vortex shedding at high Reynolds number. A numerical model for the computation of the transient flow for this configuration is developed using the ANSYS CFD simulation tool. The model is validated against published experimental and numerical data for similar flow configurations and excellent agreement is observed. A parametric study is carried out for different combinations of G and Ω for two different Reynolds numbers of 200 and 1000 to examine the development of the resulting separation bubble due to the generated APG. The mechanism of the boundary layer separation over the plane wall and the corresponding wake dynamics is investigated. Results are presented in terms of the distribution of the pressure coefficient as well as skin friction coefficient along the wall and flow patterns around and downstream of the cylinder in the proximity of the wall. The results of these computations confirm that using a rotating cylinder over a plane wall in a freestream flow is an effective technique to generate a controlled range of adverse pressure gradients.     

Transition from vortex to wall driven turbulence production in the Taylor–Couette system with a rotating inner cylinder

Axisymmetrically stable turbulent Taylor vortices between two concentric cylinders are studied with respect to the transition from vortex to wall driven turbulent production. The outer cylinder is stationary and the inner cylinder rotates. A low Reynolds number turbulence model using the k‐ω formulation, facilitates an analysis of the velocity gradients in the Taylor–Couette flow. For a fixed inner radius, three radius ratios 0.734, 0.941 and 0.985 are employed to identify the Reynolds number range at which this transition occurs. At relatively low Reynolds numbers, turbulent production is shown to be dominated by the outflowing boundary of the Taylor vortex. As the Reynolds number increases, shear driven turbulence (due to the rotating cylinder) becomes the dominating factor. For relatively small gaps turbulent flow is shown to occur at Taylor numbers lower than previously reported. Copyright © 2002 John Wiley & Sons, Ltd.     

Temporal analysis of breakup for a power law liquid jet in a swirling gas

The breakup mechanism and instability of a power law liquid jet are investigated in this study. The power law model is used to account for the non-Newtonian behavior of the liquid fluid. A new theoretical model is established to explain the breakup of a power law liquid jet with axisymmetric and asymmetric disturbances, which moves in a swirling gas. The corresponding dispersion relation is derived by a linear approximation, and it is applicable for both shear-thinning and shear-thickening liquid jets. Analysis results are calculated based on the temporal mode. The analysis includes the effects of the generalized Reynolds number, the Weber number, the power law exponent, and the air swirl strength on the breakup of the jet. Results show that the shear-thickening liquid jet is more unstable than its Newtonian and shear-thinning counterparts when the effect of the air swirl is taken into account. The axisymmetric mode can be the dominant mode on the power law jet breakup when the air swirl strength is strong enough, while the non-axisymmetric mode is the domination on the instability of the power liquid jet with a high We and a low Re _n . It is also found that the air swirl is a stabilizing factor on the breakup of the power law liquid jet. Furthermore, the instability characteristics are different for different power law exponents. The amplitude of the power law liquid jet surface on the temporal mode is also discussed under different air swirl strengths.     

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.     

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.     

Self-sustained global oscillations in a jet in crossflow

A jet in crossflow with an inflow ratio of 3, based on the maximum velocity of the parabolic jet profile, is studied numerically. The jet is modeled as an inhomogeneous boundary condition at the crossflow wall. We find two fundamental frequencies, pertaining to self-sustained oscillations in the flow, using full nonlinear direct numerical simulation (DNS) as well as a modal decomposition into global linear eigenmodes and proper orthogonal decomposition (POD) modes; a high frequency which is characteristic for the shear-layer vortices and the upright vortices in the jet wake, and a low frequency which is dominant in the region downstream of the jet orifice. Both frequencies can be related to a region of reversed flow downstream of the jet orifice. This region is observed to oscillate predominantly in the wall-normal direction with the high frequency, and in the spanwise direction with the low frequency. Moreover, the steady-state solution of the governing Navier?CStokes equations clearly shows the horseshoe vortices and the corresponding wall vortices further downstream, and the emergence of a distinct counter-rotating vortex pair high in the free stream. It is thus found that neither the inclusion of the jet pipe nor unsteadiness is necessary to generate the characteristic counter-rotating vortex pair.     

Analysis of an organized turbulent structure using a pattern recognition technique in a drag-reducing surfactant solution flow

This work presents the investigation for an organized turbulent structure in a drag-reducing flow of dilute surfactant solution by utilizing a particle image velocimetry system to perform the pattern recognition technique on a trajectory in four quadrants of streamwise and wall-normal velocity fluctuations. The pattern recognition is added to a new algorithm in order to directly capture the spatial rotation motion. The Reynolds number based on the channel height and bulk mean velocity was set to 1.5 × 10⁴. Surfactant solution with a weight concentration of 150 ppm was employed and the drag reduction rate was 65%. In the drag-reducing flow, we observe increased frequencies of occurrence of the flow events that correspond to a meandering motion in the wall-normal direction of the high-and low-speed regions. Three findings from investigation of the ensemble-averaged Reynolds shear stress and vortex structure are as follows: (i) the Reynolds shear stress in the large fluctuation range occurs in the narrow region; (ii) Size, strength, arrangement and inclination in the spatial vortex structure in the drag-reducing flow differ from those of the water; and (iii) all trajectory contributions for the wall friction coefficient decrease. Finally, we interpreted that the viscoelasticity characterizing the viscoelastic stress and relaxation time in rheological properties of the flow changes specific elementary vortex for the drag-reducing flow, and the trajectories of each flow pattern change drastically.     

Population,characteristics and kinematics of vortices in a confined rectangular jet with a co-flow

Vortex behavior and characteristics in a confined rectangular jet with a co-flow were examined using vortex swirling strength as a defining characteristic. On the left side of the jet, the positively (counterclockwise) rotating vortices are dominant, while negatively rotating vortices are dominant on the right side of the jet. The characteristics of vortices, such as population density, average size and strength, and deviation velocity, were calculated and analyzed in both the cross-stream direction and the streamwise direction. In the near-field of the jet, the population density, average size and strength of the dominant direction vortices show high values on both sides of the center stream with a small number of counter-rotating vortices produced in the small wake regions close to jet outlet. As the flow develops, the wake regions disappear, these count-rotating vortices also disappear, and the population of the dominant direction vortices increase and spread in the jet. The mean size and strength of the vortices decrease monotonically with streamwise coordinate. The signs of vortex deviation velocity indicate the vortices transfer low momentum to high-velocity region and high momentum to the low velocity region. The developing trends of these characteristics were also identified by tracing vortices using time-resolved particle image velocimetry data. Both the mean tracked vortex strength and size decrease with increasing downstream distance overall. At the locations of the left peak of turbulent kinetic energy, the two-point spatial cross-correlation of swirling strength with velocity fluctuation and concentration fluctuation were calculated. All the correlation fields contain one positively correlated region and one negatively correlated region although the orientations of the correlation fields varied, due to the flow transitioning from wake, to jet, to channel flow. Finally, linear stochastic estimation was used to calculate conditional structures. The large-scale structures in the velocity field revealed by linear stochastic estimation are spindle-shaped with a titling stream-wise major axis.     

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.     

An experimental investigation of upstream flow characteristics of viscoelastic fluids in an annular die entry

Flow behaviour of viscoelastic polymer solutions on the upstream side of an annular die entry has been experimentally investigated and compared with published results on entry flow in circular die. Stable and unstable flow patterns were observed depending on the magnitudes of Reynolds and elasticity numbers. The latter number represents the relative effects of elastic and inertial forces. The stable flow region consists of an elasticity-controlled vortex growth regime and an inertia-controlled divergent flow regime. These two flow regimes have also been observed in circular entry flow. The criteria for the onset of flow instability in an annulus, based on the maximum vortex size, agree qualitatively with various criteria proposed for polymer metls and solutions in circular entry flow. The unstable flow patterns revealed a two-stage instability with a metastable region in between. The first-stage instability is characterised by low frequency disruption of the stationary vortex; while the second-stage instability, which occurs at high Reynolds numbers, is characterised by high frequency random distortion of the flow field.     

Properties of vortices in the self-similar turbulent jet

Instantaneous, two-dimensional velocity measurements were conducted in the axial plane of a self-similar turbulent axisymmetric jet. The velocity fields were high-pass filtered to expose the vortical structures. An automated method was used to identify the radial and axial coordinates of the vortex centers and rotational sense, and to measure their size, circulation, vorticity, and energy. New insights into turbulent jets are obtained by plotting statistical distributions for vortex properties as functions of Reynolds number and radial position. While the probability of finding a vortex is uniform up to the edge of the jet, the strongest eddies in the high-pass filtered field occur near the jet axis. The average circulation is directly proportional to the vortex size. The Reynolds number strongly affects the average vorticity, circulation, and energy of the eddies. However, the normalized curves show a good collapse implying that the jet is indeed self-similar. Results for the left and right half-planes of the jet are also presented. Interestingly, we find that contrary to customary drawings of jet flows, a substantial number of both clockwise and counter-clockwise rotating eddies exist on both sides of the jet axis, with almost equal numbers of oppositely rotating vortices close to the jet axis. Further, the disparity in the number of oppositely rotating eddies in each half-plane increases with the eddy size. Nevertheless, these results are consistent with the well-known radial vorticity distribution of axisymmetric jets.     

An experimental investigation of the flow about a sphere rotating in a rivlin-ericksen fluid

Tangential and radial velocity profiles were measured for the flow about a sphere rotating slowly in a Newtonian fluid, contained in a rectangular tank. Velocities were determined from enlarged streak photographs of aluminium particles moving in a collimated “sheet” of light, at several planes through the flow field. Similar velocity profiles were measured for the flow of a 1.50% Natrosol 250 H solution about two spheres of different diameters rotating in two different sized rectangular tanks. A set of velocity distributions were also measured for a sphere rotating in a 0.9% Natrosol 250 H solution. A dye tracer study of the flow about a sphere rotating in this liquid is presented as well. Both Natrosol solutions exhibited viscoelastic behaviour. The Newtonian fluid study was carried out at a Reynolds number of 1.2 and the viscoelastic fluid studies were within the Reynolds number range of 0.05–1.24.The zero shear viscosities of the Natrosol solutions were measured using the falling-sphere method. The non-Newtonian material parameters were obtained by fitting the theoretical curves to the measured velocity data. The values of the elastic and shear thinning parameters for the two liquids obtained in the different geometrical and dynamical situations are compared.     

Experimental investigation of the unsteady structure of a transitional plane wall jet

 A laminar wall jet undergoing transition is investigated using the particle image velocimetry (PIV) technique. The plane wall jet is issued from a rectangular channel, with the jet-exit velocity profile being parabolic. The Reynolds number, based on the exit mean velocity and the channel width, is 1450. To aid the understanding of the global flow features, laser-sheet/smoke flow visualizations are performed along streamwise, spanwise, and cross-stream directions. Surface pressure measurements are made to correlate the instantaneous vorticity distribution with the surface pressure fluctuations. The instantaneous velocity and vorticity field measurements provide the basis for understanding the formation of the inner-region vortex and the subsequent interactions between the outer-region (free-shear-layer region) and inner-region (boundary-layer region) vortical structures. Results show that under the influence of the free-shear-layer vortex, the local boundary layer becomes detached from the surface and inviscidly unstable, and a vortex is formed in the inner region. Once this vortex has formed, the free-shear-layer vortex and the inner-region vortex form a vortex couple and convect downstream. The mutual interactions between these inner- and outer-region vortical structures dominate the transition process. Farther downstream, the emergence of the three-dimensional structure in the free shear layer initiates complete breakdown of the flow. Received: 8 November 1995/Accepted: 6 November 1996