GCH2545

Écoulement d’un fluide autour d’un cylindre

Mise en contexte

En mécanique des fluides, lorsqu’on considère un écoulement bidimensionnel stationnaire, on

peut exprimer les lignes de courant (i.e. les trajectoires) d’un écoulement comme une fonction de
l’espace :

ψ(x, y) = constante

où ψ(x, y) correspond à une fonction de courant, tel qu’illustré à la figure 1.

Figure 1 – Lignes et fonction de courant

Lorsqu’on connaît ψ(x, y), on peut obtenir les composantes en x et en y (soient u et v) de la

vitesse dans tout l’écoulement. Ainsi, pour un écoulement incompressible, on a :
∂ψ ∂ψ
u= v=−
∂y ∂x
Il est possible d’obtenir ψ(x, y) pour un écoulement incompressible, non-visqueux et irrotationnel
donné en résolvant l’équation de Laplace :

∇2 ψ = 0 (1)

On peut ainsi obtenir l’écoulement autour d’une géométrie d’intérêt avec les conditions frontières
appropriées.

Équations à résoudre

On s’intéresse à l’écoulement autour d’un cylindre de rayon R, tel que présenté à la figure 2.
Pour faciliter la résolution, on peut utiliser un système de coordonnées polaires tel qu’illustré à
la figure 3.

1
 Figure 2 – Domaine de calcul Figure 3 – Domaine de calcul (en bleu) en
(en bleu) en coordonnées cartésiennes coordonnées polaires

L’équation de Laplace en coordonnées polaires est donnée par :

1 ∂ψ ∂ 2 ψ 1 ∂2ψ
+ + =0 (2)
r ∂r ∂r2 r2 ∂θ2
Les conditions frontières à traiter sont :

ψ(R, θ) = ψ(r, 0) = ψ(r, 2π) = 0

!
R2
ψ(Rext , θ) = U∞ Rext sin θ 1 − 2
Rext

où U∞ est la vitesse du fluide très loin du cylindre, lorsque ce dernier n’affecte pas l’écoulement.
Les vitesses selon les axes ⃗r et θ⃗ peuvent être obtenues des équations suivantes :
1 ∂ψ ∂ψ
vr = vθ = −
r ∂θ ∂r
Finalement, on peut calculer le coefficient de pression Cp à la surface du cylindre (r = R et
θ ∈ [0, 2π]) et obtenir les coefficients de traînée Cd et de portance Cl :
2
V

Cp = 1 −
U∞
1 2π
Z
Cd = − Cp cos θdθ
2 0
1 2π
Z
Cl = − Cp sin θdθ
2 0
où V est la norme de la vitesse.
La solution exacte de l’équation 2 avec les conditions frontières spécifiées est donnée par :
!
R2
ψexacte (r, θ) = U∞ r sin θ 1 − 2 (3)
r

On travaille en unités adimensionnelles et les valeurs des paramètres du problème sont données
dans le tableau 1.

2
 Table 1 – Valeurs à utiliser des paramètres du problème
Paramètre Valeur
U∞ 10
R 1
Rext 5

Détails sur la méthode numérique

L’objectif de ce projet est d’obtenir la fonction de courant et le champ de vitesse dans le

plan cartésien (x,y) de l’écoulement autour d’un cylindre, ainsi que les coefficients de traînée et
portance. Il faut d’abord résoudre l’équation dans le système de coordonnées polaires avec la
méthode des différences finies en utilisant des approximations d’ordre 2 afin d’obtenir
la fonction de courant dans le domaine de calcul. Par la suite, les composantes de la vitesse
vr et vθ correspondantes peut être obtenues par dérivation numérique. Finalement, on peut
obtenir les coefficient de traînée Cd et de portance Cl à l’aide d’une méthode d’intégration
numérique de votre choix.
Indice : les composantes u et v de la vitesse dans le système de coordonnées cartésiennes
peuvent être obtenues à l’aide d’une matrice de rotation 2D appliquée aux composantes vr et vθ .
Soyez conscient(e)s du coût de stockage et de calcul de votre système matriciel (en fonction du
nombre de noeuds dans votre maillage). Si vous désirez simuler pour un nombre de noeuds élevé,
pensez à utilisez une matrice creuse.

Analyses demandées

Les aspects suivants doivent être abordés dans votre rapport :

— Expliquez pourquoi il est plus intéressant d’utiliser le système de coordonnées polaires
qu’un système cartésien pour résoudre le problème.
— Vérifiez votre méthode des différences finies et son implémentation à l’aide de la solution
exacte donnée par l’équation 3. En argumentant, dites si la solution est réaliste.
— Évaluez l’influence du maillage sur l’erreur commise par votre méthode des différences
finies.
— Présentez la fonction de courant dans les systèmes de coordonnées polaires et cartésiens.
Expliquez ce qu’elle représente.
— Présentez le champs de vitesse dans le plan cartésien.
— Donnez les valeurs des coefficient de traînée Cd et de portance Cl . Commentez les coefficients
obtenus. À l’aide de ces résultats, expliquez pourquoi la résolution de l’équation 2 ne peut
être utilisée pour évaluer précisément les performances d’un avion, par exemple.

Fonction Python dont l’usage est permis

Pour réaliser votre travail, vous pouvez utiliser toutes les fonctions de production de graphiques
et d’animations de Python que vous jugez pertinent. De surcroît, vous pouvez utiliser toutes
les structures de données de Python ainsi que les fonctions utilitaire de base (génération d’une
intervalle avec linspace, etc.). Vous pouvez utiliser les librairies numpy et scipy.sparse.

3
L’usage de toute autre fonction est par défaut, interdit.

Consignes

Le projet sera a réalisé en équipe de trois. La remise de votre projet devra contenir deux
composantes :
— le logiciel Python (c-à-d le code) que vous aurez programmé,
— un rapport détaillé sur votre projet.
Nous détaillons les consignes pour chacun de ces éléments ci-après.

Logiciel Python

Votre logiciel Python sera évalué tout comme votre rapport. Il devra être lisible et sa logique
devra être facilement compréhensible par la correctrice ou le correcteur. Veuillez porter une
attention particulière aux points suivants :
— Paramétrez au début de votre code les constantes physiques, les dimensions et les données
du problème afin que votre code soit le plus générique possible. L’utilisation de valeurs
numériques directement dans le coeur de votre code doit être proscrite. Employez des unités
SI tout au long de votre travail.
— Utilisez des noms de variables significatifs.
— Employez des fonctions dans votre code. Lorsqu’une opération de plusieurs lignes est répétée
plus qu’une fois, elle devrait être dans une fonction.
— Toutes vos fonctions devraient contenir des entêtes expliquant les entrées et les sorties.
Inspirez-vous des gabarits des travaux dirigés.
— Utilisez des cas de vérification pour confirmer que vous résolvez bien le problème posé.
— N’oubliez pas en Python d’initialiser tous les vecteurs ou matrices pour des raisons de
performance.

Rapport

Vous devrez remettre un rapport concis (environ 10 pages) détaillant les éléments suivants de
votre projet :
— Mise en contexte : Expliquez le contexte autour du problème que vous résolvez. Il est
intéressant ici d’aborder, par exemple, une perspective historique. Faites preuve d’originalité.
(1/2 page)
— Équations à résoudre : Détaillez les équations à résoudre. Expliquez la signification physique
des termes présents dans l’équation si possible. (1 page)
— Méthode numérique : Détaillez la méthode numérique que vous employez pour résoudre le
problème. (1-2 page)
— Vérification : Détaillez la stratégie que vous employez pour vérifier que votre implémentation
est cohérente. Présentez des résultats numériques qui démontrent la vérification
de votre code. Par exemple, démontrez que vous obtenez bien l’ordre de convergence
souhaité, etc. (2-3 pages)
— Analyse du problème : En fonction du contexte de votre projet, plusieurs analyses sont
possibles. Cette analyse peut avoir été spécifiée dans l’énoncé du projet. Justifiez le choix
des paramètres numériques (pas de temps, pas en espace, tolérance, etc.) que vous avez
pris pour la résolution du problème. (2-3 pages)

4
 — Conclusion : Concluez sur les résultats que vous avez obtenus. Identifiez les limites de votre
modèle physique ainsi que de votre modèle numérique. Discutez des défis que vous avez
rencontrés. (1/2 page)
Pour tout le rapport, il est important que les figures soient de qualité suffisante. Leur résolution
devrait être adéquate et les figures ne doivent pas être des captures d’écran. Les figures
ne doivent donc pas présenter de contours gris ou tout autre artefact qui résulte de la prise
directe de l’image. Les figures doivent comporter des légendes et des titres aux axes appropriés.
Chaque figure doit comporter une légende (caption) sous la figure tandis que chaque tableau doit
comporter une légende au-dessus du tableau.
Votre rapport doit aussi comporter une table des matières ainsi qu’une table des figures. Le
rapport peut être réalisé à l’aide de Word ou en LATEX. Si le rapport est réalisé à l’aide de LATEX,
un facteur multiplicatif de 5% boni seront donnés à l’équipe. L’équipe d’enseignement n’offrira
aucun support outre un gabarit pour la réalisation du rapport en LATEX(si vous choisissez d’écrire
le rapport en LATEX).