Les calottes polaires, ou inlandsis, sont parmi les principaux contributeurs à la montée des océans. Ces systèmes dynamiques gagnent de la masse par accumulation de neige, et en perdent par fonte au contact de l'océan et à la surface, ainsi que par le vêlage d'icebergs. Depuis plus de trois décennies, les observations ont montré que les calottes polaires de l'Antarctique et du Groenland perdent plus de masse qu'ils n'en gagnent. L'évolution des glaciers suite à ce déséquilibre de masse est devenue aujourd'hui l'une des problématiques les plus importantes des implications du changement climatique. Le Groupe d'experts intergouvernemental sur l'évolution du climat (GIEC) a identifié la contribution des glaciers comme l'un des facteurs clés d'incertitude de prédiction de l'élévation du niveau des mers. La modélisation numérique est le seul outil efficace pour répondre à cette question. Cependant, modéliser l'écoulement de glace à l'échelle du Groenland ou de l'Antarctique représente un défi à la fois scientifique et technique. Deux aspects clés de l'amélioration de la modélisation des glaciers sont abordés dans cette thèse. Le premier consiste à déterminer certaines propriétés non mesurables de la glace par méthode inverse. La friction ou la rigidité des barrières de glace, sont des paramètres qui ne peuvent être mesurés directement et doivent donc être déduits à partir d'observations par télédétection. Nous appliquons ici ces inversions pour trois modèles d'écoulement de glace de complexité croissante: le modèle bidimensionnel de MacAyeal/Morland, le modèle dit d'ordre supérieur de Blatter/Pattyn et le modèle full-Stokes. Les propriétés ainsi calculées sont ensuite utilisées pour initialiser des modèles grande-échelle et pour déterminer le degré de complexité minimum nécessaire pour reproduire correctement la dynamique des glaciers. Le second aspect abordé dans ce travail est l'amélioration de la consistance des données pour la modélisation numérique. Les données disponibles sont souvent issues de campagnes de mesures s'étalant sur plusieurs années et dont résolutions spatiales varient, ce qui rend leur utilisation pour des simulations numériques difficiles. Nous présentons ici un algorithme basé sur la conservation de la masse et les méthodes inverses pour construire des épaisseurs de glace qui sont consistantes avec les mesures de vitesse. Cette approche empêche la redistribution artificielle de masse qu'engendrent généralement les autres méthodes de cartographie de l'épaisseur de glace, ce qui améliore considérablement l'initialisation des modèles d'écoulement de glace. Les avancées présentées ici sont des étapes importantes afin de mieux caractériser de manière précise les glaciers et de modéliser leur évolution de manière réaliste.

Le réchauffement climatique actuel a une conséquence directe sur la perte de masse des calottes polaires. Reproduire les mécanismes responsables de cette perte de masse et prévoir la contribution des calottes à l’élévation du niveau des océans d’ici la fin du siècle est dès lors l’un des défis majeurs de la modélisation de l’écoulement des calottes polaires. Les modèles d’écoulement permettent de réaliser de telles prévisions mais ces simulations, à court terme, sont très sensibles à leur état initial habituellement construit à partir d’observations de terrain. Malheureusement, certains paramètres comme le frottement entre la glace et le socle rocheux ainsi que la topographie basale sont souvent méconnus à cause du manque d’observations directes ou des larges incertitudes liées à ces observations. Améliorer la connaissance de ces deux paramètres à la fois pour le Groenland et l’Antarctique est donc un pré-requis pour réaliser des projections fiables. Les méthodes d’assimilation de données et les méthodes inverses permettent alors de surmonter ce problème.Cette thèse présente deux algorithmes d’assimilation de données permettant de mieux contraindre simultanément le frottement basal et la topographie basale à partir d’observations de surface. L’un des algorithmes est entièrement basé sur la méthode adjointe tandis que le second se base sur une méthode cyclique couplant l’inversion du frottement basal avec la méthode adjointe et l’inversion de la géométrie basale à l’aide de la relaxation newtonienne. Les deux algorithmes ont été implémentés dans le modèle d’écoulement de glace éléments finis Elmer/Ice et testés dans une expérience jumelle qui montre une nette amélioration de la connaissance des deux paramètres recherchés. L’application des deux algorithmes à la région de la Terre de Wilkes réduit l’incertitude liée aux conditions basales en permettant, par exemple, d’obtenir plus de détails sur la géométrie basale en comparaison avec les modèles numériques de terrain habituels. De plus la reconstruction simultanée du frottement et de la géométrie basale permet de réduire significativement les anomalies de divergence de flux habituellement obtenues lors de l’inversion du frottement seul. Nous étudions finalement l’impact des conditions basales ainsi inversées sur des simulations pronostiques afin de comparer la capacité des deux algorithmes à mieux contraindre la contribution future des calottes polaires à l’augmentation du niveau des océans. / The current global warming has direct consequences on ice-sheet mass loss. Reproducing the responsible mechanisms and forecasting the potential ice-sheets contribution to 21st century sea level rise is one of the major challenges in ice-sheet and ice flow modelling. Ice flow models are now routinely used to forecast the potential ice-sheets contribution to sea level rise. Such short term simulations are very sensitive to model initial state, usually build from field observations. However, some parameters, such as the basal friction between icesheet and bedrock as well as the basal topography, are still badly known because of a lake of direct observations or large uncertainty on measurements. Improving the knowledge of these two parameters for Greenland and Antarctica is therefore a prerequisite for making reliable projections. Data assimilation and inverse methods have been developed in order to overcome this problem. This thesis presents two different assimilation algorithms to better constrain simulaneouslybasal friction and bedrock elevation parameters using surface observations. The first algorithm is entierly based on adjoint method while the second algorithm uses a cycling method coupling inversion of basal friction with adjoint method and inversion of bedrock topography with nudging method. Both algorithms have been implemented in the finite element ice sheet and ice flow model Elmer/Ice and tested in a twin experiment showing a clear improvement of both parameters knowledge. The application of both algorithms to regions such as the Wilkes Land in Antartica reduces the uncertainty on basal conditions, for instance providing more details to the bedrock geometry when compared to usual DEM. Moreover,the reconstruction of both bedrock elevation and basal friction significantly decreases ice flux divergence anomalies when compared to classical methods where only friction is inversed. We finaly sudy the impact of such inversion on pronostic simulation in order to compare the efficiency of the two algorithms to better constrain future ice-sheet contribution to sea level rise.

Les présents travaux s'inscrivent dans la continuité des projets soutenus depuis plusieurs années par le Campus International pour la Sécurité et l'Intermodalité dans les Transports (CISIT CPER). Une approche objective est développée dans ce mémoire afin d'identifier les causes amenant à la perte d'efficacité et à la dégradation du tribosystème de freinage. Le comportement d'une garniture de composition simplifiée, une matrice de résine phénolique renforcée de particules sphériques d'acier, est étudié dans des conditions de pression, vitesse et température de freinage automobile. Le développement de la plateforme High Speed Tribology permet la caractérisation tribologique à échelle réduite pour des vitesses de glissement de plusieurs dizaines de mètres par seconde. L'acquisition de température dans le disque est possible grâce à l'instrumentation d'un dispositif de mesure par télémétrie. Le champ de température surfacique et les flux thermiques dissipés sont identifiés à partir de techniques inverses de conduction de la chaleur. La dépendance du coefficient de frottement macroscopique vis à vis de la température de contact, de la pression moyenne et de la vitesse de glissement est évaluée pour un pion en résine pure et pour le matériau bicomposant. Les observations des surfaces témoignent de l'influence de ces paramètres sur le mode de dégradation de la zone en contact. La décohésion d'un renfort de la matrice est étudiée par le biais d'un modèle numérique à échelle mésoscopique. Les propriétés mécaniques de la résine et de l'interface sont déterminées par des essais rhéologiques dans une étuve pour différentes températures. Les simulations montrent que les hautes pressions et les hautes températures conduisent à la rupture de l'interface matrice - renfort. Des corrélations entre les essais menés sur le HST et le modèle numérique sont présentées.

La détermination des paramètres inertiels des différents segments constituants le corps humain est un problème récurrent et très actuel de la biomécanique. Les méthodes classiques déterminent ces paramètres par des équations de régression (modèles de Zatsiorsky, Dempster...) ou en modélisant les segments par des formes géométriques simples (modèles d'Hanavan, de Hatze...). Les techniques d'imagerie permettent également d'obtenir ces paramètres mais ces méthodes sont peu accessibles et très onéreuses. Des problèmes similaires d'estimation des paramètres inertiels se posent pour les robots manipulateurs. Aux méthodes classiques de C.A.O. ou de mesure, les roboticiens ont proposé des méthodes d'identification. Celles-ci consistent à calculer les paramètres inertiels à partir d'un système sur-déterminé dans lequel les paramètres inertiels apparaissent linéairement. La transposition de la méthode d'identification au système polyarticulé que constitue le corps humain nécessite de développer un « modèle de réaction externe » dans lequel la cinématique des segments et les forces externes au système interviennent. En effet, alors que pour les robots manipulateurs les efforts articulaires (entre les segments) sont connus, ils ne peuvent pas être connus directement pour le corps humain et seuls les efforts externes peuvent être mesurés directement (par une plate-forme de forces par exemple). L'originalité de la thèse est de présenter une formulation matricielle des équations du mouvement dans laquelle les efforts externes sont regroupés dans une matrice antisymétrique et les paramètres inertiels sont contenus dans une matrice (3n+1×3n+1) symétrique et définie positive (n est le nombre de segments). Un algorithme de gradient conjugué projeté est ensuite développé afin de déterminer cette matrice. L'algorithme converge vers une matrice symétrique et définie positive à partir de données provenant de la cinématique des segments et des efforts externes.<br />Plusieurs applications de cette méthode ont été réalisées et publiées. Les premières concernent la validation de la méthode d'identification à partir de mouvements simulés obtenus à l'aide du logiciel HuMAnS (INRIA). Cette application a permis de mettre en évidence l'importance de connaître précisément les centres articulaires afin de déterminer les paramètres inertiels. Les autres applications ont été menées in vivo et les paramètres inertiels des segments du membre supérieur ont ainsi été déterminés et comparés à des modélisations classiques. La méthode d'identification permet d'obtenir les paramètres inertiels personnalisés et constitue ainsi une alternative très intéressante aux modélisations classiques particulièrement dans le cas d'études incluant des populations particulières (sportifs « hors normes », obèses, handicapés) pour lesquelles aucun modèle n'est adapté.

L’objectif de cette thèse était de développer et d’appliquer une nouvelle méthode de mesure permettant de dresser des cartographies instantanées des températures et des concentrations des diverses espèces présentes dans les gaz de combustion et aider à la détermination des performances de moteurs et mieux comprendre les mécanismes de réaction.Pour cela, un moyen de mesure optique, non intrusif pour ne pas interagir sur le fonctionnement du moteur, a été mis au point, ce qui a permis ainsi une meilleure analyse des performances et donc une mise au point plus rapide des moteurs. De par ses caractéristiques, le système peut réaliser des mesures résolues temporellement et spatialement. Outre les investigations dans les écoulements réactifs des propulseurs à réaction, ce système de mesure peut être également utilisé ultérieurement pour la mise au point des moteurs cycliques à détonation, à allumage commandé (essence) ou par compression (diesel, HCCI). Les temps de réponse (inférieures à 10 μs) du moyen de mesure sont compatibles avec les phénomènes étudiés. Cette technique de mesure peut également être embarquable sur un aéronef pour effectuer des mesures en vol.Le but était donc de mettre au point un système optique permettant la mesure de ces différents paramètres à distance. La spectroscopie est un outil adapté à cette étude. C’est une méthode physique non intrusive qui, associée à un LIDAR (télédétection laser), permet la mesure précise et instantanée de la température et de la concentration des gaz de combustion et de connaître leur évolution temporelle. Il peut être envisagé par la suite de mesurer directement l’avancement de la combustion afin d’améliorer le rendement et par voie de conséquence diminuer les rejets de polluants. / The aim of the present PhD thesis was to develop and apply a new measurement technique allowing one to create maps of the temperature and t-he concentration of various species within combustion gas, to determine the performances of the engines and get a better understanding of the reaction mechanisms.In order to do so, a non-invasive optical measurement technique was setup. This has enabled us to have better an analysis of the performances and hence, a faster adjustment of the engines. Thanks to its characteristics, the system is able to do some spatial and temporal resolved measurements. Apart from reacting flows in propellers, this measurement technique can also be used for the design of detonation and spark ignition engines or compression engines. The response times (lower than 10 μs) are compatible with the phenomena under investigation. This technique can also be used in airplanes to enable on-flight measurements.The aim was hence to setup an optical system allowing the measurement of these different parameters from a distance. Spectroscopy is a well-suited tool for this work. It is a non-intrusive physical measurement which, combined with LIDAR, allows the precise and instantaneous measurement of the temperature and the concentration of combustion gas as well as the following of their temporal evolution. It would be worth considering measuring the progress of the combustion directly so as to improve the performance and hence reduce the formation of pollutants.

Les calottes polaires, ou inlandsis, sont parmi les principaux contributeurs à la montée des océans. Ces systèmes dynamiques gagnent de la masse par accumulation de neige, et en perdent par fonte au contact de l’océan et à la surface, ainsi que par le vêlage d’icebergs. Depuis plus de trois décennies, les observations ont montré que les calottes polaires de l’Antarctique et du Groenland perdent plus de masse qu’ils n’en gagnent. L’évolution des glaciers suite à ce déséquilibre de masse est devenue aujourd’hui l’une des problématiques les plus importantes des implications du changement climatique. Le Groupe d’experts intergouvernemental sur l’évolution du climat (GIEC) a identifié la contribution des glaciers comme l’un des facteurs clés d’incertitude de prédiction de l’élévation du niveau des mers. La modélisation numérique est le seul outil efficace pour répondre à cette question. Cependant, modéliser l’écoulement de glace à l’échelle du Groenland ou de l’Antarctique représente un défi à la fois scientifique et technique. Deux aspects clés de l’amélioration de la modélisation des glaciers sont abordés dans cette thèse. Le premier consiste à déterminer certaines propriétés non mesurables de la glace par méthode inverse. La friction ou la rigidité des barrières de glace, sont des paramètres qui ne peuvent être mesurés directement et doivent donc être déduits à partir d’observations par télédétection. Nous appliquons ici ces inversions pour trois modèles d’écoulement de glace de complexité croissante: le modèle bidimensionnel de MacAyeal/Morland, le modèle dit d’ordre supérieur de Blatter/Pattyn et le modèle full-Stokes. Les propriétés ainsi calculées sont ensuite utilisées pour initialiser des modèles grande-échelle et pour déterminer le degré de complexité minimum nécessaire pour reproduire correctement la dynamique des glaciers. Le second aspect abordé dans ce travail est l’amélioration de la consistance des données pour la modélisation numérique. Les données disponibles sont souvent issues de campagnes de mesures s’étalant sur plusieurs années et dont résolutions spatiales varient, ce qui rend leur utilisation pour des simulations numériques difficiles. Nous présentons ici un algorithme basé sur la conservation de la masse et les méthodes inverses pour construire des épaisseurs de glace qui sont consistantes avec les mesures de vitesse. Cette approche empêche la redistribution artificielle de masse qu’engendrent généralement les autres méthodes de cartographie de l’épaisseur de glace, ce qui améliore considérablement l’initialisation des modèles d’écoulement de glace. Les avancées présentées ici sont des étapes importantes afin de mieux caractériser de manière précise les glaciers et de modéliser leur évolution de manière réaliste. / Ice sheets are amongst the main contributors to sea level rise. They are dynamic systems; they gain mass by snow accumulation, and lose it by melting at the ice-ocean interface, surface melting and iceberg calving at the margins. Observations over the last three decades have shown that the Greenland and Antarctic ice sheets have been losing more mass than they gain. How the ice sheets respond to this negative mass imbalance has become today one of the most urgent questions in understanding the implications of global climate change. The Intergovernmental Panel on Climate Change (IPCC) has indeed identified the contribution of the ice sheets as a key uncertainty in sea level rise projections. Numerical modeling is the only effective way of addressing this problem. Yet, modeling ice sheet flow at the scale of Greenland and Antarctica remains scientifically and technically very challenging. This thesis focuses on two major aspects of improving ice sheet numerical models. The first consists of determining non-observable ice properties using inverse methods. Some parameters, such as basal friction or ice shelf hardness, are difficult to measure and must be inferred from remote sensing observations. Inversions are developed here for three ice flow models of increasing complexity: MacAyeal/Morland’s shelfy-stream model, Blatter/Pattyn’s higher order model and the full-Stokes model. The inferred parameters are then used to initialize large-scale ice sheet models and to determine the minimum level of complexity required to capture ice dynamics correctly. The second aspect addressed in this work is the improvement of dataset consistency for ice sheet modeling. Available datasets are often collected at different epochs and at varying spatial resolutions, making them not readily usable for numerical simulations. We devise here an algorithm based on the conservation of mass principle and inverse methods to construct ice thicknesses that are consistent with velocity measurements. This approach therefore avoids the artificial mass redistributions that occur in existing algorithms for mapping ice thickness, hence considerably improving ice sheet model initialization. The advances made here are important steps towards the ultimate objective of accurate characterization of ice sheets and the realistic modeling of their evolution.

La validation du modèle de réponse globale des échantillons ou des structures présente des lacunes qui peuvent être surmontés par l'utilisation de techniques de mesure de plein champ. Les techniques de mesure plein champ offrent la possibilité d'acquérir une quantité importante de données expérimentales qui pourraient être utiles dans le cadre de l'identification des paramètres de la Loi de comportement. Parmi les techniques de mesure de champ complet, le plus populaire est la corrélation d'images mono et stéréo. Les stratégies existantes pour rendre l'utilisation des données de plein champ comme la méthode de champ virtuel, la méthode écart d'équilibre, la méthode de Gap équation Constitutive et la méthode de réciprocité ont été limitées à des applications 2D. Cependant, pour un échantillon d'épaisseur finie, les effets 3D doivent être remplis. Plus important encore, dans le cas de la plasticité, triaxialité du stress joue un rôle important. Son effet doit donc être inscrits dans la modélisation du comportement de la matière. Ainsi, dans cette thèse, nous proposons une méthode pour identifier les paramètres d'une loi de comportement élasto-plastique dans lequel le modèle mécanique peut comporter cinématique 3D. La stratégie s'est avérée être le bruit robuste, presque indépendamment de la conjecture de paramètre initial et de raffinement de maillage et permet de différencier les modèles de comportement avec la même réponse globale sur la base de l'erreur de corrélation moyenne. L'identification est avérée bon pour une ou plusieurs caméras. La validation de la stratégie est faite d'acier inoxydable. La réponse globale du déplacement de charge identifié de la stratégie avec la corrélation d'image mécanique mono et stéréo est très proche des expériences. Enfin, les paramètres du matériau ont été identifiés avec l’estimation initial très différents mais tous convergent vers les mêmes valeurs finales qui montrent la solidité de la stratégie proposée. / The objective of the thesis is to formulate a strategy that gives a robust identification of constitutive law from full-field measurements taking into account 3D effects. Model validation from global response of samples or structures has shortcomings that can be overcome by the use of full-field measurement techniques. Full-field measurement techniques offer the opportunity to acquire large amount of experimental data that might be useful in the context of identification of constitutive law parameters. Among the full field measurement techniques the most popular is digital and stereo image correlation. The existing strategies to make use of full-field data like the Virtual Field Method, the Equilibrium Gap Method, the Constitutive Equation Gap Method and the Reciprocity Method were limited to 2D applications. However, for a specimen with finite thickness 3D effects must be included. Most importantly, for the case of plasticity, stress triaxiality plays an important role. Its effect must therefore be accounted for in the modelling of the constitutive behaviour of the material. Thus in this thesis we propose a method to identify the parameters of an elasto-plastic constitutive law in which the mechanical model can have 3D kinematics. The strategy has been shown to be noise robust, almost independent of initial parameter guess and mesh refinement and allows differentiating between constitutive models with same global response on the basis of mean correlation error. The identification is shown to be good for both single and multiple cameras. The strategy validation is done for stainless steel. The global identified load displacement response of the strategy for mono and stereo mechanical image correlation is very close to the experiments. Lastly, the material parameters have been identified with very different initial guess but all converge to the same final values which show the robustness of the proposed strategy.

L'évaluation et le contrôle non destructifs (E.C.N.D.) des matériaux et des structures sont une problématique industrielle très importante dans les domaines du transport, de l'aéronautique et du spatial, et dans le milieu médical. La thermographie infrarouge active est une technique d'E.C.N.D qui consiste à apporter une excitation extérieure afin d'entraîner une élévation de température dans le matériau, puis à évaluer le champ de température résultant à la surface. Cependant, les excitateurs thermiques utilisés (lampes flash, halogènes, lasers) agissent uniquement sur la surface du matériau. Plusieurs systèmes de conversion d'énergie peuvent en revanche mener à l'apparition de sources volumiques : on peut citer en particulier les phénomènes de thermo-acoustique, de thermo-induction, de thermomécanique ou de thermochimie. Par exemple, une excitation par ondes ultrasonores peut entraîner des sources thermiques volumiques si le matériau est viscoélastique ou s'il y a présence de défaut. La reconstruction de ces sources est donc la première étape permettant de remonter aux paramètres responsables de l'échauffement. Caractériser une source thermique consiste à reconstruire sa géométrie et la puissance qu'elle génère. Cependant, l'identification de sources thermiques volumiques par la mesure des champs de température de surface est un problème mathématiquement mal posé. Le caractère diffusif de la température en est le principal responsable. Dans ce travail, la reconstruction 3D des sources volumiques à partir du champ de température résultant à la surface, mesuré par InfraRouge, est étudié. Tout d'abord, une analyse du problème physique permet de spécifier les limites de la reconstruction. En particulier, un critère sur la résolution spatiale atteignable est défini et une limitation de reconstruction pour les sources en profondeur est mise en lumière. Ensuite, une méthode de reconstruction par approche probabiliste est proposée et comparée aux méthodes d'inversions existantes. Le temps d'exécution et la sensibilité au bruit de mesure sont étudiés pour chacune de ces méthodes. Des applications numériques et expérimentales seront enfin présentées pour illustrer les résultats. / Non Destructive Testing (N.D.T.) of materials and structures is a very important industrial issue in the fields of transport, aeronautics and space and in the medical domain. Active infrared thermography is a N.D.T. method that consists in providing an external excitation to cause an elevation of temperature field in the material and then to evaluate the resulting temperature field at the surface. However, thermal exciters used (flash lamps, halogen, lasers) act only on the surface of the sample. Several energy conversion systems can on the other hand lead to the generation of volumetric sources: the phenomena of thermo-acoustic, thermo-induction, thermomechanic or thermochemistry can be cited. For example, ultrasonic waves can generate volumetric heat sources if the material is viscoelastic or if there is a defect. The reconstruction of these sources is the first step for the quantification of parameters responsible of the heating. Characterizing a heat source means reconstructing its geometry and the power it generates. For example, a defect in a structure and / or the viscoelasticity of a material can be detected and quantified by this technique if it acts directly on temperature field. However, identification of volumetric heat sources from surface temperature fields is a mathematical ill-posed problem. The diffusive nature of the temperature is the main cause. In this work, the 3D reconstruction of the volumetric heat sources from the resulting surface temperature field, measured by InfraRed, is studied. First, an analysis of the physical problem enables to specify the limits of the reconstruction. In particular, a criterion on achievable spatial resolution is defined and a reconstruction limitation for in-depth sources is highlighted. Then, a probabilistic approach for the reconstruction is proposed and compared to existing inverse methods. The computation time and noise sensitivity are studied for each of these methods. Numerical and experimental applications will thus be presented to illustrate the results.

L'objectif de cette thèse est de mesurer les champs de température et de composition de la stratosphère de Saturne. Pour cela, j'ai analysé des spectres de l'émission thermique de la planète acquis au limbe par CIRS, un spectromètre infrarouge à bord de la sonde Cassini. L'analyse de ces signatures spectrales à l'aide d'un code de transfert de rayonnement et d'une méthode d'inversion permet de mesurer les profils verticaux de la température et de l'abondance des constituants. J'ai ainsi réalisé les premières cartographies, de 70° N à 80° S, et sur plusieurs échelles de hauteur, de la température et de l'abondance de cinq hydrocarbures: l'éthane, l'acétylène, le propane, le méthylacétylène et le diacétylène. La carte de température obtenue permet de quantifier la réponse de l'atmosphère aux variations saisonnières de l'ensoleillement. Les gradients mesurés sont comparés au modèle radiatif saisonnier de Greathouse et al., 2008. Je montre également la découverte d'une oscillation équatoriale de température et de vent zonal dans la stratosphère de Saturne, analogue à l'oscillation quasi-biennale terrestre, et je présente ses variations temporelles entre 2005 et 2010. L'analyse de la distribution des hydrocarbures et la comparaison à un modèle de photochimie saisonnier (Moses et al., 2005) met en évidence des indices de transport méridien dans la basse stratosphère et de transport vertical dans la haute stratosphère. L'abondance du CO2 est également déterminée à 4 latitudes et présente un enrichissement à l'équateur. Enfin, je présente un travail d'analyse des échanges radiatifs dans la stratosphère à l'aide d'une formulation en puissances nettes échangées. Ce travail servira de base à une future paramétrisation du transfert de rayonnement, étape nécessaire au développement d'un modèle de circulation générale de la stratosphère de Saturne.

La microfluidique est une technologie récente dont le principal avantage est de pouvoir miniaturiser le phénomène étudié, quelque soit le domaine (énergétique, chimie, biologie...). Cependant les techniques de mesures (vitesses, températures...) au sein des systèmes microfluidiques sont souvent complexes à mettre en oeuvre. C'est pourquoi il est important de mettre en place un nouvel outil de mesure efficace et fournissant des informations à l'échelle de ces microsystèmes. La Thermographie InfraRouge (TIR) est un moyen puissant de mesurer des champs de températures, qui peuvent ensuite être utilisés pour estimer des propriétés thermo-physiques. Le premier objectif de ce travail de thèse est donc de développer un corps de méthodes expérimentales et numériques couplant la TIR et la microfluidique dans le but d'estimer des paramètres fluidiques et thermiques. Une procédure expérimentale originale ainsi que des méthodes inverses à un paramètre ont été développées. Toutes les méthodes développées ont été testées et validées sur des champs de températures simulés numériquement. Les mesures expérimentales ont permis de cartographier et d'estimer quantitativement des paramètres tels que la vitesse, la diffusivité thermique et des sources de chaleur. Le deuxième objectif est de tester les méthodes développées sur des applications microfluidiques concrètes. Ces tests ont aboutis à d'excellents résultats tels que la mesure d'enthalpie de réaction chimique en microréacteurs, ainsi que l'étude de la cristallisation de l'eau en microcanal.