158
6 Calculs d’interactions entre solides par méthode Level Set
A l’étape initiale, il y a cinq nappes de 100 corps à distance d = 0.1 (distance
des centres des deux corps les plus proches). Les interactions entre les corps
se produisent à partir de t = 1.5.
(a) t = 0.0
(b) t = 0.75
(c) t = 1.5
(d) t = 2.25
(e) t = 3.0
(f) t = 6.9
Fig. 6.6: Simulation de 400 disques rigides soumis à la gravité (la ligne blanche
correspond à ϕ = ε). La couleur de fond représente la distance à l’union des
interfaces fluide/solides. Tiré de [90].
Comme déjà signalé, un intérêt des méthodes Level Set pour traiter les
contacts est de pouvoir prendre en compte des objets non-sphériques pour
lesquels des modèles de type force centrale ne fonctionnent pas. La figure 6.8
montre la sédimentation de 30 carrés rigides de différentes tailles.
Enfin, pour confirmer que la méthode permet aussi de traiter des objets
déformables, nous présentons un cas test de 105 vésicules biologiques évoluant
dans un écoulement de Poiseuille. Le domaine de calcul est Ω = [0, 4] × [0, 2].
La taille des vésicules étant petite, les simulations sont effectuées avec une
résolution fine, sur grille de taille (2048 × 1024).
Les résultats obtenus sont représentés dans la figure 6.9, les couleurs
indiquent la valeur de l’étiquette E 0 , les lignes blanches correspondent aux
isolignes ψ 1 = 0. Lors de l’initialisation, la région occupée par les vésicules
6 Calculs d’interactions entre solides par méthode Level Set
A l’étape initiale, il y a cinq nappes de 100 corps à distance d = 0.1 (distance
des centres des deux corps les plus proches). Les interactions entre les corps
se produisent à partir de t = 1.5.
(a) t = 0.0
(b) t = 0.75
(c) t = 1.5
(d) t = 2.25
(e) t = 3.0
(f) t = 6.9
Fig. 6.6: Simulation de 400 disques rigides soumis à la gravité (la ligne blanche
correspond à ϕ = ε). La couleur de fond représente la distance à l’union des
interfaces fluide/solides. Tiré de [90].
Comme déjà signalé, un intérêt des méthodes Level Set pour traiter les
contacts est de pouvoir prendre en compte des objets non-sphériques pour
lesquels des modèles de type force centrale ne fonctionnent pas. La figure 6.8
montre la sédimentation de 30 carrés rigides de différentes tailles.
Enfin, pour confirmer que la méthode permet aussi de traiter des objets
déformables, nous présentons un cas test de 105 vésicules biologiques évoluant
dans un écoulement de Poiseuille. Le domaine de calcul est Ω = [0, 4] × [0, 2].
La taille des vésicules étant petite, les simulations sont effectuées avec une
résolution fine, sur grille de taille (2048 × 1024).
Les résultats obtenus sont représentés dans la figure 6.9, les couleurs
indiquent la valeur de l’étiquette E 0 , les lignes blanches correspondent aux
isolignes ψ 1 = 0. Lors de l’initialisation, la région occupée par les vésicules
