Asymptotic Diffusion-Limit Accuracy of Sn Angular Differencing Schemes
In a previous paper, Morel and Montry used a Galerkin-based diffusion analysis to define a particular weighted diamond angular discretization for Sn calculations in curvilinear geometries. The weighting factors were chosen to ensure that the Galerkin diffusion approximation was preserved, which eliminated the discrete ordinates flux dip. It was also shown that the step and diamond angular differencing schemes, which both suffer from the flux dip, do not preserve the diffusion approximation in the Galerkin sense. In this paper we re-derive the Morel and Montry weighted diamond scheme using a formal asymptotic diffusion-limit analysis. The asymptotic analysis yields more information than the Galerkin analysis and demonstrates that the step and diamond schemes do in fact formally preserve the diffusion limit to leading order, while the Morel and Montry weighted diamond scheme preserves it to first order, which is required for full consistency in this limit. Nonetheless, the fact that the step and diamond differencing schemes preserve the diffusion limit to leading order suggests that the flux dip should disappear as the diffusion limit is approached for these schemes. Computational results are presented that confirm this conjecture. We further conjecture that preserving the Galerkin diffusion approximation is equivalent to preserving the asymptotic diffusion limit to first order.