Non-Newtonian

A separator is fed with a condensate/gas mixture. The condensate leaves the bottom of the separator, passes a flowmeter and is followed by a choke valve, after which the condensate is boiling. The flow can not be measured using the transit time method, due to the combination of short piping, the absence of a suitable injection point and the flow properties of the condensate, which is non-newtonian due to a high contents of wax particles The condensate can not be representatively sampled, as it boils upon depressuratioh... 

Since emulsion droplets are not rigid spheres, the coefficient of 0 is around 3-6 for many emulsion systems [3-5], More concentrated emulsions are non-Newtonian depends on shear rate and are thixotropic (ri decreasing with... 

Flow behaviour of polymer melts is still difficult to predict in detail. Here, we only mention two aspects. The viscosity of a polymer melt decreases with increasing shear rate. This phenomenon is called shear thinning [48]. Another particularity of the flow of non-Newtonian liquids is the appearance of stress nonnal to the shear direction [48]. This type of stress is responsible for the expansion of a polymer melt at the exit of a tube that it was forced tlirough. Shear thinning and nonnal stress are both due to the change of the chain confonnation under large shear. On the one hand, the compressed coil cross section leads to a smaller viscosity. On the other hand, when the stress is released, as for example at the exit of a tube, the coils fold back to their isotropic confonnation and, thus, give rise to the lateral expansion of the melt. 

Colloidal dispersions often display non-Newtonian behaviour, where the proportionality in equation (02.6.2) does not hold. This is particularly important for concentrated dispersions, which tend to be used in practice. Equation (02.6.2) can be used to define an apparent viscosity, happ, at a given shear rate. If q pp decreases witli increasing shear rate, tire dispersion is called shear tliinning (pseudoplastic) if it increases, tliis is known as shear tliickening (dilatant). The latter behaviour is typical of concentrated suspensions. If a finite shear stress has to be applied before tire suspension begins to flow, tliis is known as tire yield stress. The apparent viscosity may also change as a function of time, upon application of a fixed shear rate, related to tire fonnation or breakup of particle networks. Thixotropic dispersions show a decrease in q, pp with time, whereas an increase witli time is called rheopexy. 

Computer modelling provides powerful and convenient tools for the quantitative analysis of fluid dynamics and heat transfer in non-Newtonian polymer flow systems. Therefore these techniques arc routmely used in the modern polymer industry to design and develop better and more efficient process equipment and operations. The main steps in the development of a computer model for a physical process, such as the flow and deformation of polymeric materials, can be summarized as ... 

Non-Newtonian flow processes play a key role in many types of polymer engineering operations. Hence, formulation of mathematical models for these processes can be based on the equations of non-Newtonian fluid mechanics. The general equations of non-Newtonian fluid mechanics provide expressions in terms of velocity, pressure, stress, rate of strain and temperature in a flow domain. These equations are derived on the basis of physical laws and... 

In this chapter the general equations of laminar, non-Newtonian, non-isothermal, incompressible flow, commonly used to model polymer processing operations, are presented. Throughout this chapter, for the simplicity of presentation, vector notations are used and all of the equations are given in a fixed (stationary or Eulerian) coordinate system. 

Numerous examples of polymer flow models based on generalized Newtonian behaviour are found in non-Newtonian fluid mechanics literature. Using experimental evidence the time-independent generalized Newtonian fluids are divided into three groups. These are Bingham plastics, pseudoplastic fluids and dilatant fluids. 

The practical and computational complications encountered in obtaining solutions for the described differential or integral viscoelastic equations sometimes justifies using a heuristic approach based on an equation proposed by Criminale, Ericksen and Filbey (1958) to model polymer flows. Similar to the generalized Newtonian approach, under steady-state viscometric flow conditions components of the extra stress in the (CEF) model are given a.s explicit relationships in terms of the components of the rate of deformation tensor. However, in the (CEF) model stress components are corrected to take into account the influence of normal stresses in non-Newtonian flow behaviour. For example, in a two-dimensional planar coordinate system the components of extra stress in the (CEF) model are written as... 

Herschel, W.H. and Bulkley, R., 1927. See Rudraiah, N, and Kaloni, P.N. 1990. Flow of non-Newtonian fluids. In Encyclopaedia of Fluid Mechanics, Vol. 9, Chapter 1, Gulf Publishers, Houston. 

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. 

Kaye, A., 1962. Non-Newtonian Flow in Incompressible Fluids, CoA Note No, 134, College of Aeronautics, Cranfleld. 

Pearson,. I.R.A., 1994. Report on University of Wales Institute of Non-Newtonian Fluid Mechanics Mini Symposium on Continuum and Microstructural Modelling in Computational Rheology. /. Non-Newtonian Fluid Mech. 55, 203 -205. 

Townsend, P. and Webster, M. I- ., 1987. An algorithm for the three dimensional transient simulation of non-Newtonian fluid flow. In Pande, G. N. and Middleton, J. (eds). Transient Dynamic Analysis and Constitutive Laws for Engineering Materials Vul. 2, T12, Nijhoff-Holland, Swansea, pp. 1-11. 

Application of the weighted residual method to the solution of incompressible non-Newtonian equations of continuity and motion can be based on a variety of different schemes. Tn what follows general outlines and the formulation of the working equations of these schemes are explained. In these formulations Cauchy s equation of motion, which includes the extra stress derivatives (Equation (1.4)), is used to preseiwe the generality of the derivations. However, velocity and pressure are the only field unknowns which are obtainable from the solution of the equations of continuity and motion. The extra stress in Cauchy s equation of motion is either substituted in terms of velocity gradients or calculated via a viscoelastic constitutive equation in a separate step. 

U-V-P schemes belong to the general category of mixed finite element techniques (Zienkiewicz and Taylor, 1994). In these techniques both velocity and pressure in the governing equations of incompressible flow are regarded as primitive variables and are discretized as unknowns. The method is named after its most commonly used two-dimensional Cartesian version in which U, V and P represent velocity components and pressure, respectively. To describe this scheme we consider the governing equations of incompressible non-Newtonian flow (Equations (1.1) and (1.4), Chapter 1) expressed as... 

Level of enforcement of the incompressibility condition depends on the magnitude of the penalty parameter. If this parameter is chosen to be excessively large then the working equations of the scheme will be dominated by the incompressibility constraint and may become singular. On the other hand, if the selected penalty parameter is too small then the mass conservation will not be assured. In non-Newtonian flow problems, where shear-dependent viscosity varies locally, to enforce the continuity at the right level it is necessary to maintain a balance between the viscosity and the penalty parameter. To achieve this the penalty parameter should be related to the viscosity as A = Xorj (Nakazawa et al, 1982) where Ao is a large dimensionless parameter and tj is the local viscosity. The recommended value for Ao in typical polymer flow problems is about 10. ... 

The described application of Green s theorem which results in the derivation of the weak statements is an essential step in the formulation of robu.st U-V-P and penalty schemes for non-Newtonian flow problems. 

Bell, B.C. and Surana, K. S, 1994. p-version least squares finite element formulations for two-dimensional, incompressible, non-Newtonian isothermal and non-isothcmial fluid flow. hit. J. Numer. Methods Fluids 18, 127-162. 

Luo, X. L, and Tanner, R. L, 1989. A decoupled finite element streamline-upwind scheme for viscoelastic flow problems. J. Non-Newtonian Fluid Mech. 31, 143-162. 

Marchal, J. M. and Crochet, M.J., 1987. A new mixed finite element for calculating viscoelastic flow. J. Non-Newtonian Pluid Mech. 20, 77-114. 

Papaiiastasiou, T. C., Scriven, L. E. and Macoski, C. W., 1987. A finite element method for liquid with memory. J. Non-Newtonian Fluid Mech 22, 271-288. 

It is evident that application of Green s theorem cannot eliminate second-order derivatives of the shape functions in the set of working equations of the least-sc[uares scheme. Therefore, direct application of these equations should, in general, be in conjunction with C continuous Hermite elements (Petera and Nassehi, 1993 Petera and Pittman, 1994). However, various techniques are available that make the use of elements in these schemes possible. For example, Bell and Surana (1994) developed a method in which the flow model equations are cast into a set of auxiliary first-order differentia] equations. They used this approach to construct a least-sciuares scheme for non-Newtonian flow equations based on equal-order C° continuous, p-version hierarchical elements. 

Incorporation of viscosity variations in non-elastic generalized Newtonian flow models is based on using empirical rheological relationships such as the power law or Carreau equation, described in Chapter 1. In these relationships fluid viscosity is given as a function of shear rate and material parameters. Therefore in the application of finite element schemes to non-Newtonian flow, shear rate at the elemental level should be calculated and used to update the fluid viscosity. The shear rale is defined as the second invariant of the rate of deformation tensor as (Bird et at.., 1977)... 

---