FR 3041121 A1 20170317 FR 1655600 A 20160616 1 PROCEDE DE CONTROLE D'UN AERONEF SUIVEUR PAR RAPPORT A DES TOURBILLONS GENERES PAR UN AERONEF MENEUR DOMAINE TECHNIQUE La présente invention concerne un procédé de contrôle d'un aéronef suiveur permettant d'évaluer les positions des centres de tourbillons générés par un aéronef meneur pour permettre de commander le positionnement de l'aéronef suiveur par rapport à ces centres de tourbillons, ainsi qu'un aéronef pour lequel un tel procédé de contrôle est mis en oeuvre. ETAT DE LA TECHNIQUE ANTERIEURE Un aéronef en vol génère dans son sillage deux tourbillons, avec un tourbillon partant de chacune de ses ailes du fait de la différence de pression entre l'intrados et l'extrados de l'aile, et de la déflexion de l'écoulement d'air vers le bas qui en résulte. Ces tourbillons sont des tourbillons contrarotatifs et sont caractérisés par un champ de vent globalement ascendant à l'extérieur des tourbillons et globalement descendant entre les tourbillons. A partir des ailes, les tourbillons tendent d'abord à se rapprocher l'un de l'autre, puis à maintenir une distance plus ou moins constante entre eux tout en perdant de l'altitude par rapport à l'altitude à laquelle ils ont été générés. Du fait de cette configuration des tourbillons, il est intéressant, pour un aéronef suiveur qui suit l'aéronef meneur générant les tourbillons, de pouvoir profiter des vents ascendants afin de réduire sa consommation de carburant. EXPOSE DE L'INVENTION Un objet de la présente invention est de proposer un procédé de contrôle qui permet d'estimer la position des tourbillons générés par un aéronef meneur afin de contrôler la trajectoire d'un aéronef suiveur dans lequel est implanté un tel procédé de contrôle afin de pouvoir bénéficier des vents ascendants des tourbillons générés par l'aéronef meneur. A cet effet, est proposé un procédé de contrôle de la trajectoire d'un aéronef, dit aéronef suiveur présentant un processeur et un système de navigation, par rapport à la trajectoire d'un aéronef, dit aéronef meneur situé devant l'aéronef suiveur et générant des tourbillons, le procédé de contrôle comportant : 3041121 2 - une étape de recueil au cours de laquelle le processeur recueille des données relatives à l'aéronef meneur et à l'aéronef suiveur, ces données comprenant des données relatives à l'environnement atmosphérique, - une étape de calcul au cours de laquelle le processeur calcule, à partir des 5 données ainsi recueillies, un écart entre une position actuelle de l'aéronef suiveur et une position à atteindre pour l'aéronef suiveur, la position à atteindre correspondant à une position prédéterminée où l'aéronef suiveur utilise au mieux l'effet des tourbillons, - une étape de détermination au cours de laquelle le système de navigation to détermine des commandes à appliquer à des gouvernes de l'aéronef suiveur pour annuler l'écart ainsi calculé, - une étape de commande au cours de laquelle le système de navigation commande les gouvernes de l'aéronef suiveur selon les commandes ainsi déterminées, et 15 - un bouclage sur l'étape de calcul. Un tel procédé permet ainsi à l'aéronef suiveur de maintenir sa position en fonction des tourbillons générés par l'aéronef meneur afin de pouvoir bénéficier des vents ascendants. Avantageusement, l'étape de calcul consiste en : 20 - une première étape au cours de laquelle le processeur détermine des paramètres des tourbillons générés par l'aéronef meneur, - une deuxième étape au cours de laquelle le processeur détermine, à partir au moins des paramètres ainsi calculés, les positions des centres des tourbillons depuis leur génération par l'aéronef meneur, 25 - une troisième étape au cours de laquelle le processeur détermine la position du centre de tourbillon le plus proche de la position actuelle de l'aéronef suiveur, - une quatrième étape au cours de laquelle le processeur détermine une position cible d'un centre du tourbillon le plus proche, ainsi qu'un décalage cible par rapport à cette position cible, la position cible correspondant à une position d'un centre du 30 tourbillon proche en amont de l'aéronef suiveur pour un âge de tourbillon cible prédéterminé et le décalage cible étant un décalage entre l'aéronef suiveur et le centre du tourbillon proche et en amont, 3041121 3 - une cinquième étape au cours de laquelle le processeur détermine la position à atteindre pour l'aéronef suiveur par sommation vectorielle de la position cible du point à l'âge cible sur les centres du tourbillon proche, et le décalage cible, - une sixième étape au cours de laquelle le processeur calcule l'écart 5 correspondant à la différence entre la position à atteindre ainsi calculée et la position actuelle de l'aéronef suiveur. L'invention propose également un aéronef, dit aéronef suiveur comportant un système de navigation et un processeur, le processeur est configuré pour : - recueillir des données relatives à un aéronef, dit aéronef meneur situé devant 10 l'aéronef suiveur et générant des tourbillons, ainsi que des données relatives à l'aéronef suiveur, ces données comprenant des données relatives à l'environnement atmosphérique, - calculer, à partir des données ainsi recueillies, un écart entre une position actuelle de l'aéronef suiveur et une position à atteindre pour l'aéronef suiveur, la 15 position à atteindre correspondant à une position où l'aéronef suiveur utilise au mieux l'effet des tourbillons, - transmettre l'écart au système de navigation, et le système de navigation est configuré pour : - recevoir l'écart ainsi transmis, 20 - déterminer des commandes à appliquer à des gouvernes de l'aéronef suiveur pour annuler l'écart, - commander les gouvernes de l'aéronef suiveur selon les commandes ainsi déterminées. Avantageusement, le processeur est configuré pour : 25 - déterminer des paramètres des tourbillons générés par l'aéronef meneur, - déterminer, à partir au moins des paramètres ainsi calculés, les positions des centres des tourbillons depuis leur génération par l'aéronef meneur, - déterminer la position du centre de tourbillon le plus proche de la position actuelle de l'aéronef suiveur, 30 - déterminer une position cible d'un centre du tourbillon le plus proche, ainsi qu'un décalage cible par rapport à cette position cible, la position cible correspondant à une position d'un centre du tourbillon proche en amont de l'aéronef suiveur pour un âge de tourbillon cible prédéterminé et le décalage cible étant un décalage entre l'aéronef suiveur et le centre du tourbillon proche et en amont, 3041121 4 - déterminer la position à atteindre pour l'aéronef suiveur par sommation vectorielle de la position cible du point à l'âge cible sur les centres du tourbillon proche, et le décalage cible, - calculer l'écart correspondant à la différence entre la position à atteindre ainsi 5 calculée et la position actuelle de l'aéronef suiveur. BREVE DESCRIPTION DES DESSINS Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, apparaîtront plus clairement à la lecture de la description suivante d'un exemple de réalisation, ladite description étant faite en relation avec les dessins joints, parmi 10 lesquels : la Fig. 1 montre une vue de dessus d'un exemple de navigation entre un aéronef meneur et un aéronef suiveur selon l'invention, la Fig. 2 montre l'impact du roulis d'un aéronef meneur sur la position initiale du centre de chaque tourbillon, 15 la Fig. 3 montre l'impact du roulis d'un aéronef meneur sur la direction de la vitesse de descente du centre de chaque tourbillon à cause de l'interaction mutuelle entre les deux tourbillons, la Fig. 4 représente l'architecture d'un système de contrôle d'un aéronef mettant en oeuvre un procédé de contrôle selon l'invention, 20 la Fig. 5 est un organigramme d'un procédé de contrôle selon l'invention, et la Fig. 6 est un organigramme d'une étape de calcul de l'organigramme de contrôle de la Fig. 5. EXPOSE DETAILLE DE MODES DE REALISATION La Fig. 1 montre en vue de dessus un aéronef meneur 10 avec deux ailes 12a-b 25 qui est en vol et qui génère au niveau de chacune des ailes 12a-b, un tourbillon 14a-b (respectivement bâbord - tribord) représenté par son centre. La Fig. 1 montre également un aéronef suiveur 20 selon l'invention qui suit l'aéronef meneur 10 et qui comporte un système de contrôle 50 permettant la mise en oeuvre d'un procédé de contrôle qui est décrit ci-après et qui permet entre autres de 30 calculer, en fonction du temps depuis lequel il a été généré, la position du centre de chaque tourbillon 14a-b généré par l'aéronef meneur 10 afin de pouvoir contrôler la 3041121 5 position de l'aéronef suiveur 20 en fonction de la position ainsi calculée du centre de chaque tourbillon 14a-b et ainsi de profiter des vents ascendants des tourbillons 14a-b. Les positions de l'aéronef meneur 10 et de l'aéronef suiveur 20 sont, de préférence, déterminées dans un repère géographique (latitude, longitude, altitude) et 5 relevées par le système de navigation de l'aéronef correspondant, mais un autre repère peut également être utilisé et il faut alors mettre en oeuvre des changements de repère. Localement et par convention, on appelle X l'axe longitudinal du fuselage de l'aéronef 10, 20 orienté positivement dans le sens d'avancement de l'aéronef et passant par le centre de gravité de l'aéronef 10, 20, on appelle Y l'axe horizontal qui passe par 10 le centre de gravité de l'aéronef 10, 20 et qui est horizontal, et Z l'axe vertical qui passe par le centre de gravité de l'aéronef 10, 20 et qui est vertical, ces trois directions X, Y et Z étant orthogonales entre elles et formant un repère orthonormé ayant pour origine le centre de gravité O. La Fig. 1 montre également la trajectoire géographique 30 de l'aéronef meneur 15 10. Cette trajectoire est définie par un nombre fini de positions de l'aéronef Mo-M4, ces positions étant par exemple transmises par l'aéronef meneur 10 vers l'aéronef suiveur 20 et stockées dans le système de contrôle 50 afin d'obtenir un historique des positions. La position avec l'indice le plus petit est la position la plus récente. La Fig. 1 montre également le vecteur 40 représentatif du vent, qui est variable, 20 en ce sens que la direction et la vitesse du vent ne sont pas nécessairement constantes mais peuvent, par exemple, varier selon la position Mo-M4 de l'aéronef meneur 10. La Fig. 1 montre également l'aéronef suiveur 20 dans une « position désirée » référencée 25 et qui est la position à atteindre pour l'aéronef suiveur 20 pour qu'il bénéficie au maximum des effets des tourbillons 14a-b. 25 La Fig. 2 et la Fig. 3 montrent l'aéronef meneur 10 vu depuis son arrière et dans un repère fixé sur le centre de gravité de l'aéronef meneur 10. Les Figs. 2 et 3 montrent les positions initiales simplifiées des deux centres des tourbillons 14a-b avec l'aéronef meneur 10 en roulis avec un angle de roulis égal à I. La Fig. 2 montre les positions initiales ZBO,M, YBO,M, ZTO,M et YTO,M des deux centres des tourbillons 14a-b 30 dans ce repère. La Fig. 3 montre, par l'exemple du centre du tourbillon bâbord 14a, la décomposition de la vitesse de descente W, en fonction de l'angle de roulis (tem de l'aéronef meneur 10 de ce centre de tourbillon bâbord 14a à cause de l'interaction mutuelle des tourbillons 14a-b. 3041121 6 En référence avec la Fig. 4, le système de contrôle 50 comprend, reliés par un bus de communication 51: un processeur 52 ou CPU (« Central Processing Unit » en anglais), une mémoire vive RAM 53 (« Random Access Memory » en anglais), une mémoire morte ROM 54 (« Read Only Memory » en anglais), une unité de stockage 5 55 telle qu'un disque dur ou un lecteur de support de stockage, au moins une interface de communication 56, permettant au système de contrôle 50 de communiquer avec des modules extérieurs, tels que des modules de commande des gouvernes de l'aéronef suiveur 20, des caméras, des modules de relevés d'informations comme par exemple l'altitude, la vitesse, la température extérieure, les angles d'orientation de 10 l'aéronef.. Le processeur 52 est capable d'exécuter des instructions chargées dans la RAM 53 à partir de la ROM 54, de l'unité de stockage 55, d'une mémoire externe (non représentée), d'un support de stockage (tel qu'une carte SD), ou d'un réseau de communication. Lorsque l'équipement est mis sous tension, le processeur 52 est 15 capable de lire de la RAM 53 des instructions et de les exécuter. Ces instructions forment un programme d'ordinateur causant la mise en oeuvre, par le processeur 52, de tout ou partie des algorithmes et étapes du procédé de contrôle décrits ci-après. Tout ou partie des algorithmes et étapes du procédé de contrôle décrits ci-après peuvent être implémentés sous forme logicielle par exécution d'un ensemble 20 d'instructions par une machine programmable, par exemple un DSP (« Digital Signal Processor » en anglais) ou un microcontrôleur, ou être implémenté sous forme matérielle par une machine ou un composant dédié, par exemple un FPGA (« Field-Programmable Gate Array » en anglais) ou un ASIC (« Application-Specific Integrated Circuit » en anglais). 25 La Fig. 5 montre les étapes d'un procédé de contrôle 500 du vol de l'aéronef suiveur 20 par rapport au centre de l'un des tourbillons 14a-b générés par l'avion meneur 10 et qui comporte : - une étape de recueil 502 au cours de laquelle le processeur 52 recueille des données relatives à l'aéronef meneur 10 et à l'aéronef suiveur 20, ces données 30 comprenant des données relatives à l'environnement atmosphérique, - une étape de calcul 504 au cours de laquelle le processeur 52 calcule, à partir des données ainsi recueillies, un écart entre une position actuelle de l'aéronef suiveur 20 et une position à atteindre pour l'aéronef suiveur 20, la position à atteindre correspondant à une position prédéterminée que l'aéronef suiveur 20 doit occuper 3041121 7 pour utiliser au mieux l'effet des tourbillons 14a-b, cette position à atteindre est également appelée position désirée 25, - une étape de détermination 506 au cours de laquelle un système de navigation de l'aéronef suiveur 20 détermine les commandes à appliquer aux gouvernes de 5 l'aéronef suiveur 20 pour annuler l'écart ainsi calculé et ainsi approcher l'aéronef suiveur 20 de la position désirée 25, et - une étape de commande 508 au cours de laquelle le système de navigation commande les gouvernes de l'aéronef suiveur 20 selon les commandes ainsi déterminées. 10 Le but est d'utiliser au mieux le tourbillon 14a-b, ce qui correspond globalement à une conciliation entre une réduction de la consommation de carburant visée et les impacts non-désirés de la turbulence de sillage sur l'aéronef suiveur 20 et pouvant créer un inconfort des passagers. Après l'étape de commande 508, le procédé de contrôle 500 boucle (510) sur 15 l'étape de calcul 504 pour garantir que l'aéronef suiveur 20 suive toujours une trajectoire optimisée. De base, l'échantillonnage du procédé de contrôle 500 est adapté au taux de réception des données communiquées par l'aéronef meneur 10 à l'aéronef suiveur 20. Lorsque ce taux est trop faible et/ou irrégulier et/ou que les données sont reçues avec 20 un retard trop élevé et/ou que les données sont reçues avec trop de bruit, toutes ou une partie des données réceptionnées par l'aéronef suiveur 20 peuvent être soumises, par le processeur 52, à un échantillonnage à un taux constant et prédéfini, une correction, une stabilisation et une extrapolation, qui sont des traitements de données connues de l'homme du métier. 25 L'aéronef suiveur 20 comporte ainsi un processeur 52 configuré pour mettre en oeuvre l'étape de recueil 502 et l'étape de calcul 504 ainsi qu'à transmettre l'écart au système de navigation qui est configuré pour recevoir l'écart ainsi transmis, et mettre en oeuvre l'étape de détermination 506 et l'étape de commande 508. Le système de navigation et le processeur peuvent être deux éléments distincts 30 ou être combinés en un seul composant. Les étapes qui consistent à déterminer 506 les commandes et à commander 508 les gouvernes sont des étapes qui sont classiquement mises en oeuvre dans un aéronef. Par exemple, lorsqu'une destination à atteindre est entrée dans son système de navigation, le système de navigation calcule la différence entre la destination à 3041121 8 atteindre et la position actuelle et la trajectoire à suivre pour atteindre cette destination et donc annuler la différence, ainsi que les commandes à appliquer aux gouvernes pour diriger l'aéronef. Les données utilisées dans le procédé de contrôle 500 concernent des données 5 relatives à l'aéronef suiveur 20 comme sa position et son orientation : altitude, latitude, longitude, cap, angle de roulis, son envergure. Les données concernant l'aéronef suiveur 20 sont déjà présentes sur le système de navigation de l'aéronef suiveur 20 et sont relevées par des capteurs déjà présents dans l'aéronef suiveur 20 ou sont présentes dans son système de navigation. 10 Les données utilisées dans le procédé de contrôle 500 concernent également des données relatives à l'aéronef meneur 10 comme par exemple sa position (altitude, latitude, longitude), son orientation (cap, angle de roulis) ainsi que sa masse, la vitesse air, le facteur de charge vertical, son envergure. Les données utilisées dans le procédé de contrôle 500 concernent également des 15 données relatives à l'environnement atmosphérique comme par exemple la température statique, l'altitude-pression (pour calculer la densité de l'air au point de vol), et la vitesse du vent. Ces données sont mesurées par l'aéronef meneur 10 et/ou par l'aéronef suiveur 20. Concernant la vitesse du vent et en particulier les composantes horizontales du vent suivant les axes Nord-Sud et Est-Ouest sont à 20 recueillir, ce qui est facilité par l'utilisation de coordonnées géographiques relevées par le système de navigation. Les données utilisées dans le procédé de contrôle 500 concernent également des temps de réception par l'aéronef suiveur 20 et le temps d'échantillonnage du modèle. Les données sont datées afin d'avoir une référence commune aux deux aéronefs 10 et 25 20. Les différentes données concernant l'aéronef meneur 10 peuvent être : - estimées et/ou mesurées à bord de l'aéronef suiveur 20 par exemple à l'aide de moyens optiques tels que des caméras disposées à bord de l'aéronef suiveur 20 et orientées vers l'aéronef meneur 10 pour capturer des images de l'aéronef meneur 10, 30 et d'une unité de calcul calculant l'envergure par exemple par mesure directe de l'écart entre deux points distinctifs de la géométrie de l'appareil sur l'image (par exemple, l'écart entre deux moteurs) et l'angle de roulis cl) de l'aéronef meneur 10 à partir de l'analyse d'images capturées par les moyens optiques, et/ou 3041121 9 - estimées et/ou mesurées à bord de l'aéronef suiveur 20 par exemple à l'aide d'estimateurs de masse fonction de la classe de l'aéronef meneur 10 et d'une consommation estimée, et/ou - acquises par un système de communication du type ADS-B (position, temps de 5 réception...) depuis des données transmises par l'aéronef meneur 10, ou tout autre type de communication, comme par exemple : radio, data link dédié etc, et/ou - mesurées à bord de l'aéronef suiveur 20 en faisant l'hypothèse que les écarts entre l'aéronef meneur 10 et l'aéronef suiveur 20 sont négligeables, et qu'ils ont donc le même point de vol (par exemple la température statique, l'altitude-pression pour la 10 densité de l'air, les caractéristiques du vent...), et/ou - acquises auprès d'une base de données fonction du type d'avion qui peut être communiqué par l'aéronef meneur 10. L'angle de roulis (tem de l'aéronef meneur 10 peut être envoyé par l'aéronef meneur 10 à l'aéronef suiveur 20 par tout type de communication ou être estimé avec 15 d'autres paramètres comme le rayon R de la trajectoire horizontale de l'aéronef meneur 10, lui-même dérivé des positions historiques et actuelles de l'aéronef meneur 10, l'accélération de la pesanteur g et la vitesse air V de l'aéronef meneur 10 avec une formule comme (I) = arctan V2 R.g Le facteur de charge vertical nz de l'aéronef meneur 10 peut être estimé de 20 plusieurs manières : - soit mesuré et envoyé par l'aéronef meneur 10 à l'aéronef suiveur 20 par tout type de communication, - soit estimé avec d'autres paramètres comme l'angle de roulis en virage stabilisé (tem par exemple par la formule nz - 1 cos (1),, 25 - soit égal de manière forfaitaire à une constante (par exemple 1). La densité de l'air au point de vol de l'aéronef meneur 10 peut être estimée de plusieurs manières : - soit mesurée et envoyée par l'aéronef meneur 10 à l'aéronef suiveur 20 par tout type de communication, 30 - soit mesurée directement par le suiveur en faisant l'hypothèse que la différence entre la densité de l'air au point de vol de l'aéronef suiveur 20 et la densité de l'air au point de vol de l'aéronef meneur 10 est négligeable. Dans ce cas, l'aéronef suiveur 20 3041121 10 calcule directement cette donnée à partir de mesures de température statique et de pression statique. Les vents horizontaux (deux composantes) et verticaux rencontrés par l'aéronef meneur 10 peuvent être estimés de plusieurs manières : 5 - soit mesurés et envoyés par l'aéronef meneur 10 à l'aéronef suiveur 20 par tout type de communication, - soit calculés sur l'aéronef suiveur 20 en utilisant des données envoyées par l'aéronef meneur 10, par exemple les vecteurs des vitesses sol et air, par tout type de communication, 10 - soit mesurés directement au niveau de l'aéronef suiveur 20 par l'aéronef suiveur 20 lui-même en faisant l'hypothèse que la différence entre les vents rencontrés par l'aéronef suiveur 20 et les vents rencontrés par l'aéronef meneur 10 est négligeable. Dans ce cas, l'aéronef suiveur 20 calcule directement les vents rencontrés à partir de sa vitesse sol et sa vitesse air : 15 - soit mesurées par d'autres aéronefs passés au même endroit ou ses alentours, stockées dans une base de données au sol et envoyées par une station au sol vers l'aéronef suiveur 20, - soit par une combinaison de plusieurs des manières listées ci-dessus. La vitesse air de l'aéronef meneur 10 peut être estimée de plusieurs manières : 20 - soit mesurée et envoyée par l'aéronef meneur 10 à l'aéronef suiveur 20 par tout type de communication, - soit reconstruite à partir de la vitesse sol et des vents rencontrés par l'aéronef meneur 10, la vitesse sol étant envoyée par tout type de communication et les vents établis comme décrit ci-dessus. 25 La Fig. 6 détaille les calculs effectués au cours de l'étape de calcul 504. Les calculs sont effectués par rapport à un échantillonnage qui doit être suffisant pour pouvoir calculer un nombre suffisant de positions des centres des tourbillons 14a-b pour pouvoir évaluer ensuite une trajectoire à suivre pour l'aéronef suiveur 20. En effet, si le nombre est trop faible, les positions des centres des tourbillons 14a-b 30 calculées, c'est-à-dire les points PBO-PB4, PTO-PT4, seront très distantes les unes des autres et il sera difficile d'extrapoler une trajectoire à suivre pour l'aéronef suiveur 20. Un temps d'échantillonnage particulier est 0,2 secondes, l'échantillonnage correspond au bouclage 510. 3041121 11 A partir des données d'entrées recueillies au cours de l'étape de recueil 502, le principe du calcul consiste donc, à chaque instant de l'échantillonnage, à établir et à maintenir un historique des positions Mo-M4 de l'aéronef meneur 10 avec, en plus, un nombre des données - parmi elles le temps - associées à ces positions qui sont requises 5 ensuite pour le calcul des positions actuelles des centres des tourbillons 14a-b. Au minimum, à chaque instant de l'échantillonnage, les positions Mo-M4 de l'aéronef meneur 10 établies et stockées pendant la boucle précédente sont décalées pour que l'ancien Mo devienne le nouveau Mi, l'ancien Mi devient le nouveau M2 et ainsi de suite. Aussi, dans tous les cas, une nouvelle position Mo de l'aéronef meneur 10, dite 10 actuelle, avec les données associées doit être établie. Les données associées aux positions décalées sont soit décalées de la même façon soit recalculées. De nouveaux calculs peuvent-être privilégiés dans les cas où des améliorations de la qualité de ces données sont possibles et souhaitées, par exemple si le vent n'est pas transmis par l'aéronef meneur 10 mais établi avec des mesures faites sur l'aéronef suiveur 20. Afin 15 de ne pas excéder la mémoire disponible, les plus anciennes positions et les données associées peuvent être effacées après ce décalage des positions et données associées. Au cours d'une première étape 602, le processeur 52 détermine des paramètres des tourbillons 14a-b générés par l'aéronef meneur 10 que sont l'espacement 1), entre les centres des deux tourbillons 14a-b et la vitesse de descente \V, des centres des 20 tourbillons 14a-b, associés avec la position actuelle Mo et l'état de l'aéronef meneur 10. Pour ce faire, le processeur 52 détermine à partir de données 605 relatives à la position de l'aéronef meneur 10, la position Mo, dite position actuelle de l'aéronef meneur 10. Le cas échéant, le processeur 52 calcule aussi un rétablissement des 25 positions Mi-M4 de l'aéronef meneur 10, dites positions anciennes, par exemple par une mise à jour d'une interpolation parmi plusieurs positions réceptionnées. Le processeur 52 détermine ensuite, à partir de données 601 relatives à l'aéronef meneur 10 et de données 603 relatives à l'environnement atmosphérique, les paramètres des tourbillons 14a-b. 30 L'espacement 1), (en mètres) entre les centres des deux tourbillons 14a-b qui sont générés par l'aéronef meneur 10 dépend de son envergure 1)m et peut être calculé par exemple par la relation : b =- 4 *b (1) où : 1)m est l'envergure de l'aéronef meneur 10 (en mètres). 3041121 12 bol est issu des données transmises par l'aéronef meneur 10, ou calculé par un des autres moyens préalablement décrits. Pour calculer bv, d'autres solutions sont envisageables. Par exemple bv peut être envoyé par l'aéronef meneur 10 à l'aéronef suiveur 20, ou sur l'aéronef meneur 10, bv 5 peut être stocké comme une constante en fonction du type de l'aéronef meneur 10, ou être calculé par une fonction plus élaborée prenant en compte la position des volets et/ou la position du centre de gravité et/ou la portance, normalisée ou absolue. bv peut également être stocké dans une base des données sur l'aéronef suiveur 20 comme fonction du type d'aéronef meneur 10, où ce type d'aéronef est envoyé par 10 l'aéronef meneur 10 à l'aéronef suiveur 20. Dans les cas précédents, les différentes valeurs de bv sont soit connues par construction, soit identifiées à l'aide de valeurs expérimentales (par exemple : souffleries, méthodes de mécanique des fluides numérique, mesures effectuées en vol). 15 La vitesse de descente Wv des centres des tourbillons 14a-b dépend des données de l'aéronef meneur 10 et elle peut être calculée par la relation : m.g.nz w (2) où : m est la masse de l'aéronef meneur 10 (en kg), g est l'accélération de la pesanteur à l'altitude considérée (en m/s2), 20 p est la densité de l'air au point de vol de l'aéronef meneur 10 (en kg.m1), V est la vitesse air vraie de l'aéronef meneur 10 (en m.s-1), bv est l'espacement entre les centres des deux tourbillons 14a-b (en m), et n, est le facteur de charge vertical subi par l'aéronef meneur 10 (en g). Le cas échéant, le processeur 52 calcule également un rétablissement des 25 espacements bv entre les centres des deux tourbillons 14a-b et des vitesses de descente Wv des centres des tourbillons 14a-b associés avec les positions anciennes Mi-M4 de l'aéronef meneur 10. Un tel rétablissement peut être pertinent si la vitesse air vraie et/ou le facteur de charge et/ou autres données de l'aéronef meneur 10 sont établis par une interpolation. Dans tous les cas, le processeur 52 stocke ces paramètres des 30 tourbillons 14a-b qui sont liés aux positions Mo-M4 de l'aéronef meneur 10 associées. Dans le cadre d'une deuxième étape 604, le processeur 52 détermine ensuite, à partir de données 601 relatives à l'aéronef meneur 10 et de données 603 relatives à l'environnement atmosphérique, d'autres paramètres comme l'angle de roulis (1)m et l'angle de cap kPol de l'aéronef meneur 10, le temps T, et les composants de la vitesse 3041121 13 du vent (W,, Wy et Wz). Ces données sont établies au moins pour la position actuelle Mo de l'aéronef meneur 10. Le cas échéant, un ou plusieurs de ces paramètres qui sont associés avec des positions anciennes Mi-M4 de l'aéronef meneur 10 sont rétablis et stockés en association avec les positions anciennes correspondantes. 5 Au cours de la deuxième étape 604, le processeur 52 détermine alors, à partir des paramètres calculés et stockés ci-dessus, les positions initiales verticales et latérales des centres des tourbillons 14a-b pour chaque position Mo-M4 de l'aéronef meneur 10 préalablement établie, au moins celles associées avec la position Mo la plus récente. Ces positions initiales sont établies à partir des positions, des orientations, des 10 vitesses et autres paramètres de l'aéronef meneur 10 associées avec ces positions Mo- M4. Par exemple, sur la Fig. 1, les positions initiales associées avec la position actuelle Mo de l'aéronef meneur 10 sont les points PBO et PTO. Sur la Fig. 2, la position PBO est indiquée par ses deux composantes de position ZBO,M et YBO,M dans un repère 15 orthogonal, signalé YM, ZM et associé à l'aéronef meneur 10, et la position du point PTO est indiquée par ses deux composantes de position ZTO,M et YTO,M dans le même repère. Les positions initiales verticales et latérales (horizontales) dans le repère géodésique dépendent alors de la position du centre de gravité de l'aéronef meneur 10 20 dans ce repère, de son angle de roulis (tem, de son cap kPoi et de l'espacement 1), entre les deux centres de tourbillons 14a-b. En utilisant un repère géodésique avec des coordonnées en latitude, longitude et altitude, les positions initiales verticales et latérales peuvent être calculées par les équations : {YBO,i { 25 YTO,i XBO,i = Xcgm,i = Ycgm,i 0.5 ZBO,i = XTO,i = Xcgm,i = + 0.5 - ZTO,i + 0.5 - b',,, - cos(Dm,, - sinTm,, - M2D - b',, - cos(1)m, - cosTm,, - M2D /cosXcgm,, Zcgm,, + 0.5 - b',, - sinelm,, (3) - 0.5 - b',,, - cos(1)m,, - sinTm,, - M2D b',, - cos(1)m,, - cosTm,,- M2D /cosXcgm,, Zcgm,i - 0.5 - b',,, - sinelm,, (4) où: - i est l'indice associé avec une des positions Mo-M4 de l'aéronef meneur 10, par exemple i = 1 est associé avec la position Mi, i = 2 avec M2 et ainsi de suite, - XBO est la latitude initiale du centre du tourbillon 14a situé à bâbord (en degré), 30 - Xcgm est la latitude du centre de gravité de l'aéronef meneur 10 (en degré), -1), est l'espacement entre les centres des deux tourbillons 14a-b (en m), 3041121 14 - (1)m est l'angle de roulis de l'aéronef meneur 10 qui est compté positif aile droite en bas par convention (en radian), - Yol est le cap de l'aéronef meneur 10 qui est compté positif vers l'est par convention (en radian), 5 - M2D est la conversion entre mètre et degré, égale à 1/60/1852 = 8.9.10', - YBO est la longitude initiale du centre du tourbillon 14a situé à bâbord (en degré), - Ycgm est la longitude du centre de gravité de l'aéronef meneur 10 (en degré), - ZBO est la position verticale initiale du centre du tourbillon 14a situé à bâbord 10 (en m), - Zcgm est la position verticale du centre de gravité de l'aéronef meneur 10 (en m), - XTO est la latitude initiale du centre du tourbillon 14b situé à tribord (en m) - YTO est la longitude initiale du centre du tourbillon 14b situé à tribord (en m), 15 et - ZTO la position verticale initiale du centre du tourbillon 14b situé à tribord (en m). Les formules (3) et (4) incluent une transformation entre un repère orthogonal centré sur l'aéronef meneur 10 et le repère géodésique. Dans les formules (3) et (4), 20 cette transformation est de forme simplifiée et pourrait être remplacée par une autre transformation. Aussi, les positions peuvent être formulées dans d'autres types des repères en utilisant d'autres transformations sans que la fonctionnalité et le principe de l'invention ne soient altérés. Les formules (3) et (4) s'appliquent également à tout rétablissement des 25 positions initiales associées avec les positions anciennes Mi-M4 de l'aéronef meneur 10, comme indiqué par l'utilisation d'un indice i. Dans ce cas, les données d'entrée sont celles liées à ces points et comme préalablement établi dans la première étape 602. Toujours au cours de la deuxième étape 604, le processeur 52 calcule les 30 positions actuelles des centres des tourbillons 14a-b à partir des paramètres calculés et stockés pendant la première étape 602 et/ou les boucles précédentes et les positions initiales des centres des tourbillons 14a-b pour toutes les positions Mo-M4 de l'aéronef meneur 10 calculées comme décrit ci-dessus. D'une manière générale, au cours de la deuxième étape 604, le processeur 52 détermine les positions des centres des 3041121 15 tourbillons 14a-b depuis leur génération par l'aéronef meneur 10. C'est-à-dire la position de chaque tourbillon 14a-b selon son âge. Les positions actuelles de tous les centres des tourbillons 14a-b sont déterminées à partir de leurs positions initiales préalablement calculées en prenant en compte leurs 5 déplacements en vertu de l'influence du vent et de l'interaction aérodynamique mutuelle entre les deux tourbillons 14a-b. Les distances de ces déplacements dépendent des âges des centres des tourbillons 14a-b, où l'âge d'un centre d'un tourbillon 14a-b est le temps qui s'est écoulé depuis la génération du centre sur l'aéronef meneur 10. Ainsi sur la Fig. 1, le point PB4 représente le centre du tourbillon 10 bâbord 14b qui est le plus ancien et le point PBO représente le centre du tourbillon bâbord 14b qui est le plus récent. La même représentation est faite avec les points PT4 et PTO pour le tourbillon tribord 14a. C'est pendant ces âges que les vents et l'interaction mutuelle agissent sur les positions des centres des tourbillons 14a-b. Ces âges correspondent alors au temps qui 15 est passé au moment du calcul depuis que l'aéronef meneur 10 était aux positions correspondantes. Ce calcul se fait pour toutes les positions Mi-M4, et donne lieu à un nombre de positions des centres des tourbillons 14a-b actuels Psi à PB4 et PT1 à PT4, où, par exemple, les centres associés avec les positions PB3 et PT3 ont été créés sur l'aéronef meneur 10 quand cet aéronef était à la position M3. 20 Sur la Fig. 1, les positions des centres PBO-PB4 et PTO-PT4 des tourbillons 14a-b sont toutes des positions actuelles établies par ce processus et ensemble elles définissent les positions actuelles des deux centres des tourbillons 14a-b comme cela est montré sur la Fig. 1. Toutes les paires des centres individuels à l'instant actuel, comme par exemple PB1 et PT1 ou PB2 et PT2 ou encore PB3 et PT3, sont associées avec 25 des âges différents. Ces âges augmentent avec la distance entre les centres et la position actuelle Mo de l'aéronef meneur 10. Par exemple l'âge de PB3 est supérieur de celui de PB2 et la différence est égale au temps passé entre les moments où l'aéronef meneur 10 était aux positions M3 et M2. Les positions des centres de chaque tourbillon 14a-b à l'instant actuel ainsi 30 établies au cours de la deuxième étape de calcul 604 par le processeur 52 prennent en compte les variations dynamiques des données pertinentes, c'est-à-dire les variations influençant les déplacements et les positions des centres Psi-PB4 et PTi-PT4 de chaque tourbillon 14a-b, comme, par exemple, les variations d'état de vol de l'aéronef meneur 3041121 16 10 et les variations des conditions de l'atmosphère, rencontrés par l'aéronef meneur 10. Les positions actuelles des centres PBO-PB4 et PTO-PT4 de chaque tourbillon 14a-b peuvent être décrites dans le repère géodésique avec trois composantes comme défini 5 en suivant : XB PB, [YB ZBi (5) XT TiP -[YT ZT I, (6) où - i est l'indice associé avec une des positions Mo-M4 de l'aéronef meneur 10, par exemple i = 1 est associé avec la position Mi, i = 2 avec M2 et ainsi de suite, 10 - XB et XT sont les positions en latitude géodésique (en degré) d'un centre de tourbillon 14a-b, - YB et YT sont les positions en longitude géodésique (en degré) d'un centre de tourbillon 14a-b et - ZB et ZT sont les positions en altitude géodésique (en m) d'un centre de 15 tourbillon 14a-b. Les positions actuelles des centres PBO-PB4 et PTO-PT4 de chaque tourbillon 14a-b peuvent être calculées par les formules suivantes : { YBi = YBO,i + (Wyi - + (Wxi + Wvi ' SiTlf1),,,i - sinTni,i) - (T0 - Ti) - M2D W',,i - sin(1),,,i - cosTni,i) - (T0 - Ti) - M2D /cosX,g,,,i XBi - XBO,i { XTi = XT + i, 0,i (Wx,i + Wi ' SiT1,(1),,,,,i ' Sint lim ,i) ' (T0 - Ti) - M2D 'Ti = YTO + (wy i - ,i Wv,i - SiTlft'in,i - COSTin,i) - (T0 - Ti) - M2D /cosXcgm,i 20 où : - i est l'indice associé avec une des positions Mo-M4 de l'aéronef meneur 10, par exemple i = 1 est associé avec la position Mi, i = 2 avec M2 et ainsi de suite, - XBO est la latitude initiale du centre du tourbillon 14a situé à bâbord (en degré), - Wx est la vitesse du vent vers le nord géodésique subie par les tourbillons 14a-b (en m. s-1), 25 - W, est la vitesse de descente des centres des tourbillons 14a-b à cause de leur interaction mutuelle (en m.s-1), - (tem est l'angle de roulis de l'aéronef meneur 10 au moment de la génération du centre du tourbillon 14a-b (en degré), ZBi = ZBO,i (Wz,iW'- COS(I)m,i) - (To - Ti) ZTi = ZTO,i (Wz,iW'- COS(I)m,i) - (T0 - Ti) (7) (8) 3041121 17 - Y. est le cap de l'aéronef meneur 10 qui est compté positif vers l'est par convention au moment de la génération du centre du tourbillon 14a-b (en radian), - M2D est la conversion entre mètre et degré, égale à 1/60/1852 = 8.9.10', - T est le temps associé avec une position de l'aéronef meneur 10 (en s), par 5 exemple To est le temps actuel et Ti est le temps associé avec le moment où l'aéronef meneur 10 était à la position Mi, - YBO est la longitude initiale du centre du tourbillon 14a situé à bâbord (en degré), - Wy est la vitesse du vent vers l'est géodésique subie par les tourbillons 10 (en m. s-1), - Xcgoi est la latitude du centre de gravité de l'aéronef meneur 10 (en degré), - ZBO est la position verticale initiale du centre du tourbillon 14a situé à bâbord (en m), - Wz est la vitesse du vent vertical subie par les tourbillons 14a-b (en m.s- 15 1), - XTO est la latitude initiale du centre du tourbillon 14b situé à tribord (en m) - YTO est la longitude initiale du centre du tourbillon 14b situé à tribord (en m), et - ZTO la position verticale initiale du centre du tourbillon 14b situé à tribord (en 20 m). Les formules (7) et (8) prennent en compte le fait que la direction de la vitesse de descente W, des centres des tourbillons 14a-b à cause de leur interaction mutuelle est inclinée vis-à-vis de la verticale géodésique dans le cas où l'aéronef meneur 10 vole en virage avec un angle de roulis I. La décomposition de cette vitesse dans une 25 composante horizontale et une composante verticale est montrée sur la Fig. 3. Les formules (7) et (8) sont basées sur l'hypothèse que la vitesse et la direction du vent restent constantes pendant le déplacement de chaque centre PBO-PB4, PTO-PT4 du tourbillon 14a-b et que la vitesse de descente à cause d'interactions mutuelles des tourbillons 14a-b reste également constante pendant le déplacement de chaque centre 30 PBO-PB4, PTO-PT4 du tourbillon 14a-b. Il est également possible de prendre en compte des données variables, par exemple en fonction de l'âge du tourbillon 14a-b. La transformation entre un repère orthogonal centré sur l'aéronef meneur 10 et le repère mondial géodésique incluse dans les formules (7) et (8) est de forme simplifiée et pourrait être remplacée par une autre transformation. Aussi, les distances 3041121 18 de déplacement peuvent être formulées dans d'autres types de repères, en utilisant d'autres transformations sans que la fonctionnalité et le principe de l'invention ne soient altérés. Les positions des centres de chaque tourbillon 14a-b à l'instant actuel ainsi 5 établies lors de la deuxième étape de calcul 604 peuvent être utilisées sur l'aéronef suiveur 20 dans une étape de calcul spécifique, par exemple pour des besoins de surveillance en aide au pilotage. Le cas échéant et à la fin de la deuxième étape de calcul 604, le processeur 52 efface l'historique enregistré, c'est à dire les données qui sont associées à des positions 10 de l'aéronef meneur 10 trop vieilles. Pour-cela, les dates des données historiques sont comparées avec un âge limite. Cet âge limite est défini soit comme constant, soit comme variable dynamique, par exemple en fonction de la distance actuelle entre l'aéronef meneur 10 et l'aéronef suiveur 20 ou en fonction de l'âge actuel du vortex le plus proche à l'aéronef suiveur 20 préalablement calculé par le processeur 52 dans une 15 troisième étape de calcul 606 d'une boucle précédente. Au cours d'une troisième étape 606, le processeur 52 calcule, à partir des positions des centres des tourbillons 14a-b en état actuel calculées ci-dessus et à partir de données 607 relatives à la position de l'aéronef suiveur 20, la position du centre de tourbillon 14a-b le plus proche de la position actuelle de l'aéronef suiveur 20. 20 A partir des coordonnées établies au cours de la deuxième étape de calcul 604 et enregistrées dans la base de données, le processeur 52 détermine d'abord quel centre de tourbillon actuel, bâbord 14a ou tribord 14b est le plus proche de la position de l'aéronef suiveur 20. Pour ce faire, il calcule pour chaque tourbillon 14a-b et pour chaque position stockée des centres dudit tourbillon 14a-b, la distance entre cette 25 position et la position de l'aéronef suiveur 20 et il détermine ainsi la distance minimale de l'aéronef suiveur 20 à chaque tourbillon 14a-b. De la plus petite distance de ces deux distances minimales, le processeur 52 détermine quel tourbillon bâbord 14a ou tribord 14b est le plus proche (ici 14b sur la Fig. 1). Le processeur 52 détermine alors la distance la plus courte entre le centre du 30 tourbillon le plus proche (ici 14b sur la Fig.1) et la position de l'aéronef suiveur 20. Pour ce faire, à partir des calculs précédents, il détermine les deux positions des centres du tourbillon 14b, parmi PTO-PT4, qui sont les plus proches de la position de l'aéronef suiveur 20. 3041121 19 Pour l'exemple, on considère que l'on obtient les deux positions PT2 et PT3 de la Fig. 1. Le processeur 52 calcule ensuite la distance OHT où HT est la projection orthogonale du point O sur le segment reliant PT2 et PT3. En appelant C le centre du tourbillon 14b qui est le plus proche, l'hypothèse OC=OHT est valide si la distance 5 CHT est suffisamment petite par rapport à OHT, ce qui est généralement vérifié si le rayon de courbure entre PT2 et PT3 est suffisamment grand ou de manière équivalente, si l'historique de position des centres des tourbillons 14a est suffisamment discrétisé pour que les points PT2 et PT3 soient suffisamment proches. La distance OHT est alors identifiée comme la distance minimale actuelle entre l'aéronef suiveur 20 et le centre 10 du tourbillon le plus proche 14b. Le point HT est assimilé au centre C du tourbillon 14b le plus proche. Le processeur 52 détermine ensuite l'âge exact du centre du tourbillon le plus proche. Suivant l'exemple ici, l'âge du tourbillon 14b à ce point HT est calculé à partir des données associées aux deux points PT2 et PT3 qui sont les deux points les plus 15 proches du tourbillon 14b le plus proche. On appelle l'âge associé avec ce point HT « l'âge du tourbillon le plus proche », ici abrévié en tHT. Par exemple l'âge du tourbillon le plus proche tHT peut être calculé par l'équation suivante qui établit une simple règle de trois et est basé sur la même hypothèse déjà appliquée pour le calcul de la position HT: 1T3 - T21 20 IPT3PT2I (9) où : - To est le temps actuel associé avec la position actuelle Mo de l'aéronef meneur 10 (en s), - T2 est le temps associé avec le moment où l'aéronef meneur 10 était à la position M2 (en s), 25 - T3 est le temps associé avec le moment où l'aéronef meneur 10 était à la position M3 (en s), - CPT2 est la distance entre C et PT2 (en m), et - PT2PT3 est la distance entre PT2 et PT3 (en m). Le processeur 52 détermine alors un repère dit « repère du tourbillon le plus 30 proche ». L'origine de ce repère est le point HT préalablement calculé. Selon un exemple, l'orientation géodésique de ce repère est définie par rapport à un repère orthogonal, avec une première coordonnée X, visant globalement l'aéronef meneur 10 et passant par le point PT2 suivant l'exemple de la Fig. 1, une deuxième coordonnée Y, passant par le point HB qui est la position du centre du tourbillon bâbord 14a à l'âge tHT = To - T2 + CPT2 3041121 20 du tourbillon le plus proche et visant globalement l'aéronef suiveur 20, et une troisième coordonnée Z, complétant un repère orthogonal. Le processeur 52 détermine ensuite la position du centre de gravité O de l'aéronef suiveur 20 dans les coordonnées du repère du tourbillon le plus proche par 5 un changement de repère classique. On appelle la position de O en coordonnée Y, la « distance latérale relative au tourbillon le plus proche » et celle en coordonnée Z, la « distance verticale relative au tourbillon le plus proche ». La distance latérale relative et la distance verticale relative définissent l'écart relatif entre la position actuelle de l'aéronef suiveur 20 et la position actuelle du centre du tourbillon 14b le plus proche. 10 Les distances latérales et verticales relatives au tourbillon le plus proche, ainsi que l'âge du tourbillon le plus proche et le côté, bâbord ou tribord, des centres de tourbillons 14a-b le plus proche établis par le calcul peuvent être utilisés sur l'aéronef suiveur 20 dans une étape spécifique, par exemple pour des besoins de surveillance en aide au pilotage ou de contrôle. 15 A partir des coordonnées des centres des tourbillons 14a-b à l'instant actuel, il est alors nécessaire d'évaluer la position que l'aéronef suiveur 20 doit avoir pour bénéficier au maximum des vents ascendants générés par les tourbillons 14a-b lorsqu'il survolera successivement ces centres. Comme déjà mentionné, on appelle « position désirée », la position que doit 20 atteindre l'aéronef suiveur 20 pour bénéficier au maximum des effets des tourbillons 14a-b. Sur la Fig. 1 l'aéronef suiveur 20 se trouve à la position désirée 25. On appelle « position de l'aéronef suiveur », la position de l'aéronef suiveur 20 au moment des calculs. Au cours d'une quatrième étape 608, le processeur 52 détermine, à partir de 25 données 607 relatives à la position de l'aéronef suiveur 20, des données 601 relatives à l'aéronef meneur 10, des données 603 relatives à l'environnement atmosphérique (conditions atmosphériques), une position cible d'un centre du tourbillon le plus proche qui correspond à une position pour laquelle l'aéronef suiveur 20 utilise au mieux les effets des tourbillons 14a-b générés par l'aéronef meneur 10, en particulier, 30 il calcule la position du centre du tourbillon pour laquelle l'aéronef suiveur 20 bénéficie au maximum des effets des tourbillons 14a-b. Cette position cible correspond à une position d'un centre du tourbillon le plus proche en amont de l'aéronef suiveur 20 pour un âge de tourbillon cible prédéterminé. Sur l'exemple donné par la Fig. 1 cette position cible correspond au point Gr. 3041121 21 Comme expliqué ci-après, la position cible de l'aéronef suiveur 20 est décalée par rapport à la position cible du centre du tourbillon proche. Ce décalage est appelé « décalage cible ». La position cible du centre d'un tourbillon proche peut être définie dans le 5 même repère que les positions des aéronefs, par exemple un repère géodésique. L'âge de tourbillon cible correspond temporellement à une distance de l'aéronef suiveur 20 derrière l'aéronef meneur 10, distance qui est considérée comme optimale en prenant en compte d'une part le bénéfice de performance, et d'autre part les besoins de distance de sécurité. L'âge cible peut être une constante, par exemple 5 s. Il 10 est également possible de prendre un âge cible variable qui est une fonction d'une ou plusieurs caractéristiques d'un tourbillon et/ou une fonction d'une ou plusieurs caractéristiques de l'aéronef meneur 10 et/ou une fonction d'une ou plusieurs caractéristiques de l'aéronef suiveur 20 et/ou une fonction d'une ou plusieurs caractéristiques de l'état de l'atmosphère. Le processeur 52 choisit donc un âge cible 15 constant ou variable. A partir de l'âge cible et en utilisant les positions des centres de chaque tourbillon 14a-b à l'instant actuel avec leurs âges associés, le processeur 52 détermine, pour chaque tourbillon 14a-b, un point, centre dudit tourbillon 14a-b à l'âge cible. Sur la Fig. 1, ces points sont GB et Gr. Ces points ne font pas nécessairement partie des 20 positions des centres déjà établies. Parmi ces deux points GB et Gr, seul le point correspondant au centre du tourbillon 14b le plus proche est considéré comme la position cible du centre, ici alors Gr. Les deux points GB et Gr à l'âge cible peuvent être identifiés par règle de trois à partir des données discrètes concernant les centres des tourbillons 14a-b amont et aval, 25 ici PB1 et PB2 sur le tourbillon bâbord 14a et PT' et PT2 sur le tourbillon tribord 14b, et les âges associés avec ces quatre points. Les deux points GB et GT définissent un repère tourbillon à l'endroit et à l'âge du tourbillon cible. Ce repère peut, selon un exemple, être défini de la même façon que le repère du tourbillon le plus proche, c'est-à-dire avec l'origine à Gr, une 30 première coordonnée X, visant principalement vers l'aéronef meneur 10 et passant par le point PT1 suivant l'exemple de la Fig. 1, une deuxième coordonnée Y, passant par le point GB et visant l'aéronef suiveur 20, et une troisième coordonnée Z, complétant un repère orthogonal. Ce repère est appelé « repère du tourbillon cible ». 3041121 22 Le décalage cible est défini dans ce repère. Vu que l'aéronef suiveur 20 est censé voler globalement en parallèle du tourbillon 14a-b il est suffisant de définir deux distances dans les axes perpendiculaires à l'axe principal du tourbillon 14a-b, alors dans les directions Y, et Le décalage cible prend en compte d'un côté le 5 bénéfice de performance et de l'autre côté les influences non-désirées sur l'aéronef suiveur 20. Le décalage cible peut être constant, par exemple 50 m en Y, et 10 m en Z. Il est également possible de prendre un décalage cible variable qui est une fonction d'une ou plusieurs caractéristiques d'un tourbillon et/ou une fonction d'une ou plusieurs caractéristiques de l'aéronef meneur 10 et/ou une fonction d'une ou 10 plusieurs caractéristiques de l'aéronef suiveur 20 et/ou une fonction d'une ou plusieurs caractéristiques de l'état de l'atmosphère. Par exemple le décalage cible en Y, peut être égal à la moitié de l'envergure de l'aéronef suiveur 20 plus 15% de l'envergure de l'aéronef meneur 10. Le processeur 52 choisit donc un décalage cible prédéterminé constant ou variable. 15 Au cours d'une cinquième étape 609, le processeur 52 calcule alors la position désirée que l'aéronef suiveur 20 doit prendre par rapport à la position cible ainsi déterminée, par sommation vectorielle de la position cible, ici le point GT à l'âge cible sur les centres du tourbillon proche 14b, et le décalage cible, en prenant en compte l'orientation et définition du repère. Cette position désirée est indiquée sur la Fig. 1 20 par la position désirée 25. Au cours d'une sixième étape 610, le processeur 52 calcule l'écart, qui est la différence entre la position désirée et la position actuelle de l'aéronef suiveur 20. Les lois de pilotage et les mouvements de l'aéronef suiveur 20 ainsi commandés tendent à faire converger cet écart vers 0 avec le temps, afin que la position de l'aéronef suiveur 25 20 coïncide avec la position désirée. L'écart répond à l'équation : Ecart = Posdesiree- Poss. (10). L'écart est l'information qui est fournie en entrée de l'étape 506 et qui va servir d'entrée aux lois de pilotage qui vont asservir l'aéronef suiveur 20 afin d'annuler cet 30 écart. L'étape de calcul 504 consiste ainsi en : - une première étape 602 au cours de laquelle le processeur 52 détermine les paramètres des tourbillons 14a-b générés par l'aéronef meneur 10, 3041121 23 - une deuxième étape 604 au cours de laquelle le processeur 52 détermine, à partir au moins des paramètres ainsi calculés, les positions des centres des tourbillons 14a-b depuis leur génération par l'aéronef meneur 10, - une troisième étape 606 au cours de laquelle le processeur 52 détermine la 5 position du centre de tourbillon le plus proche de la position actuelle de l'aéronef suiveur 20, - une quatrième étape 608 au cours de laquelle le processeur 52 détermine une position cible d'un centre du tourbillon le plus proche, ainsi qu'un décalage cible par rapport à cette position cible, la position cible correspondant à une position d'un centre 10 du tourbillon proche en amont de l'aéronef suiveur 20 pour un âge de tourbillon cible prédéterminé et le décalage cible étant un décalage entre l'aéronef suiveur 20 et le centre du tourbillon proche et en amont, - une cinquième étape 609 au cours de laquelle le processeur 52 détermine la position à atteindre pour l'aéronef suiveur 20 par sommation vectorielle de la position 15 cible du point GT ainsi déterminée à l'âge cible sur les centres du tourbillon proche 14b, et le décalage cible, - une sixième étape 610 au cours de laquelle le processeur 52 calcule l'écart correspondant à la différence entre la position à atteindre ainsi calculée et la position actuelle de l'aéronef suiveur 20. 20 Le processeur 52 est ainsi configure pour : - déterminer les paramètres des tourbillons 14a-b générés par l'aéronef meneur 10, - déterminer, à partir au moins des paramètres ainsi calculés, les positions des centres des tourbillons 14a-b depuis leur génération par l'aéronef meneur 10, 25 - déterminer la position du centre de tourbillon le plus proche de la position actuelle de l'aéronef suiveur 20, - déterminer une position cible d'un centre d'un tourbillon proche, ainsi qu'un décalage cible par rapport à cette position cible, la position cible correspondant à une position d'un centre du tourbillon proche en amont de l'aéronef suiveur 20 pour un âge 30 de tourbillon cible prédéterminé et le décalage cible étant un décalage entre l'aéronef suiveur 20 et le centre du tourbillon proche et en amont, - déterminer la position à atteindre pour l'aéronef suiveur 20 par sommation vectorielle de la position cible du point GT à l'âge cible sur les centres du tourbillon proche 14b, et le décalage cible, 3041121 24 - calculer l'écart correspondant à la différence entre la position à atteindre ainsi calculée et la position actuelle de l'aéronef suiveur 20. Le système de navigation et les lois de pilotage transforment l'écart en déflections gouvernes commandées (étape 506) afin de maintenir l'aéronef suiveur à 5 proximité des tourbillons 14a-b et avec pour objectif d'annuler l'écart. Le système de navigation met (étape 508) les gouvernes aux positions correspondant aux déflections commandées par les lois de pilotage, ce qui fait bouger l'aéronef suiveur 20 autour de son centre de gravité. L'étape de bouclage 510 permet de mettre à jour régulièrement les différents 10 calculs, en particulier d'emmagasiner des données relatives aux positions des centres des tourbillons en fonction du temps, afin que le système de navigation maintienne l'aéronef suiveur 20 en position pendant que l'aéronef meneur 10 et l'aéronef suiveur 20 se déplacent dans des conditions atmosphériques variables. L'aéronef suiveur 20 peut être équipé dans son cockpit, d'un écran de 15 visualisation, sur lequel sont affichés par exemple les positions des centres des tourbillons 14a-b actuels, la position relative entre les centres des tourbillons et l'aéronef suiveur 20, l'âge actuel du tourbillon le plus proche de l'aéronef suiveur 20, la position désirée, la position de l'aéronef suiveur 20, la position de l'aéronef meneur 10, les écarts entre les différentes positions. 20 REVENDICATIONS1) Procédé de contrôle (500) de la trajectoire d'un aéronef, dit aéronef suiveur (20) présentant un processeur (52) et un système de navigation, par rapport à la trajectoire d'un aéronef, dit aéronef meneur (10) situé devant l'aéronef suiveur et générant des tourbillons, le procédé de contrôle (500) comportant : - une étape de recueil (502) au cours de laquelle le processeur (52) recueille des données relatives à l'aéronef meneur (10) et à l'aéronef suiveur (20), ces données comprenant des données relatives à l'environnement atmosphérique, - une étape de calcul (504) au cours de laquelle le processeur (52) calcule, à partir des données ainsi recueillies, un écart entre une position actuelle de l'aéronef suiveur (20) et une position à atteindre pour l'aéronef suiveur (20), la position à atteindre correspondant à une position prédéterminée où l'aéronef suiveur (20) utilise au mieux l'effet de tourbillons (14a-b), - une étape de détermination (506) au cours de laquelle le système de navigation détermine des commandes à appliquer à des gouvernes de l'aéronef suiveur (20) pour annuler l'écart ainsi calculé, - une étape de commande (508) au cours de laquelle le système de navigation commande les gouvernes de l'aéronef suiveur (20) selon les commandes ainsi déterminées, et - un bouclage (510) sur l'étape de calcul (504). 2) Procédé de contrôle (500) selon la revendication 1, caractérisé en ce que l'étape de calcul (504) consiste en : - une première étape (602) au cours de laquelle le processeur (52) détermine des paramètres des tourbillons (14a-b) générés par l'aéronef meneur (10), - une deuxième étape (604) au cours de laquelle le processeur (52) détermine, à partir au moins des paramètres ainsi calculés, les positions des centres des tourbillons (14a-b) depuis leur génération par l'aéronef meneur (10), - une troisième étape (606) au cours de laquelle le processeur (52) détermine la position du centre de tourbillon le plus proche de la position actuelle de l'aéronef suiveur (20), - une quatrième étape (608) au cours de laquelle le processeur (52) détermine une position cible d'un centre du tourbillon le plus proche, ainsi qu'un décalage cible 3041121 26 par rapport à cette position cible, la position cible correspondant à une position d'un centre du tourbillon proche en amont de l'aéronef suiveur (20) pour un âge de tourbillon cible prédéterminé et le décalage cible étant un décalage entre l'aéronef suiveur (20) et le centre du tourbillon proche et en amont, 5 - une cinquième étape (609) au cours de laquelle le processeur (52) détermine la position à atteindre pour l'aéronef suiveur (20) par sommation vectorielle de la position cible du point (Gr) à l'âge cible sur les centres du tourbillon proche (14b), et le décalage cible, - une sixième étape (610) au cours de laquelle le processeur (52) calcule l'écart 10 correspondant à la différence entre la position à atteindre ainsi calculée et la position actuelle de l'aéronef suiveur (20). 3) Aéronef, dit aéronef suiveur (20), comportant un système de navigation et un processeur (52), le processeur (52) est configuré pour : - recueillir des données relatives à un aéronef, dit aéronef meneur (10), situé 15 devant l'aéronef suiveur (20) et générant des tourbillons (14a-b), ainsi que des données relatives à l'aéronef suiveur (20), ces données comprenant des données relatives à l'environnement atmosphérique, - calculer, à partir des données ainsi recueillies, un écart entre une position actuelle de l'aéronef suiveur (20) et une position à atteindre pour l'aéronef suiveur 20 (20), la position à atteindre correspondant à une position où l'aéronef suiveur (20) utilise au mieux l'effet des tourbillons (14a-b), - transmettre l'écart au système de navigation, et le système de navigation est configuré pour : - recevoir l'écart ainsi transmis, 25 - déterminer des commandes à appliquer à des gouvernes de l'aéronef suiveur (20) pour annuler l'écart, - commander les gouvernes de l'aéronef suiveur (20) selon les commandes ainsi déterminées. 4) Aéronef suiveur (20) selon la revendication 3, caractérisé en ce que le 30 processeur (52) est configuré pour : - déterminer des paramètres des tourbillons (14a-b) générés par l'aéronef meneur (10), 3041121 27 - déterminer, à partir au moins des paramètres ainsi calculés, les positions des centres des tourbillons (14a-b) depuis leur génération par l'aéronef meneur (10), - déterminer la position du centre de tourbillon le plus proche de la position actuelle de l'aéronef suiveur (20), 5 - déterminer une position cible d'un centre du tourbillon le plus proche, ainsi qu'un décalage cible par rapport à cette position cible, la position cible correspondant à une position d'un centre du tourbillon proche en amont de l'aéronef suiveur (20) pour un âge de tourbillon cible prédéterminé et le décalage cible étant un décalage entre l'aéronef suiveur (20) et le centre du tourbillon proche et en amont, 10 - déterminer la position à atteindre pour l'aéronef suiveur (20) par sommation vectorielle de la position cible du point (GT) à l'âge cible sur les centres du tourbillon proche (14b), et le décalage cible, - calculer l'écart correspondant à la différence entre la position à atteindre ainsi calculée et la position actuelle de l'aéronef suiveur (20). 15