Search

In this paper, a unified gas kinetic particle (UGKP) method is developed for radiative transfer in both absorbing and anisotropic scattering media. This numerical method is constructed based on our theoretical work on the model reduction for an anisotropic scattering system. The macroscopic solver of this method directly solves the macroscopic anisotropic diffusion equations, eliminating the need to solve higher-order moment equations. The reconstruction of macroscopic scattering source in the microscopic solver, based on the multiscale equivalent phase function we proposed in this work, has also been simplified as one single scattering process, significantly reducing the computational costs. The proposed method has also the property of asymptotic preserving. In the optically thick regime, the proposed method solves the diffusion limit equations for an anisotropic system. In the optically thin regime, the kinetic processes of photon transport are simulated. The consistency and efficiency of the proposed method have been validated by numerical tests in a wide range of flow regimes. The novel equivalent scattering source reconstruction can be used for various transport processes, and the proposed numerical scheme is widely applicable in high-energy density engineering applications. [ABSTRACT FROM AUTHOR]

• The detailed derivations of the DSA algorithm in the r - θ - z geometry is given. • A parallel DSA algorithm is implemented in the frame of the KBA algorithm. • Based on the red-black sweep, a half-reduced PCG method has been proposed. • The effectiveness and efficiency of this parallel DSA algorithm has been verified. • The HRB(m)-PCG algorithm has well parallel scalability. A parallel DSA algorithm based on the red-black sweep PCG method is proposed to accelerate the neutron transport calculation in the 3-D cylindrical geometry. This parallel acceleration technique shares the same spatial domain decomposition strategy with the KBA method. The red-black sweep iteration is formed by moving the r - and θ -axial components of the diffusion operator into the right hand side of the diffusion equation. Using the initial values refreshed by several red-black sweep iterations simplifies the matrix-vector multiplication in the CG procedure. Through the homogeneous problems, how the times of the red-black sweep effects the DSA algorithm is analyzed. A wide range of problems including FBR benchmark, PWR Pin-by-pin calculation and PWR shielding calculation have been performed to test the effectiveness and efficiency of the new method. The solution time is reduced by 3 ~ 10 times in different problems. This new acceleration algorithm has well parallel scalability. [ABSTRACT FROM AUTHOR]