Peristaltic transport of two immiscible viscous fluids in a circular tube

Peristaltic motion of two immiscible viscous incompressible fluids in a circular tube is studied in pumping and copumping ranges under long-wavelength and low-Reynolds-number assumptions. The effect of the peripheral-layer viscosity on the time-averaged flux and the mechanical efficiency is studied. The formation and growth of the trapping zone in the core and the peripheral layer are explained. It is observed that the bolus volume in the peripheral layer increases with an increase in the viscosity ratio. The limits of the time-averaged flux (Q) over bar for trapping in the core are obtained. The trapping observed in the peripheral layer decreases in size with an increase in (Q) over bar but never disappears. The development of the complete trapping of the core fluid by the peripheral-layer fluid with an increase in the time-averaged flux is demonstrated. The effect of peripheral-layer viscosity on the reflux layer is investigated. It is also observed that the reflux occurs in the entire pumping range for all viscosity ratios and it is absent in the entire range of copumping.

Instabilities of soft elastic microtubes filled with viscous fluids: Pearls, wrinkles, and sausage strings

A linear stability analysis is presented to study the self-organized instabilities of a highly compliant elastic cylindrical shell filled with a viscous liquid and submerged in another viscous medium. The prototype closely mimics many components of micro-or nanofluidic devices and biological processes such as the budding of a string of pearls inside cells and sausage-string formation of blood vessels. The cylindrical shell is considered to be a soft linear elastic solid with small storage modulus. When the destabilizing capillary force derived from the cross-sectional curvature overcomes the stabilizing elastic and in-plane capillary forces, the microtube can spontaneously self-organize into one of several possible configurations; namely, pearling, in which the viscous fluid in the core of the elastic shell breaks up into droplets; sausage strings, in which the outer interface of the mircrotube deforms more than the inner interface; and wrinkles, in which both interfaces of the thin-walled mircrotube deform in phase with small amplitudes. This study identifies the conditions for the existence of these modes and demonstrates that the ratios of the interfacial tensions at the interfaces, the viscosities, and the thickness of the microtube play crucial roles in the mode selection and the relative amplitudes of deformations at the two interfaces. The analysis also shows asymptotically that an elastic fiber submerged in a viscous liquid is unstable for Y = gamma/(G(e)R) > 6 and an elastic microchannel filled with a viscous liquid should rupture to form spherical cavities (pearling) for Y > 2, where gamma, G(e), and R are the surface tension, elastic shear modulus, and radius, respectively, of the fiber or microchannel.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. The fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent.

Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Theoretical Analysis and Experimental Validation of the Scholte Wave Propagation in Immersed Plates for the Characterization of Viscous Fluids

This work reports the study of an attractive interfacial wave for application in ultrasonic NDE techniques for inspection and fluid characterization. This wave, called quasi-Scholte mode, is a kind of flexural wave in a plate in contact with a fluid which presents a good sensitivity to the fluid properties. In order to explore this feature, the phase velocity curve of quasi-Scholte mode is experimentally measured in a plate in contact with a viscous fluid, showing a good agreement with theory.

Damping of Vertically Excited Surface Wave in Weakly Viscous Fluid

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Variations of local heat transfer coefficient in piped flow of viscous liquids

Measurements were carried out to determine local coefficients of heat transfer in short lengths of horizontal pipe, and in the region of an discontinuity in pipe diameter. Laminar, transitional and turbulent flow regimes were investigated, and mixtures of propylene glycol and water were used in the experiments to give a range of viscous fluids. Theoretical and empirical analyses were implemented to find how the fundamental mechanism of forced convection was modified by the secondary effects of free convection, temperature dependent viscosity, and viscous dissipation. From experiments with the short tube it was possible to determine simple empirical relationships describing the axial distribution of the local 1usselt number and its dependence on the Reynolds and Prandtl numbers. Small corrections were made to account for the secondary effects mentioned above. Two different entrance configurations were investigated to demonstrate how conditions upstream could influence the heat transfer coefficients measured downstream In experiments with a sudden contraction in pipe diameter the distribution of local 1u3se1t number depended on the Prandtl number of the fluid in a complicated way. Graphical data is presented describing this dependence for a range of fluids indicating how the local Nusselt number varied with the diameter-ratio. Ratios up to 3.34:1 were considered. With a sudden divergence in pipe diameter, it was possible to derive the axial distribution of the local Nusse1t number for a range of Reynolds and Prandtl numbers in a similar way to the convergence experiments. Difficulty was encountered in explaining some of the measurements obtained at low Reynolds numbers, and flow visualization techniques wore used to determine the complex flow patterns which could lead to the anomalous results mentioned. Tests were carried out with divergences up to 1:3.34 to find the way in which the local Nusselt number varied with the diameter ratio, and a few experiments were carried out with very large ratios up .to 14.4. A limited amount of theoretical analysis of the 'divergence' system was carried out to substantiate certain explanations of the heat transfer mechanisms postulated.

A variational conjugate gradient method for determining the fluid velocity of a slow viscous flow

The problem considered is that of determining the fluid velocity for linear hydrostatics Stokes flow of slow viscous fluids from measured velocity and fluid stress force on a part of the boundary of a bounded domain. A variational conjugate gradient iterative procedure is proposed based on solving a series of mixed well-posed boundary value problems for the Stokes operator and its adjoint. In order to stabilize the Cauchy problem, the iterations are ceased according to an optimal order discrepancy principle stopping criterion. Numerical results obtained using the boundary element method confirm that the procedure produces a convergent and stable numerical solution.

Breakage of viscoelastic drops in turbulent stirred dispersions

The existing models of drop breakage in stirred turbulent dispersions are applicable only to purely viscous dispersed phases. In their present form, they are found to underpredict the diameters of the largest stable drops formed when a viscoelastic fluid is dispersed into a Newtonian liquid. In purely viscous fluids, the turbulent stresses are opposed both by the stresses due to interfacial tension and the viscous stresses generated as the drop deforms. In viscoelastic fluids, drop deformation produces additional retractive elastic stresses which also oppose turbulent stresses. As the deformation rates are large, the retractive stresses can be large in magnitude. Assuming that these additional stresses decay with time, a model of viscoelastic drop breakage in turbulent stirred dispersions has been developed. The new model quantitatively predicts the dmax of viscoelastic fluids. The model, however, does not predict the observation that when the time constant of the fluid becomes large (λ > 0.5 s), the fluid can not be dispersed into droplets up to agitator speeds of about 10 rps in our equipment.

Polyethylene stibinite phosphate esters: Novel flame-retardant plasticizers for PVC

A new series of multielement flame-retardant plasticizers containing polyethylene stibinite phosphate esters have been prepared by bulk polymerization from ethylene glycol with various antimony (III) aryloxydichlorides and arylphosphorodichloridates possessing various combinations of substituent [Cl,Br,NO2]. All the polymers are pink-coloured viscous fluids. They were characterized by inherent viscosity, density, IR, H-1, C-13 and P-31 NMR spectroscopy. The thermal behaviour of the polymers was compared by thermogravimetric analysis and correlated with their structures. The flammability studies were carried out by the limiting oxygen index test. The polymers containing P, Sb, N and Pr elements in their backbone show superior thermal-and flame-retardant characteristics than the other polymers. A comparative study was carried out with one of the synthesized polymers as a polymeric flame-retardant additive to plasticized PVC. The results showed improved LOI and mechanical properties to that of the conventional flame-retardant additive composition. (C) 1997 Elsevier Science Ltd.

Dynamics and fracture of ligaments from a droplet on a vibrating surface

A droplet residing on a vibrating surface and in the pressure antinode of an asymmetric standing wave can spread radially outward and atomize. In this work, proper orthogonal decomposition through high speed imaging is shown to predict the likelihood of atomization for various viscous fluids based on prior information in the droplet spreading phase. Capillary instabilities are seen to affect ligament rupture. Viscous dissipation plays an important role in determining the wavelength of the most unstable mode during the inception phase of the ligaments. However, the highest ligament capillary number achieved was less than 1, and the influence of viscosity in the ligament growth and breakup phases is quite minimal. It is inferred from the data that the growth of a typical ligament is governed by a balance between the inertial force obtained from the inception phase and capillary forces. By including the effect of acoustic pressure field around the droplet, the dynamics of the ligament growth phase is revealed and the ligament growth profiles for different fluids are shown to collapse on a straight line using a new characteristic time scale.

Effect of viscosity and surface tension on breakup and coalescence of bicomponent sprays

This work deals with an experimental study of the breakup characteristics of liquids with different surface tension and viscosity from a hollow cone hydraulic injector nozzle induced by pressure-swirl. The experiments were conducted at Reynolds numbers Re-p=9500-23,000. The surface tension and viscosity of the surrogate fuels were altered from 72 to 30 mN/m and 1.1 to 1.6 mN s/m(2), respectively. High speed photography and Phase Doppler Particle Anemometry were utilized to study the atomization process. Velocity and drop size measurements of the spray using PDPA in both axial and radial directions indicate a dependency on surface tension. However, these effects are dominant only at low Reynolds numbers and are negligible at high Reynolds number. Downstream of the nozzle, coalescence of droplets due to collision was also found to be significant and the diameters were compared for different liquids. For viscous fluids up to 1.6 cP, the independent effects of viscosity and injection pressure are studied. In general, the spray cone angle increases with increase in pressure. At high pressures, an increase in viscosity leads to higher drop sizes following primary and secondary breakup compared to water. This study will extend our understanding of surrogate fuel film breakup and highlight the importance of long and short wavelength instabilities. (C) 2013 Elsevier Ltd. All rights reserved

Experimental studies on the flow through soft tubes and channels

Experiments conducted in channels/tubes with height/diameter less than 1 mm with soft walls made of polymer gels show that the transition Reynolds number could be significantly lower than the corresponding value of 1200 for a rigid channel or 2100 for a rigid tube. Experiments conducted with very viscous fluids show that there could be an instability even at zero Reynolds number provided the surface is sufficiently soft. Linear stability studies show that the transition Reynolds number is linearly proportional to the wall shear modulus in the low Reynolds number limit, and it increases as the 1/2 and 3/4 power of the shear modulus for the `inviscid' and `wall mode' instabilities at high Reynolds number. While the inviscid instability is similar to that in the flow in a rigid channel, the mechanisms of the viscous and wall mode instabilities are qualitatively different. These involve the transfer of energy from the mean flow to the fluctuations due to the shear work done at the interface. The experimental results for the viscous instability mechanism are in quantitative agreement with theoretical predictions. At high Reynolds number, the instability mechanism has characteristics similar to the wall mode instability. The experimental transition Reynolds number is smaller, by a factor of about 10, than the theoretical prediction for the parabolic flow through rigid tubes and channels. However, if the modification in the tube shape due to the pressure gradient, and the consequent modification in the velocity profile and pressure gradient, are incorporated, there is quantitative agreement between theoretical predictions and experimental results. The transition has important practical consequences, since there is a significant enhancement of mixing after transition.

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.