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.     

Turbulent Air-Flow Measurement with the Aid of 3-D Particle Tracking Velocimetry in a Curved Square Bend

A three-dimensional particle tracking velocimeter (3-D PTV) was applied to air-flow measurement in a strongly curved U-bend of a square cross-section. He-filled neutral-buoyant soap bubbles were employed as a flow tracer, and turbulent statistics including all Reynolds stress components were measured. The pressure-induced secondary flow, of which magnitude reached about 30% of the bulk mean velocity, was observed. The present experimental result is mostly in good agreement with the LDA data at higher bulk-mean Reynolds number taken by Chang et al. The effect of the secondary flow on the production mechanism of turbulent kinetic energy as well as on the distributions of the invariants of stress anisotropy tensor was examined in detail. This revised version was published online in July 2006 with corrections to the Cover Date.     

Finite element Galerkin Approach for a Computational Study of Arterial Flow

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅａｓｕｒｅｍｅｎｔｏｆｓｈｅａｒｓｔｒｅｓｓｉｎａｒｔｅｒｉａｌｆｌｏｗｈａｓｂｅｅｎｂｅｓｔｏｗｅｄｂｙｎａｔｕｒｅｂｅｃａｕｓｅｏｆｔｈｅｐｏｓｓｉｂｌｅｒｅｌａｔｉｏｎｂｅｔｗｅｅｎｔｈｅｓｈｅａｒｓｔｒｅｓｓａｔｔｈｅｗａｌｌａｎｄｅｘｉｓｔｅｎｃｅｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｓ.Ｔｈｅｆｌｏｗｐｈｅｎｏｍｅｎａｉｎｌａｒｇｅａｎｄｍｅｄｉｕｍｓｉｚｅｄａｒｔｅｒｉｅｓｉｎｆｌｕｅｎｃｅｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆａｔｈｅｒｏｓｃｌｅｒｏｓｉｃｄ…     

Modeling of convective structures in a thin layer of silicone oil in air flow under heating

Numerical simulation was performed to study convective structures in a thin silicone oil layer heated from below, whose free surface is exposed to air flow generating drift flow. The basic equations are transformation to a form suitable for spectral methods. The steady flow velocity profile obtained in a laboratory experiment is calculated. It is shown that increasing the Reynolds number leads to the transition from polygonal convective cells to longitudinal rolls (elongated along the flow). The dependence of the transition Reynolds number on the temperature on the lower boundary of the layer is obtained. The calculation results are compared with experimental data.     

带扰流小槽道内单相流动阻力特性实验

采用水作为工质, 实验研究入口或出口端加入圆柱扰流的不同高宽比多槽道散热器的压降特性; 结果表明: 单位长度压降随雷诺数成线性关系. 提出了雷诺数、水力学直径和槽道高宽比的拟合准则; 拟合公式在实验数据误差范围内很好反应类似系统的流阻特性. 对于线切割小槽道, 其摩阻系数和雷诺数成反比, 且较圆管理论值偏大; 对线切割槽道阻力特性分析发现, 表面粗糙度是其最主要的影响因素, 而扰流对其影响较小; 在槽道前部设置扰流柱可以增强换热.      

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

后台阶分离流动中大涡结构演变的数值模拟

本文对后台阶分离流动中涡结构的演变进行了大涡模拟，研究了流场结构的变化规律。详细讨论了随着雷诺数的增加流场结构的典型特征的变化规律，指出流场中的涡结构随着雷诺数的增大变得十分复杂和丰富，回流区的数目、大小及其出现的位置也显著地不同。这些结果与已有的一些实验值和流场显示结果是吻合的。在此基础上，进一步研究了高雷诺数时流场中大尺度涡结构的瞬时发展和演化过程，展示了其中大涡的产生、追随、吸引、合并和破碎等过程。对于高雷诺数情况，对大涡模拟得到的数值结果进行了统计，得到的时均速度分布以及台阶后方的回流区长度与现有的其他实验结果符合得很好。本研究是针对后台阶分离流动深入开展湍流控制以及两相流动研究的基础。     

Low-diffusion preconditioning scheme for numerical simulation of low-speed flows past airfoil

The preconditioning technique can address the stiffness of a low Mach number flow, while its stability is poor. Based on the conventional preconditioning method of Roe's scheme, a low-diffusion scheme is proposed. An adjustable parameter is introduced to control numerical dissipation, especially over the dissipation in the boundary layer and extremely in a low speed region. Numerical simulations of the low Mach number and low Reynolds number flows past a cylinder and the low Mach number and high Reynolds number flows past NACA0012 and S809 airfoils are performed to validate the new scheme. Results of the three tests well agree with experimental data, showing the applicability of the proposed scheme to low Mach number flow simulations.     

后台阶流动再附着过程的大涡模拟研究

应用自主开发的大涡模拟程序数值模拟研究了后台阶流动中再附着过程的演变。在流动几何参数不变情况下，给出了再附长度随雷诺数的变化规律，并与实验进行了比较，二者相符得比较好。在此基础上，给出了三种典型雷诺数下，后台阶流动的回流区特征。在湍流情况下，研究了突扩比对再附长度的影响，与实验结果吻合的比较好。详细讨论了湍流情况下大涡拟序结构的瞬时再附着过程。 这些研究结果对具有再附着现象的流动结构的工程应用具有指导意义。     

Composite power law holdup correlations in horizontal pipes

A wide range of experimental holdup data, from different sources, are analyzed based on a theoretical model proposed in this work to evaluate the holdup in horizontal pipes. 2276 gas–liquid flow experiments in horizontal pipelines with a wide range of operational conditions and fluid properties are included in the database. The experiments are classified by mixture Reynolds number ranges and composite analytical expressions for the relationship between the liquid holdup and no-slip liquid holdup vs. the gas–liquid volumetric flow rate are obtained by fitting the data with logistic dose curves. The Reynolds number appropriate to classify the experimental data for gas–liquid flows in horizontal pipes is based on the mixture velocity and the liquid kinematic viscosity. Composite power law holdup correlations for flows sorted by flow pattern are also obtained. Error estimates for the predicted vs. measured holdup correlations together with standard deviation for each correlation are presented. The accuracy of the correlations developed in this study is compared with the accuracy of 26 previous correlations and models in the literature. Our correlations predict the liquid holdup in horizontal pipes with much greater accuracy than those presented by previous authors.     

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Slip ratio in dispersed viscous oil-water pipe flow

In this article, dispersed flow of viscous oil and water is investigated. The experimental work was performed in a 26.2-mm-i.d. 12-m-long horizontal glass pipe using water and oil (viscosity of 100 mPa s and density of 860 kg/m³) as test fluids. High-speed video recording and a new wire-mesh sensor based on capacitance (permittivity) measurements were used to characterize the flow. Furthermore, holdup data were obtained using quick-closing-valves technique (QCV). An interesting finding was the oil-water slip ratio greater than one for dispersed flow at high Reynolds number. Chordal phase fraction distribution diagrams and images of the holdup distribution over the pipe cross-section obtained via wire-mesh sensor indicated a significant amount of water near to the pipe wall for the three different dispersed flow patterns identified in this study: oil-in-water homogeneous dispersion (o/w H), oil-in-water non-homogeneous dispersion (o/w NH) and Dual continuous (Do/w & Dw/o). The phase slip might be explained by the existence of a water film surrounding the homogeneous mixture of oil-in-water in a hidrofilic-oilfobic pipe.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

A fractal analysis of subcooled flow boiling heat transfer

A fractal model for the subcooled flow boiling heat transfer is proposed in this paper. The analytical expressions for the subcooled flow boiling heat transfer are derived based on the fractal distribution of nucleation sites on boiling surfaces. The proposed fractal model for the subcooled flow boiling heat transfer is found to be a function of wall superheat, liquid subcooling, bulk velocity of fluid (or Reynolds number), fractal dimension, the minimum and maximum active cavity size, the contact angle and physical properties of fluid. No additional/new empirical constant is introduced, and the proposed model contains less empirical constants than the conventional models. The proposed model takes into account all the possible mechanisms for subcooled flow boiling heat transfer. The model predictions are compared with the existing experimental data, and fair agreement between the model predictions and experimental data is found for different bulk flow rates.     

Laser Doppler anemometry measurements of laminar flow in helical pipes

Three-dimensional laser Doppler anemometry measurements are performed on developed laminar flow in three helical pipes. The experimental observations are compared to results of numerical calculations employing the fully elliptic numerical method. Good agreement is found between measured data and numerical results. The three helical pipes, with curvature ratios of 0.0734 and 0.1374 and non-dimensional pitches of 0.0793 and 0.193, are adopted to study the effects of curvature and pitch on laminar flow in the experimental approach. The range of Reynolds numbers is 500–2000 to ensure laminar flow in the entire helical pipe. Both the profile shapes of the normal components of the secondary flow and those of the axial flow along the same centerline present not only similar patterns but also similar change when pitch, curvature ratio, and Reynolds number vary. The results demonstrate comprehensive relationships between the axial flow and the secondary flow.     

垂直湍流液-固流中大颗粒的相对速度

通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式．用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数．实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。     

An investigation of adiabatic spiral vortex flow in wide annular gaps by visualisation and digital analysis

This paper, the third in a series describing experimental investigations into spiral vortex flow, presents visual evidence illustrating adiabatic transition modes in a wide gap of radius ratio 0.848. Also, power spectra, relating to velocity fluctuations, are found to be comparable for two working fluids, oil and air. Good agreement has been found between these results and those published by the authors in which diabatic transition modes were related to the heat transfer characteristics of the flow at various axial Reynolds numbers     

A computational study of gas flow in a De‐Laval micronozzle at different throat diameters

A numerical study has been carried out to investigate the gas flows in a micronozzle using a continuum model under both slip and no‐slip boundary conditions. The governing equations were solved with a finite volume method. The numerical model was validated with available experimental data. Numerical results of exit thrust showed good agreement with experimental data except at very low Reynolds numbers. For parametric studies on the effect of geometric scaling, the nozzle throat diameter was varied from 10 to 0.1 mm, whereas throat Reynolds number was varied from 5 to 2000. A correlation has also been developed to calculate the specific impulse at specified throat diameter and Reynolds number. The effect of different gases on the specific impulse of the nozzle, such as helium, nitrogen, argon and carbon dioxide, was also examined. Copyright © 2008 John Wiley & Sons, Ltd.     

The circular cylinder in simple shear at moderate Reynolds numbers: An experimental study

The flow around a neutrally buoyant circular cylinder in simple shear was studied using particle-image velocimetry. Steady experimental results are presented for both torque-free and fixed cylinders at shear-based Reynolds numbers between 8 and 25 at two different flow confinement values. These data, which are to our knowledge the first experimental results for this Reynolds number range, show that the normalized rotation rate decays much more slowly with Reynolds number than previous theoretical predictions for unconfined flow, with the rate of decay decreasing as confinement increases. A new simple model is proposed to explain these effects. Given physical and technological limitations upon steady and unconfined flow experiments and simulations, the effect of confinement on a cylinder in simple shear merits further study. Received: 17 February 2000/Accepted: 22 June 2000