Simulation numérique en réservoir engineering 2020-2021
CHAPITRE I : Résolution des systèmes linéaires par
les méthodes directe
Introduction
Définition
Propriétés
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
I.1. Méthode de Gauss
La méthode de Gauss consiste à résoudre un système linéaire en utilisant les
propriétés d’échanges et de combinaisons linéaires vues précédemment afin d’obtenir
un système triangulaire équivalent au système initial. Précisons que cette méthode est
automatique, c’est-à-dire s’applique de la même manière à tout système. On parle de
l’algorithme de Gauss.
I.1.1 Description de la méthode :
Soit à résoudre le système suivant :
On va d’abord procéder à la triangularisation du système, cette partie contient ( n-1)
étapes de transformation.
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
I.1.1. Algorithme de la méthode de Gauss
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
I.1.3 Stratégie du choix du pivot
La méthode de Gauss repose sur la formule :
Le terme s’appelle le pivot. Il est clairement indispensable que ce coefficient ne soit
pas nul. Une implémentation de cette méthode doit donc impérativement mettre en
œuvre une stratégie de choix du pivot, consistant à rechercher le plus grand pivot en
valeur absolue ou en module si l’on utilise une matrice à valeurs complexes.
Donc il faut chercher le pivot défini par (𝑴𝒂𝒙(𝒂𝒊𝒋 )) la recherche se fait donc sur
lignes et colonnes, c’est la stratégie du pivot total.
Exemple :
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
I.2.Méthode de décomposition LU
Supposons qu’on est capable d’écrire la matrice A sous la forme d’un produit de deux
matrices triangulaires L (Inferieur) et U (Supérieur) tel que :
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Exemple :
Soit à Résoudre :
On peut écrire ce système de manière plus compacte avec la matrice augmentée :
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
I.3.Matrices tri-diagonale (Méthode de Thomas)
Dans les méthodes numériques de résolution des équations différentielles aux
dérivées partielles, on rencontre des matrices à trois diagonales (principale, sous
diagonale et sur diagonale). Ce type de matrices est dit tri diagonales. L’algorithme
de résolution de ce type de système est un cas particulier de l’élimination de GAUSS.
Soit le système à matrice tri diagonale suivant :
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Exemple :
Responsable de module : Dr. Salim GAREH
Simulation numérique en réservoir engineering 2020-2021
Responsable de module : Dr. Salim GAREH
