Brevet FR2671877A1 — publié le 19920724

I La présente invention concerne le domaine de la mesure de vitesse d'un fluide en écoulement instationnaire. Les mesures de vitesse sont des éléments essentiels à l'optimisation des processus industriels. On a déjà proposé de nombreuses méthodes de mesure de la vitesse d'un fluide en écoulement continu et permanent. A titre d'exemples, on peut citer des méthodes directes mettant en oeuvre des flotteurs, une chronophotographie, une vélocimétrie laser ou encore des moulinets et des méthodes indirectes mettant en oeuvre une mesure de la pression dymanique dans le fluide, une mesure de la masse volumique du fluide en écoulement, ou encore des méthodes à fil et film chauds. Dans le cas de fluides incompressibles, non visqueux, en écoulement continu et permanent, la méthode de mesure de vitesse fondée sur la mesure de la pression cinétique exploite l'équation de Bernouilli sous sa forme intégrée: ( 1) p + pgz + 1/2 p U 2 = cte avec p = pression statique locale dans le fluide, p = masse volumique locale du fluide, g = champ de la pesanteur, z = altitude, U = module de la vitesse. Pour des gaz, les variations de p gz peuvent être considérées comme nulles. La pression d'arrêt ou pression totale pi peut être exprimée sous la forme: ( 2) Pl = p + 1/2 p U 2. Le module de la vitesse peut donc être obtenu par différence entre la pression d'arrêt pi et la pression statique p ( 3) U = r 2 (pi p 1 I/2 _ Pl Cependant, les écoulements instationnaires rendent les opérations de mesure de vitesse particulièrement délicates Diverses mesures de vitesse utilisées pour des fluides en écoulement continu et permanent ne sont pas applicables à des fluides en écoulement instation- naire. En particulier, la mesure de vitesse fondée sur la mesure de la pression cinétique, qui a été rappelée ci-dessus, n'est pas exploitable pour des écoulements instationnaires En effet, l'équation de Bernoulli ( 1) rappelée précédemment n'est pas applicable à ceux-ci, puisque comme cela sera indiqué par la suite, les équations de la mécanique des fluides en écoulement instationnaire font notamment intervenir la dérivée temporelle de la vitesse. En conséquence, les méthodes couramment appliquées jusqu'ici pour mesurer la vitesse de fluides en écoulement instationnaire, consistent essentiellement en l'emploi de l'anémométrie à fil chaud ou bien de la vélocimétrie Laser à effet Doppler Cependant la mise en oeuvre de ces techniques coûteuses, délicates et sophistiquées, est peu aisée et particulièrement mal adpatée à une utilisation en site industriel. La présente invention a pour but d'éliminer les inconvénients de la technique antérieure. Ce but est atteint selon la présente invention grâce à un procédé comprenant les étapes qui consistent à: i) détecter les instants o l'accélération du fluide est nulle, ii) définir une mesure de la vitesse du fluide aux instants o l'accélération est nulle, et iii) déterminer les évolutions de la vitesse entre les instants d'accélération nulle par intégration numérique de la différence de pression entre deux pressions statiques. Comme cela sera précisé par la suite, le procédé de mesure conforme à la présente invention permet un contrôle industriel de débits fortement pulsés en alliant prix réduit, simplicité de mise en oeuvre et f iabilité. La présente invention concerne également un dispositif pour la mise en oeuvre du procédé précité. D'autres caractéristiques, buts et avantages de la présente invention apparaîtront à la lecture de la description détaillée qui va suivre et en regard des dessins annexés donnés à titre d'exemples non limitatifs et sur lesquels: la figure l représente une vue schématique d'un obstacle du type tube de Pitot équipé de trois prises de pression, conforme à un premier mode de réalisation de la présente invention, les figures 2 A, 2 B et 2 C représentent schématiquement différentes configurations d'autres modes de réalisation, la figure 3 représente une vue plus précise d'un obstacle du type tube de Pitot équipé de trois prises de pression conforme au premier mode de réalisation, la figure 4 représente un organigramme du procédé de mesure mis en oeuvre par ce premier mode de réalisation, la figure 5 représente les résultats de mesure obtenus à l'aide du premier mode de réalisation, les figures 6 et 7 illustrent schématiquement un autre mode de réalisation conforme à la présente invention, et la figure 8 représente schématiquement une variante de réalisation. Dans la suite de la description, on va tout d'abord expliquer les caractéristiques essentielles du procédé de mesure de vitesse conforme à la présente invention, avant de décrire des modes de réalisation particuliers mettant celui-ci en oeuvre. CARACTERISTIQUES ESSENTIELLES DE L'INVENTION Les équations de la mécanique des fluides en écoulement instationnaire font intervenir les dérivées temporelle et spatiale de la vitesse et le gradient de pression ( 4) av + p Par intégration spatiale de la relation ( 4) on voit que la pression peut être exprimée comme une intégrale de volume de la vitesse et de sa dérivée par rapport au temps. Les inventeurs ont ainsi déterminé qu'il était possible de connaître la vitesse locale et instantanée d'un écoulement instationnaire à partir de mesures de pression, en passant par une intégration temporelle et une propriété d'intégrale spatiale de la vitesse. Cette propriété intégrale peut être obtenue au moyen de l'un des théorèmes généraux de la mécanique des fluides: théorème des débits de quantités de mouvements, théorème de l'énergie cinétique, deuxième théorème de Bernoulli pour les écoulements à potentiel, Contrairement au cas des écoulements permanents, la présence de la dérivée temporelle de la vitesse dans la relation ( 4) montre qu'il n'est pas possible d'obtenir une relation locale et instantanée entre la vitesse et la pression, et qu'une intégration temporelle est donc nécessaire. De là, les inventeurs ont proposé le procédé conforme à la présente invention, qui comme indiqué précédemment, comprend les étapes consistant à: i) détecter les instants o l'accélération a V/ a t du fluide est nulle, ii) définir une mesure de la vitesse du fluide aux instants o l'accélération est nulle, et iii) déterminer les évolutions de la vitesse entre les instants d'accélération nulle par intégration numérique de l'écart entre deux pressions statiques. DESCRIPTIF GENERAL D'UN PREMIER MODE DE REALISATION Pour réaliser pratiquement un premier dispositif de mesure de vitesse instationnaire par capteurs de pression, les inventeurs proposent, comme représenté sur la figure 1 annexée, d'utiliser un obstacle 10 placé dans le fluide en écoulement (un obstacle analogue à un tube de Pitot par exemple) et sur lequel sont placées des prises de pression 14, 16 et 18. Si Vo(t) désigne la vitesse dans laquelle est placé cet obstacle, la pression en un point de l'obstacle vérifie une loi ayant la forme ( 5) ao v 2 p ( 5) lit + + -cte p dans laquelle O désigne le potentiel des vitesses et V la vitesse au point considéré, ces deux grandeurs étant proportionnelles à la vitesse Vo(t) extérieure à l'obstacle, les constantes de proportionnalité ne dépendant que de la position du point sur l'obstacle. En appliquant la relation ( 5) en deux points Ml et M 2 au même instant, on obtient une relation de la forme ( 6) d V b V 2: Pl P 2 O ( 6 dt 2 o a et b sont deux constantes qu'on peut calculer ou mesurer en écoulement permanent par les méthodes classiques de la mécanique des fluides. La connaissance de V O (t) à partir des mesures de p 1 P 2 peut ainsi être obtenue en théorie par l'intégration numérique de l'équation différentielle ( 6); si l'acquisition de la mesure de la différence des pressions p 1 P 2 est faite au moyen d'un système micro- informatique, le microprocesseur du sytème peut être programmé pour effectuer cet intégration numérique au fur et à mesure de l'acquisition des données. Cependant, l'intégration numérique de l'équation différentiel- le ( 6) entraîne des erreurs numériques qui peuvent conduire à une valeur erronée si un recalage régulier des valeurs obtenues pour V n'est pas effectué. Dans ce but, les inventeurs proposent d'exploiter le fait que la relation ( 6) fournit directement la valeur de la vitesse V O (t) au moment o l'accélération V O '(t) est nulle Plus précisément, on peut repérer directe- ment l'instant correspondant sur la mesure de p 1 P 2 pour pouvoir calculer la vitesse V 0 à cet instant de manière indépendante du processus d'intégration numérique de Vo(t). Ceci peut être réalisé, comme représenté sur la figure 1 annexée, en aménageant sur l'obstacle 10 une zone cylindrique 12 parallèle à la vitesse V O (t) du fluide en écoulement instationnaire et sur laquelle on a placé deux prises de pression statique 14, 16, distantes de 1 et entre lesquelles la différence de pression est égale à p 1 d V 0/dt Entre ces prises de pression, la relation ( 6) se réduit à ( 7) d V( Pl P 2 ( 7) dt pi De façon connue en soi, les prises de pression statique 14, 16 peuvent être formées de petits orifices pratiqués dans la partie cylindrique de l'obstacle 10 et reliés à des manomètres respectifs. De préférence, le système représenté sur la figure 1 comprend également une prise de pression 18 sensible à la pression dynamique au point d'arrêt Cette prise de pression d'arrêt 18 peut être formée d'un petit orifice pratiqué sur le nez de l'obstacle (formé de préférence d'une calotte hémisphérique), au point d'arrêt et relié à un manomètre. Dans ce cas, le procédé conforme à la présente invention comprend essentiellement les étapes qui consistent à: i) détecter les instants o l'accélération du fluide est nulle, par détection de l'égalité entre les deux pression statiques p 1 et P 2 mesurées par les prises 14, 16, ii) définir une mesure de la vitesse du fluide aux instants o l'accélération est nulle, sur la base de la pression d'arrêt mesurée par la prise 18 et de la pression statique p 1 mesurée par la prise 14, à l'aide de la relation ( 8) V 2 (parrêt P 1 et p iii) déterminer les évolutions de la vitesse entre les instants d'accélération nulle par intégration numérique de l'équation ( 6) précitée. Pour cette étape iii), on peut également effectuer l'intégra- tion directe de l'équation ( 7) qui est évidemment plus rapide que celle de ( 6) puisqu'il n'est pas nécessaire de calculer V 02 à chaque instant. DESCRIPTIF GENERAL D'UN SECOND MODE DE REALISATION Il est intéressant de choisir une prise de pression sensible à la pression dynamique au point d'arrêt comme représenté sur la figure 1. Cependant, on peut procéder d'autres manières: par exemple, en écoulement interne, on peut provoquer une variation de vitesse au moyen d'un rétrécissement de section comme représenté sur les figures 2 A, 2 B et 2 C annexées. On aperçoit sur la figure 2 A une conduite profilée 20 comportant une première partie 22 convergeante qui aboutit à une deuxième partie 24 cylindrique Deux prises de pression statique 25, 26 placées sur la partie cylindrique 24 mesurent des pressions p 1, P 2; Une autre prise de pression placée dans la première partie 22 mesure la pression p 3. On aperçoit sur la figure 2 B une autre conduite 30 comportant un rétrécissement 32 Deux prises de pression statique 35, 36 placées sur la partie cylindrique 34 de section constante de la conduite, mesurent des pressions p 1, P 2 Une troisième prise de pression statique 37 placée sur le rétrécissement 32 mesure la pression p 3. On aperçoit sur la figure 2 C une autre conduite 40 comportant un rétrécissement 42 Deux prises de pression statique 45, 46 placées sur le rétrécissement 42 mesurent des pressions P 1, P 2 Une troisième prise de pression statique 47 placée sur la partie la plus large de la conduite 40 mesure la pression p 3. On peut alors à l'aide de l'un des systèmes représentés sur les figures 2 A, 2 B ou 2 C mesurer la vitesse du fluide concerné à l'aide des étapes consistant à: i) détecter les instants o l'accélération du fluide est nulle par détection de l'égalité entre les deux pressions statiques P 1 i P 2 mesurées par les prises et 26 ou 35, 36 ou 45 et 46, ii) définir une mesure de la vitesse du fluide aux instants o l'accélération est nulle sur la base de la différence de pression p 1 p 3 mesurée par les prises correspondanes telles que 35 et 37 ou 45 et 47, à l'aide de la relation: ( 9) V = 2 (p 1 p 3), et p iii) déterminer les évolutions de la vitesse entre les instants d'accélération nulle par intégration de la différence Pl P 2. Toutes les géométries des dispositifs mesurant des pressions statiques ou totales en écoulement permanent peuvent donc être utilisées dans le cadre du procédé conforme à la présente invention, pour obtenir des vitesses instationnaires d'écoulements pulsés unidimensionnelles en temps réels, à condition toutefois que le théorème de Bernouilli permanent puisse être appliqué instantanément au moment o l'accélération est nulle: ceci suppose que le dispositif utilisé ne présente pas de décollement, ce qui exclut toutes les géométries ne présentant pas un profil aérodynamique (tels que des diaphgrammes en écoulement interne). Selon les modes de réalisation précédemment décrits, la détection des instants o l'accélération du fluide est nulle est opérée par détection de l'égalité des indications mesurées par deux prises de pression statique On peut également opérer la détection des instants o l'accéléra- tion du fluide est nulle par calcul de la vitesse en cours. PREMIER MODE DE REALISATION PARTICULIER Les inventeurs ont réalisé un premier système de mesure comprenant un obstacle en forme d'antenne de Pitot équipé de trois prises de pression, tel qu'illustré sur la figure 3 annexée Les caractéristiques géométriques de ce système de mesure respectent la Norme ISO 3966-1977 (F) des tubes de Pitot Le respect de la norme a été choisi à titre d'exemple. Plus précisément, on aperçoit sur la figure 3, un obstacle 50 similaire à celui de la figure 1 comprenant une première partie cylindrique 52 de section constante s'étendant parallèlement à l'écoulement instation- naire et prolongé à son extrémité arrière 53 par une seconde partie 54 transversale à l'écoulement instationnaire Le nez 56 de la première partie cylindrique 52 est effilé, selon une enveloppe convexe régulière Les parties cylindriques 52 et 54 possèdent de préférence le même diamètre d. La zone de transition 55 qui relie les deux parties 53 et 54 a la forme générale d'un quart de tore possédant un rayon moyen égal à 3 d. Une première prise de pression statique p 1 est prévue à distance du nez 56, sur la partie cylindrique 52 La première prise de pression 51 p 1 est prévue typiquement à une distance 8 d du nez 56. Une seconde prise de pression statique 57 P 2 est prévue en aval de la première, sur la partie cylindrique 52 La seconde prise de pression 57 est placée typiquement à une distance de la première prise 51 égale à 100 mm La seconde prise 57 est placée de préférence à une distance 8 d de l'axe de la seconde partie 54. Une troisième prise de pression 58 est ménagée axialement au sommet du nez 56 pour être sensible à la pression d'arrêt Un premier capteur de pression différentielle cl est placé entre les prises de pression 51 et 58 Il fournit un signal représentatif de la différence de pression Parrêt P 1 Un second capteur de pression différentielle c 2 est placé entre les prises de pression 51 et 57 Il fournit un signal représentatif de la différence de pression p 1 p 2 Le cas échéant ces pressions différentielles pourraient être obtenues à partir de trois capteurs de pression absolue, placés respectivement au niveau des trois prises de pression indiquées précédemment. La prise de pression 51 pl, reliée à la fois aux capteurs cl et c 2, a été doublée pour que la longueur de la liaison pneumatique allant à cl n'ait pas d'effet sur la mesure faite par c 2 et réciproquement Les capteurs de pression cl et c 2 sont reliés par l'intermédiaire de préampli- ficateurs à un micro-ordinateur ou à une petite interface électronique assurant les fonctions suivantes: Acquisition des pressions (p arrêt P 1) et (p 1 -P 2) Calcul de la vitesse de débit; Sortie d'un signal analogique proportionnel à la vitesse. Le système ainsi décrit et représenté sur la figure 3 exploite aux instants o (p 1 P 2) =O (c'est-à-dire: à chaque extremum de vitesse), le théorème de Bernouilli (loi des écoulements permantes) ( 10) v 2 soit ( 8) / Pwret Pl Pl+' + =Parret V V 2 Cela permet à la fois d'engager le processus de calcul de la vitesse et de le recaler à chaque extremum de vitesse. Après discrétisation, le schéma de calcul devient ( 11) vnl+ =vn+ t (Pl -P 2) pl quand po Pa quand (pl p 2) = o N a:( 8) v =ô 2 prê Ce schéma de calcul a été choisi à titre d'exemple Une autre mode d'intégration numérique peut être choisie pour augmenter la précision du calcul. L'organigramme du processus de calcul est représenté sur la figure 4. Après une phase de démarrage 60, le micro-ordinateur ou l'interface électronique procède à l'étape 61 à l'acquisition de pressions différentielles p 1 P 2 et Parrêt -P L'étape 61 est suivie d'une étape de test 62 correspondant à l'étape i) précitée au cours de laquelle le micro-ordinateur ou l'interface électronique détecte si les pressions p 1 et P 2 sont égales. Dans l'affirmative, l'étape de test 62 est suivie d'une étape d'initiation du calcul de la vitesse et de recalage de celle-ci à chaque extremum sur la base de la relation ( 8) L'étape 63 correspond à l'étape ii) précitée. Lorsque l'étape de test 62 est négative, elle est suivie de l'étape d'intégration 64, sur la base de l'équation ( 8) Cette étape 64 correspond à l'étape iii) précitée La vitesse est exprimée au cours de l'étape 65. La rapidité des systèmes informatiques modernes (micro-ordi- nateur type PC) permet d'effectuer ce calcul en temps réel Pendant un temps At ( 1/5000 sec), le système réalise l'acquisition des pressions, le calcul de la vitesse correspondante et délivre un signal analogique proportionnel à la vitesse d'écoulement. Le système réalisé conformément à la figure 3 a permis de mesurer la vitesse de débit d'un écoulement instationnaire formé par l'écoulement pulsé à l'admission d'un moteur automobile Cette mesure, comparée à celle d'un anémomètre à fil chaud, est tout à fait satisfaisante, et a été testée pour des fréquences de pulsations allant jusqu'à 80 Hz Le taux de modulation très important de l'écoulement ne paraît pas être une limitation. Les résultats obtenus sont donnés sur la figure 5, sur laquelle on a représenté respectivement sous la référence T Pl (pour Tube de Pitot Instationnaire) la courbe de mesure obtenue à l'aide du système conforme à la présente invention, et sous la référence AFC la courbe de mesure obtenue à l'aide d'un procédé classique d'anénomètrie à fil chaud. L'erreur moyenne de la vitesse sur la courbe présentée est de: 1,75 %. La sous-estimation moyenne du débit sur la courbe présentée est de: 0,5 %. DEUXIEME MODE DE REALISATION PARTICULIER Les inventeurs ont réalisé un deuxième système dans lequel on ne dispose plus de deux prises de pression statique placées sur une paroi plane parallèle à une vitesse unique. D'une façon générale, cette circonstance peut se rencontrer soit lorsqu'on ne dispose pas de partie cylindrique o l'écoulement est parallèle, soit lorsque celle-ci est trop courte pour donner une différence de pression Pl P 2 pouvant être mesurée avec une précision suffisante. Dans ce cas, la détection d'accélération nulle est faite grâce au calcul de la vitesse et non plus à l'aide d'une comparaison de deux pressions statiques. Les inventeurs proposent de plus d'utiliser pour ce second mode de réalisation particulier le théorème de l'énergie cinétique r V 2 2 ( 12) (p) adv p+pgh+p 2-)u N da=P -P D E o Pvr et Pv D sont respectivement la puissance fournie par viscosité sur E et la puissance dissipée par frottement dans D que l'on peut supposer nulles (voir figure 6). On peut considérer, le cas d'un écoulement entrant dans une conduite profilée 70 et un domaine D, comme représenté sur la figure 6, dont la partie E O de frontière extérieure à la conduite est dans le fluide au repos et à une pression p 0, égale à la pression atmosphérique. En désigant par p la pression statique mesurée à l'aide d'une prise de pression statique 73 sur la paroi 71 parallèle à l'écoulement, de la conduite 70, et en supposant l'entrée 72 de la conduite suffisamment bien profilée pour qu'il y ait peu de pertes de charges, on obtient l'équation: *( 13) av v 2 avec 2 pl -a = (p p 0) + ( 14) "=fv 2 +dv+: coefficient caractéristique de la répartition des D vitesses calculé une fois pour toutes comme indiqué par la suite V est une vitesse sans dimension. L'équation ( 13), comme l'équation ( 7), nécessite de fournir une valeur initiale à la vitesse et son intégration numérique est susceptible de diverger dans le temps en raison des erreurs diverses (arrondis, erreur de mesures,) Ces problèmes sont résolus comme précédemment en plaçant une prise de pression totale 74 dans l'écoulement et en appliquant le théorème de Bernouilli au moment o l'accélération s'annule. La résolution de l'équation ( 13) demande, en plus, la connaissance du coefficient c, lié à la géométrie de l'entrée de la conduite, et dont la détermination s'effectue une fois pour toutes Celle-ci peut être effectuée de différentes manières, par exemple par étalonnage comparatif; en évaluant directement la vitesse de débit d'un écoulement périodique de période T. Si la différence v(n+k T) v(n) est nulle, la valeur de cl estimée est bonne Dans le cas contraire on incrémente ou décrémente a pour converger vers la bonne valeur (assez proche de 1). Plus précisément, on aperçoit sur la figure 7, une conduite profilée 70 dont l'entrée 72 diverge Une prise de pression statique cl 73 placée à une distance 1 de l'embouchure de la conduite sur la paroi 71 de celle- ci mesure la pression statique p. Une seconde prise de pression 74 c 2 placée au centre de la conduite mesure la pression d'arrêt ou pression totale p arrêt' La structure de l'algorithme du procédé est la même que dans le cas de la première réalisation précitée. Le procédé de mesure de vitesse comprend de ce fait, pour l'essentiel, les étapes consistant à: déterminer les évolutions de la vitesse entre les instants d'accélération nulle par intégration de l'équation ( 15) n+ 1 N At p(vn)2 V =V apî ' 2 détecter les instants o l'accélération du fluide est nulle par surveillance du résultat du calcul précédent, c'est-à-dire par surveillance des instants ou ( 16) vn+ 1 =vn recalage de la valeur de la vitesse quand ( 16) vn+ 1 vn à l'aide de l'équation ( 17) v MJ 2 Pq M Pl Une étude de faisabilité, réalisée sur un tube circulaire 70 comme représenté sur la figure 7 à l'entrée profilée 72 de diamètre 56 mm, de longueur l= 70 mm (longueur entre l'atmosphère libre et le capteur de pression 73) à l'amont du système d'admission d'un moteur automobile, a montré que la méthode est capable de déterminer correctement la vitesse dans le tube. VARIANTE DE REALISATION UTILISANT LE THEOREME DES DEBITS DE QUANTITE DE MOUVENEMT On va maintenant décrire une variante de réalisation utilisant le théorème des débits de quantité de mouvement. Considérons, comme cela est schématisé sur la figure 8 annexée, un écoulement en conduite cylindrique de diamètre d (section S), comportant deux prises pl, p 2 de pression statique distantes de la longueur 1 Le domaine D est un volume de surface fermée ô, M est limitée en longueur par les deux sections S situées à l'endroit des prises de pression pi et p 2. En appliquant les approximations classiques des écoulements en conduite et en appliquant le théorème des débits de quantité de mouvement au cas considéré (en projection sur l'axe), nous obtenons la formulation suivante: pav d + Vdq xt + (pl p 2)S dv E avec 9 m = débit massique, o les forces de viscosité représentées par la somme des forces extérieures sont souvent négligeables. En régime permanent, l'équation devient U 2 (p V)d V+ 5 (Pl -O) S 1 ph)S o P et -2 sont des coefficients caractérisant le profil des vitesses, et U représente la vitesse sur l'axe de la conduite. En intégrant le premier terme sur la surface S puis sur la longueur 1, nous obtenons l'équation ci-dessous p 1 Sd U + & p rdq+ t (AI 02) Ss m(PlP 2) S o la variation de débit d U q/dt est indépendante de la répartition des vitesses Si le profil des vitesses est établi, alors,I 3 est égale à 2 Nous pouvons connaître ainsi le débit en fonction de la pression différentielle (p 1 p 2). Le procédé de mesure peut alors être ainsi mesuré i) les instants o l'accélération du fluide est nulle, c'est-à-dire les instants o les pressions détectées par les capteurs pi et p 2 sont égales, sont détectés. ii) aux instants o l'accélération est nulle, la vitesse du fluide est définie sur la base du théorème de Bernoulli comme indiqué précédemment. iii) entre les instants o l'accélération est nulle, les évolutions de la vitesse sont déterminées par intégration numérique de l'équation d U p 1 Stq (pl p 2) S Bien entendu la présente invention n'est pas limitée aux modes de réalisation particuliers qui viennent d'être décrits mais s'étend à toutes variantes conformes à son esprit. R E V E N D I C A T I O N S 1 Procédé de mesure de la vitesse d'écoulement instationnaire d'un fluide, caractérisé par le fait qu'il comprend les étapes consistant à i) détecter les instants o l'accélération du fluide est nulle, ii) définir une mesure de la vitesse du fluide aux instants o l'accélération est nulle, et iii) déterminer les évolutions de la vitesse entre les instants d'accélération nulle par intégration numérique de l'écart entre deux pressions statiques. 2 Procédé de mesure selon la revendication 1, caractérisé par le fait que l'étape i) de détection des instants o l'accélération du fluide est nulle est opérée par calcul de la vitesse en cours. 3 Procédé de mesure selon la revendication 1, caractérisé par le fait que l'étape i) de détection des instants o l'accélération du fluide est nulle et opérée par détection de l'égalité des indications fournies par deux prises de pression statique. 4 Procédé de mesure selon l'une des revendications 1 à 3, caractérisé par le fait que l'étape ii) consiste à définir la vitesse aux instants o l'accélération est nulle sur la base de la relation: ( 8) V= 2 (parrêt P) p avec V module de la vitesse p 1 pression statique P masse volumique locale du fluide. 5 Procédé de mesure selon l'une des revendications 1 à 3, caractérisé par le fait que l'étape ii) consiste à définir la vitesse aux instants o l'accélération est nulle sur la base de la relation ( 17) V avec V = module de la vitesse P= masse volumique locale du fluide. 6 Procédé de mesure selon l'une des revendications 1 à 5, caractérisé par le fait que l'étape iii) consiste à intégrer l'équation ( 6) d V O V O Pl P 2 dt 2 p avec a et b des constantes du système utilisé V O = module de la vitesse Pl et P 2 = pressions mesurées au même instant en deux points du système, et p masse volumique locale du fluide. 7 Procédé de mesure selon l'une des revendications 1 à 5, caractérisé par le fait que l'étape iii) consiste à intégrer l'équation ( 7) d V O Pl_-___ dt Pl avec V O = module de la vitesse Pl et P 2 = pressions mesurées au même instant en deux points du système, 1 représente la distance entre les deux points de mesure, et p = masse volumique locale du fluide. 8 Procédé de mesure selon l'une des revendications 1 à 5, caractérisé par le fait que l'étape iii) consiste à intégrer l'équation ( 13) av 2 apl =(P Po) + p V avec a, représentant un coefficient d'énergie cinétique, p représente la masse volumique locale du fluide, 1 représente la distance entre l'embouchure de la conduite et une prise de pression statique mesurant la pression p, p O représente la pression atmosphérique. 9 Dispositif de mesure de la vitesse d'écoulement instation- naire d'un fluide, pour la mise en oeuvre du procédé conforme à l'une des revendications 1 à 8, caractérisé par le fait qu'il comprend: des moyens aptes à détecter les instants o l'accélération du fluide est nulle, des moyens aptes à définir une mesure de la vitesse du fluide aux instants ou l'accélération est nulle, et des moyens aptes à déterminer les évolutions de la vitesse entre les instants d'accélération nulle par l'intégration numérique de l'écart entre deux pressions statiques. Dispositif de mesure selon la revendication 9, caractérisé par le fait que les moyens aptes à détecter les instants o l'accélération du fluide est nulle comprennent deux prises de pression statique ( 14, 16) et des moyens aptes à surveiller l'égalité des deux pressions mesurées. Il Dispositif de mesure selon la revendication 9, caractérisé par le fait que les moyens aptes à détecter les instants o l'accélération du fluide est nulle comprennent des moyens de surveillance de la vitesse calculée. 12 Dispositif de mesure selon l'une des revendications 9 ou , caractérisé par le fait qu'il comprend un obstacle du type tube de Pitot équipé de trois prises de pression. 13 Dispositif de mesure selon la revendication 12, caractérisé par le fait qu'il comprend deux prises de pression statique ( 51, 57) espacées sur une partie ( 52) de l'obstacle parallèle à l'écoulement, une prise de pression d'arrêt au nez de l'obstacle et deux capteurs de pression différentielle (cl, c 2), ou trois capteurs de pression absolue. 14 Dispositif de mesure selon la revendication 13, caractérisé par le fait que les moyens aptes à définir la mesure de la vitesse du fluide aux instants o l'accélération est nulle exploitent l'équation ( 8) / j 2 et, les moyens aptes à déterminer les évolutions de la vitesse exploitent l'équation ( 7) d Vo +Pl -P 2 P- 15 Dispositif de mesure selon l'une des revendications 9 ou , caractérisé par le fait qu'il comprend une conduite munie d'un rétrécissement et trois prises de pression statique. 16 Dispositif de mesure selon la revendication 15, caractérisé par le fait que les moyens de détection des instants d'accélération nulle comprennent des moyens de détection de l'égalité entre deux pressions statiques mesurées à distance sur une partie de la conduite de section constante, les moyens définissant la mesure de la vitesse du fluide aux instants d'accélération nulle exploitent la différence de pression entre l'une des pressions précitées et la troisième pression mesurée sur une zone de section différentielle de la conduite, et les moyens déterminant les évolutions de la vitesse intègrent la différence entre les deux pressions premières citées. 17 Dispositif de mesure selon la revendication 16, caractérisé par le fait que les deux pressions statiques premières citées sont prélevées sur une partie de section constante de la conduite tandis que la troisième pression est prélevée sur l'embouchure évasée de la conduite. 18 Dispositif de mesure selon la revendication 16, caractérisé par le fait que les deux pressions statiques premières citées sont prélevées à distance sur une partie de section constante de la conduite tandis que la troisième pression statique est mesurée au niveau d'un rétrécissement. 19 Dispositif de mesure selon la revendication 16, caractérisé par le fait que les deux pressions statiques premières citées sont mesurées au niveau d'un rétrécissement de section constante tandis que la troisième pression statique est mesurée à l'extérieur du rétrécissement. Dispositif de mesure selon l'une des revendications 9 ou 19, caractérisé par le fait qu'il comprend une conduite profilée ( 70), un capteur de pression statique ( 73) sur la paroi de la conduite et un capteur ( 74) de pression totale. 21 Dispositif de mesure selon la revendication 20, caractérisé par le fait que les moyens aptes à détecter les instants d'accélération nulle surveillent la vitesse calculée, les moyens définissant la vitesse aux instants o l'accélération est nulle exploitent l'équation ( 17) P p les moyens déterminant les évolutions de la vitesse entre les instants d'accélération nulle exploitent l'équation ( 15) vn' =vn At (p _ p(vo 2 apl 2 avec a représente un coefficient d'énergie cinétique, p représente la masse volumique du fluide, 1 représente la distance entre la prise de pression statique ( 73) et l'ebouchure de la conduite, p représente la pression statique, p O représente la pression atmosphérique, et Parrêt représente la pression totale.