In subroutine initialcfd.f the following fields are calculated:
- For each node in the fluid, the distance from this node to the nearest
solid surface node. This distance is stored in field yy. It is needed for
the turbulence model.
- For each solid surface node, the distance to the nearest in-flow
node. It is stored in field xsolidsurf. This quantity is also needed for the
turbulence model.
- For each node the adjacent element height. This is the shortest distance
to from this node to all nodes belonging to elements to which the node
belongs. This height is stored in field dh and is used to determine the
local increment time
\Delta t_i.
- The value of the conservative variables in all fluid nodes starting
from the physical variables. The conservative variables, stored in field
voldaux, are
\rho \epsilon_t, \rho v_i (i=1,2,3) and \rho. The physical
variables are the static temperature T, the velocity components v_i and
the static pressue p.
The fields calculated in initialcfd frequently contain distances between
nodes, which may have changed since the last call to compfluid.