Finite volume/finite element schemes for compressible two-phase flows inheterogeneous and anisotropic porous media
Volumes finis/Eléments finis pour des écoulements diphasiques compressibles en milieux poreux hétérogènes et anisotropes
Résumé
The objective of this thesis is the development and the analysis of robust and consistent numerical schemes for the approximation of compressible two-phase flow models in anisotropic and heterogeneous porous media. A particular emphasis is set on the anisotropy together with the geometric complexity of the medium. The mathematical problem is given in a system of two degenerate and coupled parabolic equations whose main variables are the nonwetting saturation and the global pressure. In view of the difficulties manifested in the considered system, its cornerstone equations are approximated with two different classes of the finite volume family. The first class consists of combining finite elements and finite volumes. Based on standard assumptions on the space discretization and on the permeability tensor, a rigorous convergence analysis of the scheme is carried out thanks to classical arguments. To dispense with the underlined assumptions on the anisotropy ratio and on the mesh, the model has to be first formulated in the factional flux formulation. Moreover, the diffusive term is discretized by a Godunov-like scheme while the convective fluxes are approximated using an upwind technique. The resulting scheme preserves the physical ranges of the computed solution and satisfies the coercivity property. Hence, the convergence investigation holds. Numerical results show a satisfactory qualitative behavior of the scheme even if the medium of interest is anisotropic. The second class allows to consider more general meshes and tensors. It is about a new positive nonlinear discrete duality finite volume method. The main point is to approximate a part of the fluxes using a non standard technique. The application of this ideato a nonlinear diffusion equation yields surprising results. Indeed,not only is the discrete maximum property fulfilled but also the convergence of the scheme is established. Practically, the proposed method shows great promises since it provides a positivity-preserving and convergent scheme with optimal convergence rates.
Cette thèse est centrée autour du développement et de l'analyse des schémas volumes finis robustes afin d'approcher les solutions du modèle diphasique compressible en milieux poreux hétérogènes et anisotropes. Le modèle à deux phases compressibles comprend deux équations paraboliques dégénérées et couplées dont les variables principales sont la saturation du gaz et la pression globale. Ce système est discrétisé à l'aide de deux méthodes différentes (CVFE et DDFV) qui font partie de la famille des volumes finis. La première classe à laquelle on s'intéresse consiste à combiner la méthode des volumes finis et celle des éléments finis. Dans un premier temps, on considère un schéma volume finis upwind pour la partie convective et un schéma de type éléments finis conformes pour la diffusion capillaire. Sous l'hypothèse que les coefficients de transmissibilités sont positifs, on montre que la saturation vérifie le principe du maximum et on établit des estimations d'énergies permettant de démontrer la convergence du schéma. Dans un second temps, on a mis en place un schéma positif qui corrige le précédent. Ce schéma est basé sur une approximation des flux diffusifs par le schéma de Godunov. L'avantage est d'établir la bornitude des solutions approchées ainsi que les estimations uniformes sur les gradients discrets sans aucune contrainte ni sur le maillage ni sur la perméabilité. En utilisant des arguments classiques de compacité, on prouve rigoureusement la converge du schéma. Chaque schéma est validé par des simulations numériques qui montrent bien le comportement attendu d'une telle solution. Concernant la deuxième classe, on s'intéressera tout d'abord à la construction et à l'étude d'un nouveau schéma de type DDFV (Discrete Duality Finite Volume) pour une équation de diffusion non linéaire dégénérée. Cette méthode permet d' avantage de prendre en compte des maillages très généraux et des perméabilités quelconques. L'idée clé de cette discrétisation est d'approcher les flux dans la direction normale par un schéma centré et d'utiliser un schéma décentré dans la direction tangentielle. Par conséquent, on démontre que la solution approchée respecte les bornes physiques et on établit aussi des estimations d'énergie. La convergence du schéma est également établie. Des résultats numériques confirment bien ceux de la théorie. Ils exhibent en outre que la méthode est presque d'ordre deux.
Origine : Version validée par le jury (STAR)
Loading...