Dans la magnétosphère de Saturne les ceintures de radiation de protons de haute énergie (de l'ordre de quelques MeV) s'avèrent être isolées de la magnétosphère moyenne et externe, et la source de ces protons de haute énergie devraient être liée aux rayons cosmiques galactiques (GCR). Pour valider cette hypothèse il est d'abord nécessaire de détenniner le flux de GCR accédant à Saturne de manière réaliste. Auparavant, seulement des tentatives théoriques ont été effectuées afm de vérifier cette idée. Dans cette thèse, pour la première est développée nne solution numérique pour la détermination de l'accès des GCR à l'atmosphère et aux anneaux de Saturne. La méthode proposée est basée sur le traçage de particules chargées et le code a été développé spécifiquement pour la magnétosphère de Saturne. Lors de la validation de la méthode les observations de Cassini MIMI1 LEMMS acquises pendant les survols de Rbéa et de Dioné ont été modélisées à l'aide du traceur et les résultats obtenus ont été comparés aux observations. TI a été découvert que le << draping >> des lignes de champ magnétique autour de ces satellites de glace, même si il produit des perturbations locales de seulement quelques pour cent du champ magnétique ambient, peut produire des changements mesurables dans la distribution spatiale et l'énergie des flux des ions énergétiques mesurés par MIMI1 LEMMS. Ces résultats sont importants pour l'interprétation correcte des données MIMI 1 LEMMS, et offrent des fonctionnalités pour l'étalonnage croisé précis en vol des instruments. Après cette validation du traceur de particules il a été appliqué pour un calcul à rebours dans le temps des GCR accédant à Saturne. L'énergie d'accès des GCR a été obtenue, les spectres des GCR ont été reconstruits et le flux intégré de GCR autour de Saturne et de ses anneaux a été calculé. Les résultats obtenus sont essentiels pour la compréhension de la formation des ceintures de radiation de protons, ainsi que pour la future investigation du processus CRAND sur Saturne, pour l'évaluation de l'intensité de la ceinture de radiation intérieure et pour d'autres projets, discutés dans cette thèse. / The MeV proton radiation belts of Saturn are isolated from the middle and outer magnetosphere and the source of these high energy protons should be related to the access of Galactic Cosmic Rays (GCRs) in the system. To validate this hypothesis it is first of all necessary to determine the realistic spectrum of GCRs at Saturn. Previously only theoretical attempts were performed in order to calculate the GCR spectra. In this thesis I provide for the first time the numerical solution for the determination of the GCR access to the upper atmosphere and rings of Saturn. The proposed method is based on the charged particle tracing technique and a code that was developed specifically for this purpose. For the validation of the code, the Cassini MIMI/LEMMS observations during the Rhea and Dione flybys were modeled using the tracer and the obtained results were compared to the observations. It was demonstrated that even a weak perturbation of the magnetic field lines can produce measurable changes in the spatial and energy distribution of fluxes measured by MIMI/LEMMS that can be accurately simulated by particle tracing. These results are important for the correct interpretation of the MIMI/LEMMS data, and offer capabilities for a precise in-flight instruments' cross-calibration besides the validation of our simulation code. After this validation the particle tracer was applied for simulating the access of the GCRs. The GCRs access to the rings and atmosphere was obtained, the GCRs spectra were reconstructed and were in part also validated using additional Cassini observations. Dependencies of the spectral parameters on the time, incidence direction, etc., were also obtained offering all necessary information for simulating the interaction of GCRs with the Saturnian system during different phases of the Cassini mission. That includes also the Proximal orbits of 2017, during which Cassini will sample for the first time the radiation belts inside the D-ring of the planet, a region which is likely populated only by GCR secondaries.

Au cours de cette thèse, le rôle joué par la turbulence de paroi dans le transport de particules a été étudié. Pour ce faire, deux types d'expériences ont été réalisés. Les premières ont porté sur un écoulement de turbulence de paroi naturelle chargé en billes de verre et de céramique, pour 4 conditions hydrauliques données. Les secondes expériences ont été réalisées en turbulence contrôlée. En effet, nous avons choisi de générer artificiellement des éjections en créant des hairpin vortex. Les premières expériences sont réalisées dans des conditions de fort transport par suspension et les secondes sont faites pour des conditions très proches du seuil de suspension, les particules étant majoritairement transportées par saltation.<br />Afin de mesurer simultanément le champ de vitesse de l'écoulement ainsi que la position et la vitesse des particules, nous avons utilisé la PIV (Vélocimétrie par Image de Particules) qui permet d'avoir sur la même paire d'image l'information liquide et solide. Nous pouvons ainsi analyser l'écoulement autour de chaque particule et mettre en évidence l'action locale du fluide sur les particules.<br />L'analyse de l'écoulement se fait dans le voisinage de chaque particule, dont la taille varie entre 5 à 10d. Tout d'abord, l'analyse par quadrants montre que les éjections sont fortement impliquées dans le transport par suspension et par saltation, notamment dans la montée des particules. Elle indique aussi que certaines particules qui descendent sont dans des éjections. Nous cherchons donc à déterminer ce qui fait qu'une éjection est capable ou non de contrer la gravité et de faire monter les particules. <br />L'étude du flux de quantité de mouvement des éjections montre que plus ce flux est important, plus les particules ont tendance à monter dans les éjections. Nous déterminons donc un seuil permettant de détecter les éjections faisant monter systématiquement les particules. Ce seuil, adimensionalisé , est compris de l'ordre de 0,5 et est quasiment constant en fonction du nombre de Reynolds de l'écoulement, de la hauteur et des caractéristiques des particules. <br />Nous observons qu'une différence importante entre la suspension et la saltation est le temps de résidence des particules dans les éjections. En effet, plus ce temps est long, plus la particule va pouvoir monter dans l'éjection et plus elle sera transportée par suspension. En revanche quand elle ne reste que peu de temps dans l'éjection, la particule n'a pas le temps de monter et elle est alors transportée par saltation. A partir de l'analyse de ces expériences, nous proposons un modèle conceptuel de transport par suspension dans les éjections en fonction du flux de quantité de mouvement autour des particules, de leur position dans les éjections et des trajectoires de suspension qui ont pues être reconstruites pour les expériences en turbulence contrôlée.

L’érosion interne est la cause principale de ruptures des ouvrages hydrauliques en terre, tels que les digues et les barrages. Les conséquences de telles ruptures sont conséquentes et couteuses. L’objectif de cette thèse est de mieux comprendre une catégorie particulière de processus d’érosion, la suffusion, en en caractérisant l’initiation et l’évolution. Cette recherche, à dominante expérimentale, s'appuie sur une modélisation physique unidimensionnelle. Les essais sont réalisés dans une conduite en plexiglas de 18 cm de diamètre intérieur, équipée de capteurs de pression, d’un débitmètre et d’un turbidimètre. Un dispositif de collecte des particules érodées est installé à l’aval de la conduite pour quantifier les particules érodées à des intervalles de temps choisis. Dans cette étude, on a réalisé une série d’essais sur des matériaux pulvérulents. Dans ce type de matériaux, l’érosion se manifeste par suffusion. Elle correspond à l’arrachement et au transport des particules fines à travers l’espace poral des particules grossières. Dans un premier temps, un protocole de reconstitution de matériaux pulvérulents est élaboré. Le principe consiste à réaliser des mélanges humides avec une teneur en eau relative à la masse des particules fines. Pour surmonter les problèmes rencontrés dans la réalisation des essais, la conduite est disposée verticalement, et les matériaux testés sont reconstitués à partir de mélanges de particules sableuses. Une étude paramétrique est menée sur l’influence de certains paramètres, tels que la granulométrie, la pression d’écoulement, la nature et la teneur en particules argileuses ainsi que la densité de mise en place. Les résultats obtenus montrent que l’érodabilité des matériaux dépend de la forme des courbes granulométrique. En effet, l’érosion est plus accentuée pour les matériaux contenant moins de particules fines. Elle est également plus forte quand la courbe granulométrique est discontinue. L’ajout de particules argileuses augmente la résistance à la suffusion des matériaux testés. Trois catégories de matériaux argileux ou fins sont testées : une illite (argile verte de Velay), une kaolinite (Speswhite), et du sable broyé (C10). Les résultats montrent que les matériaux contenant de l’illite sont plus résistants à la suffusion que les matériaux contenant de la kaolinite, tandis que les matériaux contenant du sable broyé sont les moins résistants. Cette résistance à la suffusion augmente avec la teneur en particules argileuses. L’érosion des matériaux dépend aussi de la densité initiale des matériaux testés. En effet, pour la même pression d’entrée, la quantité des particules érodées diminue en augmentant la densité de mise en place / Internal erosion is the main cause of failure of hydraulic structures such dykes and dams. The consequences of such failures are substantial and costly. The objective of this thesis is to better understand one of the erosion phenomena, the suffusion, and to characterize the initiation and the evolution of this phenomenon. A physical modeling approach is used for this study. Tests are carried out in a Plexiglas pipe of 18 cm inner diameter, equipped with pressure sensors, flowmeter and turbidimeter. A device for collecting the eroded particles, at selected time intervals, is installed in the downstream part of the device. In this study we performed a series of tests on coarse cohesionless soil. In this type of material, erosion is manifested by suffusion. It corresponds to the detachment and transport of fine particles through the pore space of the coarse particles. First, a cohesionless soil reconstitution protocol is developed. The idea is to make wet mixtures with water content depending on the fine particle content. To overcome some problems encountered in carrying out the tests, the device is arranged vertically, and the materials tested are reconstituted from mixtures of sand particles. A parametric study on the influence of several parameters such as particle size distribution, flow pressure, nature and content of clay particles and initial density are conducted. The results show that the erodibility of the material depends on the shape of the particle size distribution. Indeed erosion is more pronounced for materials containing less fine particles. This erodability is more pronounced when the curve is gap-graded. Erodibility of the tested soils increases with the hydraulic load. The addition of clay particles increases the resistance to suffusion of the soils. Three categories of clayed or fine soils are tested: illite (Argile Verte de Velay), kaolinite (Speswhite), and crushed sand (C10). The results show that materials containing illite are more resistant to suffusion than materials containing kaolinite, whereas materials containing crushed sand are less resistant. This resistance to suffusion increases with clay particle content. The erosion of materials also depends on the initial density of the soil tested. Indeed, for the same applied pressure, the quantity of eroded particles decreases with increasing the initial density

Le transport radial des particules dans les tokamaks constitue une des questions les plus cruciales pour la communauté de la fusion par con finement magnétique. En eff et, d'une part la puissance de fusion est proportionnelle au carré de la pression, d'autre part l'accumulation d'impuretés lourdes dans le coeur du plasma conduit à d'importantes pertes par rayonnement qui peuvent fi nir par causer un e ffondrement radiatif du plasma. Les dent de scie et la redistribution périodique de la température et de la densité de coeur qui lui est associée peuvent a ffecter signifi cativement le transport radial des électrons et des impuretés. Dans cette thèse, nous présentons des simulations numériques de dents de scie utilisant un code tridimensionnel non linéaire de magnétohydrodynamique appelé XTOR-2F, a n d'étudier le transport de particules pendant les dents de scie. Nous montrons que le code est capable de reproduire les structures fines de densité observées après le crash de la dent de scie avec le diagnostic de réfl ectométrie à balayage rapide sur les tokamaks Tore Supra et JET. La présence de ces structures implique la possibilité que le crash de dent de scie ne soit pas aussi effi cace que prévu pour évacuer les impuretés du coeur du plasma. Cependant, en appliquant le code aux impuretés, nous montrons que finalement le taux de redistribution est quantitativement similaire à ce qui est prévu par le modèle de Kadomtsev, un résultat inattendu a priori. Nous concluons que la dent de scie est e fficace pour évacuer les impuretés du coeur du plasma.

Le but de cette thèse est d'étudier le transport et le dépôt de particules solides dans les écoulements à travers les fractures. Dans un premier temps, l'écoulement monophasique à travers les fractures est étudié afin d'évaluer la validité de la loi cubique locale comme modèle de l'écoulement. Des canaux à parois sinusoïdales à géométrie variable sont utilisés pour représenter différents types de fractures. Un premier développement analytique montre que l'ouverture hydraulique de la fracture diffère de son ouverture moyenne lorsque la rugosité des parois est élevée. La méthode des éléments finis est ensuite utilisée pour résoudre les équations de continuité et de Navier-Stokes et comparer les solutions numériques aux prédictions théoriques de la loi cubique locale sur une gamme relativement étendue de nombres de Reynolds Re. Pour de faibles Re, typiquement inférieurs à 15, la loi cubique locale décrit raisonnablement l'écoulement, surtout lorsque la rugosité et le déphasage entre les parois sont relativement faibles. Dans un deuxième temps, les écoulements chargés de particules sont étudiés. Une approche analytique est d'abord développée pour montrer comment des particules distribuées dans un écoulement stationnaire et laminaire à travers une fracture peuvent être transportées sur de longues distances ou au contraire se déposer à l'intérieur. Plus précisément, une équation simple décrivant la trajectoire d'une particule est établie. Sur la base de cette équation, il est démontré que, quand l'inertie des particules est négligeable, leur comportement dépend directement de la géométrie de la fracture et d'un nombre adimensionnel W qui relie la vitesse de sédimentation des particules à la vitesse moyenne de l'écoulement. L'équation proposée est vérifiée en comparant ses prédictions à des simulations numériques de suivi de particules prenant en compte l'inertie des particules et résolvent complètement les équations de Navier-Stokes. Il est montré que l'équation est valide lorsque l'inertie du fluide est faible. Des diagrammes de régimes, permettant de prévoir le comportement des particules à travers la fracture sont proposés. Enfin, un appareil expérimental conçu dans le but d'effectuer une évaluation pratique du modèle analytique est présenté et les résultats préliminaires sont discutés. Les résultats expérimentaux préliminaires tendent valider le modèle analytique. De façon plus générale, les résultats obtenus à travers ce travail de thèse font progresser nos connaissances du comportement des petites particules transportées dans les écoulements de fractures. Potentiellement, ce travail devrait permettre d'améliorer notre prévision de la pollution souterraine, et peut avoir des applications dans le développement de nouvelles techniques de filtration de l'eau et de séparation des minéraux / The aim of the present thesis is to study the transport and deposition of small solid particles in fracture flows. First, single-phase fracture ow is investigated in order to assess the validity of the local cubic law for modeling ow in corrugated fractures. Channels with sinusoidal walls having different geometrical properties are considered to represent different fracture geometries. It is analytically shown that the hydraulic aperture of the fracture clearly deviates from its mean aperture when the walls roughness is relatively high. The finite element method is then used to solve the continuity and the Navier-Stokes equations and to simulate fracture ow in order to compare with the theoretical predictions of the local cubic law for Reynolds numbers Re in the range 0.067-67. The results show that for low Re, typically less than 15, the local cubic law can properly describe the fracture ow, especially when the fracture walls have small corrugation amplitudes. For Re higher than 15, the local cubic law can still be valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitude, and moderate phase lag between its walls. Second, particle-laden flows are studied. An analytical approach has been developed to show how particles sparsely distributed in steady and laminar fracture flows can be transported for long distances or conversely deposited inside the channel. More precisely, a rather simple particle trajectory equation is established. Based on this equation, it is demonstrated that when particles' inertia is negligible, their behavior is characterized by the fracture geometry and by a dimensionless number W that relates the ratio of the particles sedimentation terminal velocity to the ow mean velocity. The proposed particle trajectory equation is verified by comparing its predictions to particle tracking numerical simulations taking into account particle inertia and resolving the full Navier-Stokes equations. The equation is shown to be valid under the conditions that ow inertial effects are limited. Based on this trajectory equation, regime diagrams that can predict the behavior of particles entering closed channel flows are built. These diagrams enable to forecast if the particles entering the channel will be either deposited or transported till the channel outlet. Finally, an experimental apparatus that was designed to have a practical assessment of the analytical model is presented. Preliminary experimental results tend to verify the analytical model. Overall, the work presented in this thesis give new insights on the behavior of small particles in fracture flows, which may improve our prediction and control of underground contamination, and may have applications in the development of new water filtration and mineral separation techniques

L’algorithme Adaptive Multilevel Splitting (AMS) a récemment fait son apparition dans la littérature de mathématiques appliquées, en tant que méthode de réduction de variance pour la simulation Monte Carlo de chaı̂nes de Markov. Ce travail de thèse se propose d’implémenter cette méthode de réduction de variance adaptative dans le code Monte-Carlo de transport de particules TRIPOLI-4,dédié entre autres aux études de radioprotection et d’instrumentation nucléaire. Caractérisées par de fortes atténuations des rayonnements dans la matière, ces études entrent dans la problématique du traitement d’évènements rares. Outre son implémentation inédite dans ce domaine d’application, deux nouvelles fonctionnalités ont été développées pour l’AMS, testées puis validées. La première est une procédure d’encaissement au vol permettant d’optimiser plusieurs scores en une seule simulation AMS. La seconde est une extension de l’AMS aux processus branchants, courants dans les simulations de radioprotection, par exemple lors du transport couplé de neutrons et des photons induits par ces derniers. L’efficacité et la robustesse de l’AMS dans ce nouveau cadre applicatif ont été démontrées dans des configurations physiquement très sévères (atténuations du flux de particules de plus de 10 ordres de grandeur), mettant ainsi en évidence les avantages prometteurs de l’AMS par rapport aux méthodes de réduction de variance existantes. / The Adaptive Multilevel Splitting algorithm (AMS) has recently been introduced to the field of applied mathematics as a variance reduction scheme for Monte Carlo Markov chains simulation. This Ph.D. work intends to implement this adaptative variance reduction method in the particle transport Monte Carlo code TRIPOLI-4, dedicated among others to radiation shielding and nuclear instrumentation studies. Those studies are characterized by strong radiation attenuation in matter, so that they fall within the scope of rare events analysis. In addition to its unprecedented implementation in the field of particle transport, two new features were developed for the AMS. The first is an on-the-fly scoring procedure, designed to optimize the estimation of multiple scores in a single AMS simulation. The second is an extension of the AMS to branching processes, which are common in radiation shielding simulations. For example, in coupled neutron-photon simulations, the neutrons have to be transported alongside the photons they produce. The efficiency and robustness of AMS in this new framework have been demonstrated in physically challenging configurations (particle flux attenuations larger than 10 orders of magnitude), which highlights the promising advantages of the AMS algorithm over existing variance reduction techniques.

L’objectif de ce travail était de tester l’utilisation des radioisotopes naturels à courte période (234Th, 7Be et 210Pb) comme traceurs des transferts des particules dans le bassin versant de la Garonne-Dordogne-Gironde. Dans un premier temps, la faisabilité de la méthode basée sur le 234Th a été testée dans les rivières. Les résultats ont montré qu’une valeur de 5 mBq l-1 correspondant au seuil en 238UD est nécessaire pour l’application du 234Th en rivière. Dans un deuxième temps, le suivi sur deux années (2006-2007) des radioéléments naturels, 234Th, 7Be et 210Pb, et artificiels, 137Cs et 131I, dans le Lot, la Garonne et l’estuaire de la Gironde, a permis de distinguer les différentes étapes du transport des particules en relation avec l’hydrologie. L’atout de l’utilisation de ces traceurs est d’accéder à de premières estimations du temps de résidence des sédiments et du pourcentage des fractions anciennes et récentes. Dans le Lot et la Garonne, il existe une rétention des particules (2-89 jours) pendant l’étiage et les débits moyens, et de fortes remises en suspension des sédiments anciens à débits élevés. La moyenne du pourcentage des sédiments récents dans les rivières est d’environ 25 %, ce qui montre que les matériaux en suspension (MES) dans les rivières sont dominés par des sédiments anciens. Dans l’estuaire de la Gironde, le temps de résidence varie en fonction de l’hydrologie des rivières : pendant les périodes de forts débits, le temps de transit des particules de la Réole jusqu’à l’océan est seulement de quelque mois (~150 jours) ; alors que pendant les périodes d’étiage, l’âge des sédiments augmente dans la zone de turbidité maximum atteignant jusqu’à 400 jours. / The main objective of this study was to develop short-lived radioisotopes (234Th, 7Be and 210Pb) as tracers of particle transport in the Garonne-Dordogne-Gironde watershed. Firstly, the possibility to use 234Th/238U pair to study particle residence time in the Gironde fluvial-estuarine system was tested. The results show that dissolved 238U should be greater than 5 mBq l-1 to permit this application in fluvial systems. Secondly, the two years monitoring of natural radioisotopes, 234Th, 7Be and 210Pb, and artificial radioisotopes, 137Cs and 131I, in the Garonne and Lot rivers and in the Gironde Estuary, had allowed to observe the dynamic of particle transport, linked to the river flow. Based on this database, preliminary particle residence times and percentages of the new and old sediments were estimated. In the Lot and Garonne rivers, the particle retention ranges between 2-89 days during low and mean river discharges. The percentage of the new sediments is average 25 %, which indicates that the old resuspended sediments are dominated the suspended particulate matters (SPM). In the Gironde estuary, the particle residence time is mainly controlled by river discharge. The particle residence time during low river discharge is about 400 days. When the river discharge is increasing, the freshwaters are predominant in the estuary, the maximum turbidity zone (TMZ) is dispersed and the old sediments transported toward the ocean. During this period the particle residence time is only few months.

L'objet de cette étude est d'utiliser la décomposition en modes propres (POD) pour établir un modèle de la composante déterministe et aléatoire d'un écoulement turbulent présentant une instationnarité à grande échelle pseudo périodique.<br />Une étude expérimentale de l'écoulement bidimensionnel en moyenne en aval d'un cylindre semi-circulaire, par vélocimétrie par image de particules en deux temps (PIV2T) caractérise l'écoulement<br />Une analyse POD du champ de vitesse permet d'extraire les modes spatiaux et de définir un paramètre de phase décrivant l'instationnarité à grande échelle qui régit la partie déterministe. La modélisation de l'évolution temporelle des coefficients associés aux modes s'effectue par des fonctions soit harmoniques pour la partie déterministe, soit stochastiques pour la partie aléatoire.<br />La modélisation est en accord avec les mesures expérimentales des premiers moments statistiques en un point et des fonctions de corrélation spatio-temporelle du champ de vitesse.

Les modèles aux moments angulaires constituent des descriptions intermédiaires entre les modèles cinétiques et les modèles fluides. Dans ce manuscrit, les modèles aux moments angulaires basés sur un principe de minimisation d'entropie sont étudiés pour des applications en physique des plasmas. Ce mémoire se découpe en trois parties. La première est une contribution à la modélisation en physique des plasmas à travers le formalisme des modèles aux moments angulaires. Dans celle-ci, le domaine de validité de ces modèles est étudié en régimes non-collisionels. Il est également montré que les opérateurs de collisions proposés pour le modèle M1 permettent de retrouver des coefficients de transport plasma précis. La deuxième partie de ce document concerne la dérivation de méthodes numériques pour l'étude du transport de particules en temps long. Dans ce cadre, des schémas numériques appropriés pour le modèle M1, préservant l'asymptotique, sont construits et validés numériquement. La troisième partie représente un premier pas significatif vers la modélisation multi-espèces. Ici, le modèle aux moments angulaire M1, construit dans un référentiel mobile, est appliqué à la dynamique des gaz raréfiés. Les propriétés de ce modèle sont détaillées, un schéma numérique est proposé et une validation numérique est menée. / Angular moments models represent alternative descriptions situated in between the kinetic and the fluid models. In this work, angular moments models based on an entropy minimisation principle are considered for plasma physics applications. This manuscript is organised in three parts. The first one is a contribution to plasma physics modelling within the formalism of angular moments models. The validity domain of angular moments models in collisionless regimes is studied. It is also shown that the collisional operators proposed for the M1 angular moments model enable to recover accurate plasma transport coefficients. The second part of this document deals with the derivation of numerical methods for the long timescales particle transport. Appropriate asymptotic-preserving numerical schemes are designed for the M1 angular moments model and numerical validations are performed. The third part represents a first important step toward multi-species modelling. The M1 angular moments model in a moving frame is introduced and applied to rarefied gas dynamics. The model properties are highlighted, a numerical scheme is proposed and a numerical validation is carried out.

Le travail réalisé dans cette thèse introduit le concept de l'utilisation de la méthode Level-Set pour simuler l'interface Glace/Air au cours du temps lors du processus de givrage en vol des aéronefs. Pour cela, un code de givrage tri-dimensionnel multi-blocs et parallélisé a été implémenté au sein du solveur NSMB (Navier-Stokes-Multi-Blocks). Il comprend notamment un module de calcul des trajectoires des gouttelettes par une approche Eulérienne compatible avec l'utilisation de grilles chimères et un module thermodynamique pour le calcul des masses de glace incluant deux modèles différents : un modèle algébrique itératif et un modèle à dérivées partielles. Une attention particulière a été portée sur la vérification du code de givrage implémenté en comparant systématiquement, si possible, les résultats obtenus avec les données expérimentales et numériques existantes dans la littérature. Pour cette raison, le module de déformation de maillage existant dans NSMB a été intégré au code implémenté afin de pouvoir simuler le givrage par une méthode traditionnelle. Enfin, un nouveau principe pour le suivi de l'interface glace/air est introduit via l'utilisation d'une méthode Level-Set. Puisque dans ce travail de thèse nous nous intéressons particulièrement au concept, la méthode Level-Set développée est d'ordre un et est résolue implicitement. On montrera cependant que des résultats valides sont obtenus avec une telle approximation. / This thesis introduces the concept of the Level-Set method for simulating the evolution through time of the ice/air interface during the process of in-flight aircraft icing. For that purpose, a three-dimensionnal multi-block and parallelized icing code have been implemented in the NSMB flow solver (Navier-Stokes-Multi-Blocks). It includes a module for calculating the droplet trajectories by an Eulerian approach compatible with the use of chimera grids and a thermodynamic module to calculate the ice masses including two different models : an iterative algebraic model and a PDE model. Particular attention was paid to the validation of the icing code irnplemented by comparing results with existing experimental and numerical data in the literature. For this reason, the existing mesh deformation algorithm in NSMB was integrated into the code to simulate icing by a traditional method. Finally a new principle to track the ice/air interface is introduced using the Level-Set method. Since we are particularly interested in the concept, the Level-Set method developped is first order and solved implicitly. However it will be shown that valid results are obtained with such an approximation.