Analytical solutions of the helical flow of non newtonian fluid in eccentric annular space

ＡＮＡＬＹＴＩＣＡＬＳＯＬＵＴＩＯＮＳＯＦＴＨＥＨＥＬＩＣＡＬＦＬＯＷＯＦＮＯＮＮＥＷＴＯＮＩＡＮＦＬＵＩＤＩＮＥＣＣＥＮＴＲＩＣＡＮＮＵＬＡＲＳＰＡＣＥＺｈａｎｇＨａｉ－ｑｉａｏ（张海桥）ＷｕＪｉ－ｚｈｏｕ（吴继周）（ＴｈｅＤｅｐａｒｔｍｅｎｔｏｆ...     

A mathematical approach to study fluid-coupled vibration of eccentric annular plates

A mathematical method is proposed to study fluid-coupled vibration of axisymmetric plate structures with asymmetries due to either imperfection or practical reasons, e.g. the weight reduction of structure, natural frequency shifting, and accessibility. The suggested approach makes use of the separation of variables to determine general solutions of the partial differential equation of the plate transverse displacement, whilst defining multiple polar coordinate systems, each of which offers its own formulation of the plate deformation with respect to its coordinate system. Moreover, closed-form geometric equations and the chain rule for determining derivatives are implemented to move from one coordinate system to the other in order to satisfy boundary conditions. The mode shapes of the vibrating plate in the dry condition are determined and in turn used in the Rayleigh–Ritz method to characterize vibrational properties of the fluid-coupled plate structure. While implementing such an energy method, the fluid motion is formulated employing the velocity potential and solved using the separation of variables. Fluid–structure interaction is also taken into account satisfying the compatibility condition on the fluid–plate​ interface. The developed methodology to predict natural frequencies has been validated by comparison with results obtained by a commercial finite element program. It is also found that the eccentricity tends to reduce natural frequencies of the fluid-coupled system for the lower serial mode, but increases them for the higher serial modes regardless of the presence of liquid.     

Magnetohydrodynamic pipe flow in ducts with partial circular ring cross section and annular cross section

In this paper we use the Green function method to solve the problem of steady one dimensional flow of an incompressible viscous, electrically conducting fluid through a pipe with partial circular ring cross section and one with annular cross section, in the presence of an applied transverse uniform magnetic field. We obtain analytic solutions and carry out some numerical calculations of the velocity distribution and induced magnetic field.     

Unsteady rotative flow of non-Newtonian fluid in an annular pipe

ＵＮＳＴＥＡＤＹＲＯＴＡＴＩＶＥＦＬＯＷＯＦＮＯＮ－ＮＥＷＴＯＮＩＡＮＦＬＵＩＤＩＮＡＮＡＮＮＵＬＡＲＰＩＰＥＨｕａｎｇＪｕｎ－ｑｉ（黄军旗）（Ｄｅｐｔｏｆ＇Ｒｅｓｏｕｒ＆Ｅｎｖｉ．Ｓｃｉ．Ｂｅｉｊｉｎｇ．ＮｏｒｍａｌＵｎｉｖｅｒｓｉｔｙＢｅｉｊｉｎ...     

Stability of annular flow and slugging

In this work we propose an effective viscosity criterion for the stabilization of annular gas-liquid and liquid-particle flows and an inertial mechanism which drives waves into slugs in slugging gas-liquid flows. Annular flow is stable when the fluid having the higher effective viscosity occupies the core region and the lower viscosity fluid is in the annulus. The eddy viscosity criterion is shown to be very consistent with published work on annular flow transitions in horizontal and vertical gas-liquid flows. It also applies to a variety of liquid-solid and gas-solid flows. In the second part of the paper we propose a mechanism for explaining the growth of initially small waves and initiation of slugs in gas-liquid flow.     

A boundary integral equation method for Navier-Stokes equations—application to flow in annulus of eccentric cylinders

A novel Navier-Stokes solver based on the boundary integral equation method is presented. The solver can be used to obtain flow solutions in arbitrary 2D geometries with modest computational effort. The vorticity transport equation is modelled as a modified Helmholtz equation with the wave number dependent on the flow Reynolds number. The non-linear inertial terms partly manifest themselves as volume vorticity sources which are computed iteratively by tracking flow trajectories. The integral equation representations of the Helmholtz equation for vorticity and Poisson equation for streamfunction are solved directly for the unknown vorticity boundary conditions. Rapid computation of the flow and vorticity field in the volume at each iteration level is achieved by precomputing the influence coefficient matrices. The pressure field can be extracted from the converged streamfunction and vorticity fields. The solver is validated by considering flow in a converging channel (Hamel flow). The solver is then applied to flow in the annulus of eccentric cylinders. Results are presented for various Reynolds numbers and compared with the literature.     

Interpretation of atomization rates of the liquid film in gas-liquid annular flow

For vertical gas-liquid annular flow the fraction of the liquid in the gas is controlled by the rate of atomization of the liquid film flowing along the wall and the rate of deposition of droplets entrained in the gas. Measurements of the rate of atomization are interpreted by a Kelvin-Helmholtz mechanism. Small wavelets on the liquid film are visualized to be entrained when wave-induced variations in the gas pressure cannot be counterbalanced by surface tension effects.     

Slip flow in an annular sector duct using radial eigenfunctions

The fully developed slip flow in an annular sector duct is solved by expansions of eigenfunctions in the radial direction and boundary collocation on the straight sides. The method is efficient and accurate. The flow field for slip flow differs much from that of no-slip flow. The Poiseuille number increases with increased inner radius, opening angle, and decreases with slip.     

垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

Condensation from steam—air mixtures in a horizontal annular flow channel

Experiments were performed on the condensation of steam from steam-air mixtures in annular flow at a cooled inner tube. The range of investigation was varied for laminar and turbulent flow for 1.5 × 10³ Re 1.3 × 10⁴ and inlet concentrations 0.59 p_steam/p_total 0.95. The measurements, performed at an open test loop at p_total ≈ 0.96 bar, allowed local heat and mass transfer coefficients to be evaluated for various inlet lengths in the 2 m long annulus. The steam concentration was measured locally inside the annulus with a newly developed dew-point probe. The heat flux was measured locally using the temperature gradient in the cooled inner tube.

Near the inlet region the experiments showed a slightly higher heat flux at the bottom of the tube compared to the top, although it is expected to be smaller there owing to a thicker liquid film. Far downstream from the inlet region the heat transfer at the top was higher than at the bottom. The reasons for these effects are discussed, yielding a better understanding of the thermal and fluid processes involved in condensation from vapor-gas mixtures. The measured data allow the development of correlations for predicting the local Nusselt and Sherwood numbers in a horizontal annular-flow chanbel.     

An analytical solution and analysis of characters for viscoelastic fluid flow in annular pipe

I.Introducti0nLuiCiqunandHuangJunqiI'I(l989),ZhuWeihuiandLuiCiquri1'l(l992)sequentiallystudiedtheaxialflowofsecondorderandMaxwellfluidsandanalyzedtheflowcharactersofthesefluids.Thispaperwillstudyunsteadyrotat0ryflowofsecondorderandMaxwellfluidsinannularpi…     

Prediction of the interfacial shear-stress in vertical annular flow

Many improvements of the Wallis correlation for the interfacial friction in annular flow have been proposed in the literature. These improvements give in general a better fit to data, however, their physical basis is not always justified. In this work, we present a physical approach to predict the interfacial shear-stress, based on the theory on roughness in single-phase turbulent pipe flows. Using measured interfacial shear-stress data and measured data on roll waves, which provide most of the contribution to the liquid film roughness, we show that the interfacial shear-stress in vertical annular flow is in very close agreement with the theory. We show that the sand-grain roughness of the liquid film is not equal to four times the mean film thickness, as it is assumed in the Wallis correlation. Instead, the sand-grain roughness is proportional to the wave height, and the proportionality constant can be predicted accurately using the roughness density (or solidity). Furthermore, we show that our annular flow, which is in similar conditions to others in the literature, is fully rough. Hence, the bulk Reynolds number should not appear in the prediction of the interfacial friction coefficient, as is often done in the improvements of the Wallis correlation proposed in the literature.     

On the orbital motion of a rotating inner cylinder in annular flow

In this paper, numerical calculations have been performed to analyse the influence of the orbital motion of an inner cylinder on annular flow and the forces exerted by the fluid on the inner cylinder when it is rotating eccentrically. The flow considered is fully developed laminar flow driven by axial pressure gradient. It is shown that the drag of the annular flow decreases initially and then increases with the enhancement of orbital motion, when it has the same direction as the inner cylinder rotation. If the eccentricity and rotation speed of the inner cylinder keep unchanged (with respect to the absolute frame of reference), and the orbital motion is strong enough that the azimuthal component (with respect to the orbit of the orbital motion) of the flow‐induced force on the inner cylinder goes to zero, the flow drag nearly reaches its minimum value. When only an external torque is imposed to drive the eccentric rotation of the inner cylinder, orbital motion may occur and, in general, has the same direction as the inner cylinder rotation. Under this condition, whether the inner cylinder can have a steady motion state with force equilibrium, and even what type of motion state it can have, is related to the linear density of the inner cylinder. Copyright © 2006 John Wiley & Sons, Ltd.     

Effect of secondary flow on droplet distribution and deposition in horizontal annular pipe flow

In horizontal annular dispersed pipe flow the liquid film at the bottom is thicker and rougher than at the top of the pipe. A turbulent pipe flow experiencing a variation of roughness along the pipe wall will show a secondary flow. Such secondary flow, consisting of two counter-rotating cells in the cross-section of the tube, can change the distribution of the droplets inside the pipe and their deposition at the wall. Here, we compare the behaviour of the droplets (dispersed phase) with and without secondary flow, using large-eddy simulations. It is shown that the presence of secondary flow increases the droplet concentration in the core of the pipe and the droplet deposition-rate at the top of the pipe.     

On the role of droplets in cocurrent annular and churn-annular pipe flow

Decreasing the gas flow-rate in an initially vertical upward annular dispersed pipe-flow, will eventually lead to a down-flow of liquid. The onset of this down-flow has been related in the literature to the presence of the dispersed phase and the instability of the liquid film. Here we investigate how the dispersed-phase may influence the down-flow, performing detailed PDA-measurements in a 5 cm vertical air–water annular-flow. It is shown that the dispersed-phase does not cause the liquid down-flow, but that it delays the onset of liquid down-flow. In cocurrent annular flow the dispersed phase seems to stabilise the film flow, whereas in churn-annular flow the opposite seems to be true.     

Hamilton体系下环扇形域的Stokes流动问题

基于极坐标下Stokes流的基本方程,将径向坐标模拟为时间坐标,推导了Hamilton体系下Stokes流动问题的对偶方程,采用本征向量展开法对环扇形域Stokes流动问题进行了分析,并给出了相应的实际算例,其结果说明了本文方法的有效性。     

Heat transfer to non-Newtonian fluids in laminar flow through rectangular ducts

Heat transfer to non-newtonian fluids flowing laminarly through rectangular ducts is examined. The conservation equations of mass, momentum, and energy are solved numerically with the aid of a finite volume technique. The viscoelastic behavior of the fluid is represented by the Criminale-Ericksen-Filbey (CEF) constitutive equation. Secondary flows occur due to the elastic behavior of the fluid, and, consequently, heat transfer is strongly enhanced. It is observed that shear thinning yields negligible heat transfer enhancement effect, when compared with the secondary flow effect. Maximum heat transfer is shown to occur for some combinations of parameters. Thus, there are optimal combinations of aspect ratio and Reynolds numbers, which depend on the fluid's mechanical behavior. This result can be usefully explored in thermal designs of certain industrial processes.     

Axial annular flow of plastic fluids: dead zones and plug-free flow

In axial annular flow, the shear stress decreases from its value τ(κR) at the inner cylinder to 0 at r = λR and increases from then on to τ(R) at the outer cylinder. For plastic fluids with a yield stress τ _c, λ will be such that flow commences when τ(κR) = τ(R) = τ _c. For fluids with position-dependent yield stresses (electro- and magnetorheological fluids are examples), the situation is more complex. While it is possible that yielding and flow occur everywhere, it is also possible that flow occurs only in parts of the fluid-filled space, and a dead zone (region in which the fluid is at rest) close to one of the walls exists. In that case, the fluid will flow no matter how small the applied pressure difference is. If P is large enough, the dead zone ceases to exist and flow without any plug is possible. The fluid flows as if no yield stress exists.

Computation of a drastic flow pattern change in an annular swirling jet caused by a small decrease in inlet swirl

This paper investigates the flow pattern change in an annular jet caused by a sudden change in the level of inlet swirl. The jet geometry consists of an annular channel followed by a specially designed stepped‐conical nozzle, which allows the existence of four different flow patterns as a function of the inlet swirl number. This paper reports on the transition between two of them, called the ‘open jet flow high swirl’ and the ‘Coanda jet flow.’ It is shown that a small sudden decrease of 4% in inlet swirl results in a drastic and irreversible change in flow pattern. The objective of this paper is to reveal the underlying physical mechanisms in this transition by means of numerical simulations. The flow is simulated using the unsteady Reynolds‐averaged Navier–Stokes (URANS) approach for incompressible flow with a Reynolds stress turbulence model. The analysis of the numerical results is based on a study of different forces on a control volume, which consists of the jet boundaries. The analysis of these forces shows that the flow pattern change consists of three different regimes: an immediate response regime, a quasi‐static regime and a Coanda regime. The simulation reveals that the pressure–tangential velocity coupling during the quasi‐static regime and the Coanda effect at the nozzle outlet during the Coanda regime are the driving forces behind the flow pattern change. These physical mechanisms are validated with time‐resolved stereo‐PIV measurements, which confirm the numerical simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

A separated flow model for predicting two-phase pressure drop and evaporative heat transfer for vertical annular flow

A separated flow model has been developed that is applicable to vertical annular two-phase flow in the purely convective heat transfer regime. Conservation of mass, momentum, and energy are used to solve for the liquid film thickness, pressure drop, and heat transfer coefficient. Closure relationships are specified for the interfacial friction factor, liquid film eddy-viscosity, turbulent Prandtl number, and entrainment rate. Although separated flow models have been reported previously, their use has been limited, because they were tested over a limited range of flow and thermal conditions. The unique feature of this model is that it has been tested and calibrated against a vast array of two-phase pressure drop and heat transfer data, which include upflow, downflow, and microgravity flow conditions. The agreements between the measured and predicted pressure drops and heat transfer coefficients are, on average, better or comparable to the most reliable empirical correlations. This separated flow model is demonstrated to be a reliable and practical predictive tool for computing two-phase pressure drop and heat transfer rates. All of the datasets have been obtained from the open literature.