Cette thèse adapte le modèle GRAAL, initialement dédié à la prédiction du comportement à long terme des verres nucléaires type R7T7, à l’étude de l’altération des verres basaltiques, matériaux considérés comme des analogues naturels des verres nucléaires. Cette thèse cherche à valider les hypothèses fondatrices du modèle GRAAL et évaluer sa robustesse en l’appliquant aux verres basaltiques. Pour ce faire, deux axes de recherche sont suivis. Le premier consiste à reparamétrer le modèle GRAAL pour le verre basaltique. Le travail de description de la couche amorphe qui se forme durant l’altération du verre, a été entrepris aux moyen d’expériences de lixiviation de verres basaltiques simplifiés à 90°C et différents pH (7 ; 8,5 et 10). Des fortes analogies observées sur le comportement des éléments ont permis de conserver une partie du paramétrage du modèle GRAAL dédié aux verres nucléaires. Des expériences ont également été conduites pour affiner la compréhension des mécanismes d’altération à long terme du verre basaltique. A fort progrès de réaction, l’altération du verre basaltique n’est pas limitée par un mécanisme diffusif, comme c’est le cas pour le verre nucléaire, mais par une très lente hydrolyse du réseau vitreux entretenue par la précipitation de phases secondaires. Ainsi une modification de la loi cinétique de dissolution du verre a été réalisée.Le modèle, reparamétré avec la même méthodologie que celle employée pour le paramétrer pour le modèle de verre nucléaire parvient à rendre compte avec justesse de l’altération du verre basaltique pour les pH faiblement basiques, renforçant ainsi la confiance que l’on peut accorder au modèle sur les prédictions liées au verre nucléaire. / This thesis adapts the GRAAL model, initially dedicated to the prediction of the long-term behavior of nuclear glasses, to basaltic glasses.This thesis intends to validate the relevance of the founding hypotheses of the GRAAL model as well as to assess its robustness. To achieve this, two lines of research were carried out simultaneously. The first was to re-parameterize the GRAAL model for basaltic glass compositions. To describe the amorphous layer formed during glass alteration, simple basaltic glasses were leached at 90 ° C and various pH (7, 8.5 and 10). Strong analogies observed on the behavior of the elements allowed to preserve a part of the parameterization of the GRAAL model. Experiments have also been conducted to improve the understanding of the mechanisms of basaltic glasses alteration in silica saturated conditions. It is shown that alteration of the basaltic glass is not limited by a diffusive mechanism like nuclear glass, but by a very slow hydrolysis of the vitreous network sustained by the precipitation of secondary phases. Thus, a modification of the kinetic law of the model controlling the dissolution of the glass was carried out.The model, then, accurately accounts for the alteration of basaltic glass in slightly alkaline conditions, thus reinforcing the confidence that can be given to the model on the predictions related to nuclear glasses.

Cette thèse est dédiée à la modélisation et à la compréhension du fonctionnement de dispositifs organiques et plus particulièrement des diodes organiques. Tout d'abord, nous avons modélisé et analysé les caractéristiques courant-tension en mode statique de dispositifs planaire (diodes) ou longitudinaux (barreaux résistifs) à base de TIPS-pentacene. Nous avons pu ainsi expliquer l'origine du redressement en courant pour ce type de dispositifs à base de semi-conducteur non-dopé. L'aspect dynamique a ensuite été pris en compte, avec la proposition d'un modèle numérique et d'un modèle analytique simplifié du comportement en temporel de la diode, ce qui a permis de déterminer quels sont les paramètres physiques impactant la fréquence de coupure. Enfin, compte tenu de l'intérêt grandissant des photodiodes organiques, nous nous intéresserons également dans la dernière partie à la modélisation de la diode sous flux lumineux, pour des applications en tant que cellules solaire ou des photodétecteurs. Nous étudierons dans cette dernière partie l'impact des paramètres électriques et optiques du dispositif (épaisseur, mobilité, taux de dissociation des excitons, indice optique, etc.) sur les courants d'obscurité et sous illumination, et donc les rendements de conversion. / This report is dedicated to the modelling and understanding of organic devices functionement and more particularly to organic diodes. First, we modelled and analysed I-V caracterisations of planar (diodes) and longitudinal (resistive sticks) devices made with TIPS-pentacene. We could explain the reason of the rectifying behavior of devices based on undoped semiconductor. The dynamic aspect was then analysed with the help of a numerical model which was simplifyed in a semi-analytical one. This allowed us to find parameters who have an impact on the cutoff frequency. Finally, as there is a growing interest on organic photodiodes whose applications are organic photodetectors or solar cells, we modellised, in the last part, the impact of light on the functionnement of organic diodes. We studyed the impact of both electrical parameters and optical ones (thicknesses, mobility, dissociations rate of excitons, optical indexes, etc...) on dark current and under illumination.

En zone de montagne, et particulièrement dans les Alpes du Nord, l'amélioration des modèles de prévision de crues est un enjeu scientifique et opérationnel important, du fait notamment des changements climatiques et environnementaux qui modifient le fonctionnement des systèmes hydrologiques. Les apports liés à la fonte nivale constituent une difficulté spécifique à la modélisation hydrologique dans les Alpes. Cette thèse entend contribuer à améliorer les modèles dans ce domaine, à deux niveaux. Dans un premier temps un modèle de fonte de neige performant à l'échelle ponctuelle a été explicité, utilisant un minimum de variables météorologiques. Pour cela l'approche dite « degré jour », qui permet d'estimer en fonction de la température de l'air la lame d'eau déstockée par le manteau neigeux, a été reprise. Grâce aux données acquises par le Centre d'Etude de la Neige sur le site expérimental du Col de Porte (Chartreuse), plusieurs formulations ont été testées prenant notamment en compte le stock de froid du manteau neigeux et les effets liés au cycle saisonnier, à travers un coefficient de fonte explicité à l'aide d'une sinusoïde. Dans un deuxième temps, l'accent est mis sur la régionalisation de la donnée d'entrée fondamentale pour ce modèle degré jour : la température de l'air. La spatialisation des champs de température a été étudiée à l'aide d'outils géostatistiques, en se focalisant sur une meilleure compréhension des gradients altimétriques de température. L'objectif est de mieux cartographier les températures au pas de temps infra journalier et de déterminer un modèle d'interpolation pertinent et simple. Les perspectives sur l'utilisation des résultats obtenus dans ces deux phases de modélisation sont abordées en vue d'une utilisation future des modèles degré jour sur des bassins versants de montagne. / In mountain areas, and especially in the northern Alps, the improving of flood forecasting models is an important scientific and operational issue, due to environmental and climatic changes which modify hydrological system. Snow melt contributions are a specific difficulty in the hydrological modelling in the Alps. This thesis is concerning in improving snow melt modelling for hydrological forecasting. In a first step a snow melt model at a local scale is formulated, using a minimum of meteorological variables. For this the “degree day” model, which allows predicting the water release by snow function of air temperature, is used. Thanks to data of the Centre d'Etude de la Neige on the experimental site of Col de Porte (Chartreuse), several formulations have been tested with some interesting concepts, like the cold content of snow and the effect of seasons with a melt factor given with a sinusoidal expression. In a second step, we are interested in the regionalization of the most important parameter for the snow melt model: air temperature. The spatial interpolation of temperature in mountainous areas is done with geostatistical tools, focusing on a better understanding of the elevation lapse rate of temperature. The goal is mapping air temperatures at a sub daily time step and to propose a simple and relevant model. Prospects about these results of snowmelt modelling and temperatures interpolation are given for using them in mountainous catchments.

Le transfert radiatif joue un rôle important en combustion turbulente et doit donc êtrepris en compte dans les simulations numériques. Toutefois, à cause du fait que la combustionet le rayonnement sont deux phénomènes physiques très différents caractérisés par deséchelles de temps et d’espace également différentes, et la complexité des écoulements turbulents,l’effet du rayonnement est souvent négligé ou modélisé par des modèles très simples.Le couplage entre la combustion (LES) et le rayonnement avec l’environnement CORBAa été étudié. Dans le présent travail, quatre formulations de la méthode de Monte Carlo(méthode classique et méthode réciproque) dédiées à la résolution de l’équation de transfertradiatif ont été comparées sur un cas test de flamme 1D où l’on tient compte de l’absorptionet de l’émission du milieu en utilisant un maillage 3D. Le but de ce cas test est de valider lesolveur Monte Carlo et de choisir la méthode la plus efficace pour réaliser le couplage. Afind’améliorer la performance du code de Monte Carlo, deux techniques ont été développées.De plus, un nouveau code dédié au couplage a été proposé. Ensuite, deux solveurs radiatifs(Emission Reciprocity Monte Carlo Method et Discrete Ordinate Method), appliquésà une flamme turbulente stabilisée en aval d’un dièdre avec un modèle CK de propriétésradiatives, sont comparés non seulement en termes de description physique de la flamme,mais aussi en terme de performances de calcul (stockage, temps CPU et efficacité de laparallélisation). Enfin, l’impact de la condition limite a été discuté en prenant en comptel’émissivité et la température de paroi. / Radiative transfer plays an important role in turbulent combustion and should be incorporatedin numerical simulations. However, as combustion and radiation are characterized bydifferent time scales and different spatial and chemical treatments, and the complexity of theturbulent combustion flow, radiation effect is often neglected or roughly modelled. Couplinga large eddy simulation combustion solver and a radiation solver through a dedicated languageCORBA is investigated. Four formulations of Monte Carlo method (Forward Method,Emission Reciprocity Method, Absorption Reciprocity Method and Optimized ReciprocityMethod) employed to resolve RTE have been compared in a one-dimensional flame testcase using three-dimensional calculation grids with absorbing and emitting medium in orderto validate the Monte Carlo radiative solver and to choose the most efficient model forcoupling. In order to improve the performance of Monte Carlo solver, two techniques havebeen developed. After that, a new code dedicated to adapt the coupling work has beenproposed. Then results obtained using two different RTE solvers (Reciprocity Monte Carlomethod and Discrete Ordinate Method) applied to a three-dimensional turbulent reactingflow stabilized downstream of a triangular flame holder with a correlated-k distributionmodel describing the real gas medium spectral radiative properties are compared not onlyin terms of physical behavior of the flame but also in computational performance (storagerequirement, CPU time and parallelization efficiency). Finally, the impact of boundary conditionstaking into account the actual wall emissivity and temperature has been discussed.

La biomasse lignocellulosique est considérée comme une ressource de carbone renouvelable possédant un grand potentiel pour la une valorisation énergétique et chimique. La torréfaction de la biomasse sèche (de type bois) est un procédé de transformation thermique s’effectuant à des températures comprises entre 200°C et 300°C, et un temps de séjour compris entre quelques minutes et plusieurs heures, opérant sous pression atmosphérique et en défaut d’air. Le produit principal de la torréfaction est un combustible solide hydrophobe et stable. Cette opération génère des coproduits gazeux à haute température qui sont habituellement considérés comme des effluents pénalisants, apportant au mieux un appoint d’énergie pour le procédé. Pourtant, de nombreux constituants présents dans les condensats – récupérés par condensation des coproduits de torréfaction – pourraient être valorisés comme produits chimiques bio-sourcés. L'objectif de la thèse est de proposer un procédé de séparation-purification pour les composés condensables présents dans ces effluents gazeux. Ces condensats constituent une phase majoritairement aqueuse, contenant plus de 150 espèces organiques identifiées. Les espèces minoritaires sont présentes dans des proportions variables suivant le bois torréfié. Enfin il s’agit d’un mélange réactif et thermiquement instable, où différents équilibres chimiques sont présents. Une analyse des caractéristiques physico-chimiques des condensats a permis de proposer un milieu modèle limité à une dizaine de composés. Un modèle représentatif du comportement thermodynamique de ce mélange réactif à large spectre de masse moléculaire a été sélectionné et les paramètres d’interaction binaire identifiés. Des données expérimentales d’équilibres liquide-vapeur ont été acquises pour valider en partie ce modèle. Les composés cibles et les objectifs du procédé de valorisation ont été choisis et plusieurs stratégies de valorisation ont été élaborées et simulées sous Prosim+ sur la base de la modélisation thermodynamique. Cette étude a permis d’évaluer ces différentes stratégies en termes d’efficacité énergétique et de pureté des produits finaux pour une potentielle mise en place à l’échelle industrielle de cette filière.

Au cours d’un transitoire d’Accident de Perte de Réfrigérant Primaire (APRP), la vaporisation de l’eau dans la cuve d’un réacteur nucléaire conduit à l’assèchement et à l’augmentation de la température des crayons de combustible, à leur gonflement thermomécanique et à la rupture des gaines de ces crayons remplies de matériaux combustibles fragmentés. Ce gonflement peut être à l’origine d’un bouchage significatif d’une partie du coeur et compromettre ainsi son refroidissement lors de la phase de renoyage où l’eau injectée via les systèmes de sécurité pénètre dans le bas du coeur. Lors de cette phase de renoyage, la vaporisation importante de l’eau produit un fort débit de vapeur qui arrache et entraîne des gouttelettes dans les zones bouchées. L’enjeu de cette thèse est d’évaluer les capacités de refroidissement des crayons ballonnés par l’impact des gouttes d’eau. Pour ce faire, une étude expérimentale est réalisée au LEMTA de Nancy, afin de mesurer le flux de refroidissement d’un train de gouttes impactant une paroi chaude. L’influence des caractéristiques dynamiques des gouttes sur le flux échangé entre la goutte et la paroi est expérimentalement étudiée. Ces essais ont permis la réalisation d’une base de donnée permettant de valider un modèle de flux goutte/paroi Ce modèle d’échange goutte/paroi est réalisé à partir d’une modélisation des paramètres clefs gouvernant l’échange, à savoir, le diamètre d’étalement de la goutte sur la paroi, l’épaisseur de la couche de vapeur créée sous la goutte et le temps de séjour de la goutte sur la paroi. Ce modèle, validé expérimentalement, est intégré dans le code NEPTUNE_CFD dans le but de réaliser des simulations de refroidissement d’un assemblage de crayon, en phase de renoyage, en prenant en compte l’effet des gouttes impactantes / In a pressurized water reactor (PWR), during a Loss Of Coolant Accident (LOCA), liquid water evaporates and the fuel assemblies are not cooled anymore; as a consequence, the temperature rises to such an extent that some parts of the fuel assemblies can be deformed resulting in ’ballooned regions’. When reflooding occurs, the cooling of these partially blocked parts of the fuel assemblies will depend on the coolant flow that is a mixture of overheated vapour and undersaturated droplets. The aim of this thesis is to study the heat transfer between droplets and hot walls of the fuel rods. In this purpose, an experimental device has been designed in accordance with droplets and wall features (droplet velocity and diameter, wall temperature) representative of LOCA conditions. The cooling of a hot Nickel disk, previously heated by induction, is cooled down by a stream of monodispersed droplet. The rear face temperature profiles are measured by infrared thermography. Then, the estimation of wall heat flux is performed by an inverse conduction technique from these infrared images. The effect of droplet dynamical properties (diameter, velocity …) on the heat flux is studied. These experimental datas allow us to validate an analytical model of heat exchange between droplet and hot slab. This model is based on combined dynamical and thermal considerations. On the one hand, the droplet dynamics is considered through a spring analogy in order to evaluate the evolution of droplet features such as the spreading diameter when the droplet is squeezed over the hot surface. On the other hand, thermal parameters, such as the thickness of the vapour cushion beneath the droplet, are determined from an energy balance. In the short term, this model will be integrated in a CFD code (named NEPTUNE_CFD) to simulate the cooling of a reactor core during a LOCA, taking into account the droplet/wall heat exchange

La mixotrophie (la combinaison de la phototrophie et de l’hétérotrophie chez un même organisme) est une stratégie nutritive très courante chez les protistes planctoniques et son importance pour le fonctionnement des écosystèmes aquatiques a récemment été mise en évidence. La plupart des modèles biogéochimiques décrivant les écosystèmes aquatiques sont cependant encore basés sur le paradigme traditionnel d’une séparation stricte entre les protistes phototrophes (le phytoplancton) et les protistes hétérotrophes (le nano/microzooplancton).Dans la première partie, la thèse de doctorat propose deux modèles mathématiques mécanistiques permettant de décrire l’ensemble des stratégies nutritives déployées par les protistes planctoniques et pouvant être implémentés dans les modèles biogéochimiques. Le premier modèle décrit la capacité des protistes phototrophes à utiliser le phosphore organique dissous (DOP) lorsqu’ils sont limités en phosphate. Le second modèle décrit la mixotrophie stricte (phototrophie + phagotrophie) et permet de représenter les grands types de mixotrophies observés dans les milieux aquatiques: la mixotrophie constitutive et la mixotrophie non-constitutive.Dans la seconde partie du travail, les modèles mécanistiques sont implémentés dans le modèle MIRO (Lancelot et al. 2005) décrivant l’écosystème planctonique de la Baie Sud de la Mer du Nord, avec l’objectif d’étudier le rôle de la mixotrophie sur le fonctionnement d’un écosystème côtier eutrophisé.Les résultats montrent que lorsque les apports en nutriments par les rivières (Seine et Escaut) sont élevés et enrichis en N par rapport au P, le DOP devient une source importante en P pour les protistes phototrophes (surtout pour les colonies de l’espèce nuisible Phaeocystis) ;ne pas en tenir compte mènerait à une sous-estimation de la production primaire (jusqu’à 32%), de la sédimentation (jusqu’à 20%), de la reminéralisation en ammonium (jusqu’à 13%), des émissions de diméthylsulfure (DMS ;jusqu’à 97%), et de la capacité de l’écosystème à absorber le CO2 atmosphérique (avec, certaines années, un changement du rôle de l’écosystème, passant de source à puits de CO2). La production secondaire n’est, par contre, pas affectée par le DOP. Les mixotrophies constitutive et non-constitutive ont en général peu d’impact sur la production primaire, la production secondaire, la sédimentation, et la reminéralisation dans l’écosystème étudié. L’effet de la mixotrophie non-constitutive dépend cependant de la lumière disponible ;les printemps/été caractérisés par des intensités lumineuses élevées, la mixotrophie non-constitutive induit une augmentation de tous les processus étudiés, en particulier la production secondaire (jusqu’à 23%). / Doctorat en Sciences agronomiques et ingénierie biologique / info:eu-repo/semantics/nonPublished

Le magnésium et ses alliages sont de plus en plus utilisés dans les transports pour leur faible densité (1,7), et on cherche à les utiliser sous la forme corroyée. De structure hexagonale compact, le magnésium est très anisotrope et faiblement ductile à froid. La déformation est contrôlée par l'activation de systèmes de glissement et de maclage, seuls le glissement basal et le maclage de traction sont actifs à froid, alors qu'à chaud sont activés les glissements prismatique et pyramidal <c+a>, et le maclage de compression. Nous avons élaboré des monocristaux de magnésium pur de différentes orientations pour mesurer les contraintes de cission résolues critique (CRSS) des différents systèmes par des essais de compression plane à des températures allant de l'ambiante à 450°C. En parallèle nous avons développé un code de plasticité cristalline pour modéliser la déformation de cristaux contraints, ce programme est basée sur la loi de Schmid et le principe du travail maximum, ce qui permet de prédire les systèmes de déformation activés et la contrainte nécessaire, en fonction des CRSS. Nous avons aussi testé le modèle pour prédire le comportement de l’alliage polycristallin AZ31. / Magnesium and its alloys are increasingly used in the transport industry for their low density (1.7), and ideally could be used after plastic shaping. Due to the hcp structure, magnesium is very anisotropic and exhibits low ductility at room temperature. Strain is controlled by the activation of slip and twinning mechanisms, and only basal slip and tensile twinning operate at RT, whereas at high temperature prismatic and pyramidal <c+a> slip, and compressive twinning are active. Pure magnesium single crystals were produced with different orientations to measure the critical resolved shear stresses (CRSS) of the different deformation mechanisms, by channel-die compression from RT to 450°C. In parallel we developed a crystal plasticity model to simulate the deformation of constrained crystals, this model is based on the Schmid law and the principle of maximum work. The model enables the prediction of activated deformation mechanisms and the stress as a function of the CRSS. We have also used the model to predict the behaviour of a polycrystalline alloy AZ31.

L'étude de la planète Vénus est particulièrement intéressante car elle possède de nombreuses caractéristiques avec la Terre, des dimensions et une composition semblables. Malgré ces éléments en commun, ces planètes sont foncièrement différentes à notre époque. Vénus possède une atmosphère cent fois plus dense à sa surface que celle de la Terre. Cette atmosphère est composée essentiellement de dioxyde de carbone là où l'atmosphère terrestre est composée d'un mélange de diazote et d'oxygène moléculaire. La composition chimique de l'atmosphère vénusienne pose de nombreux problèmes aux chimistes. Vénus fait l'objet d'observations modernes depuis la Terre, de missions spatiales depuis les années 70 et l'ESA a lancé en 2005 la mission Venus Express qui est restée en orbite autour de Vénus jusqu'en décembre 2014. L'explication de la distribution et des concentrations des espèces chimiques mesurées est étudiée par une communauté de modélisateurs en parallèle des observations. Cependant, ces modélisations étaient contraintes à utiliser des modèles 1D qui contenaient de nombreuses paramétrisations. Durant ma thèse j'ai mis à jour un modèle chimique, développé au LATMOS, et développé un modèle des nuages de Vénus qui ont tout deux été introduit à un GCM 3D de Vénus développé au LMD. Ce nouveau modèle permet d'appliquer en 3D les hypothèses les plus récentes sur les cycles chimiques vénusiens, tels que la décroissance du dioxyde de soufre et de la vapeur d'eau dans la couche nuageuse, les variations en latitudes du monoxyde de carbone, la stabilité du dioxyde de carbone, la composition des gouttelettes des nuages de Vénus et la masse condensée dans les nuages. / The study of Venus is particularly interesting on many aspects. Venus shares many features with the Earth, similar size and composition. For, despite some common elements, these planets are fundamentally different in our time. Venus has an atmosphere a hundred times denser at the surface than the Earth and is composed essentially of carbon dioxide where the earth's atmosphere is composed of a mixture of molecular nitrogen and molecular oxygen. The chemical composition of the Venusian atmosphere poses many problems for chemist. Venus is the subject of modern observations from Earth and then through space missions of the 70 ESA mission launched in 2005 which remained in orbit around Venus until December 2014. the explanation of the distribution and concentrations measured chemical species is studied by a community of modelers in parallel observations. However, these models were forced to use 1D models that contained numerous parameterization. During my PhD I have updated a chemical model developed at LATMOS, and developed a model of the Venus clouds which have both been introduced to Venus 3D GCM developed at LMD. It allows to apply in 3D the latest hypothesis about the Venusian chemical cycles such as, among others, the decrease of sulfur dioxide and water vapor in the cloud layer, changes of carbon monoxide in latitudes, the stability of carbon dioxide, the composition of clouds droplets and the condensed mass in those clouds.

La leishmaniose est endémique dans plus de 98 pays avec une prévalence totale de 12 millions de cas. L’incidence annuelle de la leishmaniose viscérale est estimée à 0,3 million de cas, tandis que la leishmaniose cutanée et muco-cutanée affecte un million de personnes par année. La leishmaniose cutanée est un sérieux problème de santé publique en Tunisie, avec des milliers de cas répertoriés chaque année. Depuis sa première émergence en 1982 à Kairouan, la maladie s’est propagée vers le centre et le sud du pays. Cependant, la surveillance de la maladie est limitée à la déclaration de nombre de cas par district et les études tentant d’élucider la dynamique spatio-temporelle de la maladie sont basées uniquement sur l’analyse de séries temporaires de l’incidence de la maladie. Dans le contexte méditerranéen, la distribution des phlébotomes, vecteurs de la leishmaniose, n’a cessé de changer. Durant la dernière décade, les vecteurs de la leishmaniose ont gagné des zones géographiques plus au nord. En effet, les phlébotomes jadis présents à des latitudes au-dessous de 45° N ont été répertoriés à des latitudes au-dessus de 50° N. Des études récentes ont montré la présence de phlébotomes dans diverses régions de l’Allemagne et de la Belgique. En Europe centrale, des cas autochtones de leishmaniose sont en train d’être enregistrés dans des régions traditionnellement considérées comme non endémiques. L’objectif de ce travail, dans un premier temps, est de i) prédire la répartition géographique de cas de leishmaniose cutanée et de son principal vecteur P. papatasi, en Tunisie. ii) déterminer l’importance des facteurs environnementaux qui influencent la distribution du vecteur et de la maladie en utilisant la modélisation de la niche écologique. Par la suite, nous avons essayé de i) prédire la distribution actuelle de six phlébotomes vecteurs de la leishmaniose dans le Bassin Méditerranéen. ii) estimer le changement de la distribution géographique de ces vecteurs sous différents scénarios de changement climatique. Pour prédire la distribution géographique de P. papatasi et des cas de maladies en Tunisie la modélisation de la niche écologique a été réalisée en utilisant le modèle MaxEnt. Nous avons collecté 24 variables environnementales décrivant la topographie, l’utilisation du sol, et le climat en tant que variables explicatives. Une étude sur le terrain a été réalisée afin de déterminer les points de présence de P. papatasi et des cas de leishmaniose cutanée servant de variables d’entrées au modèle. Le modèle MaxEnt a montré une bonne performance de prédiction de la répartition du vecteur et de cas de maladie avec une aire au-dessous de la courbe ROC (Receiver Operating Characteristic) supérieure à 0,95 et un coefficient Kappa maximal supérieur à 0,7. De plus, la population sous le risque de leishmaniose cutanée a été évaluée en se basant sur la densité de population par district et la répartition du vecteur de la maladie prédite par le modèle. Afin, d’étudier l’effet du changement climatique sur la réparation des différents vecteurs de la leishmaniose dans le bassin méditerranéen, des points de présence des différents vecteurs ont été collectés à partir de la littérature scientifique entre 2006 et 2012. De plus, une série de variables climatiques décrivant la température et la précipitation pour les années 2020, 2050 et 2080 sous deux différents scénarios de changement climatique ont été obtenues à partir de la base de données mondiale WorldClim. Un modèle consensus a été dérivé à partir de six modèles individuels en se basant sur trois approches de modélisation différentes (régression, apprentissage automatique et classification). La performance du modèle consensus a été généralement bonne pour les espèces analysées avec une spécificité allant de 81 % à 96,5 % est une sensibilité allant de 88 % à 100 %. L’étude a confirmé l’hypothèse de la large propagation géographique des vecteurs de la leishmaniose sous un scénario de changement climatique. Toutes les espèces étudiées sont supposées gagner des nouvelles zones géographiques actuellement non favorables à la survie du vecteur. / Abstract : Leishmaniasis is endemic in 98 countries with an estimated global prevalence of 12 million cases. The yearly incidence of visceral leishmaniasis is estimated to be 0.3 million cases, while cutaneous leishmaniasis affects approximately one million people each year. In Tunisia, cutaneous leishmaniasis is still a serious health problem with thousands of cases reported every year. Since its first emergence as an epidemic in Kairouan in 1982, the disease has spread in several parts of the country. Control of cutaneous leishmaniasis is mainly based on surveillance of incident cases and treatment and studies carried out to elucidate the spatio-temporal dynamics of the disease are based only on time series analysis of the incidence. The geographical distribution of sandflies during the last decades has shifted northward from latitudes below 45○N in south Europe to latitudes just above 50○N. Recent studies show that some phlebotomine sandflies were recorded in several parts of Germany and Belgium. In Central Europe, some autochthone leishmaniasis cases are being recorded in region traditionally regarded as leishmaniasis-free countries. The objective of this work, in a first step is to i) predict the geographical distribution of cutaneous leishmaniasis and P. papatasi the main vector of the disease in Tunisia. ii) Assess the importance of environmental factors influencing the vector and disease distribution using ecological niche modeling. In a second step we attempted to i) predict the actual distribution of six leishmaniasis vectors in the Mediterranean basin. ii) predict the change in the geographical distribution of these vectors under different climate change scenarios. The MaxEnt model was used to predict the distribution of P. papatasi and cutaneous leishmaniasis cases. 24 environmental variables describing topography, land use and climate were collected. A filed survey was conducted to collect P. papatasi and cutaneous leishmaniasis cases occurrence points. MaxEnt showed a good performance for the prediction of the vector and the disease cases with an area under the ROC curve higher than 0.95 and a maximal Kappa coefficient higher than 0.7. Furthermore, population under the risk of cutaneous leishmaniasis was assessed based on population density by district and vector distribution predicted by the model. To elucidate the effect of climate change on the distribution of various leishmaniasis vectors in the Mediterranean basin, occurrence points of the studied species were collected from scientific literature between 2006 and 2012. Moreover, series of climatic variable describing temperature and precipitation for 2020, 2050 and 2080 under two climate change scenarios were obtained from the global climate database WorldClim. A consensus model was derived from six individual model based on three modeling techniques (statistical regression, machine learning and classification). The consensus model presented a good performance with a specificity varying from 81% to 96.5% and a sensitivity varying from 88% to 100%. The study confirmed the hypothesis of large spread of leishmaniasis vectors under climate change scenario. All studied species are supposed to invade new areas, actually not favourable to the vectors survival.