- Journal Home
- Volume 36 - 2024
- Volume 35 - 2024
- Volume 34 - 2023
- Volume 33 - 2023
- Volume 32 - 2022
- Volume 31 - 2022
- Volume 30 - 2021
- Volume 29 - 2021
- Volume 28 - 2020
- Volume 27 - 2020
- Volume 26 - 2019
- Volume 25 - 2019
- Volume 24 - 2018
- Volume 23 - 2018
- Volume 22 - 2017
- Volume 21 - 2017
- Volume 20 - 2016
- Volume 19 - 2016
- Volume 18 - 2015
- Volume 17 - 2015
- Volume 16 - 2014
- Volume 15 - 2014
- Volume 14 - 2013
- Volume 13 - 2013
- Volume 12 - 2012
- Volume 11 - 2012
- Volume 10 - 2011
- Volume 9 - 2011
- Volume 8 - 2010
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2009
- Volume 4 - 2008
- Volume 3 - 2008
- Volume 2 - 2007
- Volume 1 - 2006
Commun. Comput. Phys., 35 (2024), pp. 181-211.
Published online: 2024-01
Cited by
- BibTex
- RIS
- TXT
In this paper, we extend the unified gas kinetic particle (UGKP) method to the frequency-dependent radiative transfer equation with both absorption-emission and scattering processes. The extended UGKP method could capture the diffusion and free transport limit and provide a smooth transition in the physical and frequency space in the regime between the above two limits. The proposed scheme has the properties of asymptotic-preserving and regime-adaptive, which make it an accurate and efficient scheme in the simulation of multiscale photon transport problems. In the UGKP formulation of flux construction and distribution closure, the coefficients of the non-equilibrium free stream distribution and near-equilibrium Planck expansion are independent of the time step. Therefore, even with a large CFL number, the UGKP can preserve a physically consistent ratio of the non-equilibrium and the near-equilibrium proportion. The methodology of scheme construction is a coupled evolution of the macroscopic energy equation and the microscopic radiant intensity equation, where the numerical flux in the macroscopic energy equation and the closure in the microscopic radiant intensity equation are constructed based on the integral solution. Both numerical dissipation and computational complexity are well controlled, especially in the optically thick regime. 2D multi-thread code on a general unstructured mesh has been developed. Several numerical tests have been simulated to verify the numerical scheme and code, covering a wide range of flow regimes. The numerical scheme and code we developed are highly demanded and widely applicable in high-energy engineering applications.
}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2023-0161}, url = {http://global-sci.org/intro/article_detail/cicp/22900.html} }In this paper, we extend the unified gas kinetic particle (UGKP) method to the frequency-dependent radiative transfer equation with both absorption-emission and scattering processes. The extended UGKP method could capture the diffusion and free transport limit and provide a smooth transition in the physical and frequency space in the regime between the above two limits. The proposed scheme has the properties of asymptotic-preserving and regime-adaptive, which make it an accurate and efficient scheme in the simulation of multiscale photon transport problems. In the UGKP formulation of flux construction and distribution closure, the coefficients of the non-equilibrium free stream distribution and near-equilibrium Planck expansion are independent of the time step. Therefore, even with a large CFL number, the UGKP can preserve a physically consistent ratio of the non-equilibrium and the near-equilibrium proportion. The methodology of scheme construction is a coupled evolution of the macroscopic energy equation and the microscopic radiant intensity equation, where the numerical flux in the macroscopic energy equation and the closure in the microscopic radiant intensity equation are constructed based on the integral solution. Both numerical dissipation and computational complexity are well controlled, especially in the optically thick regime. 2D multi-thread code on a general unstructured mesh has been developed. Several numerical tests have been simulated to verify the numerical scheme and code, covering a wide range of flow regimes. The numerical scheme and code we developed are highly demanded and widely applicable in high-energy engineering applications.