Transient viscoelastic helical flows in pipes of circular and annular cross-section

Start-up helical flows for Oldroyd-B and upper-convected Maxwell fluids are studied in straight pipes of circular and annular cross-section. The differential form of the constitutive equation leads to partial differential equations which are second-order in space and time. Apart from the condition that the fluid is initially at rest another initial condition is required to complete the solution process. By comparing results derived from the integral form of the constitutive equation we show that an appropriate initial condition may be found. Numerical results for start-up rotational flow in pipes of annular cross-section are presented.     

Studies on two-phase co-current air/non-Newtonian shear-thinning fluid flows in inclined smooth pipes

In this work, co-current flow characteristics of air/non-Newtonian liquid systems in inclined smooth pipes are studied experimentally and theoretically using transparent tubes of 20, 40 and 60 mm in diameter. Each tube includes two 10 m long pipe branches connected by a U-bend that is capable of being inclined to any angle, from a completely horizontal to a fully vertical position. The flow rate of each phase is varied over a wide range. The studied flow phenomena are bubbly flow, stratified flow, plug flow, slug flow, churn flow and annular flow. These are observed and recorded by a high-speed camera over a wide range of operating conditions. The effects of the liquid phase properties, the inclination angle and the pipe diameter on two-phase flow characteristics are systematically studied. The Heywood–Charles model for horizontal flow was modified to accommodate stratified flow in inclined pipes, taking into account the average void fraction and pressure drop of the mixture flow of a gas/non-Newtonian liquid. The pressure drop gradient model of Taitel and Barnea for a gas/Newtonian liquid slug flow was extended to include liquids possessing shear-thinning flow behaviour in inclined pipes. The comparison of the predicted values with the experimental data shows that the models presented here provide a reasonable estimate of the average void fraction and the corresponding pressure drop for the mixture flow of a gas/non-Newtonian liquid.     

Gas‐assisted fluid displacement in a circular tube and a rectangular channel

In this paper the amount of liquid left inside of a circular tube and a rectangular channel when displaced by another immiscible fluid are determined by solving the full creeping‐motion equations. The exact continuity of stress on the free surface is employed with a finite difference method. In order to solve the equations, the steady‐state shape of the interface is guessed and the normal stress boundary condition is dropped. The equations based on a stream function‐vorticity formulation are solved with the aid of elliptic grid generation. The computed results are compared with the experimental results of Taylor (J. Fluid Mech. 1961; 10: 161), the theoretical results of Reinelt and Saffman (SIAM J. Sci. Stat. Comput. 1985; 6: 542) and our experimental data. The computed results are in close agreement with our experimental data and those of previous investigators. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical analysis of the transient conjugated heat transfer in a circular duct with a power-law fluid

We treat numerically in this paper, the transient analysis of a conjugated heat transfer process in the thermal entrance region of a circular tube with a fully developed laminar power-law fluid flow. We apply the quasi-steady approximation for the power-law fluid, identifying the suitable time scales of the process. Thus, the energy equation in the fluids is solved analytically using the well-known integral boundary layer technique. This solution is coupled to the transient energy equation for the solid where the transverse and longitudinal heat conduction effects are taken into account. The numerical results for the temporal evolution of the average temperature of the tube wall, _av, is plotted for different nondimensional parameters such as conduction parameter, , the aspect ratios of the tube, and ₀ and the index of power-law fluid, n.     

Secondary flows and particle centrifugation in slightly tilted rotating pipes

A theoretical analysis is presented of viscous incompressible laminar flow in a pipe which rotates around an axis held at small angle with respect to its symmetry-axis. Analogous to the results of Barua and Benton [1, 2], solutions in closed-form are given for circulatory flows in the cross-sectional plane of the pipe due to Coriolis forces in combination with Hagen-Poiseuille flow through the pipe. The solutions are used to derive analytical expressions for trajectories of solid or liquid particles entrained in the gas and being subject to centrifugation and the said secondary flows. It is shown that despite centrifugation, particles can be locked into circulatory trajectories thus remaining suspended in the gas flowing through the pipe.     

3D phase field modeling of electrohydrodynamic multiphase flows

A 3D phase field model is developed to investigate the electrohydrodynamic (EHD) two phase flows. The explicit finite difference method, enhanced by parallel computing, is employed to solve the coupled nonlinear governing equations for the electric field, the fluid flow field and free surface deformation. Numerical tests indicate that an appropriate interpolation of densities within the interface is critical in ensuring numerical stability for highly stratified flows. The 3D phase field model compares well with the Taylor theory for the deformation of a single dielectric droplet in an electric field. Computed results show that the deformation of a leaky dielectric droplet in an electric field undergoes various stages before it reaches the final oblate shape. This is caused by the free charge relaxation near the fluid–fluid interface. The coalescence of four droplets in an electric field illustrates a truly 3D deformation behavior and a complex evolving fluid flow field associated with the participating droplets. The coalescence is a result of combined actions produced by the global electric force, the circulatory flows generated by the local electrohydrodynamic stress and the electrically-induced deformation. The 3D phase field model is also applied in modeling of an electrohydrodynamic patterning process for manufacturing nanoscaled structures, in which complex 3D flow structures develop as the electrically-induced deformation evolves.     

Two‐phase pressure drop caused by sudden flow area contraction/expansion in small circular pipes

Theoretical studies have been made to determine the pressure drops caused by abrupt flow area expansion/contraction in small circular pipes for two‐phase flow of air and water mixtures at room temperature and near atmospheric pressure. Two‐phase computational fluid dynamics (CFD) calculations, using Eulerian–Eulerian model (with the air phase being compressible for pipe contraction case) are employed to calculate the pressure drop across sudden expansion and contraction. The pressure drop is determined by extrapolating the computed pressure profiles upstream and downstream of the expansion/contraction. The larger and smaller tube diameters are 1.6 and 0.84 mm, respectively. Computations have been performed with single‐phase water and air, and two‐phase mixtures in a range of Reynolds number (considering all‐liquid flow) from 1000 to 12 000 and flow quality from 1.2 × 10^?3 to 1.6 × 10^?2. The numerical results are validated against experimental data from the literature and are found to be in good agreement. The expansion and contraction loss coefficients are found to be different for single‐phase flow of air and water, and they agreed reasonably well with the commonly used theoretical predictions. Based on the numerical results as well as experimental data, correlations are developed for two‐phase flow pressure drops caused by the flow area contraction as well as expansion. Copyright © 2010 John Wiley & Sons, Ltd.     

Laminar/turbulent oscillating flow in circular pipes

A two-dimensional oscillating flow analysis was conducted simulating the gas flow inside Stirling engine heat exchangers. Both laminar and turbulent oscillating pipe flow were investigated numerically for Re_max = 1,920 (Va = 80), 10,800 (Va = 272), 19,300 (Va = 272), and 60,800 (Va = 126). The results are here compared with experimental results of previous investigators. Predictions of the flow regime on present oscillating flow conditions are also checked by comparing velocity amplitudes and phase difference with those from laminar theory and quasi-steady profile. A high Reynolds number k-ε turbulence model was used for turbulent oscillating pipe flow. Finally, the performance of the k-ε model was evaluated to explore the applicability of quasi-steady turbulent models to unsteady oscillating flow analysis.     

Numerical investigation of the influence of gravity on flutter of cantilevered pipes conveying fluid

We have carried out a numerical investigation of the three dimensional nonlinear dynamics of a cantilevered pipe conveying fluid in the presence of gravity. The pipe may be misaligned at the clamped end with respect to gravity, and the effects of this misalignment are the main objects of the present investigation. The problem has been formulated using the Cosserat rod model. First, we have computed the equilibrium solutions and used them to experimentally validate both the Cosserat model and the constitutive law. Then, we have analyzed the occurrence of flutter, via Hopf bifurcation, for critical values of the relevant parameters of the problem, such as fluid to total mass ratio, dimensionless flow rate, dimensionless gravity and misalignment angle. The influence of the equilibrium solution on flutter has been explored, and the results of the linear stability analysis show that the stabilizing or destabilizing effect of fluid flow, either in or out of the plane of the pipe, depend crucially on the misalignment. We have also computed the non-linear periodic behavior after flutter instability by two different methods: the first one is by solving the full nonlinear equations by direct integration in time and space, while the second one is by assuming the time dependence given by an appropriate ansatz. Circular periodic orbits have then been studied and found that its loss of stability via Hopf bifurcation gives rise to stable planar periodic orbits. Finally, we have also computed the multiply periodic and chaotic behaviors which take place for sufficiently large values of the flow rate.     

Periodic flows of pseudoplastic fluids in pipes

Experimental results on the flow enhancement during flow of clay slurries (n = 0.15) through oscillating pipes are compared with theoretical predictions. The agreement is fairly good, especially in the oscillating boundary layer flow regime. Flow enhancement of order 10⁸ was found with the slurries used.     

The Giesekus fluid in ω–D form for steady two-dimensional flows

Using the fact that for simple fluids the most general constitutive equation in constant stretch history flows for the extra stress tensor τ is known in an explicit form, the Giesekus fluid model is cast into this (ω–D) form for two-dimensional flows. The three material functions needed to characterize τ are listed. The explicit results for simple shear and planar elongation reveal that the parameter α should be restricted to values less than 0.5. It is demonstrated that in this explicit form the constitutive equation is free from thermodynamic objections and can thus be used as a starting point for numerical calculations of general, but steady, two-dimensional flows. Received: 9 November 1998 Accepted: 20 May 1999     

Very slow flow of Bingham viscoplastic fluid around a circular cylinder

Numerical simulations have been used to study the flow of a Bingham viscoplastic fluid around a circular cylinder in an infinite medium with negligible inertia effects. Papanastasiou's regularisation technique has been adopted to approximate the model. The case corresponding to preponderant plasticity effects has been particularly studied and convergence of the solutions examined in detail. The flow kinematics and stresses have been determined. The rigid zones have been identified and characterised. At large Oldroyd numbers, when plasticity effects become preponderant, a viscoplastic boundary layer appears around the cylinder. The characteristics of this viscoplastic boundary layer are quantified. The results are compared with existing theoretical results, concerning particularly the predictions of the viscoplastic boundary layer theory and the plasticity theory.     

Interaction between two circular cylinders in slow flow of Bingham viscoplastic fluid

The interaction between two circular cylinders was studied in the slow flow of a Bingham viscoplastic fluid in an infinite medium without any inertia effects. The configuration studied is that in which the flow direction is parallel to the centre line of the cylinders. Finite-element numerical simulations were used with an approximation by Papanastasiou's regularisation method. The case of high yield stress effect was particularly examined. The convergence of the solutions was examined in detail. Changes in the rigid zones, kinematics and stresses were determined in relation to the degree of interaction, which is a function of the distance between the cylinders and the effect of yield stress. The results compared with the case of a single cylinder show that yield stress reduces interaction effects. The transition between configurations with interacting cylinders and configuration with isolated cylinders was examined as a function of the effect of yield stress. Correlations were proposed for the drag coefficient and the stability criterion when the cylinders are interacting.     

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

考虑颗粒间碰撞的气固两相流拉格朗日模拟

在均匀,稳定的各向同性气固两相紊流场颗粒弥散的拉格朗日模拟计算方法基础上,进一步考虑了流场中颗粒之间的碰撞对于模拟计算结果的影响。与Lavieville用大涡模拟所做的计算结果进行了对比,以对本方法进行验证,并考察了颗粒间的碰撞分别对流体相和颗粒相的影响。     

Spanwise correlations in the wake of a circular cylinder and a trapezoid placed inside a circular pipe

The spanwise correlation of a circular cylinder and a trapezoidal bluff body placed inside a circular pipe in fully developed turbulent regime is studied using hotwire anemometer. The present configuration possesses complex fluid structure interaction owing to the following features: high blockage effect; low aspect ratio of the body; upstream turbulence and interaction of axisymmetric flow with a two dimensional bluff body. The spatial correlation of such configuration is seldom reported in the literature. Results are presented for Reynolds number of Re_D=1×10⁵. Three different blockage ratios (0.14, 0.19 and 0.28) are considered in the present study. Correlation coefficient is observed to improve with increase in blockage ratio. Compared to a circular cylinder, a trapezoidal bluff body possesses high correlation length. The near wall effects tend to increase the phase drift, which is reflected in low correlation coefficients close to the pipe wall. The results show that the simultaneous effect of curvature, low aspect ratio and upstream turbulence reduces the correlation coefficients significantly as compared to unconfined and confined (parallel channel) flows. The low frequency modulations with a circular cylinder are higher for lower blockage ratios. The three-dimensionality of vortex shedding for trapezoid with a blockage ratio of 0.28 was observed to be lower compared to circular cylinder and all other blockage ratios. Low frequency modulations were found to be responsible for weak vortex shedding from a circular cylinder compared to a trapezoidal bluff body. The vortex shedding is observed to be nearly two dimensional in case of a trapezoidal bluff body of blockage ratio 0.28.     

An analytical model for plug flow in microcapillaries with circular cross section

Plug flow in microcapillaries or microchannels offers significant advantages for the development of microfluidic applications and recently triggers many interests and studies. Recirculation is formed within liquid plugs due to the presence of interfaces. This paper presents an analytical model to investigate the recirculation flow and the flow resistance in microcapillaries with circular cross section. A fourth order partial differential equation is used to model the Stokes flow within the liquid plug. The results of the flow field show that the flow pattern is affected by the plug length. The flow resistance is determined through the force balance of the liquid plug. The comparison of the flow field and the flow resistance from the analytical model and the experiments shows good agreement.     

Eulerian prediction of the fluid/particle correlated motion in turbulent two-phase flows

An approximate equation governing the turbulent fluid velocity encountered along discrete particle path is used to derive the fluid/particle turbulent moments required for dispersed two-phase flows modelling. Then, closure model predictions are compared with results obtained from large-eddy simulation of particle fluctuating motion in forced isotropic fluid turbulence.