Drag Coefficient and Strouhal Number Analysis of a Rectangular Probe in a Two-Phase Cross Flow
Total Page:16
File Type:pdf, Size:1020Kb
Energy Equip. Sys./ Vol. 6/No. 1/March 2018/ 7-15 Energy Equipment and Systems http://energyequipsys.ut.ac.ir www.energyequipsys.com Drag coefficient and strouhal number analysis of a rectangular probe in a two-phase cross flow Authors ABSTRACT Erfan Kosaria* Ali Rahnamab In some case of laboratory and industrial applications, various c kind of measurement instruments must be placed in a conduit, in Mahyar Momen b which multiphase fluid flows. Vortex shedding for any immersed Pedram Hanafizadeh body in flow field is created with a frequency, which according to Mohammad Mahdi Rastegardoost b flow conditions such as flow rates, geometry of body, etc. may be a constant or variable. Failure may happen, if this frequency is close Mechanical Engineering Department, University of to one of the natural frequencies of the instruments. These flows California, Riverside, CA, USA can play a significant role in long-term reliability and safety of b School of Mechanical Engineering, University of industrial and laboratory systems. In this study, an Eulerian– Tehran, Tehran, Iran Eulerian approach is employed to simulate Air-Water two-phase flow around a rectangular probe with different volume fractions c School of Mechanical Engineering, KNT, (0.01-0.5) and Reynolds numbers (1000-3000). Two-phase flow University of Technology, Tehran, Iran characteristics around the probe have been analyzed numerically. The results show vortex shedding in all cases with distinct Strouhal number. In addition, results illustrate that shedding is intensified by increasing Reynolds number. In order to validate the results, fraction of inlet volume was set to zero, and drag coefficient and its relation with low Reynolds number Article history: (1000-3000) in single phase flow were compared to experimental Received : 31 July 2017 and numerical results in published article. The results show a Accepted : 6 January 2018 complete agreement between the simulation and available data. Keywords: Two-Phase Cross Flow, Strouhal Number, Rectangular Probe, Drag Coefficient. 1. Introduction Air/Water multiphase flows occur in a wide interest because many civil but also industrial range of natural and man-made situations. structure can be assimilate to this shape. From Examples include boiling heat transfer and an engineering point of view, the prediction cloud cavitation, bubble columns, cooling whether a structure has a potential to experience circuits of power plants, spraying of liquid fuel damaging vibrations or not is very important. In and paint, emulsions, rain, bubbles and drops addition, the study of bluff body wakes and its due to wave breaking in the oceans, and associated fluid forces and wake frequencies is explosive volcanic eruptions. of fundamental interest. In many flow systems, Attraction of multi-phase flows draws many complex conduits with flow obstacles have attempts to investigate unexplored aspects of been widely employed, and the flow in them is this field experimentally and numerically [1- often Air-Water. However, only a few studies 27]. The study of flow around bodies of of the two-phase drag experienced by such rectangular shape has a deep engineering rectangular bodies have been reported in the 8 Erfan Kosari et al. / Energy Equip. Sys. / Vol. 6/No. 1/March 2018 literature, because of its complexity and Artemiev and Kornienko, developed a heterogeneity. Knowledge of the drag numeric model to predict the manner of coefficient for a rectangular body, the volume effecting the shape of volume fraction profiles fraction distribution and bubble behavior near a on the distribution of temperature and liquid- probe subjected to two-phase cross flow is phase velocity in vertical flows [3]. Their model essential for development of mechanistic can be applied both to downward and upward understanding of fluid elastic instability of flows. They investigate how gas bubbles probes subjected to two-phase cross flow. The present in the flow improve the turbulence. The main objective of this study is to present turbulent viscosity was shown as a linear accurate and reliable data on the combined combination of two terms, one term being due effects of Reynolds number and volume to the liquid-phase turbulence, which the fraction on rectangular cylinders. Reichardt formula can be used for its calculation, and the other term modeling the 1.1. Experimental Works additional viscosity due to the relative gas motion. The second term has empirical Many works have been performed experimentally constants that making the developed numerical in the field of two phase flow. In particular, model specified in any circumstance. Mizota & Okajima [1], using a tandem type of Zaichik et al. [4], proposed a diffusion- hot-wire probe, measured successfully flow inertia model for the transport of low-inertia patterns around various rectangular cylinders particles of arbitrary density. The predicted data including a reversed flow region close to the and experimental data reported by Kamp et al. cylinders, and confirmed that the changes of the [5] were compared for gas–liquid flows in flow patterns have close correlations with those of vertical pipes (the cases of downward and drag and lift forces and Strouhal numbers in the upward flow under various conditions of the variation of the width-to-height ratio and the angle effect due to the gravity force). In principle, it of incidence. was shown to be possible to analyze bubbly Yokosawa et al. [2] worked on a single tube flows with the help of the diffusion-inertia and measured its drag coefficient under two-phase model, initially developed for gas-dispersed cross flow in the Reynolds numbers range of flows. The 푘 − 휖 model of turbulence was used 4000-300,000 and for low void fractions (0-0.1). for single-phase flow to calculate the liquid They found that the drag coefficient decreased turbulence. with increasing volume fraction for two-phase A numerical model for transport phenomena Reynolds numbers sufficiently below the single- in laminar upward bubbly flow was developed phase critical Reynolds number. However, for by Antal et al. [6]. The model was based on the Reynolds numbers above the critical value, the Eulerian two-velocity approach. The wall force drag coefficient was found to gradually increase and the lifting (Saffman) force, were the major with increasing volume fraction, although it forces acting on bubbles in the laminar flow. stayed significantly less than the value for Lopez et al. [7] reported a numerical and subcritical, single-phase flow. They also classified experimental study in which the turbulent flow patterns of the two-phase wake flow behind bubbly flow in a triangular channel was a cylinder and to investigate quantitatively the analyzed. Carrica et al. [8] and Palitano et al. change of the drag coefficient corresponding to [9] adressed the case of poly-dispersed two- the transition of the flow patterns. The former phase flows. The model makes it possible to discussed the fluctuation characteristics in surface allow for the shift of the bubble concentration pressure, lift and drag forces on a cylinder in the maximum from the near-wall zone towards the two-phase or single phase cross flow. flow core observed when the dispersed-phase size increases above some critical value. 1.2. Numerical Works Troshko and Hassan [10] developed a There are several numerical simulation studies numerical model involving the law of the wall for two-phase flows. For the most part, these for vertical mono-dispersed bubbly flow in a studies addressed the case of upward flow of pipe. two-phase gas–liquid systems. But Artemiev In Lopez et al. [7], Carrica et al. [8], Palitano and Kornienko [3] and Zaichik et al. [4] et al. [9], Troshko and Hassan [10] and Lee et reported some results for downward gas–liquid al. [11] to predict the liquid-phase turbulence, a flows. two-equation model of turbulence extended to the case of two-phase flow, was used. It can be Erfan Kosari et al. / Energy Equip. Sys. / Vol. 6/No. 1/March 2018 9 stated that, in spite of the intensive recent It is complicated to numerically simulate the efforts aimed at numerical investigation of multiphase flow using the Eulerian. It solves a downward gas–liquid flows, gained data cover set of n continuity and momentum equations for only a narrow range of parameters, and phases. Coupling through the pressure and therefore apply only to particular conditions. interphase exchange coefficients, is achieved. The purpose of the present study is a Types of involving phases, effect on the manner numerical investigation of the structure of cross in which coupling is handled; Handling of Air-Water flow over selected body. A model granular (fluid-solid) flows are different than constructed around the Eulerian representation non-granular (fluid-fluid) flows. Type of for both phases was developed to numerically mixture that is modeled, affect the momentum examine the structure of flow field and wake exchange between the phases. Particle zone behind the probe in mono-dispersed Air- suspension, fluidized beds, bubble columns, Water flow. and risers, are the applications of the Eulerian multiphase model. 2. Nomenclature 3.2. Conservation of mass D Diameter (m) F Force (N) The continuity phase averaged equation for Interphase momentum exchange phase 푞 is: K coefficient p Pressure (Pa) (1) αq ρ q . α q ρ q V q 0 t Time (s) t v Velocity (푚⁄푠) 3.3. Conservation of momentum V Volume of phase (푚3) f Frequency (Hz) Neglecting laminar stress-strain tensor, gravity ReTP Two-phase Reynolds number and other body forces, the momentum balance St Strouhal number for phase q yields: Greek symbols αq ρ qV.VV q α q ρ q q q α Volume fraction t (2) μ Dynamic viscosity (푃푎. 푠) α p. τFt ρ Density (퐾푔⁄푚3) q q pq τ Stress-strain tensor (Pa) The tilde indicates phase-averaged variables Subscribe while an over bar reflects time-averaged values.