Fast transform spectral method for Poisson equation and radiative transfer equation in cylindrical coordinate system

摘要

Direct matrix operation is extremely memory-consuming to solve the basic equations of the radiation-hydrodynamics (R-HD) problems, e.g., Poisson equation and radiative transfer equation (RTE), especially in the cases of large grid number and multidimensions. In this work, the fast transform spectral method (FTSM), which requires much less memory than the direct matrix operation, is developed to avoid large dense matrix operation. The proposed method converges monotonically for arbitrary initial value and is verified via the multidimensional Poisson equations and one-dimensional RTE in cylindrical coordinate system. Benchmarks are also introduced to demonstrate the good accuracy of this method. The performance of the FTSM has been validated by comparing with the matrix multiplication transform spectral method (MMTSM). The results show that, the presented method has good robustness and accuracy, when the directly matrix operation of MMTSM is out of memory the FTSM still works well with high accuracy when the grid number is large enough. This means that the FTSM can be a better selection for the R-HD problems when large number grid is needed and especially for multidimensions.