Abstract
The methods that are developed in this paper for differencing the discrete ordinates equations on a triangular x-y grid are based on piecewise polynomial representations of the angular flux. The first class of methods discussed here assumes continuity of the angular flux across all triangle interfaces. A second class of methods, which is shown to be superior to the first class, allows the angular flux to be discontinuous across triangle boundaries. Numerical results illustrating the accuracy and stability of these methods are presented, aad numerical comparisons between the above two classes of methods are made. The effectiveness of a fine mesh rebalance acceleration technique
is also discussed.
Users
Please
log in to take part in the discussion (add own reviews or comments).