Oldroyd equation

The Oldroyd-type differential constitutive equations for incompressible viscoelastic fluids can in general can be written as (Oldroyd, 1950)... 

A frequently used example of Oldroyd-type constitutive equations is the Oldroyd-B model. The Oldroyd-B model can be thought of as a description of the constitutive behaviour of a fluid made by the dissolution of a (UCM) fluid in a Newtonian solvent . Here, the parameter A, called the retardation time is de.fined as A = A (r s/(ri + s), where 7]s is the viscosity of the solvent. Hence the extra stress tensor in the Oldroyd-B model is made up of Maxwell and solvent contributions. The Oldroyd-B constitutive equation is written as... 

Oldroyd, J. G., 1947. A rational formulation of the equations of plastic flow for a Bingham solid. Proc. Camb. Philos. Soc. 43, 100-105. 

The first finite element schemes for differential viscoelastic models that yielded numerically stable results for non-zero Weissenberg numbers appeared less than two decades ago. These schemes were later improved and shown that for some benchmark viscoelastic problems, such as flow through a two-dimensional section with an abrupt contraction (usually a width reduction of four to one), they can generate simulations that were qualitatively comparable with the experimental evidence. A notable example was the coupled scheme developed by Marchal and Crochet (1987) for the solution of Maxwell and Oldroyd constitutive equations. To achieve stability they used element subdivision for the stress approximations and applied inconsistent streamline upwinding to the stress terms in the discretized equations. In another attempt, Luo and Tanner (1989) developed a typical decoupled scheme that started with the solution of the constitutive equation for a fixed-flow field (e.g. obtained by initially assuming non-elastic fluid behaviour). The extra stress found at this step was subsequently inserted into the equation of motion as a pseudo-body force and the flow field was updated. These authors also used inconsistent streamline upwinding to maintain the stability of the scheme. 

Many industrially important fluids cannot be described in simple terms. Viscoelastic fluids are prominent offenders. These fluids exhibit memory, flowing when subjected to a stress, but recovering part of their deformation when the stress is removed. Polymer melts and flour dough are typical examples. Both the shear stresses and the normal stresses depend on the history of the fluid. Even the simplest constitutive equations are complex, as exemplified by the Oldroyd expression for shear stress at low shear rates ... 

Schmidt et al. (102) carried out a detailed experimental study of PET blow molding with a well-instrumented machine and compared the results with theoretical predictions using FEM and an Oldroyd B constitutive equation. They measured and calculated internal gas pressure, coupled it with the thermomechanical inflation and performed experiments and computations with free parison inflation. 

The three constant Oldroyd model is a nonlinear constitutive equation of the differential corrotational type, such as the Zaremba-Fromm-Dewitt (ZFD) fluid (Eq. 3.3-11). [For details, see R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids, Second Edition, Vol. 1, Wiley, New York, 1987, Table 7.3-2.]... 

The quantities r and r] in equation (8.34) depend on the invariants of the tensor rik in accordance with equation (8.32). We ought to note that the behaviour of a non-linear viscoelastic liquid in a non-steady state would be different, if a dependence of the material parameters r and r] on the tensor velocity gradients or on the stress tensor is assumed. This is a point which is sometimes ignored. In any case, if r and r) are constant, equation (8.34) belongs to the class of equations introduced and investigated by Oldroyd (1950). 

Equations (9.6) and (9.7) make up the simplest set of constitutive equations for dilute polymer solutions, which, after excluding the internal variables j, can be written in the form of a differential equation that has the form of the two-constant contra-variant equation investigated by Oldroyd (1950) (Section 8.6). 

Solving the previous set of equations, especially with realistic boundary conditions, is a formidable task and a lot of issues are still unanswered. This is not surprising because of the complexity of the equations, and because of their recent derivation, around 1950 for the first nonlinear models, the Oldroyd models. On the other hand, the mathematical theory for the Euler and the Navier-Stokes equations for incompressible Newtonian fluids is still not complete though these equations were derived in 1755 md 1821 respectively ... 

More specific results can be obtained in some one dimensional situations, which we describe now. Following [47], we consider shearing motions of an Oldroyd fluid, such as Couette or Poiseuille flows. The dimensionless equations are easily reduced to a system for the shear component of the velocity w(2,<), the shear stress r x,t), and a linear combination of normal stresses x 1, < > 0,... 

Following [47] we restrict now the study of stability to Oldroyd models (where di = 0). It is easy to check that the steady Couette flow, solution of the steady equations corresponding to system (16)-(17), is given by... 

The constitutive equations are the Oldroyd-B model and a modified Oldroyd-B model in which the viscosity depends on the rate of strain. In [79], Laure et al. study the spectral stability of the plane Poiseuille flow of two viscoelastic fluids obeying an Oldroyd-B law in two configurations the first one is the two layer Poiseuille flow in the second case the same fluid occupies the symmetric upper and lower layers, surrounding the central fluid. (See Figure 9.)... 

Let us consider the set of equations governing the flow of an Oldroyd-B fluid ... 

In the story of numerical flow simulation, the ability to predict observed and significant viscoelastic flow phenomena of polymer melts and solutions in an abrupt contraction has been unsuccessful for many years, in relation to the incomplete rheological characterization of materials, especially in elongation. The numerical treatments have often been confined to flow of elastic fluids with constant viscosity, described by differential constitutive equations as the Upper Convected Maxwell and Oldroyd-B models. Fortunately, the recent possibility to use real elastic fluids with constant viscosity, the so-called Boger fluids [10], has narrowed the gap between experimental observation and numerical prediction [11]. 

It should be pointed out that the improvement of convergence might also be related to realistic preditions of shear and elongational viscosities by the Phan-Thien Tanner model, when compared to the Upper Convected Maxwell, Oldroyd-B and White-Metzner models. Satisfactory munerical results were also obtained with multi-mode integral constitutive equations using a spectnun of relaxation times [7, 17, 20-27], such as the K-BKZ model in the form introduced by Papanastasiou et al. [19]. 

Two equations have been selected for the viscoelastic extra-stress component a generalized Oldroyd-B model (GOB) and a multimode Phan-Thien Tanner model (mPTT). The veilues of the corresponding parameters are given in sub-section 3-2... 

---