七孔探针可压缩流场测量研究

介绍了七孔探针用于亚音速可压缩流的标定方法。作为一种可以同时获得流动速度大小、流动偏角、总压和静压的气动测量装置,七孔探针被广泛应用于各种流动测量,包括可压缩流动。但是它的校准过程周期很长,代价昂贵,影响了探针的推广。本文以数值计算为手段,对七孔探针进行亚音速可压缩流校准与测量的研究。结果表明,其校准拟合精度流动角为2%,内外区的总静压相对标准偏差都没有超过3%,高于相同状态下的实验校准精度。在实际应用中,本方法用于指导传统实验标定方法,可以节约大量的标定时间和成本,使七孔探针在亚音速可压缩流的测量变得简单可行。     

The inlet effect on the drag factor of a sphere in a tube

The creeping motion Ground a sphere situated axisymmetrically near the entrance of asemi-infinite circular cylindrical tube is analyzed using infinite series solutions for thevelocity components. pressure and the stream function. Truncating the infinite series. thecorresponding coefficients in the series are determined by a collocation technique. The dragfactor and the stress distribution on the surface of the sphere are calculated for the sphere inmotion in quiescent fluid and for the flow with uniform velocity at the entrance past a rigidlyheld sphere. The results indicate that a sphere near the entrance which has a uniformentrance velocity profile will suffer larger drag than that in infinite tube.Theconvergence of the collocation technique is tested by numerical calculation. It is shown thatthe technique has good convergence properties.     

Numerical simulation of laminar flow past a circular cylinder with slip conditions

The no‐slip condition is an assumption that cannot be derived from first principles and a growing number of literatures replace the no‐slip condition with partial‐slip condition, or Navier‐slip condition. In this study, the influence of partial‐slip boundary conditions on the laminar flow properties past a circular cylinder was examined. Shallow‐water equations are solved by using the finite element method accommodating SU/PG scheme. Four Reynolds numbers (20, 40, 80, and 100) and six slip lengths were considered in the numerical simulation to investigate the effects of slip length and Reynolds number on characteristic parameters such as wall vorticity, drag coefficient, separation angle, wake length, velocity distributions on and behind the cylinder, lift coefficient, and Strouhal number. The simulation results revealed that as the slip length increases, the drag coefficient decreases since the frictional component of drag is reduced, and the shear layer developed along the cylinder surface tends to push the separation point away toward the rear stagnation point so that it has larger separation angle than that of the no‐slip condition. The length of the wake bubble zone was shortened by the combined effects of the reduced wall vorticity and wall shear stress which caused a shift of the reattachment point closer to the cylinder. The frequency of the asymmetrical vortex formation with partial slip velocity was increased due to the intrinsic inertial effect of the Navier‐slip condition. Copyright © 2011 John Wiley & Sons, Ltd.     

跨超声速风洞扩散段流动分离控制设计参数研究

为抑制跨超声速风洞扩散段的分离,提出了一种较为完备的设计方法。由于影响扩散段性能的参数较多,完全通过试验方法进行设计的成本过高,该方法通过数值模拟,结合适当的边界条件,详细描述了扩散段角度、分流锥角度与长度、孔板开孔率对扩散段性能的影响;从数值模拟的结果可以看出,孔板开孔率和扩开角对扩散段性能有显著影响,通过比较得出较为合理的参数匹配,提高了扩散段的防分离性能,并改善了出口气流质量。数值结果与试验结果结论一致,表明本文所用的方法用于扩散段气动设计是可行的,为数值模拟方法应用于风洞部段气动设计创造了一定的条件。     

明渠交汇口分离区三维几何特性的大涡模拟研究

为更准确地把握交汇角对分离区三维几何特性的影响,建立了不同角度的交汇水槽模型并进行数值模拟.采用大涡模拟(LES)方法求解交汇区的湍流流场,并基于平衡层模型的Werner壁面函数法处理近壁区流场.模拟所得垂向流速分布及分离区尺寸等结果与实测资料吻合程度较高.以90°交汇水槽为例较详尽地分析了分离区的三维几何特性,并从流...     

Flow past a sphere in density-stratified fluid

Stratified flow past a three-dimensional obstacle such as a sphere has been a long-lasting subject of geophysical, environmental and engineering fluid dynamics. In order to investigate the effect of the stratification on the near wake, in particular, the unsteady vortex formation behind a sphere, numerical simulations of stratified flows past a sphere are conducted. The time-dependent Navier–Stokes equations are solved using a three-dimensional finite element method and a modified explicit time integration scheme. Laminar flow regime is considered, and linear stratification of density is assumed under Boussinesq approximation. The effects of stratification is implemented by density transport without diffusion. The computed results include the characteristics of the near wake as well as the effects of stratification on the separation angle. Under increased stratification, the separation on the sphere is suppressed and the wake structure behind the sphere becomes planar, resembling that behind a vertical cylinder. With further increase in stratification, the wake becomes unsteady, and consists of planar vortex shedding similar to von Karman vortex streets.     

Effect of film cooling on the aerodynamic performance of an airfoil

The effect of film cooling on the aerodynamic performance of turbine blades is becoming increasingly important as the gas turbine operating temperature is being increased in order to increase the performance. The current paper investigates the effect of blowing ratio on the aerodynamic losses of a symmetric airfoil by pressure measurements and Particle Image Velocimetry (PIV). The test model features 4 rows of holes located on the suction side at 5%, 10%, 15% and 50% of the chord length. The Reynolds number based on the airfoil chord is 1.2 × 10⁵. Experiments are performed by varying the location of air injection, the angle of attack, and the mainstream velocity. The coolant air is injected at ambient temperature and the blowing ratio is varied from 0 to 1.91. It is observed that the losses due to film cooling increase with blowing ratio of 0 to 0.48, and the wake is shifted towards the suction side. Conversely, the aerodynamic losses decrease when the blowing ratio is increased further from 0.64 to 1.91. This trend has been observed for all the experimental configurations. The effect of blowing ratio on flow separation is investigated with the time-averaged velocity fields obtained from PIV measurements. It is observed that low blowing ratios, the separation point shifts upstream and at high blowing ratios the ejected coolant energizes the flow and delays separation. The pressure field around the airfoil is reconstructed from the integration of the Poisson equation based on the PIV velocity fields. The experimental results can be used for validation of numerical models for predicting losses due to film cooling.     

Parallel DNS with Local Grid Refinement

A parallel finite volume method for unstructured grids is used for a direct numerical simulation of the flow around a sphere at Re = 5000 (based on the sphere diameter and undisturbed velocity). The observed flow structures are confirmed by visualization experiments. A quantitative analysis of the Reynolds averaged flow provides a data base for future model evaluations.     

基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

Dielectric-barrier-discharge vortex generators: characterisation and optimisation for flow separation control

We investigated the use of dielectric-barrier-discharge plasma actuators as vortex generators for flow separation control applications. Plasma actuators were placed at a yaw angle to the oncoming flow, so that they produced a spanwise wall jet. Through interaction with the oncoming boundary layer, this created a streamwise longitudinal vortex. In this experimental investigation, the effect of yaw angle, actuator length and plasma-induced velocity ratio was studied. Particular attention was given to the vortex formation mechanism and its development downstream. The DBD plasma actuators were then applied in the form of co-rotating and counter-rotating vortex arrays to control flow separation over a trailing-edge ramp. It was found that the vortex generators were successful in reducing the separation region, even at plasma-to-free-stream velocity ratios of less than 10%.     

Numerical investigation of flow structures around a sphere

A numerical investigation of flow around a sphere is performed and compared with previous studies. Here, a second-order accurate, finite volume method is used in order to predict the instantaneous and time-averaged flow characteristics using large eddy simulation (LES) on the multi-block grid system. Namely, the objectives of this article are: (i) the presentation of flow structures in the wake region downstream of the sphere with a wide variety of flow properties such as the distribution of velocity vectors, patterns of streamlines, Reynolds stress correlations, root mean square of velocity components and other time-averaged flow data in order to reveal the vortical flow structures in detail and (ii) to demonstrate the abilities of computational methods in simulation of vortical flow data. Finally, it has been concluded that there are good agreements between the experimental results and numerical predictions.     

微通道壁面突起对片上电泳分离流动特性影响数值计算研究

针对制造过程引起的通道壁面突起对微流控芯片电泳分离使用的影响进行数值计算和理论分析研究.分析了壁面突起产生机制及其对电泳分离效果的影响;阐述了毛细管电泳分离的物理模型并进行了离散化;编制了电泳分离数值计算程序,采用有限体积法计算微通道内电参数分布、缓冲溶液流场分布和样品区带分布;给出微通道壁面突起高度和突起宽度对电泳分离过程影响的计算结果.从计算结果可知:壁面突起高度是影响电渗流流速的主要因素,当壁面突起高度与微通道宽度的比值从0增加至0.2,电渗流流速变化幅度约为20%.     

Finite element analysis of the viscous flow in a vaned radial diffuser

A new finite element method was used to analyze an experimental model of a radial vaned diffuser. The new method includes a streamline upwind formulation for the advection terms in the governing equations. The streamline upwind significantly reduces numerical diffusion while maintaining the stability of the conventional upwind formulation. The new finite element method also incorporates an iterative equal-order, velocity-pressure solution method based on the well-known SIMPLER algorithm. The results of the analysis are compared to flow visualization studies of the experimental model. The flow separation point for the four blade diffuser was predicted to occur at 19, 6% of the blade length from the leading edge. The experimentally determined value was 23% of the blade length. For the eight blade diffuser model, separation was predicted to occur at 43% of the blade length from the leading edge, as compared to the experimentally observed value of 50% of the blade length. With this performance comparison, the proposed finite element method has been demonstrated to be reliable for predicting complex fluid flows.     

Numerical investigation of cone angle effect on the flow field and separation efficiency of deoiling hydrocyclones

In this study, the effect of cone angle on the flow field and separation efficiency of deoiling hydrocyclones is investigated taking advantage of large eddy simulation. The dynamic Smagorinsky is employed to determine the residual stress tensor of the continuous phase. The method of Lagrangian particle tracking with an optimized search algorithm (closest cell) is applied to evaluate the separation efficiency of deoiling hydrocyclone. Simulations are performed on a 35-mm deoiling hydrocyclone with the three different cone angles of 6, 10 and 20 degree. The numerical results revealed that the changes in the cone angle would affect the velocity and pressure distribution inside hydrocyclone, and lead to changes in the separation efficiency. However, the large cone angle increases the tangential velocity and pressure gradient inside the hydrocyclone, but reduces the separation efficiency. The reasons behind the decrease in the separation efficiency are the flow structure and reduction of oil droplets residence time in hydrocyclones with large cone angles.     

Numerical solution for the flow of viscoelastic fluids around an inclined circular cylinder.

Finite difference solutions have been obtained by the perturbation method to investigate the influence of shear thinning and elasticity on the flow around an inclined circular cylinder of finite length in a uniform flow. In this numerical analysis a generalized upper-convected Maxwell model, in which the viscosity changes according to the Cross model, has been used.The local flow over the cylinder is only slightly deflected. However, in the wake flow behind the cylinder the particle path is remarkably influenced by the axial flow and rapidly flows up parallel to the cylinder's axis. Then it gradually rejoins direction of the incoming flow. It is found that viscoelastic fluids are prone to flow axially in the vicinity of the cylinder. The numerical predictions generally agree with the flow visualization results.The numerical solutions also demonstrate that elasticity has a strong effect on the velocity profile especially around both ends of the cylinder; elasticity increases the asymmetric profiles of both circumferential velocity and axial velocity with respect to equal to 90° and decreases a difference in the circumferential velocity between the windward end and the leeward end.For non-Newtonian fluids, the length of the wake flow is influenced by not only the Reynolds number but also the cylinder diameter and it is larger for the cylinder with the smaller diameter at the same Reynolds number.Partly presented at the 9th Australasian Fluid Mechanics Conference, University of Auckland, New Zealand, 8–12 December, 1986     

Influence of velocity profile on calibration function of Lorentz force flowmeter

A Lorentz force flowmeter is a noncontact electromagnetic flow-measuring device based on exposing a flowing electrically conducting liquid to a magnetic field and measuring the force acting on the magnet system. The measured Lorentz force is proportional to the flow rate via a calibration coefficient which depends on the velocity distribution and magnetic field in liquid. In this paper, the influence of different velocity profiles on the calibration coefficient is investigated by using numerical simulations. The Lorentz forces are computed for laminar flows in closed and open rectangular channels, and the results are compared with the simplified case of a solid conductor moving at a constant velocity. The numerical computations demonstrate that calibration coefficients for solid bodies are always higher than for liquid metals. Moreover, it can be found that for some parameters the solid-body calibration coefficient is almost twice as high as for a liquid metal. These differences are explained by analyzing the patterns of the induced eddy currents and the spatial distributions of the Lorentz force density. The result provides a first step for evaluating the influence of the laminar velocity profiles on the calibration function of a Lorentz force flowmeter.     

Drag and Separation Flow Past Solid Sphere with Porous Shell at Moderate Reynolds Number

Axisymmetric viscous, two-dimensional steady and incompressible fluid flow past a solid sphere with porous shell at moderate Reynolds numbers is investigated numerically. There are two dimensionless parameters that govern the flow in this study: the Reynolds number based on the free stream fluid velocity and the diameter of the solid core, and the ratio of the porous shell thickness to the square root of its permeability. The flow in the free fluid region outside the shell is governed by the Navier–Stokes equation. The flow within the porous annulus region of the shell is governed by a Darcy model. Using a commercially available computational fluid dynamics (CFD) package, drag coefficient and separation angle have been computed for flow past a solid sphere with a porous shell for Reynolds numbers of 50, 100, and 200, and for porous parameter in the range of 0.025–2.5. In all simulation cases, the ratio of b/a was fixed at 1.5; i.e., the ratio of outer shell radius to the inner core radius. A parametric equation relating the drag coefficient and separation point with the Reynolds number and porosity parameter were obtained by multiple linear regression. In the limit of very high permeability, the computed drag coefficient as well as the separation angle approaches that for a solid sphere of radius a, as expected. In the limit of very low permeability, the computed total drag coefficient approaches that for a solid sphere of radius b, as expected. The simulation results are presented in terms of viscous drag coefficient, separation angles and total drag coefficient. It was found that the total drag coefficient around the solid sphere as well as the separation angle are strongly governed by the porous shell permeability as well as the Reynolds number. The separation point shifts toward the rear stagnation point as the shell permeability is increased. Separation angle and drag coefficient for the special case of a solid sphere of radius r = a was found to be in good agreement with previous experimental results and with the standard drag curve.     

Generalized integral transform solution for boundary layer flow over a sphere

The generalized integral transform technique is applied to the boundary layer equations for flow over a sphere in their primitive variables. Even though a diffusion-based eigenvalue problem is used, the velocity profile, shear stress and separation point have been calculated with high accuracy. Low-order approximations are shown to be accurate near the surface and the predictions of the separation point is very good. Comparison with finite difference results shows the better convergence behaviour of the integral transform method.     

Vortex dynamics of a sphere wake in proximity to a wall

Toward getting the vortex dynamics characteristics and wake structure of a sphere in proximity to a wall, the effect of a proximal flat plate on the wake of a stationary sphere is investigated via direct numerical simulation. The vortex shedding process and the significant variation of the wake structure are described in detail. The drag coefficient reduces and the wake structure of the sphere becomes complex due to the combined effect of the wake flow and the wall. A jet flow forms between the sphere and the flat plate, which suppresses the vortex separation on the bottom of the sphere. The asymmetric distributions of the coherent structures and the recirculation region behind the sphere are discussed. Besides vortex shedding patterns, the time-averaged velocity distribution, vortex dynamics, distribution regularities of turbulent kinetic energy and enstrophy are investigated.     

Mass and heat transfer during solid sphere dissolution in a non-uniform fluid flow

We studied a nonisothermal dissolution of a solvable solid spherical particle in an axisymmetric non-uniform fluid flow when the concentration level of the solute in the solvent is finite (finite dilution of solute approximation). It is shown that simultaneous heat and mass transfer during solid sphere dissolution in a uniform fluid flow, axisymmetric shear flow, shear-translational flow and flow with a parabolic velocity profile can be described by a system of generalized equations of convective diffusion and energy. Solutions of diffusion and energy equations are obtained in an exact analytical form. Using a general solution the asymptotic solutions for heat and mass transfer problem during spherical solid particle dissolution in a uniform fluid flow, axisymmetric shear flow, shear-translational flow and flow with parabolic velocity profile are derived. Theoretical results are in compliance with the available experimental data on falling urea particles dissolution in water and for solid sphere dissolution in a shear flow.