Thèse doctorat : Génie civil : La Rochelle : 2006. / Bibliogr. p. 169-175.

L'audit thermique des bâtiments existants, par la nature de l'opération, pose un certain nombre de problèmes liés à la réalisation d'outils informatiques. Un des aspects importants, concernant les informations partielles (incertaines, imprécises et lacunaires) et leur traitement dans un programme informatique, n'est jusqu'alors pas pris en compte par des outils ''classiques". On propose donc, dans cette thèse, d'établir les bases d'une nouvelle génération d'outils d'aide au diagnostic, en utilisant les techniques de programmation avancée. Après une analyse détaillée de la problématique, nous avons construit une couche d'un modèle d'acteur à partir de Scheme, et réalisé une maquette informatique. Certains concepts importants, notamment celui de la "continuation" permettant la mise en oeuvre de structures de contrôle souples, puissantes et élégantes, ont été largement exploités. La réalisation de la maquette permet de montrer les intérêts d'une telle approche.

Dus aux conditions uniques en Arctique (des températures et aux rapports de mélange de vapeur d'eau extrêmement bas, à la réflexivité élevée des surfaces de la glace de mer et de la neige, inversion de température dans la basse troposphère et à l'absence de la radiation solaire pendant des périodes prolongées), les processus macro et microphysiques contrôlant la formation des nuages sont complexes et uniques. La validation de ces paramètres atmosphériques simulés par les différents modèles numériques du climat présentement utilisés par les plus grands centres de recherche au monde avec les observations est indispensable pour mieux connaître et, par conséquent, mieux paramétriser ces processus complexes. Le nouveau Modèle Canadien Régional du Climat GEM-LAM (une version à aire limitée du modèle Global Environnemental à Multi-échelle) a été évalué pour la période de septembre 1997 à octobre 1998 au-dessus de l'océan Arctique Ouest. Cette période coincide avec la campagne de mesures du projet SHEBA (Surface Heat Budget of the Arctic Ocean). Les versions 3.2.2 et 3.3.0 du modèle sont évalués dans cette étude. La radiation solaire et terrestre vers le bas à la surface, l'albédo de surface, la vapeur d'eau dans la verticale, les contenus en eau liquide et solide et la couverture nuageuse simulés par GEM-LAM sont évalués avec les données d'observation SHEBA et comparés aux résultats des modèles participants à l'expérience d'inter-comparaison de modèles climatiques régionaux ARCMIP. Une comparaison plus poussée entre les moyennes journalières a été faite et les biais des modèles, le RMSE et le coefficient linéaire de corrélation sont calculés pour plusieurs champs radiatifs et microphysiques. Sur une base mensuelle, les modèles représentent la radiation des longues et courtes longueurs d'onde vers le bas à la surface raisonnablement bien, mais le décalage entre les simulations suivant les différents modèles est plus grand pour les SWD que pour les LWD à la surface. GEM-LAM surestime l'albédo de la surface pendant toutes les saisons dans cette expérience. De plus, le modèle GEM se comporte similairement aux autres modèles participant à cette expérience et tend à sous-estimer la quantité de vapeur d'eau intégrée verticalement pendant l'hiver, tandis que le GEM la surestime durant l'été. La majorité des modèles surestime la couverture nuageuse hivernale, par contre, ils simulent plus ou moins correctement la couverture nuageuse durant l'été. Sur une base quotidienne, tous les modèles participant à l'expérience ARCMlP simulent relativement bien la radiation des courtes longueurs d'onde vers le bas à la surface, mais les valeurs quotidiennes de l'eau liquide intégrée à la verticale sont légèrement moins bien simulées. Les deux versions du modèle GEM-LAM (v3.2.2 et v3.3.0) ont aussi de la difficulté à bien simuler la couverture nuageuse sur une base quotidienne comme les autres modèles participant à cette expérience. Les nuages optiques minces sont l'une des raisons permettant d'expliquer le grand écart entre la couverture des nuages simulée par les modèles et les observations. Pour éliminer ce problème dans notre analyse, un filtre des nuages optiquement minces a été appliqué en sortie du modèle GEM -LAM. Après filtrage des nuages optiquement minces dans le modèle GEM, pour chaque filtre sélectionné (0,5, 1, 1,5, 2), la couverture nuageuse décroit significativement durant la saison hivernale. Par contre, la couverture nuageuse est insensible au filtre durant la saison estivale. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Arctique, Modélisation, GEM-LAM, SHEBA, ARCMIP, Radiation, Couverture nuageuse.

Cette étude a pour objectifs la compréhension et la modélisation de certains effets importants du semis direct avec paillis de résidus (SDP) sur les bilans hydrique, thermique et azoté d'une culture de maïs pluvial. Elle a été menée dans le contexte d'un climat tropical semi-aride au Mexique. D'un point de vue méthodologique, le travail a consisté tout d'abord à identifier les phénomènes physiques et biologiques principaux qu'implique l'itinéraire technique SDP. En s'appuyant sur les résultats de la littérature et sur les acquis au niveau du terrain d'étude, nous avons retenu comme prédominants les effets du SDP sur i) l'évaporation et la température du sol, ii) la structure de surface du sol et l'infiltration, iii) le ruissellement, et iv) le cycle de l'azote et l'activité biologique des sols.<br />Ces effets identifiés, une lourde campagne de collecte de données de terrain a été réalisée de mai à octobre 1998, sur le site expérimental de la Tinaja (état de Jalisco, Mexique). Des mesures complémentaires ont été faites en juin et juillet 1999 sur le même site. Certaines mesures difficilement réalisables sur le terrain ont été effectuées en laboratoire, et sont venues renforcer les données du terrain. A l'issu de ce travail expérimental important, nous disposons d'une base de données relativement riche sur l'ensemble des phénomènes étudiés.<br />L'étude théorique des effets du SDP s'est déroulée en deux phases. La première s'appuie sur une stratégie d'atomisation de la problématique en volets élémentaires spécifiques, traitant les effets du SDP sur i), ii) et iii), en milieu contrôlé. Le premier volet vise à étudier et modéliser les flux d'eau et de chaleur dans le système sol-paillis-atmosphère (transferts turbulents et radiatifs au sein du paillis). Il permet de comprendre et de quantifier l'effet d'un paillis de résidus sur l'évaporation d'un sol et son amplitude thermique, en fonction du taux de couverture. Le deuxième volet consiste à caractériser les effets cumulatifs du SDP sur les propriétés hydrodynamiques d'un sol et donc sa capacité à l'infiltration. Le troisième volet est dédié à l'étude et la modélisation du ruissellement sur sol paillé éventuellement planté, en tenant compte de l'interception de l'eau de pluie, de la rétention de surface, de l'infiltration, de la canalisation du ruissellement et de son écoulement. La deuxième phase relève de la reconstruction par intégration et consiste en un volet unique. Dans ce volet les résultats des trois volets spécifiques précédents sont synthétisés et intégrés dans un modèle numérique de fonctionnement global. Ce dernier prend en compte à la fois les dynamiques de l'eau, de la chaleur et de l'azote. Il reproduit simultanément et de façon couplée tous les effets majeurs du SDP sur le système sol-paillis-plante-atmosphère.

Le compostage est un procédé aérobie de stabilisation de la matière organique (MO). Les transformations de la MO conduisent à des émissions gazeuses (NH3, N2O, CH4, CO2) et une perte d'azote utile pour la fertilisation. La variété des pratiques et les différences de nature des substrats modifient la vitesse de dégradation de la MO, la qualité finale du compost produit et la part d'émissions sous forme de polluants gazeux. L'optimisation du procédé nécessite la prédiction de ces transformations ou la réalisation d'essais empiriques. Cette thèse analyse les interactions entre les processus biologiques, biochimiques, physico-chimiques et thermodynamiques majeurs à l'origine de la stabilisation de la MO et des émissions gazeuses de CO2, H2O, NH3, N2O. Elle se focalise sur le compostage par aération passive en andain d'effluents d'élevage. La méthode s'appuie sur la modélisation dynamique, semi-empirique, du procédé et sur des expérimentations sur des andains de petite taille. Dans un premier temps, l'impact des interactions entre les facteurs clés du compostage, soit la biodégradabilité du carbone, la disponibilité de l'azote, l'humidité et la porosité, sur la cinétique et la stœchiométrie des émissions gazeuses est analysé à partir de données acquises en conditions contrôlées et sur une plateforme de compostage à la Réunion sur des périodes allant de 20 à 80 jours. Les facteurs ont été hiérarchisés : l'effet de la porosité couplé à l'effet de l'humidité joue le rôle majeur dans la régulation des transformations de la MO et des émissions gazeuses. Les effets de la porosité, de l'humidité et du type de structurant sur les transformations de la MO, observés en conditions commerciales, confirment les observations en conditions contrôlées. Alors que les écarts de répétabilité en conditions contrôlées sont faibles sur les températures et bilans de masse, ces écarts deviennent élevés en conditions de reproductibilité sur le terrain. Dans un deuxième temps, un modèle dynamique de compostage a été développé simulant la stabilisation du carbone et de l'azote ainsi que les émissions gazeuses de CO2, H2O, NH3 et N2O. Ce modèle est composé de quatre modules couplés. Le premier simule les cinétiques d'oxydation de la MO par une population microbienne hétérotrophe produisant du CO2 et de la chaleur. Le deuxième décrit les échanges thermiques se déroulant lors du compostage, à l'origine de l'aération passive de l'andain et des cinétiques de température interne et d'évaporation d'H2O. Le troisième module décrit le transfert d'oxygène au sein de l'andain modulé par la porosité et l'humidité. Le quatrième module simule les transformations et la stabilisation de l'azote ainsi que les cinétiques d'émissions d'NH3 et de N2O. Le système global représentant un andain de compostage en aération naturelle est considéré comme homogène en termes de composition et de caractéristiques thermiques. Lors de la phase thermophile, le facteur limitant le plus rapidement l'organisation de la MO est la disponibilité de l'azote. L'abaissement de la porosité entraîne une diminution des pertes gazeuses par l'augmentation de l'organisation de la MO. Le fractionnement initial de la MO ainsi que la teneur initiale en biomasse microbienne sont les facteurs clés pour prédire les cinétiques d'organisation du carbone et de l'azote. Les paramètres spécifiques des émissions d'NH3, d'H2O et de N2O sont initialisés en fonction de la nature du substrat et des caractéristiques physiques de l'andain. La calibration du modèle reste à améliorer, et par la suite le modèle devra être validé en condition industrielles. Les effets modélisés de la porosité, de l'humidité et du type de structurant sur les transformations de la MO permettent d'envisager l'usage du modèle pour optimiser le compostage d'effluents d'élevage.

Cette thèse est une étude expérimentale du comportement thermo-mécanique superélastique d'un fil nanocristallin Ti-50.9Ni at.% Ni en alliage à mémoire de forme (SMA) (diamètre 0.5 mm), après subir un cold work (CW). Les AMF sont capables d'induire des changements de température importants lorsqu'ils sont chargés mécaniquement. Ce phénomène est dû à un important couplage thermomécanique présent dans cette transformation de phase solide entre les phases Austénite (A) et Martensite (M).La chaleur latente par unité de masse (ΔH) tout au long de la transformation de phase est l'énergie responsable de cette variation de température. La détermination de ΔH est généralement effectuée par calorimétrie à balayage différentiel (DSC). Cependant, pour les SMA nanocristallins, les résultats DSC obtenus ne sont pas concluants sur la détermination de cette propriété.Dans ce travail, une méthode utilisant la corrélation d'image numérique (DIC) et les mesures de champ thermique (TFM) a été utilisée pour analyser les couplages thermomécaniques lors d'une transformation de phase induite par contrainte. Des champs cinématiques et thermiques ont été acquis lors d'essais de traction superélastiques réalisés sur des fils CW NiTi soumis à différentes températures de traitements thermiques (TTT) allant de 523 à 598 K pendant 30 min. Un tel traitement thermique à basse température favorise une boucle totalement superélastique sans plateau de contrainte et sans déformation de type Lüders. En supposant un modèle thermique uniforme, les sources de chaleur impliquées lors du chargement cyclique ont été estimées. Cette puissance thermique par unité de masse a été comparée à la puissance mécanique et intégrée au fil du temps pour obtenir l'équilibre énergétique. De plus, grâce à une analyse thermodynamique basée sur l'énergie libre de Gibbs, les valeurs de ΔH, ainsi que la fraction de martensite, ont été estimées au cours des transformations de phase A-M directe et inverse M-A. L'analyse des résultats a conduit aux conclusions suivantes: (1) Les puissances et énergies thermiques et mécaniques présentaient une dépendance significative vis-à-vis du TTT. (2) Malgré l'effet important des valeurs du TTT sur les réponses mécaniques et thermiques, les ΔH obtenues étaient très proches pour tous les TTT et dans la même gamme de valeurs fondée dans la littérature pour un alliage Ti-50.9Ni at.% Ni entièrement recuit testé par technique DSC. (3) Pour une deformation donnée, la fraction de martensite augmente avec l'augmentation de TTT. (4) Pour une contrainte imposée de 4,5%, la fraction de martensite augmente de 30% à 40% en augmentant le TTT de 523K à 598K. / This PhD thesis is an experimental study of the thermomechanical superelastic behaviour of a Ti-50.9Ni at.% Ni Shape Memory Alloy (SMA) nanocrystalline thin wire (diameter 0.5 mm), in a Cold Worked (CW) state. SMAs are capable of inducing important temperature change when they are mechanically loaded. This phenomenon is due to an important thermomechanical coupling present in this solid phase transformation between Austenite (A) and Martensite (M) phases. The latent heat per unit of mass (∆H) throughout the phase transformation is the energy responsible of this temperature variation. The determination of ∆H is generally performed by differential scanning calorimetry (DSC). However, for nanocrystalline SMAs, the obtained DSC results are non conclusive on the determination of this property.In this work, a method using digital image correlation (DIC) and thermal field measurements (TFM) was used to analyse the thermomechanical couplings during a stress induced phase transformation (SIPT). Kinematics and thermal full fields were acquired during superelastic tensile tests performed on the CW NiTi wire submitted to different heat treatments temperatures (HTT) ranging from 523 to 598 K during 30 min. Such a heat treatment at low temperature promoted a fully superelastic loop without stress plateau and no Lüders-like deformation. Assuming a uniform thermal model, the heat sources involved during the cyclic loading were estimated. This thermal power per unit of mass was compared to the mechanical one and integrated over the time to get energy balance. Further, through a thermodynamic analysis based on the Gibbs free energy, the values of ∆H, as well as the martensite fraction, were estimated during the forward A-M and reverse M-A phase transformations. The analysis of the results led to the following conclusions: (1) Thermal and mechanical powers and energies presented a significant dependence on the HTT. (2) Despite the strong effect of the values of the HTT on mechanical and thermal responses, the obtained ∆H were very close for all HTT and in the same range of values founded in the literature for a fully annealed Ti-50.9Ni at.% Ni alloy tested via DSC technique. (3) For a given strain, martensite fraction increases with increasing HTT. (4) For an imposed strain of 4.5%, the martensite fraction increases from 30% to 40% when increasing HTT from 523K to 598K.

La variété des pratiques de compostage et les différences de nature des substrats modifient la vitesse de dégradation de la matière organique (MO), la qualité finale du compost produit et la part d'émissions sous forme de polluants gazeux. L'optimisation du procédé nécessite la prédiction de ces transformations ou la réalisation d'essais empiriques. Cette thèse analyse les interactions entre les processus biologiques, biochimiques, physico-chimiques et thermodynamiques majeurs à l'origine de la stabilisation de la MO et des émissions gazeuses de CO2, H2O, NH3, N2O. Elle se focalise sur le compostage par aération passive en andain d'effluents d'élevage. La méthode s'appuie sur la modélisation dynamique et des expérimentations. L'impact des interactions entre la biodégradabilité du carbone, la disponibilité de l'azote, l'humidité et la porosité, sur la cinétique et la stœchiométrie des émissions gazeuses est analysé à partir d’observations en conditions contrôlées et sur une plateforme de compostage. Les écarts de répétabilité faibles en conditions contrôlées deviennent élevés en conditions de reproductibilité sur le terrain. Un modèle dynamique de compostage a été développé simulant la stabilisation de la MO et les émissions de CO2, H2O, NH3, N2O. Il est composé de quatre modules simulant les cinétiques d'oxydation de la MO par une population microbienne, les échanges thermiques et l'aération passive, le transfert d'oxygène, les transformations de l'azote. Lors de la phase thermophile, le facteur limitant le plus rapidement l'organisation de la MO est la disponibilité de l'azote. L'abaissement de la porosité entraîne une diminution des pertes gazeuses par l'augmentation de l'organisation de la MO. Le fractionnement initial de la MO ainsi que la teneur initiale en biomasse microbienne sont les facteurs clés pour prédire les cinétiques d'organisation de la MO. La calibration des paramètres d’initialisation du modèle reste à améliorer / The variety of the practices and the differences in nature of the substrates modify the kinetics of degradation of organic matter (OM), the final quality of the produced compost and the fraction emitted as gaseous pollutants. To optimize the composting process, it is required to predict these transformations or to do some empirical test. This thesis analyzes the interactions between the main biological, biochemical, physicochemical and thermodynamic processes which explain the OM stabilization and the gaseous emissions of CO2, H2O, NH3, N2O. Focus is done on windrow composting with passive aeration. The method is based on dynamic semi-empirical modeling of the process and experimentations. The impact of the interactions between bioavailability of carbon and nitrogen, moisture and porosity, on the kinetics and the stoechiometry of the gaseous emissions is analyzed with data gathered in controlled conditions and on a commercial composting plant. The repeatability differences are small in controlled conditions but higher in field conditions. A dynamic model of composting was developed simulating the stabilization of OM as well as the emissions of CO2, H2O, NH3, N2O. This model is composed of four coupled modules which simulate the kinetics of oxidation of OM by a heterotrophic microbial population, the heat transfers leading to the passive aeration, the transfer of oxygen, the transformations of nitrogen. During the thermopilic phase, the first factor limiting the organization of OM is the availability of nitrogen. The decrease in porosity induces a reduction in the gas losses through the increase in the organization of OM. The initial OM fractionation and the initial microbial biomass are the key factors to predict the kinetics of organization of OM. The calibration of the specific parameters used for model initialization needs a further implementation

Le joint à rainures hélicoïdales (JRH) est l'une des solutions techniques d'étanchéité sans contact utilisées dans les machines tournantes. Ce dispositif est conçu pour des applications bien particulières qui nécessitent une durée de vie et une non-tolérance aux fuites au-delà des limites que peuvent satisfaire les joints à contact. Le JRH est caractérisé par l'absence d'usure due au jeu radiale nettement supérieure aux amplitudes des aspérités et les défauts de fabrication. L'étanchéité est obtenue grâce aux rainures hélicoïdales présentes sur l'une des ses surfaces internes. Ces rainures sont à l'origine de phénomènes hydrodynamiques synthétisant un débit de pompage de même ordre que le débit de fuite.Dans ce travail, un modèle numérique pour le calcul d'étanchéité dans les JRH est proposé. Basé sur la théorie des films minces, le comportement de ce dernier est déterminé par le calcul du champ de pression et du remplissage qui satisfont l'Equation de Reynolds Modifiée (ERM). Cette dernière permet de bien gérer les frontières de rupture et de reformation du film. La résolution est faite par la méthode des éléments finis.La caractérisation du pouvoir d'étanchéité du JRH est faite par « la longueur utile » qui spécifie largeur, dans la direction axiale, de la zone occupée par le fluide lorsque l'étanchéité s'établisse. Cette étendue du domaine étant une inconnue du problème, on itère sur la longueur du joint jusqu'à l'obtention d'un débit axial nul sur le bord. Néanmoins, l'étanchéité dans le JRH dépend d'un certain nombre de paramètres géométriques et de fonctionnement. Il s'agit de la forme des rainures, leurs inclinaisons ainsi que la vitesse de rotation. Tout d'abord, une géométrie optimale de fonctionnement a été déterminée. Les aspects turbulents de l'écoulement et le comportement thermique, selon un bilan thermique global, sont également étudiés. Enfin, l'introduction des effets d'excentricité a permis de statuer sur les phénomènes dynamiques dans le JRH. / The viscoseal is one of the contacts less sealing technical solutions used in machinery. This device is designed for very specific applications requiring durability and non-tolerance leakage past the limits that can satisfy joints contact. The viscoseal is characterized by the friction absence due to radial clearance well above the asperities amplitudes and the manufacturing defects. The sealing is obtained by the helical grooves formed in one of its internal surfaces. These grooves induce a hydrodynamic phenomenon that synthesizes same pumping rate of the same order as the leakage rate.In this work, a numerical model is proposed to calculate the seal in the viscoseal. Based on the thin film theory, the behavior of the latter is determined by the calculation of the pressure field and the filling that satisfy the Modified Reynolds Equation (MRE). The latter allows managing the borders of the film breaking and reformation. Resolution is made by the finite element method.The sealing power characterization in the viscoseal is made by the "sealing length" that specifies width of the fluid full area, in the axial direction, when the sealing is established. This domain extension is unknown, it iterates over the length of the seal until a zero axial flow over the edge. However, in the sealing depends on several geometrical and operating parameters. It is about the shape of the grooves, their angle orientation of and the journal speed.First, an optimal operating geometry was determined. Turbulent aspects of flow and thermal behavior, according to a global heat balance, are also studied. Finally, the introduction of eccentricity effects allowed approving dynamic phenomena in the viscoseal.

Dans le présent projet, on cherche à évaluer le potentiel de récupération et d'intégration des déchets thermiques dans une usine de production primaire d'aluminium. Il s'agit en premier lieu de «cartographier» la chaleur et les températures en jeu; c'est-à-dire construire un inventaire détaillé des rejets thermiques, mais aussi de la consommation de chaleur dans divers bâtiments ou procédés. L'ensemble de l'usine est inclus dans l'analyse: cuves d'electrolyse, secteur de production des anodes, centre de coulée, entrepôts, bureaux, etc. En second lieu, il s'agit d'examiner les procédés qui dissipent de la chaleur pour déterminer les limites therrmodynamiques qui s'appliquent à la récupération thermique. L'analyse exergétique est employée afin d'établir pleinement les possibilités d'intégration thermoénergétiques. L'aluminerie considérée est celle d'Alcoa Deschambault Québec (ADQ). Les résultats montrent que le potentiel d'intégration réelle de la chaleur identifié en usine est faible dû i) aux équipements et technologies de récupération et intégration thermique accessibles et ii) à une demande limitée pour la chaleur de basse/moyenne température. Il n'en reste pas moins qu'il est actuellement possible thermodynamiquement de satisfaire les besoins de chauffage ambiant à l'usine d'ADQ avec les rejets thermiques contenus dans les gaz d'échappements des cuves. Aussi, -10 MW de puissance pourrait être généré en continu avec la chaleur identifiés dans les secteurs d'électrolyse, de cuisson des anodes, et de coulée, compte tenu de l'efficacité des machines thermiques disponibles sur le marché.

Surface manifestations of hydrothermal fluids such as fumaroles and hot springs provide valuable information about the level of activity of a volcano during quiescent period. Geochemical study of gas and spring waters is useful to elaborate geochemical model for magmatic-hydrothermal system. Furthermore, temporal geochemical monitoring of these fluids with time provides a better understanding in processes occurring inside the volcano and can be useful to detect any changes in the activity of the magmatic-hydrothermal system. This thesis investigates two hydrothermal systems at Kelud and Papandayan volcanoes that are located at Java Island in Indonesia. Kelud is considered as one of the most dangerous volcanoes of Java because of its frequent eruptions. After the last eruption that occurred in 1990, a new lake rapidly filled the crater of Kelud volcano. Water samples collected since 1993 are near neutral Na-K chloride fluids and are typical of aged hydrothermal system where the acidity has been completely neutralized by fluid-rock interaction and where the emission of acid magmatic gases has stopped. Two sudden increases in lake temperature in 1996 and 2001 were accompanied by rapid changes in lake water compositions and suggest the existence of two hydrothermal systems feeding the lake: a shallow hydrothermal system dominated by Ca-Mg sulfate waters and a deepest aquifer with neutral alkali chloride waters. From 2001 to 2005, measurements of CO2 emitted by the surface of the lake were performed by using the accumulation chamber method modified in order to work at the surface of a crater lake. Two statistical methods were used to process data: the graphical statistical and stochastic simulation methods. The results of graphical statistical approach showed that two different degassing processes are acting at the lake surface: one corresponding to CO2 fluxes resulting from rising bubbles and the other corresponding to equilibrium diffusion of dissolved CO2 at the water-air surface. Total CO2 emission rate estimated by stochastic simulation ranges from 105 t/day for 2001 to 32 t/day for 2005. Thermal energy released by the lake was also estimated by using an energy balance model with a new constraint using the CO2 flux. The thermal flux decreased from 200 MW (2001) to 100 MW (2002) and then remained stable. Correlation between the chemical data of waters, the fluxes of CO2 and energy show that a constant decrease in the level of activity of the volcano since 1993 occurred although the lake temperature has been stable since 2003. Since the last magmatic eruption that occurred in 1772, phreatic eruptions occur on Papandayan volcano with the last one in 2002. The volcanic material ejected during this eruption is essentially made of altered rocks from within the hydrothermal system. The interaction of acid waters with the host rocks corresponds to an advanced argilic alteration. The chemical compositions of waters from Papandayan volcano and Kelud lake waters are contrasting. Indeed, the spring waters sampled since 1994 are acid sulfate-chloride waters and acid sulfate waters. The chemical and isotopic analyses of gases and waters suggest a significant magmatic contribution in SO2, HCl and HF to the hydrothermal system. The chemical composition of waters sampled after the 2002 eruption have provided information about origin of this eruption. Decrease in chloride concentration and in delta 34S of dissolved sulfates showed that the magmatic contribution in these fluids are less important and that the waters are likely to be formed by the condensation of steam (H2O, H2S) rising from a boiling aquifer.