L'objectif de ce travail de thèse est d'enrichir la formulation du Lancer de Faisceaux Gaussiens (LFG) et de tester sa capacité à répondre à certains des besoins actuels en calculs de propagation dans le domaine du Radar terrestre. Le LFG est envisagé comme une alternative possible aux méthodes classiques (Equation Parabolique, méthodes de rayons) en environnement complexe urbanisé, en particulier en présence d'obstacles latéraux, avec une cible située en non visibilité. La méthode de LFG "de base", qui utilise des expressions analytiques obtenues par approximation paraxiale, permet des calculs de propagation rapides en environnements complexes, sans problèmes de caustiques. Elle conduit à des résultats de précision satisfaisante dans le domaine millimétrique, par exemple pour des calculs de champs intra-bâtiments. Aux fréquences plus basses comme celles utilisées en Radar terrestre, elle est limitée par une prise en compte trop approximative des effets de diffraction et par l'élargissement spatial des faisceaux gaussiens au regard des dimensions des obstacles. La théorie des frames est utilisée dans cette thèse pour dépasser ces limites. La théorie des frames fournit un cadre rigoureux pour la décomposition initiale du champ rayonné en faisceaux gaussiens, et permet de calibrer le nombre et les directions des faisceaux à lancer. Dans ce travail de thèse, l'emploi de frames de fenêtres gaussiennes pour décomposer des distributions de champs ou de sources équivalentes est généralisé aux distributions de champs incidents sur des plans ou des portions de plans, choisis en fonction des obstacles rencontrés et des distances parcourues. Les champs rayonnés à partir de ces plans sont alors obtenus par sommation des faisceaux gaussiens lancés depuis ces frames dits de "re-décomposition". Les transformations de faisceaux gaussiens par des obstacles de taille limitée sont ainsi traitées par redécomposition : les faisceaux incidents partiellement interceptés par des surfaces limitées sont "re-décomposés" successivement sur deux frames de re-décomposition, à fenêtres "étroites" puis "larges", définis dans les plans de ces surfaces. Le frame à fenêtres "étroites" permet de traiter les discontinuités physiques, tandis que le frame à fenêtres "larges" permet de propager les champs transformés sous la forme de faisceaux "collimatés". Dans cette thèse, nous présentons une formulation de ces re-décompositions permettant une mise en œuvre numériquement efficace, grâce à des expressions analytiques approchées des coefficients de frame pour la première décomposition, et des éléments de la matrice de changement de frame pour la seconde. Cette formulation est mise en œuvre numériquement, et l'influence de différents paramètres sur la précision des re-décompositions est analysée. Finalement, l'algorithme de LFG enrichi de ces re-décompositions successives est utilisé dans un scénario simplifié proche de situations rencontrées en propagation Radar terrestre.

L’objectif de ce travail de thèse est d’enrichir la formulation du Lancer de Faisceaux Gaussiens (LFG) et de tester sa capacité à répondre à certains des besoins actuels en calculs de propagation dans le domaine du Radar terrestre. Le LFG est envisagé comme une alternative possible aux méthodes classiques (Equation Parabolique, méthodes de rayons) en environnement complexe urbanisé, en particulier en présence d’obstacles latéraux, avec une cible située en non visibilité. La méthode de LFG "de base", qui utilise des expressions analytiques obtenues par approximation paraxiale, permet des calculs de propagation rapides en environnements complexes, sans problèmes de caustiques. Elle conduit à des résultats de précision satisfaisante dans le domaine millimétrique, par exemple pour des calculs de champs intra-bâtiments. Aux fréquences plus basses comme celles utilisées en Radar terrestre, elle est limitée par une prise en compte trop approximative des effets de diffraction et par l’élargissement spatial des faisceaux gaussiens au regard des dimensions des obstacles. La théorie des frames est utilisée dans cette thèse pour dépasser ces limites. La théorie des frames fournit un cadre rigoureux pour la décomposition initiale du champ rayonné en faisceaux gaussiens, et permet de calibrer le nombre et les directions des faisceaux à lancer. Dans ce travail de thèse, l’emploi de frames de fenêtres gaussiennes pour décomposer des distributions de champs ou de sources équivalentes est généralisé aux distributions de champs incidents sur des plans ou des portions de plans, choisis en fonction des obstacles rencontrés et des distances parcourues. Les champs rayonnés à partir de ces plans sont alors obtenus par sommation des faisceaux gaussiens lancés depuis ces frames dits de "re-décomposition". Les transformations de faisceaux gaussiens par des obstacles de taille limitée sont ainsi traitées par redécomposition : les faisceaux incidents partiellement interceptés par des surfaces limitées sont "re-décomposés" successivement sur deux frames de re-décomposition, à fenêtres "étroites" puis "larges", définis dans les plans de ces surfaces. Le frame à fenêtres "étroites" permet de traiter les discontinuités physiques, tandis que le frame à fenêtres "larges" permet de propager les champs transformés sous la forme de faisceaux "collimatés". Dans cette thèse, nous présentons une formulation de ces re-décompositions permettant une mise en œuvre numériquement efficace, grâce à des expressions analytiques approchées des coefficients de frame pour la première décomposition, et des éléments de la matrice de changement de frame pour la seconde. Cette formulation est mise en œuvre numériquement, et l’influence de différents paramètres sur la précision des re-décompositions est analysée. Finalement, l’algorithme de LFG enrichi de ces re-décompositions successives est utilisé dans un scénario simplifié proche de situations rencontrées en propagation Radar terrestre. / In this work the Gaussian Beam Shooting (GBS) algorithm is complemented with new original formulations, and the ability of this "augmented" GBS algorithm to address specific problems encountered in electromagnetic field computations for ground-based Radar applications is tested. GBS is considered as an alternative to methods (Parabolic Equation, ray based methods) currently used for such computations in complex urban environments, especially when lateral obstacles and Non-Line-Of-Sight (NLOS) targets are involved. The "basic" GBS algorithm makes use of analytical expressions obtained through paraxial approximations. It allows to perform fast computations in complex environments, without suffering from any caustics problems. Reasonably accurate results have been obtained with this method in the millimetric range, e.g. for indoor field calculations. At lower frequencies, such as used in ground Radar systems, "basic" GBS cannot model diffraction effects accurately enough, and Gaussian beam width with respect to obstacle dimensions becomes a problem after some propagation distance. Frame theory is used in this PhD to overcome these limitations. Frame theory provides a rigorous framework for the initial decomposition of radiated fields into a set of Gaussian beams, providing flexible rules to adjust the number and directions of the launched beams. In this thesis, frame theory is used to discretize not only the source field distribution but also incident field distributions over planes or parts of planes of interest, according to encountered obstacles and propagation distances. The radiated fields are then obtained by summation of Gaussian beams launched from these frames called "reexpansion frames". Gaussian beam transformations by finite sized obstacles are addressed by this re-expansion scheme : the incident beams partially impinging on limited areas are successively "re-expanded" on two re-expansion frames, the first one composed of "narrow" windows and the second one of "wide" windows, both defined in the plane containing the limited area. Spatially narrow window frames allow to take into account abrupt transitions in space, and spatially wide window frames radiate in the form of collimated Gaussian beams. The re-expansion formulation proposed in this work is designed for efficient numerical implementation. Approximate analytical expressions are established for expansion coefficients on narrow window frames, and for frame change matrix elements. This formulation has been implemented, and the influence of frame parameters on re-expansion accuracy is analyzed. Finally, the GBS algorithm augmented with successive re-expansions is used to compute fields in simplified scenarios similar to situations encountered in ground-based Radar propagation problems.

Une méthode pour obtenir des faisceaux avec moment angulaire orbital consiste à modifier un faisceau gaussien avec une lame à vortex, appelée aussi lame de phase en spirale (SPP). La majorité de ces lames fonctionnent par transmission. Notre technique produit une lame de phase réflectrice, qui est adéquate pour les expériences avec un faisceau à haute puissance. Dans ce mémoire nous décrivons la fabrication de lames de phase en spirale continue par réflexion avec un système de déposition. Ces lames furent utilisées pour transformer un faisceau gaussien en un faisceau de Laguerre-Gauss avec un moment angulaire de différents ordres. Ces faisceaux sont caractérisés par leur front d'onde hélicoïdal et un zéro d'intensité au centre. Les résultats furent comparés à ceux obtenus avec un élément diffractif commercial. Avec l'aide d'un axicon, ces lames furent utilisées pour la formation de faisceaux Bessel avec un moment angulaire. En focalisant un faisceau Laguerre-Gauss femtoseconde avec un axicon, nous avons produit des modifications de surface sur un échantillon de verre en BK7.

En 1969, Lit et Tremblay [1] ont démontré qu'une série d'ouvertures circulaires disposées judicieusement sur un axe peut focaliser une onde sphérique. De telles structures agissent sur les distributions spatiale et spectrale des champs électromagnétiques et, par un simple facteur d'échelle, elles sont applicables à tout le spectre électromagnétique. Nous avons vu dans ces structures plusieurs applications potentielles, dont la focalisation de rayons X et le modelage d'impulsions laser. Afin d'investiguer la faisabilité éventuelle de ces applications, nous avons effectué des simulations numériques permettant d'étudier deux effets causés par une séquence d'ouvertures circulaires disposées axialement : la focalisation de faisceaux gaussiens et le modelage d'impulsions ultrabrèves. Les résultats de ces simulations ont montré qu'il est possible de focaliser un faisceau gaussien monochromatique divergent à l'aide de trois ouvertures circulaires coaxiales. Le point focal obtenu a une intensité allant jusqu'à 16 fois l'intensité du faisceau gaussien propagé dans le vide et il a une largeur à mi-hauteur transversale d'environ deux fois la longueur d'onde. Dans le cas des impulsions, un point focal a aussi été observé et, pour les impulsions d'une durée de 20 fs et moins dont le spectre est centré sur 800 nm, l'impulsion est dédoublée temporellement. Les structures étudiées peuvent donc être adaptées à la focalisation de rayons X et au modelage fréquentiel et temporel des impulsions brèves. Les calculs présentés dans ce mémoire peuvent aussi aider à prévenir des effets indésirables de focalisation de faisceaux ou d'impulsions brèves se propageant dans certaines nanostructures.

Tableau d'honneur de la Faculté des études supérieures et postdorales, 2015-2016

Depuis plusieurs années, des travaux en collaboration avec le Département Electromagnétisme et Radar de l'ONERA et le Groupe de Recherche en Electromagnétisme du LAPLACE ont permis l'élaboration d'un modèle basé sur les faisceaux gaussiens traitant l'interaction de champs électromagnétiques avec des structures complexes de grande taille. Dans ce formalisme asymptotique le champ incident est représenté comme une somme de faisceaux élémentaires. Il permet ainsi de ramener le calcul de l'interaction globale du champ avec la structure à une somme d'interactions plus simples. En outre, la distribution gaussienne des faisceaux et la forme canonique des portions d'objets interceptées permettent l'obtention de formulations analytiques des champs rayonnés. Toutefois, ce modèle nécessite la décomposition du champ incident sur une surface courbe présentant une amplitude nulle du champ sur ses bords. Cette contrainte n'est pas limitante dans de nombreux cas sauf celui, par exemple, d'une antenne posée sur un support. En effet, le champ rasant de part et d'autre de l'antenne peut alors être non nul et même interagir fortement avec le support en excitant des ondes de surface. Cette configuration fait apparaître au niveau du modèle un problème de fermeture du domaine de décomposition et de prise en compte des ondes de surface. Ces dernières revêtent un intérêt croissant lié d'une part à l'objectif de miniaturisation des systèmes micro-ondes et d'autre part, au développement récent des métamatériaux. Dans les deux cas les ondes de surface sont soit subies soit utilisées. Cependant, la définition même des ondes de surface et de leurs variantes telles que les ondes de fuite est peu connue et apparaît de manière confuse dans la littérature. Ainsi, l'objectif de cette thèse est double. Il s'agit d'une part, de faire une synthèse la plus claire possible sur les différents types d'ondes de surface et, d'autre part, de proposer une adaptation du modèle faisceaux gaussiens permettant leur prise en compte pour le calcul du couplage antenne-structure. Dans un premier temps, nous étudions les propriétés modales et l'excitation des ondes de surface et de fuite. Puis, nous adaptons des travaux de la littérature afin de modéliser, en deux dimensions, l'excitation de ces ondes par des courants de surface équivalents. A cette occasion nous proposons une formulation analytique, valable à grande distance, de l'interaction d'un faisceau gaussien avec une lame de métamatériau sur métal. Dans la deuxième partie du travail, nous proposons, en deux dimensions, une méthode d'hybridation entre la décomposition en faisceaux gaussiens et la décomposition en courant équivalents de surface. Cette méthode est appliquée au cas d'une antenne placée sur une lame de diélectrique sur métal Nous définissons alors les paramètres pertinents pour décrire l'hybridation des deux méthodes. Après une étude de leur influence sur la qualité des résultats, nous proposons un réglage par défaut donnant de bons résultats dans le cas général.

Depuis plusieurs années, des travaux ont été menés avec succès par le laboratoire LAME de l'UPS et le DEMR de l'ONERA sur le formalisme des faisceaux gaussiens pour modéliser l'interaction d'une onde électromagnétique avec un ensemble antenne-radôme. Ce formalisme est basé sur l'utilisation de faisceaux élémentaires gaussiens et a pour principal avantage l'obtention d'une expression analytique des champs électromagnétiques. Toutefois, ce modèle devait être complété par la description du champ électromagnétique rayonné par des discontinuités éclairées par un faisceau gaussien et par le traitement des surfaces de fortes courbures. \\ <br /><br /><br />Pour traiter le cas de la diffraction, l'utilisation de la méthode spectrale de la diffraction permet d'exprimer sous forme intégrale le champ diffracté par un objet canonique, par exemple un demi-plan conducteur, lorsqu'il est éclairé par un faisceau gaussien. L'approximation de l'optique physique permet d'obtenir une expression approchée analytique des champs rayonnés en 3D par des surfaces conductrices finies rectangulaires éclairées par un faisceau gaussien. Pour les surfaces de forte courbure et éclairées sous forte incidence, les faisceaux gaussiens conformes permettent de calculer analytiquement les champs lointains rayonnés. Afin de pouvoir traiter les interactions électromagnétiques (réflexions et transmissions) entre un faisceau gaussien conforme et une paroi diélectrique, le spectre d'ondes planes d'un faisceau gaussien conforme a été formulé. <br /><br /><br />Ces techniques ont été implémentées et validées sur différents cas tests canoniques. Les résultats théoriques ont été confrontés à des mesures en chambre anéchoïque ainsi qu'à plusieurs méthodes numériques. Les résultats obtenus montrent que les faisceaux gaussiens permettent de décrire une grande variété de problèmes électromagnétiques, des interactions antennes-radômes à la propagation sur de longues distances en présentant un bon compromis entre précision et temps de calcul.

Dans certains contextes, les méthodes classiques utilisées pour le calcul de champs rayonnés ou diffractés en présence d'obstacles de grande taille par rapport à la longueur d'onde, comme l'Optique Physique ou les méthodes de rayons, ne sont pas valides ou deviennent très lourdes en temps de calcul. La théorie des frames de Gabor fournit un cadre rigoureux permettant de décomposer une distribution de sources électromagnétiques, définie dans une ouverture équivalente plane, en une somme plus ou moins redondante de fenêtres gaussiennes. Cette décomposition peut servir de base à des algorithme de lancer de faisceaux gaussiens.Jusqu'à présent cette théorie était limitée à des décompositions dans un plan (rayonnement dans un demi-espace). L'objet de cette thèse est d'utiliser cette théorie pour décomposer des champs rayonnés ou diffractés dans toutes les directions de l'espace. Ce travail de thèse commence par une étude approfondie de l'influence des paramètres utilisés pour le calcul des coefficients de frame. La mise en oeuvre numérique permet de tester l'efficacité de techniques de troncation et de compression en termes de compromis précision/temps de calcul. Le coeur de la thèse consiste en une méthode originale de partitionnement spectral, utilisant des fonctions de partition de l'unité, qui permet d'utiliser le lancer de faisceaux gaussiens à partir de frames définis dans six plans, pour un rayonnement dans tout l'espace tridimensionnel. La formulation de la méthode est présentée. Elle est appliquée à la décomposition en faisceaux gaussiens du champ rayonné par des antennes théoriques omnidirectionnelles (réseau de dipôles et dipôle demi-onde). Une antenne réaliste sert enfin de cas test pour la mise en œuvre de la décomposition à partir de données expérimentales discrètes

Dans certains contextes, les méthodes classiques utilisées pour le calcul de champs rayonnés ou diffractés en présence d'obstacles de grande taille par rapport à la longueur d'onde, comme l'Optique Physique ou les méthodes de rayons, ne sont pas valides ou deviennent très lourdes en temps de calcul. La théorie des frames de Gabor fournit un cadre rigoureux permettant de décomposer une distribution de sources électromagnétiques, définie dans une ouverture équivalente plane, en une somme plus ou moins redondante de fenêtres gaussiennes. Cette décomposition peut servir de base à des algorithme de lancer de faisceaux gaussiens.Jusqu'à présent cette théorie était limitée à des décompositions dans un plan (rayonnement dans un demi-espace). L'objet de cette thèse est d'utiliser cette théorie pour décomposer des champs rayonnés ou diffractés dans toutes les directions de l'espace. Ce travail de thèse commence par une étude approfondie de l'influence des paramètres utilisés pour le calcul des coefficients de frame. La mise en oeuvre numérique permet de tester l'efficacité de techniques de troncation et de compression en termes de compromis précision/temps de calcul. Le coeur de la thèse consiste en une méthode originale de partitionnement spectral, utilisant des fonctions de partition de l'unité, qui permet d'utiliser le lancer de faisceaux gaussiens à partir de frames définis dans six plans, pour un rayonnement dans tout l'espace tridimensionnel. La formulation de la méthode est présentée. Elle est appliquée à la décomposition en faisceaux gaussiens du champ rayonné par des antennes théoriques omnidirectionnelles (réseau de dipôles et dipôle demi-onde). Une antenne réaliste sert enfin de cas test pour la mise en œuvre de la décomposition à partir de données expérimentales discrètes / In some contexts, conventional methods used for large problems involving radiated or diffracted field computations in the presence of obstacles, such as Physical Optics and ray based methods, become really inaccurate or prohibitively time-consuming. Gabor frame theory provides a rigorous framework for the initial decomposition of equivalent source distributions into a redundant set of Gaussian windows. Frame decomposition has been introduced as a first discretization step into Gaussian Beam Shooting (GBS) algorithms. Until now, frame decomposition has essentially been restricted to planar source distributions, radiating into one half space. The main goal of this thesis is to extend the application range of this theory to radiated or diffracted field decomposition into Gaussian beams propagating into the whole space. The thesis begins with a thorough study of influence of the parameters used for frame coefficient calculation. Numerical implementation is used to test the efficiency of truncation and compression techniques in terms of accuracy / computation time balance optimization. The core of the thesis consists of an original spectral domain partitioning method involving partition of unity functions, which allows to use Gaussian beam shooting from frames defined in six planes, for radiation into the whole three-dimensional space. The formulation of the method is presented and applied to the decomposition of fields radiated by theoretical omnidirectional antennas (dipole array and half-wave dipole) into Gaussian beams. A realistic antenna is used as a test case for the implementation of decompositions based on experimental discrete initial data

La thèse, dont le sujet est "Modélisation de la propagation électromagnétique en milieux à gradient d'indice basée sur les faisceaux gaussiens - Application à la propagation en atmosphère réaliste et à la radio-occultation entre satellites" a été commencée le 2 décembre 2013, au Département ÉlectroMagnétisme et Radar (DEMR) de l'Onera de Toulouse et avec le laboratoire LAPLACE de l'Université Paul Sabatier. Elle est co-financée par l'ONERA et par la Région Midi-Pyrénées. L'encadrement a été assuré par Jérôme Sokoloff (Laplace/UPS, directeur de thèse), Alexandre Chabory (ENAC, co-directeur) et Vincent Fabbro (ONERA). L'École Doctorale est l' "École Doctorale Génie Électrique, Électronique, Télécommunications : du système au nanosystème". Le faisceau gaussien a été principalement utilisé dans la recherche scientifique afin d'étudier les systèmes optiques tels que les lasers. Des études plus rares et plus récentes ont proposé de l'utiliser pour modéliser la propagation des ondes sismiques. Ses propriétés spatiales et spectrales ont amené certains auteurs à étudier son utilisation dans des modèles de propagation atmosphériques. Cette thèse a consisté à développer un modèle, appelé GBAR (Gaussian Beam for Atmospheric Refraction), de propagation troposphérique réaliste et déterministe en utilisant le formalisme des faisceaux gaussiens. La démarche adoptée a consisté à reprendre les équations fondamentales introduites par Cerveny et Popov décrivant de façon itérative la propagation d'un faisceau gaussien en milieu inhomogène, sous hypothèse de haute fréquence (modèle asymptotique). De nouvelles équations ont été développées à partir d'elles pour obtenir une description analytique de la propagation d'un faisceau gaussien dans un milieu troposphérique décrit par les variations spatiales de l'indice de réfraction. L'hypothèse de base pour l'obtention de la formulation analytique est que le gradient de l'indice de réfraction peut être considéré vertical et constant au voisinage du faisceau. Les équations analytiques pour la description de la propagation d'un seul faisceau ont ensuite été étendues à la modélisation d'un champ quelconque dans un milieu troposphérique pouvant contenir de fortes variations du gradient d'indice, y compris des inversions de gradient. Ceci a été réalisé en couplant les équations analytiques avec la procédure de décomposition multi-faisceaux développée dans sa thèse pas Alexandre Chabory. Le modèle GBAR a été validé dans des milieux troposphériques réalistes issus de simulations du modèle météo méso-échelle WRF (Weather Research and Forecasting). Dans un troisième temps, le modèle a été utilisé pour simuler des inversions de données de radio-occultation. Des outils existent pour fournir un modèle d'interprétation de ces données pour estimer les propriétés physiques de l'atmosphère à partir des mesures en phase, amplitude, Doppler et délai des signaux GNSS transmis entre satellites en orbite autour de la Terre / The subject of this PhD thesis is " Electromagnetic propagation modeling in inhomogeneous media with refractive index gradients based on Gaussian beams - Application to realistic atmospheric propagation and radio occultation between satellites ". The study started on december 2nd, 2013 at the DEMR (Département Électromagnétisme et Radar) department of the ONERA research laboratory, in Toulouse, France. It was funded both by the ONERA and Région Midi-Pyrénées. It was supervised by Jérôme Sokoloff (LAPLACE/UPS, thesis director), Alexandre Chabory (ENAC, thesis co-director) and Vincent Fabbro (ONERA). The doctoral school was "École Doctorale Génie Électrique, Électronique, Télécommunications : du système au nanosystème ". The Gaussian beam was mostly used in scientific investigations to study optical systems such as lasers. Rarer and more recent works suggested the use of the Gaussian beam formalism in order to model the propagation of seismic waves. The properties of the Gaussian beam also led some authors to develop models for atmospheric propagation. In this thesis a model based on Gaussian beams called GBAR (Gaussian Beam for Atmospheric Refraction) was developped for tropospheric propagation in realistic and deterministic conditions. The scientific approach consisted in rewritting the fundamental equations introduced by Cerveny and Popov describing iteratively the propagation of a Gaussian beam in inhomogeneous media, under the high-frequency assumption (asymptotic model). New equations were derived from them in order to get analytical equations of the propagation of a Gaussian beam in inhomogeneous media described by the variations of the refractive index. The basic assumption under to get the analytical equations is to consider that the refractive index gradient is vertical and constant around the beam axis. The analytical equations that describe the propagation of a Gaussian beam were extended to model the propagation of an arbitrary field in a tropospheric medium with strong variations and inversions of the refractive index. This was done by coupling the analytical equations with the multibeam expansion procedure developped by Alexandre Chabory in his PhD thesis. The GBAR model was validated in tropospheric conditions, using refractive index grids from the WRF (Weather Research and Forecasting) mesoscale meteorological model. In the third and final phase, the GBAR model was used to simulate Radio Occultation data inversions. Tools exist to allow for interpretations of Radio Occultation data in order to estimate the physical properties of the atmosphere from measured phased, amplitude, Doppler shift and delay of GNSS signals transmitted between satellites orbiting around the Earth