A Free Streaming Contact Preserving Scheme for the $M_1$ Model
C. Berthon 1, J. Dubois 2, B. Dubroca 3, T.-H. Nguyen-Bui 2, R. Turpault 1*1 Universite de Nantes, Laboratoire de Mathematiques Jean Leray, 2 Rue de la Houssiniere 44322 Nantes Cedex 3, France
2 CEA, CESTA, 33114 Le Barp, France
3 Universite Bordeaux I, CELIA, 351 Cours de la liberation, 33405 Talence Cedex, France
Received 04 December 2009; Accepted (in revised version) 25 February 2010
The present work concerns the numerical approximation of the $M_1$ model for radiative transfer. The main purpose is to introduce an accurate finite volume method according to the nonlinear system of conservation laws that governs this model. We propose to derive an HLLC method which preserves the stationary contact waves. To supplement this essential property, the method is proved to be robust and to preserve the physical admissible states. Next, a relevant asymptotic preserving correction is proposed in order to obtain a method which is able to deal with all the physical regimes. The relevance of the numerical procedure is exhibited thanks to numerical simulations of physical interest.AMS subject classifications: 65M06, 85A25
Key words: Radiative transfer equation, $M_1$ model, finite volume method, Riemann solver, HLLC scheme, asymptotic preserving scheme.
Email: email@example.com (C. Berthon), Joanne.Dubois@math.u-bordeaux1.fr (J. Dubois), Bruno.Dubroca@math.u-bordeaux1.fr (B. Dubroca), Ngoc-thanh-ha.NGUYEN-BUI@CEA.FR (T.-H. Nguyen-Bui), firstname.lastname@example.org (R. Turpault)