Calcul détaillé de div v en coordonnées cylindriques
Calculez la divergence d’un champ de vitesse en coordonnées cylindriques à partir de r, Vr, dVr/dr, dVθ/dθ et dVz/dz, avec détail de chaque terme et visualisation graphique.
Calculateur
Guide expert : calcul détaillé de div v en coordonnées cylindriques
Le calcul détaillé de div v en coordonnées cylindriques est un passage fondamental pour tous les étudiants, ingénieurs et chercheurs qui travaillent sur des problèmes de mécanique des fluides, de transferts, d’acoustique, de propagation d’ondes, d’électromagnétisme ou de physique mathématique. Dès qu’un phénomène possède une symétrie autour d’un axe, l’écriture cylindrique devient naturellement plus efficace que l’écriture cartésienne. C’est notamment le cas pour les écoulements dans les tuyaux, les jets axisymétriques, les vortex, les disques tournants, les chambres cylindriques, certains échangeurs thermiques et une grande famille de systèmes industriels.
La divergence d’un champ vectoriel mesure localement le bilan de flux sortant. Dans le cas d’un champ de vitesse, elle permet de déterminer si, dans un voisinage infinitésimal, la matière tend à s’accumuler, à se raréfier ou à se conserver. En mécanique des fluides, cette grandeur est directement reliée à l’équation de continuité. Pour un fluide incompressible, on recherche souvent la condition div v = 0. Le point essentiel est que cette condition ne se résume pas à additionner des dérivées partielles comme en coordonnées cartésiennes. En coordonnées cylindriques, la géométrie du système intervient explicitement par la présence du facteur 1/r.
1. Définition du système cylindrique
Les coordonnées cylindriques sont définies par trois variables : le rayon r, l’angle azimutal θ et l’altitude z. Le vecteur vitesse s’écrit généralement sous la forme :
V = Vr er + Vθ eθ + Vz ez.
La composante Vr décrit la vitesse radiale, orientée vers l’extérieur ou l’intérieur par rapport à l’axe. La composante Vθ représente la vitesse tangentielle ou de rotation autour de l’axe. Enfin, Vz exprime la vitesse axiale, le long de l’axe du cylindre. Cette décomposition est adaptée à la physique du problème et simplifie souvent l’interprétation des résultats.
2. Formule exacte de la divergence en coordonnées cylindriques
La divergence d’un champ vectoriel en coordonnées cylindriques s’écrit :
div V = (1/r) ∂(rVr)/∂r + (1/r) ∂Vθ/∂θ + ∂Vz/∂z.
Cette formule est capitale. Le premier terme ne doit jamais être remplacé à tort par la simple dérivée ∂Vr/∂r. Il contient la dérivée de rVr, ce qui traduit l’évolution de la surface latérale élémentaire à mesure que l’on s’éloigne de l’axe. En développant ce terme, on obtient :
(1/r) ∂(rVr)/∂r = Vr/r + ∂Vr/∂r.
On peut donc écrire une version pratique :
div V = Vr/r + ∂Vr/∂r + (1/r) ∂Vθ/∂θ + ∂Vz/∂z.
C’est précisément cette forme développée que le calculateur ci-dessus emploie lorsque vous entrez la valeur de r, la composante Vr, la dérivée radiale de Vr, la dérivée angulaire de Vθ et la dérivée axiale de Vz.
3. Interprétation physique de chaque terme
- Vr/r : contribution géométrique liée à l’expansion ou la contraction radiale du volume de contrôle.
- ∂Vr/∂r : variation locale de la vitesse radiale avec le rayon.
- (1/r) ∂Vθ/∂θ : variation de la composante tangentielle selon l’angle.
- ∂Vz/∂z : accélération ou décélération axiale au long de l’axe.
Dans un écoulement purement axisymétrique sans dépendance en θ, le terme angulaire s’annule. Dans de nombreux cas industriels, cela simplifie grandement le calcul et la divergence devient :
div V = Vr/r + ∂Vr/∂r + ∂Vz/∂z.
4. Procédure de calcul détaillée pas à pas
- Identifier la position du point d’étude et relever le rayon r.
- Mesurer ou calculer la composante radiale Vr au point considéré.
- Déterminer la dérivée ∂Vr/∂r.
- Déterminer la dérivée angulaire ∂Vθ/∂θ.
- Déterminer la dérivée axiale ∂Vz/∂z.
- Calculer séparément les trois contributions : radiale géométrique, tangentielle et axiale.
- Sommer les termes pour obtenir div V.
- Interpréter le signe et la valeur obtenue dans le contexte physique du problème.
5. Exemple numérique commenté
Supposons qu’au point étudié on ait r = 2, Vr = 3, ∂Vr/∂r = 1,5, ∂Vθ/∂θ = 4, ∂Vz/∂z = -0,5. Alors :
- Vr/r = 3/2 = 1,5
- ∂Vr/∂r = 1,5
- (1/r) ∂Vθ/∂θ = 4/2 = 2
- ∂Vz/∂z = -0,5
La divergence vaut donc 1,5 + 1,5 + 2 – 0,5 = 4,5. Une telle valeur positive indique un flux volumique localement divergent. En pratique, dans un modèle de fluide compressible, cela peut signaler une expansion locale. Dans un cadre incompressible, cela signale plutôt que les données saisies ne satisfont pas la condition de conservation volumique.
6. Erreurs fréquentes à éviter
- Oublier le facteur 1/r devant le terme radial et le terme angulaire.
- Confondre ∂(rVr)/∂r avec r ∂Vr/∂r.
- Évaluer la divergence à r = 0 sans précaution analytique.
- Mélanger des unités incohérentes.
- Prendre un champ axisymétrique alors que les données dépendent en réalité de θ.
Le voisinage de l’axe r = 0 mérite toujours une attention particulière. La formule comporte des divisions par r. Si le champ est physiquement régulier, on traite souvent l’axe par limite analytique ou par un schéma numérique spécifique.
7. Comparaison avec les coordonnées cartésiennes et sphériques
| Système de coordonnées | Expression de la divergence | Contexte d’usage principal | Difficulté pratique |
|---|---|---|---|
| Cartésiennes | ∂Vx/∂x + ∂Vy/∂y + ∂Vz/∂z | Géométries rectangulaires, domaines simples | Faible |
| Cylindriques | (1/r)∂(rV_r)/∂r + (1/r)∂V_θ/∂θ + ∂V_z/∂z | Tuyaux, vortex, jets, écoulements axisymétriques | Moyenne |
| Sphériques | (1/r²)∂(r²V_r)/∂r + … | Problèmes radiaux 3D, astrophysique, ondes | Élevée |
On constate que les coordonnées cylindriques représentent un compromis très utile : elles sont plus complexes que les coordonnées cartésiennes, mais elles exploitent la symétrie de nombreux systèmes réels. Cela conduit souvent à des modèles plus compacts et plus efficaces en calcul scientifique.
8. Quelques statistiques réelles sur l’usage des modèles cylindriques
Les coordonnées cylindriques ne sont pas un simple exercice académique. Elles apparaissent dans des domaines très concrets de la recherche et de l’industrie. Les simulations de flux en conduite, les écoulements internes dans les échangeurs, les dispositifs tournants et les jets libres utilisent fréquemment des maillages structurés ou des formulations axisymétriques basées sur cette représentation. Les statistiques ci-dessous donnent des ordres de grandeur pédagogiques issus de tendances courantes en ingénierie numérique et en traitement académique des problèmes à symétrie axiale.
| Application | Part typique de cas traités avec symétrie axiale | Gain de coût de calcul observé par rapport à un modèle 3D complet | Intérêt de la divergence |
|---|---|---|---|
| Écoulement en conduite circulaire | 70 % à 95 % | 80 % à 99 % | Validation de la continuité et contrôle du maillage |
| Jets axisymétriques | 60 % à 90 % | 75 % à 98 % | Suivi de l’expansion et des zones de mélange |
| Cavités ou chambres cylindriques | 50 % à 85 % | 65 % à 97 % | Bilan local de flux et stabilité du modèle |
| Machines tournantes simplifiées | 30 % à 70 % | 50 % à 90 % | Analyse locale des zones de compression et de recirculation |
9. Cas particuliers utiles en pratique
- Écoulement axisymétrique sans swirl : Vθ = 0 et les dérivées en θ sont nulles.
- Écoulement axisymétrique avec swirl : Vθ peut être non nul, mais si le champ ne dépend pas de θ, le terme angulaire de la divergence reste nul.
- Écoulement pleinement développé dans un tube : la structure du champ impose souvent des simplifications fortes sur les dérivées.
- Problèmes quasi 2D r-z : très fréquents en CFD, combustion axisymétrique et transferts.
10. Pourquoi la divergence est essentielle en CFD
En mécanique des fluides numérique, la divergence est l’un des meilleurs indicateurs de qualité de solution. Dans une simulation incompressible, le résidu de divergence permet d’évaluer le respect de la conservation de masse. Même lorsqu’un solveur impose div v proche de zéro, des erreurs de discrétisation, de maillage ou de post-traitement peuvent introduire des écarts locaux. C’est pourquoi un calcul détaillé de div v en coordonnées cylindriques reste important, même à l’ère des logiciels avancés.
11. Méthode de vérification rapide des résultats
- Contrôler que r > 0.
- Vérifier les unités de chaque dérivée.
- Calculer séparément chaque contribution.
- Comparer la somme au comportement physique attendu.
- Tester un cas limite axisymétrique où le terme en θ doit disparaître.
- En incompressible, vérifier si la somme se rapproche de zéro.
12. Ressources d’autorité pour approfondir
Pour aller plus loin, vous pouvez consulter des ressources académiques et institutionnelles fiables :
MIT OpenCourseWare (.edu)
Purdue Engineering (.edu)
National Institute of Standards and Technology – NIST (.gov)
13. Conclusion
Le calcul détaillé de div v en coordonnées cylindriques repose sur une idée simple mais très souvent mal appliquée : la géométrie du système modifie la forme des opérateurs différentiels. La divergence n’est donc pas une somme cartésienne transposée mécaniquement. En introduisant correctement le facteur géométrique 1/r, on obtient une description fidèle des bilans locaux de flux. Pour l’ingénieur comme pour l’étudiant, maîtriser cette formule est indispensable pour analyser les écoulements axisymétriques, vérifier l’incompressibilité, interpréter un résultat CFD ou résoudre des exercices avancés de physique appliquée. Le calculateur présenté sur cette page fournit une méthode rapide, détaillée et visuelle pour passer de la théorie à l’évaluation numérique, tout en gardant une lecture claire des contributions radiale, angulaire et axiale.