气固两相分支管流动及流量分配特性的研究

在水平T型分支管道中,用压缩空气对平均粒径为0．5mm砂石进行气固两相流试验。试验结果表明,当压缩空气的流速大于33m／s时,T型分支接头处没有固相沉积,两个分支管路分配的流量几乎相同。当压缩空气的流速小于33m／s时,分支接头处出现沉积,并且沉积量和分支管路的流量分配与分支管路上阀门开度有关:开度相同时,分支接头两侧的固相沉积量和流量分配相同;开度不同时,阀门开度小的一侧分支接头处的沉积量少,其分配的流量也少。     

Stability analysis of annular flow structure, using a discrete form of the “two-fluid model”

The differential form of the “two-fluid model” for annular flow, neglecting surface tension, is ill-posed, and it is not suited for examining the stability of the steady-state solutions with respect to the average film thickness. It is shown here that a discrete (difference) representation of the two-fluid model may lead to an appropriate criterion for the stability of the steady-state solutions. Exactly the same criterion is obtained from the requirement that the kinematic waves will propagate in the downstream direction. The suggested discrete form of the “two-fluid model” is used to perform transient simulation and for examining the system response to finite disturbances.     

Flow of two stably stratified fluids in an open channel with addition of fluids along the channel length

It is shown that two stably stratified fluids flowing in an open channel have two critical flow conditions. The one at higher flowrates is equivalent to the choked flow condition of a single fluid over a broad-crested weir, when the Froude number is unity. The lower critical condition imposes restrictions, which define the system if fluids are added progressively along the channel length and the flowrates increase from low to high values. However, if the flowrate does not become sufficiently large to pass through the lower critical condition, this condition will then define a form of choking, which again determines the system.It is shown that an important special case, with the proportional flowrates of the two fluids kept constant, has an analytical solution in which the relative depths of the fluids is a constant along the channel. Other systems must be solved numerically.     

Evaporation of binary mixtures in upward annular flow

A theoretical study is presented of the evaporation of binary mixtures in upward annular flow. Heat and mass balances are written and the resulting equations solved to give axial and radial variations of concentrations, temperatures and flowrates of ethanol-water mixtures. Mass and heat transfer within the film are calculated by an extension of the Dukler-Hewitt method for heat transfer in single-component films. It is concluded from the worked examples that, for the mixture considered, the film flowrates and wall temperatures are not significantly controlled by mass transfer in either phase and can be calculated by flash vaporization methods.     

A tentative approach to the direct simulation of drag reduction by polymers

An efficient technique for drag reduction uses dilute solutions of a few p.p.m. of polymers. A possible reduction in drag of up to 80% is achieved. Several experimental observations have been made which tend to indicate that the polymers modify the turbulence structures within the buffer layer. Flow visualisations have shown that the changes consist of a weakening of the strength of the streamwise vortices. Existing literature reveals no attempts of numerical simulation of this phenomenon. In this paper an approach is pursued by using a constitutive equation which relates the elongation viscosity to the local properties of the flow. According to this model this viscosity is large in zones where the amount of strain rate is greater than the amount of vorticity, and is zero when the vorticity exceeds the strain rate. Simulations have been performed in a “minimal channel” to give good resolution with a limited number of grid points. The accuracy of the method is tested by comparison with the results of other techniques. For simulations with polymers, quantitative comparisons cannot be made, but the results reproduce the qualitative outputs of the experiments. The mean streamwise velocity is modified in the buffer layer; the peak of the streamwise turbulent intensity, in wall units, increases and its maximum moves far from the wall; and the vertical turbulent intensity is largely reduced in the wall region. An interesting outcome from both the simulation and the experiments is that the strength of the longitudinal vortices is reduced when the polymers are present.     

An experimental evaluation of the behaviour of mono and polydisperse polystyrenes in Cross-Slot flow

The behaviour of a number of mono and polydisperse polystyrenes are probed experimentally in complex extensional flow within a Cross-Slot geometry using flow-induced birefringence. Polystyrenes with similar molecular weight (M _w) and increasing polydispersity (PD) illustrated the effect of PD on the principal stress difference (PSD) pattern in extensional flow. Monodisperse materials exhibited only slight asymmetry at moderate flowrates, although increased asymmetry and cusping was observed at high flowrates. The response of monodisperse materials of different M _w at various flowrates is presented and characterised by Weissenberg numbers for both chain stretch and orientation using a theory for linear entangled polymers. The comparison of stress profiles against Weissenberg number for each process is used to determine whether the PSD pattern observed is independent of M _w and elucidate which relaxation mechanism is dominant in the flow regimes probed. For monodisperse materials, at equivalent chain orientation Weissenberg number (We _τd), different molecular weight materials were seen to exhibit similar steady state PSD patterns independent of We _τR (chain stretch We). Whilst no obvious critical Weissenberg number (We) was found for the onset of increased asymmetry and cusping, it was found to occur in the “orientating flow without chain stretch” regime.     

Dividing annular flow in a horizontal tee

The prediction of mixture composition in a branch from a manifold in which two-phase mixtures flow has been examined. A linear relationship is found to exist between the branch mass flowrates of the individual phases over a range of flow conditions. This observation is used as the basis of a correlation which contains coefficients that are functions of the manifold flow condition. 90 per cent of the data are correlated to within ±20 per cent.     

Brownian dynamics simulations of star-branched polymers in dilute solutions in extensional flow

Using Brownian dynamics (BD) simulations of FENE bead-spring models, the dynamics of star-branched polymers in dilute solutions under extensional flow have been investigated. Studies on star polymers in transient extensional flow reveal that the initial transient stress response at low strains is governed by both the number of arms and the shortest arm. On the other hand, the steady-state behavior of star polymers in elongational flow is limited by the maximum effective “contour” length of the molecules. The distribution of arm extension and birefringence of the star-branched molecule are broader and the mean is shifted to lower values, when compared to equivalent linear systems. As a result, the degree of arm extension at steady-state decreases as the number of arms in the star increases. Both an analysis of individual ensembles in Brownian dynamics simulations and a study of a simple force balance indicate that the constraint imposed on the star arms by the central branch point and contributions from “asymmetric” arm arrangements give rise to overall less extended and oriented star-branched molecules with broader arm extension and birefringence distributions. The results obtained from stress-conformation hysteresis simulation indicate that less-stretched arms exhibit more retarded relaxation, as the number of arms increases in star-branched molecules. The effect of excluded volume (EV) interactions, incorporated through the Lennard–Jones potential, on the dynamics of star polymers in extensional flow appears unimportant.     

Modelling of particle-wall collisions in confined gas-particle flows

This paper demonstrates that numerical simulations of confined particulate two-phase flows require a detailed modelling of particle—wall collisions which includes the wall surface structure and the particle shape. These effects are taken into account by “irregular bouncing” models which are based on the statistical treatment of the collision process. In this study, results obtained using various “irregular bouncing” models based on the impulse equations for a particle—wall collision are considered and compared with experimental observations. The wall roughness is simulated by assuming that the particle collides with a virtual wall which has a randomly distributed inclination with respect to the plane, smooth wall. A Gaussian distribution for this random inclination showed the best agreement with experimental results. Numerical predictions of a turbulent two—phase flow in a vertical channel, where the particle phase is treated using a Lagrangian approach, showed that the different models applied for a particle-wall collision have a strong effect on the particle velocity fluctuations and the mass flux profiles in the region of fully developed flow. The numerical simulations using the irregular bouncing models yielded considerably higher values for the particle velocity fluctuations, which also agreed better with the experimental values. This effect was most pronounced for large particles, where the distance they need to respond to the fluid flow is larger than the characteristic dimension of the confinement. On the other hand, the motion of small particles is less affected by the choice of the wall-collision model. These effects of the wall roughness on the velocity fluctuations of the dispersed phase have not been considered in previous studies using irregular bouncing models.     

Investigation of a combustion driven oscillation in a refinery flare: Part B: Visualisation of a periodic flow instability in a bifurcating duct following a contraction

A flow visualisation study was performed to investigate a periodic flow instability in a bifurcating duct within the tip of the flares at the Shell refinery in Clyde, NSW, to verify the trigger of a combustion-driven oscillation proposed in Part A of this study, and to identify its features. The model study assessed only the flow instability, uncoupled from the acoustic resonance and the combustion that are also present in the actual flare. Three strong, coupled flow oscillations were found to be present in three regions of the fuel line in the flare tip model. A periodic flow separation was found to occur within the contraction at the inlet to the tip, a coupled, periodic flow oscillation was found in the two transverse “cross-over ducts” from the central pipe to the outer annulus and an oscillating flow recirculation was present in the “end-cap” region of the central pipe. The dimensionless frequency of these oscillations in the model was found to match that measured in the full-scale plant for high fuel flow rates. This, and the strength of these flow oscillations, gives confidence that they are integral to the full-scale combustion-driven oscillation and most likely the primary trigger. The evidence indicates that the periodic flow instability is initiated by the separation and roll-up of the annular boundary layer at the start of the contraction in the fuel section of the flare tip. The separation generates an annular vortex which interacts with the blind-ended pipe downstream, leading to a pressure wave which propagates back upstream, initiating the next separation event and repeating the cycle. The study also investigated flow control devices with a view to finding a practical approach to mitigate the oscillations. The shape of these devices was constrained to allow installation without removing the tip of the flare. This aspect of the study highlighted the strength and nature of the coupled oscillation, since it proved to be very difficult to mitigate the oscillation in this way. An effective configuration is presented, comprising of three individual components, all three of which were found to be necessary to eliminate the oscillation completely.     

Macroscopic properties of a two-phase potential dispersion composed of identical unit cells

The potential flow solution for flow of fluid past dispersed objects in a “unit cell” is used to derive several macroscopic properties, including the mean pressures in the phases and on the walls, the momentum and kinetic energy density, the force function and mechanical energy flux. These properties are derived from the “resistivity” of the unit cell, which has a tensorial character in general. Various macroscopic forms of Bernoulli's equation relate the properties. Equations of motion for uniform arrays of cells are derived. Various other features, such as minimization of kinetic energy density and forces at concentration jumps, are analyzed.     

Flow of two liquids in a helix: an analogue of high pressure helical boilers

It has been shown previously (Gardner & Kubie 1976) that a two-liquid system provides a better analogue for investigating certain aspects of high pressure boiler tube hydrodynamics than air and water at atmospheric pressure. A study is presented here of the flow of $\text{n-}\text{butyl}$ acetate and water, and iso-amyl alcohol and water in helical coils and the general conclusion is that strongly stratified conditions of either completely separated phases or of drops of one phase in the other exist until velocities are high enough to produce drops sufficiently fine to be dispersed by the turbulence. The critical condition for the breakdown of the separated phases to produce drops is quantitatively correlated with the mass flowrates needed to avoid low quality dryout in both high pressure water and high pressure freon helical boilers. Some evidence is available to support two criteria indicating that low quality dryout will not be observed for sufficiently low pressures in steam-water systems and, perhaps, for small enough bore tubes.     

Transition to escape in a system of coupled oscillators

We study a Hamiltonian system of coupled oscillators derived from two forced pendulums, connected with a torsional spring. The uncoupled limit is described by two identical oscillators, each possessing a homoclinic orbit separating bounded from unbounded motion. We focus on intermediate energy levels which lead to detained motions, defined as trajectories that, though unbounded as t → ∞, oscillate within the region defined by the homoclinic orbit of the unperturbed system for a long but finite time. We analyze the existence and behavior of these motions in terms of equipotential surfaces. These curves provide bounds on the motion of the system and are shown to be closed for low energies. However, above some critical energy level the equipotential curves become open. The detained trajectories are shown to arise from the region of phase space that was, for appropriate energies, stochastic. These motions remain within this region for long times before finally “leaking out” of the opening in the equipotential curves and proceeding to infinity.     

Effects of a countercurrent gas flow on falling-film mode transitions between horizontal tubes

A liquid film falling between horizontal tubes is known to take the form of droplets, jets or sheets, depending on the liquid flow rate; the form of the flow is the so-called “falling-film mode”. Although previously neglected in studies of mode transition, a countercurrent gas flow often exists in falling-film heat exchangers, and its effect on the liquid flow might be important: it could impact the flow regime, lead to local “dryout,” and decrease the heat transfer rate. Experiments are conducted to explore the effects of a countercurrent gas flow and liquid feeding length on falling-film mode transitions for a liquid flowing over horizontal tubes. The effects on mode transition are shown to depend on fluid properties and are explained in terms of unsteadiness and film thickness. In general, transition hysteresis is reduced with an increasing gas velocity. A correlation is developed to predict the countercurrent gas flow effects on falling-film mode transitions. The liquid feeding length can affect mode transitions in quiescent surroundings and when a countercurrent gas flow imposed.     

Correlations predicting frictional pressure drop and liquid holdup during horizontal gas-liquid pipe flow with a small liquid holdup

Experimental data and correlations available in the literature for the liquid holdup ε_L and the pressure gradient ΔP_TP/L for gas-liquid pipe flow, generally, do not cover the domain 0 < ε_L < 0.06. Reliable pressure-drop correlations for this holdup range are important for calculating flow rates of natural gas, containing traces of condensate. In the present paper attention is focused on reliable measurements of ε_L and ΔP_TPIL values and on the development of a phenomenological model for the liquid-holdup range 0 < ε_L < 0.06. This model is called the “apparent rough surface” model and is referred to as the ARS model. The experimental results presented in this paper refer to air-water and air-water + ethyleneglycol systems with varying transport properties in horizontal straight smooth glass tubes under steady-state conditions. The holdup and pressure gradient values predicted with the ARS model agree satisfactorily with both our experimental results and data obtained from the literature referring to small liquid-holdup values 0 < ε_L < 0.06. Further, it has been shown that in the domain 38 < < 72 mPa m the interfacial tension of the gas-liquid system has no significant effect on the liquid holdup. The pressure gradient, however, increases slightly with decreasing surface tension values.     

Phase-resolved characterization of vortex shedding in the near wake of a square-section cylinder at incidence

The vortex formation and shedding process in the near wake region of a 2D square-section cylinder at incidence has been investigated by means of particle image velocimetry (PIV). Proper orthogonal decomposition (POD) is used to characterize the coherent large-scale flow unsteadiness that is associated with the wake vortex shedding process. A particular application of the POD analysis is to extract the vortex-shedding phase of individual velocity fields, which were acquired at asynchronous low rate with respect to the vortex shedding cycle. The phase of an individual flow field is determined from its projection on the first pair of POD modes, allowing phase averaging of the measurement data to be performed. In addition, a low-order representation of the flow, constructed from the mean and the first pair of POD modes, is found to be practically equivalent to the phase-averaged results. It is shown that this low-order representation corresponds to the basic Fourier component of the flow field ensemble with respect to the reconstructed phase. The phase-averaged flow representations reveal the dominant flow features of the vortex-shedding process and the effect of the angle of incidence upon it.     

Simulations of dispersed turbulent multiphase flow

Direct numerical simulations of homogeneous isotropic turbulence are used to investigate the effects of turbulence on the transport of particles in gas flows or bubbles in liquid flows. The inertia associated with the bubbles or the particles leads to locally strong concentrations of these in regions of instantaneously strong vorticity for bubbles or strain-rate for particles. This alters the average settling rates and other processes. If the mass-loading of the dispersed phase is significant a random “turbulent” flow is generated by the particle settling. A simple demonstration of this is given, showing the statistically axisymmetric character of this flow and how it can modify an ambient turbulent flow.     

The interplay between boundary conditions and flow geometries in shear banding: Hysteresis, band configurations, and surface transitions

We study shear banding flows in models of wormlike micelles or polymer solutions, and explore the effects of different boundary conditions for the viscoelastic stress. These are needed because the equations of motion are inherently non-local and include “diffusive” or square-gradient terms. Using the diffusive Johnson–Segalman model and a variant of the Rolie-Poly model for entangled micelles or polymer solutions, we study the interplay between different boundary conditions and the intrinsic stress gradient imposed by the flow geometry. We consider prescribed gradient (Neumann) or value (Dirichlet) of the viscoelastic stress tensor at the boundary, as well as mixed boundary conditions in which an anchoring strength competes with the gradient contribution to the stress dynamics. We find that hysteresis during shear rate sweeps is suppressed if the boundary conditions favor the state that is induced by the sweep. For example, if the boundaries favor the high shear rate phase then hysteresis is suppressed at the low shear rate edges of the stress plateau. If the boundaries favor the low shear rate state, then the high shear rate band can lie in the center of the flow cell, leading to a three-band configuration. Sufficiently strong stress gradients due to curved flow geometries, such as that of cylindrical Couette flow, can convert this to a two-band state by forcing the high shear rate phase against the wall of higher stress, and can suppress the hysteresis loop observed during a shear rate sweep.     

Application of the ultrasonic technique and high-speed filming for the study of the structure of air–water bubbly flows

Multiphase flows are very common in industry, oftentimes involving very harsh environments and fluids. Accordingly, there is a need to determine the dispersed phase holdup using noninvasive fast responding techniques; besides, knowledge of the flow structure is essential for the assessment of the transport processes involved. The ultrasonic technique fulfills these requirements and could have the capability to provide the information required. In this paper, the potential of the ultrasonic technique for application to two-phase flows was investigated by checking acoustic attenuation data against experimental data on the void fraction and flow topology of vertical, upward, air–water bubbly flows in the zero to 15% void fraction range. The ultrasonic apparatus consisted of one emitter/receiver transducer and three other receivers at different positions along the pipe circumference; simultaneous high-speed motion pictures of the flow patterns were made at 250 and 1000 fps. The attenuation data for all sensors exhibited a systematic interrelated behavior with void fraction, thereby testifying to the capability of the ultrasonic technique to measure the dispersed phase holdup. From the motion pictures, basic gas phase structures and different flows patterns were identified that corroborated several features of the acoustic attenuation data. Finally, the acoustic wave transit time was also investigated as a function of void fraction.     

Reflection of shock wave from a compression corner in a particle-laden gas region

We investigated in this paper the progression of a shock-wave reflected from a compression corner in a particle-laden gas medium using a TVD class numerical technique and a MacCormack scheme. For a gas-only flow, the numerical results agreed well with the existing experimental data, suggesting that the gas phase is correctively solved. The effect of particle size and mass fraction ratio is investigated for a dilute gas-particle flow. It has been shown that the shock-wave diffraction and the flow configuration after the shock can become remarkably different from the gas-only flow depending on the particle parameters. Relaxation phenomenon due to the momentum drag and the heat exchange between the gas and the particle phases is explained.Graduate Student of Korea Advanced Institute of Science and TechnologyThis article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.