Nitrogen/non-Newtonian fluid two-phase upward flow in non-circular microchannels

This paper presents experimental investigations on nitrogen/non-Newtonian fluid two-phase flow in vertical noncircular microchannels, which have square or triangular cross-section with the hydraulic diameters being D_h = 2.5, 2.886 and 0.866 mm, respectively, by visualization method. Three non-Newtonian aqueous solutions with typical rheological properties, i.e., 0.4% carboxymethyl cellulose (CMC), 0.2% polyacrylamide (PAM) and 0.2% xanthan gum (XG) are chosen as the working fluids. The common flow patterns are identified as slug flow, churn flow and annular flow. The dispersed bubble flow is only found in the case with nitrogen/CMC solution two-phase flow in the largest channel. A new flow pattern of nitrogen/PAM solution two-phase flow, named chained bubble/slug flow, is observed in all the test channels. The flow regime maps are also developed and the results show that the rheological properties of the non-Newtonian fluid have remarkable influence on the flow pattern transitions. The geometrical factors of the microchannel such as the cross-section shape and hydraulic diameter of the channel can also affect the flow regime map. Finally, the results obtained in this work are compared with the available flow pattern transitions.     

Phase split of liquid–liquid two-phase flow at a horizontal T-junction

This paper provides data on the split of liquid/liquid two-phase flow at a horizontal T-junction. Phase maldistribution was measured for kerosene–water flow at the T-junction with equal pipe diameters of 67.4 mm. Data were taken with both stratified flow with a mixture at interface and dispersed flows approaching the junction. The degree of phase maldistribution was not very great but preferential emergence of either phase from the side-arm was observed depending on the flow rates of the two-phases. There are similarities with the limited split data from liquid/solid flows and the degree of separation is seen to depend on the dispersed/continuous phase density ratio. The data were compared to predictions from the correlation by Seeger et al. The Seeger equation gives but reasonable agreement.     

The split of vertical annular flow at a large diameter T junction

Measurements and observations have been made of the split of gas—liquid flow at a T junction where the main pipe was vertical and the side arm horizontal. All three pipes connecting to the junction were of 0.125 m dia. The gas and liquid flow rates were chosen to ensure annular flow in the inlet pipe. The resulting data have been compared with existing models.     

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.     

Measurements and observations of the split of annular flow at a vertical T-junction

Measurements and observations have been made when annular flow divides at a vertical T. This work has extended earlier experiments in covering the entire range of take off. From the observations and measurements, three ways in which the liquid can be diverted into the side arm have been identified. A modification of an earlier model has been produced which correctly allows for two of the three phenomena.     

Numerical study on the two-phase flow distribution in a T-junction

The main objective of this paper is to investigate the ability of a two-dimensional two-fluid computer code to predict the phase separation in a T-junction. A new semi-implicit numerical scheme is developed for solving the two-fluid model equations. Special attention is directed to the modelling of the constitutive for the interfacial friction term. Detailed distribution of void fraction, pressure and velocities are obtained for an air–water mixture in a vertical tee. Good agreement was obtained between the computer code results and the experimental data for the phase separation in the T-junction.     

THE SOLUTION FOR TRANSIENT TWO-PHASE FLOW BY SPLIT FLUX VECTOR METHOD

In this paper the transient two-phase flow equations and their eigenvalues are first intreduced. The flux vector is then split into subvectors which just contain a specially signed eigenvalue. Using one-sided spatial difJerence operators finite difference equations and their solutions are obtained. Finally comparison with experiment shows the predicted results produce good agreement with experimental data.     

Fluctuation characteristics of two-phase flow splitting at a vertical impacting T-junction

In order to investigate the fluctuation characteristics of two-phase flow splitting at a T-junction, particular attention was paid on Churn flow which had the strongest fluctuation comparing with bubble flow and annular flow. The main tube of the T-junction was vertical and the two branches were horizontal. All three pipes connecting to the junction were of 15 mm inner diameter. A statistical analysis based on Root Mean Square (RMS) was applied to temporal differential pressure signals and gas flow rate signals. The Power Spectral Density (PSD) was also employed to reveal their peculiar features in frequency domain as well. The effects of the extraction flow ratio and the gas and liquid superficial velocity upstream on fluctuation characteristics of gas-liquid two-phase flow splitting at the T-junction were investigated in detail. It is found that there is a wide fluctuation in both differential pressure and gas flow rate downstream at every extraction ratio (W₃/W₁) and the fluctuation intensity increases as W₃/W₁ increasing. It is also made clear that increasing either water superficial velocity or gas superficial velocity in inlet causes fluctuation to become more intensive.     

On the case when steady converging/diverging flow of a non-Newtonian fluid in a round cone permits an exact solution

This communication considers the steady converging/diverging flow of a non-Newtonian viscous power-law fluid in a round cone. The motion is driven by a sink/source of mass at the origin. It is shown that the problem permits exact similarity solution for a particular value (n=4/3) of the fluid index. In this case a complete set of governing equations can be reduced to an ordinary differential equation, which is solved numerically for different values of the main non-dimensional parameters (the cone angle and the dimensionless sink/source intensity).     

Two-phase flow split at T junctions: Effect of side arm orientation and downstream geometry

The effect of flow pattern and geometry on the phase split of gas/liquid flows at T junctions has been examined for a horizontal main tube and horizontal and vertically upwards side arms. Important phenomena which control this split in annular and wavy stratified flow have been identified. The capability of current models to predict the split are discussed. In particular, the effect of geometry in the downstream leg of the main pipe was studied. The configurations studied had no effect in annular flow but influenced the amount of liquid taken off at high take off when stratified flow approached the junction.     

Heat transfer to non-Newtonian flows over a cylinder in cross flow

Over a range of 102<Re^*<5800, 6.5<Pr^*<79, and 0.6<n<1, circumferential wall temperatures for water and aqueous polymer (purely viscous) solution flows over a smooth cylinder were measured experimentally. The cylinder was heated by passing direct electric current through it. Aqueous solutions of Carbopol 934 and EZ1 were used as power-law non-Newtonian fluids. The peripherally averaged heat transfer coefficient for purely viscous non-Newtonian fluids, at any fixed flow rate, decreases with increasing polymer concentration. A new correlation is proposed for predicting the peripherally averaged Nusselt number for power-law fluid flows over a heated cylinder in cross flow.     

Experimental study on a two-phase critical flow with a non-condensable gas at high pressure conditions

An experimental study was performed on a two-phase critical flow with a non-condensable gas at high pressure conditions. Experimental data for the critical flow rates were generated by using sharp-edged stainless steel pipes with an inner diameter of 10.9 mm, a thickness of 3.2 mm, and a length of 1000 mm. The test conditions were varied by using the stagnation pressures of 4.0, 7.0, and 10.0 MPa, water subcoolings of 0.0, 20.0, and 50.0 °C, and nitrogen gas flow rates of 0.0–0.22 kg/s. The experimental results show that the critical mass flux decreases rapidly with an increase of the volumetric non-condensable gas fraction. Also the critical mass flux increases with an increase of the stagnation pressure and a decrease of the stagnation temperature. An empirical correlation of the non-dimensional critical mass flux, which is expressed as an exponential function of the non-condensable gas fraction of the volumetric flow, is obtained from the experimental data.     

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.     

Comparison of two numerical methods for the solution of non-Newtonian flow in ducts

Two numerical methods, the Galerkin finite element method (FEM) and the boundary-fitted co-ordinate transformation method (BFCTM), have been applied to solve inelastic non-Newtonian fluid flow in ducts of irregular cross-section. Three representative fluid models, namely the power-law, the Ellis and the Bingham models, have been analysed. The application of the FEM is straightforward, while for the BFCTM the accurate estimation of viscosity on the duct boundary and the proper mesh adjustment appear to be critical for generating convergent solutions. A detailed comparison of the two numerical methods in terms of volumetric flow rate, axial velocity, shear rate, viscosity and CPU time is given. Both methods can generate accurate solutions of velocity over a wide range of variables, but the FEM requires much less computing time to reach the same level of accuracy. Only the BFCTM can be used to approximate shear rate and viscosity with reasonable accuracy.     

Phase and velocity distributions in vertically upward high-viscosity two-phase flow

The two-phase pressure drop in vertical industrial pipes is mainly determined by gravitation and acceleration of the fluid, which means that the void fraction is key an important parameter in any model to predict pressure drops. Typically, these models are applied in industry to size pumps and, e.g., emergency relief systems. There is a shortage of void fraction data in the literature for liquids with a dynamic viscosity above 1000 mPa s. Adiabatic experiments have been performed of mixtures of nitrogen and solutions of polyvinylpyrrolidone (Luviskol^®) in water with dynamic viscosities in the range 900–7000 mPa s. Inner tube diameter was 54.5 mm. Mass flux and quality were varied in a wide range: 8–3500 kg/m²/s and 0–82%, respectively. The corresponding superficial velocities were 0.005–3.4 m/s for the liquid and 0–30 m/s for the nitrogen. For comparison, reference measurements were taken of mixtures of nitrogen with water (1 mPa s). Care has been taken to measure only well-developed flows.     

Flow regime transition criteria for two-phase flow at reduced gravity conditions

Flow regime transition criteria are of practical importance for two-phase flow analyses at reduced gravity conditions. Here, flow regime transition criteria which take the frictional pressure loss effect into account were studied in detail. Criteria at reduced gravity conditions were developed by extending an existing model from normal gravity to reduced gravity conditions. A comparison of the newly developed flow regime transition criteria model with various experimental datasets taken at microgravity conditions showed satisfactory agreement. Sample computations of the model were performed at various gravity conditions, such as 0.196, 1.62, 3.71 and 9.81 m/s² corresponding to micro-gravity and lunar, Martian and Earth surface gravity, respectively. It was found that the effect of gravity on bubbly–slug and slug–annular (churn) transitions in a two-phase flow system was more pronounced at low liquid flow conditions, whereas the gravity effect could be ignored at high mixture volumetric flux conditions. While for the annular flow transitions due to flow reversal and onset of droplet entrainment, higher superficial gas velocity was obtained at higher gravity level.     

Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

Transient MHD stagnation flow of a non-Newtonian fluid due to impulsive motion from rest

The transient boundary layer flow and heat transfer of a viscous incompressible electrically conducting non-Newtonian power-law fluid in a stagnation region of a two-dimensional body in the presence of an applied magnetic field have been studied when the motion is induced impulsively from rest. The non-linear partial differential equations governing the flow and heat transfer have been solved by the homotopy analysis method and by an implicit finite-difference scheme. For some cases, analytical or approximate solutions have also been obtained. The special interest are the effects of the power-law index, magnetic parameter and the generalized Prandtl number on the surface shear stress and heat transfer rate. In all cases, there is a smooth transition from the transient state to steady state. The shear stress and heat transfer rate at the surface are found to be significantly influenced by the power-law index N except for large time and they show opposite behaviour for steady and unsteady flows. The magnetic field strongly affects the surface shear stress, but its effect on the surface heat transfer rate is comparatively weak except for large time. On the other hand, the generalized Prandtl number exerts strong influence on the surface heat transfer. The skin friction coefficient and the Nusselt number decrease rapidly in a small interval 0<t^*<1 and reach the steady-state values for t^*≥4.