No description available.

Ce travail de thèse est une contribution à l'étude du rayonnement UVX cohérent produit par génération des harmoniques d'ordre élevé d'une impulsion laser intense focalisée dans un gaz rare. <br />Dans une première partie, nous présentons une étude approfondie de l'optimisation du flux harmonique produit soit dans une fibre creuse soit dans un jet. Dans chacun des cas, nous déterminons les facteurs limitant l'émission : accord de phase, absorption, ionisation, défocalisation. La comparaison des profils d'émission avec les simulations démontre le rôle crucial de la phase du dipôle atomique pour la génération dans une fibre. L'optimisation dans un jet conduit à des efficacités de conversion allant de 10-5 à 50nm (1010 ph/imp) à 10-7 à 15nm (5.107 ph/imp). <br />Dans une deuxième partie, nous étudions la focalisation du rayonnement harmonique par une lentille de Bragg-Fresnel, qui permet une focalisation efficace hors axe, sans aberration. Nous mesurons pour l'harmonique 39 une tache focale de l'ordre de 2µm. Compte tenu du flux élevé et de la courte durée (femtoseconde) des harmoniques, des éclairements élevés dans l'UVX peuvent être atteints.<br />Finalement, nous présentons les premières expériences de diagnostic de plasma par interférométrie UVX utilisant le rayonnement harmonique. Une technique originale est développée à partir de 2 sources harmoniques mutuellement cohérentes séparées spatialement. L'analogue temporel de l'interférométrie spatiale, qui utilise 2 impulsions harmoniques séparées en temps, est également démontré. Cette interférométrie fréquentielle a permis le premier diagnostic UVX de l'évolution temporelle d'un plasma avec une résolution femtoseconde. Une extension du schéma de 2 à 4 impulsions permet de mesurer avec une sensibilité extrême un déphasage équivalent à une demi-période harmonique, soit des écarts temporels à l'échelle attoseconde (1as=10-18s).

Lors de simulations numériques d'ordre élevé, la discrétisation sous-paramétrique du domaine de calcul peut générer des erreurs dominant l'erreur liée à la discrétisation des variables. De nombreux travaux proposent d'utiliser l'analyse isogéométrique afin de mieux représenter les géométries et de résoudre ce problème. Nous présenterons dans ce travail le couplage du schéma aux résidus distribués limité et stabilisé de Lax-Frieirichs avec l'analyse isogéométrique. En particulier, nous construirons une famille de fonctions de base permettant de représenter exactement les coniques et définies tant sur les éléments triangulaires que quadrangulaires : les fonctions de base de Bernstein rationnelles. Nous nous intéresserons ensuite à la génération de maillages précis pour l'analyse isogéométrique. Notre méthode consiste à créer un maillage courbe à partir d'un maillage linéaire par morceaux de la géométrie. Le maillage obtenu en sortie de notre procédure est non-structuré, conforme et assure la continuité de nos fonctions de base sur tout le domaine. Pour finir, nous décrirons les différentes méthodes d'adaptation de maillages développées : l'élévation d'ordre et le raffinement isotrope. Bien évidemment, la géométrie exacte du maillage courbe d'entrée est préservée au cours des processus d'adaptation.

Imager des objets non-périodiques à une échelle nanométrique et à une échelle femto seconde est un vrai challenge à notre époque. Les techniques d'imagerie " sans lentille " sont des moyens puissants pour répondre à ce besoin. En utilisant des sources ultrarapide (~fs) et cohérente (ex. laser à électron libre ou harmoniques d'ordres élevés), ces techniques nous permettent de reconstruire des objets à partir de leur figure de diffraction, remplaçant les optiques conventionnelles du système d'imagerie par un algorithme informatique. Dans ce travail de thèse, je présent des expériences d'imageries en utilisant un rayonnement extrême-UV (15~40 nm) produit par la génération d'harmoniques d'ordre élevé d'un laser infrarouge puissant. Ce manuscrit est constitué d'une introduction, un chapitre de background théorique, trois chapitres de travail de thèse et une conclusion générale avec perspectives. La première partie du travail de thèse porte sur les développements et caractérisations de la ligne de lumière avec l'objectif de générer maximum de photons harmoniques cohérents avec un front d'onde plat. La deuxième partie est consacrée aux expériences et analyses de trois techniques d'imageries " sans lentille " : Imagerie par diffraction cohérente (CDI), Holographie par la transformée de Fourier (FTH) et Holographie avec références étendues (HERALDO). Ces derniers nous permettent de reconstruire des objets avec une résolution spatiale de 78 nm dans le cas de CDI et de 112 nm dans le cas de HERALDO, tous les deux avec une résolution temporaire de 20 fs. La troisième partie est une première application physique de l'imagerie sur la ligne harmonique. Il s'agit des études statiques et dynamiques de nano-domaines magnétique avec une résolution spatiale sub-100 nm à l'échelle femto seconde. Perspective des techniques d'imagerie 3D et développement potentiel de la ligne d'harmoniques sont présentés à la fin.

Cette étude concerne le développement d'une méthode Galerkin discontinue d'ordre élevé en domaine temporel (DGTD), flexible et efficace, pour la résolution des équations de Maxwell instationnaires sur des maillages simplexes destructurés et reposant sur des schémas d'intégration en temps explicites. Les composantes du champ électromagnétique sont approximées localement par des méthodes d'interpolation polynomiale et la continuité entre éléments adjacents est renforcée de façon faible par un schéma centré pour le calcul du flux numérique à travers les interfaces du maillage. L'objectif de cette thèse est de remplir deux objectifs complémentaires. D'une part, améliorer la flexibilité de l'approximation polynomiale en vue du développement de méthodes DGTD p-adaptatives par l'étude de différentes méthodes d'interpolation polynomiale. Plusieurs aspects tels que la nature nodale ou modale de l'ensemble des fonctions de bases associées, leur éventuelle structure hiérarchique, le conditionnement des matrices élémentaires à inverser, les propriétés spectrales de l'interpolation ou la simplicité de programmation sont étudiés. D'autre part, augmenter l'efficacité de l'approximation temporelle sur des maillages localement raffinés en utilisant une stratégie de pas de temps local. Nous développerons finalement dans cette étude une méthodologie de calcul haute performance pour exploiter la localité et le parallélisme inhérents aux méthodes DGTD combinés aux capacités de calcul sur carte graphique. La combinaison de ces caractéristiques modernes résulte en une amélioration importante de l'efficacité et en une réduction significative du temps de calcul.

Cette thèse est scindée en deux parties. Dans un premier temps, nous présentons deux schémas volumes finis pour la discrétisation des problèmes de diffusion-convection-réaction sur les surfaces mouvantes. Le premier schéma présente une extension du schéma volumes finis avec flux à deux points sur les surfaces mouvantes. Le deuxième développe une méthode de type O-méthode. Cette dernière consiste à construire à partir des inconnus au centre des mailles, des solutions linéaires autour des nœuds de maillage qui intègrent la continuité des flux aux interfaces de mailles. La méthode permet aussi la construction des décentrages amont d'ordre 2 et ainsi, offre au procédé l'ordre 2 de convergence sur tout maillage non dégénéré. Ensuite, nous modélisons l'écoulement du couplage filme mince-surfactant (surface active agent) sur les surfaces mouvantes et simulons à l'aide des schémas volumes finis précédemment définis. Ici, l'utilisation du calcul tensoriel et de la méthode d'approximation par la lubrification permettent de réduire les équations de Navier-Stokes caractérisant le mouvement du filme mince en dimension 3 en un système d'équations définies sur la surface courbe mouvante dont l'inconnu est la hauteur du fluide. Le surfactant supposé insoluble est modélisé par une équation de diffusion convection à la l'interface fluid-air. Nous simulons l'ensemble en utilisant une méthode dite de capture d'interface (Interface tracking method) dérivée des volumes finis définis plus haut. Plusieurs exemples illustrent à suffisance l'efficacité et la précision des différentes méthodes.

La génération d'harmoniques d'ordre élevé en milieu gazeux est un phénomène habituellement décrit par un modèle à trois étapes : sous l'e et d'un champ laser intense, un atome (ou une molécule) est ionisé par e et tunnel. L'électron éjecté est accéléré dans le champ laser, puis il se recombine sur son ion parent en émettant un photon XUV. Ce rayonnement XUV, émis sous la forme d'impulsions attosecondes (1 as = 10^-18 s), est un outil idéal pour sonder la structure électronique des atomes ou des molécules, avec une résolution temporelle de l'ordre de l'attoseconde. Néanmoins, l'intensité de ce rayonnement n'est en général pas su sante pour induire des e ets non-linéaires (transitions à deux photons). Au cours des travaux réalisés pendant cette thèse, nous avons développé une source harmonique capable de produire un rayonnement XUV intense qui doit permettre d'accéder à la physique non-linéaire dans cette gamme de longueur d'onde. Pour parvenir à ces résultats, un travail important sur les impulsions infrarouges génératrices a été nécessaire, aussi bien dans le domaine spatial que dans le domaine temporel. Une technique de mise en forme spatiale de faisceaux laser intenses a donc été développée, ainsi qu'une technique de post compression adaptée aux impulsions laser intenses. Ce travail de thèse se divise donc en trois étapes : - Le développement de la source harmonique haute énergie et des diagnostics associés. Cette source est basée sur l'utilisation d'une chaîne laser Titane-Saphir qui délivre des impulsions de 150 mJ pour des durées de 40 fs à une cadence de 10 Hz. De bonnes conditions d'optimisation ont été obtenues, donnant lieu à des impulsions XUV dont l'énergie est de l'ordre du J lors de la génération dans l'argon. - Le développement d'une technique de mise en forme spatiale adaptée aux faisceaux laser intenses et à la génération d'harmoniques. Le dispositif est basé sur une optique en ré exion, et sur les interférences à deux faisceaux. Il permet de produire, dans la région focale, des faisceaux dont le pro l d'intensité est radialement constant (faisceaux at top) et ainsi d'apporter un contrôle supplémentaire sur la génération d'harmoniques d'ordre élevé. - Le développement d'une technique de post compression en propagation guidée basée sur l'élargissement spectral induit par ionisation. Cette technique est adaptée pour des impulsions intenses (3.5 TW) et permet de produire des impulsions de puissance crête supérieure au Térawatt dans le domaine sub-10 fs. Cette technique fournit donc une source unique pour la génération d'harmoniques d'ordre élevé. Ces deux approches ont été testées et validées pour la génération d'harmoniques d'ordre élevé, et les résultats obtenus ouvrent d'intéressantes perspectives telles que la génération d'impulsions attosecondes isolées de haute énergie (> 100 nJ).

Dans cette thèse, la dynamique de photodissociation de l'azoture de chlore (ClN3) est étudiée dans le domaine temporel par imagerie de vecteur vitesse des photofragments, spécialement du chlore et de N3. Cette imagerie résolue à l'échelle femtoseconde permet d'extraire les temps de dissociation, l'établissement temporel de la balance d'énergie de la réaction ainsi que la conservation des moments. Cette étude a permis de différencier deux domaines d'énergie: l'un menant à la formation d'un fragment N3 linéaire (étude autour de 4.5 eV d'excitation électronique) et le plus intéressant aboutissant à la formation d'un fragment N3 cyclique (autour de 6 eV). Dans une seconde étude, la dynamique de relaxation électronique du tétrathiafulvalène (C6H4S4-TTF) est étudiée autour de 4 eV par spectroscopie de masse résolue en temps ainsi que par spectroscopie de photoélectron. Les seuils d'ionisation dissociative sont extraits d'une détection en coïncidence entre les photoélectrons de seuil et les fragments ionisés réalisée sur rayonnement synchrotron. Les deux dernières expériences sont basées sur la génération d'harmoniques d'ordre élevé dans l'XUV d'une impulsion femtoseconde à 800 nm ou à 400 nm. Dans la première expérience, les harmoniques sont couplées à un imageur de vecteur vitesse en tant que rayonnement secondaire VUV. Par imagerie de photoélectron résolue en temps, nous avons révélé ainsi les dynamiques de relaxation des états de Rydberg initiée par une impulsion femtoseconde XUV à 15.5 eV dans l'argon et à 9.3 eV dans l'acétylène. Dans la seconde expérience, couramment nommée spectroscopie attoseconde, les harmoniques constituent le signal pompe sonde. Deux types de spectroscopie attoseconde ont été réalisés pour étudier la dynamique vibrationnelle de SF6: une expérience en réseau transitoire créé par deux impulsions pompe Raman avec une impulsion sonde intense générant les harmoniques à partir du réseau d'excitation et une expérience d'interférence de deux rayonnement XUV en champ lointain créés par deux impulsions sonde intenses.

Ces travaux de recherche ont pour but d’évaluer la méthode dite de la "Simulation auxEchelles Adaptées" (SAS pour Scale-Adaptive Simulation). Cette approche coïncide avec uneapproche RANS classique dans les zones pariétales attachées et adapte le niveau de viscositéturbulente dans les zones décollées pour y permettre une résolution partielle des structures turbulentes.Dans une première partie, une analyse théorique du modèle SAS original a été menéeet a permis de développer une correction visant à favoriser l’adaptation du niveau de viscositéturbulente dans les zones sièges d’instabilités de type Kelvin-Helmholtz. Le modèle ainsi corrigéest nommé SAS-αL. Les modèles SAS et SAS-αL ont été implantés dans le code de calculNavier-Stokes elsA de l’ONERA. À l’issue de cette étape, trois cas académiques d’écoulementsturbulents instationnaires, cylindre à grand nombre de Reynolds, marche descendante et cavitétranssonique, ont été simulés grâce aux trois modèles de turbulence SST, SAS et SAS-αL. Outreune comparaison aux bases de données expérimentales disponibles, une attention particulièrea été portée à l’influence de paramètres numériques tels que des schémas numériques d’ordreélevé. Enfin, afin d’étudier la viabilité de l’approche SAS dans un contexte industriel, les troismodèles de turbulence ont été testés sur une configuration issue de l’industrie aéronautique etcorrespondant à la sortie d’air chaud d’un système de dégivrage des nacelles d’avion. La comparaisondes prévisions obtenues avec les modèles SST, SAS et SAS-αL aux données expérimentalesobtenues à l’ONERA a permis de montrer un gain de précision grâce à l’emploi de l’approcheSAS et ce pour un coût de calcul compatible avec un cycle de conception industrielle. / This research work is meant to assess an upgraded URANS approach, namely the Scale-Adaptive Simulation (SAS). This method is similar to a conventional RANS approach (namelythe SSTmodel) in attached areas and is able to adapt the eddy-viscosity level in detached areas toensure the resolution, at least partially, of the turbulent structures. In a first part of this researchwork, an improvement of the SAS approach is suggestedto allowa better sensitivity of themodelto instabilities such as Kelvin-Helmholtz ones. This "improved" model is referred to as SAS-αLmodel. Both SAS and SAS-αL models were implemented in the ONERA Navier-Stokes solverelsA and both of themaswell as the SSTmodelwere tested on academic test cases : a cylinder in acrossflowat a high Reynolds number, a backward-facing step flowcorresponding to theDriver&Seegmiller experiment and the transonic flow over the M219 cavity experimentally investigatedby de Henshaw. The influence of the numerical parameters was deeply investigated and particularattention was paid to the high-order space-discretization schemes effects. The reliabilityof the SAS approach in an industrial framework was assessed on an aeronautic configurationnamely a nacelle de-icing device. Comparisons between the threemodels (SST, SAS and SAS-αL)and an experimental database available at ONERA - The French Aerospace Lab have shown thebetter accuracy of the SAS approach as well as the high potential of the SAS-αL model.

Dans ce travail de thèse, nous nous sommes intéressés à la modélisation des phénomènes de diffraction d’ondes électromagnétiques par des objets à symétrie de révolution complexes et fortement hétérogènes. La méthode que nous développons ici consiste en un couplage entre équations aux dérivées partielles (EDP) et équations intégrales (EI). Cette idée est essentiellement connue pour avoir deux avantages. Le premier est que les hétérogénéités de l’objet sont prises en compte naturellement dans la formulation du problème. Le deuxième est dû à l’utilisation des équations intégrales qui donnent une représentation exacte des solutions dans le milieu extérieur en fonction des courants surfaciques. Le domaine de simulation peut ainsi être ramené à l’objet lui-même. L’utilisation de développements en séries de Fourier combinés à la propriété d’invariance par rotation de l’objet permet alors la réduction du problème global 3D à un ensemble dénombrable de problème 2D.L’étude de ces problèmes nous a conduit à décomposer notre analyse en plusieurs parties,chacune ayant à traiter une partie du problème complet ou les méthodes d’intégrations numériques. Ces dernières étant difficiles à réaliser dans le cas des équations intégrales.Nous avons tout d’abord étudié un problème de Maxwell intérieur pour lequel nous avons développé une nouvelle méthode d’éléments finis d’ordre élevé dont nous avons montré l’efficacité et la précision sur de multiples exemples. Puis, nous avons étudié le problème de diffraction d’ondes planes pour des objets parfaitement conducteurs. La méthode d’éléments finis de frontière employée est alors construite par extension de la méthode précédente via l’opérateur de trace tangentielle. En combinant ces deux études, nous avons résolu le problème couplé en introduisant la propriété de symétrie de révolution dans une formulation variationnelle bien choisie. Par construction, les éléments finis qui y sont utilisés sont alors naturellement adaptées. L’algorithme de parallélisation de la méthode de couplage est finalement présentée et des comparaisons entre notre code AxiMax et un code 3D sont illustrées. Dans tous les cas, nous montrons que la méthode d’éléments finis d’ordre élevé permet d’obtenir des résultats d’une grande précision en fonction de la qualité des paramètres de simulation. / In this thesis, we are interested in modeling diffraction of electromagnetic waves by axisymmetric and highly heterogeneous objects. Our method consists in a coupling between partial differential equations and integral equations. This idea is mainly interesting for two reasons : heterogeneities are handled naturally in the formulation and integral equations give an analytical representation of solutions outside the object based on surface currents.These advantages allow us to limit the domain of simulation to the object itself. In addition,using Fourier series combined with the rotational invariance property of the object, the 3D problem is reduced to a countable set of 2D problems. The study of these problems is split into several parts. Each part has to deal with aspecific problem as for example the numerical integration of singular integrals which is difficult to achieve. As a first step, we study time-harmonic Maxwell’s equations in a bounded domain for which we develop a new high-order finite element method and present its efficiency and accuracy on many examples. Secondly, we consider the diffraction of plane waves by perfect electric conductors to analyse integral equations for these kind of object.The boundary finite element method applied is defined by extension of the previous one via tangential trace operator. Then, we solve the coupled problem using a well chosen formulation based on the previous studies for which our finite element method is naturally adapted by construction. In order to evaluate its efficiency, a comparison is performed between our program « AxiMax » and one based on a purely 3D model. To conclude, in the last two chapters, we present the numerical integration method and the multi-processing algorithm developed in AxiMax. In all cases, we put forward the fact that our finite element method provides accurate results depending on the quality of the simulation parameters.