Fluids viscoelastic

Depending on the method of analysis, constitutive models of viscoelastic fluids can be formulated as differential or integral equations. 

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

Other combinations of upper- and lower-convected time derivatives of the stress tensor are also used to construct constitutive equations for viscoelastic fluids. For example, Johnson and Segalman (1977) have proposed the following equation... 

Single-integral constitutive equations for viscoelastic fluids... 

Johnson, M. W. and Segalman, D., 1977. A model for viscoelastic fluid behaviour which allows non-affine deformation. J. Non-Newtonian Fluid Mech. 2, 255-270. 

Mitsoulis, E., 1990. Numerical Simulation of Viscoelastic Fluids. In Encyclopaedia of Fluid Mechanics, Vol. 9, Chapter 21, Gulf Publishers, Houston. 

In the absence of body force, the dimensionless form of the governing model equations for two-dimensional steady-state incompressible creeping flow of a viscoelastic fluid are written as... 

Crochet, M. J., 1982, Numerical simulation of die-entry and die-exit flow of a viscoelastic fluid. In Numerical Methods in Forming Processes, Pineridge Press, Swansea. 

Keunings, R., 1989. Simulation of viscoelastic fluid flow. Tn Tucker, C. L. HI (ed.), Computer Modeling for Polymer Proces.sing, Chapter 9, Hanser Publishers, Munich, pp. 403-469. 

In generalized Newtonian fluids, before derivation of the final set of the working equations, the extra stress in the expanded equations should be replaced using the components of the rate of strain tensor (note that the viscosity should also be normalized as fj = rj/p). In contrast, in the modelling of viscoelastic fluids, stress components are found at a separate step through the solution of a constitutive equation. This allows the development of a robust Taylor Galerkin/ U-V-P scheme on the basis of the described procedure in which the stress components are all found at time level n. The final working equation of this scheme can be expressed as... 

Keeping all of the flow regime conditions identical to the previous example, we now consider a finite element model based on treating silicon rubber as a viscoelastic fluid whose constitutive behaviour is defined by the following upper-convected Maxwell equation... 

Solution of the flow equations has been based on the application of the implicit 0 time-stepping/continuous penalty scheme (Chapter 4, Section 5) at a separate step from the constitutive equation. The constitutive model used in this example has been the Phan-Thien/Tanner equation for viscoelastic fluids given as Equation (1.27) in Chapter 1. Details of the finite element solution of this equation are published elsewhere and not repeated here (Hou and Nassehi, 2001). The predicted normal stress profiles along the line AB (see Figure 5.12) at five successive time steps are. shown in Figure 5.13. The predicted pattern is expected to be repeated throughout the entire domain. 

Step 5 - using updated values of viscosity and calculated vg calculate v,., 17 and p. For viscoelastic fluids also calculate the additional stress components at this step. 

For some materials the linear constitutive relation of Newtonian fluids is not accurate. Either stress depends on strain in a more complex way, or variables other than the instantaneous rate of strain must be taken into account. Such fluids are known collectively as non-Newtonian. Many different types of behavior have been observed, ranging from fluids for which the viscosity in the Navier-Stokes equation is a simple function of the shear rate to the so-called viscoelastic fluids, for which the constitutive equation is so different that the normal stresses can cause the fluid to flow in a manner opposite to that predicted for a Newtonian fluid. 

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 ... 

Many investigators beheve that the Bingham model accounts best for observations of electrorheological behavior (116,118), but other models have also been proposed (116,119). There is considerable evidence that ER materials behave as linear viscoelastic fluids while under the influence of electric field (120) thus it appears that these materials maybe thought of as elastic Bingham fluids. 

A viscoelastic material also possesses a complex dynamic viscosity, rj = rj - - iv( and it can be shown that r = G jiuj-, rj = G juj and rj = G ju), where CO is the angular frequency. The parameter Tj is useful for many viscoelastic fluids in that a plot of its absolute value Tj vs angular frequency in radians/s is often numerically similar to a plot of shear viscosity Tj vs shear rate. This correspondence is known as the Cox-Merz empirical relationship. The parameter Tj is called the dynamic viscosity and is related to G the loss modulus the parameter Tj does not deal with viscosity, but is a measure of elasticity. 

The significance of G G tan 5, Tj, and Tj is that they can be determined experimentally and used to characterize real materials. These parameters depend on frequency and temperature, and this dependence can be used to define behavior. For example, viscoelastic fluids are often characterized by log—log plots of one or more of these quantities vs the angular frequency CO, as shown in Figure 21, which illustrates the behavior of a polymer melt (149). 

Normal Stress (Weissenberg Effect). Many viscoelastic fluids flow in a direction normal (perpendicular) to the direction of shear stress in steady-state shear (21,90). Examples of the effect include flour dough climbing up a beater, polymer solutions climbing up the inner cylinder in a concentric cylinder viscometer, and paints forcing apart the cone and plate of a cone—plate viscometer. The normal stress effect has been put to practical use in certain screwless extmders designed in a cone—plate or plate—plate configuration, where the polymer enters at the periphery and exits at the axis. 

Dyna.mic Viscometer. A dynamic viscometer is a special type of rotational viscometer used for characterising viscoelastic fluids. It measures elastic as weU as viscous behavior by determining the response to both steady-state and oscillatory shear. The geometry may be cone—plate, parallel plates, or concentric cylinders parallel plates have several advantages, as noted above. 

Figure 36 is representative of creep and recovery curves for viscoelastic fluids. Such a curve is obtained when a stress is placed on the specimen and the deformation is monitored as a function of time. During the experiment the stress is removed, and the specimen, if it can, is free to recover. The slope of the linear portion of the creep curve gives the shear rate, and the viscosity is the appHed stress divided by the slope. A steep slope indicates a low viscosity, and a gradual slope a high viscosity. The recovery part of Figure 36 shows that the specimen was viscoelastic because relaxation took place and some of the strain was recovered. A purely viscous material would not have shown any recovery, as shown in Figure 16b. 

Viscoelastic fluids that are more concentrated are characteri2ed with devices that are similar to the rotational viscometers described previously. However, instead of constant rotational motion in one direction, a sinusoidal oscillatory motion is provided. Some instmments allow both viscosity and viscoelastic measurements. 

Rheometric Scientific markets several devices designed for characterizing viscoelastic fluids. These instmments measure the response of a Hquid to sinusoidal oscillatory motion to determine dynamic viscosity as well as storage and loss moduH. The Rheometric Scientific line includes a fluids spectrometer (RFS-II), a dynamic spectrometer (RDS-7700 series II), and a mechanical spectrometer (RMS-800). The fluids spectrometer is designed for fairly low viscosity materials. The dynamic spectrometer can be used to test soHds, melts, and Hquids at frequencies from 10 to 500 rad/s and as a function of strain ampHtude and temperature. It is a stripped down version of the extremely versatile mechanical spectrometer, which is both a dynamic viscometer and a dynamic mechanical testing device. The RMS-800 can carry out measurements under rotational shear, oscillatory shear, torsional motion, and tension compression, as well as normal stress measurements. Step strain, creep, and creep recovery modes are also available. It is used on a wide range of materials, including adhesives, pastes, mbber, and plastics. 

The Weissenbetg Rheogoniometer is well suited to research on homogeneous viscoelastic fluids and elastic melts. For oscillatory shear a second motor-drive mechanism is added. This allows the use of 60 frequencies in the range of 7.6 x 10 to 40 Hz at ampHtudes between 2 x 10 and 3 X 10 rad. An electronic circuit improves the precision of oscillatory measurements, particularly at frequencies neat the natural resonance frequency of the instmment itself (298). 

Fluids without any sohdlike elastic behavior do not undergo any reverse deformation when shear stress is removed, and are called purely viscous fluids. The shear stress depends only on the rate of deformation, and not on the extent of derormation (strain). Those which exhibit both viscous and elastic properties are called viscoelastic fluids. 

Steady state, fuUy developed laminar flows of viscoelastic fluids in straight, constant-diameter pipes show no effects of viscoelasticity. The viscous component of the constitutive equation may be used to develop the flow rate-pressure drop relations, which apply downstream of the entrance region after viscoelastic effects have disappeared. A similar situation exists for time-dependent fluids. 

---