Brevet FR2617971A1 — publié le 19890113

- 1 - 2617971 Procédé et dispositifs de granulométrie optique pour de larges gammes de mesure La présente invention concerne de manière générale le domaine de la granulométrie. Plus particulièrement, l'invention concerne un procédé et des' dispositifs de granulométrie pour la mesure par voie optique de tailles et de concentrations de particules en suspension dans un fluide, gazeux ou liquide. Actuellement dans le domaine de la granulométrie, les appareils de mesure les plus performants font appel à des procédés de mesure par voie optique. Les appareils existants sont de deux types. A un premier type appartiennent des appareils de comptage nécessitant la présence d'une seule particule à la fois dans un volume de mesure. A un second type appartiennent des appareils analysant par échantillonnage la diffusion lumineuse d'un volume-échantillon éclairé d'un nuage de particules afin de déterminer une fonction de distribution en taille des particules dans le nuage. Les appareils du premier type, tels que les appareils dits de KNOLLENBERG, utilisent les propriétés de diffusion de la lumière par une seule particule. Leur principal inconvénient réside dans la nécessité de la présence d'une seule particule à la fois dans le volume de mesure. Ces appareils ont de ce fait des performances limitées et ne couvrent qu'une gamme de mesure de taille bien définie. Les appareils du second type utilisent l'approximation de FRAUNHOFFER ou la théorie de LORENZ-MIE permettant de calculer la diffusion lumineuse d'un nuage de particules. Leur domaine d'utilisation est limité à une gamme de mesure de concentration volumique comprise entre 10 9 et 10-5. La gamme de mesure de taille dépend du modèle théorique de diffusion utilisé. Les principaux appareils du second type disponibles actuellement sur le marché sont les appareils commercialisés sous la marque MALVERN (marque déposée). Ces appareils font appel à l'approximation de FRAUNHOFFER. L'article de E.D. HIRLEMAN et de L.G. DODGE intitulé "Performance comparison of Malvern instruments laser driffraction drop size analysers" et publié dans la 3e "International Conférence on Liquid -2 2-6 72617971 Atomisation and Spray Systems" (I CLASS) Institut of Energy -London 1985- fait état de leurs performances. La distribution en taille des particules n'est considérée très bien estimée par les appareils MALVERN que pour des particules ayant un diamètre au-dessus de 5 pm. De plus, un inconvénient important de ces appareils est leur manque de précision dès que la fonction de distribution à mesurer n'est pas modélisable par une fonction simple de type ROSIN-RA*MMLER ou log-normal. Ce manque de précision est dû au procédé de calcul de la fonction de distribution faisant appel à une inversion matricielle et une estimation par la méthode des moindres carrés entraînant un résultat instable. Par ailleurs, les différents appareils existants ont une forme et un encombrement qui ne conviennent pas pour mesurer les tailles et les concentrations de particules dans un écoulement d'un fluide. En effet, une forme non adaptée et un encombrement trop important perturbent l'écoulement du fluide et modifient les mesures. Ces appareils sont inadaptés pour être portés par un véhicule en déplacement, par exemple par un aéronef afin de déterminer la distribution en taille de gouttes d'eau en suspension dans l'air. L'invention vise principalement à fournir un procédé et des dispositifs de granulométrie afin d'améliorer les performances obtenues par les dispositifs selon l'art antérieur. L'invention vise également à fournir des dispositifs de granulométrie adaptés pour déterminer la distribution en taille de particules dans un écoulement d'un fluide. A cette fin, le procédé selon l'invention pour évaluer la distribution en taille de particules dans un fluide à partir de la mesure du rayonnement diffusé dans un nombre entier fini P de directions prédéterminées de l'espace par un échantillon de fluide éclairé par un faisceau lumineux incident, est caractérisé en ce qu'à partir d'un vecteur de mesure de rayonnement ayant pour composantes des rayonnements moyens, dans les P directions, déterminés à la suite de plusieurs mesures successives et d'une matrice de coefficients théoriques de rayonnement, est calculée une première estimation de la distribution exprimée par une première fonction discrète de distribution contenue dans une fonction mathématique d'enveloppe de type prédéterminé et fixant respectivement des valeurs de concentration à des noeuds de taille pour lesquels ladite première fonction est définie, les noeuds de 3 26 1 2617971 taille étant répartis de manière à satisfaire à une relation récursive de type prédéterminé, et en ce que cette première fonction est prise comme donnée de départ pour déterminer, par un calcul itératif faisant intervenir des filtrages numériques, une relation récursive optimale à laquelle doivent satisfaire les noeuds d'une fonction discrète de distribution pour être répartis de telle manière que ladite fonction discrète de distribution présente une correlation maximale avec une fonction discrète d'incertitude de mesure découlant de plusieurs mesures de rayonnement successives dans les P directions, toutes les fonctions discrètes de distribution compatibles avec les incertitudes de mesure et dont les noeuds satisfont à ladite relation successive optimale étant ensuite calculées à partir du vecteur et de la matrice pour en déduire une évaluation statistique optimale de la distribution recherchée. Un dispositif pour la mise en oeuvre du procédé selon l'invention est caractérisé en ce qu'il comprend: des moyens pour produire un faisceau lumineux destiné à éclairer l'échantillon, des moyens pour mesurer le rayonnement diffusé dans les P directions de l'espace, zo des moyens pour convertir sous une forme numérique des mesures de rayonnement fournies par les moyens pour mesurer, des moyens mémorisant un programme et la matrice de coefficients de rayonnement pour traiter les mesures numériques de rayonnement et pour calculer ladite évaluation statistique optimale de la distribution recherchée, et des moyens pour visualiser l'évaluation calculée. Un dispositif pour évaluer la distribution en taille de particules dans un écoulement d'un fluide est caractérisé en ce que les moyens pour produire, pour mesurer, et pour convertir sont inclus dans une sonde de mesure du dispositif, et en ce que ladite sonde de mesure comprend également des moyens pour prélever une partie de l'écoulement et le conduire dans la sonde, et des moyens de réflexion optique pour diriger le faisceau lumineux produit par les moyens pour produire vers l'échantillon et les rayonnements diffusés par l'échantillon vers les moyens pour mesurer, les moyens pour produire et les moyens pour mesurer étant éloignés de l'échantillon, en aval de celui-ci dans le sens de 4 4 2 6 17971 l'écoulement, afin de donner a la sonde une forme prédéterminée ne perturbant pas l'écoulement è analyser au voisinage de l'échantillon. Le procédé selon l'invention autorise des performances élevées. La gamme de mesure en taille est de O,lpm à 300 pm; celle de mesure en concentration volumique est de 0,2 à 8.10 6. La fonction de distribution est déterminée en temps quasi réel avec une précision meilleure ou égale à 10 %. Les dispositifs de granulométrie selon l'invention sont particulièrement bien adaptés pour être portés par un aéronef ou pour être utilisés dans une soufflerie aérodynamique. Ils permettent la mesure dans un écoulement animé d'une vitesse allant de Mach 0 à Mach 0,9. En version étanche, un dispositif selon l'invention peut être immergé dans un liquide, tel que de l'eau afin de rechercher des diamètres et concentrations, par exemple de bulles d'air. D'autres avantages et caractéristiques de l'invention apparaîtront plus clairement la lecture de la description suivante de plusieurs réalisations préférées d'un dispositif de granulométrie selon le procédé de l'invention en référence aux dessins annexes correspondants dans lesquels: - la Fig. 1 montre sous une forme shématique fonctionnelle un dispositif de granulométrie selon-l'invention; - la Fig. 2 représente des fonctions mathématiques utilisées selon le procédé de l'invention afin de modéliser la fonction de distribution en taille recherchée; - la Fig. 3 montre l'implantation de différents éléments dans une sonde de mesure comprise dans un dispositif de granulométrie selon l'invention; et - la Fig. 4 est une vue en perspective de la sonde de mesure selon l'invention. En référence a la Fig. 1, le dispositif selon l'invention est constitué d'une sonde de mesure 1 et d'un micro-calculateur 2. La sonde de mesure 1 comprend essentiellement un générateur de faisceau laser 11, des première et seconde barrettes de mesure de rayonnement l2 et 13, et un convertisseur analogique-numérique 14. Le générateur de faisceau laser 11, par exemple de type D gaz, fournit un faisceau laser polarisé et expandé 111. -5 2 62617971 Le faisceau laser 111 est réfléchi par un miroir plan, ou prisme réflexion totale 15, afin de modifier sa direction de propagation et l'amener à traverser un fluide 16 contenant des particules en suspension de différents diamètres. Le fluide 16 est dirigé dans la sonde 1 suivant une direction d'écoulement prédéterminée 161 et est traverse perpendiculairement par le faisceau laser réfléchi 111a. L'intersection du faisceau 111 et du fluide 16 est un volume 17 sensiblement cylindrique. Le volume 17 constitue un volume de mesure prédéterminé par le choix du diamètre du faisceau laser 111 et par le dimensionnement d'un guide creux semi-circulaire 18 faisant partie d'une platine dans la sonde 1 dirigeant l'écoulement du fluide 16. Le volume de mesure 17 contient une pluralité de particules sensiblement sphériques de différents diamètres O répartis selon une fonction de distribution en taille f () dont une évaluation est l'objet de la mesure. Chacune des particules dans le volume de mesure 17 constitue un élément diffractant le faisceau laser 111a et diffuse une onde lumineuse de même longueur d'onde que le faisceau laser. La théorie rigoureuse de LORENZ-MIE permet de calculer pour une quelconque direction donnée de l'espace la quantité différentielle efficace de rayonnement lumineux diffusé dans cette direction par une particule sphérique éclairée par une onde lumineuse. Cette quantité différentielle est fonction d'un indice de réfraction complexe m de la particule et de son diamètre 0. Pour un nuage de particules, le rayonnement lumineux efficace 2Z5 diffusé par ce nuage dans une direction d de l'espace est égal, à l'infini, à la somme des quantités différentielles diffusées dans la même direction d par les différentes particules constituant le nuage. Le rayonnement lumineux diffusé par un nuage de particules dans une direction d de l'espace est donné par l'équation intégrale: b(d) =Ie maxi K(d,b,m). f(l). d mini o b(d) est le rayonnement lumineux efficace diffusé par le nuage dans la direction d, f(e) est la fonction de distribution en taille du nuage, K(d,i,m) est une fonction de rayonnement selon la théorie de LORENZ-MIE donnant pour chacune des particules du nuage, définie par son diamètre b -6- -6 - 2617971 et son indice de réfraction m, la quantité différentielle de rayonnement lumineux qu'elle diffuse dans la direction d, et o O et maxi mini maxi définissent un intervalle de diamètres dans lequel est compris le diamètre e d'une particule quelconque du nuage. Les barrettes de mesure de rayonnement 12 et 13 sont éclairées par le rayonnement diffusé par le volume de mesure 17. Des miroirs plans de renvoi, ou prismes, 18 et 19 respectivement associés aux barrettes 12 et 13 sont prévus dans la sonde 1 afin de réfléchir vers les barrettes 12 et 13 des parties du rayonnement diffusé par le volume 17. Le miroir 18 est disposé à l'opposé du miroir 15, par rapport à la direction d'écoulement 161 du fluide 16, et reçoit une partie du rayonnement diffusé par le volume de mesure 17 observé suivant la direction de propagation du faisceau laser l1la, direction de propagation perpendiculaire au plan de polarisation du faisceau. La barrette 12 est alignée suivant un axe 121 parallèle è la direction de propagation du faisceau laser Mlla et est orientée pour être éclairée par le rayonnement réfléchi par le miroir 18. Le miroir 19 est placé perpendiculairement au plan de polarisation du faisceau Mlla et reçoit le rayonnement diffusé par le volume de mesure 17 observé suivant une direction perpendiculaire à la direction de propagation du faisceau illa. La barrette 13 est alignée suivant un axe 131 incliné d'un angle prédéterminé a par rapport à la direction de propagation du faisceau laser ll1a afin de permettre un nombre plus élevé de mesures. La barrette 13 est orientée vers le miroir 19 et est éclairée par le rayonnement réfléchi par celui-ci. Les barrettes 12 et 13 comprennent chacune une pluralité de cellules photoélectriques, ou de phototransistors, alignées, 121a et 131a, constituant autant de récepteurs de mesure de rayonnement. Les barrettes 12 et 13 mesurent le rayonnement diffusé par le volume de mesure 17 suivant respectivement des premier et second groupes de directions de mesure d compris respectivement dans des premier et second plans de mesure. Le premier plan de mesure correspondant à la barrette 12 est perpendiculaire au plan de polarisation du faisceau laser Mla et passe par l'axe 121 de la barrette 12 et un axe de symétrie 181 du miroir 18. Le second plan de mesure passe par des centres de symétrie 191 et 171 respectivement du miroir 19 et du volume de mesure 17 et est incliné de l'angle a par rapport à un plan passent par les mêmes centres -7 6- - 7 - - 2617971 de symétrie 191 et 17l et perpendiculaire au plan de polarisation du faisceau laser Illa. Les barrettes de mesure 12 et 13 fournissent des vecteurs de mesure ME constitué chacun par P, typiquement égal à 32, mesures respectivement dans P directions différentes d1 ^ dp comprises dans les premier et second plans de mesure. Le convertisseur analogique-numérique 14 reçoit à des entrées parallèles 141 les vecteurs de mesure ME et les délivre sous forme numérique au micro-calculateur 2. Le micro-calculateur 2 comprend une unité centrale 21, une mémoire de coefficients 22, une mémoire de programme 23, et des moyens de visualisation 24. L'unité centrale 21 reçoit les vecteurs de mesure ME fournis par le convertisseur 14 de la sonde 1 et les traite en utilisant des coefficients mémorisés dans la mémoire 22 selon un programme de calcul et de traitement mémorisé dans la mémoire 23, afin de déterminer la fonction de distribution en taille f(4) des particules contenues dans le volume de mesure 17'. De manière générale, les fonctions de distribution en taille 2 réelles ont des formes d'onde naturelles connues et l'on peut considérer qu'une fonction de distribution en taille quelconque à déterminer est modélisable par une somme de fonctions mathématiques de type analogue. Selon le procédé de l'invention, l'on tient pour acquis qu'une fonction de distribution en taille d'un nuage de particules est modélisable par une somme de fonctions mathématiques exponentielles connues de type B - spline. Ces fonctions sont définies par des produits de convolution de fonction créneau et sont représentées par des courbes monotones avec une dérivée seconde minimale; leur valeur est nulle à l'extérieur de l'intervalle compris entre le premier et le cinquième noeud. La fonction de distribution f(4) recherchée peut être exprimée par l'égalité suivante: jW= J S4 j=J f() = z a.. S. (4), j = 1 O J o j est un indice entier variant entre 1 et une valeur maximale J. Les fonctions S.(4) sont des fonctions B-splines, et les coefficients j 8 2-6 72617971 multiplicatif a. sont des coefficients réels positifs respectivement J associés aux fonctions S(4). La mémoire 22 mémorise une matrice de coefficients T comprenant P.Q coefficients, tl à tpQ, o Q est un entier prédéterminé de valeur élevée typiquement égale à 150. Les coefficients t de la matrice T sont déterminés théoriquement par application de la théorie de LORENZ-MIE. Un coefficient t de la matrice T, o p et q sont des entiers pq compris respectivement entre 1 et P, et 1 et Q, représente un poids de rayonnement dans la direction d d'une particule ayant un diamètre q. La matrice T est calculée à l'aide d'un programme de calcul implanté à cet effet dans la mémoire de programme 23 du micro-calculateur 2, pour un indice de réfraction m donné et correspondant à la nature du nuage de particules. La matrice T représente sous une forme discrète la relation qui lie le rayonnement b diffuse dans l'espace par le volume de mesure 17 à la fonction de distribution en taille f(4) des particules incluses dans le volume de mesure 17. L'égalité suivante peut être exprimée B = T.F, o B est un vecteur de rayonnement ayant P = 32 composantes b à b i P représentant respectivement des rayonnements diffusés par le volume de mesure 17 dans les P = 32 directions d1 à dp dans les premier et second plans de mesure, et F est un vecteur représentant une forme discrète de la fonction de distribution f() et ayant Q = 150 composantes f1 à fQ, respectivement égales aux Q = 150 concentrations respectives de Q = 150 particules de diamètres différents dans le volume de mesure 17. De l'égalité ci-dessus, il découle qu'il est possible de calculer le vecteur F représentant une nouvelle forme discrète de la fonction de distribution f(p) et ayant au maximum J = 31 composantes, o J est un entier prédéterminé, connaissant le vecteur B et en inversant une forme dégénérée U de la matrice T. L'unité 21 calcule à partir de plusieurs vecteurs de mesure ME successifs fournis par la sonde 1 et à l'aide d'une table de coefficients de correction une estimation statistique du vecteur B. Chacune des cellules 121a, 131a des barrettes 12 et 13 ayant une courbe de réponse spécifique, la table de coefficient de correction contient -9 2 2617971 des coefficients à appliquer aux mesures fournies par les cellules 121, 131 afin de déterminer des rayonnements correspondants. Un vecteur moyen de rayonnement Bm et un vecteur d'écarts types correspondant aB sont calculés par l'unité 21. Les vecteurs Bm et B comprennent chacun P = 32 composantes respectivement bm1 à bu et ob à abp. A partir du vecteur moyen Bm et de la matrice U, l'unité centrale 21 détermine un vecteur Fa qualifiant une distribution et représentant les coefficients a. correspondants des J t 31 fonctions B-spline cubiques dont les noeuds 11 JN sont définis a l'origine par les relations suivantes: 1 = mini' %JN = maxi' j,n+l = C.,n, et 4j+l,n = C. j,n, o %mini et maxi sont des valeurs prédéterminées délimitant un intervalle dans lequel sont comprises toutes les valeurs des noeuds, ou n et j sont des entiers variant respectivement entre 1 et N, et 1 et J, et o C est un coefficient multiplicatif ayant une valeur C pour le o vecteur Fa. o La matrice T comprenant P.Q coefficients est transformée en la matrice U comprenant P.J coefficients u calculés à partir de la relation: UP =f jN K(d,D,m). S.(). d, Une fonction S.(ó) est montrée fig. 2; Sj($) est définie à partir d'un - nombre typique de noeuds N égal à 5. Le calcul du vecteur Fa0 à partir du vecteur Bm et de la matrice U est réalisé selon des procédés d'inversion statistique classiques dans l'art antérieur et connus par l'homme du métier. Le vecteur Fa0 peut être une solution optimale déterminée par le calcul d'un écart quadratique moyen minimum par rapport à une solution mathématique exacte égale à Bm.U '. Le vecteur de distribution modélisée Fa0 ainsi calculé est une solution instable. Le coefficient multiplicatif C entre les noeuds du vecteur Fa à une valeur calculée C0 qui peut être optimisée afin de définir le vecteur Fa0 par - 10 - 2617971 des noeuds pour lesquels le vecteur Fa0 est corrgl6 de manière maximale avec un vecteur d'incertitudes de mesure. Selon le procédé de l'invention, un filtrage numérique est réalisé pour déterminer une valeur optimale Copt du coefficient multiplicatif C. Plusieurs produits de convolution PC01, PC02,... PCoi, o I est un entier prédéterminé, sont calculés entre un vecteur d'écarts types cF = U -.aB et respectivement plusieurs vecteurs Fao01, Fao02,...FaOI. Les vecteurs Fao1 à Faoi sont déduits du vecteur Fa0 par variation du - coefficient multiplicatif C autour de la valeur CO. Aux vecteurs Fa 0 O 10 FaOi correspondent respectivement des valeurs Cll à Ci autour de la valeur C0 du coefficient multiplicatif C. Parmi les produits de convolution PC01 à PCoi correspondant respectivement aux valeurs C] à Cli du coefficient multiplicatif C, un produit de convolution a une valeur maximale. A ce produit de convolution de valeur maximale correspond respectivement une valeur C1 du coefficient multiplicatif C. La valeur C1 du coefficient multiplicatif C est une valeur optimale intermédiaire à laquelle correspond un vecteur Fa, parmi les vecteurs Fa0 à Faoi. Dans le cas o la valeur C1 est supérieure ou égale à la valeur C0 du coefficient multiplicatif C, le vecteur Fa! est légitimement interpolé du vecteur FaO, car il comprend un nombre de composantes, correspondant aux noeuds des fonctions B-splines, au plus égal à celui du vecteur Fa0 et ne porte donc que des informations contenues déj> initialement dans le vecteur FaO. En d'autres termes, le spectre de fréquence de la fonction discrète correspondant au vecteur Fa1 est inclus dans le spectre de fréquence de la fonction discrète correspondant au vecteur FaO. La valeur C1 est dans ce cas la valeur optimale recherchée Copt du coefficient multiplicatif C. Dans le cas o la valeur C1 est inférieure à la valeur C0 du coefficient multiplicatif C, le vecteur Fa1 comprend un nombre de composantes supérieur au vecteur Fa0 et n'eàt pas une interpolation légitime du vecteur FaO. Un autre vecteur de distribution modélisée Fbo, o représentant la fonction f(f) et dont les composantes correspondent à des noeuds de fonctions B-spline liées par la valeur C1 du coefficient multiplicatif C est calculé par l'unité centrale 21 à partir du vecteur moyen Bm et de la matrice U, par inversion statistique. Le vecteur Fb0 est ensuite traité de manière analogue au vecteur FaO. - il - 2617971 La détermination de la valeur optimale Copt du coefficient multiplicatif C nécessite un processus de calcul itératif, comme décrit ci-dessus. En pratique la valeur Copt est typiquement obtenue après deux ou au maximum trois itérations. La valeur optimale Copt du coefficient multiplicatif C étant &éterminée, l'unité centrale 21 calcule par inversion statistique et décalage des noeuds des fonctions B-splines tous les vecteurs de distribution F définis ?- partir de la valeur Copt du coefficient multiplicatif C et compatibles avec les vecteurs Bm et oB. Un vecteur optimum Fopt parmi tous les vecteurs F présente un écart quadratique optiu Fo-pritulsv1er pOeieu cr udaiu moyen minimum par rapport à la solution moyenne U l.Bm. Le vecteur Fot opt représente de manière optimale la fonction de distribution en taille f(%) des particules dans le volume 17. L'unité centrale 21 fournit le résultat des calculs aux moyens de visualisation 24. La fonction modélisée f(ó) peut également être visualisée sous une forme statistique calculée à partir de tous les vecteurs F. Les miroirs 15, 18 et 19 sont prévus dans la-sonde de mesure 1 afin de permettre un positionnement convenable, montré Figs. 3 et 4, du générateur de faisceau laser 11, des barettes 12 et 13, du convertisseur 14, et d' alimentations respectives, non représentées à la Fig. 1. Ces différents éléments sont positionnés en aval, selon la direction d'écoulement 161 du fluide 16, du point de mesure matérialisé par le volume 17 et ceci afin de ne pas perturber l'écoulement du fluide mesurer. Tous les éléments de la sonde de mesure sont solidaires d'une platine capotée. Cette platine est concue de manière à amener l'écoulement du fluide à- angle droit du faisceau laser par le guide creux semicirculaire 18, sans pour autant le canaliser ce qui créersit des risques de détérioration des dimensions des particules ou gouttes par modification des caractéristiques de l'écoulement. Dans un même but, un déflecteur 10, montré à la Fig. 3, équipé d'un piège à couche limite, est placé à l'avant de la sonde 1 -afin de minimiser les perturbations qu'engendrerait la présence de celle-ci au sein de l'écoulement à analyser. La sonde 1 étant destinée à fonctionner dans des conditions de températures variables, le déflecteur 10 est muni d'éléments chauffants - 12 - - 2 6 2 17971 101, tels que des résistances, assurant l'intégrité de sa forme lors de variations de température. - 13 -2617971 R e v e n d i c a t i o n s 1 - Procédé pour évaluer la distribution en taille de particules dans un fluide à partir de la mesure du rayonnement diffusé dans un nombre entier fini P de directions prédéterminées de l'espace ( d1 à dp) par un échantillon de fluide (17) éclairé par un faisceau lumineux incident (111a), caractérise en ce qu'à partir d'un vecteur de mesure de rayonnement (Bm) ayant pour composantes-des rayonnements moyens (bm bmp) dans les P directions déterminées à la suite de plusieurs mesures successives et d'une matrice de coefficients théoriques de rayonnement (U) est calculée une première estimation de la distribution exprimée par une première fonction discrète de distribution (FaO) contenue dans une fonction mathématique d'enveloppe (E a..S.(f)) de type prédéterminé et J J fixant respectivement des valeurs de concentration à des noeuds de taille (4!1 ià JN) pour lesquels ladite première fonction est définie, les noeuds de taille (11 à 4JN) étant répartis de manière à satisfaire à une relation récursive ( j,n+l = C.jn et j+ln = C.jn) de type prédéterminé, et en ce que cette première fonction est prise comme donnée de départ pour déterminer par un calcul itératif faisant intervenir des filtrages numériques une relation récursive optimale ( Copt) à laquelle doivent satisfaire les noeuds (11 aà $JN) d'une fonction discrète de distribution (F) pour être répartis de telle manière que ladite fonction discrète de distribution (F) présente une corrélation maximale avec une fonction discrète d'incertitude de mesure (oB.U b) découlant de plusieurs mesures de rayonnement successives dans les P directions, toutes les fonctions discrètes de distribution (F) compatibles avec les incertitudes de mesure (aB.U 1) et dont les noeuds (fil à JN) satisfont à ladite relation récursive optimale (Copt) étant ensuite calculées à partir du vecteur de mesure (Bm)-et de la matrice (U) pour en déduire une évaluation statistique optimale de la distribution recherchée. 2 - Procédé conforme à la revendication 1, caractérisé en ce que les coefficients de rayonnement (ull à upQ) de ladite matrice (U) sont déterminés à partir de la théorie de LORENZ-MIE. 3 - Procédé conforme à la revendication 1 ou 2, caractérisé en ce que ladite fonction mathématique d'enveloppe est une somme de fonctions mathématiques de type B-spline (E a..S.(4)). J J -14- - 14 - 2617971 4 - Procédé conforme à la revendication 1 ou 3, caractérisé en ce que ladite relation récursive lie deux noeuds successifs quelconques (jn+ j n) d'une fonction discrète de distribution (F) par un coefficient multiplicatif. 5 - Dispositif pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il comprend des moyens (11) pour produire un faisceau lumineux destiné éclairer l'échantillon (17), des moyens (12,13) pour mesurer le rayonnement diffusé dans les P directions (d1 3 dp) de l'espace, des moyens (14) pour convertir sous une forme numérique des mesures de rayonnement (ME) fournies par les moyens pour mesurer (12,13), des moyens (2) mémorisant un programme (23) et la matrice de coefficients de rayonnement (22) pour traiter les mesures numériques de rayonnement et pour calculer ladite évaluation statistique optimale de la distribution recherchée, et des moyens (24) pour visualiser l'évaluation calculée. 6 - Dispositif conforme à la revendication 5, caractérisé en ce que les moyens pour produire (11) sont un générateur de faisceau laser. 7 - Dispositif conforme à la revendication 5 ou 6, caractérisé en ce que les moyens pour mesurer (12,13) sont des barrettes de mesure constituées chacune d'une pluralité d'éléments de mesure alignés tels que des cellules photoélectriques ou des phototransistors. 8 - Dispositif conforme à la revendication 5 ou 7, caractérisé en ce que lesdits moyens pour traiter et calculer (2) sont un micro-calculateur. 9 - Dispositif conforme à l'une quelconque des revendications 5 9, caractérise en ce que les moyens pour produire (11), pour mesurer (13), et pour convertir (14) sont inclus dans une sonde de mesure (1) du dispositif. - Dispositif pour évaluer la distribution en taille de particules dans un écoulement d'un fluide (16) conforme à la revendication 9, caractérisé en ce que ladite sonde de mesure (1) comprend également des moyens (10,18) pour prélever une partie de l'écoulement et le conduire dans la sonde (1), et des moyens de réflexion optique (15, 18, 19) pour diriger le faisceau lumineux (111, - 15 - 2617971 Illa) produit par les moyens pour produire (11) vers l'échantillon (17) et les rayonnements diffusés par l'échantillon (17) vers les moyens pour mesurer (12,13), les moyens pour produire (11) et les moyens pour mesurer (12,13) étant éloignés de l'échantillon, en aval de celui-ci dans la direction de l'écoulement, afin de donner à la sonde (1) une forme prédéterminée ne perturbant pas l'écoulement à analyser au voisinage de l'échantillon (17). 11 - Dispositif conforme à la revendication 10, caractérisé en ce que les moyens pour prélever et conduire comprennent un guide creux semicirculaire 18 et un déflecteur (10) équipé d'un piège - couche limite, dans ledit déflecteur (10) étant noyées des résistances chauffantes (101) afin d'assurer l'intégrité de sa forme lors de variations de température. 12 - Dispositif conforme à la revendication 10 ou 11, caractérisé en ce que les moyens pour diriger comprennent des miroirs plans et/ou des prismes à réflexion totale (15, 18, 19).