Invariant-domain preserving IMEX schemes for the nonequilibrium Gray Radiation-Hydrodynamics equations Part I

In this work we introduce an implicit-explicit invariant-domain preserving approximation of the nonequilibrium gray radiation-hydrodynamics equations. A time and space approximation of the system is proposed using a novel split of the equations composed of three elementary subsystems, two hyperbolic and one parabolic. The approximation thus realized is proved to be consistent, conservative, invariant-domain preserving, and first-order accurate. The proposed method is a stepping stone for achieving higher-order accuracy in space and time in the forthcoming second part of this work. The method is numerically illustrated and shown to converge as advertised. This paper is dedicated to the memory of Peter Lax.