Dissertations / Theses on the topic 'Texturation de surface multi-échelles'

Une goutte d’eau roule sur une feuille de lotus sacré mais elle adhère à un pétale de rose. Ces deux surfaces présentent, aux échelles micrométrique et nanométrique, une morphologie complexe. Quel est alors le lien entre leur mouillabilité et leur morphologie multi-échelles ? L’objectif de cette thèse est de mener une approche biomimétique pour la conception de surfaces superhydrophobes tout en comprenant les stratégies mises en œuvre par le vivant. Dans une première partie, nous caractérisons des surfaces végétales dont le régime de mouillage est observé directement par microscopie confocale à balayage laser. Nous démontrons que le lotus sacré produit un mouillage en régime mixte métastable caractérisable par une profondeur d’ancrage non nulle des lignes triples à l’équilibre. De plus, nous observons un régime hiérarchique de Wenzel-Wenzel sur le pétale de rose, à rebours des hypothèses de la littérature. De ces observations, nous tirons des questionnements clés que nous confrontons aux modèles actuels de la littérature. Dans une seconde partie, nous adaptons deux approches des phénomènes capillaires à l’étude du piégeage d’air sur une topographie multi-échelles et introduisons la paramétrisation nécessaire à l’étude des régimes de mouillage mixte et de leur robustesse. Nous prédisons la profondeur d’ancrage à l’équilibre sur le lotus sacré en mettant en lumière le rôle de sa topographie nanométrique. Enfin, nous décrivons les mécanismes gouvernant les mouvements à l’avancée et à la reculée et leurs propagations par récurrence à travers les échelles topographiques constituant une surface en introduisant la notion de mouvement précurseur. Nous démontrons l’effet de la profondeur d’ancrage à l’équilibre sur l’hystérèse d’angle de contact et le rôle des sous-échelles topographiques sur la robustesse du régime de piégeage d’air. À travers l’étude expérimentale de surfaces fabriquées par photolithographie, nous confrontons ce modèle à la réalité. Enfin, dans une troisième partie, nous transposons les conclusions issues de ce modèle en un cahier des charges pour la conception de surfaces superhydrophobes robustes, déclinons la stratégie mise en œuvre par la feuille du lotus sacré et proposons deux voies de fabrication de surfaces déperlantes, par recristallisation de cire naturelle et polymérisation deux-photons
A drop of water rolls on the sacred lotus leaf but stay fiercely anchored onto a rose petal. Both surfaces display a complex morphology at the micrometric and nanometric scales. Therefore, one could ask: how are their wettability and their morphology related? The purpose of this dissertation is to carry out a biomimetic approach in order to conceive superhydrophobic surfaces and to better understand nature’s strategies. In a first part, vegetal surfaces have been characterized by directly observing the wetting state they produce with the help of confocal microscopy. We demonstrate the fact that the sacred lotus produces a metastable mixed-state wetting that is characterized by a finite equilibrium anchorage depth of triple lines. On the other hand, a Wenzel-Wenzel hierarchical wetting state is observed on the rose petal, in spite of what literature suggests. From these experiments, key questions have been highlighted and confronted to the current models available within the literature. In a second part, two approaches to capillary phenomena have been adapted to the study of a composite wetting state produced by a multiscale topography. We introduce a complete parameterization allowing us to tackle the problem of the mixed-state wetting and its stability, to predict the value of the equilibrium anchorage depth on the sacred lotus leaf and to identify the contribution of its nanoscale topography to its wetting. Then, we thoroughly describe the mechanisms underlying the advancing and receding motions of triple lines and their recursive propagation across every topographical scale constituting a surface by introducing the notion of precursor motion. We highlight the effect of the equilibrium anchorage depth on the contact angle hysteresis and the role played by topographical subscales on the robustness of the composite wetting state. Through the experimental study of model surfaces manufactured by photolithography, we compare our predictions to reality. Eventually, in a third part, the conclusions drawn from our model are transposed into technical specifications for the conception of robust superhydrophobic surfaces, the strategy of the sacred lotus leaf is thoroughly described and two promising manufacturing processes are proposed through the recrystallization of natural wax and two-photon polymerization

Nous présentons une étude détaillée, non-perturbante, de la prise du plâtre par Relaxation Magnétique Nucléaire (RMN) du proton, à bas champ magnétique. Cette technique permet de quantifier en continu le degré d'hydratation et de déduire de manière non-destructive les paramètres microstructuraux d'une pâte de plâtre en cours de durcissement. La décroissance bi-exponentielle des signaux de relaxation de l'eau dans le plâtre prouve l'existence de deux populations d'eau différentes, en échange lent. A partir des mesures de relaxométrie RMN, deux modes d'organisation de la microstructure sont identifiés en fonction du rapport de gâchage initial eau/plâtre (e/p), pour 0.4 <= e/p <= 0.6 et 0.7 <= e/p <= 1. Un modèle original d'échange entre les populations d'eau dans un milieu poreux a été établi. Une nouvelle approche multi-échelles et multi-techniques de la prise et du durcissement des pâtes de plâtre est proposée, permettant de relier les mesures locales obtenues par RMN aux propriétés mécaniques du matériau (mesurées par ultrasons). La démarche expérimentale et la modélisation présentées dans cette étude peuvent être appliquées à une large gamme de matériaux poreux en évolution.

Le palier hydrodynamique est une solution de guidage en rotation particulièrement appréciée pour ses caractéristiques d’amortissement à hautes vitesses de rotation. Cependant les performances des machines tournantes lubrifiées par un film fluide sont impactées par des effets non linéaires difficiles à analyser. La prédiction du comportement du système par la simulation nécessite une modélisation avancée de l’écoulement de lubrifiant dans le palier hydrodynamique. Enfin, l’état de surface semble avoir un impact important sur l’écoulement du fluide lubrifiant, lui-même agissant sur les caractéristiques statiques et dynamiques des parties tournantes. Cette étude vise à améliorer les modèles numériques liés à l’impact de l’état de surface des paliers hydrodynamiques sur la dynamique de ligne d’arbre. La méthode d’homogénéisation multi-échelles a été utilisée à cet effet dans un algorithme multi-physiques pour décrire l’interaction entre la structure flexible en rotation et les films fluides des supports de lubrification. Différents modèles ont été utilisés pour prendre en compte la présence de zone de rupture de film lubrifiant. Des méthodologies non-linéaires fréquentielles ont été mises en place afin de permettre l’étude paramétrique des solutions périodiques d’un tel système et de leur stabilité. Afin de confronter ce modèle complexe à la réalité, un banc d’essai miniature a également été conçu. Différents échantillons présentant des états de surface modifiés par ablation à l’aide de LASER femto-seconde ont été testés. L’étude expérimentale a permis de vérifier certaines tendances prévues par la simulation. Des améliorations des performances des paliers hydrodynamiques par rapport aux vibrations auto-entretenues du système ont été démontrées pour certaines textures. En revanche toutes les améliorations ne sont pas prédites par les algorithmes d’homogénéisation multi-échelles. La présence de recirculation dans les aspérités du motif a été mise en évidence à partir de la résolution locale des équations de Navier-Stokes. Ce résultat participe à la remise en question des hypothèses classiques utilisées en texturation, et peut justifier les améliorations obtenues expérimentalement avec les paliers texturés
The hydrodynamic bearing provides good damping properties in rotating machineries. However, the performances of rotor-bearings systems are highly impacted by nonlinear effects that are difficult to analyze. The rotordynamics prediction requires advanced models for the flow in the bearings. The surface of the bearings seems to have a strong impact on the lubricant flow, acting on the static and dynamic properties of the rotating parts. This study aims to enhance the simulation of the bearings’ surface state effect on the motion of the rotating shaft. The flexible shaft interacts with textured hydrodynamic bearings. Multi-scales homogenization is used in a multi-physics algorithm in order to describe the fluid-structure interaction. Different models are used to account for the cavitation phenomenon in the bearings. Nonlinear harmonic methods allow efficient parametric studies of periodic solutions as well as their stability. Moreover, a test rig has been designed to compare predictions to real measurements. Several textured shaft samples modified with femto-seconds LASER surface texturing are tested. In most cases the experimental study showed similar results than the simulation. Enhancements of the vibration behaviors of the rotor-bearing system have been revealed for certain texturing patterns. The self-excited vibration, also known as "oil whirl" phenomenon, is stabilized on a wide rotating frequency range. However, the simulation tool does not predict well the enhancements that are observed. Vortices in surface texturing patterns have been revealed numerically with Navier-Stokes equation resolution. These results are opposed to the classical lubrication hypothesis. It is also a possible explanation of the enhancements that are experimentally measured with textured bearings

Une voie prometteuse pour améliorer le rendement des centrales solaires à tour consiste à chauffer de l'air pressurisé à haute température afin d'alimenter un cycle thermodynamique de Brayton. Pour cela, il est indispensable de concevoir des récepteurs solaires performants,permettant de forts transferts de chaleur vers le fluide. Le développement de tels récepteurs passe par une compréhension fine de leurs écoulements internes. Il s'agit d'écoulements complexes, combinant de hauts niveaux de turbulence et un fort gradient de température entre la paroi irradiée par le flux solaire concentré et la paroi arrière isolée. On se propose dans ce travail de réaliser une étude amont numérique et expérimentale de ce type d'écoulements.D'une part, des mesures de vitesse par SPIV (vélocimétrie par images de particules stéréoscopique) sont effectuées dans une soufflerie de canal plan turbulent lisse dont la cellule de mesure est représentative d'un récepteur solaire surfacique. On observe en particulier l'influence d'un chauffage asymétrique sur les statistiques de la turbulence. Ces mesures sont d'autre part complétées par des simulations fines LES (simulation des grandes échelles)menées dans les conditions de la soufflerie. Pour finir, une simulation LES d'un canal plan texturé sur une paroi par une géométrie innovante est conduite. Cette architecture interne du récepteur combine des générateurs de tourbillon et des riblets afin d'intensifier les échanges de chaleur vers le fluide
A promising line of research to increase the efficiency of solar tower power plants consists in heating pressurized air to high temperatures in order to fuel a Brayton thermodynamic cycle. This requires to design effective solar receivers that allow for intense heat transfers toward the fluid. To develop such receivers, an in-depth understanding of their internal flows is needed. These are complex flows, combining strong turbulence and strong temperature gradient between the concentrated sun irradiated wall and the back insulated wall.The aim of this work is to investigate numerically and experimentally such flows.On one hand, velocities are measured by SPIV (Stereoscopic Particle Image Velocimetry) in a turbulent channel flow wind tunnel whom measurement cell is similar to a surface solar receiver. The influence of an asymmetric heating on the turbulence statistics are especially investigated. These measurements are supplemented by Large Eddy Simulations run under the same conditions as the wind tunnel. Finally, a Large Eddy Simulation is run in a channel flow textured on one wall by an innovative geometry. This internal receiver design combines vortex generators and riblets in order to enhance the heat transfers

Il est communément admis aujourd'hui que le contact entre deux surfaces se compose en réalité d'une multitude de contact ponctuels entre des aspérités. Cette considération amène à une surface de contact réelle significativement différente de l'aire de contact parfaite supposée dans la théorie de Hertz.De même, elle implique également la présence d'un espace libre entre les deux surfaces en contact. Dans cette situation, l'objectif principal de ces travaux de thèse est de développer des approches numériques permettant l'analyse du contact mécanique entre surfaces rugueuses dans le but de qualifier/quantifier l'étanchéité de ce contact rugueux.Deux approches différentes sont étudiées. La première consiste à analyser le contact mécanique entre une surface rugueuse et un plan rigide au moyen de la méthode des éléments finis et d'un nouveau modèle numérique. La seconde concerne l'estimation de la transmissivité d'un contact rugueux en considérant des simulations de l'écoulement d'un fluide au sein du champ des ouvertures présent entre les deux surfaces en contact.La comparaison de ces estimations numériques avec les résultats expérimentaux révèlent des écarts importants. Dans le but de comprendre ces écarts, l'influence du modèle de comportement matériau dans de telles simulations est étudiée. La plasticité cristalline, mais également l'élévation de la température dégagée par déformation plastique seront considérés. La question de la représentativité de notre problème vis-à-vis de l'approche fluide sera également discutée.

Quelques cas de fissurations de composants en acier inoxydable austénitique 304L dus à un phénomène de fatigue thermique ont été recensés, en particulier sur les circuits de Refroidissement du Réacteur à l'Arrêt (RRA) des Réacteurs à Eau Pressurisée (REP). EDF a lancé un programme de R&D afin de comprendre et d'évaluer les risques de dommages dans ces circuits. Le test développé à EDF, INTHERPOL, est un essai de fatigue thermique de type structure sur maquette tubulaire et soumis à un chargement périodique. Ces essais sont réalisés sous différents niveaux de chargement, états de surface et soudures afin de rendre compte de ces paramètres sur l'amorçage des fissures. L'objectif de cette étude est de définir une méthodologie permettant la justification d'un critère de fatigue tenant compte de l'influence de l'état de surface, à l'aide d'une approche multi-échelle. La première partie du travail consiste à caractériser le matériau (304L) et particulièrement sa surface. Les premiers résultats de la caractérisation de la surface montrent un fort gradient d'écrouissage initial du matériau sur 250 microns en profondeur. Ce gradient n'évolue pas avec le cyclage thermique. Pour cette caractérisation, des mesures de microduretéet des observations de la microstructure (MET) ont été réalisées. La deuxième partie est la modélisation macroscopique des essais INTHERPOL permettant de déterminer les champs de contrainte et déformation issus des cycles thermiques. La troisième partie du travail consiste ensuite à évaluer l'effet de la rugosité de surface et du gradient d'écrouissage à échelle plus fine à l'aide d'une modélisation à l'échelle du polycristal, basée sur la variation des densités de dislocations. Un objectif pour l'avenir est d'aider, par ces modélisations polycristallines, au choix de critères de fatigue prenant en compte l'influence de l'état de surface (critères énergétiques ou de types plan critique par exemple).

L'océan est le siège de mouvements complexes à toutes échelles spatiales et temporelles. Au sein d'une circulation moyenne et globale existe une circulation secondaire peuplée de fronts, de méandres, de jets étroits, de tourbillons, que l'on nomme circulation à mésoéchelle. L'observation spatiale permet une description et une évaluation synoptique de cette dynamique à mésoéchelle au moyen de l'altimétrie et la diffusiométrie. Cette évaluation a été le premier objectif de cette thèse et a permis de développer un produit distribué à la communauté scientifique internationale : le produit GEKCO. Cependant la description des processus submésoéchelle à plus fine résolution nécessite l'utilisation de données à super-résolution (couleur de l'eau, température de surface) qui ont la possibilité de représenter toute la complexité d'un océan en régime de turbulence pleinement développée. Une méthode à la croisée de l'océanographie physique et de la "science de la complexité" utilisant la formulation microcanonique de la cascade multiplicative, le produit GEKCO et des images de température de la mer, a fait l'objet de la seconde partie de ce manuscrit. La dynamique océanique étant la clef de voûte de tout le monde marin du vivant, la dernière partie de cette thèse s'est intéressée à l'impact de la circulation à mésoéchelle et à submésoéchelle sur la chaîne trophique marine en se focalisant sur ses deux extrémités. L'étude de la circulation à submésoéchelle a permis de montrer qu'elle joue un rôle prépondérant pour la biomasse marine ; un rôle d'activateur en océan ouvert et un rôle d'inhibiteur dans les systèmes d'upwelling de bord Est. Différentes études sur les trajets de prédateurs marins supérieurs ont démontré la nécessité de prendre en compte la dynamique océanique pour interpréter leur comportement de navigation
The ocean is governed by complex movements at all spatio-temporal scales. Within an mean and global circulation exists a secondary circulation inhabited by fronts, meanders, narrow jets, eddies, named mesoscale circulation. Observation with satellites allows a description and a synoptic evaluation of this dynamics at mesoscale by using altimetric and scatterometric data. This evaluation was the first objective of the thesis and allowed to develop a product distributed to the international scientific community : the GEKCO product. However, the description of submesoscale processes at finer resolution requires the use of super-resolution data (like ocean color or sea surface temperature) which have the ability to represent all the complexity of an ocean in fully developed turbulence. A method between physical oceanography and "science of the complexity" using the microcanonical formulation of the multiplicative cascade, the GEKCO product and sea surface temperature images was the topic of the second part of this manuscript. The ocean dynamics being the keystone of all marine ecosystems, the last part of this thesis was dedicated to the meso and submesoscale impact of the ocean dynamics to the trophic chain by focusing on its two extremities. Study of the flow at submesoscale allowed to show that it plays a preponderant role on the marine biomass ; role of activator in open ocean and role of inhibitor in the east boundary upwelling systems. Various studies on top marine predator paths demonstrated the necessity of taking into account the ocean dynamics to interpret their navigational behavior

Ce travail vise à caractériser et comparer les fonctionnements hydrodynamiques à plusieurs échelles spatiales en zones cristalline et sédimentaire de l’Ouest nigérien et par suite à proposer un modèle simple de fonctionnement hydrodynamique des bassins expérimentaux cristallins qui soit potentiellement transposable aux échelles supérieures.L’analyse qualitative du paysage montre que ces deux domaines ont des états de surface communs et spécifiques.L’étude expérimentale effectuée sur les états de surface communs (surface biologique, BIOL, d’érosion, ERO, structurale, ST-jac, et cultivée, C) montre, à l’échelle ponctuelle, que la surface ERO a la même valeur de conductivité hydraulique K dans les deux contextes. En revanche, les surfaces C et surtout ST ont des conductivités plus faibles en zone cristalline. Le suivi temporel effectué sur ces deux états de surface en zone cristalline a mis en évidence la stationnarité de la conductivité sur ST (jachère de 5 ans) durant la saison des pluies tandis qu'elle varie fortement sur la surface cultivée en fonction du calendrier cultural et du cumul de pluie. K décroit d’une valeur initiale de 170 mm/h après le sarclage à 20 mm/h (soit la conductivité mesurée sur ST) lorsque la surface sarclée reçoit un cumul de pluie de l’ordre de 70 mm. La conductivité atteint même la valeur de 10 mm/h après un cumul de pluie de 180 mm. Cette variation de la conductivité montre l’avantage à court terme du sarclage qui améliore l’infiltration mais qui à long terme tend à dégrader la surface.Les résultats obtenus à l’échelle de la surface élémentaire (10 m²) valident bien les mesures ponctuelles. La surface ERO a le même coefficient du ruissellement (Kr) en zones cristalline et sédimentaire tandis que les surfaces ST et cultivée ont un Kr plus élevé en zone cristalline. L’analyse de l’évolution temporelle du Kr de la surface cultivée du socle a montré que sur 63 événements pluvieux observés entre 2011 et 2013, 22 ont un Kr supérieur à la moyenne (qui est de 0.25) dont près de 2/3 sont observés après plus de 70 mm de pluie qui suivent le sarclage. Les résultats obtenus à ces deux échelles (ponctuelle et surface élémentaire) sont donc cohérents et montrent que la surface cultivée s’encroûte et peut ruisseler plus que la surface ST et autant que la surface ERO. A l’échelle du bassin versant (5 ha), les Kr sont plus élevés sur les bassins cristallins à cause de ces fortes valeurs de Kr des surfaces élémentaires mais aussi parce qu’ils sont composés d’autres surfaces à forte capacité ruisselante que sont la surface d’affleurement du socle altéré et la surface gravillonnaire ayant un Kr de 0.58. A ces trois échelles (ponctuelle, élémentaire, petit bassin), on note la non-dépendance du fonctionnement hydrodynamique à l'état hydrique initial.L’exploitation des résultats ponctuels obtenus sur le site cristallin (conductivité moyenne de BIOL, ERO, et ST et conductivité variable sur la surface cultivée) dans le modèle de Green et Ampt a permis de caler le potentiel de front par état de surface et décrire de manière très satisfaisante les ruissellements mesurés sur les surfaces élémentaires. Sur la base des ruissellements ainsi calculés, nous avons simulé les hydrogrammes à l’exutoire des bassins expérimentaux en assimilant le fonctionnement de ces derniers à celui des surfaces élémentaires en spatialisant leur infiltrabilité. En faisant l’hypothèse d’une ré-infiltration nulle sur les versants, nous avons appliqué une fonction de transfert simple prenant en compte la distance de chaque surface élémentaire par rapport au réseau hydrographique, une vitesse d'écoulement constante et une pluie imbibante de 3 à 4 mm devant saturer la couverture sableuse de ravine, ce qui est beaucoup moins qu'en zone sédimentaire. Finalement, les hydrogrammes simulés reproduisent assez bien les caractéristiques des hydrogrammes mesurés, ce qui offre une perspective d’application de certains principes du modèle sur de plus grands bassins
This work aims at characterizing and comparing the hydrodynamical functioning at several spatial scales within the granitic-basement and sedimentary zones of Western Niger. Then, a simple hydrological model that could be suitable for use at larger scales is proposed and tested.Qualitatively, the two geological domains have common and specific surface features.The experimental work carried out onto common surface features (biological crust, BIOL; erosion crust, ERO; fallow structural surface, ST and cultivated, C) shows that, at the point scale, ERO has the same hydraulic conductivity K value in both contexts. On the other hand, surfaces features C and especially ST have lower K values in granitic context.Monitoring of the ST and C sites along the rainy season proved the stationarity of the ST conductivity value. On the contrary, K varies widely with the amount of rain received from an initial value of 170 mm/h after weeding down to 20 mm/h (i.e. the ST measured value) after 70 mm of rain and even 10 mm/h after 180 mm of rain. This variation shows the short-term benefit of weeding onto infiltration but a degradation of the soil surface on the long term.At the plot scale (10 m2), runoff measurements are consistent with point measurements. ERO has the same runoff coefficient (Kr) in granitic and sedimentary zones while ST and C surfaces have a higher Kr in granitic context.Runoff monitoring of the granitic site cultivated plots showed that from a total of 63 rain events between 2011 and 2013, 22 had a Kr value higher than the average value (0.25) from which 2/3 are observed after the surface had received more than 70 mm rain after weeding.Results obtained at the two scales (point and 10-m2 plot) are thus consistent and show that the cultivated surface gets crusted and may produce runoff more than fallow ST sites and as much as ERO features.At the basin scale (5 ha), Kr values are higher in the granitic site, not only because of the higher Kr value for a given surface feature but also because of the specific low-infiltrating surfaces which are granite outcrops and gravel crusts (Kr = 0.58).At the three previous scales (point, plot and small basin), runoff volume was found independent of soil initial moisture.Using the previous point-scale results in a Green-Ampt infiltration model led to calibrate the wetting front pressure head for each surface feature and to satisfactorily describe runoff volumes obtained at the plot scale.By estimating runoff with the Green-Ampt infiltration model at any given point, basin-scale hydrograms were obtained by adding the contribution of all elementary surfaces. Assuming no re-infiltration of runoff water within the basins, a simple transfer function was chosen accounting for the distance of each surface to the hydrological network, a constant water velocity of 0.05 m.s-1 and a volume of 3-4 mm of water necessary to fill the kori sand cover, which is much less than that in the sedimentary context. Finally, simulated hydrograms reproduce nicely the measured ones, which offers the perspective of applying some principles of the model to larger basins

L'océan est le siège de mouvements complexes à toutes échelles spatiales et temporelles. Au sein d'une circulation moyenne et globale existe une circulation secondaire peuplée de fronts, de méandres, de jets étroits, de tourbillons, que l'on nomme circulation à mésoéchelle. L'observation spatiale permet une description et une évaluation synoptique de cette dynamique à mésoéchelle au moyen de l'altimétrie et la diffusiométrie. Cette évaluation a été le premier objectif de cette thèse et a permis de développer un produit distribué à la communauté scientifique internationale : le produit GEKCO. Cependant la description des processus submésoéchelle à plus fine résolution nécessite l'utilisation de données à super-résolution (couleur de l'eau, température de surface) qui ont la possibilité de représenter toute la complexité d'un océan en régime de turbulence pleinement développée. Une méthode à la croisée de l'océanographie physique et de la "science de la complexité" utilisant la formulation microcanonique de la cascade multiplicative, le produit GEKCO et des images de température de la mer, a fait l'objet de la seconde partie de ce manuscrit. La dynamique océanique étant la clef de voûte de tout le monde marin du vivant, la dernière partie de cette thèse s'est intéressée à l'impact de la circulation à mésoéchelle et à submésoéchelle sur la chaîne trophique marine en se focalisant sur ses deux extrémités. L'étude de la circulation à submésoéchelle a permis de montrer qu'elle joue un rôle prépondérant pour la biomasse marine ; un rôle d'activateur en océan ouvert et un rôle d'inhibiteur dans les systèmes d'upwelling de bord Est. Différentes études sur les trajets de prédateurs marins supérieurs ont démontré la nécessité de prendre en compte la dynamique océanique pour interpréter leur comportement de navigation.

Le rétrécissement continu des dispositifs microélectroniques exige la production des couches nanométriques uniformes et conformes, avec une pureté chimique et des interfaces abruptes. Le dépôt de couche atomique (ALD) est un procédé favorable à la production de tels films. Tirant ses avantages de la nature auto-limitante des réactions, ALD peut permettre un contrôle de l’épaisseur à la monocouche, produisant des films de haute pureté. Bien que l'ALD présente de nombreux avantages, des inconvénients se rencontrent lors du dépôt de films de quelques nanomètres. En particulier, la croissance initiale en îlots et la formation d'une couche interfaciale sont deux de ses limitations principales, en particulier dans le cas du dépôt d’oxyde métallique sur Si. De plus, le dépôt sur des grandes surfaces n'est pas toujours uniforme et dépend du réacteur et des conditions du procédé. Ces inconvénients doivent être supprimés afin de déposer des films d'oxydes sur Si, essentiels pour la production de transistors à effet de champ du futur. Dans cette thèse, l'ALD de Al2O3 de TMA et H2O sur Si est étudiée de manière approfondie, afin de remédier aux inconvénients ci-dessus. L'étude consiste en une approche multi-échelles numérique et expérimentale combinée. Quatre modèles numériques différents ont été développés pour traiter différentes échelles d'espace. Un ensemble de techniques de caractérisation a été utilisé, notamment l'ellipsométrie, XRR, TEM, STEM, EDX, XPS et SIMS. Dans ce cadre, les phénomènes détaillés sont illuminés, ce qui permet de comprendre le processus et l'origine des inconvénients de l'ALD. La compétition entre la désorption et les réactions de surface, s'est avérée être le facteur limitant pour le dépôt à basse température, jusqu'à 200°C. La concentration des sites réactifs en surface limite le dépôt à des températures supérieures à 300°C. Bien que l’ALD soit conçu comme un processus dépendant uniquement de la chimie de surface, l’analyse des phénomènes de transport à l’intérieur du réacteur a montré que la conception du réacteur et du processus peut affecter la distribution des réactifs et la température à l’intérieur du réacteur ALD. L'approche multi-échelles et le couplage entre les différents modèles numériques ont révélé que l'interaction entre les mécanismes de surface et les phénomènes de transport avait un effet sur l'uniformité du dépôt. En utilisant cette approche numérique, il est possible de dériver des conditions optimales garantissant une uniformité du film. Au cours des premières étapes, le dépôt du film est inhibée, ce qui a conduit à un régime de croissance en îlots. L'analyse intégrée a montré que 25 cycles d'ALD sont nécessaires pour déposer un film continu de Al2O3. Au cours de ce régime, l'oxydation interfaciale a conduit à la formation d'une couche d'oxyde interfacial d'environ 2 nm, composée de SiOx, AlOx et SixOyAl, qui altère les propriétés et donc les applications potentielles de la structure déposée. Un prétraitement in situ au plasma N2-NH3 du substrat de Si a été introduit, conduisant à la formation d'une couche de SixNyH sur la surface du substrat. Le prétraitement a augmenté la réactivité de surface, et la période d’inhibition était limitée. Une croissance linéaire est obtenue après 5 cycles. En outre, l'oxydation interfaciale du Si a été réduite, car la couche SixNyH s'est avérée servir de barrière efficace pour la diffusion de l' O et l'oxydation du Si. Le travail présenté dans cette thèse montre la nécessité de telles approches intégrées pour analyser les phénomènes impliqués dans l'ALD. Telles études permettent une compréhension approfondie des mécanismes, afin de proposer des solutions pour lutter contre les inconvénients apparus lors des premières étapes de dépôt. Cela pourrait permettre a l’ ALD de produire des couches minces nanométriques uniformes et conformes de grande pureté avec des interfaces abruptes, capables de répondre aux exigences de l’industrie électronique
The constant shrinking of microelectronic devices requires the production of conformal and uniform nanometric thin films, with a high chemical purity and abrupt interfaces. In this context, Atomic Layer Deposition (ALD) has emerged as a favorable process to produce such films. Drawing its advantages from the self-limiting nature of the surface reactions involved, ALD can yield thickness control down to the monolayer, producing conformal films of high purity. Although ALD has many advantages, drawbacks arise when depositing films of some nanometers. In particular, the initial island growth and the formation of an interfacial oxide layer are two of its main limitations, especially for the case of metal oxide ALD on Si. Moreover, the deposition on large area wafers is not always uniform, and depends on the reactor and process design. These drawbacks need to be suppressed in order to establish ALD as the adequate process for the deposition of high-k gate oxides on Si, essential for the production of field effect transistors of the future. In this thesis, the ALD of Al2O3 from TMA and H2O on Si is thoroughly investigated, in order to tackle the above drawbacks. The investigation consists of a combined multiscale computational and experimental approach. Four different numerical models were developed dealing with different space scales. A complete set of characterization techniques was used, including ellipsometry, XRR, TEM, STEM, EDX, XPS and SIMS. Using this framework, the detailed phenomena involved are illuminated, thus allowing to better understand the process and identify the factors responsible for the drawbacks of ALD. The competition between surface mechanisms, namely desorption and surface reactions, was found to be the limiting factor for deposition at low temperatures, up to 200oC. The concentration of surface reactive sites was found to limit the deposition at higher temperatures up to 300oC. Although ALD is conceived as a process depending only on surface chemistry, the analysis of the transport phenomena inside the ALD chamber showed that the reactor and process design can affect the reactant and temperature distribution inside the ALD reactor. The multiscale approach and the coupling among the different computational models revealed that the interplay between surface mechanisms and transport phenomena affects the film uniformity. Using this computational approach, it was possible to derive optimal process conditions that ensure maximum film uniformity. During the first deposition steps, the film deposition was found to be inhibited, leading to an island growth regime. The integrated analysis showed that 25 cycles are needed in order to deposit a continuous Al2O3 film. During this regime, interfacial oxidation of the Si substrate led to the formation of a ~2 nm interfacial oxide layer, consisting of SiOx, AlOx, and Al-silicates, which degrades the properties and thus the potential applications of the deposited structure. An in situ N2-NH3 plasma pretreatment of the HF-cleaned Si substrate was introduced, leading to a formation of a SixNyH layer on the substrate surface. The pretreatment was found to enhance the surface reactivity, as the inhibition period was restricted and linear ALD growth was obtained even after 5 cycles. Furthermore, interfacial Si oxidation was reduced, as the SixNyH layer was found to serve as an effective barrier for O diffusion and Si oxidation. The work presented in this thesis demonstrates the necessity of such integrated approaches to analyze the detailed phenomena involved in ALD. Such studies help in the thorough understanding of the ALD mechanisms, and consequently in elaborating solutions which restrict the drawbacks arising during the initial deposition steps. This could pave the way for the ALD process to industrially produce uniform and conformal nanometric thin films of high purity and abrupt interfaces, able to answer to the demands of the future electronic industry

Le contrôle et à la maîtrise de l’état des surfaces est un besoin majeur pour les industriels. De nombreuses études sur les interactions entre la morphologie de surface et les mécanismes physiques, chimiques ou mécaniques, ont été réalisées. Cependant une caractérisation plus précise en fonction des domaines et des besoins est nécessaire. Elle consiste à chercher les paramètres de rugosité les plus pertinents qui relient la topographie d’une surface aux phénomènes physiques qu’elle subit ou aux propriétés du matériau dont elle composé.Dans ce travail, un logiciel pour caractériser l’état de surface a été développé. Cet outil nommé « MesRug » permet de calculer des paramètres de rugosité et d’extraire les plus pertinents ainsi que de définir l’échelle la plus adéquate pour une application donnée. La recherche des paramètres les plus pertinent se fait par une approche statistique (l'analyse de la variance ‘ANOVA’ combinée avec la théorie du Bootstrap).Une caractérisation a été effectuée en utilisant des données de mesures (2D) sur des surfaces abrasives. L’influence de la forme des ondelettes discrètes et continues sur la détection de l’échelle pertinente du mécanisme d’abrasion a été testée. On déduit que la décomposition en ondelettes permet de quantifier et de localiser les échelles de l'abrasion des processus d'usinage pour tous les paramètres du processus. Cependant, la pertinence de caractériser les échelles appropriées d'abrasion ne dépend pas de la forme de l'ondelette.Dans ce travail, un nouveau paramètre de rugosité 3D est proposé pour quantifier la régularité d'une surface indépendamment de l'amplitude et des unités de longueur de balayage. L'efficacité de ce paramètre est testée sur des surfaces périodiques bruitées avec différents degrés d'anisotropie. La valeur de ce paramètre est comprise entre zéro (bruit parfait) et 100% (surface sinusoïdale parfaite). Il nous a permis de détecter les directions d'anisotropie de régularité pour une surface donnée
Monitoring and control of the state of the surfaces is a major need for industry. Numerous studies on the interactions between the surface morphology and the physical, chemical or mechanical mechanisms have been conducted. However, a more precise characterization related to industrial domains and needs is necessary. It consists in finding the most relevant roughness parameters that connect the topography of a surface with the physical phenomena which it undergoes or in the properties of the material of which it consisted.In this work, a software designed to characterize the surface condition was developed. This tool named "MesRug" allows to calculate roughness parameters then extract the most relevant ones and to define the most appropriate scale for a given application. The search for the most relevant parameters is done by a statistical approach (analysis of variance ANOVA combined with the theory of Bootstrap).A characterization was performed using (2D) data of measurement on abrasive surfaces. The influence of the form of discrete and continuous wavelet on the detection on the relevant scale mechanism of the abrasion was tested. We conclude that the wavelet decomposition allows to quantify and localize the scales of abrasion of the machining process for all process parameters. However, the relevance of appropriate scales to characterize abrasion does not depend on the shape of the wavelet.In this work, a new 3D roughness parameter is proposed to quantify the smoothness of a surface, independently of the amplitude and the scanning length units of the surface. The efficiency of this parameter is tested on noisy periodic surfaces with varying degrees of anisotropy. The value of this parameter is between zero (perfect sound) and 100 % (sine perfect surface). It enables us to identify the anisotropy directions of regularity for a given surface

Le stock d’eau, ses variations temporelles et leurs répartitions spatiales, sont des grandes inconnues du cycle hydrologique d’Afrique de l’Ouest. Les trois principaux types de gravimètres actuels (un gravimètre relatif supraconducteur –SG–, un microgravimètre relatif –CG5–, et un gravimètre absolu –FG5–) ont été déployés sur un bassin versant élémentaire de l’observatoire AMMA-CATCH de l’Ouémé supérieur (Bénin), représentatif d’un milieu cultivé de la zone soudanienne. Conjointement, un important dispositif de suivi hydrologique (piézométrie, sonde à neutrons) a été mis en place, ainsi que des mesures de géophysique de proche surface (mesures électriques, électromagnétiques et résonance magnétique des protons –RMP–) pour la caractérisation des aquifères. Ces travaux font état des apports de la gravimétrie pour 1) la détermination de la porosité de drainage en zone de socle hétérogène, 2) le suivi des variations de stock d’eau intégrées à l’échelle d’appréhension des gravimètres (parcelle de 100m de côté) et 3) l’identification de processus hydrologiques liés à la redistribution interne au sein du bassin versant et à la genèse des écoulements
Water storage, together with its temporal and spatial variations, are major unknowns of the West African hydrological cycle. The three main gravimeter types (a relative superconducting gravimeter –SG–, a relative microgravimeter –CG5–, and an absolute gravimeter –FG5–) have been set up on an elementary catchment of the Upper Ouémé (Benin) AMMA-CATCH observatory, considered as a typical cultivated environment of the Sudanian area. Hydrological monitoring (water table, neutron probe) and near surface geophysics (electrical, electromagnetic, magnetic resonance soundings –MRS–) have been developed jointly, for aquifer characterization. This work states on the inputs provided by gravimetry on 1) specific yield determination in heterogeneous basement area, 2) the monitoring of water storage changes at the gravity measurements scale (100m square plot) and 3) hydrological processes identification, linked to internal catchment redistribution and streamflow generation

Cette thèse traite du problème de l'écoulement d'un fluide dans des interfaces étroites entre des solides en contact sous un chargement normal, ce qui est important pour de nombreuses applications en tribologie, ingénierie et géophysique. Le traitement de ce problème nécessite de prévoir un couplage entre la mécanique des fluides et celle des solides. Les contraintes liées à la présence du contact, ainsi que les caractéristiques complexes de la géométrie de surface rajoutent un niveau de complexité significatif. Dans cette thèse, un solveur monolithique par éléments finis permettant la gestion du contact frottant, des écoulements visqueux incompressibles et du transfert des efforts induits par le fluide sur le solide est développé. De plus, la possibilité que le fluide se retrouve piégé dans des cavités délimitées par des zones de contact est prise en compte par l'élaboration d'un nouvel élément dit "de fluide piégé", qui utilise une loi de comportement compressible non linéaire. Le code résultant de cette méthode comprend des algorithmes d’analyse d’image permettant de distinguer les zones de contact, d’écoulement de fluide et de fluide piégé. En outre, le code convient aux approches de couplage uni- et bidirectionnel. Le cadre développé a été appliqué dans un premier temps à l'étude d'un fluide piégé entre un solide déformable présentant une surface de contact ondulée et un plan rigide. Pour un système soumis à une charge externe croissante, nous avons examiné l'évolution de la surface de contact et du coefficient de frottement global en fonction des propriétés du fluide et du solide, ainsi que de la pente du profil de surface. Nous avons ensuite étudié l’écoulement d’un fluide entre un plan rigide et un solide déformable avec une géométrie modèle ou une surface rugueuse. Nous avons obtenu une solution analytique approchée qui gouverne le flux de fluide à travers une interface de contact ondulée, et cette dernière a été comparée à nos résultats numériques. Enfin, nous avons montré pour un intervalle de paramètres physiquement pertinents, que le couplage unidirectionnel sous-estime, par rapport à une approche bidirectionnelle, la perméabilité de l’interface ainsi que la charge externe critique nécessaire à la fermeture de l’interface. Une loi phénoménologique raffinée de perméabilité macroscopique des interfaces de contact rugueuses a été proposée. Enfin, le cadre développé a été utilisé pour calculer l'évolution de la fuite de fluide à travers une interface de contact métal sur saphir en utilisant un comportement matériau élasto-plastique et des mesures réelles de la rugosité de surface
This thesis deals with the problem of a thin fluid flow in narrow interfaces between contacting solids subject to a normal loading, which is relevant for a range of tribological and engineering applications, as well as for geophysical sciences. The treatment of this problem requires coupling between fluid and solid mechanics, further complicated by contact constraints and potentially complex geometrical features of contacting surfaces. In this thesis a monolithic finite-element framework for handling frictional contact, thin incompressible viscous flow and transfer of fluid-induced tractions to the solid is developed. Additionally, we considered fluid entrapment in "pools" delimited by contact patches and formulated a novel trapped-fluid element using a non-linear compressible constitutive law. This computational framework makes use of image analysis algorithms to distinguish between contact, fluid flow and trapped fluid zones. The constructed framework is suitable for both one- and two-way coupling approaches. First, the developed framework was applied to a study of a fluid trapped between a deformable solid with a wavy surface and a rigid flat. We showed how the contact area and the global coefficient of friction evolve under increasing external load, depending on fluid and solid properties and on the slope of the surface profile. Next, we studied a thin fluid flow between a rigid flat and a deformable solid with a model geometry or random surface roughness. An approximate analytical solution for the fluid flow across a wavy contact interface was derived and compared with numerical results. We showed that for a range of physically relevant parameters, one-way coupling underestimates the interface permeability and the critical external load needed to seal the interface, compared to the two-way approach. A refined non-local phenomenological law for macroscopic permeability of rough contact interfaces was proposed. Finally, the developed framework was used to calculate the evolution of the fluid leakage through a metal-to-sapphire contact interface using an elasto-plastic material behaviour and real measurements of surface roughness

Dans les grands systèmes de réfrigération, l’utilisation de réfrigérants comme lubrifiants, à la place des huiles, simplifie la conception du système, l’allège et réduit son impact environnemental. La très faible viscosité du réfrigérant conduit à une épaisseur de film séparant les deux surfaces comparable à leur rugosité. Néanmoins, des travaux expérimentaux avec le réfrigérant R1233zd suggèrent que la lubrification est possible dans ces conditions grâce à la formation d’une couche adsorbée sur la surface d’oxyde de fer. Les analyses expérimentales in situ dans le contact sont très difficiles à cause des conditions de fort confinement et haute pression. C’est pourquoi une approche numérique multi-échelles est développée, afin d’étudier l’impact des réactions physico-chimiques à l’interface réfrigérant-surface sur les propriétés de lubrification. La théorie de la fonctionnelle de la densité est utilisée pour quantifier au niveau quantique, l’adsorption d’une molécule de réfrigérant sur une surface d’oxyde de fer. Des énergies de liaison allant de -0.92 eV à -0.22 eV sont observées et reliées à différents cas d’adsorption. Ces résultats sont exploités pour paramétrer un champ de forces interfacial, qui prédit des structures moléculaires à l’interface, différentes de celles obtenues avec des potentiels basés sur les règles de mélange classiques. Des simulations de dynamique moléculaire utilisant ce champ de forces paramétré confirment l’existence d’une couche fortement adsorbée de R1233zd sur une surface d’oxyde de fer. Avec des surfaces atomiques lisses et seulement 2 nm d’épaisseur de film de réfrigérant, les couches adsorbées résistent à des pressions allant jusqu’à 4 GPa et des vitesses de cisaillement atteignant 100 m/s. Une valeur minimale de 5 molécules de réfrigérant par nm² est nécessaire à la formation de deux couches adsorbées à 0.5 GPa. De plus, des simulations en cisaillement avec une surface rugueuse prédisent une rupture totale du film à environ 13 GPa
In large refrigeration systems, using the refrigerant as lubricant instead of oil can help to simplify the design, lighten the systems, and reduce their environmental impact. However, the very low viscosity of refrigerants leads to ultra-thin films separating the surfaces, with a thickness comparable to surface roughness. Nevertheless, experiments with the R1233zd refrigerant suggests that lubrication is still possible in that situation thanks to an adsorbed layer formed on iron oxide surfaces. Experimental in situ analysis area is very difficult because of high confinement and high pressure. That is why a multiscale numerical approach is developed here, to explore the impact of chemical reactions and physical processes at the refrigerant-surface interface on large-scale lubrication properties. Density functional theory is used to quantify the adsorption of a refrigerant molecule on an iron oxide surface at the quantum level. Binding energies ranging from -0.92 eV to -0.22 eV are measured and related to different adsorption cases. These results are then used to parametrize an interfacial force field, whose predictions of interfacial molecular structure differs from those obtained using potentials based on standard mixing rules. Large-scale molecular dynamics simulations involving this parametrized force field confirm the existence of a strongly adsorbed layer of R1233zd molecules on iron oxide surface. With atomically smooth surfaces, and a refrigerant film thickness as small as 2 nm, the adsorbed refrigerant layers resists pressures as high as 4 GPa and sliding velocities as high as 100 m/s. A minimum value of 5 refrigerant molecules per nm² is necessary to the formation of two adsorbed layers at 0.5 GPa. Moreover, sliding simulations with a rough surface reveal total film breakdown for ca. 13 GPa

Depuis quelques années, l'évolution des techniques d'acquisition a entraîné une généralisation de l'utilisation d'objets 3D très dense, représentés par des nuages de points de plusieurs millions de sommets. Au vu de la complexité de ces données, il est souvent nécessaire de les analyser pour en extraire les structures les plus pertinentes, potentiellement définies à plusieurs échelles. Parmi les nombreuses méthodes traditionnellement utilisées pour analyser des signaux numériques, l'analyse dite scale-space est aujourd'hui un standard pour l'étude des courbes et des images. Cependant, son adaptation aux données 3D pose des problèmes d'instabilité et nécessite une information de connectivité, qui n'est pas directement définie dans les cas des nuages de points. Dans cette thèse, nous présentons une suite d'outils mathématiques pour l'analyse des objets 3D, sous le nom de Growing Least Squares (GLS). Nous proposons de représenter la géométrie décrite par un nuage de points via une primitive du second ordre ajustée par une minimisation aux moindres carrés, et cela à pour plusieurs échelles. Cette description est ensuite derivée analytiquement pour extraire de manière continue les structures les plus pertinentes à la fois en espace et en échelle. Nous montrons par plusieurs exemples et comparaisons que cette représentation et les outils associés définissent une solution efficace pour l'analyse des nuages de points à plusieurs échelles. Un défi intéressant est l'analyse d'objets 3D acquis dans le cadre de l'étude du patrimoine culturel. Dans cette thèse, nous nous étudions les données générées par l'acquisition des fragments des statues entourant par le passé le Phare d'Alexandrie, Septième Merveille du Monde. Plus précisément, nous nous intéressons au réassemblage d'objets fracturés en peu de fragments (une dizaine), mais avec de nombreuses parties manquantes ou fortement dégradées par l'action du temps. Nous proposons un formalisme pour la conception de systèmes d'assemblage virtuel semi-automatiques, permettant de combiner à la fois les connaissances des archéologues et la précision des algorithmes d'assemblage. Nous présentons deux systèmes basés sur cette conception, et nous montrons leur efficacité dans des cas concrets.

Dans une première partie, nous présentons des équations de Saint-Venant. Sur le modèle proprement dit, nous remarquons tout d'abord que, suivant le lien entre la viscosité et le rapport des échelles caractéristiques, il est indispensable de conserver l'expression complète de la force de Coriolis : nous obtenons ainsi un nouveau modèle, avec un "effet cosinus". Nous montrons alors que les preuves d'existence de solutions faibles peuvent être adaptées à ce nouveau système. Des simulations numériques de certaines ondes soulignent l'importance de ce terme. Nous étudions ensuite l'influence des conditions limites (surface, fond) sur des modèles de type Saint-Venant. Nous présentons également des modèles obtenus en utilisant des échelles multiples en espace et en temps. Enfin, nous analysons théoriquement et numériquement un nouveau modèle de sédimentation puis nous donnons certains résultats pour les fluides visco-plastiques. Dans une deuxième partie, nous nous intéressons aux équations limites que sont les équations quasi-géostrophiques (QG) et les équations des lacs. L'étude numérique des équations QG 2d nous permet de voir le rôle de l'effet cosinus de la force de Coriolis. En fonction de la topographie considérée, nous montrons que celui-ci peut être non négligeable. Toujours sur les équations QG, nous donnons un schéma, basé sur des développements asymptotiques, qui permet de bien capter la couche limite mais aussi d'ajouter le terme de topographie à la solution obtenue avec fond plat, sans refaire tous les calculs. Enfin, nous expliquons l'obtention des équations des lacs avec effet cosinus, et nous prouvons que les propriétés d'existence de solutions restent valables
In a first part, we present some Shallow Water equations. About the actual model, we firstly remark that, depending on the link between the viscosity and the aspect ratio, keeping the complete Coriolis force expression is essential: this gives a new model, with a so-called "cosine effect". We then show that the proofs of existence of weak solutions can be adapted to this new system. Numerical simulations of some waves underline the fact that this term is of importance. Next we study the influence of the limit conditions (surface, bottom) on Shallow-Water type models. We also present some models obtained using multiple scales in space and time. Finally we analyze a new model of sedimentation from a theoretical and numerical point of view and then we give some results for visco-plastic fluids. In a second part, we are interested in the limit equation, namely the Quasi-Geostrophic (QG) equations and the lake equations. The numerical study of the 2d QG equations enables us to emphasize the role of the cosine effect from the Coriolis force. Depending on the topography we consider, we show that this effect can turn out to be not negligible. Still about the QG equations, we give a numerical scheme, based on asymptotic developments, which capture the boundary layer well and also give the opportunity to add a topography term to the solution for a flat bottom, without re-computing everything. Lastly we explain how to get the lake equations with cosine effect and we prove that the properties of existence of solutions to such equations are still valid

L'érosion des sols est pourtant connue comme l'un des problèmes environnementaux les plus préoccupants au monde. L'érosion des sols est particulièrement et de plus en plus entraînée par les activités anthropiques dans le cadre du changement climatique. En RDP lao, un pays tropical, l'érosion des sols est due de manière significative à une gestion inappropriée des terres sur les terrains en pente. Le Houay Pano, un bassin versant cultivé du nord de la RDP lao, est exposé à l'érosion des sols, en particulier après la conversion de la culture itinérante en plantation de teck. La mauvaise gestion des terres en défrichant le sous couvert végétal sous la plantation de teck est considérée comme une cause sous-jacente du coefficient de ruissellement (Rc) plus élevé et de l'érosion des sols. Certaines mesures d'atténuation telles que le sous couvert végétal et la végétation rivulaire sont suggérées pour atténuer l'érosion des sols. Cependant, la mesure d'atténuation de l'érosion des sols et l'effet de la gestion de l'utilisation des terres sur le ruissellement de surface (SR) et la perte en sols/rendement de sédiments (Sl) à plusieurs échelles dans la plantation de teck ne sont pas entièrement évalués. Dans ce contexte, nous émettons l'hypothèse que le sous couvert végétal et l'herbe rivulaire atténuent l'érosion du sol dans la plantation d'arbres à teck et que les plantations d'arbres à teck ont ??des impacts sur SR et Sl entraînés par des processus dominants (l'érosion en nappe, l'érosion linéaire et le dépôt de sédiment) sur diverses échelles spatiales. Par conséquent, les objectifs fixés pour ce travail sont : (1) d'évaluer l'effet de la gestion du sous couvert végétal sur le ruissellement de surface et la perte en sols dans la plantation de teck à l'échelle de la micro-parcelle ; (2) d'évaluer la capacité des zones tampons d'herbes rivulaires à atténuer SR et Sl, et d'évaluer leur efficacité de piégeage de l'eau et des sédiments dans les plantations de teck sans sous couvert végétal à l'échelle du versant ; et (3) d'évaluer l'effet de la plantation de teck sur SR et Sl à diverses échelles spatiales (échelles de micro-parcelle, de versant incluant micro-bassin versant, et de bassin versant) dans un bassin versant tropical montagneux à utilisations mixtes de terre. Dans cette étude, le village de Ban Kokngew et le bassin versant d'Houay Pano ont été sélectionnés comme zones d'étude expérimentale pendant la saison des pluies. Des micro-parcelles, des pièges Gerlach et des déversoirs ont été utilisés pour estimer SR et Sl à chaque échelle. Nous avons suivi le modèle TEST développé pour l'érosion en nappe, qui nécessite quelques paramètres, pour évaluer Sl sur la micro-parcelle et le mettre en hautes échelles spatiales pour prédire Sl à l'échelle du versant et du bassin versant. Dans une étude réalisée en 2017 dans les plantations de teck de Ban Kokngew à l'échelle micro-parcelle, nous avons montré que Rc et Sl (23%, 381 Mg·km-2, respectivement) sous teck avec sous couvert végétal étaient inférieurs à ceux sous teck sans sous couvert végétal (60 % et 5455 Mg·km-2, respectivement). Par conséquent, l'atténuation de l'érosion des sols par le maintien du sous couvert végétal sous la plantation de teck réduit Sl de 14 fois. [...]
Soil erosion is yet known as one of the most concerning problems of the environment in the world. Soil erosion is particularly and increasingly driven by anthropogenic activities under the changing climate. In Lao PDR, a tropical country, soil erosion is significantly due to inappropriate land management on the sloping land. The Houay Pano, a cultivated catchment of the northern Lao PDR, is prone to soil erosion, particularly after the conversion from shifting cultivation to teak tree plantation. Land mismanagement by clearing the understory under the teak tree plantation is considered as an underlying cause of higher runoff coefficient (Rc) and soil erosion. Some mitigations such as understory and riparian vegetation are suggested for alleviating soil erosion. However, the mitigation measure of soil erosion and the effect of land use management on surface runoff (SR) and soil loss/sediment yield (Sl) on multiple scales in the teak tree plantation are not fully assessed. In this context, we hypothesize that understory and riparian grass mitigate the soil erosion in the teak tree plantation and that teak tree plantation impacts on SR and Sl driven by dominant processes (inter rill erosion, linear erosion, and deposition) on various spatial scales. Therefore, the objectives set out for this work are: (1) to assess the effect of understory management on SR and Sl in the teak tree plantation on the microplot scale; (2) to assess the ability of riparian grass buffers to mitigate SR and Sl, and to assess their water and sediment trapping efficiencies in the teak tree plantations with no understory on the hillslope scale; and (3) to assess the effect of teak tree plantation on SR and Sl on various spatial scales (microplot, hillslope including micro-catchment, and catchment scales) in a mixed land uses mountainous tropical catchment. In this study, Ban Kokngew village and Houay Pano catchment were selected as experimental study areas during the rainy season. Microplots, Gerlach traps, and weirs were used to estimate SR and Sl on each scale. We followed the TEST model developed for inter rill erosion, which requires a few parameters, to assess Sl on the microplot and upscale it to predict Sl on the hillslope and catchment scale. In a study performed in 2017 in the teak tree plantations of Ban Kokngew on the microplot scale, we showed that Rc and Sl (23%, 381 Mg·km-2, respectively) under teak tree with understory were less than those under teak tree with no understory (60% and 5455 Mg·km-2, respectively). Hence, soil erosion mitigation by keeping the understory under teak tree plantation reduces Sl by 14 times. In a study performed in 2014 in the teak tree plantations of Houay Pano on both the microplot and the hillslope scales, we showed that leaving the riparian grass buffer of at least 6 m could limit SR and Sl discharging downstream during small storms (24-hour rainfall < 54.8 mm) with the trapping efficiency up to 88%. Lastly, in a study performed in 2014 in the teak tree plantations of Houay Pano on various scales, we showed that SR and Sl were significantly higher (p-value < 0.05) in the teak-dominated micro-catchment than in the mixed-land-use micro-catchment. SR and Sl decreased from the microplot (122 - 196 mm, 275 - 1065 Mg·km-2, respectively) to the micro-catchment (24 - 188 mm, 95 - 3635 Mg·km-2, respectively) and catchment scale (33 mm, 236 Mg·km-2, respectively), except that Sl in teak tree plantation increased from the microplot (1065 Mg·km-2) to the micro-catchment scale (3635 Mg·km-2). [...]

Dans une première partie, nous présentons des équations de Saint-Venant. Sur le modèle proprement dit, nous remarquons tout d'abord que, suivant le lien entre la viscosité et le rapport des échelles caractéristiques, il est indispensable de conserver l'expression complète de la force de Coriolis : nous obtenons ainsi un nouveau modèle, avec un "effet cosinus". Nous montrons alors que les preuves d'existence de solutions faibles peuvent être adaptées à ce nouveau système. Des simulations numériques de certaines ondes soulignent l'importance de ce terme. Nous étudions ensuite l'influence des conditions limites (surface, fond) sur des modèles de type Saint-Venant. Nous présentons également des modèles obtenus en utilisant des échelles multiples en espace et en temps. Enfin, nous analysons théoriquement et numériquement un nouveau modèle de sédimentation puis nous donnons certains résultats pour les fluides visco-plastiques.
Dans une deuxième partie, nous nous intéressons aux équations limites que sont les équations quasi-géostrophiques (QG) et les équations des lacs. L'étude numérique des équations QG 2d nous permet de voir le rôle de l'effet cosinus de la force de Coriolis. En fonction de la topographie considérée, nous montrons que celui-ci peut être non négligeable. Toujours sur les équations QG, nous donnons un schéma, basé sur des développements asymptotiques, qui permet de bien capter la couche limite mais aussi d'ajouter le terme de topographie à la solution obtenue avec fond plat, sans refaire tous les calculs. Enfin, nous expliquons l'obtention des équations des lacs avec effet cosinus, et nous prouvons que les propriétés d'existence de solutions restent valables.

La première partie de l'étude porte sur la construction de bases multi-échelles inspirées de l'analyse des ondelettes et leurs applications aux problèmes d'inversion en sismique. Plus précisément nous comparons les résultats d'identification des paramètres sismiques avec plusieurs bases multi-échelles. Dans la seconde partie l'identification des paramètres sismiques dans un milieu stratifie est remplacée par celle d'un operateur impédance équivalente. Cet operateur est approche dans une classe de fraction rationnelle. Nous présentons les résultats numériques ainsi obtenus et les difficultés rencontrées

La paramétrisation du transfert de chaleur par télédétection, basée sur le schéma SEBS, s'est déjà avérée très adaptée pour l'estimation de l'évapotranspiration (ET) sur des surfaces naturelles homogènes. Cependant, l'utilisation d'une telle méthode pour des paysages hétérogènes (e.g. régions semi-arides ou surfaces agricoles) est plus délicate, puisque le principe de la théorie de la similarité est compromis par la présence de différentes sources de chaleur et de hauteurs variées. Dans un premier temps, cette thèse a pour objectif de proposer et d'évaluer différents modèles basés sur la géométrie de la végétation qui permettent d'estimer la longueur de rugosité pour le transfert de quantité de mouvement à la surface (z0m), cette dernière étant un paramètre clé dans la caractérisation du transfert de chaleur. En revanche, une telle investigation ne peut être menée qu'à une petite échelle et à l'aide de données de télédétection très haute résolution permettant ainsi une description très détaillée de la surface. Ensuite, le second aspect de ce travail est de caractériser le transfert de chaleur dans le cas d'études régionales. Puis, la capacité de SEBS à estimer les flux de chaleur turbulents à de grandes échelles spatiales et temporelles sera évaluée. Pour ce faire, l’approche multi-échelle de SEBS (MSSEBS) a été implémentée afin de traiter une zone de 2,4 millions km2, incluant le Plateau du Tibet et l’amont des principaux fleuves d’Asie du sud-est. La combinaison de données horaires de température de surface FY-2 avec un rayonnement net journalier et des paramètres de surface avancés, permet de produire une série temporelle d’ET sur le Plateau du Tibet pour la période 2008-2010, et à une fréquence journalière
The parameterization of heat transfer by remote sensing, and based on SEBS scheme for turbulent heat fluxes retrieval, already proved to be very convenient for estimating evapotranspiration (ET) over homogeneous land surfaces. However, the use of such a method over heterogeneous landscapes (e.g. semi-arid regions or agricultural land) becomes more difficult, since the principle of similarity theory is compromised by the presence of different heat sources with various heights. This thesis aims at first to propose and evaluate some models based on vegetation geometry for retrieving the surface roughness length for momentum transfer (z0m), which is a key parameter in the characterization of heat transfer. Such an investigation can only be led at a small scale with very-high resolution remote sensing data, for a precise description of the land surface. Therefore, the second aspect of this work is to determine how to address the characterization of heat transfer for regional studies. Then, the reliability of SEBS for estimating turbulent heat fluxes at large spatial and temporal scales has been evaluated. To do so, the Multi-Scale SEBS approach (MSSEBS) has been implemented for a 2.4 million km2 area including the Tibetan Plateau and the headwaters of the major rivers of East and South Asia. The addition of gap-filled hourly FY-2 LST data to advanced daily averaged net radiation and land surface parameters, allows to compute time-series of land surface ET over the Tibetan Plateau during the period 2008-2010, and on a daily basis

Ma thématique de recherche s'intéresse à l'étude de la fatigue des matériaux, et se focalise plus particulièrement sur l'analyse de la plasticité cyclique et des endommagements induits dans les matériaux métalliques. Elle s'inscrit dans une démarche globale conjointe d'analyse expérimentale et de modélisation numérique.
Les travaux présentés se focalisent sur l'étude du comportement d'un acier inoxydable austéno-ferritique et s'articulent suivant quatre thèmes :
- Etude expérimentale du comportement mécanique en fatigue oligocyclique. Un certain nombre de résultats expérimentaux en fatigue oligocyclique sont présentés, en particulier, l'importance de caractériser le comportement mécanique sous sollicitations multiaxiales est soulignée. Deux méthodes sont comparées permettant la mesure des caractéristiques de l'écrouissage : la mesure de surfaces de plasticité en cours d'essai, et l'analyse a posteriori des boucles d'hystérésis.
- Modélisation du comportement mécanique observé. Deux approches ont été utilisées, d'une part une approche phénoménologique qui nécessite des mesures des caractéristiques de l'écrouissage, et d'autre part une approche micromécanique qui s'appuie sur l'analyse de la structure biphasée de l'acier étudié. L'identification de ce deuxième modèle a été réalisée sur une base expérimentale particulièrement réduite sans cependant amoindrir sa capacité de prédiction. Une proposition originale pour la modélisation de l'adoucissement cyclique a été faite à partir d'observations au microscope électronique à transmission.
- Etude de la déformation progressive. A partir de résultats expérimentaux biaxiaux, il est montré que les niveaux de déformation progressive obtenus peuvent s'expliquer par la direction d'écoulement plastique observée, due à la distorsion de la surface de plasticité. Un modèle simulant cette distorsion est alors identifié.
- Analyse et simulation de la plasticité à l'échelle de la microstructure. La plasticité cyclique et l'amorçage en fatigue étant fortement liés aux caractéristiques de la microstructure, une analyse locale a été entreprise. Un dispositif expérimental original a été développé afin d'analyser les mécanismes de déformation in-situ pendant un essai mécanique. Plusieurs moyens d'analyse sont utilisés conjointement (microscopie, champs cinématiques calculés à partir d'images obtenues entre deux instants, cartographie EBSD, profilométrie). Une comparaison est faite entre les champs cinématiques obtenus et un calcul d'agrégat polycristallin reprenant les lois cristallines identifiées précédemment.

La première partie du mémoire concerne l'étude de modèles multi-échelles (ou micro-macro) pour les fluides complexes. L'intérêt de ces modèles est d'éviter d'avoir à postuler des lois de comportements phénoménologiques, en couplant une description macroscopique classique (lois de conservations de la quantité de mouvement et de la masse) sur la vitesse et la pression, à des modèles microscopiques pour l'évolution des microstructures dans le fluide, à l'origine du caractère non-newtonien du fluide. La loi de comportement consiste alors à exprimer le tenseur des contraintes en fonction de la conformation des microstructures. Durant la thèse, nous avons obtenu quelques résultats d'existence et de convergence des méthodes numériques utilisées pour ces modèles. Plus récemment, nous proposons une méthode numérique pour coupler des modèles micro-macro (fins mais chers) avec des modèles macroscopiques (plus grossiers, mais beaucoup plus économiques en temps de calcul). Nous analysons par ailleurs une méthode de résolution des équations aux dérivées partielles en grande dimension qui a été proposée dans le contexte de la résolution numérique des modèles à l'échelle microscopique pour les fluides polymériques.

L'essentiel de nos travaux les plus récents sur le sujet concerne le comportement en temps long de ces modèles, avec un double objectif : théoriquement, la compréhension des modèles physiques passe souvent par l'étude de la stabilité des solutions stationnaires, et de la vitesse de convergence vers l'équilibre ; numériquement, la stabilité des schémas en temps long est cruciale, car on utilise typiquement des simulations instationnaires en temps long pour calculer des solutions stationnaires. Nous montrons comment analyser le comportement des modèles micro-macro en temps long, en utilisant des méthodes d'entropie. Cette étude a ensuite permis de comprendre le comportement en temps long de modèles macroscopiques standard (type Oldroyd-B), et de préciser sous quelles conditions les schémas numériques vérifient des propriétés similaires de stabilité en temps long.


La seconde partie du mémoire résume des travaux en simulation moléculaire, à l'échelle quantique, ou à l'échelle de la dynamique moléculaire classique. A l'échelle quantique, nous nous intéressons aux méthodes Quantum Monte Carlo. Il s'agit de méthodes numériques probabilistes pour calculer l'état fondamental d'une molécule (plus petite valeur propre de l'opérateur de Schrödinger à N corps). Essentiellement, il s'agit de donner une interprétation probabiliste du problème par des formules de Feynman-Kac, afin de pouvoir appliquer des méthodes de Monte Carlo (bien adaptées pour des problèmes de ce type, en grande dimension). Nous proposons tout d'abord une étude théorique de la méthode Diffusion Monte Carlo, et notamment
d'un biais introduit par l'interprétation probabiliste (appelé fixed node approximation). Nous analysons ensuite les méthodes numériques utilisées en Diffusion Monte Carlo, et proposons une nouvelle stratégie pour améliorer l'échantillonnage des méthodes Variational Monte Carlo.

En dynamique moléculaire, nous étudions des méthodes numériques pour le calcul de différences d'énergie libre. Les modèles consistent à décrire l'état d'un système par la position (et éventuellement la vitesse) de particules (typiquement les positions des noyaux dans un système moléculaire), qui interagissent au travers d'un potentiel (qui idéalement proviendrait d'un calcul de mécanique quantique pour déterminer l'état fondamental des électrons pour une position donnée des noyaux). L'objectif est de calculer des moyennes par rapport à la mesure de Boltzmann-Gibbs associée à ce potentiel (moyennes dans l'ensemble canonique). Mathématiquement, il s'agit d'un problème d'échantillonnage de mesures métastables (ou multi-modales), en très grande dimension. La particularité de la dynamique moléculaire est que, bien souvent, on a quelques informations sur les "directions de métastabilité" au travers de coordonnées de réaction. En utilisant cette donnée, de nombreuses méthodes ont été proposées pour permettre l'échantillonnage de la mesure de Boltzmann-Gibbs. Dans une série de travaux, nous avons analysé les méthodes basées sur des équations différentielles stochastiques avec contraintes (dont les solutions vivent sur des sous-variétés définies comme des lignes de niveaux de la coordonnée de réaction). Il s'agit en fait d'analyser des méthodes d'échantillonnage de mesures définies sur des sous-variétés de grande dimension. Plus récemment, nous avons étudié des méthodes adaptatives qui ont été proposées pour se débarrasser des métastabilités. Mathématiquement, il s'agit de méthodes d'échantillonnage préférentiel, avec une fonction d'importance qui est calculée au cours de la simulation de manière adaptative. Nous avons étudié la vitesse de convergence vers la mesure d'équilibre pour ces méthodes adaptatives, en utilisant des méthodes d'entropie. Nous avons proposé de nouvelles méthodes numériques à la communauté appliquée pour utiliser au mieux ces idées.


La troisième partie du mémoire résume des travaux issus d'une collaboration avec l'entreprise Rio Tinto (anciennement Pechiney puis Alcan), leader mondial pour la technologie des cuves d'électrolyse de l'aluminium. Cette collaboration a été entamée il y a plusieurs années par C. Le Bris, et notamment au travers de la thèse de J-F. Gerbeau. Mathématiquement, il s'agit d'analyser et de discrétiser les équations de la magnétohydrodynamique pour deux fluides incompressibles non miscibles, séparés par une interface libre. Nous expliquons le contexte industriel et la modélisation, nous résumons la méthode numérique adoptée (méthode Arbitrary Lagrangian Eulerian) et donnons quelques propriétés satisfaites par le schéma (stabilité, conservation de la masse). Nous montrons ensuite comment ce modèle permet d'étudier un phénomène (potentiellement déstabilisant) observé dans les cuves d'électrolyse : le rolling. Ces résultats ont été pour la plupart obtenus durant la thèse.

Plus récemment, dans le prolongement de l'étude industrielle, nous nous sommes intéressés à un problème de modélisation fondamentale pour les écoulements à surface (ou interface) libre: le mouvement de la ligne de contact (i.e. le bord de la surface libre qui glisse le long de la paroi). En résumé, nos travaux consistent essentiellement en deux contributions: (i) une compréhension variationnelle d'une condition aux limites permettant de modéliser correctement le mouvement de la ligne de contact (Generalized Navier Boundary Condition), et son implémentation dans un schéma Arbitrary Lagrangian Eulerian, (ii) une analyse de la stabilité du schéma obtenu.

La Corrosion Haute Température (HTC), en environnements divers et sévères, d'alliages métalliques toujours plus élaborés en termes de composition et micro-structure, est un sujet industriel et scientifique très complexe. La PhotoElectroChimie (PEC) est une technique de choix pour caractériser les propriétés physico-chimiques et électroniques des couches d'oxydation très hétérogènes formées en HTC. Sur des exemples d'alliages modèles mais représentatifs de la réalité industrielle (aciers duplex, base-Nickel 690), ce travail présente le développement et la validation d'un dispositif expérimental permettant d'appliquer pour la première fois tout l'arsenal des techniques PEC à l'échelle mésoscopique (typiquement 30 µm), ainsi que la validation d'une modélisation originale développée au SIMaP des spectres de photocourants en énergie, qui permet de les décrire et ajuster finement et d'en extraire notamment avec précision les gaps des oxydes semiconducteurs présents dans la couche thermique
High Temperature Corrosion (HTC), in various and severe atmospheres, of continually more elaborated (composition, micro–structure) metallic alloys, is a rather complex industrial and scientific topic. PhotoElectroChemistry (PEC) acquired a special place in the characterization of physico–chemical and electronic properties of the highly heterogeneous oxidation layers formed in HTC. Through studies of model but industrially representative samples (duplex stainless steel, Ni–base alloy 690), this work presents the development and validation of an experimental set–up allowing for the first time to use the whole set of PEC techniques at the mesoscopic level (typically 30 µm), as well as the validation of an original model of photocurrent energy spectra, developed at SIMaP, allowing to well describe, and accurately fit the latter spectra, and thus yielding, notably, precise bandgap values for the semiconducting components of the thermal scale

Jean M.F. Martins, CR1 CNRS, LTHE (UMR 5564), CNRS, INPG, IRD, Univ. J. Fourier, Grenoble I Résumé : L'ensemble de mes travaux est basé sur une approche pluri et interdisciplinaire développée dans le but d'améliorer nos connaissances et la compréhension du comportement biogéochimique (sorption, biotransformation et transfert) et l'impact des micropolluants dans l'environnement, afin de mieux prédire l'exposition et les impacts sur notre santé ou celle de nos écosystèmes. Au moment ou des polluants émergeants sont détectés dans les eaux et les sols dès qu'on se donne les moyens pour les détecter, et alors que la directive européenne REACH se met en place, les approches de recherche monodisciplinaires montrent de plus en plus clairement leurs limites et font ainsi ressortir le caractère indispensable de la pluridisciplinarité dans l'approche de la problématique de la biogéochimie des polluants. C'est dans ce contexte que j'ai essayé, à chaque fois que cela était nécessaire et possible, de développer une approche interdisciplinaire en rupture avec les approches classiques, en abordant notamment les phénomènes couplés ou cinétiques, dans le but d'évaluer simultanément et de modéliser le comportement biogéochimique (sorption, biotransformation et transferts) et l'impact des micropolluants dans les sols à plusieurs échelles et dans différents contextes biogéochimiques. Les approches ainsi développées ont permis de caractériser les mécanismes de sorption des polluants dans le sol et dans ses différents compartiments, minéral, organique et biotique et d'établir l'effet de ces processus sur le transfert réactifs de ces polluants en combinant des approches expérimentales et théoriques en conditions dynamiques et en faisant appel à des dispositifs technologiques spécifiquement développés. La première partie de ce mémoire concerne l'étude des processus de sorption des micropolluants organiques et métalliques dans les sols, et des divers facteurs qui les contrôlent, notamment la spéciation chimique. Cette spéciation contrôle directement l'ambiance chimique dans les sols, qui correspond à un équilibre de la chimie et de la distribution des micropolluants entre les différents compartiments du sol. Ainsi les interactions des micropolluants avec les constituants abiotiques (minéraux et matière organique) et biotiques (microorganismes) des sols et leur modélisation ont été étudiées par des approches combinées de dynamique des systèmes (colonnes), chimie de solution (batch) et spectroscopiques ou microscopiques (spectroscopie EXAFS, MEB- et MET-EDX). Le premier exemple abordé a concerné la réactivité de composés organiques (dinitrophénols) en mono- et multi-contaminations. Pour cette dernière condition, il a été montré de manière originale l'existence de processus complexes et prépondérants de compétition pour les sites de sorption dans les sols et de synergie pour la toxicité, processus largement négligés dans les études classiques mono-polluants, qui ne sont que rarement représentatives de la réalité du terrain. Ainsi, ces phénomènes induisent des modifications importantes du comportement global de ces composés en termes de mobilité (cinétiques de transfert de masse), d'impact et de dégradabilité, rendant également leur prédiction très délicate car nécessitant la prise en compte des phénomènes de synergie et de compétition précités. La réactivité des micropolluants avec une des fractions les plus réactives des sols, la biomasse vivante, a également été étudiée au travers de trois modèles bactériens contrastés (Cupriavidus metallidurans CH34, Escherichia coli, Pseudomonas putida) par une approche combinée de spectroscopie X et de chimie de solution. Ainsi nous avons contribué à la validation d'un modèle universel de réactivité de surface bactérienne aux micropolluants, notamment métalliques, basé sur la prise en compte de groupements réactifs déjà connus, phosphoesters, carboxyles, ou nouvellement identifiés, sulfhydriles, par une approche d'absorption X et de modélisation thermodynamique. Ce modèle reste toutefois peu satisfaisant pour les structures biologiques complexes telles que les biofilms, qui présentent une réactivité largement sous-estimée en raison de la mise en jeu de substances polymériques extracellulaires (EPS) complexes, particulièrement abondantes dans ces structures. Ce volet constitue une de mes principales perspectives de recherche. Dans la seconde partie de ce mémoire, des approches en conditions statiques et dynamiques et mono et multi-pollutions ont permis d'évaluer la biodégradabilité et la biotransformation de contaminants traces en sols naturels. Leur impact sur le compartiment microbien a été évalué par des méthodes d'écotoxicologie classique (bio-indicateurs, marqueurs microbiologiques globaux) et par des méthodes plus innovantes (biocapteurs, biodiversité) permettant la mise en évidence du contrôle de la toxicité des contaminants par leur biodisponibilité, elle-même largement contrôlée par leur spéciation chimique, leur réactivité et leur diffusion au contact des cellules bactériennes (accessibilité). Ces études ont contribué à montrer le rôle important (prépondérant sous certaines conditions) des microorganismes dans la rétention et dans le transport des micropolluants dans les sols, phénomènes qui ont longtemps été négligés et restent encore aujourd'hui difficiles à évaluer in situ. Cela a constitué la dernière partie de ce travail qui présente une approche de dynamique des systèmes en colonnes de sols, permettant d'identifier et modéliser les conditions biogéochimiques impliquées dans le transfert réactif des polluants dans la zone non saturée des sols. Nous avons ainsi caractérisé les principaux facteurs biotiques et abiotiques contrôlant la mobilité des polluants et des bactéries dans les sols et développé une modélisation intégrant ces facteurs dominants. En plus d'être rapidement mobilisées dans les sols, certaines cellules bactériennes peuvent, sous certaines conditions, accélérer fortement le transport de polluants, comme nous avons pu le montrer avec le pentachlorophénol ou avec le Zn, le Cd ou le Cu et le prédire par une modélisation couplée hydrodynamique et géochimique. L'approche pluridisciplinaire et multi-échelle (moléculaire à macroscopique) développée dans ces travaux, fournit des connaissances nouvelles (i) pour la compréhension des cycles biogéochimiques des micropolluants dans les sols, et (ii) pour l'évaluation des risques liés à ces pollutions (notamment les multi-pollutions) et le développement de techniques de bioremédiation de sols et d'eaux contaminés. Mots clefs : Micropolluants, biotransformations, spéciation, écotoxicité, biodisponibilité, accessibilité, réactivité de surface, transfert réactif, bioremédiation, couplages, modélisation.

Ce texte de synthèse a pour but de présenter l'´évolution de mes recherches postèrieures à ma thèse. Ce travail s'articule autour de plusieurs axes de recherche dans le cadre des équations aux dérivées partielles non linéaires et en particulier des lois de conservation.
Il s'inscrit dans l'étude des problèmes hyperboliques, des problème mixtes et des équations cinétiques. Les domaines d'application sont la mécanique des fluides ou du solide, la propagation de composants chimiques, l'électromagnétisme, l'optique.
Mon activité concerne d'abord la modélisation de phénomènes physiques ou chimiques sous forme d'équations aux dérivées partielles non linéaires telles que les équations de Bloch, Korteweg, Navier-Stokes, Saint-Venant, puis vient l'étude mathématique de ces équations à travers les
problèmes d'existence, d'unicité, de régularité avec éventuellement la mise au point de méthodes numériques de résolution.
Ce document est divisé en une introduction générale et trois chapitres qui concernent respectivement les systèmes hyperboliques avec conditions aux limites et la chromatographie, les problèmes d'analyse asymptotique et enfin les méthodes cinétiques.
Dans chaque partie, un historique et une présentation des différents résultats mathématiques sont faits et quelques problèmes ouverts sont donnés.