Analysis and performance of adjacent-cell preconditioners for accelerating multidimensional transport calculations

The formal development of the Adjacent-cell Preconditioner (AP) and its implementation in the TORT code are briefly reviewed. Based on earlier experience with diffusion type acceleration, and excellent results in slab geometry the reciprocal averaging formula is used to mix the preconditioner elements across material and mesh discontinuities. Numerical testing of the method employing the Burre Suite of Test Problems (BSTeP), a collection of 144 cases covering a wide range in parameter space, using AP, Partial Current Rebalance (PCR), and TWODANT`s Diffusion Synthetic Acceleration (DSA) is presented. While AP outperforms the other two methods for the majority of the cases included in BSTeP it consumes many more iterations than can be explained by spectral analysis of the homogeneous model problem in cases with sharp material discontinuity. In order to verify this undesirable behavior and explore potential remedies a model problem, the Periodic Horizontal Interface (PHI), is developed that permits discontinuity of nuclear properties and cell height across the interface. Fourier mode decomposition is applied to AP with the reciprocal averaging mixing formula for the PHI configuration and shown to possess a spectral radius that approaches unity as the material discontinuity gets larger. The question of whether an unconditionally stable AP exists for PHI is tackled and preliminary indications are negative. Novel preconditioners that have nontraditional cell-coupling schemes that remain stable in these regimes may have to be sought