Abstract: | A positivity-preserving, conservative and entropic numerical scheme is presented for the three-temperature grey diffusion radiation hydrodynamics model. More
precisely, the dissipation matrices of the colocalized semi-Lagrangian scheme are defined in order to enforce the entropy production on each species (electron or ion) proportionally to its mass as prescribed in 34]. A reformulation of the model is then considered to enable the derivation of a robust convex combination based scheme. This
yields the positivity-preserving property at each sub-iteration of the algorithm while
the total energy conservation is reached at convergence. Numerous pure hydrodynamics and radiation hydrodynamics test cases are carried out to assess the accuracy of the
method. The question of the stability of the scheme is also addressed. It is observed
that the present numerical method is particularly robust. |