Compressible gas-liquid mixture flow at abrupt pipe enlargements

A detailed investigation was made of the flow of compressible gas-liquid mixtures through sudden enlargements in diameter of circular pipes. One-dimensional analysis shows that the dimensionless pressure rise varies with mixture void fraction and mixture momentum, while the establishment of choking conditions at the enlargement is controlled by the length of pipe downstream in which frictional pipe flow occurs. The flows were found to exhibit two characteristic modes, jet flow and submerged flow, with intermediate flows displaying unsteady oscillation between these modes. The distance to the downstream position of maximum pressure increased steadily with mixture void fraction when the upstream pipe outlet was choked, varying from 5 to 50 times the downstream pipe diameter. If the flow was not choked, this distance was much smaller and showed discrete fixed values associated with the mode of flow.

One-dimensional analysis accurately predicted maximum pressure, but when flow was choked at the enlargement the calculation was sensitive to the pressure in the region of separated flow surrounding the central jet in the enlargement. Although analysis of maximum pressure in terms of flow expansion and normal shock gave a general indication of the maximum pressure (which was thus concluded to depend on the general flow processes expected in the enlargement), accurate prediction of maximum pressures will depend on empirical knowledge of the separated flow region pressures. The maximum pressure rise was found to be in the range extending down to 0.3 of the upstream pipe outlet pressure and reduced with void fraction; it was also influenced by the enlargement area ratio. Flows in the approach and outlet pipes were found to be compressible, frictional pipe flows of the Fanno type, with somewhat reduced friction factors occurring in the outlet pipe.     

3D numerical simulation of compressible swirling flow induced by means of tangential inlets

The objective of the present work is to predict compressible swirl flow in the nozzle of air‐jet spinning using the realizable k–ε turbulence model and discuss the effect of the nozzle pressure. The periodic change of flow patterns can be observed. The recirculation zone near the wall of the injectors upstream increases in size and moves gradually upstream, whereas the vortex breakdown in the injector downstream shifts slowly towards the nozzle outlet during the whole period. A low axial velocity in the core region moves gradually away from the centerline, and the magnitude of the center reverse flow and the area occupied by it increase with axial distance due to the vortex breakdown. From the tangential velocity profile, there is a very small free‐vortex zone. With increasing nozzle pressure, the velocity increases and the location of vortex breakdown is moved slightly downward. However, the increase in the velocity tends to decline at nozzle pressure up to a high level. Copyright © 2008 John Wiley & Sons, Ltd.     

三角翼前缘涡的某些破裂形式及特性研究

基于流动显示和PIV技术测量的实验结果，对三角翼前缘涡破裂的一些形式和破裂特性进行了分析和讨论。通过PIV测量所得到的涡量分布证实了在螺旋破裂的情况下，涡核的螺旋方向与前缘涡的旋转方向相反，及双螺旋破裂形式的存在等。进而对螺旋波的形成机理提出了与有关文献不同的看法。     

Some Observations of Vortex Breakdown in a Confined Flow with Solid Body Rotation

The occurrence of breakdown in slender vortex flows as a ``bubble'' or ``spiral'' pattern depends on the degree of radial deflection of the vortex core from its original axis as shown in [1]. A smooth transition from a bubble to a spiral-type ``mode'' can be forced by inducing a small asymmetric disturbance which led to the conclusion, that the patterns do not represent different fundamental modes of breakdown. The subject presented herein addresses the following question: how does breakdown evolve in a swirling flow in which the vortex core is forced on a straight axis? In addition, what is the effect of turbulent inflow conditions? This type of vortex conditions is achieved in a spinning tube flow. The swirl is introduced at the entrance of the rotating tube with a honeycomb package and maintained by the viscous action in the boundary layer of the spinning tube. A diffuser at the end induces an adverse pressure gradient to force the breakdown. Flow visualization experiments are carried out to characterize the nature of breakdown over a range of different flow conditions. For some selected characteristic stages, detailed velocity fields were obtained using the method of Digital Particle-Image-Velocimetry (DPIV). The results show, that for the range of parameters investigated, breakdown is initiated at Rossby-numbers below a critical value of Ro ≈ 0.6 similar to those observed in other experiments. The bursted part of the vortex has a near axi-symmetric slender conical shape containing approximately stagnant flow. Its downstream end is characterized by a jump-like contraction where the flow evolves into a jet with enhanced swirl on the axis. It is only in this region downstream of the jump-like contraction that asymmetric instabilities and wavy flow patterns could be observed. Perturbations caused by them travel upstream but do not change the near-axisymmetric shape of the bursted part of the vortex.     

可压缩流向涡与激波轴对称干扰的数值模拟

用ＮＳ方程数值模拟了可压缩流向涡和激波轴对称相互作用现象．数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较．结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征．提出了判断流向涡与运动激波相互作用中旋涡破碎的准则     

Investigation on Secondary Flow Characteristics in a Curved Annular Duct with Struts

Under the influence of duct curvature, cross-sectional area variation and internal struts, the internal flow field within a curved annular duct becomes rather complicated and contains strong secondary flow. In this paper, the secondary flow characteristics in an annular duct with struts are experimentally and numerically investigated. The results show that large pressure gradients exist on the bends of hub and shroud. Meanwhile, two counter-rotating vortex pairs appear both along the hub-side and shroud-side surfaces. The hub-side vortex pair of which the vortex cores travel downstream parallelly evolves from the horseshoe vortex which is induced by the leading edge of the upstream strut, whereas the shroud-side vortex pair originates from the strut trailing edge and the corresponding vortex cores develop in a divergent way. Additionally, the effects of the duct exit Mach number on the secondary flow characteristics are also studied. As the exit Mach number increases, the streamwise pressure gradients increase and lead to more intense vortices, higher total pressure loss and larger flow distortion.     

Effect of Volume Energy Supply on the Stability of a Subsonic Vortex Flow

Calculations of the stability of an axisymmetric vortex flow of viscous heat-conducting gas with volume energy supply are presented. The unperturbed axisymmetric vortex flow was found numerically using a quasi-cylindrical approximation of the Navier-Stokes equations under the assumption of constant peripheral-velocity circulation in the ambient co-current flow. The volume energy supply in the viscous vortex core was modeled by an additional source term in the energy equation. The stability characteristics of the viscous vortex flow in a longitudinal vortex with respect to both axisymmetric and non-axisymmetric three-dimensional waves traveling along the vortex axis and corresponding to both positive and negative values of the azimuthal wave number were found using the time-dependent formulation of the linear stability theory for compressible three-dimensional plane-parallel flows.     

Vortex structure simulation for supersonic mixing layers using nonlinear PSE method

The method of nonlinear parabolized stability equations (PSE) is applied in the simulation of vortex structures in compressible mixing layer. The spatially-evolving unstable waves, which dominate the vortex structure, are investigated through spatial marching method. The instantaneous flow field is obtained by adding the harmonic waves to basic flow. The results show that T-S waves do not keep growing exponentially as the linear evolution, the energy transfer to high order harmonic modes, and that finally all harmonic modes get saturated due to nonlinear interaction. The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones, makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up, vortex pairing, and Λ vortex.     

A method for determining valve coefficient and resistance coefficient for predicting gas flowrate

A method is introduced to determine the valve flow coefficient and resistance coefficient with the experiment of air discharging from a reservoir, and with the least squares method to fit the cumulative molar quantities discharged. The test valve is an angle-seat valve (Type 2632, Bürkert) with different apertures. At pressure difference of about 6 bar, the choked flow occurs when the valve aperture over 60%. Both the valve coefficient and resistance coefficient model can exactly predict the flowrate for the non-choked flow, while there are larger deviations for the choked flow. The modified equation for the choked flow can improve the prediction. In the resistance coefficient model, the value of resistance coefficient and the discharged cumulative molar quantities obtained with both the compressible and incompressible assumption are very close. The compressibility of air is negligible within the experimental pressure difference of about 6 bar. The additivity of the resistance coefficient makes the model more convenient to use.     

On the possible formation of a vortex flow in a supersonic inlet with three-dimensional compression

The structure of the flow in an inlet with three-dimensional compression and an adjoining channel of rectangular cross-section is experimentally investigated for the external flow parameters M_∞ = 4 and Re ≈ 52 · 10⁶ 1/m. The phenomenon of vortex generation in the air-intake and its breakdown at the channel entry, where it encounters an elevated pressure gradient, is studied     

基于速度梯度张量特征值的陷窝内旋涡分析

作为一种新型的涡流发生器,陷窝具有流动阻力小、综合传热性能高的特点,是现代高性能涡轮叶片内部冷却新技术. 旋涡的定量分析是陷窝强化传热优化设计的重要依据. 针对在不同陷窝模型下的旋涡结构、分离方式和背景压力变化引起的旋涡强度无法定量分析的问题,本文提出采用涡核速度和 涡核速度梯度张量特征值来定量分析旋涡的方法. 通过采用涡核处局部坐标系表示的速度矢量和速度梯度张量,得到了涡核的轴 向速度、径向速度、旋转角速度、轴向加速度和径向加速度,并在此基础上简化出了用最大轴向速度、最大轴向加速度和最大旋 转角速度综合表示的旋涡强度的定量分析方法. 用该方法分析了不同深宽比陷窝诱导的旋涡结构,随着深宽比的增大,最大轴向 速度、最大轴向加速度和最大旋转角速度均呈现明显的增大趋势,旋涡强度增大. 研究表明此方法具有数据处理简单、通用性强、 不受分离方式限制、不受背景压力影响的特点,且提取到的数据具有明确的物理意义,适用于各类旋涡定量分析.      

Rollup region of a turbulent trailing vortex issued from a blade with flow separation

This paper presents the results of an experimental investigation on the near field of a tip vortex generated by a blade at moderate incidence. The experiments were conducted at Re=15 000 and the boundary layer over the blade separated around midchord on the upper surface. Laser-Doppler measurements of the turbulent flow (Tu=1.5%) were performed at various stations downstream of the blade. The three components of the mean velocity field and turbulent attributes were quantified at cross-planes, characterizing both the blade wake and the tip vortex structure. This allowed the analysis of the rollup and initial stages of decay of the tip vortex in the light of known theories and models. The axial velocity defect at the center of the vortex core evolved as x⁻¹ log x, without displaying any significant outgrowth imposed by the separated flow upstream. Momentum balances were also carried out at a station downstream to the conclusion of vortex rollup. The approximate axisymmetry of the flow field in the trailing vortex was used to formulate the balances in a cylindrical coordinate system. Among other observations, it was seen that an adverse axial pressure gradient developed in the vortex core, which reinforced the tenacity of the axial velocity defect. In contrast, an area influenced by a favorable pressure gradient was found outside the core.     

The expansion of wet steam through a compressible confined vortex in a fluidic vortex diode

The pressure-flow characteristics of a Zobel-type vortex diode have been measured using a working fluid of compressible wet steam. Tests using superheated steam with inlet/outlet pressure ratios across the diode of up to 30 have shown clearly the effects of compressibility and choking on the diode characteristics. Repeating the tests using wet steam, with known dryness fractions, has shown separately the effects of wetness on the diode performance.

When the diode was installed into the pipwork in the high-resistance direction, excessive steam wetness (quality <0.93) led to a build-up of water and when this was eventually swept through to the diode the resistance was seen to fall substantially as the strong internal vortex was destroyed.     

Effusing core at the center of A potential vortex

Experiments were carried out to measure the base pressure distribution of a flow field induced by a potential vortex with its axis normal to a stationary disk. The center base pressure coefficient of the vortex, C₀(0), was found to be proportional to Reynolds number from Re = 2.0 × 10³ to Re > 2.5 × 10⁴, where Re is based on the disk radius and azimuthal velocity at the disk edge. This behavior of C₀(0) is at variance with the experimental results of Phillips (Phys. Fluids, 27, 2215, 1984) and Khoo (M. Eng. Thesis, Natl. Univ. Singapore, 1984), which showed vastly different trends depending on Re. Plausible reasons are suggested for the apparent discrepancies observed. Finally, the extent of the effusing core at the center, r₁ (taken to be the radial position where departure from the outer potential flow took place), was found to be proportional to Re^−1/2 for all Re values considered.     

可压缩流向涡与激波轴对称干扰的数值模拟

用ＮＳ方程数值模拟了可压缩流向涡和激波轴对称相互作用现象。数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较,结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征。提出了判断流向涡与运动激波相互作用中旋涡破碎的准则。     

Density and compressibility effects in turbulent subsonic jets part 1: mean velocity field

The behavior of compressible jets originated from initially turbulent pipe flows issuing in still air has been investigated at three different subsonic Mach numbers, 0.3, 0.6 and 0.9. Helium, nitrogen and krypton gases were used to generate the jet flows and investigate the additional effects of density on the flow structure. Particle image velocimetry, high-frequency response pressure transducers and thermocouples were used to obtain velocity, Mach number and total temperature measurements inside the flow field. The jets were formed at the exit of an adiabatic compressible frictional turbulent pipe flow, which was developing toward its corresponding sonic conditions inside the pipe, and continued to expand within the first four diameters distance after it exited the pipe. Theoretical considerations based on flow self-similarity were used to obtain the decay of Mach number along the centerline of the jets for the first time. It was found that this decay depends on two contributions, one from the velocity field which is inversely proportional to the distance from the exit and one from the thermal field which is proportional to this distance. As a result, a small non-linearity in the variation of the inverse Mach number with downstream distance was found. The decay of the Mach number at the centerline of the axisymmetric jets increases by increasing the initial Mach number at the exit of the flow for all jets. The decay of mean velocity at the centerline of the jets is also higher at higher exit Mach numbers. However, the velocity non-dimensionalized by the exit velocity seems to decrease faster at low exit Mach numbers, suggesting a reduced mixing with increasing exit flow Mach numbers. Helium jets were found to have the largest spreading rate among the three different gas jets used in the present investigation, while krypton jets had the lowest spreading rate. The spreading rate of each gas decreases with increasing its kinetic energy relatively to its internal energy.     

Flow regime transitions due to cavitation in the flow through an orifice

This paper presents both experimental and theoretical aspects of the flow regime transitions caused by cavitation when water is passing through an orifice. Cavitation inception marks the transition from single-phase to two-phase bubbly flow; choked cavitation marks the transition from two-phase bubbly flow to two-phase annular jet flow.

It has been found that the inception of cavitation does not necessarily require that the minimum static pressure at the vena contracta downstream of the orifice, be equal to the vapour pressure liquid. In fact, it is well above the vapour pressure at the point of inception. The cavitation number [σ = (P₃ − P_v)/(0.5 pV²); here P₃ is the downstream pressure, P_v is the vapour pressure of the liquid, ρ is the density of the liquid and V is the average liquid velocity at the orifice] at inception is independent of the liquid velocity but strongly dependent on the size of the geometry. Choked cavitation occurs when this minimum pressure approaches the vapour pressure. The cavitation number at the choked condition is a function of the ratio of the orifice diameter (d) to the pipe diameter (D) only. When super cavitation occurs, the dimensionless jet length [L/(D - d); where L is the dimensional length of the jet] can be correlated by using the cavitation number. The vaporization rate of the surface of the liquid jet in super cavitation has been evaluated based on the experiments.

Experiments have also been conducted in which air was deliberately introduced at the vena contracta to simulate the flow regime transition at choked cavitation. Correlations have been obtained to calculate the critical air flow rate required to cause the flow regime transition. By drawing an analogy with choked cavitation, where the air flow rate required to cause the transition is zero, the vapour and released gas flow rate can be predicted.     

3D tangentially injected swirling recirculating flow in a nozzle with a slotted‐tube—Effects of groove parameters

A numerical prediction for 3D swirling recirculating flow in an air‐jet spinning nozzle with a slotted‐tube is carried out with the realizable k–ε turbulence model. The effects of the groove parameters on the flow and yarn properties are investigated. The simulation results show that some factors, such as reverse flow upstream of the injector, vortex breakdown downstream of the injector, corner recirculation zone (CRZ) behind the step and vortex ring in the groove caused by the groove geometric variation, are significantly related to fluid flow, and consequently to yarn properties. With increasing groove height, the length of the CRZ increases, while the initial vortex ring in the groove decreases and a same direction rotating vortex forms in the bottom of the groove. Similarly, as the groove width increases, the extent of both vortex breakdown in downstream of the injectors and the vortex ring in the groove increases slightly, whereas the CRZ lengths in stream‐wise direction decrease. Some factors, such as the negative tangential velocities, the size of the vortex rings in the grooves and the CRZ, are constant for nozzles with different groove lengths. Copyright © 2009 John Wiley & Sons, Ltd.     

Computation of field structure and aerodynamic characteristics of delta wings at high angles of attack

A numerical investigation of the structure of the vortical flowfield over delta wings at high angles of attack in longitudinal and with small sideslip angle is presented. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics that are dominated by the effect of the breakdown of the leading-edge vortices. The methods that analyze the flowfield structure quantitatively were given by using flowfield data from the computational results. In the region before the vortex breakdown, the vortex axes are approximated as being straight line. As the angle of attack increases, the vortex axes are closer to the root chord, and farther away from the wing surface. Along the vortex axes, as the adverse pressure gradients occur, the axial velocity decreases, that is, A is negativee, so the vortex is unstable, and it is possible to breakdown. The occurrence of the breakdown results in the instability of lateral motion for a delta wing, and the lateral moment diverges after a small perturbation occurs at high angles of attack. However, after a critical angle of attack is reached the vortices breakdown completely at the wing apex, and the instability resulting from the vortex breakdown disappears.     

Numerical simulation of tangentially injected turbulent swirling flow in a divergent tube

This paper studies the properties of turbulent swirling decaying flow induced by tangential inlets in a divergent pipe using the realizable k–ε turbulence model and discusses the effects of the injector pressure and injection position. The results of transient solutions show that both the recirculation zone near the wall in upstream of the injectors and the vortex breakdown in downstream of the injectors increase in size during the whole period. A nearly axisymmetric conical breakdown is formed and its internal structure consists of two asymmetric spiral‐like vortices rotating in opposite directions. The stagnation point shifts slowly toward the pipe outlet over time. The maxima of the three velocity components, which are located near the wall, decrease gradually with streamwise direction. It can also be inferred that Mach number approaches 1.0 near the injector outlets. The velocities increase with the increasing injector pressure. However, its increasing trend is not significant. With the increase of the injection position, vortex breakdown moves in downstream direction and the pitch along the streamwise direction increases. Copyright © 2008 John Wiley & Sons, Ltd.