Construction et analyse numérique de schéma asymptotic preserving sur maillages non structurés. Application au transport linéaire et aux systèmes de Friedrichs

摘要

The transport equation in highly scattering regimes has a limit in which the dominant behavior is given by the solution of a diffusion equation. The angular discretizations like the discrete ordinate method Sn or the truncated spherical harmonic expansion Pn have the same property. For such systems it would be interesting to construct finite volume schemes on unstructured meshes which have the same dominant behavior even if the mesh is coarse (these schemes are called asymptotic preserving schemes). Indeed these models can be coupled with Lagrangian hydrodynamics codes which generate very distorted meshes. To begin we consider the lowest order angular discretization of the transport equation that is the P1 model also called the hyperbolic heat equation. After an introduction of 1D methods, we start by modify the classical edge scheme with the Jin-Levermore procedure, this scheme is not valid in the diffusion regime because the limit diffusion scheme (Two Points Flux Approximation) is not consistent on unstructured meshes. To solve this problem we propose news schemes valid on unstructured meshes. These methods are based on the nodal scheme (GLACE scheme) designed for the acoustic and dynamic gas problems, coupled with the Jin-Levermore procedure. We obtain two schemes valid on unstructured meshes which give in 1D on the Jin-Levermore scheme and the Gosse-Toscani scheme. The limit diffusion scheme obtained is a new nodal scheme. Convergence and stability proofs have been exhibited for these schemes. In a second time, these methods have been extended to higher order angular discretisation like the Pn and Sn models using a splitting strategy between the lowest order angular discretization and the higher order angular discretization. To finish we will propose to study the discretization of the absorption/emision problem in radiative transfer and a non-linear moment model called M1 model. To treat the M1 model we propose to use a formulation like a dynamic gas system coupled with a Lagrange+remap nodal scheme and the Jin-Levermore method. The numerical method obtained preserve the asymptotic limit, the maximum principle, and the entropy inequality on unstructured meshes.