Croissance et coalescence de bulles dans les magmas : analyse mathématique et simulation numérique

This thesis is devoted to mathematical and numerical study of a physical problem coming from volcanology. We look at the polydisperse modeling of bubbles growth by exsolution, decompression and coalescence. A polydisperse growth model has been proposed in literature, but it takes into accountonly the volume of bubbles, which restrict the application field, because growth by exsolution also depends on the water mass in the bubbles. In order to upgrade this model, we start with a non-dimensional monodisperse description of the bubble growth by decompression and exsolution, given by a coupled ODE system and a PDE. A numerical code is proposed to solve the monodisperse problem and is currently used. After validating this code numerically and considering several limit cases, we studied the solutions of the system and defined a flux approximation to decouple the equations system. Next, we extend the polydispers model from one to two dimensions, the volume and the water mass of bubbles. A resolution of coalescence is proposed and coupled with the polydisperse growth model. The resolution of coalescence is confronted with others numericals schemes in one and two dimensions in order to validate the proposed numerical scheme. The first numerical tests applied to a physical problem give good results.Cette thèse est consacrée à l’étude mathématiques et numérique d’un problème physique issu de la volcanologie. On s’intéresse à la modélisation polydisperse de croissance de bulles par exsolution, décompression et coalescence. Un modèle de croissance polydisperse a été proposé dans la litérature, mais ne prenait en compte que le volume des bulles, ce qui restreint le domaine d’application car la croissance par exsolution dépend également de la masse d’eau présente dans la bulle. Pour améliorer ce modèle, nous sommes parti d’une description monodisperse adimensionnelle de la croissance d’une bulle par décompression et exsolution, donnée par le couplage de deux EDO et une EDP. Un code numérique est proposé pour résoudre le problème monodisperse et est actuellement utilisé. Après avoir validé numériquement ce code et considéré plusieurs cas limites, nous avons étudié les solutions du problème et défini une approximation du flux qui nous permet de découpler le système d’équations. Ensuite, nous avons étendu le modèle polydisperse de une à deux dimensions. Une résolution de la coalescence est proposée et couplée avec le modèle de croissance polydisperse. La résolution de la coalescence est confrontée à d’autres schémas numériques en une et deux dimensions afin de valider le schéma numérique proposé. Les premiers test numériques appliqués au problème physique donnent de bon résultats