Stabilized mixed finite element method for the M1 radiation model

In this work, we present a computational approach for the numerical simulation of thermal radiation. Radiation is modeled by solving a set of two coupled partial differential equations, the so-called M1 model. A Variational Multiscale method is developed for this system, and tested on some illustrative benchmarks. The question of dealing with heterogeneous physical properties is also considered, and treated by means of an immersed method. It combines a level-set approach for representation of interfaces, mixing laws to build effective physical properties and an anisotropic mesh adaptation process, all these ingredients leading to an accurate description of the interface.