Méthode de Collocation
Méthodes de résolution d'équations différentielles
Réalisé Par : Encadré Par :
DAOUDI MOHAMED AFSSASS ABDERAHIM
[Link] REDOUANE
EL GHAZOULI ANWAR GHARBOU MOHAMED
BOUDNIBI ZINEB
2023/2024 1
Introduction 1 6 conclusion
Étapes de la Méthode de Simulation de deux
Collocation
2 PLAN 5 particules en
sédimentation
Calcul de la vitesse du fluide Distribution des Stokeslets
3 4
dans le cas axisymétrique sur les particules
2
Introduction
La méthode de collocation est une technique numérique utilisée pour résoudre des équations différentielles et des
problèmes aux limites
La méthode de collocation repose sur le principe de choisir des points spécifiques (appelés points de collocation)
dans le domaine de la solution, et de forcer l'équation différentielle à être satisfaite exactement à ces points. Cela
transforme le problème différentiel en un système d'équations algébriques qui peut être résolu numériquement.
La méthode de collocation est intégrée dans de nombreux logiciels de simulation pour résoudre les équations
différentielles en mécanique des fluides et en dynamique structurelle.
3
Histoire de méthode de collocation
Frazer, Jones et Skan (1937) l'ont développée comme une méthode générale pour résoudre des équations
différentielles ordinaires, utilisant une variété de fonctions de base et une distribution arbitraire des points de
collocation.
Le travail de Lanczos (1938) fut le premier à proposer un bon choix tant pour les fonctions de base que pour la
distribution des points de collocation
Ces études concernent l'application des polynômes de Tchebychev à des problèmes initiaux, les premières
applications de cette méthode pour des problèmes périodiques spatiaux par Kreiss et Oliger (1972) (appelée
méthode de Fourier) et Orszag (1972) (appelée méthode pseudo-spectrale).
Burgess & Mahajerin (1984) présentent des applications de la méthode de collocation pour l'écoulement interne
bidimensionnel.
Dabros (1985) et Malysa et al. (1986), par exemple, ont calculé des représentations de singularité pour plusieurs
types d'écoulement autour de particules sphériques et ellipsoïdales, ainsi que l'écoulement dû au mouvement d'une
sphère à proximité d'une autre sphère attachée à une paroi plane.
Brady & Bossis (1988) et Weinbaum et al. (1990) examinent diverses applications des méthodes de
collocation des limites, tandis que Kim & Karrila (1991) présentent des tests complets de précision.
4
Étapes de la Méthode de Collocation
Formulation du Problème :
Considérons une équation différentielle ordinaire (EDO) générale de la forme :
𝑳 [𝒖 ( 𝒙 ) ]= 𝒇 ( 𝒙 ) , 𝒙 𝝐 [ 𝒂 , 𝒃]
où est un opérateur différentiel, est la solution inconnue, et est une fonction donnée
Choix des Fonctions de Base :
On choisit un ensemble de fonctions de , souvent des polynômes, des splines ou des fonctions trigonométriques,
pour approximer la solution La solution approximative est alors exprimée comme une combinaison linéaire de ces
fonctions de base :
𝑵
𝒖(𝒙) 𝑵 = ∑ 𝒄 𝒋 𝝓 𝒋 ( 𝒙 ) Où sont les coefficients à déterminer
𝒋= 𝟎
5
Sélection des Points de Collocation :
On sélectionne un ensemble de N+1 points distincts dans l'intervalle [a,b]. Ces points peuvent être choisis de
différentes manières, telles que les points équidistants, ou d'autres points spéciaux qui améliorent la stabilité et la
précision.
Imposition des Conditions de Collocation :
On impose que l'équation différentielle soit satisfaite exactement aux points de collocation. Cela signifie que la
solution approximative doit satisfaire :
𝑳 [ 𝑢 𝑁 ( 𝒙 𝒊 ) ]= 𝒇 ( 𝒙 𝒊 ) , 𝒊=𝟎 , 𝟏 , … … . 𝑵
Cette condition entraîne un système d'équations algébriques pour les coefficients
6
Formation du Système Linéaire :
En substituant dans l'équation différentielle et en évaluant aux points de collocation, on obtient un système
d'équations linéaires :
( )
𝑵
𝑳 ∑ 𝒄 𝒋 𝝓 𝒋 ( 𝒙 𝒊 ) = 𝒇 ( 𝒙 𝒊 ) , 𝒊=𝟎 ,𝟏 ,… … . 𝑵
𝒋=𝟎
Ce système peut être écrit sous forme matricielle :
où A est une matrice de taille (N+1)×(N+1) avec les éléments ], c est le vecteur des coefficients , et est le vecteur
des valeurs
7
Calcul de la vitesse du fluide dans le cas axisymétrique
Pour le calcul du champ de vitesse dans le fluide on intégrant l’équation intégrale de surface
suivante :
En utilisant un système de coordonnées cylindrique (r, θ, z) avec , , , et intégrant
l’équation d’intégral de surface on obtient :
(1)
où , j = 1 …. Np sont les méridiens des sphères, Np est le nombre des particules
8
En raison de l’invariance par rotation du système autour de l’axe (Oz), le champ de vitesse est de révolution autour de cet axe
avec des composantes , indépendante de θ. L’équation (1) se réduit à :
(2)
Soit le rayon de la sphère la plus proche de la paroi (sphère 1) et celui de la sphère la plus éloignée de la paroi. On utilise la
vitesse v1 de la sphère 1 comme vitesse de référence. Pour rendre les équations intégrales sans dimension, on écrit la
contrainte sous la forme :
L’équation (2) devient :
est l’angle du point y dans le plan méridien de la sphère j
9
Distribution des Stokeslets sur les particules Dans le cas axisymétrique
Pour que les points de collocation sur les surfaces des deux sphères traduisent les interactions paroi-sphère et sphère-sphère,
la répartition de ces points doit se faire d’une manière appropriée :
1. Pour la moitié de la sphère 1 du côté de la paroi, la répartition des points de
collocation est proportionnelle à la distance qui sépare ce point de la paroi :
2. Pour la moitié de la sphère i (i > 1) du côté de la paroi 2, la répartition des points de
collocation est proportionnelle à la distance de ces points au plus proche plan
tangent à la sphère i − 1 et parallèle à la paroi :
Répartition des points de collocation sur deux sphères
10
3. Pour la sphère Np, la répartition des points de collocation est proportionnelle à la distance de ces points au plus
proche plan tangent à la sphère n − 1 et parallèle à la paroi (la répartition des points de collocation est comme dans
(2), avec p = M ).
Par exemple pour le cas de deux sphères identiques de rayon unité, on obtient par intégration
sur la surface S1 de la sphère 1 la répartition telle que :
et sur la surface S2 de la sphère 2 la répartition telle que :
11
Simulation Ansys de deux particules en sédimentation
Cette simulation étudie le comportement de
sédimentation de deux particules en utilisant le
logiciel Ansys. Les particules commencent à
partir de leurs positions initiales aux
coordonnées (0, -1.5, 10) et (0, 1.5, 10).
L'objectif est d'analyser les trajectoires de ces
particules et d'observer les champs de
vitesse autour d'elles.
Pour la composante de la vitesse 𝑉𝑦
12
Simulation Ansys de deux particules en sédimentation
Cette simulation étudie le comportement de
sédimentation de deux particules en utilisant le
logiciel Ansys. Les particules commencent à
partir de leurs positions initiales aux
coordonnées (0, -1.5, 10) et (0, 1.5, 10).
L'objectif est d'analyser les trajectoires de ces
particules et d'observer les champs de
vitesse autour d'elles.
Pour la composante de la vitesse 𝑉z
13
Conclusion
la méthode de collocation est une technique puissante et précise pour résoudre
une large gamme de problèmes d'équations différentielles. Elle offre une grande
précision et efficacité, surtout avec un choix judicieux des points de collocation.
Cependant, sa mise en œuvre peut être complexe, exigeante en calculs, et peut
présenter des problèmes de stabilité numérique. Elle n'est pas toujours la
meilleure option pour les problèmes avec des solutions non régulières ou des
discontinuités. Malgré ces limitations, la méthode de collocation reste une
approche précieuse pour de nombreux problèmes mathématiques et scientifiques.
14
Merci Pour Votre Attention