Brevet FR2749943A1 — publié le 19971219

SONDE VELOCIMETRIQUE OPTIQUE Le domaine de l'invention est celui des sondes vélocimétriques optiques permettant de définir la vitesse de particules en mouvement relatif par rapport à la sonde. Il existe à l'heure actuelle des dispositifs optiques comprenant deux faisceaux lasers qui interfèrent de manière à générer un volume de mesure ellipsoïdal, consistant en des franges équidistantes sombres et lumineuses. Lorsqu'une particule en mouvement traverse ce volume de mesure, le signal recueilli par un photodétecteur comprend des indications relatives à la vitesse de cette particule par rapport à la sonde optique. Il est donc possible en analysant le signal recueilli de définir la vitesse d'une particule considérée. Ce type de sonde peut être utilisé pour définir la vitesse de toute particule suffisamment petite devantles franges d'interférences créées et notamment pour définir la vitesse d'un aéronef, connaissant la vitesse relative de particules dans l'air par rapport audit aéronef équipé de ce type de sonde vélocimétrique. En effet, en raison de progrès réalisés en matière de miniaturisation et de coût au niveau des sources laser et notamment des diodes laser, il devient possible de réaliser de telles sondes vélocimétriques optiques, embarquées à bord de l'aéronef dont on cherche à mesurer la vitesse. Les performances de ce type de sonde sont directement dépendantes des moyens d'analyse du signal recueilli par un photodétecteur, pour en extraire le paramètre vitesse. C'est pourquoi I'invention a pour - objet une sonde vélocimétrique optique du type de celle décrite précédemment et comprenant des moyens de détermination de vitesse très performants. Plus précisément, I'invention a pour objet une sonde vélocimétrique optique comprenant des moyens pour éclairer un volume de mesure dans lequel peuvent passer des particules en mouvement relatif par rapport à la sonde, et des moyens de détection optique pour produire un signal électrique en réponse au passage d'une particule dans le volume éclairé, le spectre en fréquence de ce signal ayant une largeur centrée sur une fréquence dominante fo, représentative de la vitesse relative de la particule par rapport à la sonde, cette sonde comportant des moyens numériques pour établir un vecteur x(t) représentant N échantillons du signal électrique détecté et étant caractérisée en ce qu'elle comporte: des moyens pour convertir le signal électrique x(t) en un signal complexe x(t); - des moyens pour établir une fonction d'autocorrélation du signal complexe, ladite fonction étant définie de la façon suivante: N-m r;(m)= Z EJ E (t)k(t+m)} t=-N o E {x} désigne l'espérance mathématique de la variable x. - des moyens pour déterminer la phase 4(m) de la fonction r(m) reliée à la partie imaginaire de ladite fonction; - des moyens pour déterminer la fréquence fo à partir de l'ensemble des couples [mn,x(m)] - des moyens pour fournir une indication de vitesse v d'une particule par rapport à la sonde, à partir de la fréquence fo. Lorsque le signal électrique fourni par le photodétecteur résulte du rayonnement lumineux rétroréfléchi par une particule traversant un volume de mesure défini par les franges d'interférences résultant de deux faisceaux lumineux, le vecteur x(t) peut être défini comme suit: x(t) = A.e-2u2f2t2 cos(2iTft) + w(t) avec ca un coefficient fixe t est le temps w(t) un terme de bruit Les moyens employés pour transformer le signal x(t) réel en signal x(t) complexe comprennent un filtre déphaseur de -n/2 du type filtre de Hilbert. Ce type de filtre permet de définir un signal complexe x(t) associé au signal réel x(t) dont la partie réelle n'est autre que le signal de départ x(t) et la partie imaginaire est la sortie du filtre de Hilbert avec x(t) en entrée - soit encore x(t)= x(t)+j T.H. [x(t)] Si T.H. est la transformée de Hilbert correspondant au résultat du filtrage du signal par le filtre de Hilbert. Selon une variante de l'invention, les premiers et secondes moyens concernant le signal complexe x(t) et la fonction d'autocorrélation Fx(m) mis en oeuvre pour la détermination de la vitesse relative, peuvent avantageusement être regroupés de la manière suivante: - détermination du spectre discret x(co), du signal calculé par transformée de Fourier rapide (FFT), sur les N échantillons; - définition du nouveau spectre x(o) par (wo) = 2x(O))j2 pour co > 0 X(c)) = O pourro définit l'ensemble des phases cfr(m). (od On peut alors estimer la fréquence fo = fd = 2 en considérant que cette fréquence correspond à la valeur pour laquelle la fonction J(0d) est minimale avec: J(CO d) = T_ x(m)-mO d comme cela sera explicité ultérieurement. Tous ces moyens utilisés dans l'invention seront détaillés plus amplement dans la description qui va suivre. L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description qui va suivre, donnée à titre non limitatif et grâce aux figures annexées par lesquelles - la figure 1 illustre un exemple de sonde vélocimétrique selon l'invention; - la figure 2 représente le volume de mesure dans lequel passe les particules dont on veut détecter la vitesse relative. De manière générale, la sonde vélocimétrique optique de I'invention comprend deux faisceaux optiques interférant dans un volume de mesure. La figure 1 illustre un exemple de réalisation de sonde vélocimétrique montée à bord d'un aéronef, comprenant une source laser 1, pouvant être de type diode laser, en sortie de laquelle on dispose un collimateur 2, un diviseur de faisceau optique 3 créant deux faisceaux optiques L1 et L2 parallèles, dirigés au moyen d'un système de renvoi 4, vers une optique 51 en sortie de l'aéronef matérialisée par la peau de l'aéronef 6. L'optique 51 a pour fonction de faire converger les faisceaux lumineux LI et L2, dans un volume dit volume de mesure MV, extérieur compris dans l'atmosphère dans laquelle circulent des particules constitutives d'aérosols atmosphériques. Ces aérosols renvoient de la lumière rétrodiffusée selon un faisceau lumineux RL, détecté par le photodétecteur 7, pouvant typiquement être une photodiode, via des optiques 51, 52, 53. En sortie de la photodiode, la sonde vélocimétrique comprend des moyens de traitement du signal 8 comprenant des moyens de conversion du signal analogique en signal numérique, des moyens d'extraction d'informations de fréquence et des moyens de calcul de la vitesse relative des particules par rapport à l'aéronef. La figure 2 schématise le volume de mesure VM, constitué par les franges d'interférence des deux faisceaux optiques L1 et L2. La largeur 2W représente la largeur totale de ce volume, I étant la largeur d'interfrange. La sonde selon l'invention comprend les moyens 8 de traitement du signal recueilli au niveau du photodétecteur 7, permettant de déterminer avec une grande précision la vitesse d'une particule traversant le volume de mesure VM. Nous allons décrire les moyens 8 en rappelant tout d'abord les considérations suivantes: Lorsqu'une particule ayant une vitesse relative v par rapport au volume de mesure VM, pénètre dans ce réseau d'interfranges, la lumière rétrodiffusée par cette particule et recueillie par le photodétecteur 5, après échantillonnage, fournit un signal s(t) dont le modèle mathématique est connu: s(t) e-2(v/w)2t2 cos(2 v / I. t) avec t = 0, 1,... + 3w/v o v représente la vitesse de la particule; a 2w la largeur totale des franges d'interférences définissant le volume de mesure comme l'illustre la figure 2; À I la largeur d'interfrange avec aE = I/w et fd = v/I s(t) = A. e-22fd2t2 cos(27fdt) Le terme de phase de ce signal est celui qui apparaît dans la modulation en cosinus de l'amplitude, soit 27cfdt. On peut alors montrer en considérant le spectre s(t) que le signal analytique associé à celui-ci s'exprime par: s(t) = s(t) + j. T. H. [s(t)] si T.H. est la transformée de Hilbert du signal temporel s(t), résultat du filtrage de ce signal par le filtre de Hilbert h(t) s(t) = s(t)+ j[h(t)* s(t)] soit encore 1(t) e-2,2fdt2 cos(27rfdt) + je(- 2'x2fd2t2) cos 27cfdt-2 -2ca2f t2 j(27r f t) =e e Dans le cas réel de la sonde vélocimétrique le signal reçu au niveau du photodétecteur comprend un terme de bruit w(t) que l'on peut considérer comme un bruit réel blanc Gaussien de moyenne nulle et de variance a2. Ainsi le signal détecté est le signal x(t) avec x(t)=A s(t) + w(t) avec A amplitude dépendante de la taille de la particule de la puissance lumineuse incidente et du coefficient de transmission de la partie réceptrice. On a donc x(t) [A-s(t) + w(t)] + j. T. H[A. s(t) + w(t)] = A. s(t) + w(t) -22f2t2 j(2ft) =A.e d e( d) +W(t) o w(t) est le signal analytique associé au bruit réel w(t) et est un bruit blanc Gaussien complexe de moyenne nulle et de même variance &2. A partir du signal analytique x(t) on définit M échantillons de la fonction d'autocorrélation ?x(m) avec N-m rx(m)= Z EJ k*(t)x(t+m)} t__N Pour un signal échantillonné bruité, la fonction d'autocorrélation de ce signal présente un énorme avantage sur le signal temporel dans la mesure o seul l'échantillon en m=0 est bruité. On peut donc considérer à partir de 2 2 que rx(m) = A2rs(m) + (2N + 1 - m)c2 (m) Plus précisément la fonction d'autocorrélation rg(m) peut être écrite sous la forme suivante: N-m 2cL2fdt2 -j2Tcfdt e-2c2fd(t+m)2 J22tfd(t -m) rs(m)= ' e C.e2.e e t=-N N-m -2o2fdt2. 2cL2fd(t+m)2 j2nfdm Y e. e. t=-N N-m -4o2fdt(t+m) 2e2fd2m2 J2Tddm = e. e t=-N 2ox2fd2m2 j2nd = (m)e e N-m -4ct2fd2t(t+m) avec y(m)= % e t=-N La phase de la fonction d'autocorrélation re(m) est définie par óó(m) = 27fdm. On peut ainsi faire une estimation sur la fréquence recherchée 2nfd en recherchant l'argument minimal de la fonction M 2 J(fd) = = ó (m)-2Mfdm2À Ceci revient à rechercher la valeur pour laquelle la fonction dérivée première J'(fd) est nulle. Soit encore avec c)d=27fd. M 2 J(Od) = d Xx(m)-mo" d m=1 M +ú(m)-2X(m)m(l, +m2(i2 m=1 l M 2 J'(0)d)= % -2ó(m)m+2m2cd m=1 =+2 Z m[m,,,d -(m)] Mm2 M J'((d)- =O ==> Zm2ód= m=;(m) m=l m=l M On a donc la fréquence estimée d définie par M- m2 m=l lO On a doric la fréquence estimée È) d =2nfd définie par 1 M 6 M d - M 2 m=1 m(m) M(M+1)(2M+1)=1 mx(m) m=1 Les valeurs 4- (m) sont calculées à partir des valeurs d'autocorrélation r(m) avec óx(m) définie par la partie imaginaire du logarithme complexe, ou bien encore comme l'arctg lm r(m) bRe r2(m) Si Imr(m) correspond à la partie imaginaire de r9(m) et Rerx(m) correspond à la partie réelle de rx(m). Selon un mode de réalisation, la sonde vélocimétrique de l'invention comprend les moyens suivants pour directement estimer rf(m) à partir des (2N+1) échantillons prélevés sur le signal délivré par le photodétecteur dont est équipé la sonde. - des moyens pour calculer la FFT de x(t) pour N instants t=-D, O, +N; - des moyens pour calculer les valeurs du spectre discret X(co); - des moyens pour calculer les valeurs du spectre discret X(co) défini par X(co) =2IX(Co)12pourco > 0 I0 pourco REVENDICATIONS 1. Sonde vélocimétrique optique comprenant des moyens pour éclairer un volume de mesure dans lequel peuvent passer des particules en mouvement relatif par rapport à la sonde, et des moyens de détection optique pour produire un signal électrique en réponse au passage d'une particule dans le volume éclairé, le spectre en fréquence de ce signal ayant une largeur centrée sur une fréquence dominante fo, représentative de la vitesse relative de la particule par rapport à la sonde, cette sonde comportant des moyens numériques pour établir un vecteur x(t) représentant N échantillons du signal électrique détecté et étant caractérisée en ce qu'elle comporte: - des moyens pour convertir le signal électrique x(t) en un signal complexe x(t) présentant une partie réelle et une partie imaginaire, ladite partie"imaginaire comprenant un terme de phase dépendante de fd; - des moyens pour établir une fonction d'autocorrélation du signal complexe, ladite fonction étant définie de la façon suivante: N-m (m)= E{,(t),(t+m)} t=-N o E {x} désigne l'espérance mathématique de la variable de x. - des moyens pour déterminer la phase 41(m) de la fonction r(m) dépendante de la partie imaginaire de ladite fonction; - des moyens pour déterminer la fréquence fd à partir de l'ensemble des couples [m,dX(m)] - des moyens pour fournir une indication de vitesse v d'une particule par rapport à la sonde, à partir de la fréquence fd. 2. Sonde vélocimétrique optique selon la revendication 1, caractérisée en ce que le vecteur x(t) est défini par l'équation suivante: x(t) = A.e2a2fd2t2 cos(2nTft) + w(t) avec ca un coefficient fixe t est le temps w(t) un terme de bruit 3. Sonde vélocimétrique optique selon l'une des revendications 1 ou 2, caractérisée en ce que les moyens pour transformer le signal réel x(t) en signal complexe x(t) comprennent des moyens du type filtre de Hilbert. 4. Sonde vélocimétrique selon l'une des revendications 1 à 3, caractérisée en ce que les moyens pour définir la fonction complexe x(t) et les moyens pour définir la fonction -r(m) comprennent les moyens suivants: * des moyens pour définir la fonction transformée de Fourier x(o), de la fonction x(t) sur N échantillons; * des moyens pour élaborer la fonction X(co) telle que: X(o) =2,X(o)I2 si (o >0 X(co) = 0 si Co FFT-1 désigne la transformée de Fourier inverse. 5. Sonde vélocimétrique selon l'une des revendications 1 à 4, caractérisée en ce que les moyens pour définir la fréquence fd comprennent des moyens pour rechercher la valeur de fréquence pour laquelle la fonction M. 2 J(f)- 4-(m)-2 Tmf 1est minimale. 6. Sonde vélocimétrique selon la revendication 5, caractérisée en ce que la valeur fd est estimée par fd 6 M 2 dXM fd =M(M+ 1)(2M + 1) m-1 '(m) M étant le nombre de points d'autocorrélation. 7. Sonde vélocimétrique selon l'une des revendications 1 à 6, caractérisée en ce que la phase4)(m) est déterminé par l'élément [im. r(m)] arctg ave. vem) c [Re. rX(m)1 Im.r2(m): la partie imaginaire de re(m) Re.rx(m): la partie réelle de rx(m)