Deformation and breakup of viscoelastic drops in planar extensional flows

A computer-controlled four-roll mill was used to investigate the deformation and break-up of polymeric drops in the well-characterized flow of an immiscible Newtonian fluid. Aqueous polymer solutions ranging in concentration from 160 ppm to 3% by weight were examined. For zero-shear-rate viscosity ratios greater than order 1, the deformation of the drops closely followed that of Newtonian fluids, irrespective of the droplet material. However, drops with viscosity ratios less than order 1 had significantly smaller critical deformations and the critical capillary number was found to be substantially smaller. Two modes of drop break-up were discovered that differed substantially from that observed for Newtonian drops in the inclusion of cusped ends and tip streaming.     

Transient deformation of an inviscid inclusion in a viscoelastic extensional flow

The transient deformation of a bubble in a viscoelastic extentional flow is analyzed by means of a finite element algorithm for viscoelastic moving boundary problems. Using the Oldroyd-B constitutive model, we find that bubbles in a viscoelastic fluid deform to the same steady-state configurations as bubbles in a Newtonian fluid at equal values of the far-field extensional stresses (corresponding to different stretch rates). Vapor bubbles in a developed extensional flow collapse more readily in the viscoelastic liquid than bubbles in Newtonian fluids because of the large compressive stresses associated with the viscoelastic liquid.     

Effect of interfacial modifier on single drop deformation and breakup in step increasing shear flow

This work deals with in situ visualisation of deformation and breakup of a copolymer modified single Newtonian drop immersed in a Newtonian homogenous matrix. The experiments were carried out on a model system made of poly-isobutylene as the suspending fluid and two poly-dimethylsiloxanes with different molecular weights as the drop phase with viscosity ratios 0.036 and 1.13, below and above but close to unity. Three weight concentrations 0.5%, 2% and 10% of the block copolymer laying below, close to and above the critical concentration of the total drop surface coverage were examined. Single drop deformation experiments were carried out in a home-designed Couette quartz cell connected to a home-modified Paar Physica Rheometer. The variation in the length-to-diameter ratio (L/D) versus shear rate and capillary number was measured both in steady and in transient regimes till breakup. The results indicated a weaker resistance of copolymer modified drops against hydrodynamic stresses at both viscosity ratios as compared to the clean drop. However, the drop deformation was found to be complex and depends on the copolymer concentration and the viscosity ratio.     

A sharp slender cone in an incompressible flow

The range of applicability of the asymptotic solution of the problem of incompressible flow past a slender cone is studied. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 197–200, July–August, 1998.     

Heat transfer and pressure drop for purely viscous non-Newtonian fluids in turbulent flow through rectangular passages

Experimental measurements of friction factor and heat transfer for the turbulent flow of purely viscous non-Newtonian fluids in a 21 rectangular channel are compared with results previously reported for the circular tube geometry. Comparisons are also made with available analytical and empirical predictions.It is found that the rectangular duct fully established friction factor measurements are within ± 5% of the Dodge-Metzner prediction if the Kozicki generalized Reynolds number is used. A modified form of the simpler explicit equation proposed by Yoo, [i.e.f=0.079n ^0.675(Re ^*)^–0.25], is found to yield predictions for both the rectangular duct and the circular tube geometries with approximately the same accuracy as the Dodge-Metzner equation.Fully developed Stanton numbers for the rectangular duct are in good agreement with the circular tube data over a range ofn from 0.37 to 0.88 for a given Prandtl number,Pr _a , when compared at a fixed value of the Reynolds number based on the apparent viscosity evaluated at the wall shear stress. In general, the experimental data are within ± 20% of Yoo's equation,St=0.0152Re _a ^–0.155 Pr _a ^–2/3 . A new equation is proposed to bring the prediction for circular pipes as well as rectangular channels into better agreement with generally accepted Newtonian heat transfer results.

---

Wärmeübergang und Druckverlust für viskose nicht-Newtonsche Fluide in turbulenter Strömung durch rechteckige Kanäle

Zusammenfassung Es werden Messungen des Reibungsfaktors und des Wärmeübergangs bei turbulenter Strömung viskoser nicht-Newtonscher Fluide in einem rechteckigen Kanal mit dem Seitenverhältnis 21 verglichen mit früheren Ergebnissen, die an runden Rohren gewonnen wurden. Weiterhin werden Vergleiche mit aus der Literatur verfügbaren analytischen und empirischen Beziehungen gemacht.Es zeigte sich, daß die Messungen des Reibungsfaktors im rechteckigen Kanal bei vollausgebildeter Strömung auf ± 5% mit der Vorhersage von Dodge-Metzner übereinstimmen, wenn die von Kozicki verallgemeinerte Reynolds-Zahl verwendet wird. Eine modifizierte Form der einfachen von Yoo vorgeschlagenen einfachen Gleichung in explizierter Form (f=0,079n ^0,675(Re ^*)^–0,25) bewies, daß sie sowohl für den rechteckigen Kanal als auch das runde Rohr die Werte mit fast der gleichen Genauigkeit wie die Methode von Dodge-Metzner vorhersagen kann.Die Stanton-Zahlen für den rechteckigen Kanal bei vollausgebildeter Strömung sind in guter Übereinstimmung mit den Werten für das runde Rohr in einem Bereich vonn= 0,37 – 0,88 für eine gegebene Prandtl-Zahl, wenn man den Vergleich bei einem vorgegebenen Wert der Reynolds-Zahl anstellt, die auf die scheinbare Viskosität — abgeleitet aus der Wandschubspannungbezogen ist. Generell läßt sich sagen, daß die Werte auf ± 20% mit der Gleichung von Yoo (St=0,0152Re _a ^–0,155 )Pr _a ^–2/3 ) übereinstimmen. Es wird eine neue Gleichung vorgeschlagen, welche sowohl die Werte für runde Rohre als auch die für rechteckige Kanäle in bessere Übereinstimmung bringt mit den in der Literatur üblichen Ergebnissen für den Wärmeübergang an Newtonsche Fluide.

---

Nomenclature a constant in Eq. (8) - A area of cross-section of channel [m²] - b constant in Eq. (8) - c _p specific heat of test fluid [J kg^–1 K^–1] - d capillary tube diameter [m] - D _h hydraulic diameter, 4A/P[m] - f Fanning friction factor, _w/(g9 V²/2) - h axially local (spanwise averaged) heat transfer coefficient,q _w /(T_wi-T_b) [Wm^–2 K^–1] - k _f thermal conductivity of test fluid [Wm^–1K^–1] - K consistency index of power law fluid - n power law index - Nu fully established, local (spanwise averaged) Nusselt numberh D _h /k _f - P perimeter of channel [m] - Pr _a Prandtl number based on apparent viscosjity, c _p /k _f - Pr ^* defined as (Re _a Pr _a )/Re ^* - q _w wall heat flux [Wm^–2] - Re _a Reynolds number based on apparent viscosity, VD _h/ - Re Metzner's generalized Reynolds number in Eq. (2) - Re ^* Reynolds number defined in Eq. (8) - St Stanton number,h/( V c_p) - T _b local bulk temperature of the fluid [K] - T _wi local inside wall temperature [K] - T _wo local outside wall temperature [K] - V bulk flow velocity [m s^–1] - x distance from the inlet of channel along flow direction [m] Greek symbols shear rate [s^–1] - apparent viscosity [Pa s] - density of test fluid [kg m^–3] - shear stress [Pa] - _w shear stress at the wall [Pa] Dedicated to Prof. Dr.-Ing. U. Grigull's 75th birthday     

A boundary‐only approach to the deformation of a shear‐thinning drop in extensional Newtonian flow

The influence of shear thinning on drop deformation is examined through a numerical simulation. A two‐dimensional formulation within the scope of the boundary element method (BEM) is proposed for a drop driven by the ambient flow inside a channel of a general shape, with emphasis on a convergent–divergent channel. The drop is assumed to be shear thinning, obeying the Carreau–Bird model and the suspending fluid is Newtonian. The viscosity of the drop at any time is estimated on the basis of a rate‐of‐strain averaged over the region occupied by the drop. The viscosity thus changes from one time step to the next, and it is strongly influenced by drop deformation. It is found that small drops, flowing on the axis, elongate in the convergent part of the channel, then regain their spherical form in the divergent part; thus confirming experimental observations. Newtonian drops placed off‐axis are found to rotate during the flow with the period related to the initial extension, i.e. to the drop aspect ratio. This rotation is strongly prohibited by shear thinning. The formulation is validated by monitoring the local change of viscosity along the interface between the drop and the suspending fluid. It is found that the viscosity averaged over the drop compares, generally to within a few per cent, with the exact viscosity along the interface.     

Nonisothermal flow of non-Newtonian fluids in a channel

A study is made of the nonisothermal motion of a rheologically complex fluid in a plane-parallel channel in the case of boundary conditions of the third kind on the outer surfaces of the channel walls and with allowance for the dissipation of mechanical energy and temperature dependence of the coefficients. A qualitative investigation of the problem is made for arbitrary temperature dependence of the yield. Special cases are considered: flow of linear viscoplastic medium and a power-law medium with abrupt change with the temperature of the yield point and consistency. It is shown that under certain conditions several different flow regimes can be realized simultaneously in the channel and the change in the flow rate of the medium in a channel with varying temperature of the surrounding medium can exhibit hysteresis.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 3–10, May–June, 1980.     

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

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

Transition on a slender cone in hypervelocity flow

 An experimental exploratory study of laminar, transitional and turbulent boundary layers on a slender, sharp cone at zero incidence in hypervelocity flow was conducted in the shock tunnel T5, whose high operating pressure permits achieving the necessary high Reynolds numbers at high total enthalpy. Surface heat transfer rate and resonantly enhanced flow visualization were used to determine transition location and flow structure. Results suggest that the transition mechanism at the conditions tested is via the Tollmien–Schlichting instability, but since no detailed information about the noise spectrum in the facility is available, a definite conclusion is not possible. A strong, gas-dependent increase of the transition Reynolds number evaluated at the reference temperature on specific total enthalpy was observed. This effect increases monotonically with decreasing the dissociation energy of the gas. Received: 1 February 1996/Accepted: 28 June 1996     

On the breakup of slender liquid bridges: Experiments and a 1-D numerical analysis

Slender axisymmetric dielectric liquid bridges are made stable by the action of an axial electric field. In this paper, the subsequent dynamics of a slender liquid bridge after turning off the electric field is considered. The evolution in time of the bridge profile is investigated both theoretically and experimentally. A one-dimensional model is used to simulate the dynamic response of the system. Experiments are performed applying an axial electric field to a liquid bridge of 1 mm of diameter, and turning-off the electric field. The evolution of the liquid bridge is recorded using a video camera, and the digitized images are analysed. Good agreement between computations and experiments is found.     

Steady flow of a power-law non-Newtonian fluid across an unconfined square cylinder

A two-dimensional flow of a non-Newtonian power-law fluid directed normally to a horizontal cylinder with a square cross section is considered in the present paper. The problem is investigated numerically with a finite volume method by using the commercial code Ansys Fluent with a very large computational domain so that the flow could be considered unbounded. The investigation covers the power-law index from 0.1 to 2.0 and the Reynolds number range from 0.001 to 45.000. It is found that the drag coefficient for low Reynolds numbers and low power-law index (n ≤ 0.5) obeys the relationship C_D = A/Re. An equation for the quantity A as a function of the power-law index is derived. The drag coefficient becomes almost independent of the power-law index at high Reynolds numbers and the wake length changes nonlinearly with the Reynolds number and power-law index.     

Birefringence measurements in extensional flow of polyisobutylene around an expanding bubble

Birefringence in liquid polymers offers the possibility of obtaining information about stress in complex flows. In this work, this is done for extensional flows of polyisobutylene in a “breathing bubble” rheometer. In this type of rheometer, a bubble consisting of an incompressible, low-viscosity fluid (usually water) is injected into the sample with a nozzle. Expanding or collapsing the bubble by adding or removing water induces biaxial or uniaxial extension in the surrounding sample. The pressure difference between the bubble and the surroundings can be measured and compared to the predictions of constitutive equations. This measurement only gives one integral value for a complex flow history. In this paper, the birefringence around the bubble is measured in order to learn more about the flow. This is done by comparing pressure and birefringence results to those of standard constitutive equations for a polyisobutylene sample. A good agreement between the pressure and optical measurements and the theory is found with a single value of the stress-optical constant. Received: 25 June 1997 Accepted: 12 November 1997     

Spherical flow of a diatomic gas from an evaporating drop

The spherical expansion of gas from an evaporating drop is investigated on the basis of the numerical solution to a model kinetic equation for a gas with rotational degrees of freedom. Examples considered are the stationary evaporation of a drop with given temperature into the vacuum and evaporation of a drop into a gas-filled space under the condition of an energy balance on the drop surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 184–187, July–August, 1980.     

Free-surface non-Newtonian fluid flow in a round pipe

Free-surface pseudoplastic and viscoplastic fluid flows in a round pipe were studied for the case where the direction of motion coincides with the direction of gravity. Numerical modeling was performed using a technique based on a combination of the SIMPLE algorithm and the method of invariants. Three characteristic filling regimes were found to exist: a complete filling regime, a regime characterized by air-cavity formation on the solid wall, and a jet regime. Critical parameter values separating the regions of existence of these regimes were calculated. The evolution of quasisolid cores was studied for flow of a fluid with an yield point.