A second order self-consistent IMEX method for radiation hydrodynamics

We present a second order self-consistent implicit/explicit (methods that use the combination of implicit and explicit discretizations are often referred to as IMEX (implicit/explicit) methods [2], [1] and [3]) time integration technique for solving radiation hydrodynamics problems. The operators of the radiation hydrodynamics are splitted as such that the hydrodynamics equations are solved explicitly making use of the capability of well-understood explicit schemes. On the other hand, the radiation diffusion part is solved implicitly. The idea of the self-consistent IMEX method is to hybridize the implicit and explicit time discretizations in a nonlinearly consistent way to achieve second order time convergent calculations. In our self-consistent IMEX method, we solve the hydrodynamics equations inside the implicit block as part of the nonlinear function evaluation making use of the Jacobian-free Newton Krylov (JFNK) method [5], [20] and [17]. This is done to avoid order reductions in time convergence due to the operator splitting. We present results from several test calculations in order to validate the numerical order of our scheme. For each test, we have established second order time convergence.