航空发动机轴承腔润滑的气液两相均匀流研究

基于轴承腔中润滑油的两相均匀流动模型,采用湍流模式和有限差分数值方法计算轴承腔内三维定常N-S方程,对腔内润滑油的气液两相均匀流动特性进行研究,以获得气液两相均匀流条件下润滑油流场、压力场和速度场在轴承腔内的分布情况,分析转子转速、含气率和润滑油进口速度对润滑油出口压力以及润滑油与壁面之间剪切力的影响,同时对单相流和两相均匀流润滑性能差异进行比较.结果表明,转子转速、含气率和润滑油进口速度对润滑油出口压力和腔内壁面与润滑油间的平均剪切力具有不同影响,而采用2种流动模型计算出的轴承腔润滑油出口压力的差异较大,同时支持了开展航空发动机轴承腔润滑两相流动分析的必要性.     

模型燃烧室内瞬态紊流的大涡模拟

本文采用三种不同亚网格尺度模型对带有V型稳定器的模型燃烧室二维瞬态紊流流动进行了大涡模拟。并在交错网格系下用SIMPLE算法和混合差分格式求解离散方程。数值研究拟不同型式入口速度分布和不同亚网格尺度模型下模型燃烧室二维瞬态紊流流场。计算结果表明不同入口速度分布和不同亚网格尺度模型对瞬态流场和出口速度分布有一定的影响。本文通过数值模拟,揭示了V型稳定器后旋涡的产生和脱落过程。通过计算结果及实验数据的比较可知,本文采用的亚网格尺度模型可以用来模拟模型燃烧室紊流流场及稳定器后面回流区的流动情况。     

搅拌槽内气液流动大涡数值模拟研究

为优化搅拌槽体型设计、提高搅拌效率,分别利用气液两相大涡模型和标准k-ε紊流模型,结合多参考系法,对搅拌槽内的气液混合过程进行了模拟计算;控制方程的离散使用了有限体积法;速度与压力耦合求解时使用了压力隐式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法;模拟自由水面采用了VOF(Volume of Fluid)法;通过计算得到了搅拌槽在相同通气率、不同转轮转速下的气液流动速度场以及不同截面含气率的分布规律。比较两种模型,模拟结果表明:叶片旋转区域紊流的各向异性随着转速的增加而明显加强;基于各向同性假设的标准k-ε紊流模型不能描述流体流速的波动变化,而大涡模拟能够捕捉叶片附近区域流场的分布规律;在通气率一定的情况下,转轮转速的大小对下叶轮附近区域的含气率影响较小,而对上叶轮附近区域含气率的影响较大。     

定床弯道内水沙两相运动的数值模拟

在适体同位网格中采用非正交曲线坐标系下的三维k-ε-kp固液两相双流体湍流模型研究弯道内水流和悬浮泥沙运动,主要计算了试验室S型水槽内清水流动的三维流场、120°弯道内水沙两相流动中底沙与底流的运动轨迹以及S型水槽内水沙两相流动的两相流场和泥沙浓度场. 对于S型水槽内清水流动,数值结果与试验结果吻合良好. 120°弯道内水沙两相流动中固液两相的运动轨迹在弯道直线段基本重合,在弯道内泥沙轨迹逐步偏离水体轨迹,其偏离程度随泥沙粒径增大而增大. 从S型水槽内水沙两相流动计算结果中发现泥沙纵向流速在壁面附近比水流纵向速度大,在远离壁面区域比水流纵向速度小;弯道内泥沙横向流速比水流横向流速小;垂向流速在直线段和泥沙沉速相当,在弯道内受螺旋水流影响而变化;两相流速差别随泥沙粒径增大而变大;泥沙浓度呈现下浓上稀的分布,在弯道内横向断面上呈现凸岸大凹岸小的分布,泥沙浓度随泥沙粒径增大而减小.      

绕水翼超空化流动形态与速度分布

为揭示超空化流场结构特性,利用高速全流场显示技术,观察了绕hydronautics水翼的超空化流动形态,并利用数字粒子图像测速仪(DPIV)测量了其速度分布. 在测量空穴内部流速分布时,采用空化流场中的空化泡作为示踪粒子来显示流动结构. 结果表明:随着空化数的降低,超空化流动显现出了明显的阶段特征,其中水汽混合相和汽相的分布决定了空化区域的形态与流速分布;空化区和主流区的汽液交界面处存在着较大的速度梯度;低速分布区域随着空化数的降低由水翼吸力面中后部向水翼下游移动;在空化区域内部,水汽混合区的速度相对较低,而汽相区则与主流区有着相近的速度分布.关键词超空化水翼、DPIV、高速摄像、空化形态、流速分布      

不同重力下90°弯管内气液两相流流型及流动特性研究

研究了不同重力条件下90°弯管内气液两相流流型分布形态及流动特性. 通过建立90°弯管内气液两相流流动的三维数学物理模型,采用VOF 方法,对10-6g₀, 10-4g₀, 10-2g₀, 1g₀ (g₀= 9.8m/s2) 重力下的90°弯管内气液两相流流型分布特征、截面空隙率、滑速比及气相尾部最大斜向角进行了比较分析. 研究结果表明:所建立的模型能够正确模拟不同重力条件下90°弯管内气液两相流流型和截面空隙率,并得到气液两相弯管二次流与单相二次流的不同特性. 随着重力水平的提高,90°弯管对气相流型的影响作用减弱,气相整体向弯管内侧积聚靠拢,弯管对尾部的斜向作用减弱.      

激光多普勒测速技术在气液两相流中的应用

本文采用激光颗粒动态分析仪(PDA)测量了钢包底部喷吹气液两相流中气泡的直径和气泡上升速度的分布;采用激光多普勒测速仪(LDA)测量了气液两相区和液体单相循环区液体速度场的分布。测量结果表明:气泡在脱离喷嘴上浮一定距离后,其大小基本保持不变;在气液两相区中,气泡速度和液体速度的分布均服从高斯分布;液体在单相区作循环流动,在侧壁与底部交接处,存在液体流动的“死区”。     

大尺寸液桥热毛细对流失稳性地面实验研究

开展大尺寸液桥浮力-热毛细对流地面实验, 探究流场转捩的临界条件及临界状态附近的流动情况. 通过粒子图像测速方法(PIV) 获得流体速度场, 研究液桥内部定常和转捩后的流场结构以及流体运动规律;并用红外热像仪测量液桥自由面温度分布, 研究流体流动的时空演化和温度振荡. 实验发现大尺寸半浮区液桥浮力-热毛细对流临界值与几何参数有关, 在大普朗特(Prandtl) 数情况下, 流场存在由稳定态向不稳定态再到混沌的转捩过程, 在临界马兰哥尼(Marangoni) 数附近, 流场内会出现行波现象, 流动模式也会随高径比的变化而发生变化;当继续增大马兰哥尼数, 流动会进入混沌状态.      

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性;热通量辐射模型估算辐射换热,分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解,液相采用Lagrange处理,并用PSIC算法对其进行求解．通过实验结果和计算结果的比较,表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟,3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程,而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

微扩散管内幂律流体在不同外电场驱动下非稳态流动数值仿真

电场驱动下的非牛顿流体在微米级扩散管道内非稳态电渗流动特性是MEMS管设计人员关注的焦点,大部分实际液体可近似为用幂律模型描述的纯粘性流体,所以论文针对幂律流体在有限长微扩散管道内在两种不同形式的外加电场驱动下的非稳态电渗流动情况进行数值仿真.基于Ostwald-De Wael幂律模型和连续介质假说,采用高精度紧致有限差分离散二维完全Poisson-Boltzmann电势方程和Cauchy动量方程,对恒定电场及满足Maxwell方程的电场进行数值仿真,讨论了微扩散管中幂律流体在两种不同外加电场驱动下的瞬时流场分布的差异.结果表明,初始时刻固定扩散角和无量纲壁面电势,无量纲电动宽度的变化对幂律流体电渗流速度分布影响较大;在微扩散管上游等截面处,由恒定电场驱动及Maxwell电场驱动电渗流速度分布差别极小,在扩散管中下游则出现了明显的差别;由恒定电场驱动下的电渗流动在扩散管不同截面下的速度峰值相近,但Maxwell电场诱导的电渗流速度峰值则随管道半径变化出现较大差别.对于外加电场驱动的电渗流动,不同形式的外电场可使流场产生较大差别,而不同性质的流体也会形成不同的流场分布.     

Investigation of slug flow characteristics in the valley of a hilly-terrain pipeline

The objective of this study is to improve the current phenomenological understanding of slug flow characteristics over an entire hilly-terrain section, and in particular, the slug initiation mechanism at the lower dip.The experimental part of this study revealed that five possible flow behavior categories exist along a hilly-terrain section. In these categories, the flow behavior at the dip is coupled with flow conditions of the upstream downhill section. This qualitative classification was superimposed on steady-state flow pattern maps for the upstream downhill section in an attempt to relate the qualitative flow behavior at a dip to the flow pattern maps through the flow behavior in the downhill section.Statistical analyses of mean slug length, maximum slug length, slug frequency, and slug length variation across the hilly-terrain pipeline revealed that slug length distribution characteristics change across a symmetrical hilly-terrain pipeline. Physical modeling of the slug initiation mechanism and the characteristics of initiated slugs at the lower dip indicated two main mechanisms, namely, wave growth and wave coalescence initiation mechanisms. The initiated “pseudo slugs” or slug characteristics of each mechanism differ significantly with respect to frequency, length, liquid holdup and velocity. It was observed that pseudo slugs initiated by the wave coalescence mechanism have velocities less than the mixture velocity due to gas blowing through the slug body.     

Two‐dimensional modeling for stability analysis of two‐phase stratified flow

The effect of wavelength and relative velocity on the disturbed interface of two‐phase stratified regime is modeled and discussed. To analyze the stability, a small perturbation is imposed on the interface. Growth or decline of the disturbed wave, relative velocity, and surface tension with respect to time will be discussed numerically. Newly developed scheme applied to a two‐dimensional flow field and the governing Navier–Stokes equations in laminar regime are solved. Finite volume method together with non‐staggered curvilinear grid is a very effective approach to capture interface shape with time. Because of the interface shape, for any time advancement, a new grid is performed separately on each stratified field, liquid, and gas regime. The results are compared with the analytical characteristics method and one‐dimensional modeling. This comparison shows that solving the momentum equation including viscosity term leads to physically more realistic results. In addition, the newly developed method is capable of predicting two‐phase stratified flow behavior more precisely than one‐dimensional modeling. It was perceived that the surface tension has an inevitable role in dissipation of interface instability and convergence of the two‐phase flow model. Copyright © 2009 John Wiley & Sons, Ltd.     

The development of the structure of water – air bubble grid turbulence

The development of the turbulent flow field generated due to the interaction of grid turbulence with a swarm of bubbles is investigated experimentally, in a vertical channel of rectangular cross section. Void fraction and streamwise mean and rms velocity distributions have been measured at several distances from the grid, with an optical probe and Laser Doppler Velocimetry, in relation to the air flow rate ratio. The obtained results indicate that close to the grid the void fraction and velocity distributions are dictated by the bubble injectors’ location on the grid. Downstream the void fraction distribution changes to a double peak pattern. The velocity distribution is characterized by a shear layer between the wall area and the central area of the channel. The extend of this shear layer is increasing as the distance from the grid and the gas flow rate ratio are increasing, and is associated with a corresponding increase of the turbulence fluctuations. Autocorrelation and spectra measurements at the centre of the channel show a reduction of the flow scales for low void fraction. Consistently, power spectra distributions indicate that bubbles cause a redistribution of energy manifested by the relative enhancement of the intermediate scales’ energy content and a consequent reduction in the larger scales. These trends are gradually alleviated and reversed at large distances from the grid, as the air flow rate is increased.     

The use of perforated plates to control the flow emerging from a wide-angle diffuser,with application to electrostatic precipitator design

The flow through a wide-angle, pyramidal diffuser of area ratio 6.8, in which two perforated plates are used to control the emergent velocity distribution was investigated. (Wide-angle diffusers combined with perforated plates find application in electrostatic precipitator flow systems. The efficacy of these systems depends critically on the uniformity of the gas flow within the collection chamber downstream of the diffuser outlet plane.) The main results of the study are (i) the establishment of the main factor influencing the flow at the diffuser exit plane, (ii) the determination of plate characteristics which produce a uniform velocity profile in the collection chamber, and (iii) the establishment of the corresponding pressure drop characteristics of the plate—diffuser combinations. The results also extend the range of diffuser geometries for which two perforated plates provide uniform velocity profiles at exit.     

Void fraction and pressure waves in a transient horizontal slug flow

The void fraction and the pressure waves in an air–water mixture flowing in the slug regime are experimentally investigated in a horizontal line. The test section is made of a transparent Plexiglas pipe with 26 mm ID and 26.24 m long, operating at ambient temperature and pressure. The flow induced transients are made by quickly changing the air or the water inlet velocity. The test grid has four operational points. This choice allows one to create expansion and compression waves due to the changes to the gas or to the liquid. Each experimental run is repeated 100 times to extract an ensemble average capable of filtering out the intrinsic flow intermittence and disclosing the void fraction and pressure waves’ features. The slug flow properties such as the bubble nose translational velocity, the lengths of liquid film underneath the bubble and the liquid slug are also measured. The objective of the work is two-fold: access the main characteristics of the void fraction and pressure waves and disclose the mechanics of the transient slug flow as described through the changes of the slug flow properties.     

Specified discharge velocity models for numerical simulations of laminar vortex rings

We numerically and theoretically investigate the flow generated at the exit section of a piston/cylinder arrangement that is generally used in experiments to produce vortex rings. Accurate models for the velocity profile in this section (also called specified discharge velocity, SDV models) are necessary in (i) numerical simulations of laminar vortex rings that do not compute the flow inside the cylinder and (ii) in slug-models that provide a formula for the total circulation of the flow. Based on the theoretical and numerical analysis of the flow evolution in the entrance region of a pipe, we derive two new and easy to implement SDV models. A first model takes into account the unsteady evolution of the centerline velocity, while the second model also includes the time variation of the characteristics of the boundary layer at the exit plane of the vortex generator. The models are tested in axisymmetric direct numerical simulations of vortex rings. As distinguished from classical SDV model, the new models allow to accurately reproduce the characteristics of the flow. In particular, the time evolution of the total circulation is in good agreement with experimental results and previous numerical simulations including the vortex generator. The second model also provides a more realistic time evolution of the vortex ring circulation. Using the classical slug-model and the new correction for the centerline velocity, we finally derive a new and accurate analytical expression for the total circulation of the flow.     

Numerical investigation of turbulent-jet primary breakup using one-dimensional turbulence

Primary breakup to form droplets at liquid surfaces is an important fundamental process to study as it determines the initial properties of the dispersed phase, which affect mixing rates, secondary breakup, droplet collisions, and flow separation within the dispersed flow region. Primary breakup can be regarded as one of the least developed model components for simulating and predicting liquid jet breakup. However, it is of paramount importance in many technical applications, e.g. fuel injection in engines and spray painting. This paper presents a numerical investigation of primary breakup of a turbulent liquid jet in still air at standard conditions using the one-dimensional turbulence (ODT) modeling framework. ODT is a stochastic model that simulates turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. An important feature of ODT is the resolution of all relevant scales, both temporal and spatial. The restriction to one spatial dimension in ODT permits affordable high resolution of interfacial and single-phase property gradients, which is key to capturing the local behavior of the breakup process and allows simulations at high Reynolds and Weber numbers that are currently not accessible to direct numerical simulations (DNS).This paper summarizes our extensions of the ODT model to simulate geometrically simple jet breakup problems, including representations of Rayleigh wave breakup, turbulent breakup, and shear-driven breakup. Each jet breakup simulation consists of a short temporal channel section to initialize a turbulent velocity profile at the nozzle exit followed by an adjacent jet section. The simulations are carried out for jet exit Reynolds number of 11,500, 23,000, 46,000 and 92,000 while the Weber number is varied within the range 10²–10⁷. We present results on breakup statistics including spatial locations of droplet release, droplet sizes and liquid core length. The results on primary breakup are compared to experimental results and models.     

Entrainment in condensing annular flow

Experiments were carried out on low pressure, steam-water, condensing annular flow in a 38.1 mm i.d. horizontal tube. The velocity of the steam at inlet was in the range 97–186 m/s.Measurements of the liquid film flow rate at the end of the test section, which arose as a result of condensation, entrainment, and deposition, were made by extracting the film through a porous sinter bush. The liquid flow rate in the vapour core at exit was deduced from these measurements together with a heat balance on the condenser section.These results were compared with three correlations for entrainment developed from air-water studies. On the basis of the experimental data available, there was sufficient agreement in one case to warrant further investigation.     

The structure of water-air bubble grid turbulence in a square duct

Local measurements of void fraction and continuous phase velocity field in water-air bubble, grid turbulence were conducted in a channel of vertical, square test section. The measured statistics indicate that, due mainly to the interaction of mean shear with the dispersed phase, the turbulence structure of the flow is modified. The observed change is characterized by a strong spatial dependence of void fraction and liquid flow properties, and the emergence of two spatial regions controlled by different physical processes. Intensity measurements indicate significant departure from isotropy in the flow. Two distinct regimes corresponding to low and high values of void fraction have been also identified. The autocorrelation and spectra measurements indicate that for low void fraction the scales of turbulence decrease while for higher values of void fraction increase again and inverse cascade is observed.