EP2912475B1 - Procede de detection de conditions de turbulence utilisant l'interaction d'un faisceau laser avec un film mince photochromique et dispositif mettant en oeuvre ledit procede - Google Patents

Procede de detection de conditions de turbulence utilisant l'interaction d'un faisceau laser avec un film mince photochromique et dispositif mettant en oeuvre ledit procede Download PDF

Info

Publication number
EP2912475B1
EP2912475B1 EP13779888.0A EP13779888A EP2912475B1 EP 2912475 B1 EP2912475 B1 EP 2912475B1 EP 13779888 A EP13779888 A EP 13779888A EP 2912475 B1 EP2912475 B1 EP 2912475B1
Authority
EP
European Patent Office
Prior art keywords
images
turbulence
environment
assessing
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Not-in-force
Application number
EP13779888.0A
Other languages
German (de)
English (en)
Other versions
EP2912475A1 (fr
Inventor
Régis BARILLE
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Centre National de la Recherche Scientifique CNRS
Universite dAngers
Politechnika Wroclawska
Original Assignee
Centre National de la Recherche Scientifique CNRS
Universite dAngers
Politechnika Wroclawska
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Centre National de la Recherche Scientifique CNRS, Universite dAngers, Politechnika Wroclawska filed Critical Centre National de la Recherche Scientifique CNRS
Publication of EP2912475A1 publication Critical patent/EP2912475A1/fr
Application granted granted Critical
Publication of EP2912475B1 publication Critical patent/EP2912475B1/fr
Not-in-force legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/001—Full-field flow measurement, e.g. determining flow velocity and direction in a whole region at the same time, flow visualisation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/26—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/10—Selection of transformation methods according to the characteristics of the input images

Definitions

  • the field of the invention is that of methods for detecting the turbulence of a medium, exploiting the disturbances generated on the propagation of a laser beam and that of compact devices implementing such a method.
  • This type of process can also be used to measure turbulence rates (more or less strong in a given environment) than to measure the effects of turbulence on the propagation of a laser beam.
  • the turbulence analysis can find many applications, via compact and discrete means, especially in the field of the environment, detection of temperature rise, air movement in a city for example.
  • One of the methods is to simultaneously make differential measurements of arrival angles and intensity fluctuations of the laser beam from a reference using a CCD camera since the measurements of spatial fluctuations are very sensitive to the wavelength of the path traveled by the laser beam ( A. Consortini et al., Optics Communications, 216, 19-23, (2003) )). This fluctuation depends cubically on the length of the optical path traveled.
  • the tilt of the wavefront is analyzed using polynomial techniques from a shark-hartmann analysis method ( AV Sergryev et al., Appl. Opt., 50 (20), 3519-3528, (2011) )) and requires a complex and expensive laser wavefront analysis system. Fluctuations in the angle of arrival or tilt of the wavefront are used to measure the characteristic parameters of atmospheric turbulence.
  • the structure parameter of the refractive index or C n 2 is a measure of the strength of the optical turbulence along the propagation path of the laser beam. In free-space laser propagation the measurement of angular arrival fluctuations is a basic step in the study of atmospheric turbulence.
  • the subject of the present application is an easy method of implementation, opening the way to devices for evaluating the turbulence of a medium, which may notably be related to ambient air heating phenomena. , easy to implement and can be deployed massively.
  • the method comprises a step of acquiring diffraction images of the laser beam by the created surface grating, said images having diffraction spots whose characteristics are dependent on the turbulence parameters.
  • the device for evaluating the turbulence parameters of a medium comprises a set of measurement units, said measurement units being connected to said reception and processing device.
  • Another subject of the invention is the use of the device for evaluating the turbulence parameters of an ambient medium according to one of the invention, for making measurements of air movements, for example in an urban environment.
  • the invention will be described in the context of an azobenzene film, but may equally well use other types of photochromic organic materials.
  • the effect used in the present invention corresponds to a holographic recording in which the fringes build themselves during the process of molecular self-organization and print a long-term statistical average of the effect of the turbulent medium. on laser propagation.
  • Organic holographic materials have the advantage over their inorganic counterparts of being structurally versatile, easily fabricated and inexpensive.
  • This type of photochromic material such as an azopolymeric material, is capable of self-organizing on the surface under the effect of the light of a laser to produce a pattern in the form of a diffraction grating, for simultaneously measuring several characteristics of said turbulence.
  • the molecules of azobenzene have been studied for more than half a century for their particular optical properties and, even today, new phenomena related to their photoisomerization are discovered. Among these phenomena, the deformation of photo-induced material in photochromic materials containing these molecules is undoubtedly one of the most remarkable. Successive cycles of isomerization of the molecule may result in material displacement over distances of several hundred nanometers.
  • the molecules of the azobenzene family are characterized by two benzene rings linked by a double bond between two nitrogen atoms. They can exist in two isomeric forms: the trans (or E) form and the cis (or Z) form. The transition from the trans form to the cis form is excited by absorption of a photon and the reverse transition from the cis form to the trans form can be either optically excited or thermally activated.
  • This effect can be observed on a polymer thin film containing azobenzenes capable of making a topographic network on the surface of the film after exposure by an interference structure created by two laser beams.
  • this surface structure can also be achieved with a single laser beam by self-organization of the azo molecules where the direction of the grating vector of the grating depends on the polarization of the laser beam.
  • the surface networks can then be erased by using a laser beam with a circular polarization or by heat treatment.
  • the Applicant has validated the concept of the present invention by performing ambient turbulence conditions and recording the images of the surface grating created by a laser beam after propagation of said beam in the ambient environment disturbed by the turbulence created, in order to highlight the correlation between the turbulences and the structural modifications made at the surface of the film of organic photochromic material.
  • the parameters of the turbulent medium necessary for the evaluation of the present invention were first evaluated.
  • the turbulence produced is representative of atmospheric turbulence under real conditions.
  • the air flow generated presents velocities v modified by a grid having steps of different widths, said speeds being between 12.5 m / s and 7.6 m / s.
  • the intensities and the angles of arrival of the beam F 1 were evaluated in two successive stages.
  • the samples of the polymer films are made from a synthesis of a novel azobenzene-derived chromophoric monomer containing heterocyclic sulfonamide compounds.
  • the methacrylate monomer contains aliphatic spacers of different lengths between the methacrylic and chromophoric groups.
  • the monomers in question are polymerized to obtain homopolymers and copolymerized with 2-methylacrylate and 2-ethylhexyl acrylate to obtain copolymers containing several percentages of chromophoric units.
  • the thin films are prepared by spin coating from a solution diluted with THF and a concentration of 50 mg / ml.
  • the thickness of the film measured with a dektak-6M profiler is 1 ⁇ m.
  • the laser beam is first made parallel by an afocal system of two lenses (Kepler).
  • the diameter of the beam is 5.1 mm.
  • a diaphragm makes it possible to choose the central region of the wavefront and to generate a quasi-plane wave with a diameter of 3.1 mm.
  • the diameter of the diaphragm is chosen to select a region of the laser beam as the best compromise between a region of the plane wavefront and an intensity signal allowing a reasonable value of the signal-to-noise ratio.
  • the laser beam passes into the turbulent medium over a distance of 1 m and is fed onto the azopolymer thin film 2 of the present invention with a power between 100 and 300 mW / cm 2 .
  • An optical system makes it possible to collect the intensity of the diffraction order and to control over time the unfolding of the holographic inscription.
  • the diffraction intensity reaches saturation, the laser beam is cut off and the surface of the thin film on which the process is inscribed is characterized by imaging, several types of imaging may be employed.
  • the Figures 2a and 2b relate respectively to the gratings inscribed in the photochromic material film and to the diffraction images taken by the CCD camera and comprising spots given by the diffraction order, for the four turbulence conditions T1, T2, T3 and T4.
  • the spacing of the spots gives the spatial frequency of the self-organized network present on the surface of the sample.
  • the evolution of this diffraction pattern follows a power law in time until saturation of the phenomenon and depends on the power of the laser beam.
  • the change in the propagation conditions of the laser beam induces additional spatial frequencies disturbing the regularity of the pitch of the beam. These additional frequencies are identified by the appearance of intensities outside the two spots.
  • the spatial information on the direction of the additional spatial frequency wave vectors is given by their X and Y coordinates with respect to the two spots and by the intensities surface. An important surface characterizes a significant diffusion of the spatial frequencies and thus a greater turbulence.
  • the figure 2a shows in particular that turbulence creates sub-structures with steps of As sub-networks of 400 ⁇ 30 nm.
  • the polarization variation induced by the turbulence is measured by the horizontal direction of the two diffraction spots.
  • the laser beam has a linear and horizontal polarization represented by the two diffraction spots on a horizontal plane.
  • the diffraction images can be analyzed by the Radon transform, the effective calculation of the discrete Radon transform involving numerically interpolations of the image to be processed, in order to estimate it in a rotating reference, as illustrated in FIG. figure 3 .
  • the representations of the intensity as a function of the angle ⁇ as illustrated in figure 5a and 5b show the presence of a central peak representing the projection of the 2 spots given by the diffraction of the beam on the network without turbulence.
  • the peaks outside the central peak represent the contributions of the additional spatial frequencies provided by the movement of the laser beam disturbed by the turbulence.
  • the inscribed network is characterized by atomic force microscopy in contact mode.
  • the figure 6a gives the topography of the film surface after registration by the laser for the four examples of turbulence conditions T1, T2, T3 and T4.
  • the surface self-organizes under the effect of the laser to create a diffraction grating.
  • Topographies obtained for different turbulence conditions are considered as random Gaussian surfaces.
  • the observed structures are first processed by an algorithm of the present invention to recognize, quantify and analyze the turbulence-modified surface sub-areas.
  • Fractal and multifractal analysis techniques are then used to recover the scale invariance properties of all disturbed sub-zones applicable to the measurement of the turbulence effect on a laser beam.
  • the image is first and foremost integrated.
  • a Gaussian low-pass filter is used to attenuate the high spatial frequencies.
  • the average filter each points with its neighbors. It has the effect of suppressing high frequencies and increasing the low frequency signal.
  • the filter makes it possible to reduce the contrast of the zones where there is no network and makes it possible to mark the zones where the network has disappeared because of the perturbation on the laser which creates the network by the turbulence.
  • This operation amounts to drawing the contours of these Gaussian surfaces by thresholding the height of the topography with a fixed spacing ⁇ between the height of the successive contours to obtain that the sub-zones do not contain periodic structures, as illustrated in FIG. figure 6b .
  • the threshold is determined by the empirical value I Avg / 6 where I moy is the average value of the roughness of the surface.
  • Contour drawings consist of non-intersecting, closed lines in the plane that connect the points of equal heights. In this way, for each configuration of random surfaces given by the different turbid media, a configuration of a set of contours is assigned.
  • the image is divided into boxes of equal length s.
  • a line of least squares is used to adjust the box to the data (representing the trend in the box).
  • the x and y coordinates of the line segments are calculated.
  • Fluctuation trends are removed from the image by subtracting the local trend.
  • the effective value of the fluctuations of the image are then calculated.
  • the calculation is repeated for all the lengths of the boxes to characterize the local fluctuations according to the size of the boxes.
  • the figure 7 illustrates the evolution of the coefficient ⁇ (q) as a function of the number q for the different turbulences T1, T2, T3 and T4. Curves that are not linear reflect a mutltifractal behavior.
  • h (q) describes the scale behavior of boxes with large fluctuations.
  • h (q) describes the scale behavior of segments with small fluctuations.
  • T1 turbulence induces little scale behavior for small fluctuations. The trend is reversed for other turbulence conditions.
  • the width of the spectra f ( ⁇ ) measures the level of multifractality of the signal. The larger the width, the more the signal is multifractal.
  • the figure 8 illustrates these spectra for the different turbulences T1, T2, T3 and T4.
  • C n The value of C n can be evaluated according to the maximum of f ( ⁇ ).
  • the widths are very different depending on the turbulence and can be seen that with high turbulence the signal contains more multifractal fluctuations.
  • the figure 9 thus illustrates the case of low turbulence and the case of strong turbulence.
  • the method of the present invention also makes it possible to calculate the gap that evaluates the homogeneity of the distribution of the light intensities and thus to know how the object fills the space.
  • lacunarity makes it possible to measure the distribution of holes in the image: if it contains large contrasts (on the scale), then the lacunarity of the image will be high. Conversely, if the image is homogeneous, then its gap will be small.
  • Lacunarity is measured using a sliding box algorithm.
  • the algorithm uses a box of size s x s with its width placed in the upper left corner of the image. The number of pixels of the signal in the box is counted.
  • the box is dragged one pixel to the right and the number of pixels in the box is counted so on. The process is repeated until the box reaches the lower right corner of the image.
  • N (s) ( ⁇ -s + 1) 2 .
  • the width s of the box is chosen so that the numerical value is equal to 2, 4, 8, ..., /.
  • Lacunarity L (s) for a box size is defined by: N (s). Q 2 / Q 1 2
  • p (i, ⁇ ) is the number of pixels of the signal in the ith box and i is between 1 and N (s).
  • a low value of b gives a hyperbole with a broad concavity (strong concavity).
  • a low value of b describes an object with a large distribution of empty gaps.
  • the figure 10 shows the calculation of the gaps calculated from the contrasting images where the contours corresponding to the threshold of contar images ( Figure 6c ).
  • the calculation shows the different textures of the images for different turbulence conditions.
  • the values obtained on the examples of structures shown in figure 6a gives for the value of the coefficient of parabola b as a function of C n : b Cn 2.61 9.4110 -13 8.62 1.47.10 -12 10.94 3.77.10 -12 19.11 4.12.10 -12
  • the curvature of the hyperbola becomes stronger as the turbulence decreases involving a larger gap between the dark areas of the images of the Figure 6c .
  • the beam fluctuates spatially and loses its coherence reducing the structured areas on the azopolymer film. These unstructured areas represent the integration over time of all spatial fluctuations.
  • the method of the present invention can be implemented in a unit 10 comprising a laser source 20, a support 30 covered with a film of photochromic material 31 and a CCD camera type detector 40. This unit is connected to a device treatment 50.
  • the device may comprise a set of units 10, 11,..., 1 N connected to a central processing device 60 and making it possible to collect different turbulence information in different locations via the center. treatment 60.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Electrochromic Elements, Electrophoresis, Or Variable Reflection Or Absorption Elements (AREA)

Description

  • Le domaine de l'invention est celui des procédés permettant de déceler les turbulences d'un milieu, exploitant les perturbations générées sur la propagation d'un faisceau laser et celui des dispositifs compacts mettant en oeuvre un tel procédé.
  • Ce type de procédé peut aussi bien être exploité pour mesurer des taux de turbulence (plus ou moins forte dans un environnement donné) que pour mesurer les effets de turbulence sur la propagation d'un faisceau laser.
  • En effet, l'analyse de turbulence peut trouver de nombreuses applications, via des moyens compacts et discrets, notamment dans le domaine de l'environnement, de détection d'élévations de températures, de mouvements d'air dans une ville par exemple.
  • Il peut également être très intéressant de mesurer les modifications engendrées sur la propagation d'un faisceau laser en raison de turbulences ambiantes. La propagation de faisceaux laser à travers un milieu turbulent est notamment d'une grande importance dans des applications telles que les communications optiques en espace libre (mer, champ de bataille, ...), la détection de cibles, l'imagerie et la détection à longue distance. Dans beaucoup de ces applications, les lasers avec leurs propriétés de cohérence du champ électromagnétique sont utilisés.
  • Les champs cohérents qui se propagent dans des milieux aléatoires tels que la turbulence atmosphérique, sont sujets à des distorsions et scintillations spatiales et temporelles aléatoires des intensités optiques qui peuvent entraîner une dégradation considérable dans la performance du système. Ces scintillations dues à des changements rapides et continus de températures et de pression proches de la surface du sol y sont plus sévères et plus prononcées et induisent des distorsions de la structure de phase du champ d'onde par la turbulence entraînant l'onde à interférer avec elle-même.
  • Ces changements dans les données et paramètres du faisceau apportent des limitations fondamentales dans le développement des systèmes optiques de communication en espace libre. Un des effets significatifs de la turbulence atmosphérique est la fluctuation de la direction de la propagation du front d'onde référencé comme la fluctuation de l'angle d'arrivée ou le tilt du front d'onde. Cet angle est mesuré par différentes méthodes.
  • Une des méthodes consiste à faire simultanément des mesures différentielles des angles d'arrivée et des fluctuations d'intensité du faisceau laser à partir d'une référence en utilisant une caméra CCD puisque les mesures de fluctuations spatiales sont très sensibles à la longueur d'onde du trajet parcouru par le faisceau laser (A. Consortini et al., Optics Communications, 216, 19-23, (2003)). Cette fluctuation dépend cubiquement de la longueur du trajet optique parcouru.
  • D'autres méthodes consistent à utiliser une structuration du faisceau laser soit avec des franges d'interférence et en analysant leur visibilité dans le cas d'un type d'anémométrie laser Doppler (H. J. Pjeider et al., J. Opt. Soc. Am., Vol. 70(2), 167-170, (1980)) ou bien en utilisant le déplacement des franges de Moiré (S. Rasouli et al., Opt. Lett., 31 (22), 3276-3278, (2006)) résultant de la superposition avec le même pas, d'une image de franges de référence et d'une image de franges résultant de la turbulence. Les structures de Moiré donnent les fluctuations de l'angle d'arrivée.
  • Dans les deux cas les méthodes sont difficiles à mettre en oeuvre puisqu'il faut créer les interférences pour l'anémométrie doppler ou bien installer deux réseaux d'amplitude et s'assurer de leur stabilité initiale pour la technique Moiré. Le tilt du front d'onde est lui analysé en utilisant des techniques polynomiales à partir d'une méthode d'analyse de type shark-hartmann (A. V. Sergryev et al., Appl. Opt., 50(20), 3519-3528, (2011)) et nécessite un système d'analyse du front d'onde du laser complexe et coûteux. Les fluctuations de l'angle d'arrivée ou le tilt du front d'onde sont utilisées pour mesurer les paramètres caractéristiques de la turbulence atmosphérique. Le paramètre de structure de l'indice de réfraction ou Cn 2 est une mesure de la force de la turbulence optique le long du trajet de propagation du faisceau laser. Dans la propagation laser en espace libre la mesure des fluctuations de l'angle d'arrivée est une étape de base dans l'étude de la turbulence atmosphérique.
  • Dans ce contexte, la présente demande a pour objet un procédé facile de mise en oeuvre, ouvrant la voie à des dispositifs d'évaluation de la turbulence d'un milieu, pouvant notamment être lié à des phénomènes d'échauffement de l'air ambiant, faciles de mise en oeuvre et pouvant être déployés massivement.
  • Plus précisément, la présente invention a pour objet un procédé d'évaluation de paramètres de turbulence d'un milieu caractérisé en ce qu'il comprend les étapes suivantes :
    • l'illumination par un faisceau laser se propageant dans ledit milieu, d'un film de matériau organique photochromique pouvant comprendre des molécules d'azobenzène, ou de spyropyrane, ou de stylbène, ou de coumarine, ou de fulgide, ladite illumination créant un réseau de surface dans ledit film comportant des striures dont les déformations sont fonction des paramètres de turbulence ;
    • une étape d'acquisition d'images comportant des zones élémentaires de structurations liées auxdites striures, plus ou moins réduites et représentatives des paramètres de turbulence dudit milieu ;
    • une étape de traitement desdites images permettant d'évaluer les paramètres de turbulence dudit milieu.
  • Selon une variante de l'invention, le procédé comprend une étape d'acquisition d'images de diffraction du faisceau laser par le réseau de surface créé, lesdites images présentant des spots de diffraction dont les caractéristiques sont dépendantes des paramètres de turbulence.
  • Selon une variante de l'invention, l'étape de traitement desdites images comprend :
    • une opération de transformée de Radon desdites images de manière à générer des images dérivées par calcul d'images dans un repère en rotation générant depuis des zones de diffraction appartenant auxdites images, un ensemble de pics appartenant aux images dérivées,
    • une comparaison dudit ensemble de pics avec un ensemble de pics de référence de conditions de turbulence, de manière à estimer le niveau de turbulence dudit milieu.
  • Selon une variante de l'invention, le procédé comprend :
    • une étape d'acquisition d'images réalisée par microscopie AFM ;
    • une étape de traitement desdites images permettant de reconnaître, quantifier et analyser les zones élémentaires formées en surface dépendantes, des paramètres de turbulence et permettant d'obtenir des images dérivées ;
    • une étape d'analyse desdites images dérivées par analyse fractale ou multifractale, sur l'ensemble desdites zones élémentaires permettant de définir des propriétés d'invariance d'échelle et de déterminer les informations de turbulence dudit milieu.
  • Selon une variante de l'invention, l'étape de traitement desdites images comprend les opérations suivantes :
    • l'utilisation d'un filtre gaussien, moyennant chacun des points d'acquisition avec des points voisins et générant des images filtrées;
    • la construction d'un ensemble d'images dérivées par traitement desdites images filtrées, pouvant être opérée par la transformation desdites images filtrées en images filtrées binaires, puis marquage des zones exemptes de striures, de manière à augmenter le contraste entre lesdites zones élémentaires.
  • Selon une variante de l'invention, l'étape d'analyse fractale ou multifractale comprend les étapes suivantes :
    • la division de ladite image selon deux dimensions x et y, en segments N/s, définissant une partition en boîtes, N étant une dimension de l'image et le calcul de deux sommes cumulatives X(i) et Y(i) définissant un signal, telles que : X i = ∑ k = 1 i x k − < x > et Y i = ∑ k = 1 i y k − < y >
      Figure imgb0001
       avec <x> et <y> correspondant aux moyennes du signal en x et y ;
    • l'estimation dans chacune des boîtes, d'une surface définie par deux fonctions polynômiales (pv(i), qv(i)) selon les deux dimensions x et y, permettant d'approximer ledit signal et permettant de définir :
      • o une tendance dans chacune des boîtes et
      • o les variations par rapport à ladite tendance de manière à définir les paramètres Xs(i) = X(i) - pv(i) et Ys(i) = Y(i) - qv(i)
    • le calcul d'une fonction de fluctuation dépouillée desdites tendances définie par l'équation : F q s = 1 N ∑ i = 1 N s X s i Y s i q q
      Figure imgb0002
       avec q entier non nul et où : ∑ i = 1 N s X s i Y s i q ≈ s τ q
      Figure imgb0003
       permettant de définir :
      • o une loi de puissance entre la fonction Fq(s) et la taille (s) des boîtes et
      • o l'exposant d'échelle multifractale r(q) et τ (q) = qh(q) - 1;
    • la détermination d'une fonction de singularité α(q) = dτ(q) /dq et d'un spectre multifractal : f(α) = q(α).α - s[q(α)] présentant une valeur maximale f(α)max ;
    • la comparaison de la valeur f(α)max avec une table de calibration référençant des valeurs f(α)max St avec des indices de réfraction Cn et permettant de relier l'indice de réfraction Cn avec un paramètre de turbulence de référence.
  • L'invention a aussi pour objet un dispositif d'évaluation des paramètres de turbulence d'un milieu, caractérisé en ce qu'il comprend :
    • une unité de mesure comportant :
      • o un film de matériau organique photochromique pouvant comprendre des molécules d'azobenzène, ou de spyropyrane, ou de stylbène, ou de coumarine, ou de fulgide ;
      • o une source laser émettant un faisceau laser se propageant dans ledit milieu et interagissant avec ledit film de manière à créer un réseau de surface comportant des striures dont les déformations sont fonction des paramètres de turbulence ;
      • o un détecteur permettant l'acquisition d'images comportant des zones élémentaires de structurations liées auxdites striures, plus ou moins réduites et représentatives des paramètres de turbulence dudit milieu et générant des signaux relatifs auxdites images ;
    • un dispositif de réception et de traitement desdits signaux délivrant des informations d'évaluation des paramètres de turbulence dudit milieu.
  • Selon une variante de l'invention, le dispositif d'évaluation des paramètres de turbulence d'un milieu comprend un ensemble d'unités de mesure, lesdites unités de mesure étant en liaison avec ledit dispositif de réception et de traitement.
  • L'invention a encore pour objet l'utilisation du dispositif d'évaluation des paramètres de turbulence d'un milieu ambiant selon l'une d'invention, pour réaliser des mesures de mouvements d'air, par exemple en milieu urbain.
  • 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 dispositif permettant de valider le concept de la présente invention ;
    • les figures 2a et 2b illustrent des réseaux inscrits dans un film de matériau photochrome et les différentes images correspondantes de diffraction prises par une caméra CCD ;
    • la figure 3 illustre le principe de la transformée de Radon d'une image d'un réseau de diffraction inscrit dans un film de matériau photochromique ;
    • les figures 4a et 4b illustrent les images résultantes de la transformée de Radon de deux exemples d'images illustrées en figure 2a ;
    • les figures 5a et 5b illustrent les pics d'intensité issus du traitement des images illustrées en figure 4a et 4b et représentatifs du degré de turbulence du milieu analysé ;
    • les figures 6a, 6b et 6c illustrent respectivement des images prises par microscopie AFM d'un réseau de surface inscrit dans un film de matériau photochromique, les mêmes images filtrées présentant des zones définies par des contours, puis les mêmes images comportant ces zones, après marquage des intérieurs desdites zones.
    • la figure 7 illustre l'évolution de l'exposant de multifractalité en fonction du paramètre q pour différentes conditions de turbulence ;
    • la figure 8 illustre le spectre de multifractalité f(α) appelé spectre de singularités pour différentes conditions de turbulence ;
    • la figure 9 illustre l'évolution des positions maximales f(α)max en fonction d'indice de réfraction Cn ;
    • la figure 10 illustre l'évolution de la fonction de lacunarité pour différentes conditions de turbulence ;
    • la figure 11 illustre un premier exemple de dispositif utilisant le procédé de l'invention ;
    • la figure 12 illustre un second exemple de dispositif utilisant le procédé de l'invention.
  • L'invention va être décrite dans le cadre d'un film d'azobenzène, mais peut tout aussi bien utiliser d'autres types de matériaux organiques photochromiques. L'effet utilisé dans la présente invention, correspond à un enregistrement holographique dans lequel les franges se construisent d'elles-mêmes durant le processus d'auto-organisation moléculaire et impriment une moyenne statistique sur le long terme de l'effet du milieu turbulent sur la propagation laser. Les matériaux holographiques organiques ont l'avantage par rapport à leurs équivalents inorganiques d'être structurellement polyvalents, facilement fabricables et peu coûteux.
  • Ce type de matériau photochromique, tel un matériau azopolymère, est capable de s'auto-organiser en surface sous l'effet de la lumière d'un laser pour produire une structuration sous forme d'un réseau de diffraction, permettant de mesurer simultanément plusieurs caractéristiques de ladite turbulence.
  • Les molécules d'azobenzène sont étudiées depuis plus d'un demi-siècle pour leurs propriétés optiques particulières et, encore de nos jours, des phénomènes nouveaux liés à leur photoisomérisation sont découverts. Parmi ces phénomènes, la déformation de matière photo-induite dans les matériaux photochromiques contenant ces molécules est sans doute l'un des plus remarquables. Des cycles successifs d'isomérisation de la molécule, peut résulter un déplacement de matière, sur des distances de plusieurs centaines de nanomètres. Les molécules de la famille des azobenzènes sont caractérisées par deux cycles benzéniques liés par une double liaison entre deux atomes d'azote. Elles peuvent exister sous deux formes isomères : la forme trans (ou E) et la forme cis (ou Z). La transition de la forme trans vers la forme cis est excitée par absorption d'un photon et la transition inverse de la forme cis vers la forme trans peut être soit excitée optiquement soit activée thermiquement.
  • Il a par ailleurs été montré que les molécules azo du matériau photoisomérisent en répétant des isomérisations trans-cis-trans induites par l'absorption de la lumière, de l'orientation moléculaire consécutive se faisant perpendiculairement à la direction de la polarisation de la lumière.
  • Cet effet peut être observé sur un film mince en polymère contenant des azobenzènes capables de fabriquer un réseau topographique sur la surface du film après exposition par une structure d'interférences créée par deux faisceaux laser.
  • Cependant le Demandeur est parti du principe que cette structuration en surface peut se réaliser aussi avec un seul faisceau laser par auto-organisation des molécules azo où la direction du vecteur d'onde du réseau dépend de la polarisation du faisceau laser.
  • Les réseaux de surface peuvent ensuite être effacés par utilisation d'un faisceau laser avec une polarisation circulaire ou par traitement thermique.
  • Le Demandeur a validé le concept de la présente invention en réalisant des conditions de turbulences en milieu ambiant et en enregistrant les images du réseau de surface créé par un faisceau laser après propagation dudit faisceau dans le milieu ambiant perturbé par les turbulences créées, et ce afin de mettre en évidence la corrélation entre les turbulences et les modifications de structurations réalisées au niveau de la surface du film de matériau organique photochromique.
  • Les paramètres du milieu turbulent nécessaire à l'évaluation de la présente invention ont d'abord été évalués. La turbulence produite est représentative d'une turbulence atmosphérique en conditions réelles. Le montage expérimental illustré en figure 1, consiste en un faisceau laser d'un laser 1, (solide pompé par diode) à λ = 473 nm avec une puissance de 5 mW se propageant au travers d'un milieu turbulent de longueur L = 1 m produit par au moins une plaque chauffante électrique (deux plaques P1 et P2 sont représentées en figure 1) à la température de 98 ° C et générant un flux d'air chaud à la température de 42°C.
  • Le flux d'air généré présente des vitesses v modifiées par une grille ayant des pas de largeurs différentes, lesdites vitesses étant comprises entre 12.5 m/s et 7.6 m/s.
  • Les intensités et les angles d'arrivée du faisceau F1 ont été évalués en deux étapes successives.
  • On utilise pour cela une caméra CCD 3, pour mesurer l'intensité totale de la distribution de lumière détectée et les coordonnées x et y du centre du faisceau, après passage en transmission du film de matériau photochromique 2. La surface du détecteur est de 8.4 mm x 6.2 mm. Les images peuvent être enregistrées et traitées via une unité 4.
    1. 1) Les angles d'arrivée sont obtenus en mesurant les déplacements latéraux du faisceau sur la caméra par comparaison avec un faisceau de référence sans turbulence. Les angles sont comptés entre un diaphragme et une position du faisceau sur la caméra. La distance entre le diaphragme et la caméra est une zone de propagation non turbulente et est suffisante pour une bonne résolution de la mesure de déplacement.
      Les valeurs des angles d'arrivée sont données par : α x = δx / d et α y = δy / d .
      Figure imgb0004
      Le Demandeur a mesuré l'angle d'arrivée pour différentes conditions de turbulence.
      On se place dans la limite d'approximation de l'optique géométrique pour L « /o 2/λ où /o est l'échelle intérieure de la turbulence et λ la longueur d'onde du faisceau laser.
      La valeur efficace des fluctuations de l'angle d'arrivée est donnée par : < α x 2 > = < α y 2 > = 3.28 C n 2 L l o − 1 3
      Figure imgb0005
       où Cn est la constante de l'indice de réfraction du milieu (L. C. Andrews, 'Field guide to Atmospheric Optics', SPIE press book, (2004)).
    2. 2) Les fluctuations d'intensité du faisceau sont mesurées au cours du temps avec une photodiode et sont définies par : σ I 2 = < I 2 > − < I 2 > < I 2 > = 12.8 C n 2 L 3 ℓ o − 7 3 = 1.23 C n 2 k 7 6 L 11 / 6
      Figure imgb0006
       où L est le trajet du faisceau laser. On peut donc calculer l'échelle intérieure de la turbulence par : ℓ o = 1.97 < α x 2 > σ I 2
      Figure imgb0007
  • Les échantillons des films en polymère sont faits à partir d'une synthèse d'un nouveau monomère chromophorique dérivé d'azobenzène contenant des composés sulfonamide hétérocyclique. Le monomère de type methacrylate contient des espaceurs aliphatiques de différentes longueurs entre les groupes méthacrycliques et chromophoriques. Les monomères en question sont polymérisés pour obtenir des homopolymères et copolymérisés avec du 2-methylacrylate et du 2-ethylhexyl acrylate pour obtenir des copolymères contenant plusieurs pourcentages d'unités chromophoriques.
  • Les films minces sont préparés par dépôt à la tournette à partir d'une solution diluée avec du THF et une concentration de 50 mg/ml. L'épaisseur du film mesurée avec un profileur dektak-6M est de 1 µm.
  • L'absorbance maximale à λ = 488 nm est de 1.9. L'absorbance à la longueur d'onde d'expérience λ = 476 nm est de 1.6.
  • Le faisceau laser est d'abord rendu parallèle par un système afocal de deux lentilles (Kepler). Le diamètre du faisceau est de 5.1 mm. Un diaphragme permet de choisir la région centrale du front d'onde et de générer une onde quasi-plane de diamètre 3.1 mm. Le diamètre du diaphragme est choisi pour sélectionner une région du faisceau laser comme le meilleur compromis entre une région du front d'onde plane et un signal d'intensité permettant une valeur raisonnable du ratio signal sur bruit.
  • Le faisceau laser passe dans le milieu turbulent sur une distance de 1 m et est envoyé sur le film mince d'azopolymère 2, de la présente invention avec une puissance entre 100 et 300 mW/cm2.
  • Un système optique permet de collecter l'intensité de l'ordre de diffraction et de contrôler dans le temps le déroulement de l'inscription holographique. Lorsque l'intensité de diffraction arrive à saturation, le faisceau laser est coupé et la surface du film mince sur lequel le processus est inscrit est caractérisé par imagerie, plusieurs types d'imagerie pouvant être employés.
  • Les expériences ont été menées pour quatre exemples de conditions de turbulence croissante T1, T2, T3 et T4, générant des topographies de surface différentes, mais dans lesquelles, la surface s'autoorganise sous l'effet du laser pour créer un réseau de diffraction.
    • Premier mode d'analyse par analyse semi-fine :
  • On observe avec une caméra de type CCD, la diffraction de la structure auto-organisée créée en surface et obtenue par le faisceau lui-même. Elle permet de suivre au cours du temps la construction de la structure en surface. Pour un faisceau se propageant sans turbulence le réseau forme un alignement régulier de striures à la surface avec un pas du réseau A = 2λ/sinθ où θ = 90° pour une incidence du faisceau de 90°.
  • Les figures 2a et 2b sont relatives respectivement aux réseaux inscrits dans le film de matériau photochromique et aux images de diffraction prises par la caméra CCD et comportant des spots donnés par l'ordre de diffraction, pour les quatre conditions de turbulence T1, T2, T3 et T4.
  • L'écartement des spots donne la fréquence spatiale du réseau auto-organisé présent sur la surface de l'échantillon. L'évolution de cette figure de diffraction suit une loi de puissance dans le temps jusqu'à saturation du phénomène et dépend de la puissance du faisceau laser. Le changement des conditions de propagation du faisceau laser induit des fréquences spatiales additionnelles perturbant la régularité du pas du faisceau. Ces fréquences additionnelles se repèrent par l'apparition d'intensités en dehors des deux spots. L'information spatiale sur la direction des vecteurs d'onde des fréquences spatiales additionnelles est donnée par leurs coordonnées X et Y par rapport aux deux spots et par la surface des intensités. Une surface importante caractérise une diffusion importante des fréquences spatiales et donc une plus grande turbulence.
  • Avec la modification du milieu de propagation des fréquences additionnelles se créent au centre et en dehors des deux spots correspondant au réseau auto-organisé de référence. Ces fréquences spatiales sont plus basses impliquant une variation de la surface sur de grande distance et correspondent à la création de larges structures. L'apparition de fréquences sur le profil des intensités de diffraction en dehors des deux spots implique la création de sous-réseaux aux fréquences plus grandes que le réseau principal.
  • La figure 2a montre notamment que la turbulence crée des sous-structures avec des pas des sous-réseaux As de 400 ± 30 nm.
  • En monitorant la diffraction créée, on peut suivre l'évolution des sous-réseaux créés ainsi que des zones non-uniformes qui apparaissent.
  • La variation de polarisation induite par la turbulence est mesurée par la direction horizontale des deux spots de diffraction. Le faisceau laser a une polarisation linéaire et horizontale représentée par les deux spots de diffraction sur un plan horizontal.
  • Par exemple, la composante linéaire de polarisation change d'un angle θ = 7.1 ° pour l'exemple de la turbulence T4.
  • Avantageusement, les images de diffraction peuvent être analysées par la transformée de Radon, le calcul effectif de la transformée de Radon discrète impliquant numériquement des interpolations de l'image à traiter, afin de l'estimer dans un repère en rotation, comme illustré en figure 3.
  • Pour cela, on calcule la projection de l'image le long d'angles donnés. Le résultat est la somme des intensités des pixels dans chaque direction, c'est-à-dire l'intégrale de chacune des lignes. On obtient une nouvelle image R(ρ,θ) et on fait une coupe pour t = 0°. On obtient par exemple les images dérivées des figures 4a et 4b correspondant à une moyenne turbulence et à une forte turbulence.
  • Les représentations de l'intensité en fonction de l'angle θ telles qu'illustrées en figure 5a et 5b montrent la présence d'un pic central représentant la projection des 2 spots donnés par la diffraction du faisceau sur le réseau sans turbulence. Les pics en dehors du pic central représentent les contributions des fréquences spatiales additionnelles apportées par le mouvement du faisceau laser perturbé par la turbulence. En comparant les pics qui se forment de part et d'autre du pic central on peut qualitativement connaitre instantanément avec la caméra le degré de turbulence induit durant le trajet du faisceau laser
    • 2) Second mode d'analyse fine par imagerie AFM : 2.1 : Etape de traitement d'images :
  • On caractérise le réseau inscrit par microscopie à force atomique en mode contact. La figure 6a donne la topographie de la surface du film après inscription par le laser pour les quatre exemples de conditions de turbulence T1, T2, T3 et T4. La surface s'auto-organise sous l'effet du laser pour créer un réseau de diffraction.
  • Pour une faible condition de turbulence, la surface présente un réseau auto-organisé dont le pas du réseau A = 800 ± 30 nm est perturbé par les fluctuations aléatoires de l'angle d'arrivée du faisceau. Le pas du réseau donné par la théorie de la diffraction du second ordre correspond à la condition d'angle du faisceau incident A = 2λ/ sinθ.
  • On remarque que certaines zones de la structure en surface ont une structuration réduite. On considère les topographies obtenues pour différentes conditions de turbulence comme des surfaces gaussiennes aléatoires.
  • Les structurations observées sont d'abord traitées par un algorithme de la présente invention pour reconnaître, quantifier et analyser les sous-zones en surface modifiées par la turbulence.
  • On applique ensuite des techniques d'analyse fractales et multifractales pour retrouver les propriétés d'invariance d'échelles sur l'ensemble des sous-zones perturbées applicables à la mesure de l'effet de turbulence sur un faisceau laser.
  • Ces techniques permettent d'avoir une description précise de la complexité des sous-zones désorganisées.
  • Plus précisément, l'image est en premier lieu intégrée. On utilise pour cela un filtre passe-bas gaussien qui a pour but d'atténuer les hautes fréquences spatiales. Le filtre moyenne chacun des points avec ses voisins. Il a pour effet de supprimer les hautes fréquences et d'augmenter le signal de basse fréquence. Dans le cas présent, le filtre permet de réduire le contraste des zones où il y a absence de réseau et permet de marquer les zones où le réseau a disparu à cause de la perturbation sur le laser qui crée le réseau par la turbulence.
  • Cette opération revient à dessiner les contours de ces surfaces gaussiennes en seuillant la hauteur de la topographie avec un espacement fixe Δ entre la hauteur des contours successifs pour obtenir que les sous-zones ne contenant pas de structurations périodiques, comme illustré en figure 6b. Le seuil est déterminé par la valeur empirique I moy / 6
    Figure imgb0008
     où Imoy est la valeur moyenne de la rugosité de la surface. Les dessins des contours consistent en lignes non intersectantes fermées dans le plan qui connectent les points d'égales hauteurs. De cette façon, pour chaque configuration de surfaces aléatoires données par les différents milieux turbides, on assigne une configuration d'un ensemble de contours.
  • Ces contours sont ensuite remplis, comme illustré en figure 6c.
    • 2.2 : Etape de calcul sur les images :
  • A partir de l'image comportant des zones définis par des contours (tels qu'illustrés en figure 6c), contenant des taches plus ou moins importantes suivant la turbulence, on traite ces images par analyse multifractale. On cherche à mettre en évidence les processus d'auto-similarités dans les images traitées pour avoir les zones dans lesquelles le réseau de surface obtenu par auto-organisation moléculaire a disparu.
    • 2.2.1 : Première étape :
  • L'image est divisée en boîtes d'égales longueurs s. Dans chacune de ces boîtes de longueur s, une ligne de moindres carrés est utilisée pour ajuster la boîte aux données (représentant la tendance dans la boîte).
  • Les coordonnées x et y des segments de lignes sont calculés.
  • Les tendances des fluctuations sont supprimées de l'image en soustrayant la tendance locale.
  • La valeur efficace des fluctuations de l'image sont alors calculées. Le calcul est répété pour toutes les longueurs des boîtes pour caractériser les fluctuations locales en fonction de la taille des boîtes.
  • Plus précisément, la première étape consiste à diviser les profils en x et y de l'image en Ns=[N/s] segments disjoints respectivement et calculer la somme cumulative X(i) et Y(i) telle que : X i = ∑ k = 1 i x k − < x > et Y i = ∑ k = 1 i y k − < y > .
    Figure imgb0009
  • Dans chaque segment de dimension s, on approxime (« fit ») les séries du signal intégré en utilisant une fonction polynomiale pv(i) pour {xi} et qv(i) pour {yi}.
    • 2.2.2 : Seconde étape
  • On calcule la fonction de fluctuation dépouillée des tendances en soustrayant la tendance locale dans chaque boîte et on calcule la fluctuation du moment de la racine qième : F q s = 1 N ∑ i = 1 N s X s i Y s i q q
    Figure imgb0010
     où : ∑ i = 1 N s X s i Y s i q ≈ s τ q
    Figure imgb0011
     avec Xs(i) = X(i) - pv(i), Ys(i) = Y(i) - qv(i) et q entier non nul.
    • 2.2.3 : Troisième étape :
  • Une relation de loi de puissance entre la fonction Fq(s) et la taille de la boîte s indique la présence d'échelle. On obtient de cette manière les exposants de Hurst généralisés : τ (q). Ils sont liés aux différents niveaux d'amplification des petites fluctuations dans les données.
  • La dépendance de τ (q) est une fonction monotone décroissance de q pour un signal stationnaire et reflète de façon simple le fait que les faibles fluctuations surviennent plus souvent dans le signal que les plus grosses. avec τ q = qh q − 1
    Figure imgb0012
  • La figure 7 illustre l'évolution du coefficient τ (q) en fonction du nombre q pour les différentes turbulences T1, T2, T3 et T4. Les courbes qui ne sont pas linéaires reflètent un comportement mutltifractal.
  • Pour les variations positives de q, h(q) décrit le comportement d'échelle des boîtes avec de larges fluctuations.
  • Pour une forte turbulence la présente invention permet de constater un exposant d'échelle important pour les larges fluctuations. Au contraire pour les valeurs de q négatives, h(q) décrit le comportement d'échelle des segments avec de petites fluctuations.
  • On voit que la turbulence T1 induit peu de comportement d'échelle pour les faibles fluctuations. La tendance est inversée pour les autres conditions de turbulence.
  • On remarque qu'en réduisant la turbulence (turbulence T4), le signal devient de plus en plus monofractal.
    • 2.2.4. Quatrième étape :
  • On peut utiliser la fonction des singularités obtenue par une transformée de Legendre : α q = dτ q / dq
    Figure imgb0013
  • On peut aussi utiliser un spectre de multifractalité : f α = q α α − s q α
    Figure imgb0014
  • Cette dernière approche donne une représentation utile de la multifractalité.
  • La largeur des spectres f(α) mesure le niveau de multifractalité du signal. Plus la largeur est importante plus le signal est multifractal. La figure 8 illustre ces spectres pour les différentes turbulences T1, T2, T3 et T4.
  • Le maximum de la fonction f(α) en fonction de α est tracé en figure 8. On note deux zones d'évolution linéaire du maximum en fonction de la force de la turbulence.
  • La valeur de Cn peut être évaluée en fonction du maximum de f(α).
  • La largeur de f(α) à f = 0.9, caractérise le nombre de singularités multifractales c'est-à-dire que le signal a de multiples lois d'échelles et est fortement perturbé. Les largeurs sont très différentes en fonction de la turbulence et permettent de voir qu'avec une turbulence élevée le signal contient plus de fluctuations multifractales.
  • Cette dernière mesure caractérise l'évolution des sous-structures en surface de l'azopolymère qui deviennent moins régulière avec l'augmentation de la turbulence. La valeur de f(α) est liée à Cn par la relation
    • maxHdff = -0.8.Cn +0.8 pour une faible turbulence et
    • maxHdff = -0.23 .Cn +0.26 pour une forte turbulence.
    • avec f(α)max = maxHdff.
  • La figure 9 illustre ainsi le cas d'une faible turbulence et le cas d'une forte turbulence.
  • Le procédé de la présente invention permet également de calculer la lacunarité qui évalue l'homogénéité de la répartition des intensités lumineuses et donc de savoir comment l'objet rempli l'espace.
  • L'étude de la lacunarité permet de mesurer la distribution de trous dans l'image : si celle- ci comporte de grands contrastes (à l'échelle), alors la lacunarité de l'image sera élevée. A l'inverse, si l'image est homogène, alors sa lacunarité sera faible.
  • La lacunarité est mesurée en utilisant un algorithme de boîte glissante.
  • L'algorithme utilise une boîte de taille s x s avec s sa largeur placée au coin supérieur gauche de l'image. Le nombre de pixels du signal dans la boîte est compté.
  • La boîte est glissée d'un pixel sur la droite et le nombre de pixels dans la boîte est compté ainsi de suite. Le processus est répété jusqu'à ce que la boîte atteigne le coin inférieur droit de l'image.
  • Si λ est la taille de l'image le nombre de boîtes de largeur s est
  • N(s) = (λ-s+l)2. La largeur s de la boîte est choisie de telle sorte que la valeur numérique est égale à 2, 4, 8, ..., /.
  • Pour chaque boîte de largeur s, les deux quantités suivantes sont mesurées : la somme du nombre de pixels du signal dans chaque boîte
    Q1 = ∑i p (i, s) et la somme quadratique Q1 = ∑i p (i, s) 2
  • La lacunarité L(s) pour une taille de boîte est définie par :
    N(s). Q2/Q1 2
  • p(i, λ) est le nombre de pixels du signal dans la iieme boîte et i est compris entre 1 et N(s).
  • La fonction de lacunarité est évaluée à partir des mesures de la figure 10 approximée (« fitée ») avec une fonction hyperbolique :
    • L(s)= b/sa + c, où b est relatif à la concavité de l'hyperbole.
  • Une faible valeur de b donne une hyperbole avec une large concavité (forte concavité). Une faible valeur de b décrit un objet avec une forte distribution d'interstices vides.
  • La figure 10 montre le calcul des lacunarités calculées à partir des images contrastées où les contours correspondant au seuil des images contar (figure 6c). Le calcul montre les différentes texturations des images pour différentes conditions de turbulence. Les valeurs obtenues sur les exemples de structurations montrés à la figure 6a donne pour la valeur du coefficient de la parabole b en fonction de Cn :
    b Cn
    2.61 9.4110-13
    8.62 1.47.10-12
    10.94 3.77.10-12
    19.11 4.12.10-12
  • La courbure de l'hyperbole devient de plus en plus forte au fur et à mesure que la turbulence diminue impliquant un plus grand espace entre les zones foncées des images de la figure 6c. Avec la turbulence, le faisceau fluctue spatialement et perd sa cohérence réduisant les zones structurées sur le film d'azopolymère. Ces zones non structurées représentent l'intégration dans le temps de toutes les fluctuations spatiales.
  • Selon un premier exemple de dispositif représenté en figure 11, le procédé de la présente invention peut être implémenté dans une unité 10 comprenant une source laser 20, un support 30 recouvert d'un film de matériau photochromique 31 et une détecteur de type caméra CCD, 40. Cette unité est reliée à un dispositif de traitement 50.
  • Selon un second exemple de l'invention, le dispositif peut comprendre un ensemble d'unités 10, 11,..., 1 N reliées à un dispositif de traitement central 60 et permettant de collecter différentes informations de turbulences en différents endroits via le centre de traitement 60.

Claims (9)

  1. Procédé d'évaluation de paramètres de turbulence d'un milieu comprenant les étapes suivantes :
    - l'illumination par un faisceau laser se propageant dans ledit milieu, d'un film de matériau organique photochromique pouvant comprendre des molécules d'azobenzène, ou de spyropyrane, ou de stylbène, ou de coumarine, ou de fulgide, ladite illumination créant un réseau de surface dans ledit film comportant des striures dont les déformations sont fonction des paramètres de turbulence ;
    - une étape d'acquisition d'images comportant des zones élémentaires de structurations liées auxdites striures, plus ou moins réduites et représentatives des paramètres de turbulence dudit milieu ;
    - une étape de traitement desdites images permettant d'évaluer les paramètres de turbulence dudit milieu.
  2. Procédé d'évaluation de paramètres de turbulence selon la revendication 1, caractérisé en ce qu'il comprend une étape d'acquisition d'images de diffraction du faisceau laser par le réseau de surface créé, lesdites images présentant des spots de diffraction dont les caractéristiques sont dépendantes des paramètres de turbulence.
  3. Procédé d'évaluation de paramètres de turbulence selon la revendication 2, caractérisé en ce que l'étape de traitement desdites images comprend :
    - une opération de transformée de Radon desdites images de manière à générer des images dérivées par calcul d'images dans un repère en rotation générant depuis des zones de diffraction appartenant auxdites images, un ensemble de pics appartenant aux images dérivées,
    - une comparaison dudit ensemble de pics avec un ensemble de pics de référence de conditions de turbulence, de manière à estimer le niveau de turbulence dudit milieu.
  4. Procédé d'évaluation de paramètres de turbulence selon la revendication 1, caractérisé en ce qu'il comprend :
    - une étape d'acquisition d'images réalisée par microscopie AFM ;
    - une étape de traitement desdites images permettant de reconnaître, quantifier et analyser les zones élémentaires formées en surface dépendantes des paramètres de turbulence et permettant d'obtenir des images dérivées ;
    - une étape d'analyse desdites images dérivées par analyse fracale ou multifractale, sur l'ensemble desdites zones élémentaires permettant de définir des propriétés d'invariance d'échelle et de déterminer les informations de turbulence dudit milieu.
  5. Procédé d'évaluation de paramètres de turbulence selon la revendication 4, caractérisé en ce que l'étape de traitement desdites images comprend les opérations suivantes :
    - l'utilisation d'un filtre gaussien, moyennant chacun des points d'acquisition aves des points voisins et générant des images filtrées;
    - la construction d'un ensemble d'images dérivées par traitement desdites images filtrées, pouvant être opérée par la transformation desdites images filtrées en images filtrées binaires, puis marquage des zones exemptes de striures, de manière à augmenter le contraste entre lesdites zones élémentaires.
  6. Procédé d'évaluation de paramètres de turbulence selon la revendication 5, caractérisé en ce que l'étape d'analyse fractale ou multifractale comprend les étapes suivantes :
    - la division de ladite image selon deux dimensions x et y, en segments N/s, définissant une partition en boîtes, N étant une dimension de l'image et le calcul de deux sommes cumulatives X(i) et Y(i) définissant un signal, telles que : X i = ∑ k = 1 i x k − < x > et Y i = ∑ k = 1 i y k − < y >
    Figure imgb0015
     avec <x> et <y> correspondant aux moyennes du signal en x et y.
    - l'estimation dans chacune des boîtes, d'une surface définie par deux fonctions polynômiales (pv(i), qv(i)) selon les deux dimensions, permettant d'approximer ledit signal et permettant de définir une tendance dans chacune des boîtes et les variations par rapport à ladite tendance de manière à définir les paramètres Xs(i) = X(i) - pv(i) et Ys(i) = Y(i) - qv(i) ;
    - le calcul d'une fonction de fluctuation dépouillée desdites tendances définie par l'équation : F q s = 1 N ∑ i = 1 N s X s i Y s i q q
    Figure imgb0016
     avec q entier non nul et où : ∑ i = 1 N s X s i Y s i q ≈ s τ q
    Figure imgb0017
     et permettant de définir une loi de puissance entre la fonction Fq(s) et la taille (s) des boîtes et de définir l'exposant d'échelle multifractale τ(q) ;
    - la détermination d'une fonction de singularité α(q) = dτ(q) /dq et d'un spectre multifractal : f(α) = q(α).α - s[q(α)] présentant une valeur maximale f(α)max ;
    - la comparaison de la valeur f(α)max avec une table de calibration référençant des valeurs f(α)max St avec des indices de réfraction Cn et permettant de relier l'indice de réfraction Cn avec un paramètre de turbulence de référence.
  7. Dispositif d'évaluation des paramètres de turbulence d'un milieu, comprenant :
    - une unité de mesure comportant :
    ∘ un film de matériau organique photochromique pouvant comprendre des molécules d'azobenzène, ou de spyropyrane, ou de stylbène, ou de coumarine, ou de fulgide ;
    o une source laser émettant un faisceau laser se propageant dans ledit milieu et interagissant avec ledit film de manière à créer un réseau de surface comportant des striures dont les déformations sont fonction des paramètres de turbulence ;
    o un détecteur permettant l'acquisition d'images comportant des zones élémentaires de structurations liées auxdites striures, plus ou moins réduites et représentatives des paramètres de turbulence dudit milieu et générant des signaux relatifs auxdites images ;
    - un dispositif de réception et de traitement desdits signaux délivrant des informations d'évaluation des paramètres de turbulence dudit milieu.
  8. Dispositif d'évaluation des paramètres de turbulence d'un milieu selon la revendication 7, caractérisé en ce qu'il comprend un ensemble d'unités de mesure, lesdites unités de mesure étant en liaison avec ledit dispositif de réception et de traitement.
  9. Utilisation du dispositif d'évaluation des paramètres de turbulence d'un milieu ambiant selon l'une des revendications 7 ou 8, pour réaliser des mesures de mouvements d'air, par exemple en milieu urbain.
EP13779888.0A 2012-10-24 2013-10-23 Procede de detection de conditions de turbulence utilisant l'interaction d'un faisceau laser avec un film mince photochromique et dispositif mettant en oeuvre ledit procede Not-in-force EP2912475B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1260113A FR2997189B1 (fr) 2012-10-24 2012-10-24 Procede de detection de conditions de turbulence utilisant l'interaction d'un faisceau laser avec un film mince photochromique et dispositif mettant en œuvre ledit procede
PCT/EP2013/072186 WO2014064159A1 (fr) 2012-10-24 2013-10-23 Procede de detection de conditions de turbulence utilisant l'interaction d'un faisceau laser avec un film mince photochromique et dispositif mettant en œuvre ledit procede

Publications (2)

Publication Number Publication Date
EP2912475A1 EP2912475A1 (fr) 2015-09-02
EP2912475B1 true EP2912475B1 (fr) 2016-09-07

Family

ID=47989072

Family Applications (1)

Application Number Title Priority Date Filing Date
EP13779888.0A Not-in-force EP2912475B1 (fr) 2012-10-24 2013-10-23 Procede de detection de conditions de turbulence utilisant l'interaction d'un faisceau laser avec un film mince photochromique et dispositif mettant en oeuvre ledit procede

Country Status (8)

Country Link
US (1) US9684011B2 (fr)
EP (1) EP2912475B1 (fr)
CN (1) CN104871008B (fr)
BR (1) BR112015009212A2 (fr)
DK (1) DK2912475T3 (fr)
ES (1) ES2611172T3 (fr)
FR (1) FR2997189B1 (fr)
WO (1) WO2014064159A1 (fr)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9228785B2 (en) 2010-05-04 2016-01-05 Alexander Poltorak Fractal heat transfer device
US11530985B2 (en) * 2015-08-19 2022-12-20 Nanoscope Technologies, LLC Cancer diagnosis by refractive index multifractality
US10830545B2 (en) 2016-07-12 2020-11-10 Fractal Heatsink Technologies, LLC System and method for maintaining efficiency of a heat sink
US10733708B2 (en) 2018-09-24 2020-08-04 United States Of America As Represented By Secretary Of The Navy Method for estimating turbulence using turbulence parameter as a focus parameter
EP3705906A3 (fr) * 2019-03-05 2020-12-09 Visteon Global Technologies, Inc. Système de détection et de télémétrie de la lumière à couches verticales multiples, assistance au stationnement automatique et détection et maintien de la voie de vision par ordinateur
US10929962B2 (en) 2019-06-19 2021-02-23 United States Of America As Represented By The Secretary Of The Navy Imposing a priori knowledge to estimate turbulence
CN111765929B (zh) * 2020-06-22 2021-10-15 中国科学院西安光学精密机械研究所 加注管道流量图像测量方法及测量装置
CN115755559B (zh) * 2022-12-23 2024-09-17 陕西师范大学 基于ifs分形图拼接的全息信息嵌入/恢复方法及系统
CN116933053A (zh) * 2023-07-18 2023-10-24 中国人民解放军国防科技大学 平流层多尺度扰动中能量传输和物质输送的识别方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7241569B2 (en) * 1993-01-18 2007-07-10 Olympus Corporation Method and a device for the evaluation of biopolymer fitness
US7268937B1 (en) * 2005-05-27 2007-09-11 United States Of America As Represented By The Secretary Of The Air Force Holographic wavefront sensor

Also Published As

Publication number Publication date
DK2912475T3 (da) 2017-01-02
CN104871008B (zh) 2017-09-29
FR2997189B1 (fr) 2015-01-02
FR2997189A1 (fr) 2014-04-25
ES2611172T3 (es) 2017-05-05
WO2014064159A1 (fr) 2014-05-01
EP2912475A1 (fr) 2015-09-02
CN104871008A (zh) 2015-08-26
US20150293140A1 (en) 2015-10-15
BR112015009212A2 (pt) 2017-07-04
US9684011B2 (en) 2017-06-20

Similar Documents

Publication Publication Date Title
EP2912475B1 (fr) Procede de detection de conditions de turbulence utilisant l&#39;interaction d&#39;un faisceau laser avec un film mince photochromique et dispositif mettant en oeuvre ledit procede
EP3343201B1 (fr) Procédé et dispositif de comptage et de caractérisation de particules dans un fluide en mouvement
EP3274694B1 (fr) Procédé de détermination de l&#39;état d&#39;une cellule
EP3735578B1 (fr) Procede et systeme d&#39;inspection optique d&#39;un substrat
FR2880472A1 (fr) Procede et systeme de diffusometrie
Xue et al. Detection of microplastics based on spatial heterodyne Raman spectroscopy
León-Rodríguez et al. High topographical accuracy by optical shot noise reduction in digital holographic microscopy
EP1157261B1 (fr) Procede et dispositif d&#39;analyse d&#39;un front d&#39;onde a grande dynamique
Song et al. Comparison of multichannel signal deconvolution algorithms in airborne LiDAR bathymetry based on wavelet transform
Zhang et al. Non-invasive estimation of the powder size distribution from a single speckle image
CN102494614A (zh) 一种高精度数字散斑相关测量方法
FR2825150A1 (fr) Dispositif de caracterisation de reseaux optiques et procede de fabrication de reseaux optiques avec une frequence spatiale predefinie
EP4022280B1 (fr) Méthode et dispositif de caractérisation optique de particules
Zhang et al. Efficient white light interferometry for wafer‑scale high aspect ratio structures
Li et al. Two-step iterative phase demodulation method for phase-shifting interferometry
Denisov et al. Method for certification monitoring of surface inhomogeneities of optics based on frequency analysis of the surface profile
FR3100335A1 (fr) Méthode et dispositif de caractérisation optique de particules
WO2013132072A1 (fr) Procede et appareil de mesure de la structure geometrique d&#39;un composant optique
EP3995991A2 (fr) Procédé de détection d&#39;un évènement cellulaire
Du et al. Numerical refocusing algorithm for in-line holography based on sparse regularization and physical constraints
Tompkins et al. Comparison of Methods for Wavelength Determination from Laser Speckle Patterns
FR3100332A1 (fr) Méthode et dispositif de caractérisation optique de particules
FR3042293A1 (fr)
Chen et al. Deep learning-driven recovery of oversaturated interferometric signals for continuous 3D morphology measurement of complex surfaces in OFC-TS-DFT systems
Stoykova et al. Compressed dynamic speckle sensing

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20150420

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAX Request for extension of the european patent (deleted)
RIC1 Information provided on ipc code assigned before grant

Ipc: G01P 5/00 20060101AFI20160203BHEP

Ipc: G03H 1/00 20060101ALI20160203BHEP

Ipc: G01P 5/26 20060101ALI20160203BHEP

Ipc: G06T 3/00 20060101ALI20160203BHEP

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

RIN1 Information on inventor provided before grant (corrected)

Inventor name: BARILLE, REGIS

INTG Intention to grant announced

Effective date: 20160323

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: UNIVERSITE D'ANGERS

Owner name: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE

Owner name: WROCLAW UNIVERSITY OF TECHNOLOGY

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 827327

Country of ref document: AT

Kind code of ref document: T

Effective date: 20161015

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602013011215

Country of ref document: DE

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 4

REG Reference to a national code

Ref country code: SE

Ref legal event code: TRGR

REG Reference to a national code

Ref country code: DK

Ref legal event code: T3

Effective date: 20161220

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20160907

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20161207

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 827327

Country of ref document: AT

Kind code of ref document: T

Effective date: 20160907

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161031

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20161208

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2611172

Country of ref document: ES

Kind code of ref document: T3

Effective date: 20170505

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170107

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170109

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20161207

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602013011215

Country of ref document: DE

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

REG Reference to a national code

Ref country code: IE

Ref legal event code: MM4A

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161031

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161031

26N No opposition filed

Effective date: 20170608

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161023

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 5

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161023

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20161031

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20131023

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 6

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20160907

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DK

Payment date: 20200922

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20201009

Year of fee payment: 8

Ref country code: IT

Payment date: 20201009

Year of fee payment: 8

Ref country code: ES

Payment date: 20201103

Year of fee payment: 8

Ref country code: FR

Payment date: 20201030

Year of fee payment: 8

Ref country code: GB

Payment date: 20201015

Year of fee payment: 8

Ref country code: SE

Payment date: 20201013

Year of fee payment: 8

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 602013011215

Country of ref document: DE

REG Reference to a national code

Ref country code: DK

Ref legal event code: EBP

Effective date: 20211031

REG Reference to a national code

Ref country code: SE

Ref legal event code: EUG

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20211023

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211024

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211023

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20220503

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211031

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211023

Ref country code: DK

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211031

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20230217

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20211024