Dans la région de Dimbokro-Bongouanou (Centre Est de la Côte d'Ivoire), l'alimentation en eau des populations se fait essentiellement au moyen de puits et forages. Ceux-ci ont été implantés dans les aquifères du socle au cours des différentes campagnes d'hydraulique villageoise mises en place dans ce secteur. L'analyse du taux de réussite et de la productivité de ces forages montre que leur implantation n'a pas toujours été optimale. Afin d'optimiser l'implantation de nouveaux forages, il semble donc nécessaire et primordial de caractériser l'hydrodynamique du socle fracturé. A cette fin, la présente étude procède d'abord à une caractérisation géologique puis hydrogéologique de la région de Dimbokro-Bongouanou. Ensuite, sont élaborés des modèles conceptuels géologiques et hydrogéologiques, qui décrivent la structure de l'aquifère et précisent les propriétés hydrodynamiques de ses différentes couches composites. Ces modèles géologiques et hydrogéologiques permettent de simuler numériquement l'écoulement régional. La caractérisation géologique de la région de Dimbokro-Bongouanou est basée sur une description des différentes formations géologiques de la région ; des hypothèses sur leur mise en place et les cycles d'érosion et d'altération qu'elles ont subis ; et une cartographie des discontinuités régionales par photo-interprétation d'image Radar rehaussée par la technique de filtrages adaptatifs. Les discontinuités identifiées sur l'image rehaussée sont validées à partir d'informations géologiques (données de fracturation mesurées à l'affleurement et données de fracturation relevées sur carte géologique).La caractérisation hydrogéologique de la région de Dimbokro-Bongouanou est basée sur la description des aquifères de la région et l'identification du rôle des discontinuités régionales dans les écoulements souterrains. La technique proposée à cet effet consiste à analyser les corrélations entre le débit des forages et leur distance au linéament le plus proche.Il est ainsi constaté que le débit des forages est inversement proportionnel à la distance à la discontinuité la plus proche. En d'autres termes, plus le forage est proche d'une discontinuité régionale, plus son débit est important. L'interprétation des lithologs, des données géologiques à l'affleurement et des données de résistivités électriques permet de proposer un modèle conceptuel géologique représentatif du profil d'altération type pour chaque domaine géologique de la région. En particulier, le modèle proposé décrit la structuration verticale du profil d'altération pour les trois principales formations géologiques de la région: les granites, les métasédiments et les métavolcanites. La caractérisation de la structure de l'aquifère permet ensuite de proposer un modèle conceptuel hydrogéologique pour chacun des profils d'altération basé sur l'analyse statistique des arrivées d'eau et de l'interprétation des données d'essais par pompage de différents forages. Cela permet d'identifier pour chaque formation les profondeurs de fréquence, de densité et de débit des arrivées d'eau ainsi que les couches du profil associées aux zones de conductivité et de productivité (débit spécifique) élevées.L'ensemble de ces informations géologiques et hydrogéologiques permet de proposer une modélisation conceptuelle de l'hydrosystème étudié (une partie du Bassin du NZI,environ 6300km2) ainsi qu'une modélisation numériques de l'écoulement, celle-ci servant de base à une réflexion sur la pertinence de l'information disponible et d'une modélisation discrète dans ce type d'environnement au vu des objectifs affichés (modèle de gestion et d'aide à la décision pour l'implantation de nouveaux forages). / The hydraulic characterization of fractured basement in Dimbokro-Bongouanou region, Côte d'Ivoire, is an essential step to achieve an optimization of new well locations, in terms of success rates and productivity.To this end, this study first proposes conceptual models of weathering profile for the three main geological formations of the region: granites, metasediments and metavolcanic rocks. Based on in-situ observations, interpretations of lithology and cycles of erosion and alteration, the typical weathering profiles are proposed for each type geological. These are formed from top to bottom of layers alloterite, isalterite, fissured and fractured fresh basement. The fissured layer shows the maximum density of fractures in the metavolcanic rocks and métasédiments, while the fractured fresh basement concentrates the highest density of fractures in granites. Moreover, interpretation of data from pumping tests conducted in various hydraulic programs indicate that the productivity of wells is not proportional to the depth due to variation, at region scale of the thickness of the fissured layer and fractured fresh basement which totalize the maximum of water arrivals. The spatial variation of hydraulic conductivity and productivity (specific rate) is influenced by the thickness of the fissured layer, alloterite layer of the weathering profile and topography. The regional fractures of the fractured fresh basement are also areas of high productivity, as evidenced by the negative correlation between rate of productive wells and distance to these fractures. A methodology of regional fractures mapping from the use of a Radar image and hydrogeological data is suggested.All of these geological and hydrogeological informations are finally integrated into a flow model and the simulated groundwater levels are compared with measured piezometric. It thus appears that the values of hydraulic conductivities deduced from pumping tests are 103 times lower than those that reproduce the regional flow. In addition, a dual approach (equivalent continuum for alterite layer and fissured layer vs. discrete fracture of the base) confirms the role of accidents in the regional flows.

L'insuffisance cardiaque est une maladie à grande incidence dont le traitement définitif est difficile. Les pompes d'assistance ventriculaire ont été proposées comme thérapie alternative à long terme, mais la technologie est relativement jeune et selon son design, axial ou centrifuge, le dispositif favorise soit l'hémolyse, soit la stagnation de l'écoulement sanguin. Les pompes à écoulement mixte, combinant certaines propriétés des deux types, ont été proposées comme solution intermédiaire. Pour évaluer leurs performances, nous avons effectué des comparaisons numériques entre huit pompes, deux axiales, deux centrifuges, et quatre mixtes, en employant un modèle Windkessel du système cardiovasculaire avec paramètres optimisés pour l'insuffisance cardiaque résolu avec une méthode Radau IIA3, une méthode de résolution de système d'équations différentielles ordinaires L-stable appartenant à la famille des méthodes Runge-Kutta implicites. Nos résultats semblent suggérer que les pompes d'assistance mixtes ne démontrent qu'un léger avantage comparativement aux autres types en terme de performance optimale dans le cas de l'insuffisance cardiaque, mais il faudrait effectuer plus d'essais numériques avec un modèle plus complet, entre autres avec contrôles nerveux implémentés. / Heart failure is a disease with a high incidence rate that is difficult to treat definitively. Ventricular assist devices have been proposed as alternative long-term therapeutic options, but the technology is fairly recent and depending on the design, axial or centrifugal, the device will favor hemolysis or blood flow stagnation. Mixed flow devices, which combine properties of the two previous types, have therefore been proposed as intermediate solutions. To evaluate their performance, we did a numerical comparison between eight pumps (two axial, two centrifugal, and four mixed) using a Windkessel model of the cardiovascular system with parameters optimised for heart failure and solved using a Radau IIA3 method, an L-stable method used to numerically solve systems of ordinary differential equations that belongs to the family of implicit Runge-Kutta methods. Our results suggest that in terms of optimal performance, mixed flow pumps are modest improvements over the other two types, and more tests would have to be performed using a more complete model with, among other modifications, neural control systems implemented.

Nous avons étudié le transfert de chaleur et de matière au cours de l'histoire de la Terre primitive à de multiples échelles en utilisant des modèles numériques. Deux systèmes différents sont abordés. Tout d’abord, nous nous concentrons sur les premiers stades de la formation du noyau terrestre lorsque le fer se sépare des silicates et descend vers l'intérieur de la planète. Au cours de la différenciation, des interactions chimiques et thermiques se produisent entre les gouttes de fer dispersées dans des silicates fondus formant un océan de magma. Nous étudions le transport chimique des éléments traces à l'intérieur et autour des gouttes. Nous tirons quelques relations fonctionnelles dépendante du régime dynamique d'écoulement et montrons que le système tend à être en équilibre chimique extrêmement rapidement par rapport à l'échelle de temps de la descente de la goutte de fer. Peu de temps après la fin de l'accrétion de la Terre, la fusion extensive de son intérieur profond ainsi que la formation d'un océan de magma en surface a lieu. Comme le rayonnement de la chaleur dans l'espace est très efficace, les silicates fondus superficiels cristallisent très rapidement. L'histoire thermique de la couche liquide enterrée, appelée océan de magma basal (OMB), se déroule sur une longue période de temps et il est proposé que ses restes soient aujourd'hui observables sous forme de poches partiellement fondues au dessus de frontière noyau-manteau. Nous déterminons les paramètres régissant un système convectif dans lequel se produit une transition solide/liquide. Les lois d'échelle ainsi obtenues ont été appliquées à l'OMB et indiquent que la différence de température qui peut être maintenue dans les couches limites supérieure et inférieure de l'OMB est infime. Par conséquent, la température du noyau suit la température de liquidus à la base du manteau et ainsi la vitesse de refroidissement de l'OMB doit être la même que celle du noyau de la Terre. / We have studied the heat and mass transfer during the early Earth history at multiple scales and for multiple systems by means of numerical computing. Two different systems are approached. Firstly, we focus on the early stages of the Earth core formation when iron segregates from silicates and descends toward the interior of the planet. During the differentiation there are chemical and thermal interactions between dispersed iron blobs and surrounding molten silicates. We study the chemical transport of trace elements within and around the drops. We derive functional relations between critical parameters and show that the system tends to be in chemical equilibrium.During the accretion process of the Earth, extensive melting of its deep interior as well as formation of shallow magma oceans occurred.As heat radiation into space happens with high efficiency, surface molten silicates crystallize very rapidly, in about 10 My. The thermal history of the buried liquid layer, called the basal magma ocean (BMO), proceeds over a long time and it is proposed that its remnants are nowadays observable as partial melts in the core-mantle boundary region.We develop numerical models of the thermal history of the crystallizing basal magma ocean that enable to study coupling between the mantle and the core in the presence of the BMO. We derive parametrized relations for this convective system that undergoes solidification/melting. Obtained scaling equations applied to the BMO indicate that the temperature difference that can be maintained across the top and bottom boundaries of the BMO is minute. Hence, the temperature of the core follows the temperature of liquidus at the bottom of the mantle and thus the rate of the BMO cooling must be the same as that of the Earth's core.

Le concentrateur à réflecteur linéaire de Fresnel (LFR) est une technologie solaire thermodynamique en plein essor : petites applications industrielles (chaleur, froid, électricité) à centrales électriques (10-100 MWel). Ce travail de thèse établit un modèle global du procédé solaire, en régime permanent, pour un prédimensionnement du système. Le modèle comprend trois parties chaînées : (i) les transferts radiatifs dans le concentrateur optique, modélisés précisément par une méthode de Monte Carlo (environnement EDStar) ; (ii) les transferts thermiques dans le récepteur, évalués analytiquement (puissances, températures) ; (iii) le cycle thermodynamique, avec Thermoptim. L’application étudiée couple un concentrateur LFR à un moteur Ericsson. L’air est fluide caloporteur et de travail. Un prototype est en construction. L’hybridation et le stockage thermique sont des options clés. Un modèle systémique permettrait d’optimiser l’opération du procédé, en étudiant son comportement dynamique. / Linear Fresnel Reflector (LFR) is a promising Concentrating Solar Power technology. Research is booming and industrial applications are emerging. Applications range from small production units (heat, cold, electricity) to utility scale power plants (several tenths of MW). This PhD work establishes a global model of the solar process, in order to improve our knowledge of the system’s performances. It is a static model suited for a pre-design of the system. The model is chaining three parts. Radiative heat transfer in the optical concentrator is modelled by Monte Carlo statistical Method. The algorithm enables a detailed study of any geometrical configuration, especially through absorbed power flux maps on the receiver. The simulation tool is using the environment of development EDStar. The thermal model calculates analytically the useful thermal power, losses and temperature profiles along the receiver (glass cover, fluid, pipe...). The thermodynamic cycle is simulated analytically using the software Thermoptim. The studied application uses air as heat transfer and working fluid. Air directly feeds an Ericsson engine. The engine developed by LaTEP laboratory is promising for small scale cogeneration (1 to several tenths of kWel). The prototype Linear Fresnel Reflector built in Ecole des Mines d’Albi will enable experimental study of a solar process coupling an LFR concentrator and an Ericsson engine. The technology under study can feed a power plant or a cogeneration system in the industry, producing electricity and heat at 100 to 250°C. Hybridisation with an other energy source (biomass, gas...) and thermal storage (molten salt?) are key features to investigate. To optimise the operating strategy of the process, dynamic behaviour must be studied: a systemic or agent based model is a very relevant approach.

Il existe aujourd'hui un enjeu sociétal majeur à s'intéresser à la propagation du son en milieu extérieur etnotamment, dans notre contexte, à diminuer l'incertitude sur l'estimation des niveaux sonores et améliorer ainsi laprécision des diverses analyses, du bureau d'étude à l'institut de recherche. Dans le cadre de l'acoustiqueenvironnementale, l'influence des conditions météorologiques sur la propagation acoustique en milieu extérieurpeut être importante. Il est donc nécessaire d'appréhender et de quantifier les phénomènes météorologiques demicro-échelles que l'on observe dans la couche limite atmosphérique.Dans ce but, le modèle météorologique de recherche de Météo-France (Meso-NH) a été utilisé. Après avoircomparé les résultats de ce modèle à très fine résolution (de l'ordre du mètre) à l'aide des bases de données de deuxcampagnes expérimentales (Lannemezan 2005 et la Station de Long Terme), il s'est avéré nécessaire de développercet outil en intégrant la prise en compte de la force de traînée des arbres. Dès lors, les résultats issus de Meso-NH surles champs de vent, de température et d'énergie cinétique turbulente aparraissent satisfaisants. Ces informationssont par la suite utilisées en données d'entrée du modèle de propagation acoustique.Le modèle acoustique temporel utilisé est basé sur la méthode Transmission Line Matrix (TLM). Sondéveloppement a été effectué dans le but d'être appliqué à la propagation acoustique en milieu extérieur : prise encompte du relief, de différents types de sol, des conditions atmosphériques, etc. La validation numérique de laméthode TLM, par comparaison avec d'autres modèles (analytique et numérique de type Equation Parabolique), apermis de montrer la pertinence de son utilisation dans le cadre de l'acoustique environnementale.Enfin, à l'aide de ces modèles, des niveaux sonores simulés sous différentes conditions de propagation(favorables, défavorables, homogènes) ont été comparés aux mesures in-situ réalisées lors de la campagneexpérimentale de Lannemezan 2005. Les résultats se sont avérés très satisfaisants au regard de la variabilité desphénomènes observés. Cependant, l'utilisation des champs issus d'un modèle micro-météorologique de type Meso-NH reste délicate du fait de la forte sensibilité du niveau sonore aux profils verticaux de célérité du son. L'étude defaisabilité sur une expérience plus complexe (la Station de Long Terme) est encourageante et, à condition de disposerd'importants moyens de calculs, elle permet de considérer la TLM comme une nouvelle méthode de référence etainsi, d'envisager d'élargir son domaine d'utilisation à d'autres applications. / Actually, it exists a major societal issue to be interested in outdoor sound propagation and specially, in our context, toreduce the uncertainty in noise levels estimation and thus to improve the analyses accuracy, from engineers toresearch institutes. The influence of meteorology on outdoor sound propagation is significant. It is thereforenecessary to understand and quantify the micro scales phenomena into the atmospheric boundary layer.In this way, the French research meteorological model (Meso-NH) has been used. After comparing resultsof this model at very fine resolution (~1 meter) to measurements issued from the databases of two experimentalcampaigns (Lannemezan 2005 and Long Term Monitoring Station: LTMS), it appeared necessary to develop this toolin order to take into account the drag force of the trees. Finally, the Meso-NH results for the wind, temperature andturbulent kinetic energy fields are satisfactory. Then, theses informations can be used as input data for acousticmodels.Our time domain acoustic model is based on Transmission Line Matrix method (TLM). Its development wasdone in order to be applied to outdoor sound propagation: taking into account topography, soil types, meteorologicalconditions, etc. The numerical validation of the TLM method, by comparison with other models (analytical andnumerical: Parabolic Equation), has shown the relevance of its use in the context of environmental acoustics.Finally, thanks to these models, simulated noise levels in different propagation conditions (downward,upward and homogeneous refraction conditions) were compared to in situ measurements carried out during theLannemezan 2005 experimental campaign. Satisfying results were obtained regarding the observed phenomenavariability. However, using the micro-meteorological model Meso-NH is difficult because of the strong acousticsensitivity to the vertical celerity profiles. A feasibility study on a more complex experience (LTMS) is encouragingand, provided having substantial computing resources, it permits to consider the TLM as an accurate method in thecontext of environmental acoustics.

L'objectif de cette recherche est de construire un modèle réduit de petite dimension pour prévoir les réponses dynamiques dans une bande BF sur les parties rigides d'un véhicule automobile complet. Un tel modèle réduit "léger" est une aide à la phase de conception en "Avant Projet" de ces véhicules qui ont la particularité de présenter de nombreux modes élastiques locaux en BF dues à la présence de nombreuses parties flexibles et d'équipements. Pour la construction du modèle réduit, nous avons introduit une base non usuelle de l'espace admissible des déplacements globaux. La construction de cette base requiert la décomposition en sous-domaines du domaine de la structure qui peut présenter une très grande complexité géométrique et dont les modèles EF font intervenir de très nombreux types d'éléments finis. Cette décomposition en sous-domaines a été réalisée par la Fast Marching Method que nous avons due étendre pour pouvoir traiter la complexité des modèles EF des véhicules automobiles. Puis les équations matricielles du modèle EF sont projetées sur cette base. Afin de prendre en compte les incertitudes sur les paramètres du modèle, les incertitudes de modèle induites par les erreurs de modélisation et enfin les incertitudes liées à la non prise en compte des contributions locales dans le modèle réduit des déplacements globaux, un unique modèle probabiliste non paramétrique de ces trois sources d'incertitude a été implémenté sur le modèle réduit construit avec les vecteurs propres globaux. Les paramètres de dispersion de ce modèle probabiliste ont été identifiés en utilisant le principe du maximum de vraisemblance et des réponses obtenues à l'aide d'un modèle stochastique de référence qui inclut des informations expérimentales résultant de travaux précédents. Le modèle réduit stochastique, pour la prévision des déplacements globaux sur les parties rigides dans la bande BF qui a été développé, a été validé sur un modèle de structure automobile "nue" puis a été appliqué avec succès sur un modèle complet de véhicule automobile / The objective of this research is to construct a reduced-order model to predict the dynamical response, in the LF band, of the stiff parts of a complete automotive vehicle in order to facilitate the draft design. The vehicles under consideration have many elastic modes in LF due to the presence of many flexible parts and equipments. To build such a model, we introduced a non-usual basis of the admissible space of global displacements. The construction of this basis requires the decomposition of the domain of the structure. This subdomain decomposition is performed by using the Fast Marching Method that we have extended to take into account the high complexity of the mesh of an automotive vehicle. Then the matrix equations of the FE model are projected on this basis. To take into account the system parameters uncertainties, the model uncertainties induced by the modeling errors and finally, the uncertainties related to the neglecting of local contributions in the reduced-order model, a nonparametric probabilistic model of the three sources of uncertainties has been implemented on the reduced-order model constructed with the global displacements eigenvectors. The dispersion parameters of the probabilistic model are identified using the maximum likelihood method and the responses obtained from a stochastic reference model which includes experimental data resulting from previous works. This stochastic model which has been designed for the prediction of the global displacements of the rigid parts in the LF band is validated on a simple structure of an automotive model and has been successfully applied on a complete model of automotive vehicle

Le transport de la neige par le vent est une composante importante de l'interaction entre l'atmosphère et la cryosphère. En zone de montagne, il influence la distribution temporelle et spatiale de la couverture neigeuse au cours de l'hiver et a en premier lieu des conséquences sur le danger d'avalanche. La modélisation numérique de ce phénomène permet d'étudier les interactions complexes entre le manteau neigeux et le vent et d'en estimer les conséquences de manière distribuée. Dans ce contexte, cette thèse décrit le développement et l'évaluation d'un modèle couplé atmosphère/manteau neigeux dédié à l'étude du transport de la neige par le vent en zone de montagne reposant sur le modèle atmosphérique Meso-NH et le modèle détaillé de manteau neigeux Crocus. Le transport de la neige par le vent a été étudié sur le site expérimental du Col du Lac Blanc (massif des Grandes Rousses, France). Une base de données d'épisodes de transport couvrant dix hivers a tout d'abord été utilisée pour déterminer les caractéristiques principales de ces épisodes. Des simulations avec le modèle Crocus (non couplé à Meso-NH) ont ensuite montré qu'il était nécessaire de tenir compte des transformations mécaniques des grains de neige induites par le vent afin de simuler une évolution réaliste de la vitesse seuil de transport. Le site expérimental a également été le siège de deux campagnes de mesures en 2011 et 2012 visant à collecter de données de validation pour le modèle. Elles renseignent sur les conditions météorologiques près de la surface, sur les quantités de neige transportées et sur la localisation des zones d'érosion et de dépôt de la neige grâce à l'utilisation d'un laser terrestre. Le modèle de transport de neige par le vent Meso-NH/Crocus a été développé. Il intègre le transport de la neige en saltation et en suspension turbulente ainsi que la sublimation des particules de neige transportée. Un schéma à deux moments permet de simuler l'évolution spatiale et temporelle de la distribution en taille des particules. L'utilisation d'un schéma de couche limite de surface à l'interface entre Meso-NH et Crocus s'est révélé nécessaire pour représenter les forts gradients de concentration en particules de neige observés près de la surface. Meso-NH/Crocus est le premier modèle couplé atmosphère/manteau neigeux capable de simuler de manière interactive le transport de la neige par le vent en zone alpine. Meso-NH/Crocus a été évalué en relief réel grâce aux données collectées lors de la première campagne de mesure en 2011. La simulation d'un épisode de transport sans chute de neige simultanée montre que le modèle reproduit de manière satisfaisante les principales structures d'un écoulement en relief complexe ainsi que les profils verticaux de vitesse de vent et de flux de particules de neige en suspension près de la surface. En revanche, la résolution horizontale de 50 m est insuffisante pour reproduire avec précision la localisation des zones d'érosion et de dépôt autour du Col du Lac Blanc. La prise en compte de la sublimation réduit la quantité de neige déposée de l'ordre de 5%.Les techniques de descente d'échelle dynamique (grid nesting) ont ensuite été utilisées pour simuler un second épisode de transport avec chute de neige. L'augmentation de la résolution horizontale intensifie les contrastes de vitesse de vent entre versants au vent et sous le vent. En revanche, elle modifie peu les quantités et les structures spatiales des précipitations solides autour du Col du Lac Blanc. Lorsqu'il est activé, le transport devient la principale source d'hétérogénéités des accumulations neigeuses / Blowing and drifting snow are crucial components of the interaction between the cryosphere and the atmosphere. In mountainous areas, it affects the temporal and spatial distribution of snow depth throughout the winter season and influences avalanche formation. Numerical modeling offers a solution for studying the complex interaction between the snowpack and the wind field and to assess the related processes in a spatially distributed way. In this context, this PhD describes the development and the validation of a coupled snow/atmosphere model which is dedicated to the study of blowing and drifting snow in alpine terrain. The coupled model consists in the atmospheric model Meso-NH and the detailed snowpack model Crocus. Blowing and drifting snow have been monitored at the Col du Lac Blanc (Grandes Rousses range, French Alps) experimental site. A database consisting of blowing snow events observed over 10 years allowed us to identify the main features of these events. Numerical simulations using Crocus illustrated the necessity of taking the wind-dependence of snow grain characteristics into account in order to simulate satisfactorily the occurrence of blowing snow events. We also carried out two measurement campaigns at our experimental site in 2011 and 2012 in order to collect validation data for the model. This includes measurements of vertical profiles of wind speed and snow particle fluxes near the surface and the mapping of areas of erosion and deposition using terrestrial laser scanning. The coupled Meso-NH/Crocus model has been developed in order to account for blowing and drifting snow. It simulates snow transport in saltation and in turbulent suspension and includes the sublimation of suspended snow particles. In the atmosphere, a double-moment scheme allows the model to simulate the spatial and temporal evolution of the snow particle size distribution. The implementation of a surface boundary layer scheme at the interface between Meso-NH and Crocus turned out to be necessary to reproduce the strong vertical gradient of snow particle concentration near the surface. Meso-NH/Crocus is the first coupled snow-atmosphere model that can simulate snow transport in alpine terrain in an interactive way.Meso-NH/Crocus has been evaluated against data collected near Col du Lac Blanc during the first measurement campaign in 2011. The simulation of a blowing snow event without concurrent snowfall showed that the model captures the main structures of atmospheric flow in complex terrain, the vertical profile of wind speed and the snow particle fluxes. However, the horizontal resolution of 50 m is found to be insufficient to simulate the location of areas of snow erosion and deposition observed around Col du Lac Blanc. Blowing snow sublimation leads to a reduction in snow deposition of approximately 5%.We used downscaling techniques (grid nesting) to simulate a second blowing event with concurrent snowfall. The increase in horizontal resolution enhanced the contrast of wind speed between windward and leeward slopes. However, it only slightly affects the amount and the spatial pattern of snow precipitation around Col du Lac Blanc. When activated, blowing and drifting snow are the main sources of spatial variability of snow accumulation

Les volcans bouclier basaltiques se déforment sous leur propre poids, tout particulièrement en période d'éruption. À Hawaii, cette déformation co-éruptive est en majeure partie expliquée par un plan de décollement à la base des édifices, dont le glissement serait associé à la force de poussée des intrusions dans les «rift zones» (i.e. concentration de dykes subverticaux dans les zones de faiblesse de l'édifice). Cependant, cette explication n'est pas valable pour beaucoup d'autres volcans basaltiques, où l'existence d'un décollement est peu probable. Nous avons profité de l'incision intense du Piton des Neiges (La Réunion) par l'érosion pour observer la structure interne d'un volcan basaltique, et comprendre, par le terrain et la modélisation numérique, comment la déformation d'un tel édifice s'articule avec son activité magmatique. Notre étude structurale révèle que trois populations d'intrusions perpendiculaires (deux subverticales et une subhorizontale) coexistent, et que les principales zones de faiblesses sont des «sill zones» (i.e. concentrations d'intrusions subhorizontales), Parallèlement, notre étude des déformations, essentiellement cassantes, révèle que deux extensions perpendiculaires dominent dans l'édifice, avec l'apparition de régimes localement décrochants et même compressifs à proximité des sill zones. Les directions de déformations et les directions d'intrusions sont compatibles et coïncident aussi avec la direction d'écoulement des grandes avalanches de débris qui ont affectée le Piton des Neiges au cours de son histoire. À partir de ces données de terrain, nous proposons un modèle de déformation du Piton des Neiges où les intrusions dans l'édifice génèrent des cycles de permutations de contraintes. Le stade ultime de chaque cycle serait la mise en place de sills en régime compressif. En testant numériquement l'effet de telles injections sur la déformation d'un édifice, nous montrons que les sills peuvent être un facteur majeur d'instabilité, capable de conduire à des déstabilisations de flanc telles que les avalanches observées. Notre modèle conceptuel, déduit du terrain et quantifié par la modélisation, constitue ainsi une alternative au modèle de déformation hawaiien, applicable sur des édifices tels que le Piton de la Fournaise (La Réunion), Tenerife (Canaries) et Fernandina (Galapagos). Notre étude démontre enfin l'intérêt essentiel d'étudier les systèmes anciens pour mieux comprendre les volcans actifs. / Basaltic shield volcanoes deform under their own weight, especially during eruptive periods. In Hawaii, this co-eruptive deformation is mainly explained by slip events on a basal decollement plane, related to the forceful intrusion of magma into “rift zones”(i.e. high concentration of subvertical dikes in weakness areas). However, this explanation is not valid on many other basaltic shields where the existence of basal decollements is unlikely. We used the deep incision of Piton des Neiges volcano (La Réunion) to observe the internalstructure of a basaltic shield.By coupling a field work and a numerical study, we aimed at understand how deformation interacts with magmatic activity on such an edifice. Our structural study reveals that three populations of perpendicular intrusions coexist and that the main weakness areas are “sill zones” (i.e. high concentration of low-­dipping intrusions). In parallel, our study of brittle deformation structures shows that two perpendicular extensional stress fields dominate in the edifice and that strike-­slip and compressional regimes appear near sill zones. The directions of deformation are compatible with the orientations of intrusions and are also consistent with the directions of emplacement of large debris avalanches, which occurred on Piton des Neiges during its evolution. From these field data, we propose a deformation model of Piton des Neiges volcano where magma intrusion in the edifice generates cycles of stress permutations. The ultimate stage of each cycle consists in sill intrusion under a compressional regime. The numerical simulations, testing the influence of such injections on edifice deformation, reveal that sills can be major instability factors, capable of triggering large flank collapses. Our conceptual model, inferred from field work and quantified by numerical models, thus constitutes an alternative to the Hawaiian model of deformation, applicable on edifices like Piton de la Fournaise (La Réunion), Tenerife (Canary) or Fernandina (Galapagos). Our study finally demonstrates the essential interest of studying old eroded systems to understand active volcanoes.

Cette thèse porte sur l’étude de l’évolution de la qualité de l’air en Europe et en France dans les prochaines décennies à l’aide de simulations numériques. Dans les études des impacts de l’évolution du climat sur la qualité de l’air, les modèles de chimie atmosphérique utilisent des sorties de modèles climatiques globaux ou régionaux qui fournissent les « forçages », c’est-à-dire les conditions météorologiques simulées pour les périodes futures. Contrairement aux analyses météorologiques, qui représentent la variabilité jour à jour du temps, les sorties des modèles de climat doivent nécessairement être interprétées de manière statistique : elles ne représentent la météorologie que dans un sens climatologique. Afin de pouvoir commenter utilement les simulations futures de qualité de l’air, il est nécessaire d’évaluer au préalable et pour le climat présent, la qualité des simulations calculées avec des forçages climatiques par rapport aux références que constituent les simulations calculées avec des forçages analysés et, bien entendu, les observations. Trois simulations pluri-annuelles (6 ans) ont été lancées pour la période actuelle (2000-2010) et ont été comparées ; elles différent par l’utilisation d’analyses météorologiques ou de forçages de modèle de climat (pour les paramètres atmosphériques seuls et par ailleurs, pour les paramètres atmosphériques et le calcul des échanges en surface) en entrée du modèle de chimie-transport tridimensionnel de Météo-France, MOCAGE. Nous avons évalué ces différentes simulations par comparaison aux observations de la base de données européenne AirBase. Nous avons ensuite comparé les performances de ces simulations pour un grand nombre de scores quantitatifs, en analysant d’une part les effets liés aux champs météorologiques (température, vent, humidité, etc.) et d’autre part, ceux liés aux échanges en surface (comme les vitesses de dépôts, les émissions biogéniques) qui dépendent également de la météorologie. Nous avons ainsi évalué comment ces changements affectent les distributions horizontales et verticales des polluants. In fine, nous avons caractérisé la fiabilité des simulations de qualité de l’air reposant sur des forçages issus de modèles climatiques pour le climat présent : des indicateurs (biais moyens, biais moyens normalisés, RMSE, déviations standards) et des index de qualité de l’air (comme le dépassement de seuils) se distinguent et peuvent donc servir de base fiable pour l’interprétation des résultats pour les simulations du futur. Enfin dans une troisième partie, ces indicateurs considérés comme pertinents ont été utilisés pour étudier des simulations de qualité de l’air aux horizons 2030 et 2050 (5 ans). Comme attendu, l’évolution des paramètres météorologiques (température, précipitation, vent) modifie les quantités et la dispersion des polluants dans l’atmosphère, mais l’évolution des émissions en Europe et dans le reste du monde joue aussi un rôle important. Ainsi, face à l’évolution du climat et la hausse des émissions dans certains pays du monde, en Asie notamment, les effets des politiques Européennes pour réduire les émissions anthropiques sont mitigés selon les régions et les polluants, dépendant de l’influence relative des phénomènes locaux et du transport de polluant à longue distance. / This thesis aims at predicting how European and French air quality could evolve over the next decades using numerical modeling. In order to study the impacts of climate change up on regional air quality, atmospheric chemistry models rely on global or regional climate models to produce “forcings”, i.e. meteorological conditions for future periods. Unlike meteorological analyses, which can represent specifically each date and hour thanks to the assimilation of observations, climate model outputs need to be averaged and can only be interpreted in a climatological sense. And so are air quality hindcasts relying on them for their forcings. In order to properly interpret air quality simulations in a future climate, it is a pre-requisite to assess how realistic air quality hindcasts are when driven by forcings from climate models for the current period in comparison to the references, which are simulations with the same set-up but relying on meteorological analyses and observations. Three six-year simulations for the current climate (2000-2010) have been run with the three-dimensional chemistry transport model of Météo-France, MOCAGE. These simulations only differ by the meteorological forcings used, either operational meteorological analyses or outputs from climate simulations (for atmospheric parameters only ; for atmospheric parameters as well as surface exchanges, which depend also on the weather). We compared the three simulations and evaluated them against the European air quality database of the European Environment Agency, AirBase. Further, we investigated how statistical skill indicators compare in the different simulations, assessing the effects of meteorology on atmospheric fields (temperature, wind, humidity,...) and on the dependent emissions and deposition processes (such as deposition velocities, volatile organic compound emissions, ...) that depend upon meteorology. We have in particular studied how these factors affect the horizontal and vertical distributions of species. In the end, we have estimated how reliable are skill indicators for the simulations run with “climate” forcings : some indicators (mean bias, mean normalized bias, RMSE, deviation standards, number of exceedance days) are sufficiently close to the ones obtained with the reference configuration (relying on analysed meteorological forcings) to be considered reliable. They can thus be used to interpret simulations for future periods. We have run simulations of European air quality in the 2030s and 2050s (5 years for each period). They are discussed using the indicators previously indentified. As expected, the changes in meteorological parameters (temperature, precipitation, wind, ...) affect the quantities and distributions of pollutants in the atmosphere, but the future evolutions in European and global emissions also play a significant role. Faced with climate change and increased emissions in some countries in the world, as in Asia, the impacts of European policies for reducing anthropogenic emissions are mitigated, depending on the regions and the pollutants due to the respective influence of local emission and of long-range transport of pollutants.

Le secteur des transports représente 50% de la consommation mondiale en produits pétroliers. Une part significative de l'énergie est dissipée par le fonctionnement de la suspension et par la déformation du pneu. La texture de surface de chaussée et son interaction avec le pneu jouent un rôle certain dans cette dissipation d'énergie. Le travail présenté dans ce mémoire est consacré à la modélisation du contact entre un pneu et une chaussée en conditions de roulement dans le but d'évaluer la résistance au roulement. Le choc d'une sphère sur un demi-espace viscoélastique est d'abord étudié. Ensuite, le contact viscoélastique d'un objet roulant sur une surface lisse est modélisé avec une suspension attachant l'objet roulant à une masse qu'il supporte. Enfin, une méthode multipoint linéarisée est proposée pour l'étude du roulement sur une surface rugueuse sans suspension dans un premier temps, et avec suspension dans un second temps. La condition de roulement est introduite. Les surfaces rugueuses sont soit des surfaces composées d'aspérités sphériques, soit des surfaces de chaussées réelles. Les résultats numériques font apparaître une dissymétrie dans la distribution de pression de contact dans la direction du roulement. La force de résistance au roulement, définie comme le couple par rapport à l'axe de l'essieu divisé par le rayon de la roue, est calculée à partir de la distribution de pression. L'influence de différents paramètres sur cette force tels que la vitesse de roulement, la taille des aspérités ainsi que la raideur de la suspension est analysée / The transport accounts for about 50 percent of the fuel consumption. A significant proportion of energy is lost by the suspension and the deformation of the tyre. One of the important factors in the energy loss is the texture of the road and its interaction with the tyre. The present work deals with the numerical modelling of the contact between a tyre and a road in rolling conditions in order to evaluate the rolling resistance. First, the collision of a sphere and a viscoelastic halfspace is studied. Then the viscoelastic contact of a body rolling on an smooth surface is investigated with a suspension between a mass and the rolling body. Finally a linearised multipoint method is proposed for the study of the rolling contact on a rough surface first without suspension and then with a suspension. The rolling condition is introduced in the model. The rough surfaces are surfaces composed of spherical asperities or real road surfaces. Numerical results reveal an asymmetry in the distribution of the contact pressure in the rolling direction. The rolling resistance force, classically defined as the momentum opposing the rolling divided by the radius of the tyre, is calculated using the distribution of the contact pressure. The effects of different parameters on the rolling resistance force such as the rolling speed, the size of the asperities and the stiffness of the suspension are analysed