Study on the instability in separating–reattaching flow over a surface-mounted rib

A numerical study on the flow structure and instability in the separated–reattached flow over a surface-mounted rib at Re = 1000 is performed using large eddy simulation. It is found that the phenomenon of vortex pairing, which has been extensively observed in similar flows, exists in the separation zone. Based on the spectral analysis, the Kelvin–Helmholtz (K-H) instability of shear layer at St ≈ 0.361 (St ≡ fh/U₀) and its subharmonic at St ≈ 0.18 are found. It is assumed that the K-H instability reduces to its subharmonic through the vortex pairing. This process is further confirmed by the flow visualisation. The two-dimensional (2D) structures are subjected to sinusoidal undulation along the spanwise and observed to undergo helical pairing process, which is attributed to the transformation of 2D structures into 3D. However, the low frequency due to flapping of the shear layer is not found.     

Particle-turbulence interaction in a boundary layer

Particle-turbulence interaction in wall turbulent flows has been studied. A series of experiments varying particle size, particle density, particle loading and flow Re has been conducted. The results show that the larger polystyrene particles (1100 μm) cause an increase in the number of wall ejections, giving rise to an increase in the measured values of the turbulence intensities and Reynolds stresses. On the other hand, the smaller polystyrene particles (120 μm) bring about a decrease in the number of wall ejections, causing a decrease in the measured intensities and Reynolds stresses. These effects are enhanced as the particle loading is increased. It was also found that the heavier glass particles (88 μm) do not bring about any significant modulation of turbulence. In addition, measurements of the burst frequency and the mean streak-spacing show no significant change with increase in particle loading. Based on these observations, a mechanism of particle transport in wall turbulent flows has been proposed, in which the particles are transported (depending on their size, density and flow Re) by the bursting events of the wall regions.     

Two phase mixed convection Al2O3–water nanofluid flow in an annulus

Heat transfer enhancement of a mixed convection laminar Al₂O₃–water nanofluid flow in an annulus with constant heat flux boundary condition has been studied employing two phase mixture model and effective expressions of nanofluid properties. The fluid flow properties are assumed constant except for the density in the body force, which varies linearly with the temperature (Boussinesq’s hypothesis), thus the fluid flow characteristics are affected by the buoyancy force. The Brownian motions of nanoparticles have been considered to determine the effective thermal conductivity and the effective dynamic viscosity of Al₂O₃–water nanofluid, which depend on temperature. Three-dimensional Navier–Stokes, energy and volume fraction equations have been discretized using the finite volume method while the SIMPELC algorithm has been introduced to couple the velocity–pressure. Numerical simulations have been presented for the nanoparticles volume fraction (?) between 0 and 0.05 and different values of the Grashof and Reynolds numbers. The calculated results show that at a given Re and Gr, increasing nanoparticles volume fraction increases the Nusselt number at the inner and outer walls while it does not have any significant effect on the friction factor. Both the Nusselt number and the friction coefficient at the inner wall are more than their corresponding values at the outer wall.     

Free-stream turbulence effects on vortex-induced vibration of two side-by-side elastic cylinders

The effect of free-stream turbulence on vortex-induced vibration of two side-by-side elastic cylinders in a cross-flow was investigated experimentally. A turbulence generation grid was used to generate turbulent incoming flow with turbulence intensity around 10%. Cylinder displacements in the transverse direction at cylinder mid-span were measured in the reduced velocity range 1.45<U_r0<12.08, corresponding to a range of Reynolds number (Re), based on the mean free-stream velocity and the diameter of the cylinder, between Re=5000–41 000. The focus of the study is on the regime of biased gap flow, where two cylinders with pitch ratio (T/D) varying from 1.17 to 1.90 are considered. Results show that the free-stream turbulence effect is to enhance the vortex-induced force, thus to restore the large-amplitude vibration associated with the lock-in resonance. However, the enhancement is significant at a different Strouhal number (St) for different pitch ratios. When the spacing between two cylinders is relatively small (1.17<T/D<1.50), the enhancement is significant at St≈0.1. When the spacing is increased, the Strouhal number at which the enhancement is significant shifts to St≈0.16. This enlarges the range of reduced velocity to be concerned. An energy analysis showed that free-stream turbulence feeds energy to the cylinder at multiple frequencies of vortex shedding. Therefore, the lock-in region is still of main concern when the approach flow is turbulent.     

Unsteady laminar flow developing in a curved duct

The flow developing in a tightly curved U-bend of square cross section has been investigated experimentally and via numerical simulation. Both long-time averages and time histories of the longitudinal (streamwise) component of velocity were measured using a laser-Doppler velocimeter. The Reynolds number investigated was Re = 1400. The data were obtained at different bend angles, θ, and were confined to the symmetry plane of the bend. At Re = 1400, the flow entering the bend is steady, but by θ = 90° it develops an oscillatory component of motion along the outer-radius wall. Autocorrelations and energy spectra derived from the time histories yield a base frequency of approximately 0.1 Hz for these oscillations. Flow-visualization studies showed that the proximity of the outer-radius wall served to damp the amplitude of the spanwise oscillations.

Numerical simulations of the flow were performed using both steady and unsteady version of the finite-difference elliptic calculation procedure of Humphrey et al. (1977). Although the unsteadiness observed experimentally does not arise spontaneously in the calculations, numerical experiments involving the imposition of a periodic time-dependent perturbation at the inlet plane suggest that the U-bend acts upon the incoming flow so as to damp the amplitude of the imposed oscillation while altering its frequency.

The oscillations observed experimentally, and numerically as a result of the periodic perturbation, have been linked to the formation of Goertler-type vortices of the outer-radius wall in the developing flow. The vortices, which develop as a result of the centrifugal instability of the flow on the outer-radius wall, undergo a further transition to an unsteady regime at higher flow rates.     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of Knudsen number and geometry on gaseous flow and heat transfer in a constricted microchannel

A flow and heat transfer numerical simulation is performed for a 2D laminar incompressible gas flow through a constricted microchannel in the slip regime with constant wall temperature. The effects of rarefaction, creeping flow, first order slip boundary conditions and hydrodynamically/thermally developing flow are assumed. The effects of Knudsen number and geometry on thermal and hydrodynamic characteristics of flow in a constricted microchannel are explored. SIMPLE algorithm in curvilinear coordinate is used to solve the governing equations including continuity, energy and momentum with the temperature jump and velocity slip conditions at the solid walls in discretized form. The resulting velocity and temperature profiles are then utilized to obtain the microchannel C _f Re and Nusselt number as a function of Knudsen number and geometry. The results show that Knudsen number has declining effect on the C _f Re and Nusselt number in the constricted microchannel. In addition, the temperature jump on wall and slip velocity increase with increasing Knudsen number. Moreover, by decreasing the throttle area, the fluid flow characteristics experience more intense variations in the constricted region. To verify the code a comparison is carried out with available results and good agreement is achieved.     

Newtonian and Power-Law fluid flow in a T-junction of rectangular ducts

The aim of the present study is the numerical investigation of the shear-thinning and shear-thickening effects of flow in a T-junction of rectangular ducts. The employed CFD code incorporates the SIMPLE scheme in conjunction with the finite volume method with collocated arrangement of variables. The code enables multi-block computations in domains with multiple apertures, thus coping with the two-block, two-outlet layout of the current 3D computational domain. The shear-thinning and shear-thickening behaviours of the flow are covered by changing the index n of the Power-Law model within a range from 0.20 to 1.25, and the subsequent effects are investigated by means of different flow parameters namely the Reynolds (Re) number and the boundary conditions at the outlets. Results exhibit the extent of the effect of the Re number on the velocity profiles at different positions in the domain for both Newtonian and non-Newtonian cases. Similarly, the trend of the effect of shear-thinning and shear-thickening behaviours on the flow rate ratio between inlet and outlets, in the case of equal pressure imposed on outlets, is shown.     

鱼类自主游动水动力学特性的数值模拟

通过对推力和阻力进行重新定义, 从根本上解决了鱼游研究中推力和阻力无法区分的难题.在此基础上, 利用自适应网格下的ghost-cell浸没边界方法, 模拟了鱼类以鲹科模式在黏性流体(309 \le Re \le 14\,581)和无黏流体 (相当于雷诺数无穷大情形)中的二维自主游动.结果表明: (1) Strouhal数随雷诺数增大而减小,当雷诺数趋向于无穷时, Strouhal数趋向于0.25; (2)在所有雷诺数情况下, 推力主要来源于压力分量; 当Re<3000时, 阻力的压力分量小于黏性力分量, 而当Re>3000, 二者的关系就会反过来; (3)推进效率随着雷诺数的增大而增大,当雷诺数趋向于无穷大时, 推进效率最高可以达到70%, 说明鲹科模式适用于较高雷诺数下的游动.      

On the Reynolds-number independence of mixed convection in a vertical channel subjected to asymmetric wall temperatures with and without flow reversal

The problem of mixed convection in a vertical channel with asymmetric wall temperatures including situations of flow reversal is studied numerically. The SIMPLER algorithm with a staggered grid system is employed to solve the corresponding numerical equations formulated by the finite-volume method. A second-order upwind scheme is used to model the convective term, and a suitable grid distribution is introduced. The ranges of the parameters studied are 0 r_t 1, 1 Re 1000, and 0 Gr/Re 500.

The numerical results, with the streamwise coordinate scaled by the Reynolds number (Re), show that solutions for the velocity and temperature fields are independent of the Reynolds number when Re 50, even in the presence of flow reversal. These solutions, however, are dependent on r_t and Gr/Re. Subsequently, correlations are proposed for the bulk temperature distribution and the local Nusselt numbers along the hot wall and the cold wall.     

基于重叠网格与结构网格的圆柱绕流数值模拟

为研究重叠网格与结构网格在圆柱绕流数值模拟中的区别,以二维圆柱为例,利用有限元分析软件ANSYS 19.2中的DM与Mesh建立模型并划分重叠网格,利用ANSYS 19.2中的ICEM建立模型并划分结构网格。采用FLUENT 19.2中laminar模型模拟分析系统中的平均升力系数、平均阻力系数、斯特劳哈尔数St等流体动力特性。通过改变流体流速得到两种不同网格下各6组雷诺数Re,这6组雷诺数在60~160之间。结果表明:结构网格与重叠网格的St都随着Re的增加而增加,但相同雷诺数下重叠网格对应的St数值更大,St的增长速度更快;重叠网格与结构网格的平均升力系数与平均阻力系数随着Re的增加趋于稳定的速度都加快,但结构网格的平均升力系数与平均阻力系数趋于稳定的速度更快,且两种网格的平均升力系数与平均阻力系数趋于稳定速度的差距逐渐缩小,当Re=160时,两种网格的平均升力系数与平均阻力系数趋于稳定的速度几乎相同;当雷诺数在60~160之间时,采用重叠网格计算出来的斯特劳哈尔数比结构网格更加接近理论值;从升力功率谱密度分布曲线中可以看出,随着雷诺数的增加,两种网格下的频率逐渐变大,并且相同雷诺数下重叠网格的频率比结构网格大。     

Kolmogorov scales in a driven cavity flow

Turbulent flow in a three-dimensional driven cavity has been simulated directly by solving the Navier–Stokes equations. The results at Re=3200 and 10 000 compare well with the experimental data. Viscous dissipation rate has been calculated without making the assumption of isotropy. Near the top moving wall, the instantaneous dissipation rate is very high and also has high amplitude. Its frequency increases but amplitude decreases as one moves away from the wall and it becomes intermittent in the vortex core. The high Reynolds number assumption that dissipation is mainly due to the fluctuating velocity components is seen to be true in the present case except near the wall. The Kolmogorov length scale attains higher values in the core of the primary vortex due to low dissipation rate there. A value of 0.01 times the size of the cubic cavity is a good representative value at Re=10 000. Even though the present (84×84×84) grid cannot resolve this scale very well, it can resolve all the scales dynamically significant in the flow as seen from the velocity and dissipation spectra.     

Steady and transient liquid sphere heating in a convective gas stream

 A finite-difference scheme has been developed to solve the equations governing the laminar forced convection heat transfer around and inside a spherical fluid droplet moving steadily in another immiscible fluid for both steady and transient thermal conditions. For large values of the external flow Reynolds number (Re), results not available in the literature have been obtained for circulating droplets at intermediate and high interior-to-exterior viscosity ratios (μ^*). Detailed results over a wide range of viscosity ratio (μ^*) and for 200≤Re≤1000 are presented for the temperature profiles outside and inside the sphere, Nusselt number, the time required to attain a uniform surface temperature and the time required to reach the steady-state temperature. Results show that convective heating is dependent on the external flow Reynolds number (Re) and the interior-to-exterior viscosity ratio (μ^*) where increasing Re or decreasing μ^* result in increasing heat transfer rate convected to the liquid sphere. Received on 1 March 1999     

Mind the gap: a new insight into the tip leakage vortex using stereo-PIV

The tip leakage vortex (TLV), which develops in the clearance between the rotor and the stator of axial hydro turbines, has been studied for decades. Yet, many associated phenomena are still not understood. For instance, it remains unclear how the clearance size is related to the occurrence of cavitation in the vortex, which can lead to severe erosion. Experiments are here carried out on the influence of the clearance size on the tip vortex structure in a simplified case study. A NACA0009 hydrofoil is used as a generic blade in a water tunnel while the clearance between the blade tip and the wall is varied. The 3D velocity fields are measured using Stereo Particle Image Velocimetry (SPIV) in three planes located downstream of the hydrofoil for different values of the upstream velocity, the incidence angle and a large number of tip clearances. The influence of the flow conditions on the structure of the TLV is described through changes in the vortex intensity, core axial flow, vortex center position and wandering motion amplitude. Moreover, high-speed visualizations are used to highlight the vortex core trajectory and clearance flow alteration, turning into a wall jet as the tip clearance is reduced. The measurements clearly reveal the existence of a specific tip clearance for which the vortex strength is maximum and most prone to generating cavitation.     

Dependence of flow classification on the Reynolds number for a two-cylinder wake

This work aims to investigate the dependence of flow classification on the Reynolds number (Re) for the wake of two staggered cylinders. The Re examined ranges from 1.5×10³ to 2.0×10⁴. The pitch ratio, P^⁎=P/d examined is 1.2–6.0 (d is the cylinder diameter), and angle (α) is 0–90°, where P is the center-to-center spacing between two cylinders and α is the angle between the incident flow and the line through the cylinder centers. Two single hotwires were used to measure simultaneously the fluctuating streamwise velocities (u) in the vortex streets generated by the two cylinders. The power spectral density functions and the Strouhal numbers were then obtained from the u signals, based on which the flow structure pattern or mode could be determined. Over two hundred configurations of two staggered cylinders have been examined for each Re. It is found that Re has an appreciable effect on the dependence of the flow mode on P^⁎ and α. The observation is connected to the Re effect on the generic features of a two-cylinder wake such as flow separation, boundary layer thickness, gap flow deflection and vortex formation length.     

Simulation of peristaltic flow of chyme in small intestine for couple stress fluid

In this article couple stress fluid have been considered for the peristaltic flow of chyme in intestine. Problem under consideration have been formulated assuming that two non-periodic sinusoidal wave of different wavelength propagate with same speed c along the outer wall of the tube. Governing equations have been simplified under the assumptions of long wavelength and low Reynolds number approximation (such assumption are consistent that Re (Reynold number) is very small and long wavelength approximation also exists in the small intestine). Exact solutions have been evaluated for velocity and pressure rise. Physical behaviour of different parameter of couple stress fluid have been presented graphically for velocity, pressure rise, pressure gradient and frictional forces. The stream lines are also made against different parameters.     

颗粒在黏弹性流体扩张和收缩流中沉降的格子Boltzmann模拟分析

对于Oldroyd-B型黏弹性流体,本文应用格子Boltzmann方法(LBM),实现了流体在二维1:3扩张流道及3:1收缩流道中流动的数值模拟,获得了黏弹性流体在扩张和收缩流道中的流场分布．结合颗粒的受力和运动规则,基于点源颗粒模型,数值分析了颗粒在扩张流和收缩流中的沉降过程和特征,讨论了颗粒相对质量和起始位置以及雷诺数Re和威森伯格数Wi对颗粒沉降特征的影响．结果表明,颗粒相对质量和起始位置以及Re对颗粒沉降轨迹和落点影响较大,而Wi的影响则较小．     

MODELING THE UNSTEADY LIFT AND DRAG ON A FINITE-LENGTH CIRCULAR CYLINDER IN CROSS-FLOW

Semi-empirical models for unsteady lift and drag are developed to predict the spectral features of the unsteady forces on a finite-length, right circular cylinder in cross-flow. In general, the models consist of two parts; the spatial variation of r.m.s wall pressure on the cylinder, and the correlation lengths which describe the spatial extent of the correlation of the unsteady wall pressures. Experiments were conducted in a low noise wind tunnel as a function of cylinder diameter Reynolds number (19 200<Re<32 000) and the Strouhal number (0·05< St<3·33), to measure the statistics of the unsteady wall pressures on a model cylinder. These results are incorporated into the theoretical models, and predictions of the spectral characteristics of the lift and drag are made. The r.m.s. wall pressures on the cylindrical surface are found to have the largest amplitude near the cylinder end-cap, and on the rearward portion of the cylinder body. The high levels in these locations are attributed to the separated flow region over the end-cap. The circumferential and axial length-scales decrease exponentially with Strouhal number. Both length-scales exhibit maxima near the Strouhal shedding frequency of St=0·21. The axial length-scales are found to depend on the measurement reference location due to the three-dimensional flow and separated flow region near the end-cap. The unsteady lift and drag predictions using the models developed in this work agree well with previously measured unsteady force data measured on inertial hydrophones exposed to flow. The broadband unsteady lift is found to be greater than the broadband unsteady drag by nominally 3dB.     

Interaction theory of hypersonic laminar near-wake flow behind an adiabatic circular cylinder

The separation and shock wave formation on the aft-body of a hypersonic adiabatic circular cylinder were studied numerically using the open source software OpenFOAM. The simulations of laminar flow were performed over a range of Reynolds numbers (\(8\times 10^3 < Re < 8\times 10^4\)) at a free-stream Mach number of 5.9. Off-body viscous forces were isolated by controlling the wall boundary condition. It was observed that the off-body viscous forces play a dominant role compared to the boundary layer in displacement of the interaction onset in response to a change in Reynolds number. A modified free-interaction equation and correlation parameter has been presented which accounts for wall curvature effects on the interaction. The free-interaction equation was manipulated to isolate the contribution of the viscous–inviscid interaction to the overall pressure rise and shock formation. Using these equations coupled with high-quality simulation data, the underlying mechanisms resulting in Reynolds number dependence of the lip-shock formation were investigated. A constant value for the interaction parameter representing the part of the pressure rise due to viscous–inviscid interaction has been observed at separation over a wide range of Reynolds numbers. The effect of curvature has been shown to be the primary contributor to the Reynolds number dependence of the free-interaction mechanism at separation. The observations in this work have been discussed here to create a thorough analysis of the Reynolds number-dependent nature of the lip-shock.     

Dean instability and secondary flow structure in curved rectangular ducts

The pressure driven, fully developed turbulent flow of incompressible viscous fluid (water) in 120° curved ducts of rectangular cross-section is investigated experimentally and numerically. Three different types of curved duct (A-CL, B-SL and C-IL) with continuously varying curvature conform to blade profile as the inner and outer curvature walls to simplify and guide the impeller design of pumps. After validating the numerical method against Particle Image Velocimetry (PIV) measurements, the flow development in the ducts is analyzed in detail by Computational Fluid Dynamics (CFD) for a wide range of Reynolds numbers (Re = 2.4 × 10⁴–1.4 × 10⁵) and aspect ratios (Ar > 1.0, =1.0 and <1.0). The results clearly depict the existence of multiple Dean vortices along the duct: while the axial velocity profile is more related to an inner Dean vortex (called split base vortex), the wall pressure is more influenced by the Dean vortex attached to the inner curvature wall (called ICW Dean vortex). The induced multiple Dean vortices and the secondary flow patterns in the duct cannot be faithfully predicted by using traditional techniques. Therefore, a new criterion based on the vortex core velocities is proposed. With this approach, the effects of Re, Cr and Ar on the Dean instabilities in curved ducts are carefully studied. Decreasing Re promotes the generation of Dean vortices closer to the duct inlet, a trend that is as opposed to laminar flow. In addition, a new pair of vortices called entrainment Dean vortex occurs near the outlet of the curved duct with Ar = 1.0, which has not been previously reported in the literature.