L'étude présentée ici entre dans la problématique générale des arcs électriques intervenant dans des applications industrielles telles que le soudage, le découpage, le traitement des déchets. Ce travail constitue une première phase de modélisation de cette étude générale. Le plasma est décrit par un ensemble d'équations de conservation de fluide et de l'électromagnétisme, complétés par des propriétés thermodynamiques et des coefficients de transport appropriés, en formant un système d'équations non linéaires fortement couplées. Ces équations sont écrites en supposant l'équilibre thermodynamique local, une symétrie cylindrique et un écoulement laminaire stationnaire. Ce système d'équations est résolu à l'aide du logiciel commercial FLUENT de type CFD fondé sur l'approche des volumes finis. Pour pouvoir utiliser la partie solveur nous avons résolu notre modèle en utilisant les routines UDF Users-Defined-Function. Dans une première partie, nous présentons la validation du modèle à deux dimensions et à 100 A dans l'argon par des résultats de la littérature. Cette comparaison laisse apparaître un accord satisfaisant sur les profils de température dans la colonne de plasma et des différences dans les zones proches des électrodes dues aux conditions aux limites. Dans une deuxième partie, nous présentons une étude expérimentale, à l'issue de laquelle on constate que les profils de température expérimentaux sont en accord avec ceux du modèle dans la zone de colonne positive.

Au cours du dernier siècle, plusieurs séismes majeurs ont affecté l'Himalaya. Cependant, la taille maximale de ces événements et la probabilité d'occurrence de méga-séismes avec des magnitudes proches de 9 sont toujours matière à débat. L'étude de la segmentation de l'arc Himalayen est donc primordiale afin de comprendre les mécanismes qui contrôlent ces séismes ainsi que leur extension spatiale. La compréhension du cycle sismique en Himalaya est aujourd'hui essentiellement basée sur des études menées au Népal central, ce qui limite notre connaissance de son fonctionnement tri-dimensionnel. Ce travail de thèse permet d'apporter de nouvelles contraintes sur les variations latérales de la déformation crustale dans la zone comprise entre l'ouest Népal (~80°E) et l'est Bhoutan (~92°E). La première partie de cette thèse est consacrée à l'étude des variations latérales de la structure lithosphérique. Quatre campagnes gravimétriques ont été réalisées entre 2010 et 2012. Elles ont permis d'étendre le réseau gravimétrique Népalais jusqu'à l'ouest Népal et d'établir le premier réseau gravimétrique Bhoutanais. Ce nouveau jeu de données, combiné avec les données sismologiques disponibles, permet de contraindre des modèles thermo-mécaniques de la flexure de la plaque Indienne sous l'Himalaya. Les résultats suggèrent qu'il n'existe pas de variations latérales majeures du comportement mécanique de la lithosphère Indienne entre le centre et l'est Népal contrairement au Bhoutan où une rigidité flexurale plus faible est nécessaire pour expliquer les données.La deuxième partie de cette thèse est consacrée à l'étude des variations latérales des déformations récentes dans le prisme Himalayen. Les études paléosismologiques menées depuis 15 ans en Himalaya ont permis d'étendre le catalogue des événements majeurs sur le dernier millénaire. En intégrant à la fois les contraintes disponibles sur ces séismes majeurs et le chargement séculaire, nous étudions les variations spatiales et temporelles de l'état de contrainte du chevauchemlent Himalayen principal le long de l'arc. Nos calculs montrent que même dans la cas où un méga-séisme se serait produit en ~1100 en Himalaya, le chargement séculaire a quasiment compensé la chute de contrainte associée. Les modèles montrent aussi que le séisme du Shillong ne joue pas un rôle majeur sur les contraintes accumulées sur le chevauchement Himalayen depuis 1897 au niveau du Bhoutan. Enfin, nous présentons les premières contraintes sur la tectonique active au Bhoutan. L'étude morphotectonique réalisée au centre du Sud-Bhoutan a permis de montrer que la majeure partie de la déformation Holocène est accommodée au niveau du front Himalayen, comme au Népal. Nous montrons aussi des évidences de ruptures de surface avec des décalages verticaux de plusieurs mètres associés à deux séismes majeurs sur le dernier millénaire. Ces évidences remettent en cause l'interprétation d'un méga-séisme en ~1100.

This thesis aims at deciphering the respective roles of autogenic processes and allogenic forcing in the post-orogenic evolution of a coupled mountain catchment/foreland megafan system, with a focus on the Lannemezan megafan and its mountainous catchment in the northern Pyrenees (France). AFT data are consistent with previously published thermochronological data, showing (i) the main exhumation phase of the Axial Zone (AZ; ~50-30 Ma) with lateral variations in the exhumation rates, (ii) a later unroofing of the North-Pyrenean Zone (NPZ; since ~20 Ma), and (iii) a late (~20-15 Ma) exhumation phase only recorded in the western part of the range (west of the Neste River). The megafan was thus built at a late stage of the orogeny and its source catchment lays on the boundary between two tectonic domains. Petrographic data shows that the source area encompasses the NPZ and AZ but modern river sands yield more granitic material than the megafan sediments. This discrepancy suggest that (i) the foreland deposits are locally sourced, but this is not consistent with the definition of a megafan, (ii) the crystalline massifs may have been insufficiently unroofed, but this is inconsistent with the presence of granitic material in Eocene deposits of the south foreland, (iii) the main drainage divide migrated to the south around that period, but this seems inconsistent with the exhumation history of the range as inferred from thermochronological data. Comparison of my estimation of the sedimentary budget and AFT-derived exhumation rates suggests that the current watershed of the Neste is more than sufficient to provide the material for the building of the Lannemezan megafan. This implies that (i) the catchment boundaries have not necessarily changed since the Miocene and (ii) a significant portion of the sediments bypassed the megafan. I dated the abandonment of Quaternary surfaces using cosmogenic nuclides (10Be, 26Al). The abandonment of the megafan (~300 ka) and of the alluvial terraces (~100ka, ~20ka and ~1ka) suggests that incision episodes are triggered by cold-to-warm climatic shifts. The terrace slopes increase with time: I propose that the region was subjected to northward tilting by flexural isostatic rebound as a response to erosional unloading of the range during the Quaternary. From the morphometric analysis (χ-proxy, steepness index, concavity, long river profiles) I show (i) drainage reorganization near the megafan apex and (ii) no influence on the drainage network from active tectonic structures or base-level variations. The abandonment and of the Lannemezan megafan is thus linked to the capture of the Neste by the Garonne and its episodic incision is linked to Quaternary climatic cycles. Numerical modelling adequately reproduce the building of a large megafan in 15 Ma by progradation of the deposits carried by a central river that becomes dominant in the sediment routing. I suggest that, in the northern Pyrenean foreland, the Atlantic Ocean to the west and large, efficient rivers create the conditions needed for a river to be singled out and build a megafan. Negligible subsidence at during the Miocene may have encouraged progradation. I suggest that, in the northern Pyrenees, the main driving force for long-term entrenchment of the megafan is linked to autogenic processes but these are modulated by short-term climatic changes (and possibly isostatic movements). Numerical modelling also evidenced the strong coupling between the mountainous catchment and the basin and the characteristic response times (on the order of 10 ka). / Le but de cette thèse est de comprendre ce qui contrôle l’évolution géomorphologique post-orogénique d’un ensemble composé d’un bassin versant montagneux et de son piedmont à partir de l’exemple du méga-cône alluvial de Lannemezan et son bassin versant au flanc nord des Pyrénées. Les données de thermochronologie in-situ et détritiques confirment les phases d’exhumations : (1) une phase d’exhumation majeure (~50-30 Ma), enregistrée dans toute la ZA bien qu’avec des vitesses différentes à l’est (Ariège) qu’à l’ouest, (2) une excavation plus tardive de la ZNP (à partir de ~20 Ma), et (3) une phase tardive d’exhumation (vers 15-20 Ma) enregistrée seulement à l’ouest. Le méga-cône se forme donc à une période post-orogénique où la tectonique cesse de jouer un rôle dominant, mais le bassin versant d’alimentation du cône se situe à une position charnière entre deux zones qui s’exhume différemment. Les données de pétrographie des sédiments indiquent une zone source qui inclut la ZNP et la ZA. La proportion de matériel cristallin dans les sédiments Oligo-Miocène est plus faible que celle des sédiments des rivières actuellesce qui pourrait suggérer que (1) les massifs cristallins n’étaient pas suffisamment exhumés mais cela semble incohérent avec les études précédentes, (2) que la ligne de partage des eaux a migré vers le sud, mais ce n’est pas cohérent avec l’histoire d’exhumation de la chaîne. J’ai comparé le volume érodé dans le bassin pour produire le méga-cône et le flux sédimentaire estimé à partir des taux d’exhumation : il semble que le bassin versant actuel de la Neste est plus que suffisant pour produire les sédiments du cône, ce qui suggère que (1) la taille du bassin versant n’a pas nécessairement changé et (2) une partie des sédiments a été directement exportée vers le bassin marin. Les résultats des datations par nucléides cosmogéniques (10Be et 26Al) suggèrent que le cône est abandonné vers ~300 ka et les terrasses alluviales qui bordent les rivières du piémont (~100ka, ~20ka et ~1ka) pourraient suggérer un contrôle climatique sur l’incision (cycles glaciaires/interglaciaires). Pour expliquer la pente des terrasses je propose l’hypothèse d’un basculement de la région en lien avec un rebond isostatique au Quaternaire. L’analyse morphométrique du réseau de drainage montre (i) une réorganisation par capture dans la partie sud du cône et (ii) qu’il n’y a pas de structure tectonique active ou de variation du niveau marin. J’ai utilisé un modèle numérique pour comprendre les conditions de formation et d’incision d’un méga-cône et comparé les résultats avec les observations faites sur le cône de Lannemezan et son bassin versant. Les résultats suggèrent que le temps de construction d’un tel cône est de l’ordre de 15 Ma, par progradation des dépôts d’une rivière centrale qui domine les apports du piémont. La comparaison avec le cas de Lannemezan suggère que (1) l’Océan Atlantique à l’ouest et les grandes rivières, au transport efficace créent un environnement propice à la formation du méga-cône; (2) la subsidence négligeable du piémont pourrait avoir favorisé la construction du méga-cône; (3) les processus autogéniques sont suffisants pour produire au long terme une incision permanente du cône par sa rivière d’alimentation. Les forçages externes (tectoniques, climat) ne font que moduler l’incision autogénique. La réponse du paysage à un forçage climatique cyclique montre un temps caractéristique de l’ordre de 10 ka. Finalement, la modélisation a permis de mettre en évidence les interactions entre le piémont et la chaîne.

Les matériaux granulaires présentent un large spectre de propriétés mécaniques. Développer des modèles constitutifs permettant d'intégrer ces caractéristiques dans le cadre de simulations à l'échelle de l'ouvrage demeure un réel challenge scientifique. A cet égard, les approches multi-échelles constituent aujourd'hui une voie très prometteuse. Elles permettent de faire émerger des propriétés macroscopiques à partir de modèles micromécaniques calibrés à l'échelle microscopique.Parmi les modèles multi-échelles, le modèle H marque une avancée majeure pour la prise en compte des effets de la microstructure dans le comportement des matériaux granulaires. La structure du matériau granulaire est décrite par une distribution d'hexagones orientés dans l'espace. A partir d'opérations d'homogénéisation, les contraintes et les déformations incrémentales sont reliées à l'échelle de la distribution, donnant lieu à un modèle de comportement qui a la capacité à reproduire propriétés mécaniques essentielles des matériaux granulaires.Nous étudions dans un premier temps les propriétés mécaniques de l'assemblage hexagonal de grains, élément de base du modèle H, afin d'identifier les conditions menant à sa déstabilisation. Nous réalisons dans un second temps une étude de sensibilité du modèle constitutif vis-à-vis des paramètres micro-mécaniques et microstructurels. Enfin, nous démontrons les capacités opérationnelles du modèle à partir d'essais triaxiaux non drainés réalisés sur un sable lâche liquéfiable.Dans un second temps, le modèle H est implémenté en tant que loi constitutive dans un code de calcul aux différences finies. Des simulations d'essais biaxiaux non homogènes sont conduites afin d'explorer les capacités du modèle à reproduire les différents modes de rupture observés en laboratoire. L'utilisation du modèle H pour modéliser des essais biaxiaux drainés et non drainés met clairement en évidence l'influence de la microstructure sur la réponse mécanique des matériaux granulaires. Enfin, le modèle H est utilisé dans le cadre d'une simulation hydro-mécanique couplée à l'échelle de l'ouvrage pour modéliser le chargement d'une fondation superficielle et la rupture d'une digue soumise à une crue. / Granular materials generally exhibit a broad spectrum of mechanical properties. Developing constitutive models to integrate these properties in the context of simulations at the structure scale remains a real scientific challenge. In this respect, multi-scale approaches offer very promising solutions as they allow the emergence of macroscopic properties from micromechanical models calibrated on a microscopic scale.Among the multiscale models, the H model marks a major step forward in taking into account the effects of the microstructure in the behavior of granular materials. The structure of the granular material is described by an assembly of hexagons, oriented in space. From homogenization operations, stresses and incremental strains are related to the scale of the assembly, giving rise to a constitutive model that has the ability to reproduce the essential mechanical properties of granular materials.We first study the mechanical properties of the hexagonal grain assembly in order to identify the conditions leading to the triggering of its instability. We then carry out a study of the sensitivity of the constitutive model with respect to micro-mechanical and microstructural parameters. Finally, we demonstrate the operational capacities of the model from triaxial undrained tests carried out on a liquefiable loose sand.In a second step, the H model is implemented as a constitutive law in a finite difference code. Simulations of non-homogeneous biaxial tests are carried out in order to explore the model's capacities to reproduce the different failure modes observed in the laboratory. The use of the H model to model drained and undrained biaxial tests highlights the influence of the microstructure on the mechanical response of granular materials. Finally, model H is used in numerical simulations at the structure scale to model the loading of a shallow foundation and the failure of a levee subjected to a flooding event.

Dans les régions de montagne, la prévision des évènements gravitaires reste un défi pour la gestion des risques. Des méthodes de calcul telles que la méthode des éléments discrets (DEM), où les particules interagissent les unes avec les autres pour restituer un comportement global d’une masse granulaire, ont été utilisées pour aborder ce type de problématique. L’application de la DEM reste normalement limitée aux évènements de petits volumes impliquant un nombre de blocs plutôt faible, puisque les temps de calcul peuvent devenir rapidement prohibitifs avec l’augmentation du nombre de particules. Les méthodes de calcul continues sont donc une alternative intéressante car elles permettent de réduire les temps de calcul. Elles nécessitent cependant la définition d’une loi de comportement macroscopique capable de représenter correctement les principaux traits de comportement mécanique du matériau au sein de la masse. L'objectif principal du travail de thèse réside dans le développement d’un outil numérique permettant de modéliser certains aléas gravitaires tels que les écoulements en masse. Notre choix s’est porté sur une méthode Lagrangienne-Eulérienne (méthode des points matériels – MPM) capable de gérer de grandes déformations tout en bénéficiant des principaux avantages de la méthode des éléments finis (FEM). La méthode utilise une grille Eulérienne fixe sur laquelle se déplacent des points matériels pendant les simulations. Un outil numérique, nommé MPMbox (2D et 3D), a été développé entièrement durant la thèse en C++. Le code a été validé à l'aide d'une série de solutions analytiques en quasi-statique (tests géotechniques standards) ainsi que par des applications de la littérature incluant des déformations importantes et rapides (tests d'affaissement). Après validation, le code a été confronté aux prédictions d’un outil de calcul DEM (DEMbox) dans le cadre de simulations numériques impliquant l'écoulement (initiation, régimes transitoires, propagation et arrêt) d'un matériau granulaire (particules sphero-polyhédriques) sur un plan incliné. Les résultats ont été comparés en termes de distance de propagation, de forme du dépôt et d'énergies dissipées à l'interface et dans la masse pendant l'écoulement. Pour les applications qui ont suivies, des éléments discrets ont été couplés à la MPM afin qu'un bloc rigide (DEM) puisse interagir avec un sol déformable (MPM). Cette application a consisté en l'analyse (2D) de la collision entre un bloc rocheux rigide (rond ou carré) et un sol bicouche élastoplastique. Les investigations ont été largement basées sur la mesure de coefficients de restitution (rapport des énergies cinétiques avant et après impact) qui reste difficile à déterminer expérimentalement. / In mountainous regions, the prediction of gravitational phenomena remains a challenge for the management of risk. Computational methods such as the Discrete Element Method (DEM) have been used for the modeling of these types of phenomena, where particles interact with each other to give an overall behavior of the mass. Its application can be somewhat restricted to small and medium number of blocks, since the computational time can easily become too large. Continuum analyses are therefore an attractive approach, which can reduce the computational times, but that rely on a constitutive law to represent the behavior within the mass. The main objective of this PhD was to develop a numerical tool that allowed the modeling of some specific gravitational hazards, such as the flowing of mass. A Lagrangian-Eulerian method such as the Material Point Method (MPM) is able to handle large deformations, while preserving most of the capabilities of the Finite Element Method (FEM). The method uses an Eulerian grid which is only used as a numerical scratch-pad, and remains fixed during simulations. A numerical tool named MPMbox (2D and 3D) was then developed from the ground up using C++. The code was validated using a series of analytical solutions for quasi-static analysis (some standard geotechnical tests), as well as simulations including large and rather rapid deformations (slump tests). After validation, the code was first used to make a numerical comparison with the DEM. In the comparison, a parametric survey was carried out during which the flow of a granular material on a sloped surface was simulated. Results were compared in terms of run-out distance, spread of the deposit and energy dissipated at the interface and within the mass during the flow. For a second study, discrete elements were coupled with MPM so that a rigid block could interact with a deformable soil. This application consisted in the (2D) analysis of the collision between a discrete block (round and squared) and a bounded elasto-plastic double-layered soil (soft over hard layers). The investigations were largely based on the measurement of the restitution coefficient (ratio of kinetic energies before to after the impact), which cannot be easily determined experimentally.

Les éboulements rocheux sont difficilement prévisibles en raison de leur soudaineté et du peu d'informations disponibles sur la structure interne du massif rocheux. L'analyse en retour d'éboulements dans des falaises calcaires a montré qu'ils étaient initiés par la fracturation de ponts rocheux et que leur persistance était le principal paramètre à prendre en compte pour des calculs de stabilité. En 1er lieu, nous avons utilisé l'information contenue dans le bruit sismique afin d'évaluer le degré d'instabilité d'une écaille rocheuse. Les données viennent de la falaise calcaire de Chamousset (Vercors) où un réseau sismique a fonctionné 4 mois jusqu'à 2 semaines avant la chute d'une écaille de 21000 m3. Le calcul des spectres de bruit sismique a permis d'identifier la 1ère fréquence propre de l'écaille qui diminue significativement avant l'éboulement. Ce même réseau a enregistré des milliers d'événements dont les plus énergétiques ont été classés. Certains de ces signaux ont des arrivées d'ondes P et S distinctes, permettant d'analyser leur mécanisme de rupture et de calculer leur magnitude et distance hypocentrale. Des simulations favorisent l'hypothèse de ruptures en traction pour expliquer les polarisations. D'autres simulations confirment que les événements situés à l'interface écaille-massif entraînent la résonance de l'écaille, expliquant l'observation de codas mono-fréquentielle avec une amplitude sismique plus forte sur l'écaille. Enfin, l'étude de la réponse thermomécanique d'une 2nde écaille sur un an montre l'influence des paramètres météorologiques sur ses fréquences propres. Les effets réversibles et irréversibles sur ces fréquences ont été identifiés. / Rockfalls are difficult to predict events, owing to the suddenness of the phenomenon and to the little information available about the mass internal structure. Back analysis of collapses in steep limestone showed that rockfalls were initiated by intact rock bridge failures and that the persistence of joints was the main factor controlling the failure probability. First, we used the seismic noise recorded by autonomous monitoring systems to evaluate the degree of coupling of a prone-to-fall column to the rock massif. Seismic data were acquired on top of a limestone cliff in Chamousset (Vercors) during 4 months until 2 week prior the collapse of a 21,000 m3 column. The column 1st natural frequency was easily identified from ambient noise processing and exhibited a significant drop prior the column failure. Thanks to the monitoring system, thousands of seismic events were detected, out of which, the most energetic ones were classified. Some of these events exhibit clear P and S waves, making possible an analysis of the rupture mechanism as well as the calculation of their magnitude and hypocentral distance. Numerical simulations of rock bridge rupture allowed favouring traction rupture to explain the observed polarizations. Other numerical simulations showed that rupture events can trigger the column resonance, explaining the recurrent presence of specific frequencies in the late coda as well as the higher seismic motions on the column. Finally, we studied the thermomechanical response of another column. The influence of meteorological parameters on the column natural frequencies was shown. Both reversible and irreversible effects on the 1st natural frequency were identified.

Ce travail de thèse est consacré au développement d'un capteur acoustique miniature à transduction capacitive destiné à être intégré dans un système RFID afin de dépasser certaines des limitations actuelles de ce dernier. La configuration originale du capteur acoustique étudié lui offre les avantages d'avoir une bonne performance tout en gardant une forme simple qui peut être aisément réalisée avec les techniques MEMS pour une production à grande échelle. Ce transducteur est constitué d'une membrane circulaire ou carrée et une électrode arrière centrée, de même forme mais de dimensions plus petites, séparées par une très fine couche de fluide, ainsi qu'une petite cavité située à la périphérie de l'électrode, de dimensions extérieures très proches de celles de la membrane. Le comportement de ce capteur est analysé en détail. Cette étude se base sur deux approches mathématiques originales (analytique et numérique), dont les résultats convergent malgré un niveau différent des hypothèses-simplificatrices sur lesquelles reposent ces deux modèles.Finalement, une méthode de réalisation du transducteur en technologie hybride, qui associe le procédé MEMS avec les techniques classiques des circuits imprimé, est présentée. Le prototype développé est aussi caractérisé expérimentalement et les résultats obtenus correspondent bien aux caractéristiques fournies par les modèles théoriques. / This work is devoted to the development of capacitive miniature acoustic sensor to be integrated in an RFID system in order to exceed some of its current limitations.The unusual configuration of the studied sensor offers the benefits of having a significant performance while keeping a simple form that can be easily achieved with standard MEMS techniques and thus successively marketed in perspective. This transducer consists of a circular or square membrane, a centered backing electrode with the same shape, but with smaller dimensions which define a thin fluid layer between them, and a small cavity at the periphery of the backing electrode whose external dimensions are very close to the ones of the membrane. The architecture of the developed sensor is the subject of a deep theoretical study of its behavior. This investigation is based upon two original mathematical approaches (analytical and numerical ones) whose results show a high convergence despite the different-levels of simplifying assumptions on which these models are based.Finally, the transducer is made with the use of a hybrid technology, that includes the MEMS process associated with conventional printed circuit manufacturing techniques, which is presented step by step. The developed prototype is characterized experimentally and the measured results correspond to the ones predicted by theoretical models.

La chaîne Himalayenne est l'un des exemples les plus spectaculaires de déformation active à la surface de notre planète. Au cours des derniers siècles, de nombreux séismes majeurs (Mw > 7.5) ont affecté cette zone et le lourd bilan humain du séisme de Gorkha en 2015 a une nouvelle fois rappelé l’importance de parvenir à mieux estimer l’aléa sismique de cette région. Bien que les structures lithologiques et tectoniques semblent montrer au premier ordre une cylindricité le long des 2500 km de l'arc Himalayen, de nombreuses études, menées principalement au Népal, ont permis de mettre en évidence des variations latérales structurales, thermochronologiques, morphologiques, gravimétriques, sismologiques ou géodésiques. Le rôle de ces variations latérales sur la segmentation sismique reste cependant mal contraint. La taille maximale et la probabilité d’occurrence de ces méga-séismes sont donc toujours matière à débat.Ainsi, le comportement sismique de l’Himalaya du Bhoutan reste énigmatique. Pour certains les faibles taux de sismicité observés actuellement font de ce royaume une zone asismique de 350 km de long. Pour d’autres, à l’instar de l’Himalaya du Népal, la faible sismicité observée est associé à une forte accumulation de contrainte susceptible de générer des séismes majeurs. L'objectif de ce travail de thèse est d’améliorer notre connaissance de cette région en quantifiant la déformation à différentes échelles spatiales et temporelles via des études morphotectoniques et paléosismologiques.La première partie de cette thèse vise à quantifier les mouvements verticaux à l’ouest Bhutan et le long du Main Frontal Thrust, structure la plus frontale située au sud Bhoutan. Trois campagnes de terrain ont été réalisées permettant l'échantillonnage (1) de terrasses alluviales le long du front afin de quantifier et d'étudier les variations du soulèvement Holocène, (2) de bassins versants pour l'étude de la dénudation court-terme (< 20 ka) dérivée des cosmonucléides ($^{10}$Be) et (3) de terrasses alluviales dans le Moyen-Pays pour quantifier l’incision Holocène. Les vitesses verticales obtenues au front sont comparables à celle proposées le long du reste de l’arc himalayen, suggérant une cinématique relativement simple. Par contre, nos résultats indiquent une variation de la géométrie du chevauchement himalayen (Main Himalayan Thrust) entre l’est Népal et l’ouest Bhoutan.La seconde partie porte sur plusieurs études paléosismologiques le long du front ouest et centre Bhoutanais. Six sites différents ont été étudiés au cours de trois campagnes de terrain. La datation et la modélisation de charbons détritiques a permis de mettre en évidence l'occurrence d’au moins cinq séismes majeurs durant les derniers 2700 ans, faisant du Bhoutan une zone sismiquement aussi active que le Népal. A une échelle régionale, cette étude apporte donc de nouvelles contraintes et contribue au débat sur la possibilité d’occurrence d'un séisme de magnitude 9 le long de l'arc Himalayen. / The Himalayan arc is one of the most active intra-continental mountain belts in the world. Over the last centuries, several major earthquakes (Mw > 7.5) have struck this arc. The dramatic effects of the Gorkha earthquake sequence in 2015 pointed once again the crucial need to improve seismic hazard assessment of this area.Geological explorations of the Himalayas since the late 19th century have emphasized a 2500-km-long roughly cylindrical structure, with striking continuity of main units and thrust faults. However recent geophysical and geological investigations have revealed lateral variations. The relationship between these variations and earthquakes segmentation along the arc remains poorly constraint. The maximum size and the occurrence probability of such earthquakes are still a matter of debate.For instance, the seismic behavior of Bhutan remains enigmatic. The present-day low seismicity rate observed in this area can reflect two opposite fault behaviors: an aseismic creeping zone or a zone of stress accumulation for future great earthquakes as the others parts of Himalayas. The main objective of this thesis is to bring new constraints on the deformation of the Bhutan Himalayas, at different space- and time-scales, through morphotectonic and paleoseismological approaches.The first part of this thesis focuses on the vertical deformation assessment along a N-S transect in western Bhutan and along the Main Frontal Thrust, which is the southern-most thrust in southern Bhutan. In the past three years, we have carried out three fieldwork campaigns to sample (1) frontal terraces to assess Holocene uplift rates, (2) watershed basins to quantify short term denudation rates (< 20 ky) derived from cosmonuclides $^{10}$Be and (3) hinterland alluvial terraces to quantify the Holocene incision rate. Frontal Holocene uplift rates obtained in Bhutan are consistent with those obtained in the others parts of Himalayas. Furthermore, our results reveal a variation in the geometry of the Main Himalayan Thrust between eastern Nepal and western Bhutan.The second part focuses on several paleoseismic studies along the west and central bhutanese Himalayan front. Different sites were investigated during three fieldwork campaigns. Detritic charcoals sampling and modeling suggest the occurrence of at least five surface-rupturing earthquakes during the last ~2700 years. These results demonstrate that the present-day low seismicity rate observed in Bhutan is not representative of the seismic activity at longer time scale. At regional scale, they also take part of a broader discussion on the probability of occurrence of a magnitude 9 earthquake along the Himalayan arc.

La réfrigération magnétique est une alternative pertinente dans un contexte où les gaz réfrigérants sont soumis à des restrictions environnementales. Ces restrictions nécessitent l'évolution de la technologie actuelle ou bien l'émergence d'une nouvelle, d'où l'opportunité pour la réfrigération magnétique de prouver son potentiel. En effet, elle pourrait s'avérer énergiquement plus efficace et avec des densités de puissance supérieure. Ces travaux de thèse apportent des réponses sur le potentiel de la réfrigération magnétique. Dans cette logique, la thermodynamique et le magnétisme, outils indispensables à notre étude, sont développés dans le cas des matériaux à effet magnétocalorique. Puis, nous verrons que les caractérisations de ces derniers sont en mesure de fournir des modèles matériaux cohérents et réalistes, si des précautions sont prises. L'effet magnétocalorique étant limité en termes de variation de température, nous allons étudier différentes structures de réfrigération. Enfin, des modèles numériques sont développés pour permettre d'optimiser les structures à régénérations actives, qui sont les plus utilisées. Ces modèles doivent permettre de dimensionner des systèmes proches de leurs optimums. / Magnetic refrigeration is a relevant alternative in consideration of environmental restrictions of refrigerants gases. These restrictions require to improve the current technology or to pave the way for a new one, hence the opportunity for magnetic refrigeration to demonstrate its potential. Indeed, it could be energetically efficient and with higher power densities. This work aims to estimate the potential of magnetic refrigeration. Magnetism and thermodynamic, essential tools for our study, are developed in a case of magnetocaloric effect. With some care, we show that material characterizations are able to give consistence and relevant model. Magnetocaloric effect suffers of small temperature variations; therefore structures that increase the temperature span and give competitive system are studied. Finally numerical models are developed to optimize active magnetic regenerators, which are currently the most used. These models are used to calculate and design systems close to their optimum.

L'objectif de cette thèse est de bâtir un modèle prédictif d'une unité industrielle continue de dissolution de minerai de nickel/ et précipitation de fer simultanées. La méthodologie adoptée consiste dans un premier temps à réaliser des expériences en réacteur fermé ou semi-fermé pour identifier et modéliser séparément les différents phénomènes en jeu, puis dans un second temps à construire un modèle d'unité continu les rassemblant tous. Ce travail a permis certaines avancées tant sur le système retenu que sur les méthodes et modèles adoptés. Tout d'abord la thermodynamique des solutions Na-Fe(III)-H2SO4 aux alentours de 100°C a été clarifiée et la constante de solubilité de la natrojarosite Na0,84H0,16Fe2,90(SO4)2(OH)5,7 a été déduite. Ensuite les paramètres influant sur les précipitations des jarosite de sodium et potassium ont été mis en évidence et leur cinétiques de croissance ont pu être déterminées grâce à l'utilisation originale de la méthode des caractéristiques. Puis les cinétiques de dissolution des minerais ont été obtenues, en prenant en compte différentes phases du minerai et leurs granulométries. Après avoir déduit tous les paramètres nécessaires, un modèle original, permettant de simuler une cascade de réacteur de dissolution/précipitation avec recyclage a été construit puis validé. Cet outil de conception, couplé à une étude technico-économique peut permettre d'optimiser le procédé. / The goal of this thesis is to build a predictive model for a continuous industrial unit combining simultaneously nickel ore dissolution and iron precipitation. A two steps method was used. First experiments in batch or semi-batch reactors were done in order to understand and model separately the different phenomenon that take place. Then all the equations and associated constants were used to build a model. This work allowed some advances both on the system studied and on the methods used. First thermodynamics of Na-Fe(III)-H2SO4 solutions in the 70-100°C temperature range was clarified and natrojarosite solubility constant Na0,84H0,16Fe2,90(SO4)2(OH)5,7 was deduced. Secondly parameters acting on sodium and potassium jarosite precipitation were highlighted and their growth kinetics were deduced through an original use of caracterisctics method. Then ore dissolution kinetics were found, taking into account different ores phases and their granulometry. After deducing all the necessary parameters, an original model allowing to simulate a cascade of dissolution/precipitation reactors with recyling loop was build and validated. This conception tool, coupled with a technico-economic study allows the optimisation of such a process.