Les récentes applications des théories de l'inversion à la physique du Globe ont été adaptées à l'échelle du Génie-Civil. Des mesures de durées de propagation entre puits permettent d'obtenir la répartition des vitesses sismiques dans le plan contenant les deux forages. Le problème de l'inversion en facteur de quaiité peut se poser dans les mêmes termes que ceux de l'inversion en vitesse. Des mesures d'amplitude~ permettent d'obtenir la répartition des facteurs de qualité. La mise au point des programmes d'inversion a nécessité une élaboration préalable de programmes de simulation. Les données synthétiques obtenues ont permis de juger de l'influence des divers paramètres sur les résultats des inversions ainsi que de tester les critères de contrôle des résultats. Par la suite, un matériel spécifique a été mis au point afin d'obtenir des données réelles. Les géométries des premières Images calculées, en vitesse et en facteur de qualité sont identiques. Ceci apporte une crédibilité importante aux résultats de l'inversion en facteur de qualité. L'application particulière décrite ici est le contrôle des injections. Des mesures identiques avant et après injection ont permis de déterminer clairement les conséquences diverses de l'injection dans les sous-sols étudiés.

Les moteurs thermoacoustiques sont des oscillateurs autonomes constitués d'un résonateur acoustique partiellement occupé par un matériau poreux (stack) soumis à un important gradient de température grâce à un apport de chaleur externe. Lorsque le gradient de température imposé le long du stack devient supérieur à un certain gradient critique, appelé seuil de déclenchement, l'interaction fluide-parois se traduit par l'amplification d'une onde acoustique auto-entretenue de fort niveau sur le mode le plus instable du résonateur. L'objet des travaux présentés dans ce mémoire est double. D'une part, il est de proposer un formalisme pour la description du fonctionnement de moteurs thermoacoustiques facilement généralisable à l'ensemble de ces systèmes, qu'ils soient à ondes stationnaires ou à ondes progressives. D'autre part, il est de proposer une approche expérimentale pour la caractérisation du noyau thermoacoustique (incluant le stack et la portion de guide inhomogène en température), qui permette de décrire le comportement de systèmes thermoacoustiques sans formuler d'hypothèses sur la forme du champ de température ou la géométrie du stack.Une modélisation analytique des conditions marginales de stabilité et du taux d'amplification de l'onde est tout d'abord proposée, basée sur l'écriture des matrices de transfert des différents éléments qui constituent le moteur. Ces matrices de transfert associées aux conditions aux limites du système étudié conduisent à une équation caractéristique dont la forme dépend de la géométrie de moteur considérée. La solution de cette équation est une pulsation acoustique complexe dont la partie imaginaire correspond au coefficient d'amplification thermoacoustique.La mesure de la matrice de transfert du noyau thermoacoustique constitue la partie expérimentale des travaux exposés. Elle est réalisée pour différentes conditions de chauffage au moyen d'une méthode à quatre microphones. Dans un premier temps, les résultats expérimentaux sont introduits dans le modèle développé précedemment pour prédire le seuil de déclenchement de divers moteurs thermoacoustiques équipés de ce noyau. Les résultats obtenus grâce à cette méthode sont très proches des observations expérimentales, validant ainsi le banc de mesure et le modèle décrivant les conditions de stabilité. Dans un second temps, les données expérimentales sont utilisées pour affiner un modèle analytique décrivant les mécanismes couplés de propagation acoustique et de transport de la chaleur dans le noyau thermoacoustique : ceci permet notamment d'ajuster les valeurs de paramètres acoustiques et thermiques au moyen d'une méthode inverse.Au-delà du seuil de déclenchement, l'amplification et la saturation de l'onde résultent pour une bonne part du transport de chaleur thermoacoustique et de la convection forcée liée à la génération d'un écoulement redressé (vent acoustique), ces deux mécanismes étant généralement responsables d'une dynamique d'évolution complexe de l'amplitude de pression acoustique au cours du régime transitoire. La dernière partie de ces travaux est consacrée à l'introduction de ces deux effets dans le modèle décrit ci-avant, donnant ainsi accès à la description du régime transitoire de l'onde. Une modélisation simplifiée des transports de chaleur associés au vent acoustique de Rayleigh est notamment proposée, qui permet de montrer que cet effet joue vraisemblablement un rôle important dans la dynamique des régimes transitoires observés expérimentalement dans un générateur thermoacoustique quart d'onde.

L'objectif de la thèse est de proposer une stratégie robuste permettant l'identification en dynamique transitoire de loi de comportement gouvernant la réponse du matériau jusqu'à la rupture. Il s'agit d'être à même d'identifier les paramètres matériaux supposés gouverner la rupture dans un contexte où les données expérimentales sont fortement incertaines. Faisant suite à un premier travail où la stratégie avait été élaborée dans un cadre élastique, le travail s'est concentré sur l'extension de la méthode d'identification pour les cas non linéaires tout d'abord la viscoplasticité puis les modéles d'endommagement à taux limités.<br />Les difficultés rencontrées dans ces cas résident dans la non-linéarité et le caractère instable du problème de minimisation sous contraintes non linéaires auquel la formulation nous amène. Une extension de la méthode LATIN aux problèmes mal posés a été proposée et développée afin de permettre la résolution itérative de ce type de problèmes d'optimisation. La résolution de ces derniers fait appel à une méthode de traitement robuste issue du contrôle optimal et basée sur l'équation de Riccati.<br />Une fois ces difficultés résolues et dans les cas simples unidimensionnels traités pour le moment, la stratégie d'identification proposée s'avère très robuste face aux perturbations des mesures même dans le cas très sévère de la localisation et de la rupture.

Ce travail de thèse est consacré au développement de modèles et outils numériques permettant la simulation dans la gamme des fréquences 50 Hz¡100 kHz des phénomènes induits par des champs magnétiques dans le corps humain. Pour cela, il existe deux problématiques principales : - Déterminer la répartition du champ rayonné (ou champ de fuite) par un appareil électrique . - Calculer les courants induits qui en résultent dans le corps humain. Après avoir étudié l'influence de divers paramètres sur la répartition des champs de fuite générés par les systèmes électrotechniques, nous avons recensé les méthodes numériques aptes à calculer ce champ pour des systèmes connus. Cette étude montre que ces méthodes en l'état ne sont pas adaptées au calcul des champs dans l'air, car trop coûteuses. Nous présentons donc un modèle 3D ne prenant en compte que les aspects essentiels de la géométrie du système, beaucoup moins coûteux et suffisamment précis. Dans ce modèle, le champ de fuites est calculé à partir d'une représentation en termes de charges magnétiques fictives. Ce modèle a été validé sur des structures géométriques simples, de type bobine + circuit magnétique + entrefer. L'exposition aux champs en milieu industriel ou domestique est également due à des sources dont la structure est inconnue. Nous avons donc développé un système expérimental, permettant de caractériser le champ de fuites à partir d'un nombre minimum de mesures localisées. Pour cela, nous avons développé divers modèles de la source rayonnante s'appuyant sur les notions de dipôle et multipôles. Nous nous sommes consacrés au divers problèmes d'estimation (influence du choix des mesures, prise en compte du bruit). Le problème étant mal conditionné, des techniques de régularisation ont été mises en oeuvre. La faisabilité a été démontrée sur des structures rayonnantes simples. Enfin, nous présentons une formulation 3D adaptée pour calculer les courants induits dans le corps humain, à partir de la connaissance des champs de fuites. Cette formulation a été implémentée dans un code aux éléments finis, dans lequel le domaine de résolution est limité au seul corps humain : ce code permet de modéliser des structures anatomiques réalistes où les valeurs des conductivités des divers organes sont prises en compte. Différentes situations d'exposition du corps aux rayonnements issus de sources simples filaires, et à des sources réelles ont pu être ainsi simulées, et les applications notamment dans des contextes normatif et industriel sont particulièrement importantes.

L'objectif de ce mémoire de thèse est l'identification et la mise au point de techniques permettant la mesure de la permittivité complexe des bétons dans le domaine de fréquence de fonctionnement des radar géophysiques. Deux techniques types ont été développées : une cellule en transition linge coaxiale/guide cylindrique (50 MHz – 1,6 GHz) et deux sondes en guide d'onde rectangulaire (bandes spectrales retenues : [4 – 7 GHz] et [7 – 13 GHz]). <br /><br />Outre leurs caractères large bande et peu destructif, des simulations numériques réalisées à l'aide d'un logiciel électromagnétique 3D ont montré que ces techniques présentent les avantages d'avoir une grande profondeur de pénétration (supérieure à 6 cm), ainsi qu'une bonne sensibilité de mesure. Elles permettent alors de satisfaire le critère du Volume Élémentaire Représentatif (VER - taille minimale du matériau pour avoir des mesures représentatives).<br /><br />L'étude théorique des sondes/cellule a nécessité l'analyse des problèmes direct et inverse. L'étude du problème direct permet de calculer le coefficient de réflexion en fonction des propriétés diélectriques du matériau à caractériser. Le calcul de la permittivité à partir du coefficient de réflexion mesuré définit le problème inverse. Après la conception des dispositifs de mesure, les logiciels de dépouillement des données et d'étalonnage ont été validés sur des matériaux étalons.<br /><br />Enfin, Les mesures réalisées sur différents mélanges de béton ont montré que la nature des granulats joue un rôle prépondérant sur la variation de la permittivité, que la taille maximale des granulats influence la permittivité dans le domaine centimétrique, et que cette dernière varie linéairement en fonction de la compacité.

Dans le domaine du raffinage, les mesures morphologiques de particules sont devenues indispensables pour caractériser les supports de catalyseurs. À travers ces paramètres, on peut remonter aux spécificités physico-chimiques des matériaux étudiés. Une des techniques d'acquisition utilisées est la tomographie électronique (ou nanotomographie). Des volumes 3D sont reconstruits à partir de séries de projections sous différents angles obtenues par microscopie électronique en transmission (MET). Cette technique permet d'acquérir une réelle information tridimensionnelle à l'échelle du nanomètre. Les projections sont obtenues en utilisant les possibilités d'inclinaison du porte objet d'un MET. À cause des limitations mécaniques de ce porte objet (proximité avec les lentilles magnétiques et déplacement nanométrique), on ne peut acquérir qu'un nombre assez restreint de projections, celles-ci étant limitées à un intervalle angulaire fixe. D'autre part, l'inclinaison du porte objet est accompagnée d'un déplacement mécanique nanométrique non parfaitement contrôlé. Ces déplacements doivent être corrigés après l'acquisition par un alignement des projections suivant le même axe 3D de rotation. Cette étape est un pré-requis à la reconstruction tomographique. Nous suggérons d'utiliser une méthode de reconstruction tomographique par une approche de type "problème inverse". Cette méthode permet d'aligner des projections et de corriger les lacunes de l'acquisition de l'objet observé en introduisant de façon pertinente des informations a priori. Ces informations sont donc basées à la fois sur la physique de l'acquisition (nature physique des images MET, géométrie et limitation spécifique de l'acquisition des projections, etc...) et sur la nature des objets à reconstruire (nombre et répartition des phases, critères morphologiques de type de connexité, etc...). L'algorithme proposé permet de réaliser la reconstruction nanotomographique avec une grande précision et un temps de calculs réduit considérablement par rapport à la technique classique. Nous avons testé avec succès notre méthode pour les projections réelles de différents supports de catalyseur

Les travaux présentés dans cette thèse se focalisent sur l'étude de l'effet de dégradations des composants de puissance, plus particulièrement au niveau de l'environnement proche des puces (métallisations, connexions, brasures puces/DCB), sur le comportement électrique et thermique des puces ainsi que de leur assemblage. Pour ce faire nous avons cherché à étudier la répartition des courants et des températures à la surface de la puce à l'aide d'un modèle électrothermique 2D distribué. Nous avons aussi évalué l'effet de la dégradation des brasures dans le volume de l'assemblage, à l'aide cette fois d'un modèle thermique relié à la constitution de l'assemblage. La première partie de cette thèse consiste à mettre en place un modèle électrothermique distribué de puce MOSFET, qui tient compte à la fois du caractère distribué de la dissipation de la puissance et du couplage électrothermique en régime transitoire. Ce modèle électrothermique s'appuie sur un modèle électrique aux variables d'états et un modèle thermique par éléments finis couplé au modèle électrique. Les modèles électriques et thermiques ont été développés respectivement sous Matlab et sous CAST3M, et le couplage des deux modèles a été fait sous Simulink. Dans une deuxième partie, pour la validation des résultats des températures et pour l'analyse de l'effet du vieillissement et des dégradations (sur la distribution et la dynamique de température de la surface supérieure de la puce), une méthodologie de mesure rapide de température et un banc expérimental pour thermographie infrarouge ont été mis en place. Les difficultés rencontrées lors des mesures thermiques IR sous variation rapide de la température nous ont poussé à envisager d'autres méthodes d'analyse thermique. Enfin, nous avons cherché à évaluer la réponse impulsionnelle du composant testé en estimant, par des simulations thermiques, la fonction de transfert dans le domaine fréquentiel à l'aide du logiciel COMSOL Multiphysics. Nous avons également étudié la pertinence de modèles RC équivalents (réseau RC de Cauer). Ces modèles ont ensuite été utilisés pour rendre compte de différents modes de dégradation notamment cette fois au niveau des couches de brasures entre puce et DCB et entre DCB et semelle. Mots clef : Modules de puissance à semi-conducteur, Vieillissement, Métallisation, Modélisation électrothermique, Court-circuit, Distribution de courant et de température, Problème inverse, Caméra IR, Réseaux de Cauer.

Ce travail de thèse apporte des éléments méthodologiques pour l'interprétation rapide de données électromagnétiques transitoires (TEM) aéroportées. Dans un premier temps, nous nous sommes intéressés au problème de l'inversion 1D qui est, à ce jour, un traitement standard. Nous nous sommes ensuite focalisés sur le problème de l'inversion 2D rapide. En nous basant sur les résultats de modélisations numériques, nous proposons un modèle empirique de sensibilité pour effectuer une inversion 2D rapide. Les tests effectués sur des données synthétiques et réelles s'avèrent très prometteurs. Ensuite, nous étudions analytiquement l'effet de la topographie. Pour un demi-espace homogène incliné, nous montrons que la direction des courants induits dans le sous-sol dépend de l'angle de l'interface air/sol. Nous proposons alors une méthode pour prendre en compte cet effet. Enfin, nous testons nos méthodes de traitement sur des données VTEM acquises dans le bassin de Franceville au Gabon.

L'imagerie moléculaire occupe une place majeure dans le domaine de la recherche préclinique. Parmi les modalités existantes, les techniques optiques fondées sur la détection d'un rayonnement visible ou proche infrarouge sont les plus récentes et sont principalement représentées par les méthodes de tomographie optique de luminescence. Ces méthodes permettent une caractérisation 3D d'un milieu biologique par la reconstruction de cartes de concentration ou la localisation de marqueurs luminescents sensibles à des processus biologiques et chimiques se déroulant à l'échelle de la cellule ou de la molécule. La tomographie de luminescence se fonde sur un modèle de propagation de la lumière dans les tissus, un protocole d'acquisition du signal en surface du milieu et une procédure numérique d'inversion de ces mesures afin de reconstruire les paramètres d'intérêts. Ce travail de thèse s'articule donc autour de ces trois axes et apporte un élément de réponse à chacun des problèmes. L'objectif principal de cette étude est d'introduire et de présenter des outils d'évaluation des performances théoriques d'une méthode de tomographie optique. L'un des aboutissements majeurs est la réalisation de reconstructions tomographiques expérimentales à partir d'images acquises par un imageur optique conçu pour l'imagerie planaire 2D et développé par la société Quidd. Dans un premier temps nous abordons la théorie du transport en milieu diffusant afin de poser les concepts et outils sur lesquels vont s'appuyer l'ensemble des travaux. Nous présentons particulièrement deux modèles de propagation différents ainsi que les méthodes de résolution et les difficultés théoriques qui leur sont liées. Dans une deuxième partie nous introduisons les outils statistiques utilisés pour caractériser les systèmes tomographiques et leur résolution potentielle. Nous définissons une procédure et nous l'appliquons à l'étude de quelques situations simples en tomographie de luminescence. La dernière partie de ce travail présente la mise au point d'une procédure d'inversion. Après avoir présenté le cadre théorique dans lequel cette procédure s'inscrit nous la validons à partir de données numériques avant de l'appliquer avec succès à des mesures expérimentales.

La tomographie des temps de première arrivée cherche à estimer un modèle de vitesse de propagation des ondes sismiques à partir des temps de première arrivée pointés sur les sismogrammes. Le modèle de vitesse obtenu peut alors permettre une interprétation structurale du milieu ou bien servir de modèle initial pour d'autres traitements de l'imagerie sismique. Les domaines d'application de cette méthode s'étendent, à des échelles différentes, de la géotechnique à la sismologie en passant par la géophysique pétrolière. Le savoir-faire du géophysicien joue un rôle important dans la difficile résolution du problème tomographique non-linéaire et mal posé. De nombreuses recherches ont entrepris de faciliter et d'améliorer cette résolution par des approches mathématique ou physique. Dans le cadre de ce travail, nous souhaitons développer une approche pragmatique, c'est-à-dire que nous considérons que le problème tomographique doit être résolu par un algorithme interactif dont les paramètres de réglage sont clairement définis. L'aspect interactif de l'algorithme facilite l'acquisition du savoir-faire tomographique car il permet de réaliser, dans un temps raisonnable, de nombreuses simulations pour des paramétrisations différentes. Le but poursuivi dans cette thèse est de définir, pour le cas spécifique de la tomographie des temps de première arrivée, un algorithme qui réponde au mieux à ces critères. Les algorithmes de tomographie des temps de première arrivée classiquement mis en oeuvre aujourd'hui ne répondent pas à nos critères d'une approche pragmatique. En effet, leur implémentation ne permet pas d'exploiter l'architecture parallèle des supercalculateurs actuels pour réduire les temps de calcul. De plus, leur mise en oeuvre nécessite une paramétrisation rendue complexe du fait de la résolution du système linéaire tomographique. Toutes ces limitations pratiques sont liées à la formulation même de l'algorithme à partir de la méthode de Gauss-Newton. Cette thèse repose sur l'idée de formuler la résolution du problème tomographique à partir de la méthode de plus grande descente pour s'affranchir de ces limitations. L'étape clé de cette formulation réside dans le calcul du gradient de la fonction coût par rapport aux paramètres du modèle. Nous utilisons la méthode de l'état adjoint et une méthode définie à partir d'un tracé de rais a posteriori pour calculer ce gradient. Ces deux méthodes se distinguent par leur formulation, respectivement non-linéaire et linéarisée, et par leur mise en oeuvre pratique. Nous définissons ensuite clairement la paramétrisation du nouvel algorithme de tomographie et validons sur un supercalculateur ses propriétés pratiques : une parallélisation directe et efficace, une occupation mémoire indépendante du nombre de données observées et une mise en oeuvre simple. Finalement, nous présentons des résultats de tomographie pour des acquisitions de type sismique réfraction, 2-D et 3-D, synthétiques et réelles, marines et terrestres, qui valident le bon comportement de l'algorithme, en termes de résultats obtenus et de stabilité. La réalisation d'un grand nombre de simulations a été rendue possible par la rapidité d'exécution de l'algorithme, de l'ordre de quelques minutes en 2-D.