Brevet FR2749942A1 — 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 o10 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 devant les 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 l'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 de mesure éclairé, le spectre en fréquence de ce signal a une largeur AF centrée sur une fréquence dominante Fo représentative de la vitesse relative de la particule par rapport à la sonde, on a: Fo = i, i étant l'interfrange. Cette sonde comporte des moyens numériques pour établir un vecteur x(t) représentant N échantillons numériques du signal électrique détecté et étant caractérisée en ce qu'elle comporte en outre: - des moyens pour calculer la densité de probabilité p(x,O) du vecteur signal x, constitué des échantillons de signal {x(t)}t=0+ +D, en fonction d'un vecteur paramètre e=[A,fd,c2] caractérisant complètement x; - des moyens pour déterminer la fréquence fo pour laquelle la probabilité est maximale; des moyens pour fournir une indication de vitesse v d'une particule par rapport à la sonde, à partir de la fréquence fo. Le vecteur signal x = {x(t)}t=0,+ l,... +D est défini par x(t) = A e2'x2f2t2 cos(2rft)+co(t) o ce est un coefficient fixe, t est le temps, co(t) est un bruit gaussien d'énergie c2. Le vecteur 0 est quant à lui défini comme 0 = [A, fd, a2] Les moyens utilisés dans l'invention pour déterminer la fréquence f permettent d'atteindre une précision sur ladite fréquence supérieure à celle obtenue par une méthode classique de calcul utilisant la transformée de Fourier. En effet, le signal recueilli au niveau de la photodiode étant d'une durée déterminée, le spectre en fréquence obtenu par ce moyen d'analyse présente une certaine largeur spectrale en fréquence entraînant une imprécision sur la détermination de ladite fréquence. 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 parmi lesquelles: - la figure 1 illustre un exemple de sonde vélocimétrique selon I'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 l'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 I'aéronef 6. L'optique 51 a pour fonction de faire converger les faisceaux lumineux L1 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 phototiode, 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 fournit un signal i(t) dont le modèle mathématique est connu: i(t) = A. s(t) = A. e-2(v/w)2t2 cos(27 v /.t) + co (t) avec t = O, +1,... + w/v o I'amplitude A dépend de la taille de la particule, de la puissance lumineuse incidente et du coefficient de transmission de la partie réceptrice; v représente la vitesse de la particule; 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; * o(t) bruit blanc supposé gaussien d'énergie o2. avec ca = I/w et fd = vil on obtient i(t) = A. s(t)+ o(t) -= A. e2(2fd2t2 cos(2ltfdt)+wo(t) - la sonde selon l'invention comprend des moyens pour prélever N échantillons du signal, à des instants t = -D,.. ., O,..., +D, I'ensemble de ces échantillons permet de définir le vecteur x tel que: x = A.s + c); - elle comprend également des moyens pour définir un estimateur du maximum de vraisemblance, permettant d'obtenir une estimée fo de fd qui rende les plus vraisemblables possibles les mesures effectuées sur le signal et les valeurs obtenues à partir du modèle théorique pour la fréquence fo; - pour définir cette fonction de probabilité on peut déterminer un vecteur 0 (relatif au modèle mathématique) dépendant des paramètres (A, f, ). Le bruit co(t) peut être estimé égal à une valeur moyenne c2. 0 ne dépend alors que des paramètres (A, f). La fonction de densité de probabilité p(x,0) peut alors être définie comme suit. Le vecteur s contient les données théoriques correspondant à 0 alors que x contient les données mesurées. p(x, O) = 2 -2î/.exp lx - As 112}(1) p(x, ') = (2na 2)(2D+ l)/2.exp. 11 2 chercher à maximiser la fonction p(x,O) est équivalent à chercher à maximiser le logarithme de p(x,0) soit la fonction A(x,0) avec: A(xO) (2) A(x, 0)=-2D+l) 2 log2 -i2Dé+ ld 2r 2lix Asa 2 En définissant une estimée du paramètre A avec sT.x sT.s sT étant la transposée de s. L'équation (2) devient l'équation suivante (3): _2 12 +1 2 i T.x2 A(xO)z(+)1g7-( 2 XsT. (3) log27i -22 2 sT.s On est donc ramené à rechercher le minimum de la fonction sT.x 2 avec x = [x(-D),..., x(o), ., x(D)]T s = [s(-D),..., s(o),..., s(D)]T et s(t) = e- 2a2f2t2.cos(2 rft) Comme la fonction J(f) est une fonction non linéaire de f, il n'existe pas de solution analytique pour déterminer la valeur fo pour laquelle cette fonction J(f) est minimale, c'est-à-dire possédant une dérivée première..DTD: nulle J'(f) et une dérivée seconde positive J"(f). Une solution proposée dans le cadre de la sonde selon l'invention est de procéder de manière itérative à partir d'une première valeur fo1 définie, par transformée de Fourier du signal recueilli au niveau du photodétecteur. Plusieurs types d'itérations peuvent être utilisés pour converger vers la valeur estimée fo recherchée. On peut notamment utiliser un algorithme de type Gauss-Newton dans lequel l'incrément de fréquence estimée est déterminée par l'équation (4) suivante: f(n+l) = (f) = dJ(f) d2 j(f) (4) fl "If(n) df d2f On détermine une valeur estimée fo = f(n) lorsque la différence f(n+1) - f(n) est inférieure à un seuil prédefini sT x en considérant s = x-- s, x c, S on peut noter J(f) = |s 2 et l'on a J'(f) = 2s'T.s si s' d-(f df la procédure de Gauss Newton repose également sur l'approximation suivante donnée par l'équation (5): J"(f) = 2s"T.' + 2 - 2 s'T s' (5) on peut montrer que: [(XT.s'f)(sT.s)- 2(XT.s)(sT.sf)]s + (XT.s)(sT. s)sf (sT. s)2 avec s'f = ds(f) df Dans l'estimateur de vraisemblance utilisé dans les moyens de traitement du signal de la sonde vélocimétrique, on peut avantageusement approximer les différents produits scalaires sT.s _ _ 4oLf sT.sf sf 8o.f2 -T (3cO2 + 22) s'.s'f-O3 f 'f 16o3f3 En assimilant les sommes discrètes sur 2T+1 points à des sommes discrètes infinies, puis en assimilant, de même, ces sommes discrètes infinies à des intégrales continues de + oo a - oo. En reportant dans les dérivées premières et secondes, on obtient les expressions suivantes J'(f) - {2f(xT s'f) + (xT.s)}(xT. s) JL(f) 2 {t2e2f(XTsu f (XTce. s) +2(c2 + d2)(xTe s) La formule d'iteration sur la fréquence devenant: f(n+l) -f (n)f [(-22f). {2f.(x s'f)+ (x.s)}(xT. s)l a 2[2f. (XT S f) + (XT s)] + 2(0r2 + x2) (XT' S)2 f(n) Chaque iteration nécessite ainsi le calcul de xT.s et xT.s'f, soit 2(2D+1) multiplications et additions plus environ 7 multiplications et 2 additions sur des scalaires. L'algorithme utilisé comprend ainsi les trois phases: initialisation de l'algorithme avec f01= fFFT (fFFT: fréquence obtenue par la transformée de Fourier) - actualisation de la fréquence estimée par la formule d'iteration - I'arrêt du calcul est fixé par la validation d'un critère de précision sur l'estimation. Ce critère dépend de la valeur d'un paramètre 6 défini de la manière suivante: |f (n+l) _ f (n) que l'on peut avoir dans l'estimation de la fréquence. De manière générale, pour tout estimateur que l'on peut construire sur le signal, il existe une mesure de l'efficacité - en terme de précision des résultats qu'il fournit - de cet estimateur, constituée par les bornes de Rao-Cramer (correspondant à des bornes associées à un modèle de bruit, ici supposé bruit blanc gaussien). Les bornes de Rao-Cramer (usuellement notées CRB) constituent en termes de variance, la performance ultime de tout estimateur que l'on pourrait construire sur le signal. C'est-à-dire qu'elles donnent la limite inférieure de la variance des paramètres estimés du signal par rapport à son modèle théorique. La méthode d'estimation proposée dans l'invention appelée Estimateur du Maximum de vraisemblance est en terme de précision de mesures, le meilleur estimateur possible, atteignant les bornes de Rao- Cramer. Ainsi on peut considérer que dans 99 % des cas, la mesure effectuée sur f est comprise dans l'intervalle: [f0 - / RB(f);fo + CRB(fo0)] Le calcul de la borne CRB(fo) est effectué de la manière suivante: A partir de l'équation de probabilité (2) on définit la matrice FIM (Fisher Information Matrix) dont les éléments sont donnés par l'équation: Fij = -E ' J J avec O{ est la ième composante du vecteur 0 et les bornes CRB(O) sont définies par l'équation CRB(O) = F-1 on peut notamment montrer que CRB(fo) répond à l'approximation suivante: o2 8oLE3f3 CRB(fo) d O A2 -JC(OE 72)+n 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 de mesure éclairé, 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 numériques du signal électrique détecté et étant caractérisée en ce qu'elle comporte en outre: - des moyens pour calculer la densité de probabilité p(x,O) du vecteur signal x, constitué des échantillons de signal {x(t)}t=0+ _D, en fonction d'un vecteur paramètre 0=[A,fd,a2] caractérisant complètement x; - des moyens pour déterminer la fréquence fo pour laquelle la probabilité est maximale; - des moyens pour fournir une indication de vitesse v d'une particule par rapport à la sonde, à partir de la fréquence fo. 2. Sonde vélocimétrique optique selon la revendication 1, caractérisée en ce que le vecteur signal x = {x(t)}t=0 ++lD est défini par x(t) = A.e-2(x2f2t2 cos(2nft)+0)(t o c est un coefficient fixe, t est le temps, c)(t) est un bruit gaussien d'énergie a2. 3. Sonde vélocimétrique optique selon l'une des revendications 1 ou 2, caractérisée par l'expression suivante de la densité de probabilité des données mesurées en fonction du vecteur paramètre 0: log P(x,O) = (2D+1)/2. log 27c - (2D+1)/2. Ioga2 - 1/2c211lx-A.sll o s(t) = e2(xc2f2t2cos(2 ft) et o -D,..., 0,..., +D correspondent aux instants o sont prélevés les N échantillons. 4. Sonde vélocimétrique optique selon l'une des revendications 1 à 3, caractérisée en ce que l'amplitude A est estimée à la valeur dont la valeur est donnée par la formule suivante = sT.x/sT.s o sT est la transposée du vecteur s. 5. Sonde vélocimétrique optique selon la revendication 4, caractérisée en ce que les valeurs f possibles pour calculer la probabilité P(x,O) sont une succession de valeur fn o n est un indice et la valeur f(n+1) se déduit de la valeur f(n) selon une itération de type Gauss-Newton. 6. Sonde vélocimétrique optique selon la revendication 5, caractérisée en ce que f(n+l) = f(n) _ J i(n) avec J' et J" les fonctions dérivées première et seconde de la sT. fonction J1 (f) définie par J(f) = llsll2 avec s = - s-.s sT.s 7. Sonde vélocimétrique selon la revendication 6, caractérisée en ce que J'(f) est calculé par le produit scalaire 2s'T.., F'T est le valeur transposé de s'. 8. Sonde vélocimétrique selon la revendication 7, caractérisée en ce que la dérivée seconde J"(f) est calculée par le produit scalaire 2'". s. 9. Sonde vélocimétrique selon l'une des revendications 7 ou 8, caractérisée en ce que les produits scalaires sT.s, sT.s'f, s'fT.s'f sont approximées par les valeurs suivantes: sT.s./ sT.sf - - 'f T.s 'f -(3c2 + 2n2) - 4c.f ' - 8f2 - 16.cc3f3