  We present an asymptotic-preserving (AP) numerical method for solving the
three-temperature radiative transfer model, which holds significant importance
in inertial confinement fusion. A carefully designedsplitting method is
developed that can provide a general framework of extending AP schemes for the
gray radiative transport equation to the more complex three-temperature
radiative transfer model. The proposed scheme captures two important limiting
models: the three-temperature radiation diffusion equation (3TRDE) when opacity
approaches infinity and the two-temperature limit when the ion-electron
coupling coefficient goes to infinity. We have rigorously demonstrated the AP
property and energy conservation characteristics of the proposed scheme and its
efficiency has been validated through a series of benchmark tests in the
numerical part.
