Weissenberg effect and its dependence upon the experimental geometry

The flow of a second-order fluid with a free surface between two coaxially mounted cylinders of finite length, the inner one of which is rotating, is being studied. In the case of slow flow and small shear rates the flow can be divided into a primary flow in the plane perpendicular to the axis of rotation and a secondary flow in the meridional plane. These flow components are numerically calculated and the results are compared with the analytical results for the semi-infinite cylinder approximation. The influence of the finiteness of the cylinders (end effect) upon the free surface deformation is analysed. The numerical results for the secondary flow are compared with results obtained by flow visualisation.     

基于PIV测量的超声波流量计内流场特性研究

采用PIV(Particle Image Velocimetry)测量手段,考察了小口径超声波流量计的流动特性。首先针对前端安装直管段时,不同流量条件下的流场特性建立基本认识,实验结果表明,在低流量条件下,流量计内流场存在明显的不稳定演变和非定常流动特征。进一步以上游前端安装球阀为典型案例,考察了安装条件对超声波流量计响应特性和测量偏差的影响。结合直管段的实验观测结果,发现此种结构超声波流量计的适应性与其流场非定常性的关系具有很好的一致性,即流场结构稳定则适应性强。此外,综合多参数的实验结果表明,雷诺数是判断小口径超声波流量计测量准确性的重要无量纲参数。     

Experimental investigation of the flow characteristics within a vortex tube with different configurations

The energy separation in a vortex tube is a combined result of different factors and its explanation remains debatable. As a classical fluid mechanics phenomenon, understanding of the complex helical flow mechanism within a vortex tube is a necessary foundation. The small scale of an industrial vortex tube and the extremely complex flow conditions are the two main challenges in obtaining the internal flow properties. This paper reports the results of an experimental investigation on the flow behaviour within a confined cylindrical system having different configurations corresponding to the actual flow field in a vortex tube at different conditions. Transparent devices were used to enable flow visualisation and Particle Image Velocimetry (PIV) measurement. The results of the flow visualisation and PIV experiments show that a precessing vortex core is significant only in a specific range of swirling strength. A good agreement between the observed flow characteristics and previously published results was observed.     

固体火箭燃气射流驱动液柱过程的CFD分析

固体火箭燃气射流驱动液柱过程会产生一个复杂的非稳态多相流场,为了研究液柱对固体火箭发动机工作过程中射流流场的降温效果,并揭示燃气冲击液柱的流动演化和气水之间的相互作用,利用FLUENT软件中耦合了液态水汽化方程的VOF多相流计算模型对燃气与液柱之间的耦合流动及相变过程进行了数值模拟,并与无液柱情况下射流流场的计算结果进行了对比分析。计算结果表明,当有液柱平衡体时射流流场中的压力、温度、速度波动幅度均减小,减弱了射流流场中的湍流脉动强度;液柱与燃气之间的汽化以及液柱的阻碍作用减小了射流流场的轴向发展位移,尾管后的完全发展射流流场核心区域内的压力峰值降低了0.9 MPa,温度峰值降低了503 K,速度峰值降低了291 m/s,验证了实验中液柱对燃气射流流场的降温效果。     

Forces on oscillating bodies in viscous fluid

This paper describes a method for determining the fluid forces on oscillating bodies in viscous fluid when the corresponding flow problem has been solved using the finite element method. These forces are characterized by the concept of added mass, added damping and added force. Numerical results are obtained for several example body shapes. Comparison is made with exact analytical results and other finite element results for the limiting cases of Stoke's flow and inviscid flow, and good agreement is obtained. The results for finite values of the body amplitude parameter β show the appearance of added force from the steady streaming component of the flow for asymmetric bodies. Results are also obtained for the associated flow where the fluid remote from a fixed body is oscillating.     

气固耦合振动叶栅非定常流动分析研究

引入结构动力学方程建立了二维N-S/结构振动耦合方程组,采用双时间法建立了气固耦合方程组的非定常数值求解体系,研究了叶栅间的二维非定常粘性流动及叶栅振动特性。对两种叶型分别计算了不同折合振动频率下的流场,振动叶栅位移随时间变化的曲线表明,采用气固耦合得到的叶栅振动频率与非耦合自振频率相比均有所下降;振动位移-时间曲线在不同振动折合频率下有显著差别。在气固耦合情况下叶栅振动规律及其稳定性与非耦合情形差异较大,因此研究叶栅振动稳定性应当考虑气动/结构的耦合。     

Thermocapillary flow without return flow–linear flow

The experimental realization of thermocapillary flow without return flow is reported. This type of flow (linear flow) was proposed and analyzed theoretically by Smith and Davis (J. Fluid Mech., 132:119–144, 1983). We suppressed the return flow by providing channels and side channels with lower flow resistance compared to that of the return flow. Cooling the layer with linear flow from above results in the Marangoni instability of longitudinal rolls as the most dangerous mode. Strong linear flow stabilizes the system against longitudinal rolls. We report preliminary results on the threshold and on the wavelength of the longitudinal rolls.     

Prediction of the slug-to-churn flow transition in vertical two-phase flow

An assessment is made of the various viewpoints on the slug-to-churn flow transition in vertical upward flow in the light of recent experimental results obtained at Harwell Laboratory. It is found that the flooding model of McQuillan & Whalley and the bubble entrainment model of Barnea & Brauner give satisfactory results at low and high liquid flow rates, respectively. An improved model for flooding, which takes account of the effect of the falling film, has been proposed. It is shown that this new model is in good agreement with experimental results at both low and high liquid flow rates.     

Correlations of flow maldistribution parameters in an air cooled heat exchanger

The present paper provides correlations of flow maldistribution parameters in air‐cooled heat exchangers. The flow field in the inlet header was obtained through the numerical solution of the governing partial differential equations including the conservation equations of mass and momentum in addition to the equations of the turbulence model. The results were obtained for different number of nozzles of 2–4, different inlet flow velocities of 1–2.5m/s and different nozzle geometries in addition to incorporation of a second header. The results are presented in terms of mass flow rate distributions in the tubes of the heat exchanger and their standard deviations. The results indicate that the inlet flow velocity has insignificant influence on maldistribution while the nozzle geometry shape has a slight effect. Also, the results indicate that reducing the nozzle diameter results in an increase in the flow maldistribution. A 25% increase is obtained in the standard deviation as a result of decreasing the diameter by 25%. Increasing the number of nozzles has a significant influence on the maldistribution. A reduction of 62.5% in the standard deviation of the mass flow rate inside the tubes is achieved by increasing the number of nozzles from 2 to 4. The results indicate that incorporating a second header results in a significant reduction in the flow maldistribution. A 50% decrease in the standard deviation is achieved as a result of incorporation of a second header of seven holes. It is also found that the hole‐diameter distribution at the exit of the second header has a slight influence on the flow maldistribution. Correlations of the flow maldistribution in terms of the investigated parameters are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

表面粗糙度对微细管内气体流动特性的影响

采用了表面粗糙度粘性系数模型对微细管内的气体流动进行数值模拟，以研究微管内壁表面粗糙度对微管内气体流动的影响。运用本文改进的表面粗糙度粘性系数模型，数值模拟与实验数据十分吻合。计算结果表明，进出口压力一定时，表面粗糙度对流场的压力、密度及温度分布的影响不大，但是对速度场的影响十分显著，表面粗糙度使气体流动速度减小，并使壁面附近的速度梯度减小，从而使通过管道的气体质量流量减小，在微管内的气体流动中，表面粗糙度的影响是不能被忽略的。     

The role of shear and elongation in the flow of solutions of semi-flexible polymers through porous media

The rheological behavior of hydrophobically modified hydroxyethyl cellulose (HMHEC) and xanthan gum solutions has been characterized in simple shear flow, opposed-jets flow, and flow through porous media. Both polymers exhibit shear thinning in simple shear flow and apparent shear thinning in flow through porous media. Analysis of the results shows there is a direct correspondence between shear viscosities determined in simple shear experiments and apparent viscosities in porous media flow at relatively low shear rates. At high shear rates the extensional component of the flow in porous media appreciably increases the apparent viscosity over the simple shear values. This increase is shown to correlate with results obtained in opposed-jets experiments, and is attributed to formation of transient entanglements.     

Numerical simulation of flow in Hartmann resonance tube and flow in ultrasonic gas atomizer

The gas flow in the Hartmann resonance tube is numerically investigated by the finite volume method based on the Roe solver.The oscillation of the flow is studied with the presence of a needle actuator set along the nozzle axis.Numerical results agree well with the theoretical and experimental results available.Numerical results indicate that the resonance mode of the resonance tube will switch by means of removing or adding the actuator.The gas flow in the ultrasonic gas atomization (USGA) nozzle is also studied by the same numerical methods.Oscillation caused by the Hartmann resonance tube structure,coupled with a secondary resonator,in the USGA nozzle is investigated.Effects of the variation of parameters on the oscillation are studied.The mechanism of the transition of subsonic flow to supersonic flow in the USGA nozzle is also discussed based on numerical results.     

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

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

Mass flow rate measurements in gas micro flows

The main objective of this experimental investigation on the gas flow slip regime is to measure the mass flow rate in isothermal steady flows through cylindrical micro tubes. Two technical procedures devoted to mass flow rate measurements are compared, and the measured values are also compared with the results yielded by different approximated analytical solutions of the gas dynamics continuum equations. Satisfactory results are obtained and the way is clearly opened to measuring mass flow rates for higher Knudsen numbers, over all the micro flow transitional regime.     

Three‐dimensional simulation of planar contraction viscoelastic flow by penalty finite element method

The planar contraction flow is a benchmark problem for the numerical investigation of viscoelastic flow. The mathematical model of three‐dimensional viscoelastic fluids flow is established and the numerical simulation of its planar contraction flow is conducted by using the penalty finite element method with a differential Phan‐Thien–Tanner constitutive model. The discrete elastic viscous split stress formulation in cooperating with the inconsistent streamline upwind scheme is employed to improve the computation stability. The distributions of velocity and stress obtained by simulation are compared with that of Quinzani's experimental results detected by laser–doppler velocimetry and flow‐induced birefringence technologies. It shows that the numerical results agree well with the experimental results. The numerical methods proposed in the study can be well used to predict complex flow patterns of viscoelastic fluids. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

小型气驱动式U形振荡水槽

本文介绍气驱动式U形振荡水楷(U形管)的设计特点、流场品质和初步实验结果。在结构设计上采用了理论计算的弯管型线、步进电机带动蝶阀的控制系统和吸气驱动等较先进的技术措施,使得该设备具有流场品质好、结构简单、调节方便等优点,是开展振荡流实验研究的基本设备。     

Some Aspects of the Lift Force on a Spherical Bubble

Several situations in which a spherical bubble experiences a lift force are examined, especially through the use of computational results obtained by solving the full Navier–Stokes equations. The lift force is computed over a wide range of Reynolds number for the case of pure shear flow, pure strain and solid body rotation. Using these results, the validity of asymptotic solutions derived in the limit of low Reynolds number or inviscid flow is discussed. A general expression of the lift force valid for low to moderate shears is proposed. It is shown that for such shears, the lift force in a complex flow can be predicted by superposing the results obtained in pure strain flow and solid body rotation flow. Finally, the interaction force experienced by two bubbles rising side-by-side is studied. The computational results reveal that, at variance with the predictions of potential theory, the sign of this force changes when the Reynolds number or the separation distance between the bubbles decreases below a critical value. All these results are discussed in terms of vorticity. The respective role played by the vorticity generated at the bubble surface and by the one that is eventually present in the unperturbed flow is emphasized.     

On the Electric Potential Induced by a Homogeneous Magnetic Field in a Turbulent Pipe Flow

In this paper we study a turbulent pipe flow of a weakly electrical conducting fluid subjected to a homogeneous magnetic field which is applied perpendicular to the flow. This configuration forms the basis of a so-called electromagnetic induction flow meter. When the Hartmann number is small so that modification of flow by the Lorenz force can be neglected, the influence of the magnetic field results only in a spatially and temporally varying electric potential. The magnitude of the potential difference across the pipe is then proportional to the flow rate and this constitutes the principle of the flow meter. In this study the flow and electric potential are computed with help of a numerical flow simulation called Large-Eddy Simulation (LES) to which we have added an equation for the electrical potential. The results of the LES have been compared with experiments in which the electric potential is measured as a function of time at several positions on the circumference of the pipe. Both the experimental and numerical results for the mean potential at the pipe wall agree very well with an exact solution that can be obtained in this particular case of a homogeneous magnetic field. Furthermore, it is found that fluctuations in the electric potential due to the turbulence, are small compared to the velocity fluctuations. Based on the results we conclude that electrical-magnetic effects in pipe flow can be accurately computed with LES.     

Pressure-driven and electroosmotic non-Newtonian flows through microporous media via lattice Boltzmann method

The lattice Boltzmann method is developed to simulate the pressure-driven flow and electroosmotic flow of non-Newtonian fluids in porous media based on the representative elementary volume scale. The flow through porous media was simulated by including the porosity into the equilibrium distribution function and adding a non-Newtonian force term to the evolution equation. The non-Newtonian behavior is considered based on the Herschel–Bulkley model. The velocity results for pressure-driven non-Newtonian flow agree well with the analytical solutions. For the electroosmotic flow, the influences of porosity, solid particle diameter, power law exponent, yield stress and electric parameters are investigated. The results demonstrate that the present lattice Boltzmann model is capable of modeling non-Newtonian flow through porous media.