La prévention des risques industriels nécessite de simuler la dispersion turbulente de polluants. Cependant, les outils majoritairement utilisés à ce jour ne permettent pas de traiter les champs proches dans le cas de géométries complexes, et il est nécessaire d'utiliser les outils de CFD (“ Computational Fluid Dynamics ”) plus adaptés, mais plus coûteux. Afin de simuler les écoulements atmosphériques avec dispersion de polluants, les modèles CFD doivent modéliser correctement d'une part, les effets de flottabilité, et d'autre part les effets de la turbulence. Plusieurs approches existent, notamment dans la prise en compte des effets de flottabilité et la modélisation de la turbulence, et nécessitent des méthodes numériques adaptées aux spécificités mathématiques de chacune d'entre elles, ainsi que des schémas numériques précis pour ne pas polluer la modélisation. Une formulation d'ordre élevé en volumes finis, sur maillages non structurés, parallélisée, est proposée pour simuler les écoulements atmosphériques avec dispersion de polluants. L'utilisation de schémas d'ordre élevé doit permettre d'une part de réduire le nombre de cellules et diminuer les temps de simulation pour atteindre une précision donnée, et d'autre part de mieux contrôler la viscosité numérique des schémas en vue de simulations LES (Large Eddy Simulation), pour lesquelles la viscosité numérique des schémas peut masquer les effets de la modélisation. Deux schémas d'ordre élevé ont été étudiés et implémentés dans un solveur 3D Navier Stokes incompressible sur des maillages volumes finis non structurés. Nous avons développé un premier schéma d'ordre élevé, correspondant à un schéma Padé volumes finis, et nous avons étendu le schéma de reconstruction polynomiale de Carpentier (2000) aux écoulements incompressibles. Les propriétés numériques des différents schémas implémentés dans le même code de calcul sont étudiées sur différents cas tests bi-dimensionnels (calcul de flux convectifs et diffusifs sur une solution a-priori, convection d'une tâche gaussienne, décroissance d'un vortex de Taylor et cavité entraînée) et tri-dimensionnel (écoulement autour d'un obstacle cubique). Une attention particulière a été portée à l'étude de la précision et du traitement des conditions limites. L'implémentation proposée du schéma polynomial permet d'approcher, pour un maillage identique, les temps de simulation obtenus avec un schéma décentré classique d'ordre 2, mais avec une précision supérieure. Le schéma compact donne la meilleure précision. En utilisant une méthode de Jacobi sans calcul implicite de la matrice pour calculer le gradient, le temps de simulation devient intéressant uniquement lorsque la précision requise est importante. Une alternative est la résolution du système linéaire par une méthode multigrille algébrique. Cette méthode diminue considérablement le temps de calcul du gradient et le schéma Padé devient performant même pour des maillages grossiers. Enfin, pour réduire les temps de simulation, la parallélisation des schémas d'ordre élevé est réalisée par une décomposition en sous domaines. L'assemblage des flux s'effectue naturellement et différents solveurs proposés par les librairies PETSC et HYPRE (solveur multigrille algébrique et méthode de Krylov préconditionnée) permettent de résoudre les systèmes linéaires issus de notre problème. / The prevention of industrial risks requires simulating turbulent dispersion of pollutants. However, the tools mostly used so far do not allow near fields treated in the case of complex geometries, and it is necessary to utilize the tools of CFD (Computational Fluid Dynamics ") more suitable but more expensive. To simulate atmospheric flows with dispersion of pollutants, the CFD models must correctly model the one hand, the effects of buoyancy, and secondly the effects of turbulence. Several approaches exist, including taking into account the effects of buoyancy and turbulence modeling, and require numerical methods adapted to the specific mathematics of each, and accurate numerical schemes to avoid pollution modeling. A formulation of high order finite volume on unstructured meshes, parallelized, is proposed to simulate the atmospheric flows with dispersion of pollutants. The use of high order schemes allow one hand to reduce the number of cells and decrease the simulation time to achieve a given accuracy, and secondly to better control the viscosity numerical schemes for simulation LES (Large Eddy Simulation), for which the numerical viscosity patterns may mask the effects of modeling. Two high-order schemes have been studied and implemented in a 3D Navier Stokes solver on unstructured mesh finite volume. We developed the first high-order scheme, corresponding to a Padé finite volume scheme, and we have extended the scheme of reconstruction polynomial Carpentier (2000) for incompressible flows. The numerical properties of the various schemes implemented in the same computer code are studied different two-dimensional test cases (calculation of diffusive and convective flow on a solution a priori, a task Gaussian convection, decay of a vortex of Taylor and driven cavity) and tri-dimensional (flow past an obstacle cubic). Particular attention has been paid to the study of the accuracy and treatment of boundary conditions. The implementation of the polynomial allows to obtain quasi identical simulation time compared to a classical upwind scheme of order 2, but with higher accuracy. The compact layout gives the best accuracy. Using a Jacobi method without calculation implied matrix to calculate the gradient, the simulation time becomes interesting only when the required accuracy is important. An alternative is the resolution of linear system by an algebraic multigrid method. This method significantly reduces the computation time of the gradient and the Padé scheme is effective even for coarse meshes. Finally, to reduce simulation time, the parallelization schemes of high order is achieved by a decomposition into subdomains. The assembly flow occurs naturally and different solvers provided by PETSc libraries and HYORE (algebraic multigrid solver and preconditioned Krylov method) used to solve linear systems from our problem. The work was to identify and determine the parameters that lead to lowest time resolution simulation. Various tests of speed-up and scale-up were used to determine the most effective and optimal parameters for solving linear systems in parallel from our problem. The results of this work have been the subject of a communication in an international conference "Parallel CFD 2008" and an article submitted to "International Journal for Numerical Methods in Fluids" (Analysis of high-order finite volume schemes for the incompressible Navier Stokes equations)

These thesis review some recent results on the construction of very high order multidimensional upwind schemes for the solution of steady and unsteady conservation laws on unstructured triangular grids.<p>We also consider the extension to the approximation of solutions to conservation laws containing second order dissipative terms. To build this high order schemes we use a subtriangulation of the triangular Pk elements where we apply the distribution used for a P1 element.<p>This manuscript is divided in two parts. The first part is dedicated to the design of the high order schemes for scalar equations and focus more on the theoretical design of the schemes. The second part deals with the extension to system of equations, in particular we will compare the performances of 2nd, 3rd and 4th order schemes.<p><p>The first part is subdivided in four chapters:<p>The aim of the second chapter is to present the multidimensional upwind residual distributive schemes and to explain what was the status of their development at the beginning of this work.<p>The third chapter is dedicated to the first contribution: the design of 3rd and 4th order quasi non-oscillatory schemes.<p>The fourth chapter is composed of two parts: we start by understanding the non-uniformity of the accuracy of the 2nd order schemes for advection-diffusion problem. To solve this issue we use a Finite Element hybridisation.<p>This deep study of the 2nd order scheme is used as a basis to design a 3rd order scheme for advection-diffusion.<p>Finally, in the fifth chapter we extend the high order quasi non-oscillatory schemes to unsteady problems.<p>In the second part, we extend the schemes of the first part to systems of equations as follows:<p>The sixth chapter deals with the extension to steady systems of hyperbolic equations. In particular, we discuss how to solve some issues such as boundary conditions and the discretisation of curved geometries.<p>Then, we look at the performance of 2nd and 3rd order schemes on viscous flow.<p>Finally, we test the space-time schemes on several test cases. In particular, we will test the monotonicity of the space-time non-oscillatory schemes and we apply residual distributive schemes to acoustic problems. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

De nos jours, la simulation des écoulements diphasiques a de plus en plus d’importance dans les chambres de combustion aéronautiques en tant qu’un des éléments requis pour analyser et maîtriser le processus complet de combustion, afin d’améliorer la performance du moteur et de mieux prédire les émissions polluantes. Dans les applications industrielles, la modélisation du combustible liquide trouvé en aval de l’injecteur sous forme de brouillard de gouttes polydisperse, appelé spray, est de préférence faite à l’aide de méthodes Eulériennes. Ce choix s’explique par les avantages qu’offrent ces méthodes par rapport aux méthodes Lagrangiennes, notamment la convergence statistique intrinsèque, le couplage aisé avec la phase gazeuse ainsi que l’efficacité pour le calcul haute performance. Dans la présente thèse, on utilise une approche Eulérienne basée sur une fermeture au niveau cinétique de type distribution Gaussienne Anisotrope (AG). L’AG résout des moments de vitesse jusqu’au deuxième ordre et permet de capter les croisements des trajectoires (PTC) à petite échelle de manière statistique. Le système d’équations obtenu est hyperbolique, le problème est bien-posé et satisfait les conditions de réalisabilité. L’AG est comparé au modèle monocinétique (MK) d’ordre 1 en vitesse. Il est approprié pour la description des particules faiblement inertielles. Il mène à un système faiblement hyperbolique qui peut générer des singularités. Plusieurs schémas numériques, utilisés pour résoudre les systèmes hyperboliques et faible- ment hyperboliques, sont évalués. Ces schémas sont classifiés selon leur capacité à traiter les singularités naturellement présentes dans les modèles Eulériens, sans perdre l’ordre global de la méthode ni rompre les conditions de réalisabilité. L’AG est testé sur un champ turbulent 3D chargé de particules dans des simulations numériques directes. Le code ASPHODELE est utilisé pour la phase gazeuse et l’AG est implémenté dans le code MUSES3D pour le spray. Les résultats sont comparés aux de simulations Lagrangiennes de référence et aux modèle MK. L’AG est validé pour des gouttes modérément inertielles à travers des résultats qualitatifs et quantitatifs. Il s’avère prometteur pour les applications complexes comprenant des PTC à petite échelle. Finalement, l’AG est étendu à la simulation aux grandes échelles nécessaire dans les cas réels turbulents dans le domaine industriel en se basant sur un filtrage au niveau cinétique. Cette stratégie aide à garantir les conditions de réalisabilités. Des résultats préliminaires sont évalués en 2D pour tester la sensibilité des résultats LES sur les paramètres des modèles de fermetures de sous mailles. / In aeronautical combustion chambers, the ability to simulate two-phase flows gains increasing importance nowadays since it is one of the elements needed for the full understanding and prediction of the combustion process. This matter is motivated by the objective of improving the engine performance and better predicting the pollutant emissions. On the industrial scale, the description of the fuel spray found downstream of the injector is preferably done through Eulerian methods. This is due to the intrinsic statistical convergence of these methods, their natural coupling to the gas phase and their efficiency in terms of High Performance Computing compared to Lagrangian methods. In this thesis, the use of Kinetic-Based Moment Method with an Anisotropic Gaussian (AG) closure is investigated. By solving all velocity moments up to second order, this model reproduces statistically the main features of small scale Particles Trajectories Crossing (PTC). The resulting hyperbolic system of equations is mathematically well-posed and satisfies the realizability properties. This model is compared to the first order model in the KBMM hierarchy, the monokinetic model MK which is suitable of low inertia particles. The latter leads to a weakly hyperbolic system that can generate δ-shocks. Several schemes are compared for the resolution of the hyperbolic and weakly hyperbolic system of equations. These methods are assessed based on their ability to handle the naturally en- countered singularities due to the moment closures, especially without globally degenerating to lower order or violating the realizability constraints. The AG is evaluated for the Direct Numerical Simulation of 3D turbulent particle-laden flows by using ASPHODELE solver for the gas phase, and MUSES3D solver for the Eulerian spray in which the new model is implemented. The results are compared to the reference Lagrangian simulation as well as the MK results. Through the qualitative and quantitative results, the AG is found to be a predictive method for the description of moderately inertial particles and is a good candidate for complex simulations in realistic configurations where small scale PTC occurs. Finally, within the framework of industrial turbulence simulations a fully kinetic Large Eddy Simulation formalism is derived based on the AG model. This strategy of directly applying the filter on the kinetic level is helpful to devise realizability conditions. Preliminary results for the AG-LES model are evaluated in 2D, in order to investigate the sensitivity of the LES result on the subgrid closures.

L'objectif principal de ce travail est d'améliorer la précision et l'efficacité des modèles LES au moyen des méthodes Galerkine discontinues (DG). Deux thématiques principales ont été étudiées: les stratégies d'adaptation spatiale et les modèles LES pour les méthodes d'ordre élevé.Concernant le premier thème, dans le cadre des méthodes DG la résolution spatiale peut être efficacement adaptée en modifiant localement soit le maillage (adaptation-h) soit le degré polynômial de la solution (adaptation-p). L'adaptation automatique de la résolution nécessite l'estimation des erreurs pour analyser la qualité de la solution locale et les exigences de résolution. L'efficacité de différentes stratégies de la littérature est comparée en effectuant des simulations h- et p-adaptatives.Sur la base de cette étude comparative, des algorithmes statiques et dynamiques p-adaptatifs pour la simulation des écoulements instationnaires sont ensuite développés et analysés. Les simulations numériques réalisées montrent que les algorithmes proposés peuvent réduire le coût de calcul des simulations des écoulements transitoires et statistiquement stationnaires.Un nouvel estimateur d'erreur est ensuite proposé. Il est local, car n'exige que des informations de l'élément et de ses voisins directs, et peut être calculé en cours de simulation pour un coût limité. Il est démontré que l'algorithme statique p-adaptatif basé sur cet estimateur d'erreur peut être utilisé pour améliorer la précision des simulations LES sur des écoulements turbulents statistiquement stationnaires.Concernant le second thème, une nouvelle méthode, consistante avec la discrétisation DG, est développée pour l'analyse a-priori des modèles DG-LES à partir des données DNS. Elle permet d'identifier le transfert d'énergie idéal entre les échelles résolues et non résolues. Cette méthode est appliquée à l'analyse de l'approche VMS (Variational Multiscale). Il est démontré que pour les résolutions fines, l'approche DG-VMS est capable de reproduire le transfert d'énergie idéal. Cependant, pour les résolutions grossières, typique de la LES à nombres de Reynolds élevés, un meilleur accord peut être obtenu en utilisant un modèle mixte Smagorinsky-VMS. / The main goal of this work is to improve the accuracy and computational efficiency of Large Eddy Simulations (LES) by means of discontinuous Galerkin (DG) methods. To this end, two main research topics have been investigated: resolution adaptation strategies and LES models for high-order methods.As regards the first topic, in the framework of DG methods the spatial resolution can be efficiently adapted by modifying either the local mesh size (h-adaptation) or the degree of the polynomial representation of the solution (p-adaptation).The automatic resolution adaptation requires the definition of an error estimation strategy to analyse the local solution quality and resolution requirements.The efficiency of several strategies derived from the literature are compared by performing p- and h-adaptive simulations. Based on this comparative study a suitable error indicator for the adaptive scale-resolving simulations is selected.Both static and dynamic p-adaptive algorithms for the simulation of unsteady flows are then developed and analysed. It is demonstrated by numerical simulations that the proposed algorithms can provide a reduction of the computational cost for the simulation of both transient and statistically steady flows.A novel error estimation strategy is then introduced. It is local, requiring only information from the element and direct neighbours, and can be computed at run-time with limited overhead. It is shown that the static p-adaptive algorithm based on this error estimator can be employed to improve the accuracy for LES of statistically steady turbulent flows.As regards the second topic, a novel framework consistent with the DG discretization is developed for the a-priori analysis of DG-LES models from DNS databases. It allows to identify the ideal energy transfer mechanism between resolved and unresolved scales.This approach is applied for the analysis of the DG Variational Multiscale (VMS) approach. It is shown that, for fine resolutions, the DG-VMS approach is able to replicate the ideal energy transfer mechanism.However, for coarse resolutions, typical of LES at high Reynolds numbers, a more accurate agreement is obtained by a mixed Smagorinsky-VMS model.

Cette thèse porte sur différents aspects des dynamiques ultra-rapides d’atomes et de molécules soumises à des impulsions laser infrarouges courtes et intenses. Nous étudions des processus fortement non linéaires tels que l’ionisation tunnel, la génération d’harmoniques d’ordre élevé ou l’ionisation au-dessus du seuil. Deux approches différentes sont utilisées. D’un côté nous mettons au point des modèles analytiques approchés qui nous permettent de construire des interprétations physiques de ces processus. D’autre part nous appuyons les interprétations données par ces modèles avec les résultats obtenus par des simulations numériques qui résolvent explicitement l’équation de Schrödinger dépendante du temps en dimension réduite. Nous étudions également une méthode numérique basée sur l’interaction de configuration dépendante du temps afin de pouvoir des décrire des systèmes à plusieurs électrons plus gros et plus complexes. / In this thesis we study different aspects of the ultrafast dynamics of atoms and molecules triggered by intense and short infrared laser pulses. Highly non-linear processes like tunnel ionization, high order harmonic generation and above threshold ionization are investigated. Two different and complementary approaches are used. On the one hand we construct approximate analytical models to get physical insight on these processes. On the other hand, these models are supported by the results of accurate numerical simulations that explicitly solve the time dependent Schrödinger equation for simple benchmark models in reduced dimensions. A numerical method based on time dependent configuration interaction is investigated to describe larger and more more complex systems with several electrons

La génération d’harmoniques d’ordre élevé est un processus d’interaction lumière-matière hautement non-linéaire permettant la synthèse d’impulsions sub-femtosecondes, dites attosecondes (1 as = 10⁻¹⁸ s). Mes travaux de thèse portent sur l’étude du transfert de moment angulaire lors de ce processus, afin de contrôler les caractéristiques spatiales et de polarisation du rayonnement émis dans l’extrême ultraviolet. Comme pour la matière, le moment angulaire de la lumière peut être séparé en une composante de spin, associée à l’état de polarisation du faisceau, et une composante orbitale, reliée à la forme du front d’onde. La maitrise complète du moment angulaire des harmoniques nécessite de recourir à des schémas de génération à deux faisceaux non-colinéaires, créant un réseau de diffraction dans le milieu générateur. Nous avons montré que, bien que les règles de transfert obéissent à des lois de conservation du moment angulaire, la description fine du phénomène requiert une analyse précise du champ laser dans le milieu de génération. Ces travaux ouvrent des perspectives de mise en forme avancée des impulsions attosecondes. / High-order harmonic generation is a highly nonlinear laser-matter interaction process which allows the synthesis of sub-femtosecond pulses, also called attosecond (1 as = 10⁻¹⁸ s) pulses. My PhD is centered around the study of angular momentum transfer during this process, in order to control spatial and polarization features of the radiation which is emitted in the extreme ultraviolet. As for matter, the angular momentum of light can be divided into a spin component, associated with the beam’s polarization, and an orbital component, related to the shape of the wavefront. The control of high harmonics’ angular momentum requires generating schemes involving two crossing beams, thus creating a diffraction grating in the generating medium.We have shown that, although the transfer rules obey conservation laws of the angular momentum, the fine description of the phenomenon requires an accurate analysis of the laser field in the generation medium. This work opens the road for advanced shaping of attosecond pulses.

Les impulsions ultra-brèves de rayonnement ultraviolet extrême (UVX) ont un grand champ d’application dans les domaines tels que le diagnostic de plasmas, la spectroscopie ou l’étude de la dynamique ultrarapide dans les atomes et les molécules.Aujourd’hui, il existe trois sources délivrant ce genre d’impulsions. Les harmoniques d’ordre élevé (HHG, en anglais) générés dans les gaz rares ou sur les solides peuvent fournir des impulsions attosecondes. Cependant, leur énergie, le plus souvent de l’ordre du nanojoule, limite les applications. L’amplification des impulsions harmoniques dans les plasmas créés par laser (SXRL, en anglais) a démontré pouvoir fournir des énergies de plusieurs dizaines de microjoules. Des énergies plus élevées peuvent être obtenues avec les lasers à électrons libres (LEL) UVX injectés, mais ce sont des Très Grandes Infrastructures ayant un accès limité.Ces dernières années, des progrès significatifs ont été réalisé avec chacune des ces sources, avec pour objectif la génération d’impulsions plus brèves. Il est devenu nécessaire de développer des nouvelles techniques de métrologie temporelle des impulsions UVX ultra-brèves. De plus, beaucoup d’expériences, comme ceux impliquant des phénomènes non-linéaires, nécessitent de hautes intensités UVX. La focalisation efficace des impulsions de faibles énergies peut significativement augmenter le domaine d’application. De bons fronts d’onde sont nécessaires pour focaliser les impulsions UVX à haute intensité, et les optiques doivent aussi être de bonne qualité et alignées avec précision.Dans cette thèse, les propriétés spatiales des harmoniques d’ordre élevé ont été extensivement étudiées grâce à un senseur de front d’onde UVX. Cet appareil couplé à une source HHG a démontré être utile pour la caractérisation de table et à la longueur d’onde ainsi que pour l’optimisation de systèmes optiques UVX.Le problème de la mise en place de la complète caractérisation temporelle d’impulsions UVX est aussi discuté en détail, et deux nouveaux schémas pour la reconstruction d’impulsions de LEL injectés et de lasers X à plasma sont présentés. Finalement, la première implantation d’un système d’amplification à dérive de fréquence (CPA, en anglais) sur un LEL UVX est présentée et son implantation pour les lasers X à plasmas est aussi discutée. / Ultrashort pulses of extreme-ultraviolet (XUV) radiation have a wide range of applications in fields such as plasma probing, spectroscopy, or the study of ultrafast dynamics in atoms and molecules.Nowadays, there are three main sources of such pulses. High-order harmonic generation (HHG) in rare gases or solid surfaces is able to provide attosecond pulses. However, their limited energy, of the order of nanojoules, limits its number of applications. The amplification of high-harmonic pulses in laser-driven plasmas (SXRL) has been demonstrated to provide energies of tens of microjules. Higher pulse energies can be obtained from seeded XUV free-electron lasers (FELs), large-scale facilities with more limited accessibility.In recent years, significant progress has been made with each of these sources towards the generation of shorter pulses. It is thus necessary to develop new techniques for full temporal metrology of ultrashort XUV pulses. Additionally, many experiments, such as those involving nonlinear phenomena, require high XUV intensities. Efficient focusing of low-energy pulses can significantly increase their range of application. Good wavefronts are required in order to focus XUV pulses to high intensities, and the optics must be of high quality and precisely aligned.In this thesis, the spatial properties of high-harmonic pulses are extensively explored thanks to the use of an XUV Hartmann wavefront sensor. This device is also proven here to be useful for tabletop, at-wavelength characterization and optimization of XUV optical systems with HHG sources.The problem of performing full temporal characterization of XUV pulses is also discussed in detail, and two new schemes for complete pulse reconstruction for seeded XUV FELs and seeded SXRLs are presented. Finally, the first implementation of chirped pulse amplification (CPA) in a seeded XUV FEL is reported, and its implementation in seeded SXRLs is discussed as well.

Depuis sa première observation, la génération d'harmoniques d'ordre élevé (GHOE) dans les gaz a demontré son importance, ouvrant la voie à la science attoseconde. Cette technique produit un rayonnement impulsionnel XUV qui s'étend dans le domaine spectral intermédiaire entre l'ultraviolet et les rayons X. Ces impulsions attosecondes donnent accès à des résolutions temporelles extrêmes, permettant ainsi d'observer des dynamiques électroniques dans des atomes ou des molécules. En effet le processus de généneration d'harmonique repose sur l'oscillation de paquets d'électrons attosecondes issus des molécules, accélérés par le champ de laser intense et se recombinant radiativement avec leurs ions moléculaires parents. Ainsi, le rayonnement harmonique émis lors de la recombinaison permet d'encoder l'information structurale sur le ou les orbitales impliquées avec une résolution spatiale de l'ordre l'Angström et temporelle femtoseconde ou attoseconde. La génération d'harmonique peut être utilisée comme signal de sonde dans des expériences de spectroscopie pompe-sonde résolue en temps. Ces expériences de spectroscopie harmoniques permettent d'étudier la structure des orbitales et les dynamiques moléculaires ultra-rapides. L'objectif de cette thèse est d'utiliser le processus de la GHOE, pour sonder les processus fondamentaux qui interviennent dans les atomes, les molécules et la matière condensée. Tout d'abord, pour comprendre comment extraire des informations dynamiques ou structurelles sur les orbitales à partir du signal harmonique nous avons étudié un système simple et connu: l'argon. Une nouvelle approche théorique développée par Fabre et Pons a permis de reproduire fidèlement l'expérience. Nous avons continué à étudier la structure et la dynamique moléculaire dans N2 et CO2. Les molécules issues d'un jet supersonique Even-Lavie qui permettait d'obtenir des températures rotationelles de moins de 10K ont été alignées par laser avec un fort degré d'alignement. Ce type de jet permet d'améliorer la sensibilité à la structure des orbitales impliquées et d'identifier la contribution de plusieurs orbitales. Ensuite nous avons utilisé la sensibilité de la génération des harmoniques d'ordre élevé à la structure des orbitales moléculaires pour sonder la dynamique complexe du NO2 excité autour d'une intersection conique. Nous avons appliqué la méthode du réseau d'excitation transitoire qui permet d'améliorer la sensibilité aux molécules excitées. Nous avons donc mené une étude dans les agrégats. A l'aide d'une étude différentielle en température et d'une méthode de cartographie spectrale et spatiale, nous avons pu isoler la contibution des grands agrégats. Notre analyse suggère un nouveau mécanisme de génération par des agrégats et permet même une estimation de la longeur de corrélation des électrons dans les agrégats. Ce manuscrit se termine avec la présentation d'une ligne de lumière XUV. Cette technique consiste à utiliser le rayonnement XUV fs produit par la GHOE comme impulsion sonde pour ioniser des fragments de dissociation moléculaire à l'aide d'une transition à un photon. / Since the first observation of high-order harmonic spectra in gases, high harmonic generation (HHG) has demonstrated its importance, opening a door to the field of attosecond sience. The bandwidth of the emitted spectrum reaches up to the XUV. The attosecond pules reach a very high time resolution, allowing the study of electron dynamics in atoms or molecules. The generation mechanism of HHG is based on the oscillation of the attosecond electron wavepacket emitted by the atoms/molecules, accelerated by the laser field. The electron wavepacket finally recombines radiatively with its parent ion. Thus the structural information of the probed orbital is encoded in the high harmonic spectrum with a spatial resolution of one Angtröm and a temporal resolution of few femtoseconds. HHG can be used as a probe signal resolved for pump-probe spectroscopy. High harmonic spectroscopy allows the study of the orbital structure and ultra-fast molecular dynamics.In this thesis the fundamental mechanisms playing a role in atoms, molecules and condensed matter are probed using HHG. In order to understand how to extract dynamical and structural information of orbitals from a harmonic signal, we have studied an easy and well known systems: the argon atom. A new theoretical approach developped by Fabre and Pons allowed us to reproduce the experimental results in good agreement. We continued with a study of the molecular structure and dynamics of N2 and CO2. A supersonic Even-Lavie jet permitted to reach rotational temperatures lower than 10K with an excellent alignment distribution. Owing to the good alignment in such gas jet, we were able to resolve the orbital structure with a higher sensitivity and to identify the contribution of several orbitals. In the next step we used the sensitivity of HHG towards the structure of molecular orbitals in order to probe the complex dynamics of NO2 in the vicinity of a conical intersection. We applied HHG combined with transient grating spectroscopy which leads to a higher sensitivity of the excited molecules. We then continued with studying cluster. We were able to disentangle the contribution of large clusters to the harmonic signal due to a 2D spatio-spectral representation of a temperature dependent differential measurement. Our analysis suggests a new generation mechanism in clusters and allows an estimation of the electron correlation length in clusters. This thesis ends with the presentation of a XUV beamline. This technique uses the emitted fs-XUV radiation, provided by HHG, as a probe pulse for ionizing the photofragments by a one photon transition.

Cette thèse est consacrée à l’étude théorique et numérique de plusieurs équations aux dérivées partielles non linéaires qui apparaissent dans la modélisation de la séparation de phase et des micro-systèmes électro-mécaniques (MSEM). Dans la première partie, nous étudions des modèles d’ordre élevé en séparation de phase pour lesquels nous obtenons le caractère bien posé et la dissipativité, ainsi que l’existence de l’attracteur global et, dans certains cas, des simulations numériques. De manière plus précise, nous considérons dans cette première partie des modèles de type Allen-Cahn et Cahn-Hilliard d’ordre élevé avec un potentiel régulier et des modèles de type Allen-Cahn d’ordre élevé avec un potentiel logarithmique. En outre, nous étudions des modèles anisotropes d’ordre élevé et des généralisations d’ordre élevé de l’équation de Cahn-Hilliard avec des applications en biologie, traitement d’images, etc. Nous étudions également la relaxation hyperbolique d’équations de Cahn-Hilliard anisotropes d’ordre élevé. Dans la seconde partie, nous proposons des schémas semi-discrets semi-implicites et implicites et totalement discrétisés afin de résoudre l’équation aux dérivées partielles non linéaire décrivant à la fois les effets élastiques et électrostatiques de condensateurs MSEM. Nous faisons une analyse théorique de ces schémas et de la convergence sous certaines conditions. De plus, plusieurs simulations numériques illustrent et appuient les résultats théoriques. / This thesis is devoted to the theoretical and numerical study of several nonlinear partial differential equations, which occur in the mathematical modeling of phase separation and micro-electromechanical system (MEMS). In the first part, we study higher-order phase separation models for which we obtain well-posedness and dissipativity results, together with the existence of global attractors and, in certain cases, numerical simulations. More precisely, we consider in this first part higher-order Allen-Cahn and Cahn-Hilliard equations with a regular potential and higher-order Allen-Cahn equation with a logarithmic potential. Moreover, we study higher-order anisotropic models and higher-order generalized Cahn-Hilliard equations, which have applications in biology, image processing, etc. We also consider the hyperbolic relaxation of higher-order anisotropic Cahn-Hilliard equations. In the second part, we develop semi-implicit and implicit semi-discrete, as well as fully discrete, schemes for solving the nonlinear partial differential equation, which describes both the elastic and electrostatic effects in an idealized MEMS capacitor. We analyze theoretically the stability of these schemes and the convergence under certain assumptions. Furthermore, several numerical simulations illustrate and support the theoretical results.

No description available.