Stability of the interface in two-layer couette flow of upper convected maxwell liquids

The linear stability of a two-layer Couette flow of upper convected Maxwell liquids is considered. The fluids have different densities, viscosities, and elasticities, with surface tension at the interface. At low speeds, the interfacial mode may become unstable, while other modes stay stable. The shortwave asymptotics of the interfacial mode is analyzed. It is found that an elasticity difference can stabilize or destabilize the flow even in the absence of a viscosity difference. As the viscosity difference increases, the range of elasticities for which there is shortwave stability widens. A linearly stable arrangement can be achieved by placing the less viscous fluid in a thin layer to stabilize longwaves and using elasticities to stabilize shortwaves. Such an arrangement can be stable even when the density stratification is adverse.     

Weakly nonlinear instability of the flow of two immiscible liquids with different viscosities in a pipe

The weakly nonlinear stability due to an axisymmetric disturbance of two immiscible liquids with different viscosities in a pipe is worked out. The most important conclusions drawn from this work are the independence of the Reynolds number in the stability criterion, the existence of stable waves with sawtooth profile when the more viscous fluid is located at the core and it occupies most of the pipe, and the destabilizing effect of the surface tension, effect that can be overtaken by the stabilizing effect due to the difference in viscosity of the two fluids, when the viscosity of the ring fluid is much smaller than the viscosity of the core fluid and the thickness of the ring flow is smaller than the thickness of the core flow.     

Precise viscosity measurements of Newtonian liquids withv < 1 mm2/s for the selection of suitable standards

The application of long capillary Ubbelohde viscometers in the determination of low viscosities is described. Corrections and measurement uncertainties are discussed. Viscosity and density data are given for 18 Newtonian liquids of commercially available purity, with viscosities below about 1 mm²/s. For several liquids the commonly used purity specifications were found to be sufficient to qualify these liquids as viscosity standards, which can be used for viscometer calibration.     

Two-dimensional study of drop deformation under simple shear for Oldroyd-B liquids

The effect of viscoelasticity on the deformation of a circular drop suspended in a second liquid in shear is investigated with direct numerical simulations. A numerical algorithm based on the volume-of-fluid method for interface tracking is implemented in two dimensions with the Oldroyd-B constitutive model for viscoelastic liquids. The code is verified against a normal mode analysis for the stability of two-layer flow in a channel; theoretical growth rates are reproduced for the interface height, velocity and stress components. Drop simulations are performed for drop and matrix liquids of different viscosities and elasticities. A new feature is found for the case of equal viscosity, when the matrix liquid is highly elastic and surface tension is low; hook-like structures form at the drop tips. This is due to the growth of first normal stress differences that occur slightly above the front tip and below the back tip as the matrix elasticity increases above a threshold value.     

Magnetohydrodynamics instability of interfacial waves between two immiscible incompressible cylindrical fluids

The problem of nonlinear instability of interfacial waves between two immiscible conducting cylindrical fluids of a weak Oldroyd 3-constant kind is studied. The system is assumed to be influenced by an axial magnetic field, where the effect of surface tension is taken into account. The analysis, based on the method of multiple scale in both space and time, includes the linear as well as the nonlinear effects. This scheme leads to imposing of two levels of the solvability conditions, which are used to construct like-nonlinear Schr6dinger equations （1-NLS） with complex coefficients. These equations generally describe the competition between nonlinearity and dispersion. The stability criteria are theoret- ically discussed and thereby stability diagrams are obtained for different sets of physical parameters. Proceeding to the nonlinear step of the problem, the results show the appearance of dual role of some physical parameters. Moreover, these effects depend on the wave kind, short or long, except for the ordinary viscosity parameter. The effect of the field on the system stability depends on the existence of viscosity and differs in the linear case of the problem from the nonlinear one. There is an obvious difference between the effect of the three Oldroyd constants on the system stability. New instability regions in the parameter space, which appear due to nonlinear effects, are shown.     

An efficient and accurate fully discrete finite element method for unsteady incompressible Oldroyd fluids with large time step

This paper proposes a second‐order accuracy in time fully discrete finite element method for the Oldroyd fluids of order one. This new approach is based on a finite element approximation for the space discretization, the Crank–Nicolson/Adams–Bashforth scheme for the time discretization and the trapezoid rule for the integral term discretization. It reduces the nonlinear equations to almost unconditionally stable and convergent systems of linear equations that can be solved efficiently and accurately. Here, the numerical simulations for L², H¹ error estimates of the velocity and L² error estimates of the pressure at different values of viscoelastic viscosities α, different values of relaxation time λ₁, different values of null viscosity coefficient μ₀ are shown. In addition, two benchmark problems of Oldroyd fluids with different solvent viscosity μ and different relaxation time λ₁ are simulated. All numerical results perfectly match with the theoretical analysis and show that the developed approach gives a high accuracy to simulate the Oldroyd fluids under a large time step. Furthermore, the difference and the connection between the Newton fluids and the viscoelastic Oldroyd fluids are displayed. Copyright © 2015 John Wiley & Sons, Ltd.     

A microtube viscometer with a thermostat

The viscometer presented in this paper is suitable for measuring the viscosity of liquids in micro-litre quantities. It consists of a micro-flow experimental system with a thermostat. Using the measurements of the flow rates and pressure drops of a liquid passing through a microtube, the liquids viscosity can be calculated from on Hagen-Poiseuille theory. After calibration, the viscometer was used to measure viscosities of deionized water and ethyl alcohol at temperatures ranging from 0 to 40 °C. For both test liquids, the relative deviation of the measured values from those quoted in the literature (obtained using other viscometers) was less than 2.6%. The relative uncertainty of the experimental system was reduced to ±1.8% using the relative measuring method. Due to the micro-scale of the test section, only a micro-litre quantity of liquid is needed for a test; this is a potential advantage for measurement of bio-liquid viscosities.     

On the viscosity of magnetic fluid with low and moderate solid fraction

The design of a pressurized capillary rheometer operating at prescribed temperature is described to measure the viscosity of magnetic fluids (MFs) containing Fe3O4 magnetic nanoparticles (MNPs). The equipment constant of the rheometer was obtained using liquids with predetermined viscosities. Experimentally measured viscosities were used to evaluate different equations for suspension viscosities. Deviation of measured suspension viscosities from the Einstein equation was found to be basically due to the influence of spatial distribution and aggregation of Fe3O4 MNPs. By taking account of the coating layer on MNPs and the aggregation of MNPs in MFs, a modified Einstein equation was proposed to fit the experimental data. Moreover, the influence of external magnetic field on viscosity was also taken into account. Viscosities thus predicted are in good agreement with experimental data. Temperature effect on suspension viscosity was shown experimentally to be due to the shear-thinning behavior of the MFs.     

Purely elastic interfacial instabilities in superposed flow of polymeric fluids

Purely elastic interfacial stability of superposed plane Poiseuille flow of polymeric liquids has been investigated utilizing both asymptotic and numerical techniques. It is shown that these instabilities are caused by an unfavorable jump in the first normal stress difference across the fluid interface. To determine the significance of these instabilities in finite experimental geometries, a comparison between the maximum growth rates of purely elastic instabilities with instabilities driven primarily by a viscosity or a combined viscosity and elasticity difference is made. Based on this comparison, it is shown that purely elastic interfacial instabilities can play a major role in superposed flow of polymeric liquids in finite experimental geometries.     

Stability of Rheometric Viscoelastic Fluid Flow with Respect to Shear Perturbations

Flow between two plates is considered for a fluid obeying the DeWitt rheological equation of state with the Jaumann derivative. It is found analytically that the steady-state Couette flow is stable or unstable with respect to plane shear perturbations when the Weissenberg numbers are less or greater than unity, respectively. The flow acceleration stage is studied analytically and numerically, a comparison with the case of an Oldroyd fluid is carried out, and the neutral stability curves are constructed. The fundamental role of perturbations of the type considered among the set of instability types which can act on the fluid in such a flow is noted.     

The stability of two stratified ”power-law“ liquids in couette flow

Consideration is given to the stability of the flow of two ”power-law“ liquids between two infinite parallel planes when one of the planes moves with constant velocity in its own plane. It is found that the ratios of the power-law parameters for each layer have a dramatic effect and can be chosen to destabilize the flow.     

On the viscosity of magnetic fluid with low and moderate solid fraction

The design of a pressurized capillary rheometer operating at prescribed temperature is described to measure the viscosity of magnetic fluids （MFs） containing Fe3O4 magnetic nanoparticles （MNPs）. The equipment constant of the rheometer was obtained using liquids with predetermined viscosities. Experimentally measured viscosities were used to evaluate different equations for suspension viscosities. Deviation of measured suspension viscosities from the Einstein equation was found to be basically due to the influence of spatial distribution and aggregation of Fe3O4 MNPs. By taking account of the coating layer on MNPs and the aggregation of MNPs in MFs, a modified Einstein equation was proposed to fit the experimental data. Moreover, the influence of external magnetic field on viscosity was also taken into account. Viscosities thus predicted are in good agreement with experimental data. Temperature effect on suspension viscosity was shown experimentally to be due to the shear-thinning behavior of the MFs.     

Squeeze film flow of ideal elastic liquids

An exact solution is presented for the squeeze film flow of an Oldroyd B. fluid. The solution demonstrates that the flow kinematics is similar to the Newtonian (or Maxwellian) one. Theoretical predictions for constant velocity squeezing are compared to experimental observation for well characterized non-shear thinning elastic fluids. It is shown both theoretically and experimentally that the effect of elasticity in a constant velocity squeeze film flow is to always reduce the load relative to the inelastic (Newtonian) prediction and that this load reduction falls between the upper and lower asymptote prediction by the exact solution for the Oldroyd B fluid. The upper load asymptote is given by the Stefan solution for the viscosity of the polymer solution and the lower asymptote is given by the Stefan solution for the viscosity of the solvent. Experimental observations agree with the theoretical prediction for the Oldroyd B fluid at low shear rates where it is shown that the steady and dynamic flow properties of the test fluids used in the experimental program are well represented by the Oldroyd B constitutive equation. With the exception of the work of Lee et al. [6] for constant load squeezing of a Maxwell fluid, this work represents one of the few cases where experimental observation of large effects due to elasticity are indeed predicted with a constitutive equation which actually describes the steady and dynamic shear properties of the fluids used in the experimental program.     

Electroviscous potential flow in nonlinear analysis of capillary instability

We study the nonlinear stability of electrohydrodynamic of a cylindrical interface separating two conducting fluids of circular cross section in the absence of gravity using electroviscous potential flow analysis. The analysis leads to an explicit nonlinear dispersion relation in which the effects of surface tension, viscosity and electricity on the normal stress are not neglected, but the effect of shear stresses is neglected. Formulas for the growth rates and neutral stability curve are given in general. In the nonlinear theory, it is shown that the evolution of the amplitude is governed by a Ginzburg–Landau equation. When the viscosities are neglected, the cubic nonlinear Schrödinger equation is obtained. Further, it is shown that, near the marginal state, a nonlinear diffusion equation is obtained in the presence of viscosities. The various stability criteria are discussed both analytically and numerically and stability diagrams are obtained. It is also shown that, the viscosity has effect on the nonlinear stability criterion of the system, contrary to previous belief.     

Viscosity effect on the flow patterns in T-type micromixers

Flow patterns and mixing of liquids with different viscosities in T-type micromixers are numerically investigated on the Reynolds number range from 1 to 250. The viscosity ratio of the mixing media varied from 1 to 2; its effect on the flow structure and the mixing is studied. The dependences of the mixing efficiency and the pressure difference in the channel on the viscosity ratio and the Reynolds number are obtained. It is shown that the viscosity ratio has a considerable effect on the flow structure before and after transition from the symmetric to the asymmetric flow pattern. The self-similar behavior of the asymmetric flow pattern is established.     

Linear rheology of viscoelastic emulsions with interfacial tension

Emulsions of incompressible viscoelastic materials are considered, in which the addition of an interfacial agent causes the interfacial tension to depend on shear deformation and variation of area. The average complex shear modulus of the medium accounts for the mechanical interactions between inclusions by a self consistent treatment similar to the Lorentz sphere method in electricity. The resulting expression of the average modulus includes as special cases the Kerner formula for incompressible elastic materials and the Oldroyd expression of the complex viscosity of emulsions of Newtonian liquids in time-dependent flow.     

Effects of a transverse flows with variable viscosity in micropolar fluids

A regular perturbation analysis is presented for three laminar natural convection flows in micropolar fluids in liquids with temperature dependent viscosity: a freely-rising plane plume, the flow above a horizontal line source on an adiabatic surface (a plane wall plume) and the flow adjacent to a vertical uniform flux surface. While these flows have well-known power-low similarity solutions when the fluid viscosity is taken to be constant, they are non-similar when the viscosity is considered to a function of temperature. A single similar flow, that adjacent to a vertical isothermal surface, is also analysed for comparison in order to estimate the extent of validity of perturbation analysis. The formulation used here provides a unified treatment of variable viscosity effects on those four flows. Computed first-order perturbation quantities are presented for all four flows. Numerical results for velocity, angular velocity and thermal functions has been shown graphically or tabulated for different values of micropolar parameters. Received on 20 October 1997     

液晶高分子各向异性粘弹性流体本构方程理论

将液晶高分子－各向异性流体的本构方程,建立在Oldroyd随体导数观点基础上。推广上随机Oldroyd B流体模型,提出共转OldroydB流体模型,同时将微观结构的影响通过宏观参数表示出来,使在宏观理论中包含微观结构的贡献,即引入取向物质函数,非线性各向异性黏度函数和各向异性松弛时间及推迟时间等,表征取向运动对黏度和松弛及推迟现象的影响,在此基础上开展了一类新的液晶高分子－Oldroyd型本构方程理论,由该类型本构方程得出的物质函数,液晶高分子流体的第一、第二法向应力差与实验结果一致,解释了液晶高分子溶液的第一、第二法向应力差的特殊流变学行为。     

Interface stabilization in two-layer channel flow by surface heating or cooling

The linear stability of plane Poiseuille flow of two immiscible Newtonian liquids in a differentially heated channel is considered. The equations of motion and energy are fully coupled via temperature-dependent fluid-viscosity coefficients. A long-wave asymptotic formulation of the stability problem is presented together with numerical results for disturbances of arbitrary wavelength. Two combinations of immiscible liquids are analyzed: silicone/water and oil/water (water at the bottom layer in both cases). It is shown that an imposed wall temperature difference can be stabilizing or destabilizing depending on the disturbance wavenumber and layer thickness ratio. Interfacial tension has a stabilizing effect on the interface. Stabilizing influence of interfacial tension is observed at intermediate and large wavenumbers. Most importantly, for certain ranges of the controlling dimensionless parameters, stable interfaces at all disturbance wavelengths can be attained.     

A new method to determine viscosity of liquids using vibration principles

A new method for determining viscosity of liquids is examined. The method employs the principles of vibration and measures the viscous damping due to the motion of a liquid placed in a cylindrical tube. The apparatus and the test liquid are treated as a dynamic system and the measured mechanical impedances are used to calculate energy dissipation due to the viscous damping. The newly designed apparatus is able to generate shear deformations in the liquid without using moving solid surfaces. A harmonic varying force with a frequency close to the resonance frequency of the system is applied through a piston and the resulting velocities of the oscillations generated in the system are measured. Liquids with higher viscosities result in lower velocities due to the higher damping. Analytical equations are provided to relate the viscous damping of the dynamic system to the viscosity of the liquids. The viscosities obtained from the proposed method are in good agreement with the ones obtained from standard rotational viscometry using a cone and plate geometry.