A half-moment model for radiative transfer in a 3D gray medium and its reduction to a moment model for hot, opaque sources

We present in this paper a new 3D half-moment model for radiative transfer in a gray medium, called the M11/2 model, which uses maximum entropy closure. This model is a generalization to 3D of the 1D version recently proposed in (J. Comp. Phys. 180 (2002) 584). The direction space Ω is divided into two pieces, Ω+ and Ω-, in a dynamical way by the plane perpendicular to the total radiative flux, and the half moments are defined from these subspaces. The model closure and the integrations of the radiative transfer equation performed on the moving Ω± spaces are detailed. 1D planar results, which have motivated the extension of the model of (J. Comp. Phys. 180 (2002) 584) to multi-dimensions, are shown. These results are very good. The model is thereafter derived for 3D spherically symmetric geometry, where the correctness of the non-trivial border terms can be checked. Two 3D spherically symmetric problems are numerically solved in order to show the accuracy of the closure and the role of the border terms. Once again, compared to the solution obtained with a ray tracing solver, results are very good. From the 3D half-moment model, a new moment model, called M1+, is derived for the particular case of a 3D hot and opaque source radiating into a cold medium, for applications such as simulations of stellar atmospheres and fires. Two-dimensional numerical results are presented and compared to those obtained solving the RTE and with other moment models. They demonstrate the very good accuracy of the M1+ model, its good convergence properties, and better prediction compared to all other existing moment models in its domain of applicability.