Element schemes

In a first discretization step, we apply a suitable spatial discretization to Schrodinger s equation, e.g., based on pseudospectral collocation [15] or finite element schemes. Prom now on, we consider tjj, T, V and H as denoting the corresponding vector and matrix representations, respectively. The total... 

Development of weighted residual finite element schemes that can yield stable solutions for hyperbolic partial differential equations has been the subject of a considerable amount of research. The most successful outcome of these attempts is the development of the streamline upwinding technique by Brooks and Hughes (1982). The basic concept in the streamline upwinding is to modify the weighting function in the Galerkin scheme as... 

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 following section representative examples of the development of finite element schemes for most commonly used differential and integral viscoelastic models are described. 

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. 

Finite element schemes for the integral constitutive models... 

Therefore the viscoelastic extra stress acting on a fluid particle is found via an integral in terms of velocities and velocity gradients evalua ted upstream along the streamline passing through its current position. This expression is used by Papanastasiou et al. (1987) to develop a finite element scheme for viscoelastic flow modelling. 

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

Nguen, N. and Reynen, J., 1984. A space-time least-squares finite element scheme for advection-diffusion equations. Cornput. Methods Appl Mech. Eng. 42, 331- 342. 

In this section the governing Stokes flow equations in Cartesian, polar and axisymmetric coordinate systems are presented. The equations given in two-dimensional Cartesian coordinate systems are used to outline the derivation of the elemental stiffness equations (i.e. the working equations) of various finite element schemes. 

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

WORKING EQUATIONS OF THE FINITE ELEMENT SCHEMES hence the stress terms in Equation (4,138) can be ignored to obtain... 

Temperature variations are found by the solution of the energy equation. I he finite element scheme used in this example is based on the implicit 0 time-stepping/continuous penalty scheme described in detail in Chapter 4, Section 5. 

The required working equations are derived by application of the following finite element schemes to the described governing model ... 

With this approach, when an element becomes severly distorted, it is eliminated from the computational grid and becomes a free mass point. Clearly, care must be taken to avoid eliminating elements that could potentially influence the problem at some later time. An example of a three-dimensional Lagrangian calculation that uses the eroding element scheme is presented in the next section. 

Recently, silicon-tethered diastereoselective ISOC reactions have been reported, in which effective control of remote acyclic asymmetry can be achieved (Eq. 8.91).144 Whereas ISOC occur stereoselectively, INOC proceeds with significantly lower levels of diastereoselection. The reaction pathways presented in Scheme 8.28 suggest a plausible hypo thesis for the observed difference of stereocontrol. The enhanced selectivity in reactions of silyl nitronates may he due to 1,3-allylie strain. The near-linear geometry of nitrile oxides precludes such differentiating elements (Scheme 8.28). 

These schemes require the calculations of the second and mixed derivatives, which normally result in poor accuracy when the computations are performed on discrete data. For noisy data, computed values of H and K depend on the finite element scheme used to calculate the first, second, and mixed derivatives. 

---