An investigation into the interfacial shear stress contribution in two-phase stratified flow

This paper provides a combined theoretical and experimental investigation into the contribution of interfacial shear stress in certain co and counter-current flows in circular pipes. Based on momentum balance two equations were developed for such flows then two fluid systems of significantly different density ratio were experimentally tested to quantify these equations.     

Stable two-layer flows at all Re; visco-plastic lubrication of shear-thinning and viscoelastic fluids

Multi-fluid flows are frequently thought of as being less stable than single phase flows. Consideration of different non-Newtonian models can give rise to different types of hydrodynamic instability. Here we show that with careful choice of fluid rheologies and flow paradigm, one can achieve multi-layer flows that are linearly stable for Re = ∞. The basic methodology consists of two steps. First we eliminate interfacial instabilities by using a yield stress fluid in one fluid layer and ensuring that for the base flow configurations studied we maintain an unyielded plug region at the interface. Secondly we eliminate linear shear instabilities by ensuring a strong enough Couette component in the second fluid layer, imposed via the moving interface. We show that this technique can be applied to both shear-thinning and visco-elastic fluids.     

Weakly nonlinear stability analysis of a falling film with countercurrent gas flow

Weakly nonlinear stability analysis of a falling film with countercurrent gas–liquid flow has been investigated. A normal mode approach and the method of multiple scales are employed to carry out the linear and nonlinear stability solutions for the film flow system. The results show that both supercritical stability and subcritical instability are possible for a film flow system when the gas flows in the countercurrent direction. The stability characteristics of the film flow system are strongly influenced by the effects of interfacial shear stress when the gas flows in the countercurrent direction. The effect of countercurrent gas flow in a falling film is to stabilize the film flow system.     

The instability mechanism of single and multilayer Newtonian and viscoelastic flows down an inclined plane

The instability mechanism of single and multilayer flow of Newtonian and viscoelastic fluids down an inclined plane has been examined based on a rigorous energy analysis as well as careful examination of the eigenfunctions. These analyses demonstrate that the free surface instability in single and multilayer flows in the limit of longwave disturbances (i.e., the most dangerous disturbances) arise due to the perturbation shear stresses at the free surface. Specifically, for viscoelastic flows, the elastic forces are destabilizing and the main driving force for the instability is the coupling between the base flow and the perturbation velocity and stresses and their gradient at the free surface. For Newtonian flows at finite Re, the driving force for the interfacial instability in the limit of longwaves depends on the placement of the less viscous fluid. If the less viscous fluid is adjacent to the solid surface then the main driving force for the instability is interfacial friction, otherwise the bulk contribution of Reynolds stresses drives the instability. For viscoelastic fluids in the limit of vanishingly small Re, the driving force for the instability is the coupling of the base flow and perturbation velocity and stresses and their gradients across the interface. In the limit of shortwaves the interfacial stability mechanism of flow down inclined plane is the same as plane Poiseuille flows (Ganpule and Khomami 1998, 1999a, b). Received: 20 October 2000/Accepted: 11 January 2001     

Two-fluid modeling in analyzing the interfacial stability of liquid film flows

The area-averaged two-fluid model formulation of a separated two-phase flow system is used to investigate interfacial stability of liquid film flows. The analysis takes into account the effects of phase change at the interface as well as the dynamic effects of the adjacent vapor flow on the interfacial stability. Wave formation and instability criteria are established in terms of the generalized fluid and flow parameters. The criteria are applied to investigate the stability of laminar liquid film flow with interfacial shear and phase change. The influence of various dimensionless parameters characterizing film thickness, gravity, phase change and interfacial shear are studied with respect to the neutral stability, temporal growth factor and the wave propagation velocity. The results of the present study indicate that the interfacial stability analysis developed within the frame of the two-fluid model formulation proves to be quite accurate as judged by comparing its results with the available experimental data and with the results of much longer and more complex analytical investigations which are valid only for the liquid film free of interfacial shear.     

Exact solutions of core-annular laminar inclined flows

An exact solution for laminar two-phase eccentric core-annular flows (CAF) in inclined pipes is derived. This solution complements the exact solutions that were obtained for inclined stratified flows with curved interfaces as to provide a set of solutions for two-phase laminar separated flows. A unified set of three dimensionless parameters for separated flows is defined and used to explore the effects of the system parameters and separated flow configurations on the velocity profiles and the resulting holdup, pressure gradient and pumping power requirement in horizontal and inclined concurrent and countercurrent flows. It is shown that similarly to stratified flows, also in CAF multiple solutions for the holdup and the associated flow characteristics can be obtained in inclined flows. The boundaries of the multiple solution regions are mapped and the effect of the core eccentricity and system parameters boundaries are demonstrated and discussed.The benefits of adding a lubricating phase for transportation of a viscous fluid in inclined CAFs is investigated. An adverse effect of the upward pipe inclination on the power savings in all of the separate flow configurations is demonstrated. Independently of the density of the lubricant, namely, whether it is lighter or heavier than the viscous fluid, the effect of hydrostatic pressure gradient always hinders the possibility of reducing the pumping requirement for transporting the viscous phase. However, surprisingly, a heavier lubricant is preferable form the view point of power saving. The implications of turbulent flow of the lubricating phase and the susceptibility to Ledinegg instability on the potential power savings are also considered and discussed. The application of the model for the analysis of experimental data of the holdup and pressure drop obtained in horizontal and inclined CAF is also demonstrated.     

Flow pattern transition,pressure gradient,hold-up predictions in gas/non-Newtonian power-law fluid stratified flow

This work focuses on gas/non-Newtonian power-law fluid stratified pipe flow. Two different theoretical approaches to obtain pressure gradient and hold-up predictions are presented: the steady fully developed two-fluid model and the pre-integrated model. The theoretical predictions are compared with experimental data available for horizontal and for slightly downward inclined air/shear thinning fluid stratified flow taken from literature. The predictions of the pre-integrated model are validated showing a good agreement when compared with experimental data. The criteria for the transition from the stratified flow pattern are applied to gas/non-Newtonian stratified flow. The neutral stability analysis (smooth/wavy stratified flow) and the well-posedness (existence region of stratified flow) of governing equations are carry out. The predicted transition boundaries are obtained using the steady fully developed two-fluid model and the pre-integrated model, where the shape factors and their derivatives are accounted for. A comparison between the predicted boundaries and experimental flow pattern maps is presented and shows a good agreement. A comment on the shear stress modeling by the pre-integrated model is provided.     

Analysis of stability of a thin layer of granular material moving on an inclined plane

The stability of a layer of a granular medium on an inclined plane has been studied within the framework of the model of a non-Newtonian fluid with an index of 2, which ensures the experimentally found quadratic dependence of the shear stress on the shear rate. It is shown analytically and numerically that these flows are stable or unstable depending on the value of the generalized Reynolds number relative to the critical value equal to5 cot α. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 6, pp. 113–117, November–December, 1998.     

Gas/shear-thinning liquid flows through pipes: Modeling and experiments

In chemical and oil industry gas/shear-thinning liquid two-phase flows are frequently encountered. In this work, we investigate experimentally the flow characteristics of air/shear-thinning liquid systems in horizontal and slightly inclined smooth pipes. The experiments are performed in a 9-m-long glass pipe using air and three different carboxymethyl cellulose (CMC) solutions as test fluids. Flow pattern maps are built by visual observation using a high-speed camera. The observed flow patterns are stratified, plug, and slug flow. The effects of the pipe inclination and the rheology of the shear-thinning fluid in terms of flow pattern maps are presented. The predicted existence region of the stratified flow regime is compared with the experimental observation showing a good agreement. A mechanistic model valid for air/power-law slug flow is proposed and model predictions are compared to the experimental data showing a good agreement. Slug flow characteristics are investigated by the analysis of the signals of a capacitance probe: slug velocity, slug frequency, and slug lengths are measured. A new correlation for the slug frequency is proposed and the results are promising.     

Cohesive zone modeling of interfacial stresses in plated beams

The elastic analysis of interfacial stresses in plated beams has been the subject of several investigations. These studies provided both first-order and higher-order solutions for the distributions of interfacial shear and normal stresses close to the plate end in the elastic range. The notable attention devoted to this topic was driven by the need to develop predictive models for plate end debonding mechanisms, as the early models of this type adopted debonding criteria based on interfacial stresses. Currently, approaches based on fracture mechanics are becoming increasingly established. Cohesive zone modeling bridges the gap between the stress- and energy-based approaches. While several cohesive zone analyses of bonded joints subjected to mode-II loading are available, limited studies have been conducted on cohesive zone modeling of interfacial stresses in plated beams. Moreover, the few available studies present complex formulations for which no closed-form solutions can be found. This paper presents an analytical cohesive zone model for the determination of interfacial stresses in plated beams. A first-order analysis is conducted, leading to closed-form solutions for the interfacial shear stresses. The mode-II cohesive law is taken as bilinear, as this simple shape is able to capture the essential properties of the interface. A closed-form expression for the debonding load is proposed, and the comparison between cohesive zone modeling and linear-elastic fracture mechanics predictions is discussed. Analytical predictions are also compared with results of a numerical finite element model where the interface is described with zero-thickness contact elements, using the node-to-segment strategy and incorporating decohesion and contact within a unified framework.     

The roles of inertia and shear-thinning in flow of an inelastic liquid through an axisymmetric sudden contraction

The roles of luid inertia and shear-rate dependent viscosity in determining the velocity field in an axisymmetric sudden contraction are assessed by finite-element analysis for a generalized Newtonian fluid with viscosity function given by a Carreau equation. Acting alone, either increasing shear-thinning of the viscosity or increasing fluid inertia suppresses the upstream vortex that surrounds the opening to the small tube. For creeping flows, shear thinning does not produce concavities and off-centered maxima in the axial velocity profile just inside the small tube, even at high Carreau numbers where the velocity field approaches the limiting form for a power-law fluid. Peaks in the axial velocity away from the center of the tube were found only for moderate and high Reynolds numbers and were enhanced by shear thinning, which decreased the viscosity and consequently increased the “local” Reynolds number near the wall of the small tube. The effect of steep velocity gradients near this surface on the accuracy of the finite-element approximations is discussed.     

Solutions of xanthan gum/guar gum mixtures: shear rheology,porous media flow,and solids transport in annular flow

Mixtures of xanthan and guar gum in aqueous solution were studied in two flow situations: simple shear and porous media. In addition, solids transport in vertical annular flow of sand suspensions was explored. The zero shear rate viscosity of the solutions displayed a pronounced synergy: the viscosity of the mixture is higher than that of the polymer solutions in a wide range of relative concentrations of the two polymers, in agreement with previous literature. However, at relatively high shear rates, the viscosity approaches the value of the more viscous xanthan gum solutions at mass fractions of xanthan gum between 0.1 and 0.15, and the degree of synergy substantially decreases. Stress relaxation experiments in simple shear indicate that the polymer mixtures exhibit a well-defined yield stress after relaxation that is absent in solutions of pure polymers. In porous media flow experiments, a synergistic behavior mimicking the shear flow results was obtained for the polymer mixtures at low shear rates. However, at a critical shear rate, the apparent viscosity in porous media flows exceeds the shear viscosity due to the elongational nature of flow in the pores. The solids transport capacity in annular flows is well-represented by trends in shear viscosity and stress relaxation behavior. However, the lack of viscosity synergy at high shear rates limits the applicability of the mixtures as a way to improve solids suspension capacity in annular flows.     

Resistance of rough plates in a rotating fluid

Generalizing Navier’s partial slip condition, the flow due to a rough or striated plate moving in a rotating fluid is studied. It is found that the motion of the plate, the fluid surface velocity, and the shear stress are in general not in the same direction. The solution is extended to the case of finite depth, or Couette slip flow in a rotating system. In this case an optimum depth for minimum drag is found. The solutions are also closed form exact solutions of the Navier–Stokes equations. The results are fundamental to flows with Coriolis effects.     

A low cost laser-speckle photographic technique for velocity measurement in slow flows

A low cost, low power laser-speckle photographic technique has been developed and is duscussed for the measurement of point velocities in slow laminar flows. The technique is particularly suitable for axisymmetric flows where the two velocity components can be easily measured. The accuracy of the technique is established by measurement of the velocity distribution for Poiseuille flow and from data obtained for acceleration of an inelastic Newtonian fluid through a four-to-one circular contraction. Preliminary results are also presented in the contracting flow field for a non shear-thinning highly elastic fluid. These data are particularly significant for verification of finite element numerical solutions currently being developed for viscoelastic fluids in circular entry flows.     

Non-modal instabilities of two-dimensional disturbances in plane Couette flow of a power-law fluid

Instabilities of fluid flows have traditionally been investigated by normal mode analysis, i.e. by linearizing the equations of flow and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In this paper we study the instabilities of two-dimensional Couette flow of a polymeric fluid in the framework of non-modal stability theory rather than normal mode analysis. A power-law model is used to describe the polymeric liquid. We focus on the response to external excitations and initial conditions by examining the pseudospectra structures and the transient energy growths. For both Newtonian and non-Newtonian flows, the results show that there can be a rather large transient growth even though the linear operator of Couette flow has no unstable eigenvalue. The effects of non-Newtonian viscosity on the transient behaviors are examined in this study. The results show that the “shear-thinning/shear-thickening” effect increases/decreases the amplitude of responses to external excitations and initial conditions.     

INTERFACIAL BEHAVIOR OF PIPE JOINTS UNDER TORSION LOADS

The interfacial behavior of pipe joints is studied in this paper. Firstly, through nonlinear fracture mechanics, the analytical expressions of interfacial shear stress and the loaddisplacement relationship at loaded end of pipe joints under torsion loads are obtained. Thus the shear stress propagation and the debonding process of the whole interface for different bond lengths can be predicted. Secondly, through the analytical solutions, the influences of different bond lengths on the load-displacement curve and the ultimate load are studied. The stress transfer mechanism, the interface crack propagation and the ductility behavior of the joints can be explained.     

Determination of yield stress fluid behaviour from inclined plane test

The aim of this paper is to determine precisely under which conditions an inclined plane can be used as a rheometer, which could represent a practical and rapid technique for various types of industrial or natural viscoplastic coarse suspensions. We first examine its efficiency and relevancy for determining fluid yield stress in a straight way by measuring the deepest fluid layer able to stay on the inclined plane. We have made experiments with different materials (clay-water suspensions) whose yield stress ranged from 35 to 90 Pa, using 1 m long open rectangular channels with a slope ranging from 10 to 30° and a width ranging from 5 to 25 cm. Our procedure involved measuring the final fluid depth far from edges a long time after the end of the slow gravity-induced emptying of a dam placed upstream. The fluid yield stress was also estimated independently by fitting a Herschel-Bulkley model to simple shear rheometry data obtained within a relatively wide shear rate range. A good agreement between inclined rectangular channel tests and independent usual rheometrical tests is obtained even for aspect ratios (flow depth to channel width ratio) as large as 1 when one assumes that, when the fluid has stopped, the side and bottom wall shear stresses are equal to the fluid yield stress. These results prove the efficiency of the inclined plane test for determining yield stress when appropriate experimental precautions are taken for both tests. In addition we examine the possibility of determining the simple shear flow curve of a mud suspension from fluid depth, velocity and discharge measurements of different steady flows in a wide open channel (8 m long; 60 cm wide) equipped with a recirculating system. The results obtained from inclined plane tests are in good agreement with independent rheometrical data (with torsional geometries). However it is technically difficult to cover a wide shear rate range from the inclined plane technique since this requires a rather wide channel flow rate range.     

Interaction of shear flows of an ideal incompressible fluid in a channel

The problem of the decay of an arbitrary discontinuity for the equations describing plane-parallel shear flows of an ideal fluid in a narrow channel is considered. The class of particular solutions corresponding to fluid flows with piecewise constant vorticity is studied. In this class, the existence of self-similar solutions describing all possible unsteady wave configurations resulting from the nonlinear interaction of the specified shear flows is established. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 34–47, November–December, 2006.     

Closure relations effects on the prediction of the stratified two-phase flow stability via the two-fluid model

The possibility of predicting the exact long wave linear stability boundary via the two-fluid (TF) model for horizontal and inclined stratified two-phase flow is examined. The application of the TF model requires the introduction of empirical closure relations for the velocity profile shape factors and for the wave induced wall and interfacial shear stresses. The latter are recognized as the problematic closure laws. In order to explore the closure relations effects and to suggest the necessary modifications that can improve the stability predictions of the TF model, the results are compared with the exact long wave solution of the Orr–Sommerfeld equations for the two-plate geometry. It is demonstrated that with the shape factors corrections and the inclusion of wave induced stresses effects, the TF model is able to fully reproduce the exact long wave neutral stability curves. The wave induced shear stresses in phase with the wave slope, which give rise to the so called “sheltering force”, were found to have a remarkable destabilizing effect in many cases of horizontal and inclined flows. In such cases, the sheltering effects must be included in the TF model, otherwise the region of smooth stratified flow would be significantly over predicted. Based on the results of the exact analysis, a simple closure relation for the sheltering term in the TF model is provided.