This thesis develops a computationally efficient way of employing Gaussian wave packets to study laser-induced electron dynamics in atomic and molecular systems by directly solving the time-dependent Schrödinger equation (TDSE). First, we investigate charge migration (treating the nuclei classically), high-order harmonic generation (HHG), and single-isolated attosecond pulse generation in the Hydrogen molecular ion subjected to intense laser fields in a different range of frequencies with a basis of static coherent states (SCS). Then, seeking for a smarter way of constructing and guiding a minimal set of time-dependent basis functions, we introduce a fast and accurate approach for optimizing s-type Gaussian type orbitals (GTOs) and apply it to calculate electronic states of different 1D and 3D time-independent systems. Finally, we apply our optimization approach to time-dependent problems. With our approach we obtain excellent agreement with the exact results for HHG spectra of the 1D Hydrogen atom and molecular ion exposed to intense laser fields, which is not possible even with a much larger basis of static s-type GTOs. / Diese Arbeit sucht nach einem computereffizienten Ansatz für die Verwendung von Gaußschen Wellenpaketen zur Untersuchung der Quantenelektronendynamik in atomaren und molekularen laserinduzierten Systemen durch direkte Lösung der zeitabhängigen Schrödingergleichung (TDSE). Beginnend mit statischen kohärenten Zuständen (SCS) untersuchen wir die Ladungsmigration (wobei wir die Kerne klassisch behandeln), die Erzeugung von Oberwellen höherer Ordnung (HHG) und die Erzeugung von isolierten Attosekundenimpulsen im 3D-Wasserstoffmolekül-Ion \ih, das intensiven Laserfeldern in einem unterschiedlichen Frequenzbereich ausgesetzt ist. Auf der Suche nach einer intelligenteren Methode zur Konstruktion und Führung eines minimalen Satzes von Basisfunktionen stellen wir einen schnellen und genauen Ansatz zur Optimierung von Gauß-Orbitalen (GTOs) vom s-Typ vor und wenden ihn erfolgreich zur Berechnung gewünschter elektronischer Zustände verschiedener 1D- und 3D-Quantensysteme an. Letztendlich erweitern wir unseren Optimierungsansatz auf zeitabhängige Szenarien. Wir demonstrieren, wie diese Methode eine ausgezeichnete Übereinstimmung mit den exakten Ergebnissen in den HHG-Spektren des 1D-Wasserstoffatoms und des 1D-\ih, die intensiven Laserfeldern ausgesetzt sind, erzielt, wo die nicht optimierten s-Typ GTOs nicht übereinstimmen, selbst nach den ersten paar Harmonien im Plateaubereich.

Magnetism is a familiar phenomenon, and is applied in a variety of devices, from simple ones, like compasses, to the more sophisticated magnetic hard disk drives. It is also known that the magnetization of a material can change, for example by heating or by exposure to an external magnetic field. The dynamics of transitions between different magnetizations, however, is largely unknown, particularly in complex materials. To further the understanding of such dynamics, this thesis presents an observational study of the dynamics of laser induced demagnetization of permalloy (Ni0.8Fe0.2). Dynamics were studied with element-specificity, i.e. Fe and Ni were studied simultaneously, but separately, rather than studying the overall material. The study was conducted at the HELIOS laboratory at Uppsala University, which features equipment for the study of magnetization dynamics. Important concepts like high-order harmonic generation (HHG) and the transverse magneto-optic Kerr effect (T-MOKE) are discussed. HHG is the laser induced generation of high energy photons, and T-MOKE relates the intensity of reflected light to the magnetization of the reflecting material. The study revealed a very short demagnetization time, and maximum demagnetization of both elements was achieved within 1 picosecond. An onset delay in the demagnetizations of Fe and Ni of about 25 femtoseconds was also observed. Both phenomena have been previously reported. The results further imply that the magnetizations diverge over a 10 picosecond time interval after the onset of demagnetization, which has not been previously reported. The apparent divergence may be due to an unknown transient setup-related issue. The short demagnetization times, as well as the onset delay could potentially contribute to the development of a more complete theory of magnetization dynamics. / Magnetism är ett bekant fenomen, som utnyttjas i allt från enkla tillämpningar, som kompasser,till mer sofistikerade sådana, som hårddiskar. Det är också känt att magnetiseringen i ett material kan ändras, t.ex. genom upphettning eller genom att det utsätts för ett yttremagnetiskt fält. Dynamiken vid övergångar mellan olika magnetiseringstillstånd är dock ett relativt okänt ämne, i synnerhet när det kommer till komplexa material. För främjandet av en större förståelse inom ämnet, presenteras här en observationell studie av dynamiken för laserinducerad avmagnetisering av permalloy (Ni0.8Fe0.2). Dynamiken undersöktes med grundämnesspecificitet, d.v.s. Fe och Ni undersöktes samtidigt, men var för sig, snarare än att materialet undersöktes som helhet. Undersökningen gjordes vid HELIOS-laboratoriet vid Uppsala Universitet, som tillhandahåller utrustning för undersökning av magnetiseringsdynamik. Viktiga koncept diskuteras, såsom övertonsgenerering och den transversella magnetooptiska Kerr-effekten (T-MOKE). Övertonsgenerering innebär laserinducerad generering av högenergifotoner, och T-MOKE relaterar reflekterad intensitet till magnetiseringen i det reflekterande materialet. Undersökningen påvisade en mycket kort avmagnetiseringstid, och maximal avmagnetisering nåddes inom en pikosekund. En relativ tidsförskjutning mellan avmagnetiseringsförloppen för Fe och Ni om ungefär 25 femtosekunder observerades också. Båda fenomen har rapporterats tidigare. Resultatet visar även en divergens mellan magnetiseringsförloppen under ett tidsspann på 10 pikosekunder efter avmagnetiseringens början, vilket inte har rapporterats förr. Den skenbara divergensen kan bero på ett okänt, tillfälligt problem i uppställningen. Den korta avmagnetiseringstiden och den relativa tidsförskjutningen skulle kunna bidra till utvecklingen av en mer komplett teori för magnetiseringsdynamik.

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.

La simulation numérique est aujourd'hui un outils majeur dans la conception des objets aérodynamiques, que ce soit dans l'aéronautique, l'automobile, l'industrie navale, etc... Un des défis majeurs pour repousser les limites des codes de simulation est d'améliorer leur précision, tout en utilisant une quantité fixe de ressources (puissance et/ou temps de calcul). Cet objectif peut être atteint par deux approches différentes, soit en construisant une discrétisation fournissant sur un maillage donné une solution d'ordre très élevé, soit en construisant un schéma compact et massivement parallèlisable, de manière à minimiser le temps de calcul en distribuant le problème sur un grand nombre de processeurs. Dans cette thèse, nous tentons de rassembler ces deux approches par le développement et l'implémentation de Schéma Distribuant le Résidu (RDS) d'ordre très élevé et de compacité maximale. Ce manuscrit commence par un rappel des principaux résultats mathématiques concernant les Lois de Conservation hyperboliques (CLs). Le but de cette première partie est de mettre en évidence les propriétés des solutions analytiques que nous cherchons à approcher, de manière à injecter ces propriétés dans celles de la solution discrète recherchée. Nous décrivons ensuite les trois étapes principales de la construction d'un schéma RD d'ordre très élevé : - la représentation polynomiale d'ordre très élevé de la solution sur des polygones et des polyèdres; - la description de méthodes distribuant le résidu de faible ordre, compactes et conservatives, consistantes avec une représentation polynomiale des données de très haut degré. Parmi elles, une attention particulière est donnée à la plus simple, issue d'une généralisation du schéma de Lax-Friedrichs (\LxF); - la mise en place d'une procédure préservant la positivité qui transforme tout schéma stable et linéaire, en un schéma non linéaire d'ordre très élevé, capturant les chocs de manière non oscillante. Dans le manuscrit, nous montrons que les schémas obtenus par cette procédure sont consistants avec la CL considérée, qu'ils sont stables en norme $\L^{\infty}$ et qu'ils ont la bonne erreur de troncature. Même si tous ces développements théoriques ne sont démontrés que dans le cas de CLs scalaires, des remarques au sujet des problèmes vectoriels sont faites dès que cela est possible. Malheureusement, lorsqu'on considère le schéma \LxF, le problème algébrique non linéaire associé à la recherche de la solution stationnaire est en général mal posé. En particulier, on observe l'apparition de modes parasites de haute fréquence dans les régions de faible gradient. Ceux-ci sont éliminés grâce à un terme supplémentaire de stabilisation dont les effets et l'évaluation numérique sont précisément détaillés. Enfin, nous nous intéressons à une discrétisation correcte des conditions limites pour le schéma d'ordre élevé proposé. Cette théorie est ensuite illustrée sur des cas test scalaires bidimensionnels simples. Afin de montrer la généralité de notre approche, des maillages composés uniquement de triangles et des maillages hybrides, composés de triangles et de quandrangles, sont utilisés. Les résultats obtenus par ces tests confirment ce qui est attendu par la théorie et mettent en avant certains avantages des maillages hybrides. Nous considérons ensuite des solutions bidimensionnelles des équations d'Euler de la dynamique des gaz. Les résultats sont assez bons, mais on perd les pentes de convergence attendues dès que des conditions limite de paroi sont utilisées. Ce problème nécessite encore d'être étudié. Nous présentons alors l'implémentation parallèle du schéma. Celle-ci est analysée et illustrée à travers des cas test tridimensionnel de grande taille. / Numerical simulations are nowadays a major tool in aerodynamic design in aeronautic, automotive, naval industry etc... One of the main challenges to push further the limits of the simulation codes is to increase their accuracy within a fixed set of resources (computational power and/or time). Two possible approaches to deal with this issue are either to contruct discretizations yielding, on a given mesh, very high order accurate solutions, or to construct compact, massively parallelizable schemes to minimize the computational time by means of a high performance parallel implementation. In this thesis, we try to combine both approaches by investigating the contruction and implementation of very high order Residual Distribution Schemes (RDS) with the most possible compact stencil. The manuscript starts with a review of the mathematical theory of hyperbolic Conservation Laws (CLs). The aim of this initial part is to highlight the properties of the analytical solutions we are trying to approximate, in order to be able to link these properties with the ones of the sought discrete solutions. Next, we describe the three main steps toward the construction of a very high order RDS: - The definition of higher order polynomial representations of the solution over polygons and polyhedra; - The design of low order compact conservative RD schemes consistent with a given (high degree) polynomial representation. Among these, particular accest is put on the simplest, given by a generalization of the Lax-Friedrich's (\LxF) scheme; - The design of a positivity preserving nonlinear transformation, mapping first-order linear schemes onto nonlinear very high order schemes. In the manuscript, we show formally that the schemes obtained following this procedure are consistent with the initial CL, that they are stable in $L^{\infty}$ norm, and that they have the proper truncation error. Even though all the theoretical developments are carried out for scalar CLs, remarks on the extension to systems are given whenever possible. Unortunately, when employing the first order \LxF scheme as a basis for the construction of the nonlinear discretization, the final nonlinear algebraic equation is not well-posed in general. In particular, for smoothly varying solutions one observes the appearance of high frequency spurious modes. In order to kill these modes, a streamline dissipation term is added to the scheme. The analytical implications of this modifications, as well as its practical computation, are thouroughly studied. Lastly, we focus on a correct discretization of the boundary conditions for the very high order RDS proposed. The theory is then extensively verified on a variety of scalar two dimensional test cases. Both triangular, and hybrid triangular-quadrilateral meshes are used to show the generality of the approach. The results obtained in these tests confirm all the theoretical expectations in terms of accuracy and stability and underline some advantages of the hybrid grids. Next, we consider solutions of the two dimensional Euler equations of gas dynamics. The results obtained are quite satisfactory and yet, we are not able to obtain the desired convergence rates on problems involving solid wall boundaries. Further investigation of this problem is under way. We then discuss the parallel implementation of the schemes, and analyze and illustrate the performance of this implementation on large three dimensional problems. Due to the preliminary character and the complexity of these three dimensional problems, a rather qualitative discussion is made for these tests cases: the overall behavior seems to be the correct one, but more work is necessary to assess the properties of the schemes in three dimensions.

Ce manuscrit présente des expériences d’imagerie par diffraction réalisées en utilisant une source de rayonnement cohérent basée sur la génération d’harmoniques d’ordre élevé d’un laser Ti:Sa. Elles démontrent que cette source extrême-UV de laboratoire produit un nombre suffisant de photons par impulsion pour enregistrer une figure de diffraction d’objets tests en « simple tirs ». Le signal ainsi enregistré permet l’obtention d’une image de l’objet avec une résolution d’une centaine de nanomètres. Deux schémas sont utilisés pour reconstruire l’objet : le premier utilise un algorithme itératif de reconstruction de la phase perdue pendant la détection de la figure de diffraction ; le second utilise une configuration holographique par transformée de Fourier. Les travaux réalisés comportent deux parties. La première concerne l’optimisation de la source harmonique et inclut une étude expérimentale d’un dispositif de filtrage spatial du faisceau laser de génération par propagation dans une fibre creuse. La seconde partie présente les expériences d’imagerie par diffraction, et notamment une démonstration du schéma holographique HERALDO qui est une extension de l’holographie par transformée de Fourier à des références en forme de polygones. L’utilisation de ces références « étendues » a pour avantage d’optimiser l’enregistrement holographique tout en conservant une reconstruction directe et sans ambigüité de l’objet. Une analyse signal-sur-bruit ainsi qu’une comparaison des reconstructions d’hologramme pour différentes formes de références sont effectuées. / This manuscript presents diffraction imaging experiments performed using a source of coherent radiation based on high order harmonics generation of a Ti:Sa laser. They demonstrate that this laboratory size XUV source produces a number of photons per pulse sufficient to record the diffraction pattern of test objects in « single shot ». The signal thus recorded allows obtaining an image of the object with a resolution of around 100 nanometers. Two schemes are used to reconstruct the object: the first one uses an iterative algorithm to retrieve the phase lost during the detection of the diffraction pattern; the second uses a configuration of Fourier transform holography. The work presented here is separated in two parts. The first one concerns the optimization of the harmonic source, including an experimental study of a spatial filtering device for laser beams by propagation in a hollow core fiber. The second part deals with the diffraction imaging experiments. In particular, I present a demonstration of the holographic scheme HERALDO, which is an extension of the Fourier transform holography with polygonal references. The use of these « extended » references allows the optimization of the holographic recording while maintaining a direct and non-ambiguous reconstruction of the object. An analysis of signal-to-noise ratio and a comparison of hologram reconstructions for different types of references are performed.

L’imagerie du sous-sol par équations d’onde est une application de l’ingénierie pétrolière qui mobilise des ressources de calcul très importantes. On dispose aujourd’hui de calculateurs puissants qui rendent accessible l’imagerie de régions complexes mais des progrès sont encore nécessaires pour réduire les coûts de calcul et améliorer la qualité des simulations. Les méthodes utilisées aujourd’hui ne permettent toujours pas d’imager correctement des régions très hétérogènes 3D parce qu’elles sont trop coûteuses et /ou pas assez précises. Les méthodes d’éléments finis sont reconnues pour leur efficacité à produire des simulations de qualité dans des milieux hétérogènes. Dans cette thèse, on a fait le choix d’utiliser une méthode de Galerkine discontinue (DG) d’ordre élevé à flux centrés pour résoudre l’équation des ondes acoustiques et on développe un schéma d’ordre élevé pour l’intégration en temps qui peut se coupler avec la technique de discrétisation en espace, sans générer des coûts de calcul plus élevés qu’avec le schéma d’ordre deux Leap-Frog qui est le plus couramment employé. Le nouveau schéma est comparé au schéma d’ordre élevé ADER qui s’avère plus coûteux car il requiert un plus grand nombre d’opérations pour un niveau de précision fixé. De plus, le schéma ADER utilise plus de mémoire, ce qui joue aussi en faveur du nouveau schéma car la production d’images du sous-sol consomme beaucoup de mémoire et justifie de développer des méthodes numériques qui utilisent la mémoire au minimum. On analyse également la précision des deux schémas intégrés dans un code industriel et appliqués à des cas test réalistes. On met en évidence des phénomènes de pollution numériques liés à la mise en oeuvre d'une source ponctuelle dans le schéma DG et on montre qu'on peut éliminer ces ondes parasites en introduisant un terme de pénalisation non dissipatif dans la formulation DG. On finit cette thèse en discutant les difficultés engendrées par l'utilisation de schémas numériques dans un contexte industriel, et en particulier l'effet des calculs en simple précision. / Oil engineering uses a wide variety of technologies including imaging wave equation which involves very large computing resources. Very powerful computers are now available that make imaging of complex areas possible, but further progress is needed both to reduce the computational cost and improve the simulation accuracy. The current methods still do not allow to image properly heterogeneous 3D regions because they are too expensive and / or not accurate enough. Finite element methods turn out to be efficient for producing good simulations in heterogeneous media. In this thesis, we thus chose to use a high order Discontinuous Galerkin (DG) method based upon centered fluxes to solve the acoustic wave equation and developed a high-order scheme for time integration which can be coupled with the space discretization technique, without generating higher computational cost than the second-order Leap Frog scheme which is the most widely used . The new scheme is compared to the high order ADER scheme which is more expensive because it requires a larger number of computations for a fixed level of accuracy. In addition, the ADER scheme uses more memory, which also works in favor of the new scheme since producing subsurface images consumes lots of memory and justifies the development of low-memory numerical methods. The accuracy of both schemes is then analyzed when they are included in an industrial code and applied to realistic problems. The comparison highlights the phenomena of numerical pollution that occur when injecting a point source in the DG scheme and shows that spurious waves can be eliminated by introducing a non-dissipative penalty term in the DG formulation. This work ends by discussing the difficulties induced by using numerical methods in an industrial framework, and in particular the effect of single precision calculations.

Le sujet de ce travail de thèse est l’application de la méthode de caractéristiques longues (MOC) pour résoudre l’équation du transport des neutrons pour des géométries à trois dimensions extrudées. Les avantages du MOC sont sa précision et son adaptabilité, le point faible était la quantité de ressources de calcul requises. Ce problème est même plus important pour des géométries à trois dimensions ou le nombre d’inconnues du problème est de l’ordre de la centaine de millions pour des calculs d’assemblage.La première partie de la recherche a été dédiée au développement des techniques optimisées pour le traçage et la reconstruction à-la-volé des trajectoires. Ces méthodes profitent des régularités des géométries extrudées et ont permis une forte réduction de l’empreinte mémoire et une réduction des temps de calcul. La convergence du schéma itératif a été accélérée par un opérateur de transport dégradé (DPN) qui est utilisé pour initialiser les inconnues de l’algorithme itératif and pour la solution du problème synthétique au cours des itérations MOC. Les algorithmes pour la construction et la solution des opérateurs MOC et DPN ont été accélérés en utilisant des méthodes de parallélisation à mémoire partagée qui sont le plus adaptés pour des machines de bureau et pour des clusters de calcul. Une partie importante de cette recherche a été dédiée à l’implémentation des méthodes d’équilibrage la charge pour améliorer l’efficacité du parallélisme. La convergence des formules de quadrature pour des cas 3D extrudé a aussi été explorée. Certaines formules profitent de couts négligeables du traitement des directions azimutales et de la direction verticale pour accélérer l’algorithme. La validation de l’algorithme du MOC a été faite par des comparaisons avec une solution de référence calculée par un solveur Monte Carlo avec traitement continu de l’énergie. Pour cette comparaison on propose un couplage entre le MOC et la méthode des Sous-Groupes pour prendre en compte les effets des résonances des sections efficaces. Le calcul complet d’un assemblage de réacteur rapide avec interface fertile/fissile nécessite 2 heures d’exécution avec des erreurs de quelque pcm par rapport à la solution de référence.On propose aussi une approximation d’ordre supérieur du MOC basée sur une expansion axiale polynomiale du flux dans chaque maille. Cette méthode permet une réduction du nombre de mailles (et d’inconnues) tout en gardant la même précision.Toutes les méthodes développées dans ce travail de thèse ont été implémentées dans la version APOLLO3 du solveur de transport TDT. / The topic of our research is the application of the Method of Long Characteristics (MOC) to solve the Neutron Transport Equation in three-dimensional axial geometries. The strength of the MOC is in its precision and versatility. As a drawback, it requires a large amount of computational resources. This problem is even more severe in three-dimensional geometries, for which unknowns reach the order of tens of billions for assembly-level calculations.The first part of the research has dealt with the development of optimized tracking and reconstruction techniques which take advantage of the regularities of three-dimensional axial geometries. These methods have allowed a strong reduction of the memory requirements and a reduction of the execution time of the MOC calculation.The convergence of the iterative scheme has been accelerated with a lower-order transport operator (DPN) which is used for the initialization of the solution and for solving the synthetic problem during MOC iterations.The algorithms for the construction and solution of the MOC and DPN operators have been accelerated by using shared-memory parallel paradigms which are more suitable for standard desktop working stations. An important part of this research has been devoted to the implementation of scheduling techniques to improve the parallel efficiency.The convergence of the angular quadrature formula for three-dimensional cases is also studied. Some of these formulas take advantage of the reduced computational costs of the treatment of planar directions and the vertical direction to speed up the algorithm.The verification of the MOC solver has been done by comparing results with continuous-in-energy Monte Carlo calculations. For this purpose a coupling of the 3D MOC solver with the Subgroup method is proposed to take into account the effects of cross sections resonances. The full calculation of a FBR assembly requires about 2 hours of execution time with differences of few PCM with respect to the reference results.We also propose a higher order scheme of the MOC solver based on an axial polynomial expansion of the unknown within each mesh. This method allows the reduction of the meshes (and unknowns) by keeping the same precision.All the methods developed in this thesis have been implemented in the APOLLO3 version of the neutron transport solver TDT.

In this work we review 12 years of developments in the field of residual based discretizations for hyperbolic problems and their application to the solution of the shallow water equations. Fundamental concepts related to the topic are recalled and he construction of second and higher order schemes for steady problems is presented. The generalization to time dependent problems by means of multi-step implicit time integration, space-time, and genuinely explicit techniques is thoroughly discussed. Finally, the issues of C-property, super consistency, and wetting/drying are analyzed in this framework showing the power of the residual based approach.

No description available.

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.