Les études présentées dans ce mémoire portent sur la simulation numérique et l'analyse physique des milieux granulaires immergés dans un fluide. Des développements numériques ont été réalisés pour coupler la méthode Lattice Boltzmann pour la dynamique du fluide avec la méthode Contact Dynamics en 2D et avec la méthode Molecular Dynamics en 3D pour la dynamique des grains. Ces outils numériques ont été utilisés pour étudier l'initiation des avalanches sur un plan incliné en fonction de la compacité initiale et de l'angle d'inclinaison en 3D. Les résultats sont en bon accord quantitatif avec les expériences et ont permis de mettre en évidence la stabilisation de la pente granulaire par une pression négative du fluide interstitielle induite par la dilatance, et l'évolution spatiotemporelle des grandeurs telles que la compacité et la déformation de cisaillement. Ces évolutions dans la phase de fluage qui précède la rupture de pente ont pu être mises à l'échelle par un modèle théorique incorporant la loi de Darcy et l'effet de la dilatance sur l'angle de frottement interne. L'analyse de la texture granulaire a révélé la distortion du réseau des contacts pendant le fluage et la saturation de l'anisotropie comme un critère de rupture. La propagation des avalanches granulaires a été étudiée dans une configuration 2D pour deux géométries différentes : 1) l'effondrement et l'étalement d'une colonne sous son propre poids, 2) l'étalement d'une pente sous l'effet d'une énergie cinétique injectée. Nous avons en particulier montré que la distance et la durée d'étalement obéissent à des lois de puissance en fonction du rapport d'aspect initial ou de l'énergie injectée. Le fluide exerce deux effets contradictoires : réduire les temps de relaxation et lubrifier les contacts. Ces effets ont été analysés dans le régime visqueux en fonction des conditions initiales et la viscosité du fluide.

La spectroscopie ro-vibrationelle de haute résolution de l’ion moleculaire H2+ par REMPD requiert une source sélective en état interne. Dans cette thèse nous présentons notre travail de conception et de réalisation d’une telle source utilisant la photo-ionisation multiphotonique résonante (REMPI) de H2 à l’aide d’un laser pulsé à 303 nm.Dans un deuxième temps nous présentons nos simulations numériques de refroidissement sympathique dans un piège de Paul linéaire, avec pour application principale le projet GBAR qui implique le refroidissement sympathique d’un ion d’antimatière H̄+ par des ions Be+ refroidis par laser. Nous avons dévelopé un code GPU utilisant un pas de temps variable permettant de décrire les interactions coulombiennes de façon efficace. Nous discutons de l’influence du chauffage RF et de lois d’échelles entre le temps de capture, l’énergie initiale et le nombre d’ions dans le cristal. Nous montrons que le refroidissement sympathique de H̄ requis pour GBAR pourrait fonctionner avec un cristal dissymétrique de Be+/HD+ qui semble plus efficace que le Be+ seul. Nous montrons qu’avec un tel cristal la capture du H̄+ pourrait être détectée expérimentalement par analyse de Fourier des données de fluorescence. / The high-resolution ro-vibrational spectroscopy of the H2+ molecular ion by resonance-enhanced multiphoton dissociation (REMPD) requires a state-selective source of H2+. In this thesis we present work on a functional state selective H2+ ion source using resonance enhanced multiphoton ionisation (REMPI) with a 303 nm pulsed laser.The second part of the thesis presents numerical simulations of sympathetic cooling in linear RF traps, whose main application is the GBAR project (Gravitational Behaviour of Antihydrogen at Rest) which involves sympathetic cooling of an antimatter H̄+ ion by laser-cooled Be+ ions. We developed a GPU code using a variable timestep technique enabling a fast description of Coulomb interactions. We discuss the influence of RF heating and scaling laws between cooling times, initial energy and ion numbers in the cooling crystal. We show that the H̄ sympathetic cooling step of GBAR could be feasible using a rotationally asymmetric two-component Be+/HD+ crystal which appears more effective than a single-component Be+ crystal. We find that the H̄+ ion’s capture by this crystal could be detected experimentally by Fourier analysis of the fluorescence data.

Dans le contexte du calcul scientique, l'équilibrage de la charge est un problème crucial qui conditionne la performance des simulations numériques parallèles. L'objectif est de répartir la charge de travail entre un nombre de processeurs donné, afin de minimiser le temps global d'exécution. Une stratégie populaire pour résoudre ce problème consiste à modéliser la simulation à l'aide d'un graphe et à appliquer des algorithmes de partitionnement. En outre, les simulations numériques tendent à se complexifier, notamment en mixant plusieurs codes représentant des physiques différentes ou des échelles différentes. On parle alors de couplage de codes multi-physiques ou multi-échelles. Dans ce contexte, le problème de l'équilibrage de charge devient également plus difficile, car il ne s'agit plus d'équilibrer chacun des codes séparément, mais l'ensemble de ces codes pris dans leur globalité. Dans ce travail, on propose de resoudre ce problème en utilisant le modèle de partitionnement à sommets fixes qui pourrait représenter efficacement les contraintes supplémentaires imposées par les codes couplés (co-partitionnement). Nous avons donc développé un algorithme direct de partitionnement de graphe qui gère des sommets fixes. L'algorithme a été implémenté dans le partitionneur Scotch et une série d'expériences ont été menées sur la collection des graphes DIMACS. Ensuite nous avons proposé trois algorithmes de co-partitionnement qui respectent les contraintes issues des codes couplés respectifs. Nous avons egalement validé nos algorithmes par une étude expérimentale en comparant nos méthodes aux strategies actuelles sur des cas artificiels ainsi que sur des codes réels couplés. / Load balancing is an important step conditioning the performance of parallel applications. The goal is to distribute roughly equal amounts of computational load across a number of processors, while minimising interprocessor communication. A common approach to model the problem is based on graph structures and graph partitioning algorithms. Moreover, new challenges involve the simulation of more complex physical phenomena, where different parts of the computational domain exhibit different physical behavior. Such simulations follow the paradigm of multi-physics or multi-scale modeling approaches. Combining such different models in massively parallel computations is still a challenge to reach high performance. Additionally, traditional load balancing algorithms are often inadequate, and more sophisticated solutions should be explored. In this thesis, we propose new graph partitioning algorithms that balance the load of such simulations, refered to as co-partitioning. We formulate this problem with the use of graph partitioning with initially fixed vertices which we believe represents efficiently the additional constraints of coupled simulations. We have therefore developed a direct algorithm for graph partitioning that manages successfully problems with fixed vertices. The algorithm is implemented inside Scotch partitioner and a series of experiments were carried out on the DIMACS graph collection. Moreover we proposed three copartitioning algorithms that respect the constraints of the respective coupled codes. We finally validated our algorithms by an experimental study comparing our methods with current strategies on artificial cases and on real-life coupled simulations.

Ce travail de thèse est consacré à l'étude de la corrélation entre la microstructure des matériaux granulaires cimentés, la morphologie des fissures qui peuvent y apparaître et leur perméabilité apparente. Une approche numérique par discrétisation sur réseau a permis de prendre en compte la structure polyphasique hétérogène de ces matériaux et d'étudier les mécanismes qui contrôlent leur comportement et leur rupture. Trois régimes de rupture ont été identifiés en fonction de la fraction volumique de la phase cimentaire et de l'adhésion aux interfaces inclusions/matrice cimentaire. Dans le cadre de la sûreté nucléaire, et dans une optique d'analyse d'étanchéité sur des bétons fissurés, une méthodologie de génération de microstructures représentatives du béton, combinant les fractions volumiques des phases et les distributions des tailles des grains, a été proposée. Cette méthodologie permet d'analyser le réseau de fissures qui peut y apparaitre sous chargement en fonction de la microstructure. En particulier, la tortuosité des fissures est analysée en fonction de la fraction volumique des inclusions et de la distribution de leurs diamètres. Enfin, la perméabilité apparente des des fissures a été étudiée par la simulation de l'écoulement d'un liquide par la methode Lattice Boltzmann. Une corrélation microstructure-tortuosité-perméabilité a été ainsi obtenue. Les outils de modélisation et d'analyse proposés sont suffisamment génériques pour permettre de prendre en compte la complexité microstructurelle d'autres types de matériaux polyphasiques et leurs évolutions. / In this work, we investigate the relationships between the microstructure of cemented granular materials, morphology of cracks that may appear under tensile loading and their apparent permeability. In order to simulate the complexity of the multiphase heterogeneous structure of these materials, a numerical approach based on a lattice discretization was developed and used to study the mechanisms that control their behavior and rupture. Three distinct regimes of crack propagation were evidenced in terms of the combined influence of the matrix volume fraction and particle/matrix adherence. In the context of nuclear safety and in view of analyzing cracked concrete toughness, a methodology is proposed for generating representative microstructures with controlled phase volume fractions and particle size distributions. The cracks obtained under tensile loading are analyzed and a relationship is obtained between the microstructure of concrete and the tortuosity of the cracks. The permeability of cracked numerical microstructures was obtained by the simulation of liquid flow through the cracks by means of the Lattice Boltzmann method. Finally, a microstructure-tortuosity-permeability relation is proposed. The modeling and analysis tools developed in this work are generic enough to be applied to other complex multiphase heterogeneous materials.

Cette thèse s'inscrit dans le cadre des énergies renouvelables au sein du programme HARVEST du laboratoire LEGI, qui consiste à développer un concept d'hydrolienne de type Darrieus. L'hydrolienne peut être équipée d'un dispositif appelé carénage afin de transformer une portion plus grande de l'énergie cinétique contenue dans le courant d'eau en électricité. Les travaux présentés se sont focalisés sur ces systèmes de carénage, autour de trois axes : l'explication du principe de fonctionnement hydrodynamique du carénage, la quantification des performances de l'hydrolienne carénée et la mise en évidence des grandeurs géométriques clés du carénage permettant d'améliorer ou d'optimiser la performance du système. Toutes les études ont été réalisées à l'aide des calculs RANS 2D et des données expérimentales mises à disposition et, comparées aux résultats obtenus pour une hydrolienne non-carénée. / The general context of the present thesis is renewable energies within the HARVEST program initialized at LEGI laboratory, which consists in developing a Darrieus-type water current turbine (WCT). The WTC can be equipped with a channelling device which allows transforming a bigger amount of the kinetic energy contained in the flowstream into electricity. The present work is focused on the channelling devices. Studies concern three main topics: the explanation of the channelling device hydrodynamic functioning, the evaluation of the performance of the shrouded WCT and the revealing of the system geometrical parameters which allow its improvement or optimisation. All studies have been carried out by 2D RANS calculations and available experimental data, and have been compared to bare WTC results.

Les vitesses particulières de galaxies fournissent des informations importantes sur notre mouvement par rapport aux grandes structures observées et peuvent être utilisées afin de contraindre la distribution de matière noire sous-jacente. En ce basant sur cette approche, la technique des réalisations contraintes permet de calculer des simulations numériques qui ressemblent à l'Univers Local observé. Ceci fourni un laboratoire numérique puissant pour étudier la dynamique, la formation et l' évolution des structures bien connues de l'Univers Local. L' étape cruciale est de générer, à partir des données observationnelles, des conditions initiales appropriées. Nous présentons ici des améliorations de la technique des simulations contraintes, accompagnées d'un code numérique nouvellement développé et hautement optimisé, qui peut ingérer les énormes ensemble de données observationnelles. Les galaxies évoluent à partir du collapse gravitationnel des surdensités primordiales de l'univers très jeune; leurs mouvements créent un champ de déplacement à grande échelle. Une source majeure d'erreurs systématiques dans les simulations contraintes est produite si l'on ne tient pas compte de cet effet. Afin de dépasser cette limitation, nous avons développé la méthode de reconstruction par approximation inverse de Zeldovich (RZA). La RZA permet de reconstruire les déplacements et les positions initiales des galaxies observées et de générer une estimation significativement meilleure des conditions initiales de l'Univers Local. Cette méthode est intensivement testée sur des données de simulations. Nous étudions aussi l'influence de la qualité des données et de diverses erreurs observationnelles et/ou systématiques. Nous démontrons qu'avec la technique RZA, la qualité de la reconstruction des champs de densité et de vitesse est drastiquement améliorée. La position des objets dans les simulations contraintes évoluées sont retrouvées plus précisément et sur un plus grand intervalle de masses / Galaxy peculiar velocities provide valuable information about our motion with respect to theobserved large-scale structure and can be used to constrain the underlying dark matter distribution.Based on this approach, the technique of constrained realisations allows us to run numericalsimulations that resemble the observed Local Universe. This provides a powerful numerical laboratoryto study the dynamics, formation and evolution of structure in the Local Universe. Thecrucial step is to generate appropriate initial conditions from the observational data.We present here improvements on the technique of constrained simulations, along with anewly developed highly optimised numerical code that can handle the upcoming large observationaldatasets. Galaxies evolve from the gravitational collapse of primordial overdensities in theearly Universe; their motion leads to a large-scale displacement field. A major source of systematicerrors in constrained simulations arises by not accounting for this effect. To overcome thislimitation, we develop the Reverse Zeldovich Approximation (RZA) reconstruction method. TheRZA allows to reconstruct displacements and initial positions of observed galaxies and generatea significantly better estimate of the initial conditions of the Local Universe. This method isextensively tested on simulation data. We also study the influence of data quality and variousobservational and systematic errors. We show that with the RZA technique, the reconstructionquality of the density and velocity fields improves significantly. The positions of objects in theevolved constrained simulations are recovered more accurately and over a wider range of masses. / Die Pekuliargeschwindigkeiten von Galaxien liefern wertvolle Informationen über ihre Bewegungrelativ zu den beobachteten großräumigen Strukturen und können verwendet werden, um diezugrundeliegende Verteilung der dunklen Materie abzuschätzen. Auf dieser Grundlage könnenmit der Methode der “Constrained Realisations” numerische Simulationen durchgeführt werden,die die Materieverteilung im Lokalen Universum widerspiegeln und somit leistungsfähige numerischeExperimente ermoeglichen, um die Dynamik, Entstehung und Evolution von Strukturenim lokalen Universum zu untersuchen. Eine wichtige Voraussetzung für diese Experimente ist es,korrekte Anfangsbedingungen aus den Beobachtungsdaten abzuleiten.In dieser Arbeit werden Verbesserungen der “Constrained Simulations”-Technik sowie ein fürdiesen Zweck neu entwickelter, hochgradig optimierter numerischer Code vorgestellt, welcheres ermöglicht, die kommenden, sehr umfangreichen Beobachtungsdaten zu verarbeiten. Galaxienentwickeln sich aus dem Gravitationskollaps von primordialen Überdichten im frühen Universum.Deren Bewegung führt zu einer großräumigen Dislokation, deren Vernachlässigung die größtesystematische Fehlerquelle für “Constrained Simulations” darstellt. An dieser Stelle knüpft dieReverse-Zeldovich-Näherung (RZA) an, die hier vorgestellt wird. Mithilfe der RZA können wirdiese Dislokation abschätzen und somit die ursprünglichen Positionen der beobachteten Galaxienrekonstruieren. Dadurch erhält man wesentlich genauere Anfangsbedingungen fuer die Entwicklungdes Lokalen Universums. Diese Methode wird ausführlich an Simulationsdaten getestet.Der Einfluß der Datenqualität und verschiedener Beobachtungs- und systematischer Fehler wirdeingehend untersucht. Die Rekonstruktionsgenauigkeit der Dichte- und Geschwindigkeitsfelderkann durch den Einsatz der RZA signifikant erhöht werden. In unseren Tests stimmen die Positionenvon Objekten in den mit RZA erzeugten “Constrained Simulations” wesentlich besser mitden ursprünglichen überein, als ohne die RZA, und das für einen deutlich größeren Bereich vonMassen

Le milieu interstellaire (MIS) est un système extrêmement complexe. Il correspond à une échelle intermédiaire entre les étoiles et les galaxies. Le gaz interstellaire est présent dans toute la galaxie, remplissant l’espace entre les étoiles. Une grande diversité de processus couplés, comme la gravité, le champs magnétiques, la turbulence et la chimie, participe à son évolution, faisant de la modélisation du MIS un problème ardu. Une description correcte du MIS nécessite un bon traitement des équations de la magnetohydrodynamique (MHD), de la gravité, du bilan thermique et de l’évolution chimique à l’intérieur du nuage moléculaire.L’objectif de ce travail de thèse est une meilleure compréhension de la formation et de l’évolution des nuages moléculaires, et plus particulièrement de la transition du gaz atomique en gaz moléculaire. Nous avons réalisé des simulations numériques de la formation des nuages moléculaires et de la formation de l’hydrogène moléculaire sous l’influence de la gravité et de la turbulence MHD, en utilisant des estimations précises de l’écrantage par les poussières et de l’auto-écrantage par la molécule H2. Ceci a été calculé grâce à une méthode en arbre, à même de fournir une rapide estimation des densités de colonne.Nous avons trouvé que l’hydrogène moléculaire se forme plus rapidement que prévu par les estimations classiques du fait de l’augmentation de densité locale provoquée par les fluctuations turbulentes du gaz. L’hydrogène moléculaire, formé à des densités plus élevées, peut alors migrer vers les régions plus chaudes et moins denses.Les densités de colonne totale d’hydrogène moléculaire montrent que la transition HI-H2 se produit à des densités de colonne de quelques 10^20 cm−2. Nous avons calculé les populations des niveaux rotationnels de H2 à l’équilibre thermique et intégré le long de plusieurs lignes de visée. Ces résultats reproduisent bien les valeurs observées par Copernicus et FUSE, suggérant que la transition observée et les populations excitées pourraient être une conséquence de la structure multi-phasique des nuages moléculaires. Comme la formation de H2 précède la formation des autres molécules, le H2 chaud pourrait permettre le développement d’espèces endothermiques et éventuellement expliquer certains aspects de la richesse moléculaire observée dans l’ISM. / The interstellar medium (ISM) is a highly complex system. It corresponds to an intermediate scale between stars and galaxies. The interstellar gas is present throughout the galaxy, filling the volume between stars. A wide variety of coupled processes, such as gravity, magnetic fields, turbulence and chemistry, participate in its evolution, making the modeling of the ISM a challenging problem. A correct description of the ISM requires a good treatment of the magnetohydrodynamics (MHD) equations, gravity, thermal balance, and chemical evolution within the molecular clouds.This thesis work aims at a better understanding of the formation and evolution of molecular clouds, specially how they become "molecular", paying particular attention to the transition HI-to-H2. We have performed ideal MHD simulations of the formation of molecular clouds and the formation of molecular hydrogen under the influence of gravity and turbulence, using accurate estimates for the shielding effects from dust and the self-shielding for H2, calculated with a Tree-based method, able to provide fast estimates of column densities.We find that H2 is formed faster than predicted by the usual estimates due to local density enhancements created by the gas turbulent motions. Molecular hydrogen, formed at higher densities, could then migrate toward low density warmer regions.Total H2 column densities show that the HI-to-H2 transition occurs at total column densities of a few 10^20 cm−2. We have calculated the populations of rotational levels of H2 at thermal equilibrium, and integrated along several lines of sight. These two results reproduce quite well the values observed by Copernicus and FUSE, suggesting that the observed transition and the excited populations could arise as a consequence of the multi-phase structure of molecular clouds. As H2 formation is prior to further molecule formation, warm H2 could possibly allow the development of a warm chemistry, and eventually explain some aspects of the molecular richness observed in the ISM.

Les fluides à seuil sont des matériaux utilisés dans de nombreux procédés industriels au cours desquels ils sont mis en écoulement via leur interaction avec des outils solides. Afin de mieux comprendre comment ils s'étalent sur des surfaces, nous avons étudié d'une part leur enduction verticale sur une plaque mince et d'autre part leur étalement horizontal à l'aide d'une lame. Ces deux études se sont appuyées sur des mesures macroscopiques permettant de suivre la force appliquée à l'outil et des mesures microscopiques de PIV pour déterminer les champs de vitesse dans le matériau. Après confirmation de leur validité par comparaison avec nos résultats expérimentaux, des simulations numériques basées sur la programmation conique ont permis de préciser les écoulements en jeu. Pour l'enduction par trempage, nous observons que dans notre gamme d'étude, le dépôt de fluide à seuil sur la plaque est millimétrique et d'épaisseur constante, sauf aux extrémités. Selon le rapport des forces visqueuses et plastiques, l'enduction est régie par un équilibre "seuillo-gravitaire" ou "visco-gravitaire". Nous avons ensuite caractérisé l'écoulement généré dans le bain par le déplacement de la plaque afin d'expliciter l'origine du phénomène d'enduction. Les simulations numériques précisent la forme de cet écoulement. L'étalement horizontal dans un canal à l'aide d'une lame mince conduit au déplacement d'un amas de fluide faiblement cisaillé par rapport à une région de fluide au repos, par l'intermédiaire d'une couche cisaillée. Nous montrons qu'un modèle simple permet de prédire la dynamique de croissance de l'amas et la relie à la force normale à la lame / Yield stress fluids are used in various industrial processes in which solid tools make them flow. To have a better understanding on how they spread on surfaces, we first studied the vertical coating on a thin plate and then the horizontal blade-coating. Both studies are based on macroscopic measures recording the force applied on the plate and on the microscopic determination on the velocity fields in the fluid by PIV. Numerical simulations based on cone programming, which validity had been confirmed in comparison to our experimental results, specified the flows at stake. In dip-coating, we observed than within our framework, the yield stress fluid deposit on the plate is millimetric and of constant thickness, except on the tips. According to the value of the viscous over plastic forces ratio, the coating is led by a "yield-gravity" or a "visco-gravity" balance. Then, we characterized the flow generated in the bath by the plate displacement in order to clarify the origin of the coating phenomenon. The numerical simulations specified the flow characteristics. The horizontal blade-coating in a channel involves the displacement of a weakly sheared cluster of fluid in relation to a part of fluid at rest, through a sheared layer. We showed that the growing dynamic of the cluster is described with a simple model which links it to the normal force to the plate

Dans cette thèse nous nous intéressons au problème de l’impact d’une goutte sur une surface solide. Nous proposons pour cela de nous placer dans un cadre plus général en utilisant les analogies existantes avec d’autres problèmes d’impact. Dans la première partie de ce manuscrit nous proposons de revisiter le problème de l’impact de goutte pour les temps courts à la lumière de son problème dual à savoir l’impact d’un objet solide dans un bain liquide lorsque l’inertie est l’effet dominant. De cette analogie est déduit un modèle reposant sur la théorie des écoulements potentiels. L’analyse asymptotique nous permet de dégager à l’ordre dominant les mécanismes essentiels de ce problème puis nous mettons en évidence la structure autosimilaire des champs de pression et de vitesse induits par l’impact. La structure de la couche limite est également étudiée. Les prédictions théoriques issues de ce modèle sont comparées à des solutions numériques obtenues à l’aide d’un solveur des équations de Navier-Stokes. Nous étudions ensuite les temps intermédiaires de l’impact, correspondants au moment où la solution autosimilaire cesse d’être valide et nous déterminons les causes de cette transition. Dans la troisième partie nous étudions un cas particulier d’évolution aux temps longs en revisitant le problème de l’impact d’une goutte sur un disque de même taille. Nous obtenons les solutions analytiques pour les champs de pression et de vitesse à l’instant initial et nous proposons ensuite différentes directions de recherche pour l’étude de l’évolution de la nappe liquide induite par l’impact. Nous finissons ce manuscrit par une brève introduction aux impacts de goutte de fluides à seuil. / In this thesis we consider the problem of drop impact onto a solid surface. In order to study this phenomenon we consider a more general framework by using analogies with some other impact problems which are a priori very different. In the first part of the thesis we propose to revisit the inertia-dominated drop impact problem for short times at the light of the dual problem defined by the impact of a solid object onto a liquid bath. We deduce from this analogy a model based on potential flow theory. Then asymptotic analysis is used to determine the essential mechanisms of the problem at leading order. This approach reveal a self-similar structure both for the velocity field and the pressure field induced by the impact. The structure of the boundary layer is also studied. Theoretical predictions deduced from this model are compared with numerical solutions obtained with the Navier-Stokes multiphase flow solver Gerris. Then we study the impact for intermediates times which correspond to the period of the breakdown of the self-similar solution. The origin of the transition is determined by using new numerical experiments. In a third part we propose to study a particular case of long time evolution by revisiting the problem of drop impact onto a solid target matching its own size. We obtain analytical solutions for pressure and velocity fields at initial time by using pressure impulse theory and we propose few research directions for the study of the evolution of the liquid sheet induced by the impact. This thesis ends with a brief introduction to drop impact of Bingham fluids.

Ce travail porte sur le développement et la vérification d’une nouvelle approche de coordonnées alignées FCI (Flux-Coordinate Independent), qui tire partie de l’anisotropie du transport dans un plasma immergé dans un fort champ magnétique. Sa prise en compte dans les codes numériques permet de réduire grandement le coût de calcul nécessaire pour une précision donnée. Une particularité de l’approche nouvellement développée dans ce manuscrit est en particulier sa capacité à traiter, pour la première fois, des configurations avec point X. Toutes ces analyses ont été conduites avec FENICIA, code modulaire entièrement développé dans le cadre de cette thèse, et permettant la résolution d’une classe de modèles génériques. En résumé, la méthode développée dans ce travail est validée. Elle peut s’avérer pertinente pour un large champ d’application dans le contexte de la fusion magnétique. Il est montré dans cette thèse que cette technique devrait pouvoir s’appliquer aussi bien aux modèles fluides que gyrocinétiques de turbulence, et qu’elle permet notamment de surmonter un des problèmes fondamentaux des techniques actuelles, qui peinent à traiter de manière précise la traversée de la séparatrice. / The primary thrust of this work is the development and implementation of a new approach to the problem of field-aligned coordinates in magnetized plasma turbulence simulations called the FCI approach (Flux-Coordinate Independent). The method exploits the elongated nature of micro-instability driven turbulence which typically has perpendicular scales on the order of a few ion gyro-radii, and parallel scales on the order of the machine size. Mathematically speaking, it relies on local transformations that align a suitable coordinate to the magnetic field to allow efficient computation of the parallel derivative. However, it does not rely on flux coordinates, which permits discretizing any given field on a regular grid in the natural coordinates such as (x, y, z) in the cylindrical limit. The new method has a number of advantages over methods constructed starting from flux coordinates, allowing for more flexible coding in a variety of situations including X-point configurations. In light of these findings, a plasma simulation code FENICIA has been developed based on the FCI approach with the ability to tackle a wide class of physical models. The code has been verified on several 3D test models. The accuracy of the approach is tested in particular with respect to the question of spurious radial transport. Tests on 3D models of the drift wave propagation and of the Ion Temperature Gradient (ITG) instability in cylindrical geometry in the linear regime demonstrate again the high quality of the numerical method. Finally, the FCI approach is shown to be able to deal with an X-point configuration such as one with a magnetic island with good convergence and conservation properties.