Les surfaces de subdivision fournissent une méthode alternative prometteuse dans la modélisation géométrique, et ont des avantages sur la représentation classique de trimmed-NURBS, en particulier dans la modélisation de surfaces lisses par morceaux. Dans ce mémoire, nous considérons le problème des opérations géométriques sur les surfaces de subdivision, avec l'exigence stricte de forme topologique correcte. Puisque ce problème peut être mal conditionné, nous proposons une approche pour la gestion de l'incertitude qui existe dans le calcul géométrique. Nous exigeons l'exactitude des informations topologiques lorsque l'on considère la nature de robustesse du problème des opérations géométriques sur les modèles de solides, et il devient clair que le problème peut être mal conditionné en présence de l'incertitude qui est omniprésente dans les données. Nous proposons donc une approche interactive de gestion de l'incertitude des opérations géométriques, dans le cadre d'un calcul basé sur la norme IEEE arithmétique et la modélisation en surfaces de subdivision. Un algorithme pour le problème planar-cut est alors présenté qui a comme but de satisfaire à l'exigence topologique mentionnée ci-dessus. / Subdivision surfaces are a promising alternative method for geometric modelling, and have some important advantages over the classical representation of trimmed-NURBS, especially in modelling piecewise smooth surfaces. In this thesis, we consider the problem of geometric operations on subdivision surfaces with the strict requirement of correct topological form, and since this problem may be ill-conditioned, we propose an approach for managing uncertainty that exists inherently in geometric computation. We take into account the requirement of the correctness of topological information when considering the nature of robustness for the problem of geometric operations on solid models, and it becomes clear that the problem may be ill-conditioned in the presence of uncertainty that is ubiquitous in the data. Starting from this point, we propose an interactive approach of managing uncertainty of geometric operations, in the context of computation using the standard IEEE arithmetic and modelling using a subdivision-surface representation. An algorithm for the planar-cut problem is then presented, which has as its goal the satisfaction of the topological requirement mentioned above.

Cette thèse est consacrée aux écoulements homogènes de bulles, ainsi qu'à leur couplage avec le transport d'un scalaire et la turbulence. Elle s'intéresse plus spécifiquement aux effets de taille finie, des interactions hydrodynamiques et de la microstructure de la suspension qui sont étudiés à l'aide de simulations numériques directes à l'échelle d'une seule bulle. La dynamique d'une suspension laminaire de bulles induite par la seule gravité est d'abord revisitée. L'influence de la fraction volumique sur la vitesse de dérive des bulles est établie analytiquement et numériquement pour une suspension parfaitement ordonnée, puis des ressemblances entre suspensions ordonnées et suspensions désordonnées sont mises en évidence. Ces résultats sont ensuite mis à profit pour la modélisation du transport d'un scalaire passif au sein d'une suspension laminaire, tel que décrit par une diffusivité effective tensorielle, et des différences essentielles entre systèmes ordonnés et systèmes désordonnés concernant le transport de scalaire sont mises en exergue. Enfin, la turbulence est prise en compte dans les simulations et son interaction avec une bulle de taille finie est caractérisée. Il est montré que le comportement dynamique d'une bulle de taille comparable à la microéchelle de Taylor ressemble qualitativement à celui d'une microbulle, avec, notamment, une préférence pour certaines régions caractéristiques de l'écoulement. Une définition de l'écoulement vu par la bulle compatible avec les modèles standards de masse ajoutée et de portance est finalement proposée / This thesis is devoted to the study of homogeneous bubbly flows and their coupling with scalar transport and turbulence. It focuses on the effects of finite size, hydrodynamic interactions, and suspension microstructure, which are investigated using direct numerical simulations at the bubble scale. The dynamics of laminar buoyancy-driven bubbly suspensions is first revisited. More specifically, the effect of volume fraction on the bubble drift velocity is clarified by connecting numerical results to theory for dilute ordered systems, and similarities between perfectly ordered and free disordered suspensions are evidenced. These results are then used for the modeling of passive scalar transport in laminar suspensions as described by an effective diffusivity tensor, and crucial differences between ordered and disordered systems with respect to scalar transport are highlighted. Lastly, turbulence is included in the simulations, and its interaction with a finite-size bubble is characterized. The behavior of a bubble as large as Taylor microscale is shown to share a number of common features with that of a microbubble, most notably, the flow sampled by the bubble is biased. A definition of the liquid flow seen by the bubble, as it enters in usual models for the added mass and the lift forces, is finally proposed

Ce travail de thèse porte sur l’étude de l’accélération de particules fluides et inertielles en déplacement dans une turbulence soumise à un gradient de vitesse moyen. L’objectif est de récupérer des données de référence afin de développer des modèles LES stochastiques pour la prédiction de l’accélération de sous-maille et l’accélération de particules inertielles dans des conditions inhomogènes. La modélisation de l’accélération de sous-maille est effectuée à l’aide de l’approche LES-SSAM introduite par Sabel’nikov, Chtab et Gorokhovski[EPJB 80:177]. L’accélération est modélisée à l’aide de deux modèles stochastiques indépendants : un processus log-normal d’Ornstein-Uhlenbeck pour la norme d’accélération et un processus stochastique Ornstein-Uhlenbeck basé sur le calcul de Stratonovich pour les composantes du vecteur d’orientation de l’accélération. L’approche est utilisée pour la simulation de particules fluides et inertielles dans le cas d’une turbulence homogène isotrope et dans un cisaillement homogène. Les résultats montrent une amélioration des statistiques à petites échelles par rapport aux LES classiques. La modélisation de l’accélération des particules inertielles dans le cisaillement homogène est effectuée avec l’approche LES-STRIP introduite par Gorokhovski et Zamansky[PRF 3:034602] et est modélisée avec deux modèles stochastiques indépendants de manière similaire à l’accélération de sous-maille. Nos calculs montrent une amélioration de l’accélération et de la vitesse des particules lorsque le modèle STRIP est utilisé. Enfin dans une dernière partie, nous présentons une équation pour décrire la dynamique de particules ponctuelles de taille supérieure à l’échelle de Kolmogorov dans une turbulence homogène isotrope calculée par DNS. Les résultats sont comparés avec l’expérience et montrent que cette description reproduit bien les propriétés dynamiques des particules. / The main objective of this thesis is to study the acceleration of fluid and inertial particles moving in a turbulent flow under the influence of a homogeneous shear in order to develop LES stochastic models that predict subgrid acceleration of the flow and acceleration of inertial particles. Subgrid acceleration modelisation is done in the framework of the LES-SSAM approach which was introduced by Sabel’nikov, Chtab and Gorokhovski[EPJB 80:177]. Acceleration is predicted with two independant stochastic models : a log-normal Ornstein-Uhlenbeck process for the norm of acceleration and an Ornstein-Uhlenbeck process expressed in the sense of Stratonovich calculus for the components of the acceleration orientation vector. The approach is used to simulate fluid and inertial particles moving in a homogeneous isotropic turbulence and in a homogeneous sheared turbulence. Our results show that small scales statistics of particles are better predicted in comparison with classical LES approach. Modelling of inertial particles acceleration is done in the framework of the LES-STRIP which was introduced by Gorokhovski and Zamansky[PRF 3:034602] with two independant stochastic models in a similar way to the subgrid fluid acceleration. Computations of inertial particles in the homogeneous shear flow present good predicitons of the particles acceleration and velocity when STRIP model is used. In the last chapter, we present an equation to describe the dynamic of point-like particles which size is larger than the Kolmogorov scale moving in a homogeneous isotropic turbulence computed by direct numerical simulation. Results are compared with experiments and indicate that this description reproduces well the properties of the particles dynamic.

Dans un système de machines rotatives comme un moteur à turbine, les paliers d'éléments roulants à grande vitesse jouent un rôle important dans le support de l'arbre ou du rotor rotatif, et ont besoin d'une lubrification pour assurer leur fonction. Sauf qu'une petite quantité d'huile est nécessaire pour former le film lubrifiant élastohydrodynamique dans la zone de contact, la plus grande partie du lubrifiant reste en suspension dans l'air, formant un mélange huile/air. Ce phénomène entraîne des pertes hydrauliques parasitaires excessives lorsque les éléments roulants se translatent et tournent dans l'environnement fluide, ce qui peut constituer une partie relativement importante de la perte de puissance totale du roulement, appelée traînée d'enroulement et pertes de barattage. Pour une vitesse de rotation jusqu'à 3× 106 Ndm, la contribution de la traînée/dérive au total peut atteindre 50%. Cependant, jusqu'à présent, il existe peu d'approches utilisées directement pour l'estimation des pertes par traînage, qui ne pouvait fournir qu'une approximation plutôt grossière. Dans cette thèse, la méthode CFD est utilisée pour analyser d'abord l'écoulement autour d'un cylindre circulaire de longueur finie avec deux extrémités libres dans un espace ouvert. Ensuite, le modèle est remplacé par plusieurs cylindres circulaires en ligne pris en sandwich par deux parois plates, ce qui représente une approche simplifiée. Le fluide est ici considéré comme incompressible, représentant un fluide monophasé équivalent pour l'écoulement diphasique huile/air à l'intérieur de la cavité de palier avec des propriétés de fluide spécifiées. Les résultats indiquent que l'écoulement autour de l'élément de rouleau de longueur finie est perturbé par ses deux extrémités libres, les anneaux environnants, la cage et d'autres éléments roulants. Il est proposé une relation entre le coefficient de traînée et le nombre de Reynolds approprié pour un cylindre circulaire dans les roulements à rouleaux (1<L/D<6), ainsi qu'une formulation pour la prévision des pertes de barattage. L'écoulement diphasique huile/air à l'intérieur de la cavité de palier avec lubrification sous la course est également étudié dans ce travail. Le volume couplé de niveau de fluide (CLS-VOF) est utilisé pour démontrer la distribution du lubrifiant le long de la circonférence du palier. L'effet de divers facteurs est étudié, par ex. la vitesse d'injection d'huile, le diamètre de la buse, les propriétés de l'huile et l'angle d'injection de l'huile. La vitesse de rotation de tous les composants du palier est étudiée en particulier pour quantifier leur influence sur la fraction du volume d'huile à l'intérieur de la cavité du palier. Les résultats démontrent que non seulement la vitesse de rotation relative de l'anneau interne, mais la vitesse de la cage elle-même pourrait changer la distribution d'huile. / In a rotating machinery system like turbine engine, high speed rolling element bearings play an important role in supporting the rotating shaft or rotor, and need lubrication to insure their function. Except a small quantity of oil is needed to form the elastohydrodynamic lubricant film in the contact zone, most of lubricant remains in suspension in air, forming an oil/air mixture. This phenomenon leads to excessive parasitic hydraulic losses when rolling elements translate and rotate into the fluid environment, which may constitute a relatively large portion of the bearing's total power loss, named windage drag and churning losses. For high speed applications, i.e. for rotational speed up to 3× 10^6 Ndm, the contribution of drag/windage loss to the total one may reach up to 50%. However, so far there are few approaches used directly for drag and churning losses estimation, which could only provide a rather gross approximation. In this thesis, the Computational Fluid Dynamics (CFD) method is employed to analyze first the flow around one finite-length circular cylinder with two free ends in an open space. Then the model is changed to several in-line circular cylinders sandwiched by two flat walls, which represents a simplified approach. The fluid here is regarded as incompressible, representing an equivalent one-phase fluid for the oil/air two-phase flow inside the bearing cavity with specified fluid properties. The results indicate that the flow around the finite length roller element is perturbed by its two free ends, the surrounding rings, the cage and other rolling elements. A relationship between the drag coefficient and the Reynolds number suitable for circular cylinder in roller bearings (1<L/D<6) is proposed, as well as a formulation for churning losses prediction. The oil/air two phase flow inside the bearing cavity with under-race lubrication is also studied in this work. The coupled level-set volume of fluid (CLS-VOF) method is employed to demonstrate the lubricant distribution along the bearing circumference. The effect of various factors is studied, e.g. the oil injection velocity, the nozzle diameter, the oil properties, and the oil injecting angle. Rotational speed of all the bearing components are studied particularly to quantify their influence to the oil volume fraction inside the bearing cavity. The results demonstrate that not only the inner-ring relative rotational speed, but the cage speed itself could change the oil distribution. The results can be used for the precise lubrication design to optimate the oil distribution inside the bearing.

Dans ce mémoire, nous présentons une stratégie basée sur une approche hybridedans le domaine temporel, couplant une méthode de résolution des équations de Maxwelldans le domaine 3D (FDTD) avec une méthode de résolution des équations de ligne detransmission, afin de pouvoir simuler des problèmes électromagnétiques de grande échelle. Lemémoire donne les éléments d’hybridation pour deux cadres d’utilisation de cette approche :une approche multi-domaine et une approche multi-résolution ou d’échelle.L’approche multi-domaine est une extension de la méthode FDTD 3D à plusieurs sousdomainesreliés par des structures filaires sur lesquelles on résout une équation de lignes detransmission par un formalisme FDTD 1D. La difficulté est d’abord d’avoir une définitionimplicite du champ électromagnétique dans la théorie des lignes de transmission, et d’autrepart de prendre en compte les effets du sol sur les courants induits au niveau des lignes etsur les champs électromagnétiques.L’approche multi-résolution ou d’échelle est conçue pour étendre les capacités de la méthodeFDTD au traitement du routage de câbles complexes ayant une section plus petite quela taille de la cellule. Ce mémoire présente différentes techniques pour évaluer les paramètresde la ligne, basées sur la résolution d’un problème de Laplace 2D, ainsi qu’une méthode decouplage champs/câbles basée sur le courant de mode commun.L’ensemble de ce travail nous a permis de proposer une méthode numérique efficace pourcalculer les effets électromagnétiques induits par une source (type onde plane ou dipolaire)sur des sites de grande dimension, composés de plusieurs bâtiments reliés entre eux par unréseau de câbles. Dans ce cadre une application à la foudre a été réalisée. / In this thesis, we present a strategy based on a hybrid approach in the timedomain, by coupling 3D method (FDTD) with a multi-conductors transmission line (MTL)method, in order to simulate complex large scale electromagnetic problems. This reportgives the theoretical and numerical elements for coupling these approaches for two kindof problems, which are the multi domains approach and the multi scale approach. Themultiple domains approach is an extension of the classical FDTD method taking into accountseveral 3D subdomains, interconnected by a wire network, on which a 1D transmission lineformalism is used. The main issues are, on one hand to have an implicit expression ofthe electromagnetic field in the transmission line approach, and on the other hand to beable to take into account the ground effects on the induced currents, on the transmissionline parameters and on the electromagnetic field. The multi scale approach is developed toextend the capabilities of FDTD to deal with complex cables routing. We assume that thecross section of the cables are smallest than the cell size, and in these problems, the 1Dtransmission line problem is physically included in the 3D global computational domain.The work done in this thesis leaded to a new field to transmission line coupling based onthe common mode current, and an evaluation of the transmission. line parameters basedon a Laplace equation resolution in 2D. In this work, we have elaborated and proposedefficient numerical strategies for the computation of electromagnetic induced effects on largeand complex sites, composed of several interconnected distant buildings. An application tolightning problems have been done.

Les étoiles à neutrons sont parmi les objets les plus extrêmes dans l'univers. Elles sont des étoiles compactes, nées à la suite d'une explosion de supernova gravitationnelle, au point final de l'évolution stellaire. Le champ gravitationnel y est très fort, et la matière à l'intérieur atteint des densités extrêmement élevées. Elles sont donc des "laboratoires" prometteurs, non seulement pour tester le régime de champ fort en relativité générale, mais aussi pour en apprendre davantage sur la physique nucléaire à haute densité, qui actuellement ne peut pas être reproduit avec des expériences terrestres. Ainsi, les étoiles à neutrons nous permettent d'adresser des questions telles que l'existence éventuelle de particules autres que nucléons à haute-densité. À cause de la naissance violente de ces objets, les étoiles à neutrons très jeunes, que l'on appelle proto-étoiles à neutrons, sont également très chaudes, et souvent en rotation différentielle rapide. Dans cette thèse nous avons pour but de développer un modèle stationnaire d'une telle proto-étoile à neutrons.Ainsi, nous présentons une nouvelle méthode pour calculer numériquement les équations d'équilibre d'un fluide parfait relativiste, axisymétrique et stationnaire, en rotation différentielle et à température finie, valable pour une équation d'état réaliste. Nous présentons en détail le code (accessible au public) développé. Nous avons appliqué ce code avec des nouvelles équations d'état réalistes à température finie, basée sur une théorie relativiste du champ moyen, en incluant tous les degrés de liberté hyperoniques. Nous avons calculé des modèles relativistes stationnaires de proto-étoiles à neutrons en rotation différentielle rapide. Nous allons discuter les applications de nos modèles pour explorer plus en détail la physique de ces objets. / Neutron stars are among the most extreme objects in the universe. They are compact stars born as the aftermath of a core-collapse supernova explosion, at the endpoint of stellar evolution, with a strong gravitational field, and extremely high densities. They are therefore promising 'laboratories', not only to test the strong-field regime of general relativity, but also to learn about nuclear physics in the high density regime, which presently is not accessible in earth based experiments. This allows to address questions such as the possible existence of particles other than nucleons at high-densities. As a consequence of the violent birth of these objects, new-born (proto-)neutron stars are extremely hot and, in general, rapidly rotating, which raises interesting problems in the task of developing a stationary model of such objects.In this thesis, we present a new self-consistent method to numerically compute the equilibrium equations of stationary axisymmetric relativistic (differentially) rotating perfect fluids at finite temperature, with a realistic equation of state. We introduce in detail the (publicly available) code in which we implemented the described numerical scheme. We use newly developed realistic equations of state with finite temperature, which are based on density dependent relativistic mean field theory, and in which all hyperonic degrees of freedom are included, to compute realistic stationary relativistic models of rapidly differentially rotating proto-neutron stars. We discuss future applications of our code for further exploring the physics of proto-neutron stars.

Les nano-objets unidimensionnels alliés présentent des propriétés physiques spécifiques qui résultent à la fois de leur morphologie, de leur taille et de la répartition chimique des atomes. Nous exploitons un modèle d’Ising sur réseau qui rend compte en particulier des effets de ségrégation au sein de nanofils bimétalliques pour obtenir une compréhension fine des effets gouvernant cette répartition à l’équilibre.Dans une première section, nous détaillons l’équilibre d’un nanofil en fonction de sa taille et de sa composition, de manière à mettre en évidence le rôle des effets de taille finie sur la thermodynamique d’équilibre d’objets bimétalliques 1D. Contrairement aux systèmes infinis, l’équilibre dépend de l’ensemble statistique considéré. Ainsi la ségrégation est plus marquée dans l’ensemble canonique, où la concentration du nanofil est imposée, que dans l’ensemble pseudo-Grand Canonique (p-GC) où le nanofil est en équilibre avec un réservoir qui fixe la différence de potentiel chimique entre les espèces. De même, la contrainte de composition dans l’ensemble canonique induit des corrélations chimiques d’occupation des sites qui favorisent davantage les paires hétéroatomiques. Nous montrons que l’écart observé entre les isothermes des deux ensembles croît avec la courbure de l’isotherme canonique et avec l’amplitude des fluctuations de la concentration nominale dans l’ensemble p-GC. Ces fluctuations diminuant avec la taille du nanofil considéré, l’écart entre les ensembles s’annule à la limite thermodynamique. Les effets de taille finie se traduisent par ailleurs par l’apparition, à basse température et pour de petits nanofils, d’une coexistence d’un mode pur en l’espèce ségrégeante et d’un mode de faible concentration nominale constitué principalement de configurations de type cœur-coquille et Janus. Nous développons alors un formalisme permettant de caractériser cette bimodalité.Alors que les résultats évoqués précédemment concernent un nanofil considéré seul, nous étudions dans la deuxième section l’équilibre de l’ensemble des nanofils formant un co-dépôt unidimensionnel inférieur à la mono-couche. Nous montrons que la distribution en taille de ces nanofils varie globalement selon une loi de puissance, quelle que soit la composition du codépôt, de sorte que la ségrégation n’a que peu d’influence sur la microstructure observée. Par contre, en raison du rapport surface/volume et des corrélations chimiques dans ces objets, la composition des nanofils du co-dépôt varie très fortement selon leur taille, les petits nanofils étant plus riches en l’espèce ségrégeante que les plus grands. Enfin, nous étendons le diagramme de bimodalité d’un nanofil seul à l’ensemble des nanofils du co-dépôt et montrons que cette bimodalité est difficilement observable car elle ne concerne que des amas de petite taille qui sont très minoritaires du fait de la cohésion atomique. / The chemical configuration and the specific shape of 1D bimetallic nano-objects endow them with physical properties (such as magnetic ones) that strongly differ from their bulk counterparts. To get a deep insight of the parameters that govern the equilibrium configuration, we consider a rigid lattice-gas Ising model that accounts for segregation effects within bimetallic nanowires that decorate a step edge.In a first section, we detail the equilibrium of a nanowire as a function of both its size and composition in order to specify the role of finite-size effects onto the equilibrium thermodynamics of 1D bimetallic objects. Contrary to infinite systems, this equilibrium depends on the statistical ensemble to be considered. The segregation profile is indeed stiffer in the canonical ensemble where the nanowire concentration is imposed, than in the semi-grand-canonical ensemble (s-GC) where the nanowire is in equibrium with a reservoir that sets the difference of chemical potentials between the species. Moreover, the composition constraint in the canonical ensemble yields chemical correlations between occupation sites that favor heteroatomic pairs. We show that the deviation observed between the isotherms related to the two ensembles increases with the curvature of the canonical isotherm and with the amplitude of the fluctuations of the nominal concentration within the s-GC ensemble. As these fluctuations decrease with the nanowire size, the deviation between ensembles vanishes at the thermodynamical limit. The finite-size effects also imply at low temperature for small nanowires, that a pure mode of the segregating species coexists with a low-concentration mode that mainly corresponds to core-shell and Janus configurations. We develop a framework to characterize the resulting two-mode density of compositions.While the abovementioned results deal with a fixed-size nanowire, we study in the second section the equilibrium of the set of nanowires that forms a submonolayer 1D-codeposit. We show that the size distribution of these nanowires globally varies as a power law, whatever the codeposit composition, so that segregation has a slight influence onto the observed microstructure. However, due to the surface/volume ratio and chemical correlations within these objects, the composition of the nanowires of the codeposit varies strongly with their size, the smaller the richer in the segregating species. Finally we extend the two-mode diagram of the single nanowire to the set of nanowires forming the codeposit and show that this two-mode distribution is hardly visible as it concerns only short nanowires which are very rare, mainly due to atomic cohesion that is reinforced at low temperature.

Cette thèse porte sur l'étude de milieux granulaires constitués de particules molles. Elle s'appuie d'une part sur la méthode de corrélation d'images numériques (DIC) et d'autre part sur des simulations associant la méthode des éléments finis (FEM) et la méthode de Dynamique des Contacts (CD) permettant de tenir compte de la déformabilité élevée des particules et des interactions entre particules. Pour une large gamme de matériaux complexes (élastique, plastique et mousse solide), l'étude préliminaire de particules uniques comprimées radicalement valide la technique de DIC pour les mesures des champs de déformation, et pour la détermination de l'énergie de déformation pour des matériaux élastiques. Des grandeurs macroscopiques et micro-structurales tels que la compacité, la coordinence, l'énergie élastique, la densité de la probabilité de force de contact ainsi que de la densité d'énergie sont mesurées, et comparées entre les résultats numériques et expérimentaux en particulier dans le régime où les déformations des particules sont d'amplitudes finies. L'accord quantitatif entre ces approches permet de valider à la fois la méthode expérimentale pour l'étude de milieux granulaires à particules déformables, et d'autre part l'approche numérique. Forts de ces résultats, nous avons mené une étude numérique de la compression uniaxiale d'un assemblage de cylindres Néo-Hookéens. L'effet du frottement sur les paramètres macroscopiques et mico-structuraux est déterminé, y compris lorsque la compacité du système est proche de 100% / This thesis deals with the study of granular media composed of soft particles. It relies on the one hand on the method of Digital Image Correlation(DIC) and on the other hand, on simulations coupling the Finite Element Method and the Contact Dynamics method, taking into account the high deformability of particles and interactions between particles. For a wide range of complex materials (elastic, plastic and solid foam), the preliminary study of radially compressed single particles validates the DIC technique for strain field measurements, and for the determination of the strain energy for elastic materials. Macroscopic and micro-structural quantities such as packing fraction, coordination, elastic energy, probability density of the contact force as well as strain energy density, are measured and compared between the numerical and experimental results in particular in the regime where the deformations of the particles are large. The quantitative agreement between these approaches makes it possible to validate both the experimental method for the study of granular media with deformable particles, and the numerical approach. Based on these results, we conducted a numerical study of the uniaxial compression of a Neo-Hookean cylinder assembly. The effect of friction on the macroscopic and mico-structural parameters is determined, even when the packing fraction of the system is close to 100%.

Ce travail de thèse est surtout consacré à l'étude du ferromagnétisme dans les systèmes désordonnés, en particulier les semiconducteurs magnétiques dilués (DMS). Le formalisme utilisé, basé sur les fonctions de Green à température finie, est très général pour l'étude du magnétisme et du transport dans ces systèmes. Dans un premier temps, il est montré que l'approche champ moyen - RKKY (MF-RKKY) souvent utilisée n'est pas appropriée pour décrire les propriétés magnétiques des systèmes dilués. Un meilleur traitement du Hamiltonien de Heisenberg RKKY dans le cadre de la théorie RPA locale auto cohérente (SC-LRPA) a permis de montrer en particulier que l'approche MF-RKKY surestime largement les températures critiques ainsi que l'étendue de la zone de stabilité du ferromagnétisme. Dans un second temps on étudie le modèle non perturbatif « V-Jpd » par diagonalisation exacte pour chaque configuration de désordre et on calcule explicitement les échanges magnétiques. En dehors de la limite perturbative, ces échanges n'ont pas le caractère RKKY. Ensuite, le Hamiltonien de Heisenberg effectif est traité dans le cadre de la théorie SC-LRPA. Cette approche en 2 étapes montre en particulier (i) l'importance du désordre et des fluctuations thermiques et transverses et (ii) que le potentiel coulombien V joue un rôle crucial pour comprendre l'origine du ferromagnétisme dans les DMS. Ce modèle minimal tient compte de la percolation et des diffusions multiples des porteurs itinérants sur les impuretés et permet d'unifier la description des DMS. Enfin, une étude numérique des effets de taille finie et de l'importance de l'échantillonnage statistique a permis de montrer les insuffisances sévères du traitement Monté Carlo « complet » du modèle dilué « V-Jpd ». Ce modèle microscopique permet de combler le fossé entre les approches modèles trop simplistes et celles basées sur les calculs ab initio.

La taille finie des nanoparticules leur confère des caractéristiques spécifiques, tant en surface qu'en leur coeur. Dans ce contexte, nous avons étudié la thermodynamique des agrégats bimétalliques à l'aide d'un modèle d'Ising traité par simulations Monte Carlo et également dans l'approximation de champ moyen. Nous avons considéré le système Cu-Ag, prototype des systèmes à forte tendance à la démixtion. Pour un cuboctaèdre de quelques milliers d'atomes, les simulations Monte Carlo ont montré que les différentes facettes de la surface sont le lieu d'un équilibre dynamique à suffisamment basse température. Ainsi, au sein d'un même agrégat, nous prédisons la coexistence de facettes (100) riches en cuivre et riches en argent, ce résultat s'étendant aux facettes (111) et mêmes aux arêtes. Un équilibre dynamique similaire est également observé pour les couches de coeur, équivalent de la démixtion pour le système infini. L'étude de la chaîne linéaire nous a permis de déterminer les facteurs régissant la compétition entre le mouillage (enrichissement en Ag partant de la surface et gagnant progressivement le coeur) et l'établissement d'un équilibre dynamique initié au coeur. L'influence de la taille de la chaîne sur cette compétition a été spécialement étudiée, ce qui conduit à l'ébauche d'un diagramme de phase de milieu fini intégrant la ségrégation superficielle. Nous avons complété cette étude sur réseau par la prise en compte des relaxations atomiques, via un potentiel interatomique à N-corps, dans le cas de l'icosaèdre. Nous illustrons les effets des contraintes propres à cette structure, qui se manifestent par de fortes hétérogénéités de compositions intra et inter couches.