Taylor-Galerkin time stepping

The described continuous penaltyf) time-stepping scheme may yield unstable results in some problems. Therefore we consider an alternative scheme which provides better numerical stability under a wide range of conditions. This scheme is based on the U-V-P method for the slightly compressible continuity equation, described in Chapter 3, Section 1.2, in conjunction with the Taylor-Galerkin time-stepping (see Chapter 2, Section 2.5). The governing equations used in this scheme are as follows... 

Extension of the streamline Petrov -Galerkin method to transient heat transport problems by a space-time least-squares procedure is reported by Nguen and Reynen (1984). The close relationship between SUPG and the least-squares finite element discretizations is discussed in Chapter 4. An analogous transient upwinding scheme, based on the previously described 0 time-stepping technique, can also be developed (Zienkiewicz and Taylor, 1994). 

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

---