Incompressible

The total fugacity, if the liquid is considered to be incompressible, isj calculated as a function of the vagor pressure by the expression ... 

In 1873, van der Waals [2] first used these ideas to account for the deviation of real gases from the ideal gas law P V= RT in which P, Tand T are the pressure, molar volume and temperature of the gas and R is the gas constant. Fie argried that the incompressible molecules occupied a volume b leaving only the volume V- b free for the molecules to move in. Fie further argried that the attractive forces between the molecules reduced the pressure they exerted on the container by a/V thus the pressure appropriate for the gas law isP + a/V rather than P. These ideas led him to the van der Waals equation of state ... 

Geometrically, Liouville s theorem means that if one follows the motion of a small phase volume in Y space, it may change its shape but its volume is invariant. In other words the motion of this volume in T space is like that of an incompressible fluid. Liouville s theorem, being a restatement of mechanics, is an important ingredient in the fomuilation of the theory of statistical ensembles, which is considered next. 

Surface waves at an interface between two innniscible fluids involve effects due to gravity (g) and surface tension (a) forces. (In this section, o denotes surface tension and a denotes the stress tensor. The two should not be coiifiised with one another.) In a hydrodynamic approach, the interface is treated as a sharp boundary and the two bulk phases as incompressible. The Navier-Stokes equations for the two bulk phases (balance of macroscopic forces is the mgredient) along with the boundary condition at the interface (surface tension o enters here) are solved for possible hamionic oscillations of the interface of the fomi, exp [-(iu + s)t + i V-.r], where m is the frequency, is the damping coefficient, s tlie 2-d wavevector of the periodic oscillation and. ra 2-d vector parallel to the surface. For a liquid-vapour interface which we consider, away from the critical point, the vapour density is negligible compared to the liquid density and one obtains the hydrodynamic dispersion relation for surface waves + s>tf. The temi gq in the dispersion relation arises from... 

These are the two components of the Navier-Stokes equation including fluctuations s., which obey the fluctuation dissipation theorem, valid for incompressible, classical fluids ... 

In two classic papers [18, 46], Calm and Flilliard developed a field theoretic extension of early theories of micleation by considering a spatially inliomogeneous system. Their free energy fiinctional, equations (A3.3.52). has already been discussed at length in section A3.3.3. They considered a two-component incompressible fluid. The square gradient approximation implied a slow variation of the concentration on the... 

Let us consider die scattered intensity from a binary incompressible mixture of two species (containing molecules of particle 1 and molecules of particle 2) as in (B 1.9.112) we can rewrite the relationship as... 

Continuum theory has also been applied to analyse tire dynamics of flow of nematics [77, 80, 81 and 82]. The equations provide tire time-dependent velocity, director and pressure fields. These can be detennined from equations for tire fluid acceleration (in tenns of tire total stress tensor split into reversible and viscous parts), tire rate of change of director in tenns of tire velocity gradients and tire molecular field and tire incompressibility condition [20]. 

This latter modified midpoint method does work well, however, for the long time integration of Hamiltonian systems which are not highly oscillatory. Note that conservation of any other first integral can be enforced in a similar manner. To our knowledge, this method has not been considered in the literature before in the context of Hamiltonian systems, although it is standard among methods for incompressible Navier-Stokes (where its time-reversibility is not an issue, however). 

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. 

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

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

Brooks, A. N, and Hughes, T. J.R., 1982. Streamline-upwind/Petrov Galerldn formulations for convection dominated hows with particular emphasis on the incompressible Navier -Stokes equations. Cornput. Methods Appl Meek Eng. 32, 199-259. 

Weighted residual finite element methods described in Chapter 2 provide effective solution schemes for incompressible flow problems. The main characteristics of these schemes and their application to polymer flow models are described in the present chapter. 

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

Further details of the BB, sometimes referred to as Ladyzhenskaya-Babuska-Brezi (LBB) condition and its importance in the numerical solution of incompressible flow equations can be found in textbooks dealing with the theoretical aspects of the finite element method (e.g. see Reddy, 1986), In practice, the instability (or checker-boarding) of pressure in the U-V-P method can be avoided using a variety of strategies. 

The main strategies for obtaining stable results by the U -V -P scheme for incompressible flow are as follows ... 

Algorithms based on the last approach usually provide more flexible schemes than the other two methods and hence are briefly discussed in here. Hughes et al. (1986) and de Sampaio (1991) developed Petrov-Galerkin schemes based on equal order interpolations of field variables that used specially modified weight functions to generate stable finite element computations in incompressible flow. These schemes are shown to be the special cases of the method described in the following section developed by Zienkiewicz and Wu (1991). 

As already explained the necessity to satisfy the BB stability condition restricts the types of available elements in the modelling of incompressible flow problems by the U-V P method. To eliminate this restriction the continuity equation representing the incompressible flow is replaced by an equation corresponding to slightly compressible fluids, given as... 

Using different types of time-stepping techniques Zienkiewicz and Wu (1991) showed that equation set (3.5) generates naturally stable schemes for incompressible flows. This resolves the problem of mixed interpolation in the U-V-P formulations and schemes that utilise equal order shape functions for pressure and velocity components can be developed. Steady-state solutions are also obtainable from this scheme using iteration cycles. This may, however, increase computational cost of the solutions in comparison to direct simulation of steady-state problems. 

The penalty method is based on the expression of pressure in terms of the incompressibility condition (i.e. the continuity equation) 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. ... 

Elimination of the pressure term from the equation of motion does not automatically yield a robust scheme for incompressible flow and it is still necessary to satisfy the BB stability condition by a suitable technique in both forms of the penalty method. 

The use of selectively reduced integration to obtain accurate non-trivial solutions for incompressible flow problems by the continuous penalty method is not robust and failure may occur. An alternative method called the discrete penalty technique was therefore developed. In this technique separate discretizations for the equation of motion and the penalty relationship (3.6) are first obtained and then the pressure in the equation of motion is substituted using these discretized forms. Finite elements used in conjunction with the discrete penalty scheme must provide appropriate interpolation orders for velocity and pressure to satisfy the BB condition. This is in contrast to the continuous penalty method in which the satisfaction of the stability condition is achieved indirectly through... 

The basic procedure for the derivation of a least squares finite element scheme is described in Chapter 2, Section 2.4. Using this procedure the working equations of the least-squares finite element scheme for an incompressible flow are derived as follows ... 

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

Using a known solution at the inlet. To provide an example for tins option, let us consider the finite element scheme described in Section 2.1. Assuming a fully developed flow at the inlet to the domain shown in Figure 3.3, v, (dvy/dy) = 0 and by the incompressibility condition (dvx/dx) - 0, x derivatives of all stress components are also zero. Therefore at the inlet the components of the equation of motion (3.25) are reduced to... 

The right-hand side of Equation (3,87) is set to zero considering that DA//Dt, DFIDt and the divergence of the velocity field in incompressible fluids are all equal to zero. Therefore, after integration Equation (3.87) yields... 

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. 

Gresho, P. M., Lee, R. L. and Sani, R. L., 1980. On the time-dependent solution of the incompressible Navier-Stokes equations in two and three dimensions. In Recent Advances in Numerical Methods in fluids, Ch. 2, Pineridge Press, Swansea, pp. 27-75. 

Kheshgi, H. S. and Scriven, L. E., 1985. Variable penalty method for finite element analysis of incompressible flow. Int. J. Numer. Methods Fluids 5, 785-803. 

Tanner, R.I. 2000. Engineering Rheology, 2nd edti, Oxford University Press, Oxford. Taylor, C., Ranee, J. and Medwell, J. O., 1985. A note on the imposition of traction boundary conditions when using FEM for solving incompressible flow problems. Comnmn. Appl. Numer. Methods 1, 113-121. 

---