Computing and
, defined in Equation 17–39, involves evaluations of spatial gradients for various quantities
at face centers. Applying Green-Gauss theorem the gradient of a scalar
at face center
is computed as,
(17–48) |
where the subscript represents the center of each edge of a film face,
the area of the face and
the length vector (whose magnitude is the
length of the edge and whose direction is normal to the edge), as
in the sketch below,
Clearly, at the center of the gradient calculation is how to
obtain . The following describes the steps taken
in the current implementation.
(1) Primary gradient computation
A primary gradient is calculated using Green-Gauss theorem by simply estimating edge center value as,
(17–49) |
(17–50) |
(2) Reconstruction of edge-center values
The edge-center values are then computed using either the first-order upwind or the second-order upwind schemes,
1st-order upwind:
2nd-order upwind:
is the value
of
at the UPWIND face center;
is the distance vector from the UPWIND face center
to the edge center.
(3) Final gradient calculation
Finally the gradient at face center is computed as,
(17–51) |