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 ...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.展开更多
We present a new numerical method to approximate the solutions of an Euler-Poisson model,which is inherent to astrophysical flows where gravity plays an important role.We propose a discretization of gravity which ensu...We present a new numerical method to approximate the solutions of an Euler-Poisson model,which is inherent to astrophysical flows where gravity plays an important role.We propose a discretization of gravity which ensures adequate coupling of the Poisson and Euler equations,paying particular attention to the gravity source term involved in the latter equations.In order to approximate this source term,its discretization is introduced into the approximate Riemann solver used for the Euler equations.A relaxation scheme is involved and its robustness is established.The method has been implemented in the software HERACLES[29]and several numerical experiments involving gravitational flows for astrophysics highlight the scheme.展开更多
文摘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.
基金supported by the A.N.R.(Agence Nationale de la Recherche)through the projects SiNeRGHY(ANR-06-CIS6-009-01)and Anemos(ANR-11-MONU002).
文摘We present a new numerical method to approximate the solutions of an Euler-Poisson model,which is inherent to astrophysical flows where gravity plays an important role.We propose a discretization of gravity which ensures adequate coupling of the Poisson and Euler equations,paying particular attention to the gravity source term involved in the latter equations.In order to approximate this source term,its discretization is introduced into the approximate Riemann solver used for the Euler equations.A relaxation scheme is involved and its robustness is established.The method has been implemented in the software HERACLES[29]and several numerical experiments involving gravitational flows for astrophysics highlight the scheme.
