Calcul de champ de vitesse PIVlab
Estimez rapidement les composantes de vitesse, la norme du vecteur, l’angle d’écoulement et la résolution spatiale d’un champ PIV à partir du déplacement mesuré entre deux images.
Calculateur
Renseignez les paramètres puis cliquez sur Calculer le champ de vitesse.
Visualisation
Guide expert du calcul de champ de vitesse avec PIVlab
Le calcul de champ de vitesse avec PIVlab repose sur une idée simple : mesurer le déplacement apparent de particules traceuses entre deux images consécutives, puis convertir ce déplacement en vitesse réelle à l’aide d’un étalonnage spatial et d’un pas de temps connu. En pratique, cette simplicité apparente cache plusieurs points critiques : choix du temps entre images, taille des fenêtres d’interrogation, niveau de recouvrement, qualité de l’illumination laser, homogénéité de l’ensemencement et traitement des vecteurs aberrants. Si vous cherchez un outil de calcul de champ de vitesse PIVlab, l’objectif n’est pas seulement d’obtenir un nombre, mais de comprendre comment ce nombre est construit, quelle est sa sensibilité aux hypothèses de départ, et comment l’interpréter correctement dans un contexte de mécanique des fluides.
Dans PIVlab, la séquence de base consiste à charger une paire d’images, définir la calibration spatiale, renseigner le temps entre les deux acquisitions, puis lancer une corrélation croisée sur des fenêtres d’interrogation. La corrélation identifie le décalage le plus probable du motif de particules. Ce décalage est généralement fourni en pixels sur les axes x et y. Le calculateur ci-dessus reprend précisément cette logique. Il vous demande le déplacement horizontal dx, le déplacement vertical dy, l’étalonnage image, le temps entre images Δt, ainsi que deux paramètres utiles pour l’interprétation du maillage : la taille de fenêtre d’interrogation et le recouvrement.
La formule de base utilisée dans le calcul
Le principe mathématique est direct :
- Convertir le déplacement en pixels en déplacement réel grâce à la calibration spatiale.
- Convertir le temps entre images en secondes.
- Diviser le déplacement réel par le temps pour obtenir une vitesse.
Si votre calibration est connue en px/mm, alors :
mm/px = 1 / (px/mm)
Ensuite :
- Vx = dx × (mm/px) / Δt
- Vy = dy × (mm/px) / Δt
- |V| = √(Vx² + Vy²)
- Angle = atan2(Vy, Vx)
Ce calcul est exactement celui qu’on exploite ensuite pour reconstruire un champ vectoriel complet. À l’échelle d’une image entière, PIVlab n’effectue pas le calcul sur un seul point, mais sur un ensemble de fenêtres réparties spatialement. Le résultat est un maillage de vecteurs vitesse qui permet de visualiser jets, couches de cisaillement, recirculations, tourbillons et gradients locaux.
Pourquoi le temps entre images est si important
Le paramètre Δt détermine directement l’amplitude du déplacement mesuré. S’il est trop court, les particules se déplacent de moins d’un pixel et le rapport signal sur bruit de la corrélation devient médiocre. S’il est trop long, les particules quittent la fenêtre ou la déformation locale du motif devient trop forte. En pratique, de nombreux expérimentateurs ciblent un déplacement moyen de quelques pixels jusqu’à environ 10 ou 15 pixels dans la fenêtre finale, selon le niveau de cisaillement et la qualité de l’image.
Interprétation de la fenêtre d’interrogation et du recouvrement
La fenêtre d’interrogation contrôle le compromis entre résolution spatiale et robustesse du calcul. Une grande fenêtre contient plus de particules et améliore la corrélation, mais elle lisse les structures fines. Une petite fenêtre révèle mieux les gradients locaux, mais augmente le risque d’erreurs. Le recouvrement agit quant à lui sur la densité du champ final. Avec 50 % de recouvrement, le pas de grille correspond approximativement à la moitié de la taille de fenêtre. Avec 75 %, on obtient un champ plus dense, mais avec une dépendance plus forte entre vecteurs voisins.
Le calculateur estime un pas de grille physique par la relation suivante :
pas ≈ fenêtre × (1 – recouvrement) × mm/px
Ce n’est pas une mesure de vitesse, mais une information très utile pour savoir si votre maillage est suffisamment fin pour capturer la physique d’intérêt. Par exemple, une couche de cisaillement étroite ou un petit vortex ne sera pas correctement représenté si votre pas de grille est plus grand que la taille caractéristique de la structure.
Exemples de réglages typiques en PIV 2D
| Cas expérimental | Fenêtre finale | Recouvrement | Déplacement cible | Objectif principal |
|---|---|---|---|---|
| Jet libre subsonique | 32 × 32 px | 50 % | 6 à 12 px | Bon compromis robustesse / résolution |
| Couche de cisaillement | 32 × 32 px | 75 % | 4 à 10 px | Maillage plus dense pour gradients élevés |
| Microfluidique | 16 × 16 px | 50 à 75 % | 2 à 6 px | Capturer de petites structures |
| Sillage derrière obstacle | 32 ou 64 px | 50 % | 5 à 15 px | Stabilité des vecteurs en turbulence modérée |
Les valeurs du tableau ne sont pas des règles absolues, mais elles reflètent des pratiques courantes en laboratoire. En recherche, on affine souvent la stratégie grâce à des passes successives, par exemple 64 px puis 32 px, avec déformation de fenêtre et validation des vecteurs aberrants. Le gain est particulièrement net lorsque l’écoulement présente une forte rotation ou une variation rapide de vitesse dans l’espace.
Qualité d’image, densité de particules et précision
Un calcul de champ de vitesse PIVlab fiable dépend au moins autant de la qualité des images que des formules de conversion. Une image PIV utile doit contenir des particules suffisamment nombreuses, mais pas au point d’engendrer une saturation visuelle ou une perte de contraste. Le diamètre apparent des particules sur l’image est également important. Dans beaucoup de configurations PIV, un diamètre image de l’ordre de 2 à 4 pixels est recherché pour un bon compromis entre intensité et précision de localisation. Si les particules sont trop petites, le signal devient bruité. Si elles sont trop grosses, la précision sub-pixel peut se dégrader.
Il faut aussi tenir compte de la profondeur de champ, de l’épaisseur de la nappe laser et du mouvement hors plan. En PIV 2D classique, tout déplacement hors du plan lumineux tend à dégrader la corrélation. Cela ne signifie pas forcément que les résultats sont inutilisables, mais simplement qu’il faut rester prudent dans l’interprétation, notamment près des zones où la tridimensionnalité est forte.
Statistiques utiles pour interpréter un champ PIV
Une fois le calcul effectué, plusieurs grandeurs peuvent être extraites :
- la vitesse moyenne locale,
- la fluctuation de vitesse,
- les profils de vitesse sur une ligne,
- la vorticité,
- le taux de cisaillement,
- la localisation des centres de vortex,
- les grandeurs intégrales comme le débit ou l’épaisseur de mélange.
Pour beaucoup d’utilisateurs de PIVlab, le premier calcul important n’est pas la carte complète, mais la conversion correcte d’un déplacement pixel en vitesse physique. C’est pourquoi ce calculateur met en avant les composantes Vx et Vy, la norme |V|, l’angle d’écoulement et une estimation du pas spatial du maillage.
| Paramètre | Valeur couramment admise | Impact si trop faible | Impact si trop élevé |
|---|---|---|---|
| Déplacement inter-image | Environ 5 à 15 px | Signal de corrélation faible, bruit relatif élevé | Perte de corrélation, biais en fort cisaillement |
| Recouvrement | 50 à 75 % | Champ moins dense | Vecteurs voisins plus corrélés entre eux |
| Diamètre image des particules | Environ 2 à 4 px | Faible intensité, précision limitée | Pic de corrélation élargi, sous-résolution |
| Fenêtre finale | 16 à 64 px selon le cas | Corrélation fragile | Lissage excessif des structures fines |
Ces statistiques synthétiques correspondent aux tendances généralement observées dans la littérature expérimentale et dans l’usage courant des logiciels de PIV. Elles permettent de poser un diagnostic rapide lorsque les résultats semblent incohérents. Si vous obtenez des vitesses très élevées alors que les vecteurs sont visiblement instables, le problème vient souvent d’un Δt mal saisi, d’une calibration erronée, ou d’une confusion entre px/mm et mm/px.
Erreurs fréquentes dans le calcul de champ de vitesse PIVlab
- Confondre px/mm et mm/px : c’est l’erreur d’échelle la plus courante. Une inversion entraîne un facteur d’erreur énorme sur les vitesses calculées.
- Saisir Δt dans la mauvaise unité : beaucoup d’expériences utilisent des microsecondes. Si vous les traitez comme des millisecondes, le résultat sera faux d’un facteur 1000.
- Choisir une fenêtre trop petite : le champ paraît détaillé, mais devient souvent plus bruité et plus sensible aux vecteurs aberrants.
- Utiliser un recouvrement insuffisant pour des structures fines : les gradients sont alors mal résolus spatialement.
- Négliger la validation des vecteurs : il faut toujours contrôler les outliers avant toute conclusion physique.
Comment valider vos résultats
La validation d’un champ PIV ne se limite pas à l’inspection visuelle. Il est recommandé de comparer les ordres de grandeur obtenus avec une estimation indépendante : débit imposé, vitesse moyenne théorique dans un canal, mesure par anémométrie, conservation de la masse, ou profils attendus d’un écoulement canonique. Un bon champ PIV doit être à la fois cohérent localement et plausible globalement. Les résultats peuvent aussi être confrontés à des références de méthodes expérimentales diffusées par des organismes reconnus comme NIST, aux ressources en mécanique des fluides de NASA Glenn Research Center, ou à des supports universitaires comme ceux diffusés sur des sites .edu.
Dans un workflow sérieux, on examine également :
- le pourcentage de vecteurs rejetés puis remplacés,
- la continuité des lignes de courant,
- la sensibilité des résultats à un léger changement de Δt,
- la stabilité du champ moyen sur plusieurs paires d’images,
- l’accord entre la physique attendue et l’orientation des vecteurs.
Quand utiliser ce calculateur
Ce calculateur est particulièrement utile dans trois situations. D’abord, lors de la préparation d’une campagne d’essais, pour estimer rapidement quel Δt donnera un déplacement adapté. Ensuite, pendant le traitement, pour vérifier qu’un vecteur moyen observé dans PIVlab correspond bien à l’ordre de grandeur physique attendu. Enfin, après le traitement, pour documenter clairement la chaîne de conversion entre déplacement pixel et vitesse réelle dans un rapport d’essai, une note technique ou un article scientifique.
En résumé, le calcul de champ de vitesse PIVlab n’est pas qu’une simple conversion numérique. C’est une étape clé qui relie l’image expérimentale à l’interprétation physique de l’écoulement. En maîtrisant les unités, l’étalonnage, le temps inter-image, la taille des fenêtres et le recouvrement, vous obtenez des champs de vitesse plus crédibles, plus exploitables et plus faciles à défendre scientifiquement.
Conseil pratique : si les vitesses affichées semblent irréalistes, revérifiez en priorité la calibration, l’unité de Δt et l’ordre de grandeur du déplacement inter-image. Ce sont les trois sources d’erreur les plus fréquentes dans les premiers traitements PIV.