Laser-Doppler measurements of laminar and turbulent flow in a pipe bend

Laser-Doppler measurements are reported for laminar and turbulent flow through a 90° bend of circular cross-section with mean radius of curvature equal to 2.8 times the diameter. The measurements were made in cross-stream planes 0.58 diameters upstream of the bend inlet plane, in 30, 60 and 75° planes in the bend and in planes one and six diameters downstream of the exit plane. Three sets of data were obtained: for laminar flow at Reynolds numbers of 500 and 1093 and for turbulent flow at the maximum obtainable Reynolds number of 43 000. The results show the development of strong pressure-driven secondary flows in the form of a pair of counter-rotating vortices in the streamwise direction. The strength and character of the secondary flows were found to depend on the thickness and nature of the inlet boundary layers, inlet conditions which could not be varied independently of Reynolds number. The quantitative anemometer measurements are supported by flow visualization studies. Refractive index matching at the fluid-wall interface was not used; the measurements consist, therefore, of streamwise components of mean and fluctuating velocities only, supplemented by wall pressure measurements for the turbulent flow. The displacement of the laser measurement volume due to refraction is allowed for in simple geometrical calculations. The results are intenden for use as benchmark data for calibrating flow calculation methods.     

Effect of the pattern of initial perturbations on the steady pipe flow regime

The influence of the inlet flow formation mode on the steady flow regime in a circular pipe has been investigated experimentally. For a given inlet flow formation mode the Reynolds number Re* at which the transition from laminar to turbulent steady flow occurred was determined. With decrease in the Reynolds number the difference between the resistance coefficients for laminar and turbulent flows decreases. At a Reynolds number approximately equal to 1000 the resistance coefficients calculated from the Hagen-Poiseuille formula for laminar steady flow and from the Prandtl formula for turbulent steady flow are equal. Therefore, we may assume that at Re > 1000 steady pipe flow can only be laminar and in this case it is meaningless to speak of a transition from one steady pipe flow regime to the other. The previously published results [1–9] show that the Reynolds number at which laminar goes over into turbulent steady flow decreases with increase in the intensity of the inlet pulsations. However, at the highest inlet pulsation intensities realized experimentally, turbulent flow was observed only at Reynolds numbers higher than a certain value, which in different experiments varied over the range 1900–2320 [10]. In spite of this scatter, it has been assumed that in the experiments a so-called lower critical Reynolds number was determined, such that at higher Reynolds numbers turbulent flow can be observed and at lower Reynolds numbers for any inlet perturbations only steady laminar flow can be realized. In contrast to the lower critical Reynolds number, the Re* values obtained in the present study, were determined for given (not arbitrary) inlet flow formation modes. In this study, it is experimentally shown that the Re* values depend not only on the pipe inlet pulsation intensity but also on the pulsation flow pattern. This result suggests that in the previous experiments the Re* values were determined and that their scatter is related with the different pulsation flow patterns at the pipe inlet. The experimental data so far obtained are insufficient either to determine the lower critical Reynolds number or even to assert that this number exists for a pipe at all.     

Large Eddy Simulations of convergent–divergent channel flows at moderate Reynolds numbers

The paper is focused on the study of fully turbulent channel flows, using Large Eddy Simulations (LES), in order to address the effects of adverse pressure gradient regions. Analyses of the effects of streak instabilities, which have been shown to be relevant in such regions, are extended to moderate Reynolds numbers. The work considers two different channel geometries in order to further separate influences from wall curvature, flow separation and adverse pressure gradients. Turbulent kinetic energy and Reynolds stress budgets are investigated at separation and re-attachment points. The numerical approach used in the present work is based on the incompressible Navier–Stokes equations, which are solved by a pseudo-spectral methodology for structured grids. Wall-resolved LES calculations are performed using the WALE subgrid scale model. The study shows that the streak instability mechanism persists at higher Reynolds numbers with and without wall curvature in the adverse pressure gradient regions. Moreover, the observed effects are also present regardless of the existence of flow separation regions. Finally, the study of turbulent kinetic energy budgets indicates that, independently of the flow condition, there are well-defined patterns for such turbulent properties at separation and re-attachment points.     

Large eddy simulation using a dynamic mixing length subgrid‐scale model

A novel dynamic mixing length (DML) subgrid‐scale model for large eddy simulations is proposed in this work to improve the cutoff length of the Smagorinsky model. The characteristic mixing length (or the characteristic wave number) is dynamically estimated for the subgrid‐scale fluctuation of turbulence by the cutoff wave‐number, k_c, and the dissipation wave‐number, k_d. The dissipation wave number is derived from the kinetic energy spectrum equation and the dissipation spectrum equation. To prove the promise of the DML model, this model is used to simulate the lid‐driven cubical cavity with max‐velocity‐based Reynolds numbers 8850 and 12,000, the channel flows with friction‐velocity‐based Reynolds numbers 180, 395, 590, and 950, and the turbulent flow past a square cylinder at the higher Reynolds number 21,400, respectively, compared with the Smagorinsky model and Germano et al.'s dynamic Smagorinsky model. Different numerical experiments with different Reynolds numbers show that the DML model can be used in simulations of flows with a wide range of Reynolds numbers without the occurrence of singular values. The DML model can alleviate the dissipation of the Smagorinsky model without the loss of its robustness. The DML model shows some advantages over Germano et al.'s dynamic Smagorinsky model in its high stability and simplicity of calculation because the coefficient of the DML model always stays positive. The characteristic mixing length in the DML model reflects the subgrid‐scale fluctuation of turbulence in nature and thus the characteristic mixing length has a spatial and temporal distribution in turbulent flow. Copyright © 2011 John Wiley & Sons, Ltd.     

SCALE ADAPTIVE SIMULATION OF FLOWS AROUND A CIRCULAR CYLINDER AT HIGH SUB-CRITICAL REYNOLDS NUMBER

借助γ-Re_θ转捩模型,实现了高亚临界雷诺数（Re=1.4×10⁵）下圆柱层流分离流动的尺度自适应模拟.统计平均结果表明数值计算和实验测量较为接近,尤其在圆柱后半段的分离区中,压力系数和实验符合得很好,误差主要源于分离点预测的不准确. 瞬态流动则显示,层流分离的剪切层中出现了展向不稳定,且在向下游的输运过程中不断增强,最后转捩为完全湍流. 在湍流分离模拟中,由于缺乏剪切层失稳的非定常性,SST-SAS 模型的尺度分辨能力变弱,因此在分离区以及下游尾迹中求解出的湍流尺度要明显较层流分离时大.     

High‐order discontinuous Galerkin spectral element methods for transitional and turbulent flow simulations

In this paper, we investigate the accuracy and efficiency of discontinuous Galerkin spectral method simulations of under‐resolved transitional and turbulent flows at moderate Reynolds numbers, where the accurate prediction of closely coupled laminar regions, transition and developed turbulence presents a great challenge to large eddy simulation modelling. We take full advantage of the low numerical errors and associated superior scale resolving capabilities of high‐order spectral methods by using high‐order ansatz functions up to 12th order. We employ polynomial de‐aliasing techniques to prevent instabilities arising from inexact quadrature of nonlinearities. Without the need for any additional filtering, explicit or implicit modelling, or artificial dissipation, our high‐order schemes capture the turbulent flow at the considered Reynolds number range very well. Three classical large eddy simulation benchmark problems are considered: a circular cylinder flow at Re_D=3900, a confined periodic hill flow at Re_h=2800 and the transitional flow over a SD7003 airfoil at Re_c=60,000. For all computations, the total number of degrees of freedom used for the discontinuous Galerkin spectral method simulations is chosen to be equal or considerably less than the reported data in literature. In all three cases, we achieve an equal or better match to direct numerical simulation results, compared with other schemes of lower order with explicitly or implicitly added subgrid scale models. Copyright © 2014 John Wiley & Sons, Ltd.     

Direct Numerical Simulation of Self-Similar Turbulent Boundary Layers in Adverse Pressure Gradients

Direct numerical simulations of the Navier–Stokes equations have been carried out with the objective of studying turbulent boundary layers in adverse pressure gradients. The boundary layer flows concerned are of the equilibrium type which makes the analysis simpler and the results can be compared with earlier experiments and simulations. This type of turbulent boundary layers also permits an analysis of the equation of motion to predict separation. The linear analysis based on the assumption of asymptotically high Reynolds number gives results that are not applicable to finite Reynolds number flows. A different non-linear approach is presented to obtain a useful relation between the freestream variation and other mean flow parameters. Comparison of turbulent statistics from the zero pressure gradient case and two adverse pressure gradient cases shows the development of an outer peak in the turbulent energy in agreement with experiment. The turbulent flows have also been investigated using a differential Reynolds stress model. Profiles for velocity and turbulence quantities obtained from the direct numerical simulations were used as initial data. The initial transients in the model predictions vanished rapidly. The model predictions are compared with the direct simulations and low Reynolds number effects are investigated.     

Numerical study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

Parametric study of unsteady laminar flows

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

Flow fields in a simulated axisymmetric terminal aneurysm with symmetric and asymmetric outflows

Laser-Doppler velocimetric measurements and flow visualization were performed in a glass axisymmetric aneurysm model with symmetric and asymmetric outflows through the branches. The bifurcation angle was fixed at 140°, and the Reynolds number based on the steady bulk average velocity and diameter of the affarent conduit was 500. The flow characteristics such as flow separation in the afferent conduit and flow activity inside the aneurysm for the symmetric and asymmetric outflow cases were compared in detail, and the case that is susceptible to thrombosis was identified. In addition, the onset of transition from laminar to turbulent flow inside the aneurysm was evidenced by the presence of vortex breakdown and the steep increase in the fluctuation level. Finally, the effect of pulsation on the flow pattern in the aneurysm was examined.     

Large eddy simulations of wall‐bounded flows using a simplified immersed boundary method and high‐order compact schemes

This paper presents a solution algorithm based on an immersed boundary (IB) method that can be easily implemented in high‐order codes for incompressible flows. The time integration is performed using a predictor‐corrector approach, and the projection method is used for pressure‐velocity coupling. Spatial discretization is based on compact difference schemes and is performed on half‐staggered meshes. A basic algorithm for body‐fitted meshes using the aforementioned solution method was developed by A. Tyliszczak (see article “A high‐order compact difference algorithm for half‐staggered grids for laminar and turbulent incompressible flows” in Journal of Computational Physics) and proved to be very accurate. In this paper, the formulated algorithm is adapted for use with the IB method in the framework of large eddy simulations. The IB method is implemented using its simplified variant without the interpolation (stepwise approach). The computations are performed for a laminar flow around a 2D cylinder, a turbulent flow in a channel with a wavy wall, and around a sphere. Comparisons with literature data confirm that the proposed method can be successfully applied for complex flow problems. The results are verified using the classical approach with body‐fitted meshes and show very good agreement both in laminar and turbulent regimes. The mean (velocity and turbulent kinetic energy profiles and drag coefficients) and time‐dependent (Strouhal number based on the drag coefficient) quantities are analyzed, and they agree well with reference solutions. Two subfilter models are compared, ie, the model of Vreman (see article “An eddy‐viscosity subgrid‐scale model for turbulent shear flow: algebraic theory and applications” in Physics and Fluids) and σ model (Nicoud et al, see article “Using singular values to build a subgrid‐scale model for large eddy simulations” in Physics and Fluids). The tests did not reveal evident advantages of any of these models, and from the point of view of solution accuracy, the quality of the computational meshes turned out to be much more important than the subfilter modeling.     

Numerical simulation of two‐dimensional laminar slot‐jet impingement flows confined by a parallel wall

Two‐dimensional laminar incompressible impinging slot‐jet is simulated numerically to gain insight into flow characteristics.Computations are done for vertically downward‐directed slot‐jets impinging on a plate at the bottom and confined by a parallel surface on top. The behaviour of the jet with respect to aspect ratio (AR) and Reynolds number (Re) are described in detail. The computed flow patterns for various AR (2–5) and for a range of jet‐exit Reynolds numbers (100–500) are analysed to understand the flow characteristics. The transient development of the flow is also simulated for AR = 4 and Re = 300. It is found that the reattachment length is dependent on both AR and Reynolds number for the range considered. The correlation for reattachment length is suggested. The maximum resultant velocity V_rmax and its trajectory is reported. A detailed study of horizontal velocity profile at different downstream locations is reported. It is found that the effect of Reynolds number and AR is significant to the bottom wall vorticity in the impingement and wall jet regions. Copyright © 2007 John Wiley & Sons, Ltd.     

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

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

NUMERICAL STUDIES OF TRANSITIONAL TURBULENT PULSATILE FLOW IN PIPES WITH RING–TYPE CONSTRICTIONS

Pulsatile flows in the vicinity of mechanical ring-type constrictions in pipes were studied for transitional turbulent flow with a Reynolds number (Re) of the order of 10⁴. The Womersley number (Nw) is in the range 30–50, with a corresponding Strouhal number (St) range of 0·0143–0·0398. The pulsatile flows considered are a pure sinusoidal flow, a physiological flow and an experimental pulsatile flow profile for mechanical aortic valve flow simulations. Transitional laminar and turbulent flow characteristics in an alternating manner within the pulsatile flow fields were studied numerically. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity and turbulent shear stress are smaller during the acceleration phase than during the deceleration period. Various parametric equations have been formulated through numerical experimentation to better describe the relationships between the instantaneous flow rate (Q), the pressure loss (ΔP), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max) for turbulent pulsatile flow in the vicinity of constrictions in the vascular tube. An elliptic relationship has been found to exist between the instantaneous flow rate and the instantaneous pressure gradient. Other linear and quadratic relations between various flow parameters were also obtained.     

Buoyancy modelling with incompressible SPH for laminar and turbulent flows

This work aims to model buoyant, laminar or turbulent flows, using a two‐dimensional incompressible smoothed particle hydrodynamics model with accurate wall boundary conditions. The buoyancy effects are modelled through the Boussinesq approximation coupled to a heat equation, which makes it possible to apply an incompressible algorithm to compute the pressure field from a Poisson equation. Based on our previous work [1], we extend the unified semi‐analytical wall boundary conditions to the present model. The latter is also combined to a Reynolds‐averaged Navier–Stokes approach to treat turbulent flows. The k ? ? turbulence model is used, where buoyancy is modelled through an additional term in the k ? ? equations like in mesh‐based methods. We propose a unified framework to prescribe isothermal (Dirichlet) or to impose heat flux (Neumann) wall boundary conditions in incompressible smoothed particle hydrodynamics. To illustrate this, a theoretical case is presented (laminar heated Poiseuille flow), where excellent agreement with the theoretical solution is obtained. Several benchmark cases are then proposed: a lock‐exchange flow, two laminar and one turbulent flow in differentially heated cavities, and finally a turbulent heated Poiseuille flow. Comparisons are provided with a finite volume approach using an open‐source industrial code. Copyright © 2015 John Wiley & Sons, Ltd.     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

高亚临界雷诺数圆柱绕流的尺度自适应模拟

借助γ-Re_θ转捩模型,实现了高亚临界雷诺数（Re=1.4×105）下圆柱层流分离流动的尺度自适应模拟.统计平均结果表明数值计算和实验测量较为接近,尤其在圆柱后半段的分离区中,压力系数和实验符合得很好,误差主要源于分离点预测的不准确. 瞬态流动则显示,层流分离的剪切层中出现了展向不稳定,且在向下游的输运过程中不断增强,最后转捩为完全湍流. 在湍流分离模拟中,由于缺乏剪切层失稳的非定常性,SST-SAS 模型的尺度分辨能力变弱,因此在分离区以及下游尾迹中求解出的湍流尺度要明显较层流分离时大.      

Continuation methods applied to the 2D Navier–Stokes equations at high Reynolds numbers

Nonlinearities arise in aerodynamic flows as a function of various parameters, such as angle of attack, Mach number and Reynolds number. These nonlinearities can cause the change from steady to unsteady flow or give rise to static hysteresis. Understanding these nonlinearities is important for safety validation and performance enhancement of modern aircraft. A continuation method has been developed to study nonlinear steady state solutions with respect to changes in parameters for two‐dimensional compressible turbulent flows at high Reynolds numbers. This is the first time that such flows have been analysed with this approach. Continuation methods allow the stable and unstable solutions to be traced as flow parameters are changed. Continuation has been carried out on two‐dimensional aerofoils for several parameters: angle of attack, Mach number, Reynolds number, aerofoil thickness and turbulent inflow as well as levels of dissipation applied to the models. A range of results are presented. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical analysis of axisymmetric and planar sudden expansion flows for laminar regime

In this study, the Navier–Stokes equations are solved numerically for axisymmetric and planar sudden expansion flows. The flow is considered laminar and steady state, and the fluid is incompressible. Finite difference equations are obtained using a control volume method in a non‐staggered grid arrangement, and solved by line‐by‐line TDMA technique using the SIMPLEM (SIMPLEM‐Modified) algorithm. Calculations are performed for higher expansion ratios, β, ranging from 1.5 to 10 and Reynolds numbers from 0.1 to 500. Results are presented in terms of streamlines, relative eddy intensity, location of the eddy center, and the eddy reattachment length depending on Re number and β values for both axisymmetric and planar sudden expansions. It is aimed to provide a picture of the effects of high values of expansion ratio and Reynolds number on the sudden expansion flow. As a result, it is found that the flow characteristics keep their structure for both higher expansion ratios and higher Reynolds numbers. Further, correlations were developed for the nondimensional eddy reattachment length, location of the eddy center and the relative eddy intensity, which have agreeable results to the computed results available from the literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Heat transfer enhancement in grooved channels due to flow bifurcations

The flow bifurcation scenario and heat transfer characteristics in grooved channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element methods for increasing Reynolds numbers in the laminar and transitional regimes. The Eulerian flow characteristics show a transition scenario of two Hopf bifurcations when the flow evolves from a laminar to a time-dependent periodic and then to a quasi-periodic flow. The first Hopf bifurcation occurs to a critical Reynolds number Rec₁ that is significantly lower than the critical Reynolds number for a plane-channel flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁ and m· ω₁+n·ω₂, respectively, with m and n integers. Friction factor and pumping power evaluations demonstrate that these parameters are much higher than the plane channel values. The time-average mean Nusselt number remains mostly constant in the laminar regime and continuously increases in the transitional regime. The rate of increase of this Nusselt number is higher for a quasi-periodic than for a periodic flow regime. This higher rate originates because better flow mixing develops in quasi-periodic flow regimes. The flow bifurcation scenario occurring in grooved channels is similar to the Ruelle-Takens-Newhouse transition scenario of Eulerian chaos, observed in symmetric and asymmetric wavy channels.