Slow motions of a rigid sphere in incompressible viscoelastic fluids

Problems concerned with the force of resistance and the moment of forces acting from the side of a viscoelastic fluid on a sphere moving with acceleration are considered in a linear formulation. Fairly simple relations are obtained for a fluid with a single relaxation time or a single after-action time. A discussion of the asymptotic expressions is given for a fluid with a large number of times.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 9–16, March–April, 1976.     

Finite amplitude acceleration waves in incompressible BKZ viscoelastic fluids— The propagation condition

A general propagation condition is derived which permits the calculation of the speed of propagation of a second-order acceleration wave passing through a particular non-linear incompressible viscoelastic fluid. The viscoelastic fluid is taken to obey the Bernstein, Kearsley and Zapas single integral constitutive model. The analysis is valid for arbitrary finite amplitude waves propagating through a medium undergoing an arbitrary large deformation. Three examples, rest history, steady simple shearing flow and steady simple extensional flow, are given to demonstrate the utility of the propagation condition.     

On the fully developed tube flow of a class of non-linear viscoelastic fluids

The fully developed pipe flow of a class of non-linear viscoelastic fluids is investigated. Analytical expressions are derived for the stress components, the friction factor and the velocity field. The friction factor which depends on the Deborah and Reynolds numbers is substantially smaller than the corresponding value for the Newtonian flow field with implications concerning the volume flow rate. We show that non-affine models in the class of constitutive equations considered such as Johnson-Segalman and some versions of the Phan-Thien-Tanner models are not representative of physically realistic flow fields for all Deborah numbers. For a fixed value of the slippage factor they predict physically admissible flow fields only for a limited range of Deborah numbers smaller than a critical Deborah number. The latter is a function of the slippage.     

Observability of viscoelastic fluids

We apply the observability rank condition to study the observability of various viscoelastic fluids under imposed shear or extensional flows. In this paper the observability means the ability of determining the viscoelastic stress from the time history of the observations of the first normal stress difference. We consider four viscoelastic models: the upper convected Maxwell (UCM) model, the Phan–Thien–Tanner (PTT) model, the Johnson–Segalman (JS) model and the Giesekus model. Our study reveals that all of the four models have observability for all stress components almost everywhere under shear flow whereas under extensional flow most of the models have no observability for the shear stress component. More specifically, for UCM and JS models under imposed shear flow, the observations of the first normal stress difference allow the reconstruction of all components of viscoelastic stress. For UCM and JS models under extensional flow, the two normal stress components can be determined from the measurements of the first normal stress difference; the shear stress component does not affect the evolution of the normal stress components and consequently it cannot be extracted from the observations. Under shear flow, the PTT and Giesekus models have observability almost everywhere. That is, all components of the viscoelastic stress can be determined from the observations when the vector formed by the components of viscoelastic stress does not lie on a certain surface. Under extensional flow, the PTT model has observability almost everywhere for normal stress components whereas the Giesekus model has observability almost everywhere for all stress components. We also run simulations using the unscented Kalman filter (UKF) to reconstruct the viscoelastic stress from observations without and with noises. The UKF yields accurate and robust estimates for the viscoelastic stress both in the absence and in the presence of observation noises.     

Minimum principles for incompressible viscoelastic solids

Summary Four minimum principles are established, governing the quasi-static deformations of linear, isotropic, incompressible viscoelastic solids.

---

Sommario Si stabiliscono quattro principi di minimo che governano le deformazioni quasi statiche di solidi viscoelastici, lineari, isotropici e incompressibili.

     

Unsteady flows of inhomogeneous incompressible fluids

In this paper, we study the unsteady motion of inhomogeneous incompressible viscous fluids. We present the results corresponding to Stokes' second problem and for the flow between two parallel plates where one is oscillating.     

Nonlinear waves in incompressible viscoelastic Maxwell medium

In this paper we study two-dimensional flows of incompressible viscoelastic Maxwell media with Jaumann corotational derivative in the rheological constitutive law. In the general case, due to the incompressibility condition, the equations of motion have both real and complex characteristics. Group properties of this system are studied. On this basis, two submodels of the Maxwell model are selected, which can be reduced to hyperbolic ones. More precisely, we consider plane shear flow between two parallel planes and Couette type flow caused by the inertial cylinder rotation. As a result, we obtain the closed systems of three equations of mixed type, which describe nonlinear transverse waves in an incompressible Maxwell fluid. It is demonstrated that discontinuities can develop in elastic media even from smooth initial data. Stability of shocks in the Maxwell fluid with and without retardation time is discussed.     

Mathematical model of an incompressible viscoelastic Maxwell medium

Nonstationary motions of incompressible viscoelastic Maxwell continuum with a constant relaxation time are considered. Because in an incompressible continuous medium, pressure is not a thermodynamic variable but coincides with the stress-tensor trace to within a factor, it follows that, separating the spherical part from this tensor, one can assume that the remaining part of the stress tensor has zero trace. In the case of an incompressible medium, the equations for the velocity, pressure, and stress tensor form a closed system of first-order equations which has both real and complex characteristics, which complicates the formulation of the initial-boundary-value problem. Nevertheless, the resolvability of the Cauchy problem can be proved in the class of analytic functions. Unique resolvability of the linearized problem was established in the classes of functions of finite smoothness. The class of effectively one-dimensional motions for which the subsystem of three equations is a hyperbolic one was studied. The results of an asymptotic analysis of the latter imply the possible formation of discontinuities during the evolution of the solution. The general system of equations of motion admits an infinite-dimensional Lie pseudo-group which contains an extended Galilean group. The theorem of the invariance of the conditions on the a priori unknown free boundary was proved to obtain exact solutions of free-boundary problems. The problem of deformation of a viscoelastic strip subjected to tangential stresses applied to the free boundary is considered as an example of application of this theorem. In this problem, a scale effect of short-wave instability caused by the absence of diagonal dominance of the stress tensor deviator was found.     

Nonlinear Taylor stability of viscoelastic fluids

Using Stuart's energy method, the torque on the inner cylinder, for a second order fluid, in the supercritical regime is calculated. It is found that when the second normal stress difference is negative, the flow is more stable than for a Newtonian fluid and the torque is reduced. If the second normal stress difference is positive, then the flow is more stable and there is no torque reduction. Experimental data related to the present work are discussed.Nomenclature a amplitude of the fundamentals - A _ij ⁽¹⁾ , A _ij ⁽²⁾ first and second Rivlin-Ericksen tensors - d r ₂–r ₁ - D d/dx - E - F - g _ij metric tensor - G torque on the inner cylinder in the supercritical regime - h height of the cylinders - k ₀ /d ² - k ₁ /d ² - I ₁ - I ₂ - I ₃ - I ₄ - r ₁, r ₂ radii of inner and outer cylinders respectively - r ₀ 1/2(r ₁+r ₂) - R Reynolds number ₁ r ₁ d/ ₀ - R _c critical Reynolds number - T Taylor number r _{1 1} ² d ^{3 2}/ ₀ ² *) - T _c critical Taylor number - u ₁, v ₁, w ₁ Fundamentals of the disturbance - u _i , v _i , w _i , (i>1) harmonics - mean velocity (not laminar velocity) - u –u ₁/ar _{1 1} - v v ₁/Rar _{1 1} - x (r–r ₀)/d - , material constants - 0 viscosity - wave number d - density - ₁ angular velocity of inner cylinder - tilde denotes complex conjugate     

Numerical simulation of a developed shear turbulence

Results of two-dimensional numerical studies of turbulence that arises at the interface of two flows of poorly compressible gases are described. The results were obtained using a MAKh software system. The interrelation between spatial and time problems on the development of a turbulent zone induced by shear instability is analyzed. A constant that characterizes the degree of turbulent shear mixing is calculated. The effect of the density difference of the mixing fluids on the growth rate of the turbulence zone is studied. For all density differences considered, the coefficient of heterogeneity of the resultant mixture is evaluated. Institute for Technical Physics, Snezhinsk 456770. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 1, pp. 77–83, January–February, 2000.     

Stokes' first problem for viscoelastic fluids

The theory given in this paper is based on a generalization of Boltzmann's equation of linear viscoelasticity in which the presence of a Newtonian viscosity is acknowledged. The solution of Stokes' first problem for this kind of fluid, with a viscosity and a relaxation kernel, are derived here for the first time. The formulas given in this paper form a basis for the numerical interpretation of the idea of an effective viscosity and relaxation modulus.     

Spherical bubble collapse in viscoelastic fluids

The collapse of a spherical bubble in an infinite expanse of viscoelastic fluid is considered. For a range of viscoelastic models, the problem is formulated in terms of a generalized Bernoulli equation for a velocity potential, under the assumptions of incompressibility and irrotationality. The boundary element method is used to determine the velocity potential and viscoelastic effects are incorporated into the model through the normal stress balance across the surface of the bubble. In the case of the Maxwell constitutive equation, the model predicts phenomena such as the damped oscillation of the bubble radius in time, the almost elastic oscillations in the large Deborah number limit and the rebound limit at large values of the Deborah number. A rebound condition in terms of ReDe is derived theoretically for the Maxwell model by solving the Rayleigh–Plesset equation. A range of other viscoelastic models such as the Jeffreys model, the Rouse model and the Doi-Edwards model are amenable to solution using the same technique. Increasing the solvent viscosity in the Jeffreys model is shown to lead to increasingly damped oscillations of the bubble radius.     

Analysis of dynamical behavior of viscoelastic fluids

Summary A theoretical study is undertaken on the dynamical behaviour of viscoelastic fluids in simple shearing flow. In the present study differential-type constitutive equations are chosen in attempting to describe the dynamical behaviour of polymer solutions experimentally observed byPhilippoff. The dynamical system studied here takes into account shear rate and normal stresses as the dependent variables with shear stress kept constant during the transient period. It is found that the 3-constantOldroyd model describes qualitatively the experimentally observed behaviour.

---

Zusammenfassung Es wird theoretisch das dynamische Verhalten viskoelastischer Flüssigkeiten in einer einfachen Scherströmung untersucht. Als Versuch zur Beschreibung des dynamischen Verhaltens von Polymerlösungen, die experimentell vonPhilippoff untersucht worden waren, werden hier Zustandsgieichungen vom Differentialtyp benutzt. Das hier untersuchte dynamische System berücksichtigt die Schergeschwindigkeit und die Normalspannungen als abhängige Variable von der Schubspannung, die während der Anlaufzeit konstant gehalten wird. Es zeigt sich, daß das Dreikonstantenmodell vonOldroyd das experimentell beobachtete Verhalten qualitativ beschreibt.

     

Detecting instabilities in flows of viscoelastic fluids

There is a growing interest in developing numerical tools to investigate the onset of physical instabilities observed in experiments involving viscoelastic flows, which is a difficult and challenging task as the simulations are very sensitive to numerical instabilities. Following a recent linear stability analysis carried out in order to better understand qualitatively the origin of numerical instabilities occurring in the simulation of flows viscoelastic fluids, the present paper considers a possible extension for more complex flows. This promising method could be applied to track instabilities in complex (i.e. essentially non‐parallel) flows. In addition, results related to transient growth mechanism indicate that it might be responsible for the development of numerical instabilities in the simulation of viscoelastic fluids. Copyright © 2003 John Wiley & Sons, Ltd.     

On the dynamic response of viscoelastic fluids

Summary To define the dynamic shear behaviour of a viscoelastic fluid we require two functions — one elastic and one viscous. There are two simple alternatives based on the Voigt and Maxwell concepts which are mathematically interrelated. Current practice interprets the dynamic response of fluids as a function of angular frequency () in terms of the storage (G) and loss (G) moduli: the loss function is commonly converted to a viscosity = G/. As is well known the parameters andG are the elements of a Voigt model whereas it is near universal practice to interpret steady flow in terms of the Maxwell model. This paper shows how the interpretation of dynamic experiments on fluids in terms of the apparent Maxwell parameters is more simple, more sensitive, more consistent with steady flow behaviour, and physically more realistic.

---

Zusammenfassung Zur Festlegung des dynamischen Scherverhaltens einer viskoelastischen Flüssigkeit benötigt man zwei Funktionen, eine für das elastische und eine für das viskose Verhalten. Dafür gibt es zwei einfache, mathematisch miteinander verknüpfte Alternativen, die entweder auf dem Voigtschen oder dem Maxwellschen Konzept gründen. Die übliche Praxis beschreibt das dynamische Verhalten der Flüssigkeiten als Funktion der Winkelgeschwindigkeit durch den SpeichermodulG und den VerlustmodulG bzw. die zugeordnete Viskosität = G/. Die Parameter undG sind bekanntlich die Elemente eines Voigt-Modells, wohingegen es nahezu allgemein üblich ist, das stationäre Fließen durch ein Maxwell-Modell zu beschreiben. Diese Arbeit zeigt nun, daß die Beschreibung dynamischer Versuche mit Hilfe von apparenten Maxwell-Parametern einfacher, empfindlicher, mit dem stationären Verhalten konsistenter und in physikalischer Hinsicht realistischer ist.

---

---

Paper presented at British Society of Rheology conference. Rheometry: Methods of measurement and analysis of results, Shrivenham, April 8–11, 1975.

---

With 13 figures and 2 tables     

Flow of viscoelastic fluids between rotating disks

Few boundary-value problems in fluid mechanics can match the attention that has been accorded to the flow of fluids, Newtonian and non-Newtonian, between parallel rotating disks rotating about a common axis or about distinct axes. An interesting feature which has been recently observed is the existence of solutions that are not axially symmetric even in the case of flow due to the rotation of disks about a common axis. In this article we review the recent efforts that have been expended in the study of both symmetric and asymmetric solutions in the case of both the classical linearly viscous fluid and viscoelastic fluids.The support of the Air Force Office of Scientific Research is gratefully acknowledged.