Ce travail de thèse s'inscrit dans le cadre de recherches sur le retraitement futur du combustible nucléaire usé et notamment sur la gestion poussée des radionucléides à vie longue tels que les actinides mineurs. Il concerne la fabrication de pastilles céramiques de Couvertures Chargées en Actinides Mineurs (CCAM) dédiées à la transmutation pour les réacteurs à neutrons rapides. Les pastilles céramiques utilisées dans le milieu nucléaire sont classiquement fabriquées en utilisant les procédés issus de la métallurgie des poudres. En raison de la pulvérulence des précurseurs poudres utilisés et de la forte radioactivité des actinides mineurs et notamment de l’américium, un procédé « sans poudre » innovant a été proposé. Ce procédé prévoit la fabrication de pastilles céramiques en utilisant des précurseurs oxydes non plus sous forme de poudre mais sous forme de microsphères (Calcined Resin Microspheres Pelletization, CRMP). Le principe de ce procédé consiste à compacter sous forme de pastilles des microsphères d’oxyde puis à fritter le comprimé obtenu.La première partie de cette étude concerne la synthèse et la caractérisation de précurseurs oxyde sous la forme de microsphères sub-millimétriques d’oxyde de lanthanides (simulants les actinides) par le procédé aux résines. Différents lots de microsphères ont pu être synthétisés afin de mieux appréhender l’influence de certains paramètres de synthèse, tel que la température de calcination, le tri en taille de la résine ou encore la teneur en lanthanides dans le cas des oxydes mixtes ont été investigués dans le but de déterminer leur impact sur les propriétés microstructurales et mécaniques de microsphères d’oxyde.L'étude aborde dans un deuxième temps la mise en forme des microsphères d’oxyde à travers une approche à la fois expérimentale et numérique. L’approche numérique utilise la méthode éléments discrets (Discrete Element Method, DEM), bien adaptée pour ces milieux granulaires. Les microsphères d’oxyde, prises individuellement, sont caractérisées mécaniquement notamment à travers la mesure de leur résistance à l'écrasement. Une mise en relation des conditions de synthèse et des tenues mécaniques des microsphères a été entreprise afin de comprendre l'impact de ces paramètres de synthèse sur le comportement mécanique du matériau. La compaction en matrice de lots de ces microsphères sous forme de pastille a été étudiée. En particulier, la compressibilité d'un certain nombre de microsphères a été analysée expérimentalement puis simulée par la DEM en mesurant la densification de la pastille en cru en fonction de la contrainte axiale appliquée et en décrivant l’évolution de sa microstructure. / One option envisioned for the future management of high level nuclear waste is the transmutation of minor actinides into short-lived fission products in sodium fast reactor. This route requires the development of pellet fabrication processes to prepare Minor Actinide Bearing Blanket (MABB) for the transmutation of americium.Currently, those ceramic pellets are produced by powder metallurgy processes involving numerous grinding and milling steps that generate very fine and highly contaminating and irradiating particles. A viable option for reducing the amount of those fine particles would be to develop a dustless process by working on much coarser particles. In this context, this study is concerned with the pelletization of porous and spherical lanthanides oxide precursors (surrogates of actinides). The present work uses both experimental data and numerical simulations to optimize the pelletization step. The final aim is to obtain, after sintering, homogeneous, dense and undistorted ceramic pellets.Firstly, this study concerns the synthesis and characterisation of these oxide microspheres precursors by the Weak Acid Resin process, which consists in loading beads of ion exchange resin with lanthanides cations and mineralizing the metal loaded resin leads into sub-millimetric-sized oxide microspheres. Comprehensive characterisations of the microstructure and mechanical properties of oxide microspheres have been carried out to better understand their behaviour into the matrix when producing pellets.Secondly, the mechanical properties of a single microsphere were investigated in order to better understand its behaviour during compaction steps. They were also analysed using multi-scale simulations based on the Discrete Element Method (DEM), which is well suited for such particulate materials. In a second approach, compaction studies were carried out in a three parts die to characterize the mechanical behaviour under pressure of a large number of oxide precursors. The behaviour of several microspheres in the matrix was finally simulated using DEM in order to describe interactions between microspheres and to have a better understanding of their evolution during pressing.

Les pulsars sont des étoiles à neutron hautement magnétisées en rotation rapide produisant un rayonnement pulsé. Cette thèse est dédiée à leur magnétosphère, c'est à dire la zone proche de l'étoile à neutron, remplie d'un plasma entraîné par la rotation rapide de celle-ci. Il a été montré dès 1969 que la magnétosphère doit avoir des zones très peu denses arborant des champs électriques intenses capables d'accélérer le plasma raréfié de ces régions à des énergies très élevée le long du champ magnétique. La courbure des lignes de champ, couplé avec la rotation d'une particule autour du champ, cause un rayonnement dit de « synchro-courbure ». L'énergie est rayonnée essentiellement en photons gamma (g). Ces photons peuvent ensuite être convertis par interaction quantique photon γ-champ magnétique ou γ-γ en une paire électron-positron e+e- dont chaque composante rayonne à son tour, résultant en une cascade qui alimente la magnétosphère en plasma. Cette thèse traite particulièrement de deux phénomènes clefs de ces cascades : le rayonnement de synchro-courbure et la création de paires par interaction γγ.La théorie quantique du rayonnement de synchro-courbure est développée pour la première fois à partir des principes de base de l'électrodynamique quantique. Les paramètres compatibles avec les approximations du calcul correspondent à une large gamme de conditions physiques typiques des magnétosphères de pulsars. Les transitions quantiques sont considérées dans l'approximation continue lorsqu'elles impliquent un saut de l'impulsion de la particule dans la direction parallèle au champ, et discrète dans la direction perpendiculaire. Il en résulte un spectre tendant asymptotiquement vers les descriptions classiques des rayonnement de courbure et de synchro-courbure mais présentant des déviations très importantes lorsque les transitions discrètes dominent le rayonnement.L’interaction γγ→e+e- a été étudiée dans le cas où un gamma réagit sur un fond de photons de basse énergie. Ce mécanisme est considéré comme potentiellement important lorsque le champ magnétique n'est pas assez fort pour produire des paires par le mécanisme γ-champ magnétique. Tout indique que le fond est anisotrope, c'est pourquoi nous avons développé un formalisme permettant de prendre en compte arbitrairement les anisotropies et de produire les spectres des particules produites. Appliqué à un modèle simple d'étoile rayonnant thermiquement en X, il en résulte une dépendance forte du taux de réaction sur la direction du photon gamma.Cette thèse comprend également un modèle de chronométrage du pulsar milliseconde dans un système triple J0337+1715. Ce pulsar orbite avec deux étoiles naines blanches dont les interactions mutuelles ne sont pas négligeables. Une intégration numérique, à l'ordre newtonien et post-newtonien, a été développée pour déterminer les orbites. Un modèle complet incluant le calcul des retards du système du pulsars au télescope a été réalisé. Le modèle s'ajuste aux données de chronométrage provenant du radiotélescope de Nançay avec des résidus d'écart-type inférieur à 2 µs. Un tel système permet en principe le test du principe d'équivalence fort gravitationnel par une technique similaire à celle employée lors des expériences de laser-lune, mais avec une précision sans précédent en régime de champ fort. Ce test requiert une évaluation rigoureuse des incertitudes sur chaque paramètre, échantillonnées grâce à un code MCMC. La validation du code et l'évaluation des incertitudes sont en cours. / Pulsars are highly magnetized fast rotating neutron stars producing a pulsed radiation. This thesis is dedicated to their magnetosphere, namely the zone surrounding the star and filled with a plasma dragged by the rotation of the star. It was shown as soon as 1969 that the magnetosphere must have vacuum gaps, where intense electric fields develop that are capable of accelerating the rarefied plasma to very high energies along the magnetic field. The curvature of the field lines, together with the rotation around the magnetic field, results in the so-called «  synchrocurvature » radiation. The energy is mostly radiated in gamma photons (γ). These photons may then be converted by the quantum processes γ photon-magnetic field or γ-γ in an electron-positron pair e+e-, each component of which then radiates at its turn which results in a cascade that provides plasma to the magnetosphere. This thesis particularly deals with two key phenomena of these cascades : synchrocurvature radiation and γγ pairs.The quantum theory of synchrocurvature radiation is developed for the first time from the first principles of quantum electrodynamics. The range of parameters compatible with the approximations of the derivation covers a wide range of physical conditions typical of pulsar magnetospheres. Quantum transitions are considered in the continuous limit when they imply a jump of the particle impulsion parallel to the magnetic field, and discrete when the jump is in the perpendicular direction. It results in a spectrum that asymptotically tends to the classical descriptions of curvature and synchrocurvature radiations but that presents very important deviations when the discrete transitions dominate the radiation.The γγ→e+e- process was studied in the case of the reaction of a gamma photon on a soft photon background. This mechanism is considered as potentially important when the magnetic field is nopt strong enough for the γ-magnetic field process to efficiently produce pairs. The soft background is most likely anisotropic, and that is why we developed a formalism allowing to arbitrarily take into account anisotropies, as well as produce the spectra of the outgoing particles so as to be able to feed the subsequent cascade consistently. Applied to a simple model of a star radiating thermal X rays, it results in a strong dependence of the reaction rate on the direction of the gamma photon.This thesis also includes a timing model of the millisecond pulsar in a triple system J0337+1715. This pulsar orbits with two white-dwarf stars, and their mutual interactions are not negligible. It follows that a numerical integration of the orbits was developed at Newtonian and first post-Newtonian orders. A complete model including the computation of delays from the star to the telescope was realized. This model is able to fit the timing data from the Nançay (France) radiotelecope with a standard deviation of less than 2µs. In principle, such a system allows to test the strong equivalence principle by a technique similar to that employed in Lunar-laser-ranging experiments, but with an unprecedented accuracy in the strong-field regime. This test demands a careful estimate of the uncertainties on each parameter, which we sample using a MCMC code. The validation of the code and the evaluation of the uncertainties are ongoing.

Cette thèse présente les résultats de deux études : la première concerne la convection naturelle turbulente dans une cavité rectangulaire chauffée uniformément par le bas et remplie d’un nanofluide et la seconde concerne l’investigation du transfert de chaleur conjugué dans un échangeur de chaleur à tubes ailetés.L’enceinte de la première étude a un faible rapport d’aspect. Ses parois gauche, droite et supérieure sont maintenues à une température relativement basse. Le fluide de travail est un nanofluide constitué d’eau et de nanoparticules, soient d’alumine (Al2O3), ou de cuivre (Cu) ou d’oxyde de cuivre (CuO). L’influence des paramètres tels que le nombre de Rayleigh (basé sur la hauteur H de la cavité et la densité de flux de chaleur), le type de nanofluide et la fraction volumique des nanoparticules sur la performance de refroidissement est présentée. Les équations de Navier-Stokes et les équations de conservation de la masse et de l'énergie sont résolues pour une géométrie bidimensionnelle par la méthode numérique des volumes finis. L'algorithme SIMPLE est utilisé pour le couplage pression-vitesse. La discrétisation des termes convectifs est faite avec le schéma QUICK. Le modèle de turbulence k-epsilon standard est utilisé. Le maillage du domaine simulé est généré par le code Gambit. Les résultats montrent que pour toutes les valeurs de Ra, le nombre de Nusselt moyen augmente d’une façon linéaire et monotone avec l’augmentation de la concentration des nanoparticules dans le fluide de base. Le flux de chaleur moyen prend des valeurs qui diminuent en fonction de l’ordre suivant : Cu, CuO et Al2O3.La deuxième étude est une investigation numérique des caractéristiques dynamique et thermique d'un échangeur de chaleur. Les calculs supposent un transfert de chaleur et un écoulement en régime permanent. Le nombre de Nusselt et le coefficient de frottement qui caractérisent l'échangeur de chaleur sont déterminés pour différentes valeurs du nombre de Reynolds. L’équation de conservation de l'énergie dans le fluide et l’équation de conduction de la chaleur dans le solide en trois dimensions ont été résolues avec les équations de la conservation de la masse et de la quantité de mouvement afin de déterminer ces caractéristiques. Les deux régimes d'écoulement laminaire et turbulent sont considérés. L'effet de la modélisation de la turbulence a été étudié en utilisant trois modèles différents (Spalart-Allmaras modèle de turbulence à une équation, le modèle k-epsilon ; standard et le modèle RSM). La validation du modèle a été effectuée en comparant les facteurs de frottement, f, et le facteur, j, de Colburn avec les données expérimentales trouvées dans la littérature. Les résultats tracés ont montré un bon accord qualitatif entre les résultats numériques et les données expérimentales. Les résultats montrent également que le plus simple des trois modèles de turbulence testés (à savoir. Spalart-Allmaras) donne les valeurs les plus proches des données expérimentales. / This thesis presents the results of two studies: the first concerns natural turbulent convection in a rectangular cavity heated from the bottom wall and filled with a nanofluid and the second relates to the investigation of conjugate heat transfer in a fin-and-tube heat exchanger.The cavity of the first study is tall and has a heat source embedded on its bottom wall, while its left, right and top walls are maintained at a relatively low temperature. The working fluid is a water based nanofluid having three nanoparticle types: alumina, copper and copper oxide. The influence of pertinent parameters such as the Rayleigh number, the type of nanofluid and solid volume fraction of nanoparticles on the cooling performance is studied. Steady forms of twodimensional Reynolds-Averaged-Navier-Stokes equations and conservation equations of mass and energy, coupled with the Boussinesq approximation, are solved by the volume control based on the discretisation method employing the SIMPLE algorithm for pressure-velocity coupling. Turbulence is modeled by using the standard k-epsilon model. The Rayleigh number, Ra, is varied from 2.49xE09 to 2.49xE11. The volume fractions of nanoparticles where varied too. Stream lines, isotherms, velocity profiles and temperature profiles are presented for various combinations of Ra, the type of nanofluid and solid volume fraction of nanoparticles. The results are reported in the form of average Nusselt number on the heated wall. It is shown that for all values of Ra, the average heat transfer rate from the heat source increases almost linearly and monotonically as the solid volume fraction increases. Finally the average heat transfer rate takes on values that decrease according to the ordering Cu, CuO and Al2O3.In the second study We determined the heat transfer and friction characteristics of a realistic fin-and-tube heat exchanger. The computations assume steady-state heat transfer and fluid flow. Nusselt number and friction factor characteristics of the heat exchanger are presented for various values of Reynolds numbers. The energy conservation and the heat conduction equations in 3 dimensions have been solved in the fluid and the solid respectivelyalong with the mass and momentum conservation equations in order to determine these characteristics. Both laminar and turbulent flow regimes are considered. The effect of turbulence modeling was investigated using three different models (the one equation Spalart-Allmaras turbulence model, the standard k-epsilon; model and the RSM model). The computations allowed the determination of the dynamic and thermal fields. Model validation was carried out by comparing the calculated friction factor f and Colburn j-factor to experimental results found in the literature. The plotted results showed a qualitatively good agreement between numerical results and experimental data. The results obtained also showed that the simplest of the three turbulence models tested(i.e. Spalart-Allmaras) gives the closest values to the experimental data.

La fabrication des pièces composites structurales aéronautiques nécessite l'utilisation de renforts complexes et notamment des renforts tridimensionnels. Le procédé d'infusion permettrait la fabrication de ces pièces tout en maîtrisant leur qualité et les coûts. L'objectif de la thèse est de développer un modèle, à différentes échelles, de l'écoulement de résine à travers le réseau fibreux du renfort et de proposer une simulation numérique plus fine du procédé d'infusion. L'influence des structures de renfort 2D et 3D, sera ainsi étudié aux différentes échelles : micro (intra-mèches), méso (inter-mèches)et macro(renfort, pièce). / The manufacture of composite aeronautical structural parts requires the use of more complex reinforcements, especially the 3D reinforcements. Resin infusion is the key process of Liquid Composites Molding (LCM) which will determine the quality and cost of the parts. The target of this thesis is to develop models at different scales of fibrous network of the reinforcement and to propose a finer numerical simulation of the infusion process. The influence of the 2D and 3D reinforcement structures on the infusion process will be studied at the different scales: micro (intra-tow), meso (inter-tows) and macro (reinforcement, part).

Ces travaux de thèse s’inscrivent dans un thème de recherche portant sur l’optimisation de la phase de conception et la préparation de la fabrication par estampage de pièces de formes complexes. Le préformage des renforts tissés secs est un enjeu important pour plusieurs procédés de production de pièces en matériaux composites tel que, par exemple, le procédé RTM (Resin Transfer Molding). Au cours de cette phase, la préforme est soumise à des déformations importantes. La connaissance du comportement du tissu sec est alors un enjeu majeur en vue de l’optimisation des procédés de mise en forme. Pour mettre en œuvre les renforts tissés, il est nécessaire de tenir compte de leurs caractéristiques intrinsèques aux différentes échelles, de leurs très grandes déformations en cisaillement et du comportement fortement orienté de ces matériaux. Un point important réside dans la détermination des orientations des renforts après formage. Face à la complexité de mise au point expérimentale de la mise en forme des renforts tissés, la simulation est un outil important pour l’optimisation de conception de pièces composites. Dans ce travail, une nouvelle approche hybride discrète non linéaire, basée sur l’association d’éléments continus hypoélastiques (comportement en cisaillement non linéaire) avec des connecteurs spécifiques de comportement non linéaire a été abordée. Elle permet de prédire les contraintes au niveau des fibres et de déterminer avec plus de précision, les angles de cisaillement en se basant sur la modification de l’orientation en grande déformation. En outre, elle permet d’analyser et de prévoir le comportement global du tissu à partir de sa structure interne. Le nombre de paramètres à identifier est faible et le temps de calcul est raisonnable. Cette approche a été programmée via une routine VUMAT et implémentée dans le code de calcul élément fini ABAQUS/Explicit. L’identification et la validation du modèle ont été effectuées en utilisant des essais de caractérisation standard des tissus. Les résultats de mise en forme des renforts tissés ont été comparés à des résultats expérimentaux.

Le ballottement dans les réservoirs d'ergols est une des perturbations les plus importantes de la stabilité d'un satellite en orbite. En considérant des manœuvres faiblement inertielles, il n'existe pas de modèle analytique et l'expérimentation nécessite de longues périodes de temps en micro-gravité. Nous nous proposons donc, dans cette thèse, de réaliser des simulations numériques de ce phénomène. L'étude est basée sur le solveur DIVA résolvant les équations de Navier-Stokes diphasiques avec les méthodes level-set et Ghost Fluid. Deux technologies de réservoirs sont à l'étude : les réservoirs classiques, ne contenant que l'ergol liquide et le gaz pressurisant, et les réservoirs à membrane, pour lesquels une membrane hyperélastique sépare les deux fluides. Dans le premier cas, une étude paramétrique complète sur les effets du ballottement lors d'une manœuvre de rotation est menée et les différents régimes d'écoulement obtenus sont décrits par rapport aux paramètres d'études. Par la suite, les données de l'expérience FLUIDICS, envoyée à bord de la Station Spatiale Internationale, sont comparées aux résultats numériques et montrent un très bon accord. Par rapport au second cas, un modèle d'interaction fluides-membrane est développé en s'inspirant de travaux sur la déformation de cellules biologiques. Les déformations et contraintes propres à la membrane sont suivies de façon Eulérienne, les efforts exercés par la membrane sur les fluides environnants en sont déduits et intégrés au solveur diphasique. Les résultats obtenus sont validés par comparaison à des cas-tests de la littérature. / Propellant sloshing in tanks is one of the most important disturbances of satellite stability in orbit. Considering low-inertial manoeuvres, there is no analytical model and experimental facilities require long time period of microgravity conditions. Thus, this PhD thesis aims to predict this phenomenon by numerical simulations. The study is based on the DIVA code which solves the Navier-Stokes equations for two-phase flows with the level-set method and the Ghost Fluid method. Two tank technologies are studied: simple tanks, which only contain the liquid propellant and the gas maintaining the pressure, and diaphragm tanks, for which a hyperelastic membrane separates both fluids. In the first case, a parametric study on the sloshing effects is done considering rotational manoeuvres and the different behaviours observed are described in relation to the study parameters. Thereafter, the data from the FLUIDICS experiment, sent to the International Space Station, are compared to the numerical results and exhibit good agreement. In the second case, a fluids-membrane interaction model inspired from works on the deformation of biological cells is developed. The membrane strains and stresses are computed in an Eulerian way, from which the force exerted by the membrane on the surrounding fluids is deduced and integrated in the two-phase flows solver. The numerical results are validated by comparison with benchmarks from the literature.

L'équilibrage de charge est une étape importante conditionnant les performances des applications parallèles. Dans le cas où la charge varie au cours de la simulation, il est important de redistribuer régulièrement la charge entre les différents processeurs. Dans ce contexte, il peut s'avérer pertinent d'adapter le nombre de processeurs au cours d'une simulation afin d'obtenir une meilleure efficacité, ou de continuer l'exécution quand toute la mémoire des ressources courantes est utilisée. Contrairement au cas où le nombre de processeurs ne varie pas, le rééquilibrage dynamique avec un nombre variable de processeurs est un problème peu étudié que nous abordons ici.Cette thèse propose différentes méthodes basées sur le repartitionnement de graphe pour rééquilibrer la charge tout en changeant le nombre de processeurs. Nous appelons ce problème « repartitionnement M x N ». Ces méthodes se décomposent en deux grandes étapes. Dans un premier temps, nous étudions la phase de migration et nous construisons une « bonne » matrice de migration minimisant plusieurs critères objectifs comme le volume total de migration et le nombre total de messages échangés. Puis, dans un second temps, nous utilisons des heuristiques de partitionnement de graphe pour calculer une nouvelle distribution optimisant la migration en s'appuyant sur les résultats de l'étape précédente. En outre, nous proposons un algorithme de partitionnement k-aire direct permettant d'améliorer le partitionnement biaisé. Finalement, nous validons cette thèse par une étude expérimentale en comparant nos méthodes aux partitionneursactuels. / Load balancing is an important step conditioning the performance of parallel programs. If the workload varies drastically during the simulation, the load must be redistributed regularly among the processors. Dynamic load balancing is a well studied subject but most studies are limited to an initially fixed number of processors. Adjusting the number of processors at runtime allows to preserve the parallel code efficiency or to keep running the simulation when the memory of the current resources is exceeded.In this thesis, we propose some methods based on graph repartitioning in order to rebalance the load while changing the number of processors. We call this problem \M x N repartitioning". These methods are split in two main steps. Firstly, we study the migration phase and we build a \good" migration matrix minimizing several metrics like the migration volume or the number of exchanged messages. Secondly, we use graph partitioning heuristics to compute a new distribution optimizing the migration according to the previous step results. Besides, we propose a direct k-way partitioning algorithm that allows us to improve our biased partitioning. Finally, an experimental study validates our algorithms against state-of-the-art partitioning tools.

La pollution d’aquifères par des hydrocarbures est très persistante et difficile à traiter, devenant un enjeu majeur pour l’environnement en raison de l’effet négatif sur la santé humaine. La génération de mousse in situ combinée au lavage de sols est une technique de dépollution innovante, permettant un meilleur contrôle de la mobilité des fluides dans des formations hétérogènes. Le but de cette méthode est de bloquer l’écoulement dans les couches à forte perméabilité pour améliorer l’efficacité de balayage des strates peu perméables.Le mémoire se divise en deux parties principales décrivant une étude expérimentale et une étude numérique.Une sélection préalable de tensioactifs respectueux de l’environnement à été réalisée sur la base de leur moussabilité. Le sucrose laurate (émulsifiant alimentaire) a été retenu et comparé avec un tensioactif dit classique. Des tests sur colonnes 1-D ont permis d’étudier l’influence de la perméabilité, du type de tensioactif et de la qualité de la mousse sur le facteur de résistance à la mobilité du gaz. Les résultats, interprétés à l’aide de la pression capillaire, révèlent deux régimes d’écoulement, correspondants aux deux textures de mousse : « weak » et « strong ». Des expériences sur pilotes 2-D hétérogènes ont permis de visualiser les différentes phases du processus de remédiation et de comparer l’efficacité de balayage avec et sans une région bloquée par de la mousse.La partie modélisation comprend une simulation numérique des tests sur colonne 1-D en utilisant le code UTCHEM et une étude de sensibilité des paramètres les plus importants du modèle, ainsi qu’une modélisation des tests de traçage sur pilote 2-D avec le code MODFLOW. / Aquifer pollution with hydrocarbons is very persistent and difficult to treat, becoming a major issue for the environment because of the negative effect on human health. In situ foam generation combined with soil washing is an innovative remediation technology that allows a better fluid mobility control in heterogeneous formations. The purpose of this method is to block the flow through the high permeability layers in order to improve the sweep efficiency of low permeability strata.The thesis is divided into two main parts describing an experimental study and a numerical study.A preliminary selection of environmentally-friendly surfactants was carried out on the basis of their foamability. Sucrose laurate (emulsifier from food industry) was selected and compared with a conventional surfactant. 1-D column tests were performed to study the influence of permeability, surfactant type and foam quality on the resistance factor. The results, interpreted using the capillary pressure, show two flow regimes, corresponding to the two foam textures: "weak" and "strong". Experiments on 2-D heterogeneous pilot allowed visualizing the different steps of the remediation process and comparing the sweep efficiency with and without a region blocked by foam.The numerical part includes modeling 1-D column tests using UTCHEM, a sensitivity study of the most important parameters of the model, and a simulation of tracer tests on 2-D pilot using MODFLOW.

La propagation acoustique en milieu externe est fortement influencée par l'environnement. Les effets liés à la géométrie, comme la topographie ou la présence d'obstacles, sont principalement les réflexions et les diffractions. Concernant l'effet de l'atmosphère, les gradients moyens génèrent des réfractions tandis que la turbulence provoque des fluctuations aléatoires et une perte de cohérence du signal. La plupart de ces effets sont généralement bien décrits de manière théorique, cependant dans les configurations réelles le cumul de tous ces effets rend l'utilisation des expressions analytiques très difficile. Les études expérimentales présentent également des limites liées à la difficulté de connaître l'environnement parfaitement et d'isoler un effet physique particulier. Dans cette perspective, la simulation numérique est une alternative pratique et complémentaire à la théorie et l'expérimentation. Parmi les modèles numériques de propagation existants, ceux basés sur une résolution par différences finies dans le domaine temporel (FDTD pour Finite-Difference Time-Domain) des équations d'Euler linéarisées sont récents et particulièrement prometteurs. Cependant comme pour tout modèle nouveau, il reste à montrer qu'effectivement l'ensemble des phénomènes physiques d'intérêt sont retranscrits.Dans le cadre de ses études sur la propagation acoustique extérieure, l’Institut franco-allemand de recherches de Saint-Louis (ISL) a implémenté un tel modèle de propagation. Cette implémentation est ci-après appelée ITM, pour ISL FDTD Model. L'objectif de cette thèse, proposée par l'ISL en collaboration avec le Laboratoire de Mécanique des Fluides et d'Acoustique (LMFA), est de poursuivre le développement et les validations de cette implémentation. Une part importante du travail consiste également à illustrer les potentialités du code ITM pour des applications de propagation de signaux acoustiques complexes dans un environnement complexe. [...] / Outdoor sound propagation is strongly influenced by the environment. The geometry, such as topography and the presence of obstacles, alters the sound through reflexions and diffractions. Regarding atmosphere-related effects, the mean gradients produce refractions while turbulence cause random fluctuations and signal coherence loss. Most of those effects are generally well described theoretically. Still, in real configurations, the accumulation of those effects makes the use of analytical expressions difficult. Experimental studies are also limited because of difficulties in perfectly determining the environment or in separating a precise physical effect. In that perspective, numerical simulation is a convenient and complementary alternative approach to theory and experimentation. Among the existing numerical propagation models, those based on a Finite-Difference resolution in the Time-Domain (FDTD) of the linearized Euler equations are recent and particularly promising. However as for every new model, it remains to show that indeed the physical phenomena of interest are reproduced. In the framework of its studies on outdoor sound propagation, the french-german research Institute of Saint-Louis (ISL) has implemented such a propagation model. This implementation is hereafter called ITM, for ISL FDTD Model. The objectives of the thesis, proposed by ISL in collaboration with the Laboratory of Fluid Mechanics and Acoustics (LMFA), are to pursue the developments and validations of this implementation. An important part of the work is also given on the illustration of the potentialitiesof the ITM code in propagating complex acoustic signals in complex environments. […]

In the present work, CO2 corrosion of API 5L X65 pipeline steel in the presence of acetic acid (HAc) was investigated. The objective of this research was to study a qualitative behaviour frequently found in the case of the so-called top-of-line corrosion (TLC). That is, increasing corrosion rates and steep shrinkage after certain - not predictable - time delay. The focus was on the role of HAc in the growth and stifling of localized attacks in CO2-containing media. The kinetic behaviours of carbon steel in such brines indicated that the overall effect of HAc is a balance between an enhanced cathodic reaction rate and a slightly inhibited anodic reaction rate. Results also showed that the active dissolution is not directly related to the presence of acetic acid but rather to what is generically referred to as CO2 corrosion. Zero resistance ammeter (ZRA) measurements on artificial pit electrode assembly, coupled to numerical simulation results, indicated that local HAc concentration gradients sustain the pit growth mechanism to a certain critical depth. Beyond this critical depth, the coupling current sharply drops off and the pit undergoes a process of stifling. However, a substantial effect of the purely ohmic drop on the pit stabilization process appeared unlikely in such conditions. Instead, initiated pits were shown to propagate and stifle mainly according to the counteracting depletions of HAc and CO2. On the basis of the obtained results, a complete scenario was therefore proposed for the morphological trend of localized attacks, which is in good agreement with field occurring TLC cases.