Nodal

With the density fiinctional theory, the first step in the constmction of a pseudopotential is to consider the solution for an isolated atom [27]. If the atomic wavefiinctions are known, tire pseudo-wavefiinction can be constmcted by removing the nodal stmcture of the wavefiinction. For example, if one considers a valence... 

IlyperCl hem can display molecular orbitals and the electron density ol each molecular orbital as contour plots, showing the nodal structure and electron distribution in the molecular orbitals. 

Equation (2.1) provides an approximate interpolated value for / at position x in terms of its nodal values and two geometrical functions. The geometrical functions in Equation (2.1) are called the shape functions. A simple inspection shows that (a) each function is equal to 1 at its associated node and is 0 at the other node, and (b) the sum of the shape functions is equal to 1. These functions, shown in Figure 2.3, are written according to their associated nodes as Aa and Ab-... 

By the insertion of the nodal coordinates into Equation (2.5) nodal values of / can be found. This is shown as... 

The described direct derivation of shape functions by the formulation and solution of algebraic equations in terms of nodal coordinates and nodal degrees of freedom is tedious and becomes impractical for higher-order elements. Furthermore, the existence of a solution for these equations (i.e. existence of an inverse for the coefficients matrix in them) is only guaranteed if the elemental interpolations are based on complete polynomials. Important families of useful finite elements do not provide interpolation models that correspond to complete polynomial expansions. Therefore, in practice, indirect methods are employed to derive the shape functions associated with the elements that belong to these families. 

In this element the velocity and pressure fields are approximated using biquadratic and bi-linear shape functions, respectively, this corresponds to a total of 22 degrees of freedom consisting of 18 nodal velocity components (corner, mid-side and centre nodes) and four nodal pressures (corner nodes). 

In conjunction with the use of isoparametric elements it is necessary to express the derivatives of nodal functions in terms of local coordinates. This is a straightforward procedure for elements with C continuity and can be described as follows Using the chain rule for differentiation of functions of multiple variables, the derivative of a function in terms of local variables ij) can be expressed as... 

After the insertion of the boundaiw conditions the solution of the system of algebraic equations in this case gives the required nodal values of 7 (i.e. T2 to 7io) as... 

Lagrangian-Eulerian (ALE) method. In the ALE technique the finite element mesh used in the simulation is moved, in each time step, according to a predetermined pattern. In this procedure the element and node numbers and nodal connectivity remain constant but the shape and/or position of the elements change from one time step to the next. Therefore the solution mesh appears to move with a velocity which is different from the flow velocity. Components of the mesh velocity are time derivatives of nodal coordinate displacements expressed in a two-dimensional Cartesian system as... 

The unknowns in this equation are the local coordinates of the foot (i.e. and 7]). After insertion of the global coordinates of the foot found at step 6 in the left-hand side, and the global coordinates of the nodal points in a given element in the right-hand side of this equation, it is solved using the Newton-Raphson method. If the foot is actually inside the selected element then for a quadrilateral element its local coordinates must be between -1 and +1 (a suitable criteria should be used in other types of elements). If the search is not successful then another element is selected and the procedure is repeated. 

After identification of the elements that contain feet of particle trajectories the old time step values of F at the feet are found by interpolating (or extrapolating for boundary nodes) its old time step nodal values. In the example shown in Figure 3.6 the old time value of Fat the foot of the trajectory passing through A is found by interpolating its old nodal values within element (e). 

Step 2 an initial configuration representing the partially filled discretized domain is considered and an array consisting of the appropriate values of F - 1, 0.5 and 0 for nodes containing fluid, free surface boundary and air, respectively, is prepared. The sets of initial values for the nodal velocity, pressure and temperature fields in the solution domain are assumed and stored as input arrays. An array containing the boundary conditions along the external boundaries of the solution domain is prepared and stored. 

Step 8 - the new values of the nodal velocities found at the end of step 7 are used as input and the free surface equation is solved. 

In Equation (5,14), (77j ) is found by interpolating existing nodal values at the old time step and then transforming the found value to the convccted coordinate system. Calculation of the componenrs of 7 " and (/7y ) depends on the evaluation of first-order derivahves of the transformed coordinates (e.g, as seen in Equation (5.9). This gives the measure of deformation experienced by the fluid between time steps of n and + 1. Using the I line-independent local coordinates of a fluid particle (, ri) we have... 

Step 1 - assume a Newtonian flow and obtain the nodal pressures. 

The described algorithm may not yield a converged solution in particular for values of power law index less than 0.5. To ensure convergence, in the iteration cycle (h + 1) for updating of the nodal pressures, an initial value found by... 

It is important to note that finite element computations on multi-block grids involving a discontinuous interface are not straightforward and special arrangements for the transformation of nodal data across the internal boundaries are required. 

STRESS. Applies the variational recovery method to calculate nodal values of pressure and, components of the stress. A mass lumping routine is called by STRESS to diagonalize the coefficient matrix in the equations to eliminate the... 

GETNOD Reads and echo prints nodal coordinates fonnatting should match the output generated by the pre-processor. 

SETPRM Rearranges numbers of nodal degrees of freedom to make them compatible with the velocity components at each node. For example, in a niiie-noded element allocated degree of freedom numbers for v i and vj at node n are X and X +9, respectively. 

ELEt-TENT CONIIECTIVITY ARRAY NODAL COORDINATES ARRAY... 

---