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. / During high order simulations, the approximation error may be dominated by the errors linked to the sub-parametric discretization used for the geometry representation. Many works propose to use an isogeometric analysis approach to better represent the geometry and hence solve this problem. In this work, we will present the coupling between the limited stabilized Lax-Friedrichs residual distributed scheme and the isogeometric analysis. Especially, we will build a family of basis functions defined on both triangular and quadrangular elements and allowing the exact representation of conics : the rational Bernstein basis functions. We will then focus in how to generate accurate meshes for isogeometric analysis. Our idea is to create a curved mesh from a classical piecewise-linear mesh of the geometry. We obtain a conforming unstructured mesh which ensures the continuity of the basis functions over the entire mesh. Last, we will detail the curved mesh adaptation methods developed : the order elevation and the isotropic mesh refinement. Of course, the adaptation processes preserve the exact geometry of the initial curved mesh.

Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests. / The construction of compact high-order Residual Distribution schemes for the discretizationof steady multidimensional advection-diffusion problems on unstructuredgrids is presented. Linear and non-linear scheme are considered. A piecewise continuouspolynomial approximation of the solution is adopted and a gradient reconstructionprocedure is used in order to have a continuous representation of both thenumerical solution and its gradient. It is shown that the gradient must be reconstructedwith the same accuracy of the solution, otherwise the formal accuracy ofthe numerical scheme is lost in applications in which diffusive effects prevail overthe advective ones, and when advection and diffusion are equally important. Thenthe method is extended to systems of equations, with particular emphasis on theNavier-Stokes and RANS equations. The accuracy, efficiency, and robustness of theimplicit RD solver is demonstrated using a variety of challenging aerodynamic testproblems.

Cette thèse présente la construction d’un schéma aux Résidus Distribués p-adaptatif pour la discrétisation des équations d’Euler ainsi qu’un schéma aux Résidus Distribués hp-adaptatif pour les équations de Navier- Stokes pénalisées. On rappelle tout d’abord les équations d’Euler et de Navier-Stokes ainsi que leurs versions non dimensionnelles. Les définitions et propriétés de base des schémas aux Résidus Distribués sont ensuite présentées. On décrit alors la construction d’un schéma aux Résidus Distribués p-adaptatif pour les équations d’Euler. La construction du schéma p-adaptatif est basée sur la possibilité d’exprimer le résidu total d’un élément K de degré k (au sens où l’élément fini (K; P; Sigma ) est un élément fini de degré k) comme une somme pondérée des résidus totaux de ses sous-éléments de degré 1. La solution discrète ainsi obtenue est en général discontinue à l’interface entre un élément subdivisé et un élément non subdivisé. Ceci contredit l’hypothèse de continuité de la solution qui est utilisée pour démontrer le théorème de Lax-Wendroff discret pour les schémas aux Résidus Distribués. Cependant, on montre que cette hypothèse peut être assouplie. La conséquence pratique est que si l’on emploie des quadratures particulières dans l’implémentation numérique, on peut quand même démontrer le théorème de Lax-Wendroff discret, ce qui garantit la convergence du schéma numérique vers une solution faible des équations d’origine. Les formules qui permettent d’exprimer le résidu total comme une somme pondérée des résidus totaux des sous-éléments sont à la base de la méthode de p-adaptation présentée ici. Dans le cas quadratique, la formule est obtenue avec les classiques fonctions de base de Lagrange en dimension deux et avec des fonctions de base de Bézier en dimension trois. Ces deux formules sont ensuite généralisées à des degrés polynomiaux quelconques en dimension deux et trois avec des fonctions de base de Bézier. Dans la deuxième partie de la thèse, on présente l’application du schéma p-adaptatif aux équations pénalisées de Navier-Stokes avec adaptation de maillage anisotrope. . En pratique, on combine le schéma p-adaptatif avec la méthode IBM-LS-AUM (Immersed Boundary Method with Level Sets and Adapted Unstructured Meshes). La méthode IBM-LS-AUM permet d’imposer les conditions aux bords grâce à la méthode de pénalisation et l’adaptation anisotrope du maillage à la solution numérique et à la level-set augmente la précision de la solution et de la représentation de la surface. Une fois la méthode IBM-LS-AUM combinée avec le schéma p-adaptatif, il est alors possible d’utiliser des éléments d’ordre élevés en-dehors de la zone où la pénalisation est appliquée. La méthode est robuste comme le montrent les diverses expérimentations numériques à des vitesses faibles à élevées et à différents nombres de Reynolds. / This thesis presents the construction of a p-adaptive Residual Distribution scheme for the steady Euler equations and a hp-adaptive Residual Distribution scheme for the steady penalized Navier-Stokes equations in dimension two and three. The Euler and Navier-Stokes equations are recalled along with their non dimensional versions. The basis definitions and properties of the steady Residual Distribution schemes are presented. Then, the construction of a p-adaptive Residual Distribution scheme for the Euler equations is considered. The construction of the p-adaptive scheme is based upon the expression of the total residual of an element of a given degree k (in the Finite Element sense) into the total residuals of its linear sub-elements. The discrete solution obtained with the p-adaptive scheme is then a one degree polynomial in the divided elements and a k-th degree polynomial in the undivided ones. Therefore, the discrete solution is in general discontinuous at the interface between a divided element and an undivided one. This is in apparent contradiction with the continuity assumption used in general to demonstrate the discrete Lax-Wendroff theorem for Residual Distribution schemes. However, as we show in this work, this constrain can be relaxed. The consequence is that if special quadrature formulas are employed in the numerical implementation, the discrete Lax-Wendroff theorem can still be proved, which guaranties the convergence of the p-adaptive scheme to a weak solution of the governing equations. The formulas that express the total residual into the combination of the total residuals of the sub-elements are central to the method. In dimension two, the formula is obtained with the classical Lagrange basis in the quadratic case and with the Bézier basis in dimension three. These two formulas are then generalized to arbitrary polynomial degrees in dimension two and three with a Bézier basis. In the second part of the thesis the application of the p-adaptive scheme to the penalized Navier-Stokes equations with anisotropic mesh adaptation is presented. In practice, the p-adaptive scheme is used with the IBM-LS-AUM (Immersed Boundary Method with Level Sets and Adapted Unstructured Meshes) method. The IBM-LS-AUM allows to impose the boundary conditions with the penalization method and the mesh adaptation to the solution and to the level-set increases the accuracy of the representation of the surface and the solution around walls. When the IBM-LSAUM is combined with the p-adaptive scheme, it is possible to use high-order elements outside the zone where the penalization is applied. The method is robust as shown by the numerical applications at low to large Mach numbers and at different Reynolds in dimension two and three.

Le bruit de cavité est un phénomène très fréquent dans le domaine des transports aériens.Il survient notamment lors de l’approche à l’atterrissage, où des interactions entre la cellule de l’aéronef et l’écoulement sont à l’origine de fortes émissions tonales. Il devient dès lors une source de pollution acoustique non-négligeable pour les populations résidant à proximité de zones aéroportuaires. Les études numériques et expérimentales décrites jusqu’à présent dans la littérature abordent essentiellement le cas des cavités rectangulaires. Pourtant, les cavités rencontrées en pratique dans l’industrie aéronautique impliquent des géométries souvent plus complexes. Lorsque ces cavités sont soumises à une excitation de nature aérodynamique, leur spécificité géométrique conduit le plus souvent à des réponses acoustiques assez éloignées des estimations issues de modèles académiques construits sur l’observation de cavités rectangulaires. Quelques travaux seulement abordent le cas des cavités cylindriques.Ce travail est consacré à l’étude aéroacoustique des cavités cylindriques, à l’initiative d’Airbus. Il s’inscrit dans le cadre du projet AEROCAV soutenu par la Fondation de Recherche pour l’Aéronautique & l’espace (FRAE). Son objectif est de déterminer les mécanismes impliqués dans les émissions acoustiques intenses et tonales pour les configurations étudiées.Une première partie présente les résultats expérimentaux issus des campagnes de mesures menées dans la soufflerie anéchoïque du Centre Acoustique du LMFA et de l’école Centrale de Lyon. Un modèle semi-empirique, reposant sur l’hypothèse d’une résonance acoustique pilotée par les instabilités présentes dans la couche de cisaillement à l’ouverture de la cavité,est construit à partir du modèle d’Elder (1978). Le modèle permet d’estimer les fréquences susceptibles de dominer l’acoustique rayonnée en champ lointain à partir de la donnée du champ moyen de vitesse longitudinale, que l’on mesure dans le plan de l’écoulement par Vélocimétrie par Imagerie des Particules (PIV).Une seconde partie est destinée au calcul direct du bruit rayonné par un écoulement laminaire ou turbulent affleurant une cavité cylindrique de référence. Il consiste à calculer le champ acoustique directement à partir de la résolution des équations tridimensionnelles de la mécanique des fluides. Le solver Alesia est présenté dans une version modifiée et adaptée à la mise en oeuvre d’une approche multidomaine d’ordre élevé faisant intervenir plusieurs maillages se recouvrant. Des techniques d’interpolation sont spécifiquement développées en vue d’assurer une communication bidirectionnelle entre les différents maillages, malgré des contraintes géométriques fortes. Un modèle d’excitation de l’écoulement est aussi développé afin de disposer de fluctuations dans l’écoulement incident, pour le cas turbulent. Ces deux points font l’originalité des calculs réalisés.Les simulations, menées sur une cavité de rapport d’aspect géométrique égal à 1 et soumise à un écoulement incident à Mach 0.2, montrent que le rayonnement acoustique peut être fidèlement reproduit numériquement. La couche de cisaillement est caractérisée par la présence de deux larges structures tourbillonnaires s’amplifiant lors de leur convection. Leur présence s’accompagne de fortes fluctuations de vitesse à l’origine d’un débit aérodynamique de fluide à l’ouverture qui excite la cavité acoustiquement. Une résonance forcée s’établit dans celle-ci, excitant la couche de mélange au voisinage du point de séparation. Ce couplage auto-entretenu est à l’origine du rayonnement acoustique intense et fortement tonal de la cavité. Il s’établit à une fréquence proche de la fréquence prédite par le modèle semi-empirique développé. / Cavity noise is a very frequent phenomenon in air transport. It occurs in particular during landing approaches, where airframe–flow interactions are responsible for strong tonal emissions. Accordingly, it turns to be a non negligible source of acoustic pollution for populations living near airport areas. Numerical and experimental studies reported in the literature tackle essentially the case of rectangular cavities. Nevertheless, cavities may often exhibit more complex shapes in practice. When subject to aerodynamic excitations, and because of their geometrical specificity, these cavities may have acoustic responses which can be rather far from estimations resulting from academic models designed for rectangular cavities. Only asmall number of studies tackle the case of cylindrical cavities.The present work requested by Airbus is dedicated to the study of aeroacoustics in cylindrical cavities. This work was been supported by the Fondation de Recherche pour l’Aéronautique& l’Espace (FRAE) under contract reference AEROCAV. It aims at discerning the mechanisms responsible for strong and tonal acoustic emissions for the studied configurations.Experimental data resulting from measurements performed in the anechoic wind-tunnel of the Centre Acoustique at ´Ecole Centrale de Lyon are presented in a first part. A semi-empirical model based on the hypothesis of a shear-layer driven acoustic resonance is constructed from the Elder model (1978). The model provides an estimation of the frequences which are likely to be predominant in the far-field acoustics, given the mean streamwise velocity field, currently measured in the flow plane by Particle Image Velocimetry (PIV).A second part deals with the direct computation of the noise radiated by a laminar or turbulent grazing flow over a standard cylindrical cavity. The method consists in the calculationof the acoustic field directly from the resolution of the tridimensional Navier–Stokes equations. The Alesia solver is presented in a modified form, adapted to the implementationof a high-order chimera method involving several overlapping grids. Interpolation techniques have been specifically developed to achieve a bidirectional communication between the meshes in spite of strong geometrical constraints. A flow excitation model has also been constructed in order to obtain fluctuations into the incoming flow in the turbulent case. These two last points make the present computations original. The simulations, which are performed on a cavity of geometric ratio taken as 1 and subject to a grazing flow of Mach 0.2, reveal that it is possible to retrieve the radiated noise numerically with high fidelity. They indicate the presence of two large amplifying vortices in the shearlayer. These vortices go with strong velocity fluctuations giving rise to an inflow of fluid at the cavity mouth which excites the cavity acoustically. A forced acoustic resonance occurs into the cavity, then destabilises the shear layer near the separation point. This self-sustained coupling is responsible for strong tonal radiations from the cavity. The frequency of the radiated noise is close to the one predicted by the semi-empirical model.

La conception optimale de turbomachines repose usuellement sur des méthodes itératives avec des évaluations soit expérimentales, soit numériques qui peuvent conduire à des coûts élevés en raison des nombreuses manipulations ou de l’utilisation intensive de CPU. Afin de limiter ces coûts et de raccourcir les temps de développement, le présent travail propose d’intégrer une méthode de paramétrisation et de métamodélisation dans un cycle de conception d’une turbomachine axiale basse vitesse. La paramétrisation, réalisée par l’étude de sensibilité d’ordre élevé des équations de Navier-Stokes, permet de construire une base de données paramétrée qui contient non seulement les résultats d’évaluations, mais aussi les dérivées simples et les dérivées croisées des objectifs en fonction des paramètres. La plus grande quantité d’informations apportée par les dérivées est avantageusement utilisée lors de la construction de métamodèles, en particulier avec une méthode de Co-Krigeage employée pour coupler plusieurs bases de données. L’intérêt économique de la méthode par rapport à une méthode classique sans dérivée réside dans l’utilisation d’un nombre réduit de points d’évaluation. Lorsque ce nombre de points est véritablement faible, il peut arriver qu’une seule valeur de référence soit disponible pour une ou plusieurs dimensions, et nécessite une hypothèse de répartition d’erreur. Pour ces dimensions, le Co-Krigeage fonctionne comme une extrapolation de Taylor à partir d’un point et de ses dérivées. Cette approche a été expérimentée avec la construction d’un méta-modèle pour une hélice présentant un moyeu conique. La méthodologie fait appel à un couplage de bases de données issues de deux géométries et deux points de fonctionnement. La précision de la surface de réponse a permis de conduire une optimisation avec un algorithme génétique NSGA-2, et les deux optima sélectionnés répondent pour l’un à une maximisation du rendement, et pour l’autre à un élargissement de la plage de fonctionnement. Les résultats d’optimisation sont finalement validés par des simulations numériques supplémentaires. / The turbomachinery optimal design usually relies on some iterative methods with either experimental or numerical evaluations that can lead to high cost due to numerous manipulations and intensive usage of CPU. In order to limit the cost and shorten the development time, the present thesis work proposes to integrate a parameterization method and the meta-modelization method in an optimal design cycle of an axial low speed turbomachine. The parameterization, realized by the high order sensitivity study of Navier-Stokes equations, allows to construct a parameterized database that contains not only the evaluations results, but also the simple and cross derivatives of objectives as a function of parameters. Enriched information brought by the derivatives are utilized during the meta-model construction, particularly by the Co-Kriging method employed to couple several databases. Compared to classical methods that are without derivatives, the economic benefit of the proposed method lies in the use of less reference points. Provided the number of reference points is small, chances are a unique point presenting at one or several dimensions, which requires a hypothesis on the error distribution. For those dimensions, the Co-Kriging works like a Taylor extrapolation from the reference point making the most of its derivatives. This approach has been experimented on the construction of a meta-model for a conic hub fan. The methodology recalls the coupling of databases based on two fan geometries and two operating points. The precision of the meta-model allows to perform an optimization with help of NSGA-2, one of the optima selected reaches the maximum efficiency, and another covers a large operating range. The optimization results are eventually validated by further numerical simulations.

Le but de ce travail est d’étudier le potentiel d’un bitume naturel (BN) extrait en Albanie pour la production des bitumes durs (BD) et des enrobés à module élevé (EME). Pour la production des BD, différentes techniques de raffinage du pétrole existent. Néanmoins, avec certains bruts pétroliers il est impossible de les fabriquer. De plus, les BD de raffinerie comportent souvent des défauts qui limitent leurs applications. Ces raisons, couplées à des questions d’approvisionnement, conduisent à chercher des méthodes de substitution. Dans ce contexte, il est intéressant d’utiliser des BN. Cette étude a montré que le BN d’Albanie donne des BD et des EME en conformité avec la Normalisation Européenne, résistants au vieillissement et performants aux basses températures. Un dosage adéquat permet de formuler des matériaux d’une dureté désirée en réponse des exigences techniques des applications routières visées tout en facilitant la gestion des stocks de bitumes dans les centrales d’enrobages. / The aim of this thesis is to study the potential of a natural bitumen (NB) mined in Albania in the production of hard bitumens (HB) and that of high modulus asphalt concrete (HMAC). Various petroleum refining techniques are available for the production of HB. Nevertheless, this is impossible with some crude oils. In addition, HB from refineries often present shortcomings which limit their applications. These reasons, coupled with practical issues related to HB supply, motivate the research of alternative methods. In this industrial context it is interesting to use NB. This study has shown that the Albanian NB provides HB and HMAC in accordance with European Standards, resistant to aging and relatively performant at low temperatures. The proper dosage of this NB allows the formulation of materials of desired properties, in response to the technical requirements of the considered road applications, facilitating in this way the use of bitumen stocks in asphalt concrete production plants.

Les travaux de cette thèse concerne l’étude d’une méthode élément finis d’ordre spatial élevé que l’on peut assimilé à une extension du schéma de Yee. On parle alors de méthode différences finies d’ordre élevé. Après avoir donné, dans un premier chapitre, un historique non exhaustif des principales méthodes utilisées pour résoudre les équations de Maxwell dans le cadre de problèmes de CEM et montré l’ intérêt de disposer d’un solveur de type "différences finies d’ ordre élevé", nous présentons dans un deuxième chapitre le principe de la méthode. Nous donnons pour cela les caractéristiques du schéma spatial et temporel en précisant les conditions de stabilité de la méthode. En outre, dans une étude purement numérique, nous étudions la convergence du schéma. On se focalise ensuite sur la possibilité d’utiliser des ordres spatiaux variable par cellules dans chaque direction de l’espace. Des comparaisons avec le schéma de Yee et un schéma de Galerkin Discontinu particulier sont ensuite effectuées pour montrer les gains en coûts calcul et mémoire et donc l’intérêt de notre approche par rapport aux deux autres. Dans un troisième chapitre, nous nous intéressons à l’étude de modèles physiques indispensable au traitement d’un problème de CEM. Pour cela, nous nous focalisons particulièrement sur un modèle de fil mince oblique, des modèles de matériaux volumiques et minces et enfin sur la prise en compte de sol parfaitement métallique dans une agression de type onde plane. Chaque modèle est détaillé et validé par comparaison avec des solutions analytiques ou résultant de la littérature, sur des exemples canoniques. Le quatrième chapitre est dédié à une technique d’hybridation entre notre méthode et une approche Galerkin Discontinu en vue de traiter des géométries possédant des courbures. Nous donnons pour cela une stratégie d’hybridation basée sur l’échange de flux qui garantie au niveau continue la conservation d’une énergie. Nous présentons ensuite quelques exemples montrant la validité de notre approche dans une stratégie multi-domaines/multi-méthodes que nous précisons. Enfin le dernier chapitre de cette thèse concerne l’exploitation de notre méthode sur des cas industriels en comparaisons avec d’autres méthodes ou des résultats expérimentaux. / This thesis is about the study of a high spatial finite element method whichcan be assimilated at an extension of the Yee schema. In the next, this method is also called high order finite difference method. In the first chapter, we give a non exhaustive recall of the major methods used to treat EMC problems and we show the necessity to have this kind of schema to simulate efficiently some EMC configurations. In the second chapter, the principle of the numerical method is presented and a stability condition is given. A numerical study analysis of the schema convergence is also done. Next, we show the interest to have the possibility to use local spatial order by cell in each direction of the computational domain. Some canonic examples are given to show the advantages interms of CPU time and memory storage of the method by comparison with Yee’s scheme and DG approach. In the third chapter, we define and validate on several examples,some physical models as thin wire, materials and perfectly metallic ground in presence of a plane wave, to have the possibility to treat EMC problems. The fourth chapter is about a hybridization strategy between our high order FDTD method and a DG schema.We focalize our study on a hybrid method which provides an energy conservation of the continuous problem. A numerical example is given to validate the method. Finally, in the last chapter, we present some simulations on industrial problems to show the possibility of the method to treat realistic EMC problems.

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'effet d'un champ laser intense, un atome (ou une molécule) est ionisé par effet 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 suffisante pour induire des effets 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éflexion et sur les interférences à deux faisceaux. Il permet de produire, dans la région focale, des faisceaux dont le profil d'intensité est radialement constant (faisceaux flat 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). / The generation of high order harmonics in a gaseous medium is a phenomenon conveniently described by a three steps model : subject to a strong laser field irradiation, an atom (or molecule) can undergo a tunneling ionization. The ejected electron is accelerated in the laser field and recombine on its parent ion leading to the emission of an XUV photon. The XUV radiation can be emitted as attosecond pulses (1 as = 10-18 s), and it is then an ideal tool to probe the electronic structure of atoms or molecules which require the highest time resolution. However, the intensity of this radiation is usually not sufficient to induce non-linear processes (two-photon transitions).In the frame of this work, we have developed a harmonic source capable of producing an intense XUV radiation to access non-linear physics in this wavelength domain.To achieve these results, significant work on the infrared generating pulses was necessary, both in the spatial and temporal domain. We have developed a technique for spatial shaping of intense laser beams, and a post compression technique fitted to high energy pulses.This thesis is therefore divided into three parts:- The development of an high energy harmonic source and related diagnostics. We use a Ti: sapphire laser system for this source which delivers 40-fs pulses up to an energy of 150 mJ at 10 Hz repetition rate. Good optimization conditions were obtained, leading to XUV pulse energies of the order of μJ in the case of generation in argon.- The development of a spatial shaping technique adapted to intense laser beams and to harmonic generation. The device is based on reflection optics and the interferences of two beams. It can produce, in the focal region, beams with a radially constant intensity over a large volume (flat top beams) and thus provide additional control of the harmonics generating process.- The development of a post compression technique in guided geometry based on the ionization induced spectral broadening. This technique is suitable for intense pulses (3.5 TW) and produces pulses above the terawatt level in the 10-fs range. This technique therefore provides a unique source for harmonic generation.These two approaches have been tested and validated for high order harmonics generation, and the results open interesting perspectives such as the generation of isolated attosecond pulses of high energy (> 100 nJ).

Lors d’un lancement spatial, le bruit des jets supersoniques chauds, générés par les moteurs-fusées au décollage et en interaction avec le pas de tir, est dommageable pour le lanceur et en particulier sa charge utile. Par conséquent, les acteurs du spatial cherchent à renforcer leur compréhension et leur maîtrise de cette ambiance acoustique, entre autres grâce à des méthodes et outils numériques. Toutefois, ils ne disposent pas d’une approche numérique globale capable de prendre en compte simultanément la génération fidèle du bruit, la propagation acoustique non-linéaire, les effets d’installation complexes et les géométries réalistes, pourtant inhérents aux applications spatiales. Dans cette optique, cette étude consiste à mettre en place et valider une méthodologie de simulation numérique par couplage fort Navier-Stokes − Euler, puis à l’appliquer à des cas réalistes de bruit de jet supersonique. L’objectif est d’affiner les capacités de prévision et de contribuer à la compréhension des mécanismes de génération de bruit dans de tels jets. Le solveur Navier-Stokes repose sur une méthode LES sur maillage non-structuré et le solveur acoustique sur une méthode de Galerkine discontinue d’ordre élevé sur maillage non-structuré. La méthodologie est tout d’abord évaluée sur des cas académiques visant à valider la simulation par couplage fort. Après des calculs préliminaires, la méthodologie est appliquée à la simulation du bruit d’un jet libre supersonique à Mach 3.1. Une méthode de déclenchement géométrique de la turbulence est implémentée sous la forme d’une marche à la paroi de la tuyère. La simulation aboutit à des estimations du bruit très proches des mesures réalisées au banc MARTEL et met en évidence des effets non-linéaires significatifs ainsi qu’un mécanisme singulier de rayonnement des ondes de Mach. Dans une démarche de progression vers des cas toujours plus réalistes, l’ensemble de l’approche numérique est finalement adaptée avec succès à la simulation du bruit d’un jet en présence d’un carneau. À terme, elle pourra être étendue à des configurations multi-jets réactifs, avec injection d’eau, voire à l’échelle 1. / During a space launch, the noise from hot supersonic jets, generated by rocket engines at liftoff and interacting with the launch pad, is harmful to the launcher and in particular its payload. Consequently, space actors are seeking to strengthen their understanding and control of this acoustic environment through numerical methods and tools, among the others. However, they do not dispose of a comprehensive numerical strategy that can simultaneously take into account accurate noise generation, nonlinear acoustic propagation, complex installation effects and realistic geometries, which are inherent to space applications. For this purpose, the present study consists in setting up and validating a numerical simulation methodology using a Navier-Stokes − Euler two-way coupling approach, then applying it to realistic cases of supersonic jet noise in order to improve prediction capabilities and contribute to the understanding of the noise generation mechanisms in such jets. The Navier-Stokes solver is based on an LES method on unstructured mesh and the acoustic solver on a high-order discontinuous Galerkin method on unstructured mesh. The methodology is first assessed on academic cases to validate the use of the two-way coupling. After preliminary computations, the methodology is applied to the simulation of the noise from a supersonic free jet at Mach 3.1. A geometric turbulence tripping method is implemented via a step at the nozzle wall. The computation leads to noise predictions very close to the experimental measurements performed at the MARTEL test bench and highlights significant nonlinear effects as well as a quite particular Mach waves radiation mechanism. Targeting even more realistic cases, the entire numerical approach is finally successfully adapted to the simulation of the noise from a supersonic jet configuration including a flame trench. In the future, it may be extended to configurations with clustered reactive jets, water injection devices or even at full scale.

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. / Ultrafast imaging of isolated objects with nanometric spatial resolution is a great challenge in our time. The lensless imaging techniques have shown great potential to answer this challenge. In lensless imaging, one can reconstruct sample images from their diffraction patterns with computational algorithms, which replace the conventional lens systems. Using ultrafast and coherent light sources, such as free electron laser and high order harmonics, one can investigate dynamic phenomena at the femtosecond time scale. In this thesis work, I present the lenless imaging experiments using XUV radiation provided by a laser driven high order harmonic beamline. The manuscript is composed of an introduction, a chapter of theoretical background, three chapters of main research work and a general conclusion with perspectives. The first part of this work concerns the development of the harmonic beamline to optimize the illumination condition for lensless imaging. The second part concentrates on the imaging techniques: the Coherent Diffraction Imaging (CDI), the Fourier Transform Holography (FTH) and the Holography using extended references (HERALDO). The reconstructions have achieved 78 nm spatial resolution in case of CDI and 112 nm resolution in case of HERALDO, both in single-shot regime corresponding to a temporal resolution of 20 fs. The third part presents the first physical application on the harmonic beamline using the lensless imaging. Samples with magnetic nano-domains have been studied with sub-100 nm spatial resolution, which paves the way for ultrafast magnetic dynamic studies. At the end, single-shot 3D imaging and further beamline development have been discussed.