Stabilizing the Consistent Quasidiffusion Method with Linear Prolongation

TL;DR

This paper introduces the lpCQD method, enhancing the stability of the consistent quasidiffusion approach for neutron transport calculations by incorporating linear prolongation techniques.

Contribution

It proposes the lpCQD method, combining the CQD approach with linear prolongation to improve stability in optically thick spatial cells.

Findings

Fourier analysis shows spectral radius similar to NDA

lpCQD improves stability over CQD in thick cells

Method maintains scalar flux consistency

Abstract

The quasidiffusion (QD) method, also known as the Variable Eddington Factor (VEF) method in the astrophysical community, is an established iterative method for accelerating source iterations in SN calculations. A great advantage of the QD method is that the diffusion equation that accelerates the SN source iterations can be discretized in any valid discretization without concern for consistency with the transport discretization. QD has comparable effectiveness with diffusion synthetic acceleration (DSA), but the converged scalar flux of the diffusion equation will differ from the transport solution by the spatial truncation errors. Larsen et al. introduced a new consistent QD method (CQD), which includes a straightforwardly defined transport consistency factor closely related to the well-known coarse mesh finite difference (CMFD) and DSA methods. The CQD method preserves the discretized scalar flux solution of the SN equations, and it is stable for problems with optically thin spatial cells, but just like nonlinear diffusion acceleration (NDA), it degrades in performance and eventually becomes unstable when the spatial cells become greater than about one mean free path thick. In this paper, we performed a formal Fourier analysis of the CQD method to show that its theoretical spectral radius is essentially the same as that of the NDA method. To improve the stability of CQD, we introduce the lpCQD method, which adopts the idea of the linear prolongation CMFD (lpCMFD) method.