Brevet FR3067764A1 — publié le 20181221

- 1 - Dispositif de contrôle actif du recollement d'un écoulement sur un profil L'invention concerne le domaine de la conception des véhicules et en particulier de l'aménagement des zones de surfaces en contact avec les fluides dans lesquels les véhicules circulent, appelées profils. Il est bien connu que tout véhicule qui se déplace dans un fluide, gaz ou liquide, très généralement newtonien comme l'air atmosphérique ou l'eau, induit par son déplacement une force de traînée qui va à l'encontre du déplacement du véhicule et qui est issue de l'interaction entre le profil du véhicule et le fluide dans lequel il circule. Cette force de tramée devient une véritable cause de surconsommation d'énergie et de perte de performances à vitesses relatives élevées, et à partir de 50km/h. De manière générale, la traînée est causée, lors de l'écoulement du fluide autour du profil, aussi appelé flux, par les frottements entre ces derniers ainsi que par les décollements de l'écoulement autour du profil, qui forment des tourbillons aussi appelés vortex, induisant des dépressions en aval de l'écoulement. Bien qu'il soit très difficile de réduire les frottements qui sont inhérents au mouvement relatif du profil et du fluide, il est possible d'adapter le profil pour réduire le plus possible les décollements du flux. Les travaux de diminution de la tramée se sont principalement axés sur la recherche de solutions passives, c'est-à-dire consistant à faire varier la forme du profil dans le but de garder l'écoulement collé au profil. - 2 - Cependant, ces solutions passives semblent aujourd'hui avoir atteint leurs limites. Plus récemment, l'intérêt s'est porté sur des solutions actives de réduction du décollement, faisant intervenir des contrôleurs divers, 5 comme le décrit Brunton and Noack (2015). Parmi ces contrôleurs, on retrouve généralement l'utilisation d'injecteurs de fluides, comme le décrit Volino (2003). Ces injecteurs débouchent à la surface du profil et y injectent un fluide, de même type que celui s'écoulant, possédant une 10 énergie assez élevée pour recoller l'écoulement sur le profil et éviter l'apparition de vortex. Toutefois se pose le problème de contrôler convenablement ces injecteurs afin d'obtenir un résultat fiable et prouvé de recollement, tout en économisant le plus possible l'injection de fluide qui 15 est consommatrice d'énergie. Une des solutions communément utilisées est d'appliquer un contrôle des injecteurs en boucle ouverte. La conception et l'application d'un tel contrôle est 20 simple, cependant il n'est pas robuste et le moindre changement du profil ou de son environnement peut provoquer des résultats inattendus et potentiellement dangereux. On trouve également dans la littérature des contrôles boucle fermée qui utilisent une mesure de taux de décollement de 25 l'écoulement en aval du profil comme entrée de boucle de contrôle. On pourra ici citer Tian et al. (2006). Ceux-ci sont dans certains cas basés sur les équations Navier Stokes mais leur complexité empêche une utilisation en temps réel. Le plus souvent, ils sont basés sur des modèles 30 linéaires ou sont des contrôleurs dits sans modèle. On peut citer, notamment, le PID, pour contrôleur Proportionnel Intégral Dérivé, qui est un contrôle très classique de l'automatique bien qu'il ne soit pas bien adapté pour - 3 - contrôler des relais. Les résultats pour le contrôle des écoulements détachés sont assez mitigés et sa robustesse est limitée. On trouve également depuis quelques années du contrôle fait par Machine Learning. Il s'agit de faire 5 apprendre à l'algorithme par de très nombreux tests quel est le meilleur contrôle à appliquer. Les résultats sont bons et une certaine robustesse est observée expérimentalement, cependant ce type d'algorithmes manque de preuves mathématiques et leurs limites sont difficiles à 10 déterminer. L'invention vient proposer une solution pour pallier les problèmes ci-dessus évoqués et recoller, avec un contrôle fin, efficace et robuste, l'écoulement du fluide 15 sur le profil. A cet effet, l'invention propose un dispositif de pénétration dans un fluide, comprenant au moins un profil, un injecteur de fluide en périphérie du profil, au moins un 20 moyen de mesure de décollement d'un écoulement le long du profil en aval de l'injecteur par rapport à un sens d'écoulement prévu, et un contrôleur d'activité de l'injecteur. Le contrôleur est agencé pour effectuer les opérations suivantes : 25 - recevoir des paramètres ai ; bl ; cil ; hl ; Ti ; Tfi, tels que o le système N1 N3/ N2 ji(t) = (aiy(t - Ti)) + b1 + (city(t - T'1)) Kt - i=1 /=1 j=1 y(S) = 0, u(s) = 0, pour S 0, V / = 1, , N3 b o pour ymax = N3 N2 Ymax (1 + EN1 a. T N3 1=1 1=1 CiIT'i) I1=1(b1h1) > y* 1 + N3 ENià_ aiTi -FENI=31E7=21Cii-C -ENI=31E7=21Ciihi - N3 . . -N3 -N2 10 1=1(b1111)+11=12cjihiy N E a- V y* E MYrrtaxE N3 . -N3 s-,N2 i=1 1=1(b1n1)+2;1=12cilY o les zéros des quasi-polynômes N1 N3 N2 MS) = aie-sri + cite -s'el - S. i=1 1=11=1 Ni (22(s) =aie-sri -s i.i soient à partie réelle négative , 15 o le système {7 Iv1 ii(t) = raizi(t - Ti)) + b1 + (citzi(t - T'1)) Z2(t - h1) + N3 N2 N3 N3 N2 1=11=1 1 i=1 1=1 \ j=1 1Vr1 1V3 7 1V2 \ z1(t) bi Z2 (t) (Z2 (t) 1)Z1(t - h1) = MS) = (I0 S E [-Tmax, 0] 1=1 1=1 1=1 soit localement asymptotiquement stable, - recevoir, et mémoriser pendant une durée minimum Tmax = max(ci, ..., hN3, T'Ar2), des mesures y ( t ) au - 5 - cours du temps t du moyen de mesure du décollement de l'écoulement le long du profil, - envoyer, et mémoriser pendant une durée minimum Tm= = max(ci, , , hN3, , TIN2), des instruction u(t)= 0 5 ou 1 au cours du temps t respectivement d'inactivité ou d'activité de l'injecteur, - calculer une grandeur a(t) au cours du temps t comme étant au moins la somme des termes y(t) ; ft_ri y(s)ds ; a2 ftt_T2 y(s)ds ; b1 ft_hiu(s)ds ; frC1110, (s)ds ; J+hi y. cil ft_hi y(s - + hi)u(s)ds , puis la comparer à un paramètre a* comme étant au moins la somme des termes y* ; y*alTi ; y*a2T2 ; ; y*a2h1 ; ; -y*cnhi , y* étant une valeur de consigne cible du moyen de mesure du décollement de 15 l'écoulement le long du profil, telle que y* E 10, N , et j=1C N3 N2 o si a(t) est plus petit que a*, envoyer un ordre d'activité de l'injecteur, o si a(t) est plus grand que a*, envoyer un ordre 20 d'inactivité de l'injecteur. La loi de commande selon l'invention est particulièrement fiable et efficace car elle est issue d'une définition du système physique décrit par l'invention 25 au moyen d'un modèle bilinéaire à retards inspiré des équations Navier-Stokes qui décrivent la physique des fluides. Ce modèle décrit dans Feingesicht et al. (2016) est comme suit : - 6 - N3 / N2 3./(t) = (cti,y(t - Ti)) + b1 + (cjiy(t - /.1 1=1 \ 1=1 N1 u(t - h1) avec y correspondant aux mesures de décollement l'écoulement au cours du temps t, avec u étant égal de à 0 pour une inactivité de l'injecteur, et égal à 1 pour une 5 activité de l'injecteur, avec ai ; bl ; cil des coefficients, 10 avec hl ; Ti ; Tfi des retards. En considérant que le modèle recoupait suffisamment la réalité d'un système physique d'écoulement autour du profil, la demanderesse a prouvé selon l'annexe A du présent mémoire que si les paramètres ai ; bl ; cil ; hl ; Ti ; Tfi, sont tels que o le système N1 N3/ N2 = (Ctiy(t (city(t - T'1)) 21.(t - i=1 1=1 j=1 15 y(s) = 0, u(s) = 0, pour s 0,V 1 = 1, ... , N3 o pour ymax = Ni N3 N2 j=ic Ymax (1 + EN a.T - + N3 N2 V y* E j=1CiIT' i)> y* 1 + EN' ,N3 ai T N3 N3 N2 Ei./(bito cjihiy* N rN3 VN2 ' VN3 VN2 +Et.iEj. E. a. z-,1=1L-,./.1 Li/ 1-111e1 -1=1 1 11 Y,N3 (b h -FE1=1E j=iCity* 1-1 1 5 ]°,Ymax[ o les zéros des quasi-polynômes Ni N3 N2 = a Q1(s) -srf - - s 1=1 1=1 Q2(S) = ae - s i=1 soient à partie réelle négative o le système N1 N3/ N2 {7 N1 \i=1 ai + C. 1=1 N3 N2 \ 11 z1(t) -F bi Z2 (t) 1=1 I 1=1 N3 1=1 1=1 N3 N2 ±1(t) = (aizi(t - Ti)) + b1 + (citzift - T' 1)) Z2(t - h1) i=1 1=1 \ j=1 Zi(S) = cf) , s E [- Tmax,O] 10 soit localement asymptotiquement stable, une loi de contrôle d'activation ou non de l'injecteur selon l'invention, qui est de type par modes glissants, permet de tendre vers une valeur cible y* du décollement mesuré y(t) en un temps fini. Cette valeur cible y* doit 15 être défini telle qu'elle fasse partie des valeurs atteignables par y(t) selon le modèle défini, donc être comprise entre 0 et une grandeur Ll=ib1 dépendant des paramètres choisis. Ymax = N N3 N2 j=i cil 20 Cette loi de contrôle est décrite dans l'invention dans son mode le plus minimaliste, c'est à dire se reposant au moins sur le modèle tel que N1=2 ; N2=1 et N3=1. Pour un - 8 - contrôle gagnant en précision, mais gagnant en complexité, on peut se reposer sur un modèle tel que N1>2 ; N2>1 ou N3>1. Dans ce cas, les sommes permettant de déterminer a(t) et a*, en vue d'effectuer le contrôle, gagneront des termes 5 supplémentaires en plus de ceux décrits a minima dans l'invention. En effet, dans le cas général, c'est à dire pour tout N, on exprime : N1 N3 t t = y(t) +(ai f y(s)ds)+(bi f u(s)ds) t-ri 1=1 t-h1 N3 N2 N3 N2 t t + Cii f Y(S - + 111)U(S)dS + Cil f y(s)ds 1=11=1 t-h1 1=1 j=1 t-T'j+111 10 et N1 N3 N2 V N3 V N2 h ( Lit j.1 N2 cjlY*) \ 0-* = 1 +aiTi (Ti] - h1) r N3 r N2 ( bl ' i=1 1=11=1 L'l=1L-V=1 N2 Ce contrôle par modes glissants qui découle spécifiquement du modèle, si les critères présentés ci- 15 avant sont respectés, asservit l'injection d'air de manière que la valeur y(t) de mesure du décollement de l'écoulement tende vers la valeur de consigne ciblée y* de décollement de l'écoulement n un temps fini. Les résultats expérimentaux ont bien démontré que ce contrôle était assez 20 robuste. L'annexe A décrit également les conditions à respecter pour les paramètres dans le cas plus simple dans lequel N1=2 ; N2=1 et N3=1. Si ces conditions sont respectées, 25 alors le dispositif permet de réguler la mesure y(t) de décollement de l'écoulement à une valeur cible y*. - 9 - Il est possible de simplifier le modèle en choisissant certains paramètres égaux, tels que les retards Ti et h par exemple. Dans ce cas, il en découle pour le contrôleur qu'il peut être configuré pour ne recevoir qu'un seul 5 paramètre au lieu de deux, et simplifier ainsi la loi de contrôle. Il est donc évident que ce cas particulier, qui nécessite au contrôleur de recevoir un paramètre en moins, est en fait couvert par l'invention décrite en tant que cas généralisé. 10 Afin de calculer la grandeur a(t) qui fait intervenir des intégrales, le contrôleur effectue des approximations en appliquant une méthode de calcul numérique d'intégrales de grandeurs discrètes, telles que la méthode des 15 rectangles, la méthode des trapèzes, la méthode Simpson, etc. Ainsi, grâce à l'utilisation d'un modèle prenant en compte les réalités physiques du système, et un contrôle 20 robuste et spécifiquement adapté, l'invention permet de contrôler finement le décollement de l'écoulement autour du profil par l'actionnement des injecteurs, comme le montrent les résultats expérimentaux. Le respect des conditions énoncées dans l'invention apporte l'assurance prouvée 25 mathématiquement que le dispositif pourra toujours être régulé en décollement d'écoulement en un temps fini, ce qui est un réel avantage en sécurité et en performances. De plus, la loi de contrôle est assez simple et économe en moyens de calcul pour être effectuée en temps réel. 30 Avantageusement, les actions qui consistent à - 10 - - recevoir, et mémoriser pendant une durée minimum tin', les mesures y(t), - envoyer, et mémoriser pendant une durée minimum tin', les instruction u(t), 5 - calculer la grandeur a(t), puis la comparer au paramètre a* et o si a(t) est plus petit que a*, envoyer un ordre d'activité de l'injecteur, o si a(t) est plus grand que a*, envoyer un ordre 10 d'inactivité de l'injecteur, sont répétées à une fréquence d'au moins 20Hz. Cette fréquence permet que les opérations effectuées par le contrôleur s'effectuent à un rythme suffisant pour que la 15 mesure y(t) de décollement ne s'éloigne jamais de la valeur cible y* selon un écart indésirable. Préférentiellement, le moyen de mesure du décollement de l'écoulement est un capteur par film chaud disposé à la 20 surface du profil dont la tension de sortie sert à calculer la mesure y. Ce film chaud peut être macrométrique ou micrométrique. Selon le principe de fonctionnement de ce 25 capteur, un film est chauffé par résistance à un courant électrique le parcourant, et l'on mesure son refroidissement causé par son environnement. Dans le cas ici présenté, lorsque l'écoulement est collé au profil, la couche limite d'écoulement reste en contact avec la surface 30 du profil, donc le film chaud. La couche limite d'écoulement présente une viscosité telle qu'elle possède un mouvement relatif faible voir nul au contact du profil. Il y a donc peu d'échanges de chaleur entre le fluide et le film chaud, qui reste donc chaud. Au contraire, lorsque l'écoulement se décolle, il se crée des vortex en mouvement par rapport au film chaud et qui échangent de la chaleur 5 avec ce dernier, le refroidissant ainsi. Un capteur film chaud est donc particulièrement adapté pour mesurer le taux du décollement du flux autour du profil. On trouve également des capteurs fil chaud qui ont un 10 principe de fonctionnement équivalent à celui des capteurs film chaud, bien qu'un fil soit utilisé à la place d'un film. Bien que le présent mémoire décrive principalement un cas d'usage de l'invention avec des capteurs film chaud, il faut considérer que l'alternative d'usage de capteurs fil 15 chaud est en tous points équivalente. On peut également imaginer utiliser un capteur de pression comme moyen de mesure du décollement de l'écoulement. On pourra disposer une pluralité de moyens de mesure 20 du décollement de l'écoulement sur le profil, de préférence alignés perpendiculairement au sens d'écoulement du fluide, et recevoir alors des mesures y(t) moyennées sur l'ensemble des moyens de mesure. Le décollement mesuré n'est ainsi pas soumis à certains effets locaux dans l'écoulement, et 25 représente mieux le décollement général de l'écoulement autour du profil. Afin d'améliorer la convergence de y(t) vers y*, et d'avoir donc un contrôle plus précis du décollement de 30 l'écoulement, les paramètres ai ; bl ; cil ; hl ; Ti ; Tfi sont choisis à partir d'une méthode d'identification entre - un modèle bilinéaire à retard - 12 - N1 N3/ N2 3i(t) = (ctiy(t - Ti)) + b1 +(city(t - T'1)) U(t i=1 1=1 \ 1=1 de comportement de la mesure y de décollement de l'écoulement dans le temps t en fonction d'états préalables de la mesure y de décollement et de 5 l'injection u, l'injection u prenant la valeur 0 pour l'état inactif, et la valeur 1 pour l'état actif, - et des enregistrements dans des conditions expérimentales d'écoulement autour du dispositif de la mesure y et de l'injection u. 10 Une méthode d'identification permet de rapprocher le modèle utilisé le plus possible des conditions particulières liées au profil et au type d'écoulement auquel il fait face. La méthode d'identification produit 15 des paramètres optimisés qui permettent une convergence plus rapide des mesures de décollement vers la valeur cible y*, au moyen de la loi de contrôle qui s'appuie sur ce modèle amélioré. 20 Lorsque on utilise la méthode d'identification ici décrite pour choisir les paramètres ai ; bl ; cil ; hl ; Ti ; on obtient d'excellents résultats en termes de précision par rapport aux enregistrements expérimentaux, ce qui prouve que le modèle est de très bonne qualité. En 25 effet, on a pu obtenir un critère d'erreur --11Jv - expérimental - cov(31 expérimentabY simulé) de corrélation p- eYexpérimentaleYsimulé 0,99, et un indicateur FIT=(1, Y expérimental-Ysimulél1L2 ) X100% > IlY expérimental-Y expérimentcallg 87,10%. Les paramètres issus d'une telle identification peuvent donc être considérés comme permettant au contrôleur Ysimutélliz 80% et plus élevée de 30% par rapport au modèle de l'art antérieur décrit dans Dandois et al. (2013). 5 Dans un mode de réalisation préféré de l'invention, la valeur de consigne cible y* du moyen de mesure du décollement de l'écoulement le long du profil varie dans le temps en fonction des performances recherchées. 10 En effet, il peut être avantageux dans certains cas de recherche un recollement maximal pour favoriser les effets d'aérodynamisme tel que la portance ou la suppression d'un sillage, quitte à sur-solliciter les injecteurs et perdre en performances. Dans d'autres cas, on pourra vouloir 15 minimiser l'énergie totale consommée et ne pas sur-solliciter les injecteurs. Dans un mode de fonctionnement avantageux de l'invention, les paramètres ai ; bl ; cil ; hl ; Ti ; Tfi 20 varient dans le temps en fonction de variations intrinsèques au profil et/ou dans l'écoulement du fluide autour du profil. Il est en effet connu, pour des véhicules, de posséder 25 des profils qui peuvent varier de formes au moyen de dispositifs mécaniques. On peut notamment citer les volets des avions ou les ailerons qui se déplacent en fonction des commandes de guidage. Il est en effet courant que les véhicules qui se déplacent dans un fluide rencontrent des 30 conditions d'écoulement variables, en fonction par exemple de la densité du fluide, de sa composition, de sa vitesse d'écoulement, de sa température, etc. Ces conditions varient donc au cours d'une même utilisation du véhicule. - 15 - Il est donc avantageux de disposer de paramètres de modèle particulièrement adaptés à ces différentes conditions rencontrées, et de les utiliser en temps voulu dans le contrôleur afin de conserver des performances de contrôle 5 homogènes quelles que soient ces conditions. Selon un mode de réalisation de l'invention, l'injecteur de fluide comprend une source de gaz comprimé reliée à un tube dont un orifice débouche à la surface du 10 profil. La source d'air comprimé peut être une chambre sous pression ou un moteur de mise sous pression. Le gaz peut être de l'air ambiant. De préférence, le tube débouche à la surface du profil selon un angle vers l'aval par rapport à l'écoulement. L'orifice possède typiquement un diamètre 15 compris entre 0,3mm et 2mm et peut posséder par exemple une forme ronde ou en fente. De préférence, l'injecteur débouche juste en amont d'une cassure de forme dans le profil. C'est à cet endroit 20 qu'il est le plus utile pour recoller l'écoulement au niveau de la cassure en question. Il peut y avoir une pluralité d'injecteur le long du profil, de préférence alignés perpendiculairement par rapport au sens d'écoulement afin d'obtenir une homogénéité de recollement 25 de l'écoulement sur l'ensemble du profil. L'invention concerne également un aéronef comprenant un dispositif tel que décrit précédemment. 30 L'invention concerne aussi un procédé de contrôle d'un injecteur de fluide en périphérie d'un profil, au moins un moyen de mesure du décollement d'un écoulement le long du profil étant disposé en aval de l'injecteur par rapport au - 16 - sens d'écoulement prévu. Le procédé de contrôle comprend les étapes suivantes : - recevoir des paramètres ai ; bl ; cil ; hl ; Ti ; Tfj, tels que 5 o le système N3/ N2 3./(t) = (12,,y(t - T3) b1+(cj1Y(t - j)) 1-1 1-1 j=1 y(S) = 0, U(S) = 0, pour s 0, / = , N3 N3 b 15 o pour ymax = N3 N2 j=1 Cil u(t- h1) Ymax(1 vNi 1h1) > Y* z, a. y+ VN3 vN2 i L-1/=1L]=1CjiT'j) EN3 (b 20 EN' a.- c- c. T +EN3 N2 ' - EN3 c- h 1-1 1-1 Ji I - 1=1(b1111)+2: 24=11j=iC jihlY +ENi_31E iN=21c V Y* e MYmax[f Ei.1 ai ENi_31(bito o les zéros des quasi-polynômes N1 N3 N2 Q1 (S) = aie-sri + ciie-"' i - s i=1 1=11=1 Ni (22(s) = aie-sri -s i.i soient à partie réelle négative, o le système - 17 - {7 Ni N3 N2 N3 N3 N2 \ C.11 (t) bi Z2 (t) (Z2 (t) 1)Z1(t - h1) = 1=1 j=1 1=1 1=1 j=1-a, + MS) = (I0 S E [-Tmax, 0] Ni 1V3 7 IV2 ii(t) = (aizi(T - Ti)) + b1 + (cjizi(t - T'1)) i=1 1=1 \ 1=1 z2(t - h1) soit localement asymptotiquement stable, - recevoir, et mémoriser pendant une durée minimum Tmax = maxeci, -cAti, ... hN3, des mesures y(t) au 5 cours du temps t du moyen de mesure du décollement de l'écoulement le long du profil, - envoyer, et mémoriser pendant une durée minimum Tmax= max(Ti, hN3, ...,T'N2), des instruction u(t)= 0 ou 1 au cours du temps t respectivement d'inactivité 10 ou d'activité de l'injecteur, - calculer une grandeur a(t) au cours du temps t comme étant au moins la somme des termes y(t) ; - ft_ri y(s)ds ; a2t y(s)ds ; b1 ft-h1u(s)ds ; - it_r, -En y(s)ds ; cn ft_hiY(s + hi)u(s)ds , puis la 15 comparer à un paramètre a* comme étant au moins la somme des termes y* ; ; y*a2T2 ; ; y*a2h1 ; ; -y*cntil , y* étant une valeur de consigne cible du moyen de mesure du décollement de l'écoulement le long du profil, telle que y* E N3 20 10, N N N2 cil , et o si a(t) est plus petit que a*, envoyer un ordre d'activité de l'injecteur, o si a(t) est plus grand que a*, envoyer un ordre d'inactivité de l'injecteur. 25 - 18 - D'autres caractéristiques et avantages de l'invention apparaîtront mieux à la lecture de la description qui suit, tirée d'exemples donnés à titre illustratif et non limitatif, tirés des dessins sur lesquels : 5 la figure 1 représente un dispositif de pénétration dans un fluide selon l'invention vu de côté par rapport au sens d'écoulement du fluide, en mode passif, la figure 2 représente un dispositif de pénétration 10 dans un fluide selon l'invention vu de côté par rapport au sens d'écoulement du fluide, en mode actif, la figure 3 représente l'algorithme de fonctionnement du contrôleur selon l'invention, 15 la figure 4 représente des mesures expérimentales de décollement de l'écoulement grâce à l'invention. L'annexe B rassemble les références de publications citées dans le présent mémoire. 20 Les dessins et la description ci-après contiennent, pour l'essentiel, des éléments de caractère certain. Ils pourront donc non seulement servir à mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant. 25 La présente description est de nature à faire intervenir des éléments susceptibles de protection par le droit d'auteur et/ou le copyright. Le titulaire des droits n'a pas d'objection à la reproduction à l'identique par 30 quiconque du présent document de brevet ou de sa description, telle qu'elle apparaît dans les dossiers - 19 - officiels. Pour le reste, il réserve intégralement ses droits. En outre, la description détaillée est augmentée de 5 l'annexe A, qui donne la formulation de preuves mathématiques mises en oeuvre dans le cadre de l'invention. Cette Annexe est mise à part dans un but de clarification, et pour faciliter les renvois. Elle est partie intégrante de la description, et pourra donc non seulement servir à 10 mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant. Dans la figure 1, on représente un profil 10 de pénétration dans un fluide. Le profil 10 est prévu pour 15 pénétrer dans un fluide selon un sens général prédéfini, ici vers la gauche. Le profil 10 possède une zone centrale élargie depuis son extrémité gauche, puis, plus à droite, une rupture 20 de profil 10 depuis laquelle le profil 10 présente un amincissement abrupt vers la droite, jusqu'à 20 son extrémité droite aussi appelée bord de fuite. Ce profil 10 correspond ici à celui d'une aile d'avion qui permet un effet aérodynamique de portance, quand il pénètre dans l'air vers la gauche. 25 On représente par une flèche 18 l'écoulement d'un fluide, ici de l'air, autour du profil. L'origine de l'écoulement correspond à l'amont et la destination de l'écoulement correspond à l'aval. 30 On dispose sur le profil, juste à gauche, ou en amont, de la rupture 20 de profil, des orifices chacun reliés à un tube connectant l'orifice à une source d'air comprimé. L'ensemble source d'air comprimé, tube et orifice forme un - 20 - injecteur 14. Les orifices des injecteurs 14 possèdent un diamètre d'environ 1 mm, et pouvant aller de 0,2 mm à 1 cm, sont circulaires, bien qu'ils puissent aussi avoir la forme de fentes ou autres, et peuvent être orientés selon des 5 angles variables par rapport à la normale du profil. Alternativement, et selon des modes également fonctionnels, on peut disposer les orifices des injecteurs 14 au niveau de la rupture 20 de profil, ou bien juste en aval de la rupture 20 de profil. 10 Dans l'exemple ici décrit, les injecteurs 14 d'air comprennent une chambre d'air comprimé. Il est également possible de manière fonctionnelle d'utiliser des moteurs de compression d'air. Dans des cas de profils de pénétration 15 dans un liquide, on utilisera de manière équivalente des injecteurs de liquide. Il est également possible d'utiliser des injecteurs de gaz pour un profil de pénétration dans un liquide, ou des 20 injecteurs d'eau pour un profil de pénétration dans un gaz. Les injecteurs 14 comprennent des moyens d'activation ou non de l'injection, représentés sur la figure 1 par des clapets d'ouverture ou non de la chambre sous pression. 25 Dans le cas d'utilisation de moteurs de compression d'air, il pourra s'agir de commutateur de mise en marche des moteurs. Les injecteurs 14 injectent de l'air à un débit 30 massique à partir de 10 g/s, autour de 25g/s de préférence, ou à une vitesse à partir de 15 m/s. - 21 - On dispose des moyens de mesure de décollement de l'écoulement du fluide le long du profil. Ces moyens de mesure correspondent ici à des capteurs 16 film chaud disposés à la surface du profil 10 à droite de la rupture 5 20 de profil, au niveau de l'amincissement du profil. Il est également possible d'utiliser des capteurs fil chaud, des capteurs de pression, etc., comme moyens de mesure de décollement. 10 Un capteur film chaud permet de mesurer le refroidissement d'un film métallique chauffé en continu par asservissement de manière connue. Le film chaud peut posséder plusieurs dimensions, de l'ordre du nanomètre au millimètre. 15 Selon les lois d'écoulement d'un fluide, lorsque ce dernier adopte un comportement laminaire le long d'un profil selon lequel l'écoulement est recollé au profil, ce qui est recherché pour les effets aérodynamiques tels que 20 la diminution de traînée, l'écoulement présente une couche limite en contact du profil dans laquelle le fluide présente une viscosité telle que le fluide possède une vitesse relative au profil quasi-nul. De ce fait, il y a peu d'échange d'énergie avec le profil et le film chaud du 25 capteur reste chaud. Au contraire, lorsque l'écoulement présente des tourbillons ou vortex en périphérie du profil, caractéristiques d'un décollement de l'écoulement autour du profil et néfastes pour les effets aérodynamiques généralement recherchés, le fluide possède une vitesse 30 relative au profil élevée, et il y a beaucoup d'échange de chaleur avec le film chaud du capteur, refroidissant ainsi ce dernier. - 22 - Le film chaud est donc bien adapté à la mesure de décollement de l'écoulement le long du profil. Toutefois, il est aussi connu et fonctionnel d'utiliser d'autres types de capteurs, comme des capteurs de pression en surface du 5 profil, des capteurs optiques, des capteurs de débit, vitesse, du fluide, ou autres, pour mesurer ce décollement. Un contrôleur 12 est agencé dans le dispositif. Il est capable de recevoir l'information de mesure des capteurs 10 16, et de contrôler la mise en activité ou non des injecteurs 14. Le contrôleur 12 possède une mémoire 22 de stockage de données. Dans le cadre de l'invention, la mémoire 22 peut 15 être tout type de stockage de données propre à recevoir des données numériques : disque dur, disque dur à mémoire flash (SSD en anglais), mémoire flash sous toute forme, mémoire vive, disque magnétique, stockage distribué localement ou dans le cloud, etc. Les données calculées par le dispositif 20 peuvent être stockées sur tout type de mémoire similaire à la mémoire, ou sur celle-ci. Ces données peuvent être effacées après que le dispositif a effectué ses tâches ou conservées. 25 Le contrôleur 12 possède également un calculateur 24 accédant directement ou indirectement à la mémoire 22. Il peut être réalisé sous la forme d'un code informatique approprié exécuté sur un ou plusieurs processeurs. Par processeurs, il doit être compris tout processeur adapté 30 aux calculs décrits plus bas. Un tel processeur peut être réalisé de toute manière connue, sous la forme d'un microprocesseur pour ordinateur personnel, d'une puce dédiée de type FPGA ou SoC (« system on chip » en anglais), - 23 - d'une ressource de calcul sur une grille, d'un microcontrôleur, ou de toute autre forme propre à fournir la puissance de calcul nécessaire à la réalisation décrite plus bas. Un ou plusieurs de ces éléments peuvent également 5 être réalisés sous la forme de circuits électroniques spécialisés tel un ASIC. Une combinaison de processeur et de circuits électroniques peut également être envisagée. Dans le cadre de l'exemple ici décrit, on utilise 10 comme contrôleur 12 un système connu sous le nom de Arduino® remarquable par sa simplicité et sa faible utilisation de ressources. Le contrôleur 12 selon l'invention nécessite en effet un processeur fonctionnant à une fréquence de seulement 4MHz, et une mémoire 22 d'une 15 taille totale de seulement quelques kilooctets, ce qui implique une très faible consommation d'énergie et un très faible coût d'implémentation. Comme montré en figure 3, le contrôleur 12 est 20 programmé pour stocker à une étape 110 dans sa mémoire 22 des paramètres al ; a2 ; bl ; cl ; hl ; Tl ; T2 ; Tfl, ainsi que pendant une durée Tma x égale à max(h,TifT2,T') les mesures y(t) de décollement de l'écoulement et les commandes u(t) égales à 1 pour une activation des 25 injecteurs 14 et égales à 0 pour une désactivation des injecteurs 14. Pour -t.0, il est considéré que y(t)=0 et u(t)=0. Une valeur de consigne cible y*, qui correspond à la valeur de mesure y(t) de décollement que l'on souhaite atteindre, est également stockée en mémoire 22. 30 Lors d'une étape 112 qui suit, le contrôleur 12 est configuré pour calculer deux grandeurs a(t) et a* et les - 24 - comparer entre elles. La grandeur a* est calculée à partir des paramètres stockés en mémoire 22 et est telle que o-* = y* + y*ai-ci + y*a2T2 - y*ailii - y*a2hi + y*citi - y*cilii 5 La grandeur a(t) est calculée à partir des paramètres, des mesures y(t) et des commandes u(t) stockés en mémoire 22 et est telle que t t r a(t) = y(t) + a2 ft_r2 y(s)ds + c1 f ,t-rl-Fik y(s)ds + al ft_ri y(s)ds + b1 ftt_h1 u(s)ds + cl ft y(s - -C' + 11)1,t(S)dS 10 Si a(t) a*, alors à une étape 114 on attribue 1 à u(t) et on positionne les injecteurs 14 15 dans un état activé à une étape 118. Puis, on incrémente t selon la fréquence de fonctionnement du contrôleur 12, de préférence supérieure à 25Hz. Les valeurs de y(t) et u(t) pendant la durée 20 d'incrémentation sont stockées en mémoire 22 en répétant l'étape 110, et l'algorithme de contrôle est répété. Nous allons maintenant expliquer pourquoi cette configuration du contrôleur 12 permet de réguler la mesure 25 y(t) de décollement de l'écoulement à la valeur de consigne ciblée y* avant une durée bornée. Selon l' invention, les paramètres al ; a2 ; b1 ; cl ; hl ; 'cl ; 'C2 ; 'Cf 1 sont choisis tels que 30 o le système - 25 - N1 N3/ N2 3i(t) = (ctiy(t - Ti)) b1+ (cily(t - T'1)) - h1) i=1 1=1 \ 1=1 y(S) = 0, U(S) = 0, pour s 0, / = 1, ..., N3 b 10 o pour ymax = N3 N2 j=1 Cil Ymax(1+ EN' a.T yN3 rN2 N3 N1 > Y* 1 + Ei=laiTi N3 . s-,N3 vN2 ' v N3 VN2 Y.N3 chl Et=i(bih1)+2,1=12L EN.1 a. /-1 Lit=i -/-1 E (b010+2, N3 . -N3 vN2 1=1 1=1 1 '11Y V Y* e YmaxL o les zéros des quasi-polynômes N1 N3 N2 15 Q1(s) = ai e + c-ie-sr' - - s i=1 1=1 j=1 (22 (s) = aie - S soient à partie réelle négative, o le système 0 - 26 - { N1 N3/ N2 z2(t - h1) ±1(t) = (t - Ti)) + + (cjizift - i=1 1=1 \ 1=1 7 N1 N3 N2 \ N3 N3 N2 a, + Cil z1 (t) + bi Z2 (t) + Cii (Z2 (t) - 1)Zi (t t'j = 1=1 1=1 j=1 z1(s) = (P, s E [-Tmax, 0] soit localement asymptotiquement stable. Alors, comme prouvé dans l'annexe A, la loi de 5 contrôle décrite ci-avant permet bien de réguler la mesure y(t) de décollement à la valeur cible y* en un temps fini. Afin de mieux choisir les paramètres al ; a2 ; bl ; cl ; hl ; Ti ; T2 ; Tfl, on a en plus procédé à une identification 10 du modèle utilisé 3i(t) = aiy(t - + a2y(t - -c2) + (b1 + clly(t - -C1))u(t - h1) où interviennent ces paramètres, par rapport à des situations expérimentales du profil 10 et de l'écoulement dans lesquelles sont enregistrées les mesures y(t) de 15 décollement de l'écoulement et les états u(t)=0 ou 1 inactifs ou actifs des injecteurs 14. Plus précisément, afin d'effectuer cette identification, on propose des retards hl ; Ti ; Tfi. En 20 intégrant alors les valeurs y(t) et u(t) expérimentales au modèle, on détermine les coefficients ai ; bl ; cil par la méthode des moindres carrés. On calcule alors soit le critère d'erreur P1/ -V expérimental simulélliz soit le critère de cov(YexpérimentalYsimuié) corrélation p- , soit l'indicateur FIT corrélation eY expérimentale)/ simulé I1Y expérimental-Y simulél1L2 25 x100%. Y expérimental-Y expérimentalll - 27 - Un algorithme d'optimisation globale, NOMAD ou bien l'algorithme Genetic, est utilisé pour répéter en boucle ces étapes. Il est configuré pour proposer à chaque 5 itération de nouveaux retards hl ; Ti ; Tfi qui, après identification des coefficients ai ; bl ; cil par la méthode des moindres carrés, découleront soit sur une minimisation du critère d'erreur, soit sur une maximisation du critère de corrélation, soit sur une maximisation de l'indicateur 10 FIT. L'algorithme d'optimisation globale est configuré pour s'arrêter quand les critères atteignent une valeur cible ou bien ne varient plus suffisamment à chaque itération. L'algorithme Genetic est décrit dans Davis, L. (1991), 15 Handbook of genetic algorithms. L'algorithme NOMAD est connu et décrit par exemple dans Audet et al. (2009), Le Digabel (2011). Dans la figure 1, le profil 10 se déplace à une 20 vitesse d'environ 40 m/s par rapport à l'air. Le contrôleur 12 est désactivé. Les injecteurs 14 ne fonctionnent donc jamais, et il apparaît dans ces conditions des vortex 20 le long de la zone d'amincissement du profil 10 en aval de la rupture 20. Ces vortex 20 apparaissent à partir d'une 25 vitesse d'écoulement d'environ 10 m/s. Ils créent une dépression en aval du profil 10 qui induit une force de traction du profil 10 vers l'aval, ou traînée. Dans la figure 2, le contrôleur 12 est activé de 30 manière à réguler la mesure y(t) de décollement de l'écoulement à une valeur cible y*=0,25 , en agissant sur les injecteurs 14. Pour les capteurs 16 film chaud, cette valeur correspond à un bon recollement de l'écoulement. Les - 28 - injecteurs 14 sont donc activés périodiquement suivant la loi de contrôle décrite précédemment. Grâce à ce contrôle actif des injecteurs 14, l'écoulement 18 reste collé au profil 10 en aval de la rupture 20 de profil, et on ne voit 5 pas apparaître de vortex 20. La figure 4 fournit les mesures relevées lors des cas expérimentaux représentés en figures 1 et 2. La courbe 210 indique le relevé des mesures y(t) de décollement selon les 10 capteurs 16 film chaud en ordonnée dans le temps t en secondes indiqué en abscisse. La courbe 212 indique le relevé des injections u(t) en ordonnée dans le temps t en secondes indiqué en abscisse. Les deux courbes sont synchronisées dans le temps. De t=7 secondes à t=40 15 secondes, on active le contrôleur 12 selon l'invention. Le reste du temps, on désactive le contrôleur 12. On voit bien qu'à l'activation du contrôleur 12, la mesure y(t) de décollement passe d'une valeur moyenne de 0, significative d'un mauvais recollement de l'écoulement, à une valeur 20 moyenne régulée de 0,25, significative d'un bon recollement de l'écoulement, et qui correspond à la valeur de consigne cible y* recherchée. De plus, on observe sur la zone agrandie 214 de la courbe u(t) que les injecteurs 14 sont activés à intervalles variables, ce qui évite de 25 surconsommer en moyens d'injection. L'invention est donc bien fonctionnelle et permet bien de réguler le décollement de l'écoulement à une valeur voulue. - 29 - Annexe A I Proposition On note : N1 N3 N3 N2 i=1 1=1 1=1 j=1 N1 N3 N3 N2 N3 N2 ; cr, = ciir'i ; cil= C-iiii I i=1 1=1 1=1 j=1 1=1 j=1 Ymax - - a + c Tmax = max(ri, TATI, hi, ..., hN3, r'i, ...,r'N2) Si : les coefficients du modèle sont tels que : à+ E O,V/= 1, , N3 (3) 761 + * + Ey* a) V Y* E MYmax[- - Condition I : le système 5 10 15 ( (4) N3/ 25 5/(t) = (aiy(t - Ti)) + b1 + (citY(t - i=1 1=1 \ j=1 y(s) = 0, u(s) = 0, pour s c11(z(t) - 1)z1ft 1=1 j=1 z1(s) = cp, s E [--c,,,,,, 0] est localement asymptotiquement stable ; alors la loi de contrôle u(t) = {1 si 0-(t) a* (6) 20 avec - 31 - N3 0-(t) = y(t) +(a, f y(s) ds) +(b1 f u(s)ds) 1=1 t-ri 1=1 t-h1 et N3 N2 Cil 1=11=1 N3 N2 Cil 1=11=1 y(s - -c' + hi)u(s)ds fr y(s)ds Jr-T' 1F111 ft-h1 (7) a b +y* ) et y* E Ymax[ - garantit y(t) -> y* lorsque t -> 00 . De plus, 0 0 et tout U E ,GOE) (Agarwal 2012, Lemme 2.1) . Remarque : La Condition 1 est vérifiée si (voir 15 Agarwal 2012, Theorème 2.6) il existe une constante À > 0 telle que : N1 N3 N2 (CtiFe) 0 : 6-(t) = (à + E)y(t) + 5 10 Or, puisque à +c 0 . - Si o-(t) > o-* , alors u(t) = 0 et il faut d(t) O. Supposons l' inverse : o-(t) O. Cela 5 signifie que u(t) = 1 pour tout t > O. Alors : N3 N3 N2 3> (T) =(ctiY(T - Ti)) +b1 +cji,Y(T - i) i=1 1=1 1=1 1=1 Si les zéros de Q1(s) sont à partie réelle négative 10 alors ce système est globalement asymptotiquement stable et converge vers son unique point d' équilibre Ymax tel que IIY(t) Ymaxiic 0 - t->+. On peut alors écrire Y(t) = Ymax + E(t) avec E(t) O. t->+. 15 On obtient : N1 ft Jt-Ti (Ymax + E(s)) ds) a(t) = Ymax + E(t) + ( a 1-1 N3 N2 C11 ..I (YmaxaX E(S))dS t-T1' 1=111 20 G(t) = (Ymax c-C-rYmax + N3 N2 E(t) (ai f E(s) cis) + cji f E(s)cis \ i.1 t-ri / 1=1 j=1 t-rf Puisque E(t) -> 0 alors (E(t) EA:à (a, ftt_ri E(t) ds) + t->+. ENià EiN=21 c11 ft_r, E(t)ds) 0 et a(t)-> Gmax - Ymax(1 + + cr,) t->+co t_)+00 25 bh . Grâce à l'equation (5), on obtient que a -max > ci* pour - 34 - tout Y* e Yn[Puisque o (t) t->+co ()Max , alors il existe un temps t* > 0 tel que o-(t) > o-* pour tout Y* e MYmax[ En utilisant la méthode du contrôle équivalent (V. 5 Utkin, Sliding modes in gptimization and control problems. Springer Verlag, New York, 1992. ) , on obtient : ENt=31E7-21 citY(t + hl) - (à + î')y(t) ueq(t) = +IN1=31I I=21Ci1Y(t j + hl) Le système obtenu en remplaçant u(t) par ueg(t) dans 10 (1) doit être asymptotiquement stable pour garantir y (t) -> y* lorsque t -> 00. , ce qui est donné par la Condition 3 de la Proposition. En remplaçant u par Ueq et y par y dans a, on 15 obtient a* : N3 N2 0-* (t) = y* (t) + (ai t y* (s) cls) f y* (s)ds 1=1 -Ti 1=1 j=1 t-TI ft-h1 N3 N2 Ueq(S) dS) C.11 f Y* hjueq(s)ds 1=1 j=1 t-h1 1=1 N3 N2 Ci* (t) = y* (t) (ai ft-Tiy* (s) ds) i+hiy* (s)ds i=1 1=1 j=1 f (-N9 + * T'i hi)) Ueq (s) ds t-h1 1=1 j=1 20 N3 N2 t - 35 - N3 N2 rt (t) = y* (t) +(ai f y* (s) ds) cjiJ y* (s)ds i=1 t-ri J-Erti 1=1 j=1 bi / N3 N2 Cpciy(S N3 N2 t CitY*(S + hi) 1 Y 2 ft-hi vN3 (LM city(s Tin + hm) \q=1 p=1 1=1 j=1 z-Im=1z-In=1 N2 - + hq) - (à + E)y(t) N1 N3 N2 t t Ci* (t) = y* (t) + (ai f y*() ds) +c'11 t-ri f t-rf i+hiy*(s)ds 1.1 j=1 bi + i3 i2 f t - + Cj/Y*(S - T'i + hi) N2 1=1 1=1 j=1 L vim.1Z-In=1 (e2 + CilY(S - T'il + N3 V b ds 5 7 N3 N2 hm)) \q=1 p=1 Cpcy hq) ds bi N3 N2 t + cily* - + hi) v2 - ((à 1=1 j=1 ft-hi rN3 VN2 (-brn c -iy(s - T'n + hm)) z-Im=1z-In=1 N2 10 + E)y(t)) ds Si v est une constante : 15 a b +y* ) Si N3 = 1 : 5 III Cas particuliers Si N12,N21,N31 : 10 ji(t) = aiy(t - T1) + a2y(t - T2) + (b1 + clly(t - - hl) N2 Cf* = y* 1 + (ri - hi) a2 (x2 - hi) + - hi) a2 0 C11 15 Ymax bi a2 20 Tmax = max(ci, T2, h1, ti1) (1 + cti(ci - hi) + a2(T2 - hi) + - hi)) > 0 a(t) = y(t) + y(s) ds + a2 f y(s) ds + b1 f u(s) ds ft-T1 t-T2 t-111 y(s)ds ci* = Y* (1 + ai (Ti - h1) + a2 (x2 - h1) + (el - hl)) + cll f y(s - t'l hi)U(S)dS + Cl t-h1 ft-rf - 37 - - Condition I : le système (1) est positif : pour tout signal d' entrée non trivial u E L', on a 0 G y(t) . - Condition 2 : les zéros de Q1(s) et Q2(s) sont à partie réelle négative : Q1(s) = aie-srl + a2e-sr2 + clie-sr'l - S. Q2(S) = aie-sri + a2e-sr2 - S - Condition 3 : pour tout signal d' entrée u E LOE), la solution du système (1) vérifie 15 (a1 + a2 + cii)y(t) G clly(t --Cl+ hl) - Condition 4 : le système til(t) = aizi(t - T1) + azzi(t - T2) + (b1 + clizi(t - T'1))z2(t - hl) (a1 + a2 + cii)zi(t) + b1z2(t) + c11(z2(t) - 1)z1(t - -C1 + h1) = 0 z1(s) = cf), s E [-Tmax, 0] 20 est localement asymptotiquement stable. 5 10 - 38 - Annexe B Liste des références : 5 Brunton SL, Noack BR. Closed-loop turbulence control: progress and challenges. Applied Mechanics Reviews 2015 Volino RJ. Separation Control on Low-Pressure Turbine Airfoils Using Synthetic Vortex Generator Jets, In ASME 10 Turbo Expo 2003, collocated with the 2003 International Joint Power Generation Conference., Atlanta, Georgia, USA, 2003:845-859 Tian Y, Cattafesta L, Mittal R. Adaptive control of 15 separated flow, 2006. American Institute of Aeronautics and Astronautics, Audet C, Jr J. ED, Digabel SL. Globalization strategies for Mesh Adaptive Direct Search. Computational Optimization and 20 Applications 2009; 46(2):193-215 Le Digabel S. Algorithm 909: NOMAD: Nonlinear optimization with the MADS algorithm. ACM Transactions on Mathematical Software (TOMS) 2011; 37(4):1-15 25 Feingesicht M, Raibaudo C, Polyakov A, Kerherve F, Richard J-P. A bilinear input-output model with statedependent delay for separated flow control, 2016 30 Dandois J, Garnier E, Pamart P-Y. NARX modelling of unsteady separation control. Experiments in Fluids 2013; 54(2):1-17 - 39 - Agarwal RP, Berezansky L, Braverman E, Domoshnitsky A. Nonoscillation Theory of Functional Differential Equations with Applications. Springer New York : New York, NY, REVENDICATIONS1. Dispositif de pénétration dans un fluide, comprenant au moins un profil (10) , un injecteur (14) de fluide en périphérie du profil, au moins un moyen de mesure (16) de décollement d' un écoulement (18) le long du profil en aval de l' injecteur par rapport à un sens d' écoulement prévu, et un contrôleur (12) d' activité de l' injecteur, 10 caractérisé en ce que le contrôleur est agencé pour effectuer les opérations suivantes : - recevoir des paramètres ai ; b1 ; cii ; ; ; , tels que o le système N1 N3/ N2 15 3i(t) = (ctiy(t - Ti)) + b1 + (cily(t - T '1)) - h1) i=1 1=1 \ 1=1 y(S) = 0, U(S) = 0, pour S 0,V / = 1, , N3 b o pour ymax = N3 N2 j=ic - 41 - Ymax (1 + EN1 a.TN3 CIT') ,N3 jj L1.1(b1h1) > y* 1 + ENià_ aiTi -FENI=31E7=21CiiT'i -ENI=31E7=21Ciihi- 'N2 V Y* E MYmax[r -ini)-N3 2,1=12ciiii-IY E; ai N3 , 1=1(b1n1)+11=11j=1C1131 o les zéros des quasi-polynômes N1 N3 N2 5 Q1(s) = ae +Cile-STI S 1=1 1=11=1 N1 (22(s) = aie-sri - s soient à partie réelle négative, o le système N1 N3 N2 il(t) = (aizi(t - -ci)) + b1 + (ciizi(t - '1)) Z2(t h1) i=1 1=1 \ j=1 N1 N3 N2 \ N3 N3 N2 ai +c11 z1(t) Z2(t) -FCii(Z2(t) - 1)Z1(t - + h1) = \i=1 1=1 j=1 1 1=1 1=1 j=1 z1(s) = (i0 S E [-Tmax, 0] soit localement asymptotiquement stable, - recevoir, et mémoriser pendant une durée minimum Tmax = max(cl, TATI, ..., hN3, T'Ar2), des mesures y (t) au cours du temps t du moyen de mesure du décollement de l' écoulement le long du profil, - envoyer, et mémoriser pendant une durée minimum Tmax = max(cl, ..., hN3, T'Ar2), des instruction u (t) = 0 ou 1 au cours du temps t respectivement d'inactivité ou d'activité de l'injecteur, - calculer une grandeur o(t) au cours du temps t comme étant au moins la somme des termes y (t) ; ft_ri y(s)ds ; a2t y(s)ds ; b1 ft-h1u(s)ds - 42 - rtt cil J-rfi-Ehi Y(S)ds ; cll fY(s - + hi)u(s)ds , puis la comparer à un paramètre a* comme étant au moins la somme des termes y* ; Y*a]:«ci ; y*a2T2 ; ; y*a2h1 ; ; -y*cntil , y* étant une valeur de 5 consigne cible du moyen de mesure du décollement de l'écoulement le long du profil, telle que y* E 10, N , et j=1C il N3 N2 o si a(t) est plus petit que a*, envoyer un ordre d'activité de l'injecteur, 10 o si a(t) est plus grand que a*, envoyer un ordre d'inactivité de l'injecteur. 2. Dispositif selon la revendication 1, caractérisé en ce que le contrôleur (12) est agencé pour 15 - recevoir, et mémoriser pendant une durée minimum Tmaxl les mesures y(t), - envoyer, et mémoriser pendant une durée minimum Tmaxl les instruction u(t), - calculer la grandeur a(t), puis la comparer au 20 paramètre a* et o si a(t) est plus petit que a*, envoyer un ordre d'activité de l'injecteur (14), o si a(t) est plus grand que a*, envoyer un ordre d'inactivité de l'injecteur, à une fréquence d'au moins 20Hz. 3. Dispositif selon l'une des revendications 1 ou 2, caractérisé en ce que le moyen de mesure (16) du décollement de l'écoulement (18) est un capteur par film - 43 - chaud disposé à la surface du profil (10) dont la tension de sortie sert à calculer la mesure y. 4. Dispositif selon l'une des revendication 1 à 3, 5 caractérisé en ce que les paramètres ai ; bl ; cil ; hl ; Ti ; -CF sont choisis à partir d'une méthode d'identification entre - un modèle bilinéaire à retard N1 N3/ N2 3i(t) = (ctiy(t - Ti)) + b1 + (city(t - T'1)) U(t i=1 1=1 \ 1=1 10 de comportement de la mesure y de décollement de l'écoulement dans le temps t en fonction d'états préalables de la mesure y de décollement et de l'injection u, l'injection u prenant la valeur 0 pour l'état inactif, et la valeur 1 pour l'état actif, 15 - et des enregistrements dans des conditions expérimentales d'écoulement autour du dispositif de la mesure y et de l'injection u. 5. Dispositif selon la revendication 4, caractérisé en 20 ce que la méthode d'identification est réalisée au moyen de la méthode des moindres carrés et d'un algorithme d'optimisation globale, tel que l'algorithme NOMAD ou l'algorithme Genetic. 25 6. Dispositif selon la revendication 5, caractérisé en ce que l'algorithme d'optimisation globale est configuré pour optimiser les paramètres hl ; Ti ; Tfi dits retards, en proposant par itération des jeux de retards soit minimisant le critère d'erreur --11Jv - expérimental simulélliz soit maximisant - 44 - cov(YexpérimentalYsimulé) le critère de corrélation P- soit -Y expérimentale)/ simulé maximisant l'indicateur FIT=(1 IlY expérimental-Y simulél X100%, I 1L2 , Il Yexpérimental-YexpérimentallIg) un jeu de paramètres ai ; bl ; cil , dits coefficients, étant calculé à chaque itération au moyen de la méthode des 5 moindres carrés afin de déterminer ysimuié. 7. Dispositif selon l'une des revendications précédentes, caractérisé en ce que le contrôleur (12) est agencé pour faire varier la valeur de consigne cible y* du 10 moyen de mesure du décollement de l'écoulement (18) le long du profil (10) dans le temps en fonction des performances recherchées. 8. Dispositif selon l'une des revendications 15 précédentes, caractérisé en ce que le contrôleur (12) est agencé pour faire varier les paramètres ai ; bl ; cil ; hi ; Ti ; Tfi dans le temps en fonction de variations intrinsèques au profil (10) et/ou dans l'écoulement du fluide autour du profil. 20 9. Dispositif selon l'une des revendications précédentes, caractérisé en ce que l'injecteur (14) de fluide comprend une source de gaz comprimé reliée à un tube dont un orifice débouche à la surface du profil. 25 10. Aéronef comprenant le dispositif selon l'une des revendications précédentes.