Depuis quelques années, les codes de calculs Monte Carlo évoluant qui couplent un code Monte Carlo, pour simuler le transport des neutrons, à un solveur déterministe, qui traite l'évolution des milieux dû à l'irradiation sous le flux neutronique, sont apparus. Ces codes permettent de résoudre les équations de Boltzmann et de Bateman dans des configurations complexes en trois dimensions et de s'affranchir des hypothèses multi-groupes utilisées par les solveurs déterministes. En contrepartie, l'utilisation du code Monte Carlo à chaque pas de temps requiert un temps de calcul prohibitif.Dans ce manuscrit, nous présentons une méthodologie originale évitant la répétition des simulations Monte Carlo coûteuses en temps et en les remplaçant par des perturbations. En effet, les différentes simulations Monte Carlo successives peuvent être vues comme des perturbations des concentrations isotopiques de la première simulation. Dans une première partie, nous présenterons donc cette méthode, ainsi que la méthode de perturbation utilisée: l'échantillonnage corrélé. Dans un second temps, nous mettrons en place un modèle théorique permettant d'étudier les caractéristiques de la méthode des échantillons corrélés afin de comprendre ses effets durant les calculs en évolution. Enfin, dans la troisième partie nous discuterons de l'implémentation de cette méthode dans TRIPOLI-4® en apportant quelques précisions sur le schéma de calcul qui apportera une accélération importante aux calculs en évolution. Nous commencerons par valider et optimiser le schéma de perturbation à travers l'étude de l'évolution d'une cellule de combustible de type REP. Puis cette technique sera utilisée sur un calcul d'un assemblage de type REP en début de cycle. Après avoir validé la méthode avec un calcul de référence, nous montrerons qu'elle peut accélérer les codes Monte Carlo évoluant standard de presque un ordre de grandeur. / For several years, Monte Carlo burnup/depletion codes have appeared, which couple Monte Carlo codes to simulate the neutron transport to deterministic methods, which handle the medium depletion due to the neutron flux. Solving Boltzmann and Bateman equations in such a way allows to track fine 3-dimensional effects and to get rid of multi-group hypotheses done by deterministic solvers. The counterpart is the prohibitive calculation time due to the Monte Carlo solver called at each time step.In this document we present an original methodology to avoid the repetitive and time-expensive Monte Carlo simulations, and to replace them by perturbation calculations: indeed the different burnup steps may be seen as perturbations of the isotopic concentration of an initial Monte Carlo simulation. In a first time we will present this method, and provide details on the perturbative technique used, namely the correlated sampling. In a second time we develop a theoretical model to study the features of the correlated sampling method to understand its effects on depletion calculations. In a third time the implementation of this method in the TRIPOLI-4® code will be discussed, as well as the precise calculation scheme a meme to bring important speed-up of the depletion calculation. We will begin to validate and optimize the perturbed depletion scheme with the calculation of a REP-like fuel cell depletion. Then this technique will be used to calculate the depletion of a REP-like assembly, studied at beginning of its cycle. After having validated the method with a reference calculation we will show that it can speed-up by nearly an order of magnitude standard Monte-Carlo depletion codes.

Cette thèse a pour but de contribuer à la qualification du formulaire de calcul neutronique NARVAL, dédié aux coeurs de réacteurs à plaques. En particulier, l’objectif est de développer des méthodes innovantes permettant d’utiliser les données expérimentales inédites du programme HIPPOCAMPE pour évaluer la précision du profil de puissance calculé. La complexité provient de la localisation de l’instrumentation (chambres à fission placées entre les assemblages) et des hétérogénéités caractéristiques de ce type de coeurs (géométrie à plaques, poisons consommables et de contrôle solides). Pour aborder ce problème deux voies ont été mises en oeuvre : la première voie consiste à « combiner puis extrapoler » les écarts C/E observés afin de déterminer les incertitudes associées aux facteurs de puissance. Nous avons utilisé, pour ce faire, la méthode « P/A », traditionnellement employée dans les REP électrogènes mais jamais appliquée aux coeurs à plaques à ce jour. La deuxième voie passe, en revanche, par la reconstruction d’une nappe de puissance à utiliser comme référence (comparaison calcul/« expérience-reconstruite ») : nous avons focalisé notre travail sur des techniques géostatistiques. Après avoir constaté que les deux méthodes conduisent à des résultats satisfaisants (erreur comparable à l’incertitude expérimentale cible) nous avons continué notre recherche, en explorant les possibles développements et en introduisant en particulier une nouvelle méthode hybride (associant les techniques géostatistiques à la méthode P/A) qui permet d'améliorer ultérieurement la qualification du profil de puissance (écart-type des écarts C/E cohérent avec la constatation expérimentale). / The aim of this doctoral thesis work is to contribute to the experimental validation of a neutron physic code, called NARVAL, devoted to the analysis of slab cores reactors. The primary objective is to develop some innovative methods in order to validate the computed power map starting from the original experimental data, provided by the HIPPOCAMPE campaign. The particular position of the instrumentation (fission chambers located between the assemblies) and the strong heterogeneities, characterising this specific core design (slab geometry, burnable and control neutron absorbers in solid state) represent the main challenge of this work. Two different approaches are investigated : the first one consists in “combining and extrapolating” the observed calculated/experimental results in order to evaluate the uncertainty of power coefficients. Among different solutions, the “P/A” method is chosen : it is usually employed to perform conventional PWR plant analysis and has never been applied before to slab cores. The latter aims to reconstruct a power map that could be used as a direct reference for code validation : in this case the geostatistical techniques are selected. These methods provide satisfactory results as estimated errors are in good agreement with the experimental uncertainty target. Nevertheless, in this work a new hybrid method, applying the geostatistical technics to the P/A scheme, is investigated and developed. The good agreement between the experimental and the estimated validations of the computed power map attests the noteworthy performance of this new method.

Nous nous sommes attaché à étudier les liens entre propriétés cristallographiques et magnétiques de phases intermétalliques R M-B où R est un élément de terres rares, l'yttrium ou le thorium et M un métal de transition 3d. La présence de bore dans l'environnement atomique des sites de métaux de transition, notamment, agit fortement sur les propriétés magnétiques tant à l'échelle macroscopique que microscopique.<br />Nous apportons une contribution à l'étude des propriétés magnétiques particulières du composé YCo4B, composé modèle et de référence de la famille Rn+mCo3n+5mB2n. L'anisotropie magnétocristalline originale de ce composé est caractérisée à la fois par une réorientation de spin et un processus d'aimantation du premier ordre. Elle s'avère également fort sensible à la pression comme l'ont montré nos mesures magnétiques. La substitution du fer au cobalt dans les composés RCo4B est possible et donne lieu à un schéma préférentiel comme l'illustre nos mesures de diffraction des neutrons. En combinant ces résultats à des mesures magnétiques diverses et de spectroscopie Mössbauer, les propriétés magnétiques de chaque site cristallographique de métaux de transition 3d ont pu être particularisées.<br />Nous avons étudié la solution solide Y1-xThxCo4B et montré l'existence d'une substitution préférentielle du thorium à l'yttrium sur un site cristallographique. Ces deux éléments R non magnétiques ont des effets différents sur les propriétés magnétique des composés. Par comparaison à des études magnétiques sous pression, il a été possible de montrer que ces différences sont principalement dues aux différences de configuration électronique entre le thorium et l'yttrium et peu aux évolutions structurales induites par la substitution.<br />De nouveaux composés de la famille structurale Rn+mCo3n+5mB2n ont été mis à jour récemment. Nous présentons des études des conditions de synthèses et des structures magnétiques de quelques phases R3Co13B2 et R5Co19B6. Le composé Pr5Co19B6 présente notamment des anomalies magnétiques interprétables comme une réorientation de spin à basse température.<br />Nous présentons des résultats d'absorption des rayons X réalisés sur plusieurs compositions de la solution solide CeCo12-xFexB6. Dans cette structure, le cérium présente un état de valence intermédiaire qui évolue notablement en fonction de la teneur en fer. <br />Enfin, nous avons précisé les conditions de synthèse de phases métastables Nd2Fe23B3 et Y3Fe62B14, très riches en fer. Ces structures originales présentent de nombreux types d'environnements atomiques du fer ce qui en fait des composés aux propriétés magnétiques intéressantes. Nous avons notamment montré les évolutions inhabituelles de ces propriétés magnétiques entre les états amorphe et cristallisé de ces alliages.

Le fil de traîne de l'araignée est un matériau remarquable, un biopolymère aux exceptionnelles propriétés mécaniques. La soie d'araignée a été au centre d'une intense recherche utilisant une large gamme de techniques expérimentales et de modélisations théoriques. Néanmoins, ses propriétés macroscopiques n'ont pas encore été correctement reliées à sa structure et sa dynamique microscopiques.<br />Dans ce travail de thèse, les techniques de diffusion neutronique ont été pour la première fois utilisées pour l'étude de la soie d'araignée. La forte contribution de l'hydrogène en diffusion des neutrons ainsi que la différence des longueurs de diffusion de l'hydrogène et du deutérium ont permis de porter un regard nouveau sur les propriétés structurales et dynamiques des soies d'araignée. Ainsi les résultats appuient un modèle hiérarchique à trois phases de nanofibrilles composées de domaines cristallins et d'ordre à courte portée, contenues dans une matrice amorphe. Des expériences complémentaires de diffusion du rayonnement synchrotron suggèrent que l'eau absorbée par la matrice amorphe forme une glace amorphe à basses températures. Des expériences de diffraction de neutrons (abréviation anglaise : WANS) ont montré un pic méridional hors réseau -non observé par les expériences en rayons X (abréviation anglaise : WAXS)- attribué à une structure smectique de feuillets beta dans les domaines d'ordre à courte portée. L'échange de H2O contre D2O pour les expériences de diffusion de neutrons aux petits angles (abréviation anglaise : SANS) a permis d'observer la variation de contraste à l'intérieur des nanofibrilles et entre les nanofibrilles et la matrice. La mobilité moléculaire a été sondée par des techniques de diffusion inélastique et quasiélastique des neutrons. Il semble qu'une hiérarchie de phénomènes de relaxations décrive la soie hydratée, tandis que la soie native a un comportement vitreux à température ambiante.

Ce travail concerne la synthèse des composés à anions mixtes (O, S, F) à base de terre rare et leurs propriétés d'absorption des radiations UV-visible. De nouveaux composés appartenant à la classe des (oxy)fluorosulfures de terres rares ont été isolés. Les déterminations structurales ont été réalisées à l'aide des différentes techniques de diffraction (X, neutrons, électrons). La plupart de ces réseaux dérivent soit de la matlockite PbFCl soit du type structural a-La2O3. Une de ces formulations se révèle être une structure modulée incommensurable dérivant de La2O2S par un mécanisme de cisaillement. Parmi ces matériaux, l'origine de l'absorption optique des fluorosulfures a-LnSF a été déterminée à partir d'une étude couplant calcul de structures électroniques et mesure XPS (bande de valence). Les propriétés optiques d'absorption des différents composés ont été corrélées à leurs caractéristiques structurales.

Le recours à l'énergie électronucléaire pour diminuer les émissions anthropiques de CO2 nécessite des avancées technologiques majeures. Les réacteurs nucléaires de IVe génération doivent répondre à plusieurs contraintes, telles qu'une sûreté améliorée, la régénération du combustible et la minimisation de la production de déchets radioactifs. De ce point de vue, l'utilisation du Cycle Thorium en Réacteurs à Sel Fondu semble prometteuse. Cet axe de recherche, étudié dans le passé, avait cependant débouché sur un concept dont les défauts ont empêché la réalisation. Une nouvelle réflexion est menée afin de trouver des solutions et d'aboutir au concept de Thorium Molten Salt Reactor. Le couplage d'un code de transport de neutrons avec un code d'évolution des matériaux permet de simuler le comportement d'un coeur nucléaire, et de suivre son évolution tout au long de sa vie. Par cette méthode, nous avons étudié un large éventail de configurations de réacteurs. Les performances de ces systèmes ont été évaluées, grâce à un jeu de contraintes qu'ils doivent satisfaire au mieux. Ce travail a permis de comprendre bon nombre de phénomènes physiques régissant le comportement de ces réacteurs. Grâce à cette nouvelle compréhension, la recherche de configurations acceptables a pu aboutir à des solutions satisfaisantes, apportant un souffle nouveau dans le domaine des Réacteurs à Sel Fondu.

Dans l'éventualité d'une contribution significative du nucléaire aux besoins énergétiques mondiaux, le problème de la gestion des déchets actuels remet en question la pérennité de la liere REP. En complément des recherches sur le stockage et l'incinération de ces déchets, il convient donc d'envisager des systemes innovants dédiés à une production d'énergie nucléaire durable, aussi sobre, propre et sûre que possible. Nous nous intéressons ici au concept de réacteur à sels fondus, dont le combustible est liquide. Cette particularité autorise un retraitement en ligne pyrochimique, qui peut permettre de s'affranchir de certaines limites neutroniques. Le projet MSBR (Molten Salt Breeder Reactor) d'un réacteur à sels fluorures et à modérateur graphite a ainsi montré à la fin des années 1960 que la surgénération en cycle thorium et en spectre thermique est accessible avec un retraitement en ligne approprié.<br> A l'aide d'outils de simulation développés autour du code Monte Carlo MCNP, nous réévaluons dans un premier temps les performances d'un système de référence inspiré du projet MSBR. L'étude complète de la phase de mise à l'équilibre de ce réacteur d'une puissance de 2500 MWth, démarré en 232Th/233U, nous permet de valider nos choix de référence. L'équilibre obtenu montre une réduction importante des inventaires et des radiotoxicités induites par rapport aux autres cycles possibles. Le retraitement en ligne associé est suffisamment efficace pour rendre le systeme surgénérateur, avec un temps de doublement d'une trentaine d'années à l'équilibre.<br> A partir du système de référence, nous testons ensuite différentes options en termes d'économie de neutrons, de transmutation et de contrôle de la réactivité. Il en ressort que c'est le retraitement en ligne qui apporte le plus de flexibilité à ce système particulièrement bien adapté à la production d'énergie en cycle thorium. L'étude de scenarios de transition à cette filière quantifie les limites d'un éventuel déploiement à partir du parc français actuel, et montre enfin qu'une gestion raisonnée du plutonium disponible serait nécessaire dans tous les cas.

Au cours de ce travail nous avons dosé par activation neutronique aux neutrons thermiques le chlore, le sodium et le manganèse dans une centaine d'échantillons de névé prélevés en Antarctique de l'Est. Les prélèvements ont été effectués entre la surface et 90 cm de profondeur le long d'un axe côte -intérieur de 830 km. Nous exposons en détail les techniques de prélèvement et d'analyse utilisées ainsi que les méthodes que nous avons employées pour évaluer les contaminations et contrôler la validité de nos résultats. Ce point est particulièrement important étant donné le bas niveau des concentrations que nous avons mesurées (de 700 ppb à 40 ppb pour le chlore qui est un des éléments les plus abondants en Antarctique, et de 2 ppb à 2. 10- 2 ppb pour le manganèse). Après avoir discuté la représentativité des échantillons analysés, nous avons étudié l'influence de différents paramètres (distance à la côte, accumulation, altitude) sur les concentrations en éléments traces. Nous avons ensuite déterminé l'origine des trois éléments que nous avons dosés; nous avons alors accordé une place importante au chlore qui semble exister à l'état particulaire mais aussi en phase gazeuse au dessus de l'Antarctique. Nous avons enfin calculé les taux de déposition du chlore, du sodium et du manganèse au-dessus de l'Antarctique à partir de la relation existant entre l'accumulation et les quantités déposées de ces trois éléments. En utilisant ces résultats et la composition chimique des différentes sources d'aérosols en Antarctique nous avons estimé le bilan des retombées d'origines marine et continentale; pour faire ce calcul, nous avons traité séparément la frange côtière et la partie intérieure du continent. Le s valeurs obtenues (90.10 3 t. an- 1 pour l'apport marin dans la zone côtière, 140.10 3 t. an- 1 pour l'apport marin ainsi que pour l'apport continental dans la zone intérieure) sont du même ordre de grandeur que les rares résultats publiés jusqu'à présent. Un appendice est consacré au dosage de Cl, Mn, Na, Mg et K par activation neutronique dans des échantllions de glace alpine.

Etude des micro-contraintes dans les matériaux texturés hétérogènes par diffraction et modèles de comportement L'objectif de ce travail est le développement d'une méthodologie d'analyse des contraintes utilisant des modèles théoriques pour décrire le comportement élasto-plastique des matériaux polycristallins. L'étude vise d'abord l'interprétation de résultats expérimentaux par des modèles de déformation qui décrivent la création des champs de contrainte dans les matériaux polycristallins déformés. Une attention particulière est portée à l'explication des phénomènes physiques à l'origine des contraintes résiduelles et à la prédiction de leur évolution et de leur influence sur les propriétés du matériau. Dans le premier chapitre, la méthode classique, dite des sin2y, d'analyse des contraintes est présentée. Ensuite, la nouvelle méthode d'analyse, méthode de multiréflexions, basée sur les mesures de déformation en utilisant plusieurs réflexions hkl est introduite. Dans cette méthode, tous les pics de diffraction sont analysés simultanément et la distance interréticulaire dhkl est remplacée par une distance équivalente a. Aussi, sont présentées les méthodes de calcul des constantes élastiques radiocristallographiques qui jouent un rôle crucial dans la détermination des contraintes. La détermination de ces constantes est indispensable pour l'interprétation des différents résultats expérimentaux. De nouvelles méthodes de calculs des constantes élastiques radiocristallographiques utilisant le modèle autocohérent ont été développées et testées. Une attention particulière a été portée au calcul par ce nouveau modèle autocohérent dans le cas des couches superficielles (surface libre). Dans ce modèle, le calcul des forces et contraintes normales à la surface est effectué selon le modèle de Reuss et pour les deux autres directions, c'est le modèle auto-cohérent qui est utilisé. Cette méthode de calcul est particulièrement adaptée au cas de la diffraction des rayons X où seulement une couche superficielle du matériau est examinée (généralement de quelques mm d'épaisseur). Dans le chapitre suivant, deux modèles de déformation ont été développés et utilisés pour déterminer l'évolution des contraintes et analyser les propriétés du matériau. Le premier modèle (LF) est basé sur les formulations de Leffers (Lefers 1968) qui ont été reprises et développées par Wierzbanowski (Wierzbanowski 1978, 1982). Le second est le modèle autocohérent (SC) (Hutchinson 1964, Berveiller et Zaoui 1979). Dans ce travail, le calcul est réalisé à partir de l'algotithme développé par Lipinski et Berveiller (Lipinski et Berveiller 1989). Dans cette approche, le tenseur de Green est utilisé pour décrire les interactions entre les grains. Les grains du polycristal sont considérés comme des inclusions ellipsoïdales (en 3D) dans une matrice homogène. Ces deux modèles de déformation elasto-plastique (LW et SC) sont des outils très utiles pour l'étude des propriétés mécaniques des matériaux polycristallins. Ils permettent la prédiction des propriétés macroscopiques du matériau (texture, courbes contrainte-déformation, surfaces d'écoulement plastique, densité des dislocations, état final des contraintes résiduelles, etc.) à partir de ses caractéristiques microsructurales (systèmes de glissement, loi d'écrouissage, texture initiale, état initial des contraintes résiduelles, etc.) (Wierzbanowski 1978). Des résultats typiques: de texture, écrouissage et énergie stockée, obtenus par ces modèles, ont été comparés aux résultats expérimentaux. Le chapitre 3 est consacré principalement à l'explication des origines physiques des contraintes et de la prédiction de leur évolution, ainsi qu'à leur influence sur les propriétés du matériau. Les contraintes internes sont classées en trois types selon l'échelle : contraintes d'ordre I, II ou III. Une attention particulière est portée aux contraintes d'ordre I et II car ce sont les seules qui sont déterminées à partir de la position des pics de diffraction. Les modèles de déformation ont été utilisés pour l'analyse des contraintes à l'échelle des grains (contraintes du second ordre). L'évaluation quantitative de ce type de contraintes ne peut pas être effectuée directement par des mesures mais elle est possible grâce aux modèles. Les matériaux multi-phasés ont été également étudiés. Pour ces matériaux, l'interprétation des données expérimentales est plus complexe que celle du cas des matériaux monophasés en raison de la nécessité de prendre en compte l'interaction entre les phases. C'est pourquoi, une nouvelle méthode adaptée aux matériaux multi-phasés a été développée et appliquée au cas des aciers inoxydables austéno-ferritiques (aciers Duplex). Les paramètres de déformation plastique ( ph c t - scission critique résolu et Η ph - paramètre d'écrouissage) de chacune des phases ont pu être déterminés. Lors de la déformation plastique, l'évolution des contraintes dans les phases et la création de contraintes d'incompatibilité de second ordre, sont observées et l'influence de la texture cristallographique et de l'anisotropie élastique est étudiée. La méthodologie développée et utilisée dans ce travail a, donc, permis de déterminer quantitativement les contraintes du premier et du second ordre, pour chaque phase. Il a été montré qu'une bonne corrélation entre les déformations déterminées expérimentalement et les résultats théoriques, n'est obtenue que si l'influence des contraintes du second ordre est prise en compte. Aussi, le meilleur lissage des courbes expérimentales est obtenu quand les calculs intègrent les constantes d'élasticité anisotropiques et la texture réelle initiale de l'échantillon. Les méthodes de détermination des contraintes du premier et du second ordre, présentées au troisième chapitre, sont employées pour l'étude des contraintes résiduelles dans des alliages écrouis par laminage croisé (Chapitre 4). Le laminage croisé a été retenu pour ajouter une symétrie de la texture cristallographique et, donc, de réduire l'anisotropie de la pièce (comparé au laminage uniaxial). Les résultats sont présentés pour des séries d'éprouvettes en acier et en alliage de cuivre. Dans le cas de l'alliage de cuivre, les résultats montrent de très faibles niveaux de contraintes d'incompatibilité de second ordre qui peuvent être négligées. Par contre, dans le cas de la ferrite, il faut en tenir compte car leur niveau s'avère important. Les oscillations observées sur les courbes des sin2y peuvent être expliquées, dans ce cas, principalement par la présence de contraintes du second ordre. Enfin, au chapitre 5, une nouvelle méthode d'analyse des contraintes utilisant un faisceau de rayons X avec un angle d'incidence faible et constant (méthode de diffraction en incidence rasante GID-sin2y). Cette méthode présente l'avantage d'une profondeur de pénétration des rayons X constante, contrairement à la méthode des sin2y classique qui présente l'inconvénient d'une forte variation de la pénétration avec l'angle y. C'est pour cette raison que la méthode classique des sin2y est mal adaptée pour l'étude des matériaux à forts gradients de contraintes. Moyennant un choix optimisé des angles d'incidence et du type de rayonnement, la nouvelle méthode s'avère efficace pour l'étude des matériaux à forts gradients de contraintes, en permettant des mesures dans différentes couches proches de la surface. L'incertitude des mesures a été évaluée et le rôle de l'absorption, de l'indice de réfraction et des facteurs de Lorentz-polarisation et de diffusion atomique ont été étudiés. A partir de mesures sur des poudres de référence, l'influence de chacun de ces paramètres a été évaluée et prise en compte dans la détermination de la position des pics de diffraction. Les analyses effectuées ont confirmé la faible influence de l'absorption et des facteurs de Lorentz-polarisation et de diffusion atomique sur la contrainte déterminée. Par contre, ils ont révélé un effet important de l'indice de réfraction, en particulier aux petits angles d'incidence. Pour des angles d'incidence a£100, les corrections sont importantes et modifient les résultats des contraintes d'une manière significative (la correction peut atteindre 70 MPa dans le cas de la poudre). Cet effet et, donc, la correction nécessaire décroît quand l'angle d'incidence augmente.

Pour la conception des cœurs de réacteurs de 4ème génération, une précision accrue est requise pour les calculs des différents paramètres neutroniques. Les ressources mémoire et le temps de calcul étant limités, une solution consiste à utiliser des méthodes de raffinement de maillage afin de résoudre l'équation de transport des neutrons. Le flux neutronique, solution de cette équation, dépend de l'énergie, l'angle et l'espace. Les différentes variables sont discrétisées de manière successive. L'énergie avec une approche multigroupe, considérant les différentes grandeurs constantes sur chaque groupe, l'angle par une méthode de collocation, dite approximation Sn. Après discrétisation énergétique et angulaire, un système d'équations hyperboliques couplées ne dépendant plus que de la variable d'espace doit être résolu. Des éléments finis discontinus sont alors utilisés afin de permettre la mise en place de méthodes de raffinement dite hp. La précision de la solution peut alors être améliorée via un raffinement en espace (h-raffinement), consistant à subdiviser une cellule en sous-cellules, ou en ordre (p-raffinement) en augmentant l'ordre de la base de polynômes utilisée.Dans cette thèse, les propriétés de ces méthodes sont analysées et montrent l'importance de la régularité de la solution dans le choix du type de raffinement. Ainsi deux estimateurs d'erreurs permettant de mener le raffinement ont été utilisés. Le premier, suppose des hypothèses de régularité très fortes (solution analytique) alors que le second utilise seulement le fait que la solution est à variations bornées. La comparaison de ces deux estimateurs est faite sur des benchmarks dont on connaît la solution exacte grâce à des méthodes de solutions manufacturées. On peut ainsi analyser le comportement des estimateurs au regard de la régularité de la solution. Grâce à cette étude, une stratégie de raffinement hp utilisant ces deux estimateurs est proposée et comparée à d'autres méthodes rencontrées dans la littérature. L'ensemble des comparaisons est réalisé tant sur des cas simplifiés où l'on connaît la solution exacte que sur des cas réalistes issus de la physique des réacteurs.Ces méthodes adaptatives permettent de réduire considérablement l'empreinte mémoire et le temps de calcul. Afin d'essayer d'améliorer encore ces deux aspects, on propose d'utiliser des maillages différents par groupe d'énergie. En effet, l'allure spatiale du flux étant très dépendante du domaine énergétique, il n'y a a priori aucune raison d'utiliser la même décomposition spatiale. Une telle approche nous oblige à modifier les estimateurs initiaux afin de prendre en compte le couplage entre les différentes énergies. L'étude de ce couplage est réalisé de manière théorique et des solutions numériques sont proposées puis testées.