Vorticity in plane parallel, axially symmetrical or spherical flows of viscous fluid

For viscous (barotropic or incompressible) fluids it is shown that, if the vorticity and the viscous force are orthogonal, vortex lines are convected by a vector field which fits with the velocity field when viscosity vanishes (extension of Helmholtz theorem); it is also found that energy remains constant along the field lines of this vector field (extension of Bernoulli theorem).If, moreover, vorticity and velocity are orthogonal too, the magnitude of the vorticity then behaves as the density of a fluid which flows along streamsheets according to this very same vector field. These properties are mainly encountered for plane parallel flows, axially symmetrical flows, spherical flows, but also for some other miscellaneous flow geometries such as unidirectional or radial flows. The set of the former three flows can even be characterized by these properties; that enhances this set of important flow geometries, avails a general view on vorticity behavior, and explains the great simplicity of vorticity equations in these cases. Numerous examples and comments are given for illustrating.     

Electrohydrodynamic instability of two superposed Walters B´ viscoelastic fluids in relative motion through porous medium

Summary  The electrohydrodynamic Kelvin–Helmholtz instability of the interface between two uniform superposed viscoelastic (B′ model) dielectric fluids streaming through a porous medium is investigated. The considered system is influenced by applied electric fields acting normally to the interface between the two media, at which there are no surface charges present. In the absence of surface tension, perturbations transverse to the direction of streaming are found to be unaffected by either streaming and applied electric fields for the potentially unstable configuration, or streaming only for the potentially stable configuration, as long as perturbations in the direction of streaming are ignored. For perturbations in all other directions, there exists instability for a certain wavenumber range. The instability of this system can be enhanced (increased) by normal electric fields. In the presence of surface tension, it is found also that the normal electric fields have destabilizing effects, and that the surface tension is able to suppress the Kelvin–Helmholtz instability for small wavelength perturbations, and the medium porosity reduces the stability range given in terms of the velocities difference and the electric fields effect. Finally, it is shown that the presence of surface tension enhances the stabilizing effect played by the fluid velocities, and that the kinematic viscoelasticity has a stabilizing as well as a destabilizing effect on the considered system under certain conditions. Graphics have been plotted by giving numerical values to the parameters, to depict the stability characteristics. Received 27 March 2000; accepted for publication 3 May 2001     

Espaces d'écoulements dits « universels », suite 2Spaces of universal flows,part 2

At small velocities, every universal motion can be performed in a fluid whenever the constitutive equation approaches a most general linear model. Both the celebrate theorems of Appell, Cauchy, Helmholtz, Kelvin hold true when dealing with the universal motions. Here the theorem of Bernoulli admits four generalisations. At last, the Helmholtz–Rayleigh theorem about dissipation may be regarded as defining the universal motions. To cite this article: M. Bouthier, C. R. Mecanique 332 (2004).     

On co-rotational and other rate type constitutive equations

We derive expressions for the dilatational properties of suspensions of gas bubbles in incompressible fluids, using a cell model for the suspension. A cell, consisting of a gas bubble centered in a spherical shell of incompressible fluid, is subjected to a purely dilatational boundary motion and the resulting stress at the cell boundary is obtained. The same dilatational boundary motion is prescribed at the boundary of an “equivalent” cell composed of a one-phase, uniformly compressible fluid with unknown dilatational properties. By specifying that the stress at the boundary of the one-phase cell is equal to the stress at the boundary of the two-phase suspension cell, we obtain expressions for the unknown dilatational properties as a function of observable properties of the suspension. The dilatational viscosity of a suspension with a Newtonian continuous phase and the analogous properties for suspensions with non-Newtonian continuous phases are obtained as functions of the boundary motion, volume fraction of gas, and properties of the incompressible continuous phase. Results are presented for continuous phases which are Newtonian fluids, second-order fluids, and Goddard—Miller model fluids.     

Vorticity addition method

It is shown that in the general case it is not possible to propose the Lagrangian viewpoint on the vorticity evolution, which would be unique for the entire flow, using the existing analogs of the Helmholtz theorems. This is related to the fact that, as distinct from the Helmholtz theorem oneself, these analogs are valid only for nonzero vorticity zones. New analogs of the Helmholtz theorems are proposed for the general case of flows (from incompressible fluid to viscous gas). They describe the vorticity evolution at all the points including the points of nonzero vorticity.     

Processing the capillary viscometry data of fluids with yield stress

The capillary viscometer is used to measure the shear stress-shear rate relationship of a wide range of purely viscous fluids. It is however not considered as an appropriate instrument for obtaining the yield stress and the post-yield behaviour of fluids that have a yield stress. This is partly because conventional methods of processing the capillary viscometry data of purely viscous fluids cannot be applied to similar data of fluids with yield stress. The unavoidable experimental noise in the capillary data, particularly at low shear rates, also makes it difficult to obtain a reliable estimate of the yield stress from capillary data. In this investigation the problem of converting the capillary viscometry data of yield stress fluids into a shear stress-shear rate curve and a yield stress is formulated as a Volterra integral equation of the first kind. This is an ill-posed problem i.e. noise in the data will be amplified by inappropriate methods of data processing. A method, based on Tikhonov regularisation that takes into account the ill-posed nature of the problem, is then developed to solve this problem for fluids with yield stress. The performance of this method is assessed by applying it to a set of “synthetic” capillary viscometry data with added random noise and to a set of experimental data for a concentrated suspension of TiO₂ taken from the literature. In both cases Tikhonov regularisation was able to extract the complete shear properties of these fluids from capillary viscometry data alone. Received: 22 November 1999/Accepted: 17 December 1999     

The dynamics and dissolution of gas bubbles in a viscoelastic fluid

This paper reports an experimental study of the motion of dissolving and non-dissolving gas bubbles in a quiescent viscoelastic fluid. The objective of the investigation was to determine the influence of the abrupt transition in bubble velocity, which had been observed at a critical radius of approx. on the rate of mass transfer. Thus, a range of bubble sizes from an equivalent (spherical) radius of 0.2–0.4 cm was employed using CO₂ gas, and five different fluids, including one Newtonion glycerine/water solution and four viscoelastic solutions of Separan AP30 in water (0.1, 0.5, 1% by weight) and in a water/glycerine mixture.The experimental data on bubble velocity shows that the discontinuous increase with bubble volume observed previously for air bubbles in viscoelastic fluids, does not occur for dissolving CO₂ bubbles—presumably due to the continuous decrease in bubble volume. Instead, a very steep but definitely continuous transition is found. Mass transfer rates are found to be significantly enhanced by viscoelasticity, and comparison with available theoretical results shows that the increase is greater than expected for purely viscous, power-law fluids. We conclude that a fully viscoelastic constitutive model would be necessary for a successful analysis of the dissolution of a gas bubble which is translating through a (high molecular weight) polymer solution.     

Film thickness in blade coating of viscous and viscoelastic liquids

Coating of viscous and viscoelastic liquids is examined both theoretically and experimentally. A rigid blade, accurately positioned over a rotating roll, provides an experimental system in which coating thickness is measured as a function of geometric parameters. A perturbation solution to the Navier—Stokes equations yields a lubrication theory which shows agreement with the data to an extent depending on the specific geometry.The effect of a non-Newtonian viscosity is explored by adopting a purely viscous power-law model. The lubrication equations are solved by the method of Horowitz and Steidler [1], and predict an increase in coating thickness relative to the Newtonian case. Data for viscoelastic fluids show both an increase and a decrease in coating thickness compared with Newtonian liquids depending on the relative magnitude of shear thinning and elastic effects.     

Instability of the interface in two-layer flows with large viscosity contrast at small Reynolds numbers

The Kelvin–Helmholtz instability is believed to be the dominant instability mechanism for free shear flows at large Reynolds numbers. At small Reynolds numbers, a new instability mode is identified when the temporal instability of parallel viscous two fluid mixing layers is extended to current-fluid mud systems by considering a composite error function velocity profile. The new mode is caused by the large viscosity difference between the two fluids. This interfacial mode exists when the fluid mud boundary layer is sufficiently thin. Its performance is different from that of the Kelvin–Helmholtz mode. This mode has not yet been reported for interface instability problems with large viscosity contrasts.These results are essential for further stability analysis of flows relevant to the breaking up of this type of interface.     

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.     

Spatially periodic suspensions of convex particles in linear shear flows. II. Rheology

Following the purely kinematical developments of Part 1, a rigorous analysis is presented of the “almost” time-periodic low Reynolds number hydrodynamics of a spatially periodic suspension of identical convex particles in a Newtonian liquid undergoing a macroscopically homogeneous linear shear flow. By considering the case of a single particle within a unit cell of the instantaneous spatially periodic configuration, the quasistatic dynamical analysis of this infinite-particle system is effected in much the same way as for a single particle suspended in an unbounded fluid. This is accomplished via the introduction of a partitioned hydrodynamic Stokes resistance matrix, linearly relating the force, couple and stresslet on the particle in the unit cell to the translational and rotational particle-(mean) suspension slip velocities and the mean rate-of-strain of the suspension. In contrast with the unbounded fluid case for a given geometry of the individual particles, the (purely geometric) elements of the resistance matrix depend upon the instantaneous lattice configuration.These dynamic quasistatic calculations for a given instantaneous lattice conformation, in particular that for the stresslet, are then appropriately averaged over both space and time for the class of almost time-periodic, lattice-reproducing, flows discussed in Part I. (In actually performing the time average, an important distinction is drawn between the ergodic and deterministic shear processes whose kinematical basis was laid in Part I.) In turn, this averaged dynamical information is translated into knowledge of the rheological properties of the macroscopically homogeneous suspension.A rigorous asymptotic, lubrication-theory analysis is performed during the course of an illustrative calculation of the rheological properties of a concentrated suspension of almost-touching spheres in a simple shear flow. Contrary to the findings of a previous heuristic treatment of this same lubrication-theory problem—one that ignores evolutionary variations in the instantaneous geometrical configuration of the spatially periodic suspension as the shear proceeds—the time-average properties of the suspension are found to be nonsingular in the limit.Finally, brief comments are offered on potential extensions of the scheme to include nonlinear phenomena, such as nonNewtonian fluids and inertial effects.     

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations.     

Uniqueness and drag for fluids of second grade in steady motion

For incompressible fluids of second grade that are compatible with the Clausius-Duhem inequality, non-uniqueness of steady flows with small Reynolds number (i.e. creeping flows) is possible provided the ‘absorption number’ is also small. We discuss this uniqueness question, generalize a well-known theorem of Tanner concerning how solutions of the Stokes equations may be used to generate solutions of the creeping flow equations for fluids of second grade, and give a new uniqueness theorem appropriate to a class of problems for the steady creeping flow of fluids of second grade. Under the conditions for uniqueness, we obtain a simple formula for the drag force on a fixed body which is immersed in a fluid of second grade which is undergoing uniform creeping flow. For bodies with certain geometric symmetries, the non-Newtonian nature of the fluid has no effect upon the drag.     

About Clapeyron's Theorem in Linear Elasticity

We examine some elementary interpretations of the classical theorem of Clapeyron in linear elasticity theory. As we show, a straightforward application of this theorem in the purely mechanical setting leads to an apparent paradox which can be resolved by referring either to dynamics or to thermodynamics. These richer theories play an essential part in understanding the physical significance of this theorem.     

Volume integrals of ellipsoids associated with the inhomogeneous Helmholtz equation

Problems of wave phenomena in the fields of acoustics, electromagnetics and elasticity are often reduced to an integration of the inhomogeneous Helmholtz equation. Results are presented for volume integrals associated with the Helmholtz operator, ?² + α², for an ellipsoidal region. By using appropriate Taylor series expansions and the multinomial theorem, these volume integrals are obtained in series form regions r > r′ and r < r′, where r and r′ are the distances from the origin to the point of observation and the source. Derivatives of these integrals are easily evaluated. When the wavenumber approaches zero the results reduce directly to the potentials of ellipsoids of variable densities.     

Coaxial axisymmetric vortex rings: 150 years after Helmholtz

This article addresses the fascinating 150 years history of the classical Helmholtz paper that laid the foundation of the vortex dynamics. Among general theorems on vortex motion, this memoir contains the special section on circular vortex filaments and axisymmetric vortex rings, in particular. The objective of this article is both to clarify some purely mathematical questions connected with the Dyson model of coaxial vortex rings in inviscid incompressible fluid and to provide a historical overview of achievements in experimental, analytical, and numerical studies of vortex rings interactions. The model is illustrated by several examples both of regular and chaotic motion of several vortex rings in an unbounded fluid.     

Instabilités à l'interface entre fluides miscibles par forçage oscillant horizontalInstabilities at the interface between miscible fluids under a horizontal oscillating forcing

Experiments were performed in order to study the instability at the interface between miscible fluids placed in a closed tank when submitted to a horizontal and sinusoidal oscillating forcing. For a given amplitude a, the frequency f was progressively increased untill a threshold frequency above which a stationnary relief can be observed in the reference frame of the tank. The difference of density involves an oscillating shear which lead to a Kelvin–Helmholtz type instability. An experimental stability curve was determined and confirms the physical mecanism. The capillar effects which could be expected at the interface between such fluids can be neglicted. To cite this article: M. Legendre et al., C. R. Mecanique 331 (2003).     

A Helmholtz free-energy function for a Cu–Al–Ni shape memory alloy

In this paper we develop a free energy function for a Cu–Al–Ni alloy undergoing cubic–orthorhombic phase transitions. We use the irreducible Lagrangian strain polynomial invariants of the cubic austenite parent phase to construct a polynomial expansion for the Helmholtz free energy. Our expansion retains quartic terms of the strain components in order to describe the two-phase material—the cubic austenite phase and six variants of the orthorhombic martensite phase. The coefficients of the free energy polynomial function are temperature dependent and are fitted to appropriate experimental data for austenite and martensite phases from literature. The resulting Helmholtz free energy function is given by (3.20), (4.11). We examine the response predicted by the model for shear in the twinning directions.