Diffusion terms

The general form of the diffusion term involves evaluating the derivative of a field variable at an integration point using shape functions. For a given node, the flow solver calculates the field variable using shape functions, which are polynomials determined by the flow solver.

The general form of the diffusion term is:

The derivatives at an integration point ip are evaluated using shape functions.

Given the (s, t, u) coordinates of the nodes within an element, the value of the field variable, ϕ, is calculated using the shape functions, as follows:

where NNODE is the number of nodes in the element.

The flow solver determines the shape functions for an element, indentifying the (s, t, u) coordinates of the element nodes. An appropriate polynomial in (s, t, u) coordinates and that the shape function corresponding to a particular node is unity at that node and zero at all other element nodes. The derivatives of the field variables are then evaluated using ϕ nodal values and derivatives of the shape functions as follows: