Brevet FR2929399A1 — publié le 20091002

Débitmètre instationnaire temps réel DOMAINE DE L'INVENTION La présente invention concerne le domaine de la mesure du débit instantané d'un fluide en écoulement instationnaire. Elle trouve notamment, mais non limitativement, application dans les domaines du génie des procédés, et de l'industrie automobile. Plus particulièrement, elle trouve avantageusement application dans le contrôle et la régulation des moteurs thermiques, les centres d'essais et laboratoires de recherche, ainsi que toutes les applications actuelles des débitmètres de l'art antérieur. Plus particulièrement, elle trouve avantageusement application dans toutes les situations où l'écoulement étant instationnaire, la connaissance du débit impose une intégration dans le temps et dans l'espace de la vitesse. ETAT DE LA TECHNIQUE Les mesures de débits sont des éléments indispensables à la mise en oeuvre et à l'optimisation des processus industriels. A l'heure actuelle, pour connaître le débit d'un fluide en écoulement instationnaire à un instant donné, il est généralement nécessaire de mettre en oeuvre des techniques de mesure permettant de déterminer la vitesse de ce fluide. Une fois que l'on a déterminé un ensemble de vitesses pour différents instants, on intègre spatialement cet ensemble de vitesses afin d'avoir accès à l'évolution du débit en fonction du temps. Il a déjà été proposé un certain nombre de systèmes et méthodes permettant l'obtention du débit d'un fluide en écoulement instationnaire par la détermination de la vitesse de ce fluide. A titre d'exemple, on peut citer les méthodes dites de Vélocimétrie par Images de Particules (PIV) ou les méthodes dites de Vélocimétrie Doppler Laser (LDV), ainsi que les systèmes qui leurs sont associés. Ces méthodes sont basées sur la mesure du déplacement de petites particules contenue dans le fluide en écoulement, ce déplacement étant mesuré par une technique de traitement d'image d'une part, ou par traitement de fréquences Doppler d'autre part. Toutefois ces méthodes et systèmes présentent un certain nombre d'inconvénients. En particulier, les méthodes PIV et LDV ne permettent pas de déterminer le débit d'un fluide en temps réel, puisque ces techniques imposent un post-traitement des images acquises afin d'obtenir la vitesse du fluide. De plus, la mise en oeuvre de telles techniques est lourde et coûteuse. Ces méthodes sont également mal résolues en temps et possèdent donc une faible bande passante (généralement moins de 10 Hz pour la méthode PIV). Une autre méthode permettant d'obtenir le débit d'un fluide en écoulement instationnaire à partir de sa vitesse utilise le tube de Pitot (ou sonde de Prandtl). Le tube de Pitot est un instrument destiné à être plongé dans une conduite pour déterminer la vitesse d'un fluide en écoulement grâce à la mesure d'une différence de pression. Toutefois, la direction et le sens de l'écoulement doivent être connus et constants pour mettre en place le tube de Pitot. En outre, lorsque le sens de l'écoulement varie, la mesure de pression totale est impossible, ce qui induit une erreur dans l'obtention de la mesure de la vitesse. Enfin, ce type de système est intrusif et perturbe donc l'écoulement. Le débit d'un fluide en écoulement instationnaire peut également être obtenu à partir de sa vitesse par l'utilisation d'un système à fil chaud ou à film chaud placé dans le fluide en écoulement. Néanmoins, une telle méthode est ponctuelle et insensible au sens de l'écoulement du fluide. En outre, sa mise en oeuvre est complexe et coûteuse, en particulier parce que le fil ou film chauffant est très fragile et vieillit rapidement ce qui nécessite une maintenance régulière. Enfin, c'est une méthode intrusive qui modifie donc l'écoulement. Pour remédier aux inconvénients des méthodes présentées ci-dessus, il a été développé un système de mesure particulier permettant de mesurer en temps réel le débit instantané d'un fluide en écoulement et prenant en compte le sens de l'écoulement du fluide. Ce système et la méthode de mesure associée sont plus particulièrement décrits dans la demande PCT/FR2005/000352 publiée le 1 er septembre 2005 sous la référence WO2005/080924, et que l'on incorpore ici par référence. Un tel système est particulièrement performant pour mesurer le débit instantané de fluide dont l'écoulement varie à une fréquence élevée de plusieurs dizaines de Hertz. Il présente toutefois quelques inconvénients liés notamment à l'encombrement du système proposé.35 Un but de la présente invention est donc de fournir un système de mesure en temps réel du débit d'un fluide en écoulement instationnaire, permettant de résoudre la plupart des inconvénients précités. EXPOSE DE L'INVENTION A cette fin, on propose un système de mesure en temps réel du débit instantané d'un fluide en écoulement stationnaire ou instationnaire dans une conduite, comprenant un organe d'écoulement du fluide muni d'au moins deux prises de pression en paroi, un moyen de mesure d'une différence de pression couplé aux deux prises de pression, et un moyen de calcul pour calculer le débit en temps réel par la résolution d'une équation différentielle ordinaire non-linéaire reliant le débit instantané à la différence de pression, la différence de pression dans ladite formule étant positive ou négative en fonction de la variation de la vitesse d'écoulement du fluide dans la conduite et/ou du sens de l'écoulement du fluide, caractérisé en ce que l'organe d'écoulement comprend un filtre placé entre les deux prises de pression pour augmenter la perte de charge. Des aspects préférés mais non limitatifs de ce système de mesure sont les suivants : - l'organe d'écoulement a une géométrie cylindrique avec une section circulaire de diamètre constant ; - l'organe d'écoulement a une géométrie tronconique de type Venturi ; - l'organe d'écoulement comprend deux filtres supplémentaires agencés respectivement en amont et en aval de l'organe d'écoulement par rapport à l'écoulement du fluide de manière à conditionner l'écoulement ; - le ou les filtres sont agencés sensiblement perpendiculaires à l'axe de l'organe d'écoulement ; - le ou les filtres sont des grilles ; - le ou les filtres ont une structure en nid d'abeille ; - le ou les filtres sont formés dans un matériau poreux ; - l'organe d'écoulement comprend en outre une ou plusieurs sondes de mesure de température couplée au moyen de calcul ; - l'organe d'écoulement comprend en outre une sonde de mesure de pression statique couplée au moyen de calcul.35 DESCRIPTION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des dessins annexés, sur lesquels : - la figure 1 est une représentation schématique d'un système de mesure du débit instantané d'un fluide selon l'invention ; - la figure 2 est un graphique permettant de comparer le débit mesuré par un système selon l'invention au débit mesuré par un système à fil chaud. DESCRIPTION DETAILLEE DE L'INVENTION La figure 1 représente schématiquement le système proposé, capable de mesurer en temps réel le débit d'un fluide dans une conduite à partir d'une différence de pression. Ce système est un débitmètre et se compose essentiellement d'un organe d'écoulement 1 destiné à s'insérer dans la conduite, d'un moyen 2 de mesure d'une différence de pression, et d'un moyen 3 de calcul adapté pour calculer le débit en temps réel à partir d'une différence de pression mesurée dans l'organe d'écoulement. L'organe d'écoulement 1 a une forme particulièrement simple adaptée pour être facilement inséré, ou interposé, dans la conduite d'écoulement du fluide dont on veut mesurer le débit instantané. Plus particulièrement, l'organe d'écoulement 1 comprend une section intérieure S quelconque, l'organe pouvant avoir une géométrie cylindrique. De manière préférée, la forme générale de l'organe d'écoulement 1 est cylindrique avec une section circulaire de diamètre constant, avec les mêmes caractéristiques que la conduite où il s'insère. Il peut également avoir une forme tronconique, formant ainsi un Venturi. Deux prises de pression statique A et B sont placées sur la paroi de cet organe d'écoulement 1. Ces prises de pression statique A et B sont couplées au moyen de mesure 2 d'une différence de pression du système. Le moyen 2 permet donc de mesurer la différence entre les deux pressions statiques P1 et P2 prises respectivement en A et B. On utilise par exemple un capteur de pression différentielle. Le moyen de calcul 3 est un calculateur électronique. Ce moyen de calcul 3 est programmé pour mettre en oeuvre un algorithme de résolution d'une formule que l'on explicitera plus loin, qui est une équation différentielle ordinaire non- linéaire permettant le calcul du débit instantané, le signe de ce débit rendant compte du sens de l'écoulement. Le calculateur électronique 3 est couplé au capteur 2 de pression différentielle puisque ce calculateur 3 calcule le débit instantané à partir de la différence de pression mesurée dans l'organe d'écoulement 1 qui est elle aussi instantanée. Comme on l'a indiqué plus haut, l'organe d'écoulement 1 présente comme première caractéristique d'avoir une forme cylindrique particulièrement simple à mettre en oeuvre. L'organe d'écoulement peut par exemple avoir une section circulaire avec un diamètre constant, ce qui simplifie d'autant sa mise en oeuvre. En outre, cet organe d'écoulement 1 intègre un filtre 4 disposé dans le passage du fluide à l'intérieur de l'organe d'écoulement 1, et placé entre les deux prises de pression A et B. Ce filtre 4 est prévu pour induire une perte de charge dans l'organe d'écoulement 1, plus particulièrement entre la prise de pression A et la prise de pression B, de manière à créer une différence de pression entre ces deux prises de pressions A et B. Ce filtre 4 est donc prévu pour introduire dans l'organe d'écoulement 1 une perte de charge supplémentaire par rapport à la perte de charge régulière due à la géométrie, et plus particulièrement aux parois, de l'organe d'écoulement 1. Comme il est décrit en détail plus loin, il existe une relation entre la pression et le débit instantané pour un fluide en écoulement instationnaire dans une conduite, ou plus précisément une relation entre d'une part la pression, et d'autre part l'accélération du fluide et son énergie cinétique (qui est directement liée aux pertes de charge). La perte de charge supplémentaire introduite par le filtre 4 permet d'équilibrer l'importance de l'énergie cinétique de l'écoulement par rapport à l'accélération du fluide, de manière à pouvoir calculer plus précisément le débit instantané du fluide dans la conduite à partir de la différence de pression mesurée. Le filtre 4 permet donc d'adapter la perte de charge à la dynamique particulière de l'écoulement du fluide. Le filtre 4 peut avoir un agencement, une forme, et une constitution 30 quelconques, tant qu'il permet d'introduire une perte de charge homogène entre les deux prises de pression A et B. De préférence, on utilise un filtre 4 dont la section correspond sensiblement à la section intérieure S de l'organe d'écoulement 1, de manière à couvrir la quasi-totalité de cette section intérieure S. En outre, on dispose de préférence le filtre 4 35 de façon sensiblement perpendiculaire à l'écoulement du fluide dans l'organe d'écoulement 1, cet écoulement du fluide étant représenté schématiquement dans la figure 1 par les flèches noires. Le filtre 4 peut par exemple prendre la forme d'une grille. Ce filtre 4 peut avoir une structure en forme de nid d'abeilles. On peut également prévoir de former ce filtre 4 dans un matériau poreux. Le système de mesure de débit présenté est donc particulièrement simple à concevoir et à mettre en oeuvre. Ceci est en particulier le cas lorsque l'organe d'écoulement 1 consiste en une portion cylindrique avec une section circulaire de diamètre constant, formant ainsi un tube, dont les dimensions correspondent sensiblement aux dimensions de la conduite à l'intérieur de laquelle circule le fluide à mesurer, et qu'un simple filtre 4 est inséré à l'intérieur de cette portion cylindrique. En outre, une telle configuration présente d'autres avantages notables par rapport aux dispositifs de l'art antérieur. En particulier, ce système a une configuration symétrique ce qui permet en premier lieu de simplifier son insertion dans la conduite d'écoulement du fluide. D'autre part, cela permet d'avoir un fonctionnement et une fiabilité indépendants du sens de l'écoulement du fluide dans la conduite. En effet, cette symétrie dans la géométrie de l'organe d'écoulement 1 implique que le système de mesure est très peu sensible aux conditions amont et aval puisque le filtre 4 a exactement les mêmes effets en termes de perte de charge quelque soit le sens considéré d'écoulement du fluide. La précision du système de mesure est donc indépendante du sens de l'écoulement, ce qui est particulièrement avantageux dans certains domaines d'application comme celui de l'automobile, pour lesquels certains écoulements, comme par exemple les écoulements à l'aspiration des moteurs thermiques, peuvent changer périodiquement de sens. Un autre avantage induit par cette configuration particulière de l'organe d'écoulement 1 réside dans le faible encombrement du système de mesure correspondant. En effet, la longueur de l'organe d'écoulement 1 peut être très sensiblement réduite, par rapport par exemple à des systèmes reposant sur des profils de vitesses particuliers (de type venturi par exemple) qui imposent de concevoir des organes d'écoulement avec une géométrie ad hoc, avec des conduites, à l'amont et à l'aval de la zone de mesure, de longueurs suffisantes pour conditionner l'écoulement.. L'utilisation d'un filtre 4 dans l'organe d'écoulement 1 présente en outre l'avantage de réduire les perturbations du fluide en écoulement, ce qui fiabilise les mesures de pressions effectuées. L'organe d'écoulement 1 peut en outre comprendre des filtres supplémentaires 5 et 6 disposés respectivement en amont et en aval de l'organe d'écoulement 1 par rapport au sens principal d'écoulement du fluide. De tels filtres 5 et 6 peuvent avoir les mêmes formes, structures et agencements que le filtre 4 disposé entre les prises de pression A et B. Ces filtres amont 5 et aval 6 permettent de s'affranchir des éventuelles structures de grande échelle de l'écoulement du fluide. En effet, si l'écoulement du fluide contient des structures tourbillonnaires de tailles importantes, alors les filtres 5 et 6 placés aux deux extrémités de l'organe d'écoulement 1 permettent de conditionner l'écoulement du fluide, c'est à dire de stabiliser et régulariser l'écoulement du fluide à l'intérieur de l'organe d'écoulement 1 pour que les mesures de pression soient encore plus fiables. Cette solution permet donc d'effectuer des mesures fiables de débit instantané de l'écoulement du fluide quelque soit la structure générale d'écoulement dudit fluide, tout en conservant un très faible encombrement. En effet, la solution proposée n'impose pas d'augmenter la longueur de l'organe d'écoulement 1 pour régulariser l'écoulement. Le principe de fonctionnement du système de mesure de la figure 1 est le suivant. Le fluide s'écoule dans une conduite et passe dans l'organe d'écoulement 1. Le moyen 2 de mesure d'une différence de pression mesure la différence entre la pression statique P1 acquise au niveau de la prise de pression A, et la pression statique P2 acquise au niveau de la prise de pression B. La différence de pression mesurée par le moyen 2 de mesure est transmise en entrée du moyen 3 de calcul. Le moyen 3 de calcul calcule le débit du fluide en temps réel à partir de la différence de pression reçue en entrée. En effet, comme on l'a dit, le moyen 3 de calcul est adapté pour résoudre la formule reliant le débit à la différence de pression, cette formule permettant le calcul du débit en temps réel et donnant en outre le sens de l'écoulement du fluide. A la sortie du moyen 3 de calcul, on obtient en temps réel le débit du fluide, même en écoulement instationnaire. Le dispositif de la figure 1 permet donc de mesurer en temps réel le débit instationnaire d'un fluide dans une conduite. La mesure de ce débit ne nécessite que deux mesures de pression statique, la différence de pression correspondante étant par exemple obtenue par un unique capteur de pression différentielle 2 couplé aux prises de pression, ou par deux capteurs de pression relative (ou absolue). Comme il a été indiqué plus haut, le calcul du débit en temps réel de l'écoulement du fluide est réalisé à partir d'une relation établie entre le débit instantané et la pression pour un fluide en écoulement instationnaire dans une conduite. Rappelons que la valeur de ce débit peut être positive ou négative, de sorte que le signe du débit calculé indique le sens de l'écoulement du fluide. Une telle relation avait déjà été établie pour le débitmètre instationnaire temps réel décrit dans la demande PCT publiée sous la référence WO 2005/080924, à laquelle on se référera utilement pour une explication détaillée. La formule utilisée dans le cadre du présent système de mesure de débit est donnée par une équation différentielle ordinaire non-linéaire de la forme générale suivante : K1.dq(t)/dt + K2.f(q(t)) = ÈP(t) [Equation 1] où: - q(t) représente la forme générale du débit instantané recherché, - KI et K2 sont des constantes, - dq(t)/dt représente la dérivée par rapport au temps du débit recherché, - f(q(t)) est une fonction non linéaire du débit volumique instantané, - ÈP(t) est la différence de pression mesurée. En conséquence, si l'on note qm(t) le débit massique instantané recherché, et qv(t) le débit volumique instantané recherché, alors l'équation 1 peut s'écrire : soit ki.dqm(t)/dt + k2.f(gm(t))/rho = ÈP(t) soit ki.dgv(t)/dt + k2.f(gv(t)) = AP(t)/rho où : - kl et k2 sont des constantes déterminées par étalonnage. - rho est la masse volumique du fluide, avec qm(t)=rho. qv(t) Pour f(q(t)) qui est une fonction non linéaire du débit volumique instantané, 35 on pourra par exemple prendre : f(q(t))=Iq(t)I.q(t) [Equation 4]) [Equation 2] [Equation 3] La présence dans les formules de termes dont le signe varie en fonction de la variation de la vitesse d'écoulement du fluide dans la conduite et/ou du sens de l'écoulement (termes dq(t)/dt et ÈP(t)), et de terme dont le signe varie en fonction du sens de l'écoulement du fluide dans la conduite (terme f(q(t))) permet de calculer un débit positif ou négatif, dont le signe rend compte du sens de l'écoulement du fluide. Notamment, dans les formules décrites ci-dessus, le terme ÈP(t) est pris en compte tel quel, sans valeur absolue. Le signe attribué à la différence de pression ÈP(t) varie notamment en fonction de la variation de vitesse d'écoulement du fluide et/ou du sens de l'écoulement du fluide. La différence de pression ÈP(t) est donc positive ou négative, par exemple lorsque le fluide s'écoulant dans la conduite accélère ou décélère, et/ou lorsque le fluide s'écoule dans la conduite dans un sens ou dans le sens opposé. Par ailleurs, dans les formules décrites ci-dessus, les termes dq(t)/dt et f(q(t)) sont également pris en compte tels quels, sans valeur absolue. Le signe du débit calculé permet de rendre compte du sens de l'écoulement, et les inversions de signe du débit calculé permettent de rendre compte des inversions de sens d'écoulement du fluide. Les équations différentielles ordinaires non-linéaires formulées possèdent la propriété remarquable de toujours converger vers une solution bornée pour peu que la condition initiale soit convenablement choisie de l'ordre de grandeur du débit à trouver. Préférentiellement, on choisira comme condition initiale q(t=0)=0. En pratique, lorsque la condition initiale est correctement choisie, le signal de sortie du calculateur 3 converge vers la valeur instantanée du débit recherché. L'algorithme de résolution peut être mis en oeuvre aussi bien de façon numérique qu'analogique, et permet d'obtenir la valeur du débit instantané en temps réel. Le débit mesuré grâce à cette méthode correspond au débit q(t) d'un fluide incompressible en écoulement instationnaire dans une conduite. La masse volumique du fluide est une donnée qui doit être connue pour résoudre les équations 1, 2 ou 3 ci-dessus. Il convient de noter qu'elle est elle-même fonction d'un certain nombre d'autres grandeurs physiques comme par exemple la température, la pression statique et éventuellement la masse molaire du fluide. Dans le cas où la compressibilité du fluide ne peut pas être négligée, on intègre au système de mesure une sonde de mesure de la température et/ou une sonde de mesure de la pression statique régnant dans l'écoulement, ce qui permet de prendre en compte les fluctuations de masse volumique dans l'écoulement de fluide considéré. S'il est possible par un moyen tierce de déterminer la valeur instantanée de la valeur absolue du débit volumique IgvI alors l'instrument proposé permet d'en déduire le débit massique instantané qrr, (mais cette fois-ci avec son sens), indépendamment des éventuelles fluctuations de la masse volumique rho, qui n'a pas besoin d'être déterminée dans ce cas. C'est ce que montre l'équation 5 tirée des équations 2 et 4 : ki.dqm(t)Idt + k2.gm(t),Igv(t)I = ÈP(t) [Equation 5] De même s'il est possible par un moyen tierce de déterminer une grandeur statistique liée à la valeur absolue du débit volumique IgvI, comme la valeur moyenne (ou la valeur efficace) du débit volumique IgvI, alors l'instrument proposé permet d'en déduire le débit massique instantané qm (avec son sens), indépendamment des éventuelles fluctuations lentes de la masse volumique rho, qui n'a donc plus besoin d'être déterminée dans ce cas. C'est ce que montre l'équation 6 tirée des équations 2 et 4 : ki.dqm(t)Idt + k2.gm(t),Igm(t)I. / = ÈP(t) [Equation 6] où <> représente l'opérateur moyenne arithmétique ; représente donc la valeur moyenne de la valeur absolue du débit pendant un certain laps de 25 temps. Le graphique de la Figure 2 met en avant les résultats de mesure de débit instantané obtenus en utilisant le système de la figure 1. Plus précisément, ce graphique représente l'évolution du débit pulsé dans une conduite d'aspiration 30 d'un moteur thermique monocylindre mesuré par le système proposé, et comparée avec le débit reconstitué à partir de la vitesse du fluide déterminée par un système à fil chaud. La première courbe Cl et la deuxième courbe C2 de la figure 2 présentent l'évolution du débit dans une conduite en fonction du temps. Les traits pleins 35 (courbe Cl) représentent les résultats obtenus avec le système proposé, tandis que les traits pointillés (courbe C2) représentent les résultats obtenus avec le système de mesure par fil chaud. On peut remarquer sur la courbe C2 que le système de mesure par fil chaud, qui est incapable par conception de donner le sens de l'écoulement, donne toujours un débit positif, même lorsque celui-ci s'inverse. Ce n'est pas le cas du système de la présente invention qui suit parfaitement les inversions de sens de l'écoulement comme on peut le voir sur la courbe Cl. Ceci est dû au fait que le système de mesure de la présente invention tient compte du sens de l'écoulement. Le système proposé possède une très bonne résolution en temps (grande bande passante) puisqu'il a été testé avec succès jusqu'à un régime de 3000 tours/minute. Sur la courbe Cl est reportée l'évolution du débit pour un régime de ralenti de 900 tours/minute. Ainsi, le système proposé a une résolution temporelle importante, comparable aux systèmes à fil ou films chauds, sans présenter leurs nombreux inconvénients. Le lecteur aura compris que de nombreuses modifications peuvent être apportées sans sortir matériellement des nouveaux enseignements et des avantages décrits ici. Par conséquent, toutes les modifications de ce type sont destinées à être incorporées à l'intérieur de la portée du système de mesure selon l'invention.20 REVENDICATIONS1. Système de mesure en temps réel du débit instantané d'un fluide en écoulement stationnaire ou instationnaire dans une conduite, comprenant un organe d'écoulement (1) du fluide muni d'au moins deux prises de pression (A, B) en paroi, un moyen de mesure (2) d'une différence de pression couplé aux deux prises de pression (A, B), et un moyen de calcul (3) pour calculer le débit en temps réel par la résolution d'une équation différentielle ordinaire non-linéaire reliant le débit instantané à la différence de pression, la différence de pression dans ladite formule étant positive ou négative en fonction de la variation de la vitesse d'écoulement du fluide dans la conduite et/ou du sens de l'écoulement du fluide, caractérisé en ce que l'organe d'écoulement (1) comprend un filtre (4) placé entre les deux prises de pression (A, B) pour augmenter la perte de charge. 2. Système selon la revendication 1, caractérisé en ce que l'organe d'écoulement (1) a une géométrie cylindrique avec une section (S) circulaire de diamètre constant. 3. Système selon la revendication 1, caractérisé en ce que l'organe d'écoulement (1) a une géométrie tronconique de type Venturi. 4. Système selon l'une quelconque des revendications 1 à 3, caractérisé en ce que l'organe d'écoulement (1) comprend deux filtres (5,6) supplémentaires agencés respectivement en amont et en aval de l'organe d'écoulement (1) par rapport à l'écoulement du fluide de manière à conditionner l'écoulement. 5. Système selon l'une quelconque des revendications 1 à 4, caractérisé en ce que le ou les filtres (4,5,6) sont agencés sensiblement perpendiculaires à l'axe de l'organe d'écoulement. 6. Système selon l'une quelconque des revendications 1 à 5, caractérisé en ce que le ou les filtres (4,5,6) sont des grilles. 7. Système selon l'une quelconque des revendications 1 à 6, caractérisé en ce que le ou les filtres (4,5,6) ont une structure en nid d'abeille. 8. Système selon l'une quelconque des revendications 1 à 7, caractérisé en ce 30 que le ou les filtres (4,5,6) sont formés dans un matériau poreux. 9. Système selon l'une quelconque des revendications 1 à 8, caractérisé en ce que l'organe d'écoulement (1) comprend en outre une ou plusieurs sondes de mesure de température couplée au moyen de calcul (3). 10. Système selon l'une quelconque des revendications 1 à 9, caractérisé en ce 35 que l'organe d'écoulement (1) comprend en outre une sonde de mesure de pression statique couplée au moyen de calcul (3).