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

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

Flow and heat transfer performance of a mini-channel radiator with cylinder disturbed flow

Mini-channel heat sinks have relatively low Nusselt number due to small Reynolds number. For heat transfer enhancement purpose, a mini-channel radiator with cylinder disturbed flow was proposed. The disturbed flow was created by a circular cylinder placed horizontally in front of channels entrance. The performance of heat transfer and pressure drop with/without disturbed flow was studied experimentally. It was found that the friction factor of mini-channel flow was larger than that of the macro-channel flow due to larger surface roughness, and the pressure drop caused by cylinder disturbed flow was less than 5%. It also concluded that the average Nusselt number increases with augment of Reynolds and Prandtl number. The Nusselt number correlations as the function of the Reynolds and Prandtl number were given for evaluation the heat removal performance of similar heat radiators. There is an inflexion point in the empirical formulas when the channel length equals to the thermal entrance length. For the mini-channels heat radiators with disturbed flow, the inflexion Reynolds number is larger than that of without disturbed flow. Due to the flow pulsing caused by circular cylinder placed in front of channels entrance, the thermal entrance length increases. On the other hand, for both mini-channels with or without disturbed flow, the thermal resistance increases with the decrease of pressure drop.     

Transition from laminar to turbulent flow in liquid filled microtubes

The transition to turbulent flow is studied for liquids of different polarities in glass microtubes having diameters between 50 and 247 µm. The onset of transition occurs at Reynolds numbers of ~1,800–2,000, as indicated by greater-than-laminar pressure drop and micro-PIV measurements of mean velocity and rms velocity fluctuations at the centerline. Transition at anomalously low values of Reynolds number was never observed. Additionally, the results of more than 1,500 measurements of pressure drop versus flow rate confirm the macroscopic Poiseuille flow result for laminar flow resistance to within –1% systematic and ±2.5% rms random error for Reynolds numbers less than 1,800.An erratum to this article can be found at     

Dominant parameters for maximum velocity induced by body-force models for plasma actuators

This study investigates the relationship between body-force fields and maximum velocity induced in quiescent air for development of a simple body-force model of a plasma actuator. Numerical simulations are conducted with the body force near a wall. The spatial distribution and temporal variation of the body force are a Gaussian distribution and steady actuation, respectively. The dimensional analysis is performed to derive a reference velocity and Reynolds number based on the body-force distribution. It is found that the derived Reynolds number correlates well with the nondimensional maximum velocity induced in quiescent conditions when the center of the Gaussian distribution is fixed at the wall. Additionally, two flow regimes are identified in terms of the Reynolds number. Considering the variation of the center of gravity of force fields, another Reynolds number is defined by introducing a new reference length. The nondimensional maximum velocity is found to be scaled with the latter Reynolds number, i.e., the maximum induced velocity in quiescent conditions is determined from three key parameters of the force field: the total induced momentum per unit time, the height of the center of gravity, and the standard deviation from it. This scaling turns out to be applicable to existing body-force models of the plasma actuator, despite the force distributions different from the Gaussian distribution. Comparisons of velocity profiles with experimental data validate the results and show that the flow induced by a plasma actuator can be simulated with simple force distributions by adjustment of the key body-force parameters.     

Slip and induced magnetic field effects on peristaltic transport of Johnson-Segalman fluid

The peristaltic flow of a Johnson-Segalman fluid in a planar channel is investigated in an induced magnetic field with the slip condition.The symmetric nature of the flow in a channel is utilized.The velocity slip condition in terms of shear stresses is considered.The mathematical formulation is presented,and the equations are solved under long wavelength and low Reynolds number approximations.The perturbation solutions are established for the pressure,the axial velocity,the micro-rotation component,the stream function,the magnetic-force function,the axial induced magnetic field,and the current distribution across the channel.The solution expressions for small Weissenberg numbers are derived.The flow quantities of interest are sketched and analyzed.     

Crown behavior and bubble entrainment during a drop impact on a liquid film

Physical and mathematical models are established to simulate a single liquid drop impinging onto a liquid film using the coupled level set and volume of fluid method. The crown liquid sheet after impact is obtained, which coincides well with the experimental results in literatures. Influence of Weber number, Reynolds number and the dimensionless film thickness on the crown diameter and height is discussed quantitatively. Results indicate that the crown diameter is independent of the two non-dimensional numbers, while it can be increased by reducing the dimensionless film thickness. The crown height increases with the increasing of Weber number, but Reynolds number has small effect on it. Mechanism of the jet formation process is revealed by analyzing pressure distribution and velocity field in the liquid. It is found that both pressure difference in the neck region and velocity discontinuity can greatly affect the jet formation. Besides, the bubble entrainment phenomenon during a liquid drop impact on a liquid film is successfully captured with this numerical method. It is found that the increase in both impact Weber number and the drop diameter contributes to the emerging of bubble rings.     

Analysis of Laminar Flow and Forced Convection Heat Transfer in a Porous Medium

The flow of an incompressible Newtonian fluid confined in a planar geometry with different wall temperatures filled with a homogenous and isotropic porous medium is analyzed in terms of determining the unsteady state and steady state velocities, the temperature and the entropy generation rate as function of the pressure drop, the Darcy number, and the Brinkman number. The one-dimensional approximate equation in the rectangular Cartesian coordinates governing the flow of a Newtonian fluid through porous medium is derived by accounting for the order of magnitude of terms as well as accompanying approximations to the full-blown three-dimensional equations by using scaling arguments. The one-dimensional approximate energy and the entropy equations with the viscous dissipation consisting of the velocity gradient and the square of velocity are derived by following the same procedure used in the derivation of velocity expressions. The one-dimensional approximate equations for the velocity, the temperature, and the entropy generation rate are analytically solved to determine the velocity, the temperature, and the entropy distributions in the saturated porous medium as functions of the effective process parameters. It is found that the pressure drop, the Darcy number, and the Brinkman number affect the temperature distribution in the similar way, and besides the above parameters, the irreversibility distribution ratio also affects the entropy generation rate in the similar way.     

Effects of two-dimensional V-shaped grooves on turbulent channel flow

Detailed Laser Doppler velocimeter (LDV) measurements have been carried out in a turbulent rectangular channel flow with one rough wall. The roughness elements of two-dimensional spanwise 120° V-shaped grooves are periodically arranged with different depths and pitches. The Reynolds number based on the centerline velocity, and the channel half height ranges from 2,740 to 20,000. The comparisons of turbulence statistics over smooth and rough walls indicate that the present roughness leads to a significant change in the turbulence both in the inner and in the outer flow. Particularly, the distribution density of the grooves is a key parameter to evaluate the effect of roughness. The low-Reynolds-number dependence of turbulence statistics is also observed. The rough walls with the same pitch-to-depth ratio exhibit the equivalent roughness function under the corresponding Reynolds numbers. The disagreement of velocity defect profiles between smooth and rough walls challenges the defect universal law. The variations of the turbulence stresses and Reynolds shear stress decomposition in the outer layer suggest that the turbulent motions may be modified by the present grooves. The importance of sweep events for the present groove-roughened walls is reflected by the differences in relative contribution to Reynolds shear stress from each quadrant and the higher-order moments over smooth and rough walls.     

Flow of Particulate-Fluid Suspension in a Channel with Porous Walls

The problem of the dispersed particulate-fluid two-phase flow in a channel with permeable walls under the effect of the Beavers and Joseph slip boundary condition is concerned in this paper. The analytical solution has been derived for the longitude pressure difference, stream functions, and the velocity distribution with the perturbation method based on a small width to length ratio of the channel. The graphical results for pressure, velocity, and stream function are presented and the effects of geometrical coefficients, the slip parameter and the volume fraction density on the pressure variation, the streamline structure and the velocity distribution are evaluated numerically and discussed. It is shown that the sinusoidal channel, accompanied by a higher friction factor, has higher pressure drop than that of the parallel-plate channel under fully developed flow conditions due to the wall-induced curvature effect. The increment of the channel’s width to the length ratio will remarkably increase the flow rate because of the enlargement of the flow area in the channel. At low Reynolds number ranging from 0 to 65, the fluids move forward smoothly following the shape of the channel. Moreover, the slip boundary condition will notably increase the fluid velocity and the decrease of the slip parameter leads to the increment of the velocity magnitude across the channel. The fluid-phase axial velocity decreases with the increment of the volume fraction density.     

Pressure Drop for Low Reynolds-Number Flows Through Regular and Random Screens

Creeping flow through both regular and irregular screens was simulated by the lattice-Boltzmann method, and the dependence on screen porosity and Reynolds number of the pressure drop across the screen was analyzed. Regular structures were planar arrays of straight fibers or woven one-layer structures. The irregular planar structures were composed of randomly located and oriented fibers of finite length. A simple function of screen porosity based on partly numerical scaling arguments was found to describe accurately the simulated pressure drop across all regular screens. Due to their bigger surface area, the flow resistance of woven screens was found to be about 15% larger than that of regular planar screens with the same porosity. The pressure drop across irregular planar screens was found to be described by the same screen-porosity function with a slightly different ‘scaling’ exponent which thus appears to be dependent on the structure of the screen. The flow resistance of irregular structures was found to be clearly smaller than that of regular structures because of channelling of the flow through very few largest pores.     

Pressure losses in transitions between square and rectangular ducts of the same cross-sectional area

Experimental results are presented for the pressure loss in transitions between square and rectangular ducts where the two ends have the same cross-sectional area. The aspect ratios at the rectangular end ranged from 0.3 to 0.625, and the transition length from 1 to 2 times the hydraulic diameter. Reynolds numbers ranged from 50 000 to 125 000. The pressure drop may be divided into components arising from friction and velocity profile distortion. The friction component, which may be evaluated by normal pipe flow methods, accounts for the observed variation with Reynolds number. The velocity profile component increases as the aspect ratio of the rectangular end falls, and is significantly higher for rectangular to square than for square to rectangular transitions. There is an optimum length to hydraulic diameter ratio, for which the pressure loss is a minimum; it has not been found exactly, but is less than 2 and probably below 1.     

Experiments on elongational flow of dilute polymer solutions Part I: Jet reaction and excess pressure drop for the flow through small apertures

Experiments have been made with dilute polymer solutions on the reaction of jets issuing from small orifices and the excess pressure drop for orifice and capillary flows.Under the flow conditions with vortices occurring upstream of the aperture, the jet reaction is nearly zero below some mean velocity for PEO solutions and similarly zero below some generalized Reynolds number for Separan solutions. The normalized jet reactions, when they possess positive values, are correlated with the generalized Reynolds number irrespective of the aperture diameters for both kinds of solution.In most cases, the pressure is higher than in the corresponding water flow, but for some flows with no vortex it is lower. For the vortex flow of PEO solutions the normalized excess pressure drop is inversely proportional to the Reynolds number for both orifices and capillaries, while for Separan solutions this quantity is not correlated with the generalized Reynolds number for orifice flow but is correlated with it for capillary flow.     

Laminar flow development and heat transfer in converging plane-walled channels

The characteristics of flow development and heat transfer in converging plane-walled channels are studied by the finite difference method. The velocity and temperature profiles in both angular and radial directions, the average Nusselt number and the pressure drop are calculated for three different taper angles. The results show that the transport process is governed by three parameters: the inlet Reynolds number, the Péclet number and the taper angle. The increase of the taper angle yields an increase of the Nusselt number and a decrease of the pressure drop.     

Laminar flow through channels with porous walls and with an applied transverse magnetic field

Summary The problem of two-dimensional steady laminar flow of a viscous incompressible and electrically conducting fluid through a channel with two equally porous walls in the presence of a transverse magnetic field has been extended to include all values of Hartmann number and small suction velocity at the walls. Expressions for the velocity components, the pressure and the wall friction in terms of the Hartmann number and the suction Reynolds number are given. It is found that the pressure drop in the major flow direction and the wall friction decrease with the increase in suction and increase with the increase in the strength of the magnetic field.     

Turbulent Drag Reduction by Uniform Blowing Over a Two-dimensional Roughness

Direct numerical simulation (DNS) of turbulent channel flow over a two-dimensional irregular rough wall with uniform blowing (UB) was performed. The main objective is to investigate the drag reduction effectiveness of UB on a rough-wall turbulent boundary layer toward its practical application. The DNS was performed under a constant flow rate at the bulk Reynolds number values of 5600 and 14000, which correspond to the friction Reynolds numbers of about 180 and 400 in the smooth-wall case, respectively. Based upon the decomposition of drag into the friction and pressure contributions, the present flow is considered to belong to the transitionally-rough regime. Unlike recent experimental results, it turns out that the drag reduction effect of UB on the present two-dimensional rough wall is similar to that for a smooth wall. The friction drag is reduced similarly to the smooth-wall case by the displacement of the mean velocity profile. Besides, the pressure drag, which does not exist in the smooth-wall case, is also reduced; namely, UB makes the rough wall aerodynamically smoother. Examination of turbulence statistics suggests that the effects of roughness and UB are relatively independent to each other in the outer layer, which suggests that Stevenson’s formula can be modified so as to account for the roughness effect by simply adding the roughness function term.     

Large Eddy Simulation of boundary layer transition over an isolated ramp-type micro roughness element

Boundary layer transition over an isolated surface roughness element is investigated by means of numerical simulation. Large Eddy Simulation (LES) flow-modeling approach is employed to study flow characteristics and transition phenomenon past a roughness element immersed within an incoming developing boundary layer, at a height-based Reynolds number of 1170. LES numerical results are compared to experimental data from literature showing the time-averaged velocity distribution, the velocity fluctuation statistics and the instantaneous flow topology.Despite slight difference in the intensity of streamwise velocity fluctuations, the present LES results and experimental data show very good agreement. The mean flow visualization shows streamwise counter-rotating vortices pairs formation downstream of the obstacle. The primary pair induces an upwash motion and a momentum deficit that creates a Kelvin-Helmholtz type flow instability. The instantaneous flow topology reveals the formation of coherent K-H vortices downstream that produce turbulent fluctuations in the wake of the roughness element. These vortices are streched and lifted up when moving downstream. The velocity fluctuations results show that the onset of the turbulence is dominated by the energy transfer of large-scale vortices.     

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.      

平屋盖风压分布的数值模拟

基于Reynolds时均N-S方程和RSM模型对平屋面的风压分布进行了数值模拟,在此基础上系统研究了风向角、跨高比、地面粗糙度、风速等因素对屋面风压分布的影响,探讨了结构周围流场的绕流特性,最后根据屋面的结构形式及风压分布特点将屋面进行分区,给出了屋面在不同风向角下的分区风载体型系数以供工程设计参考。     

Friction factor improved correlations for laminar and turbulent gas–liquid flow in horizontal pipelines

We develop improved correlations for two-phase flow friction factor that consider the effect of the relative velocity of the phases, based on a database that includes 2560 gas–liquid flow experiments in horizontal pipes. The database includes a wide range of operational conditions and fluid properties for two-phase friction factor correlations. We classify the experiments by liquid holdup ranges to obtain composite analytical expressions for two-phase friction factor vs. the Reynolds number by fitting logistic dose curves to the experimental data with. We compute the liquid holdup values used to classify the experimental data using correlations proposed previously. The Reynolds number is based on the mixture velocity and the liquid kinematic viscosity. The Fanning friction factor for gas–liquid is defined in term of the mixture velocity and density. Additionally, we sort the experimental data by flow regime and obtain the two-phase friction factor improved correlations for dispersed bubble, slug, stratified and annular flow for different holdup ranges. We report error estimates for the predicted vs. measured friction factor together with standard deviation for each correlation. The accuracy of the correlations developed in this study is compared with that of other 21 correlations and models widely available in the specialized literature. Since different authors use different definitions for friction factors and Reynolds numbers, we present comparisons of the predicted pressure drop for each and every data point in the database. In most cases our correlations predict the pressure drop with much greater accuracy than those presented by previous authors.     

激波诱导高速气流中液滴的初期变形

基于激波管平台和高速摄影方法对平面激波诱导高速气流中液滴的早期变形现象进行实验研究。研究发现在相近的We数或Re数下,实验参数的改变可导致液滴形态发展出现显著差异。这种差异主要体现在背风面的脊状环形突起、褶皱区以及后驻点区的凹凸形态。对刚性圆球外流的数值模拟显示,液滴变形早期形态与外流场结构和表面气动力分布之间存在鲜明的对应关系。最后采用简化理论推导出一组估测液滴早期变形的表达式。将数值模拟所得气动力数据代入计算发现:导致液滴变形的主要驱动力是液滴表面不均匀压力的挤压效应,而不是界面剪切摩擦所引起的切向流动堆积效应,前者高出后者约2个数量级;此外,采用压力作用理论计算所得液滴外形在主要变形特征和变形量级上均可与实验图像很好地吻合。