Depuis plus d’un siècle, la radiographie sur film est utilisée pour le contrôle non destructif (CND) de pièces industrielles. Avec l’introduction de méthodes numériques dans le domaine médical, la communauté du CND industriel a commencé à considérer également les techniques numériques alternatives au film. La radiographie numérique (en anglais Computed radiography -CR) utilisant les écrans photostimulables (en anglais imaging plate -IP) est une voie intéressante à la fois du point de vue coût et facilité d’implémentation. Le détecteur (IP) utilisé se rapproche du film car il est flexible et réutilisable. L’exposition de l’IP aux rayons X génère une image latente qui est ensuite lue et numérisée grâce à un système de balayage optique par laser. A basse énergie, les performances du système CR sont bonnes ce qui explique son utilisation importante dans le domaine médical. A haute énergie par contre, les performances du système CR se dégradent à la fois à cause de la mauvaise absorption de l’IP mais également de la présence de rayonnement diffusé par la pièce qui, étant d’énergie plus faible, est préférentiellement absorbée par l’IP. Les normes internationales préconisent l’utilisation d’écrans métalliques pour améliorer la réponse des systèmes CR à haute énergie. Néanmoins, la nature et l’épaisseur de ces écrans n’est pas clairement définie et la gamme des configurations possibles est large. La simulation est un outil utile pour prévoir les performances d’une expérience et déterminer les meilleures conditions opératoires. Les méthodes Monte Carlo sont communément admises comme étant les plus précises pour simuler les phénomènes de transport de rayonnement, et ainsi comprendre les phénomènes physiques en jeu. Cependant, le caractère probabiliste de ces méthodes implique des temps de calcul importants, voire prohibitifs pour des géométries complexes. Les méthodes déterministes au contraire, peuvent prendre en compte des géométries complexes avec des temps de calcul raisonnables, mais l’estimation du rayonnement diffusé est plus difficile. Dans ce travail de thèse, nous avons tout d’abord mené une étude de simulation Monte Carlo afin de comprendre le fonctionnement des IP avec écrans métalliques à haute énergie pour le contrôle de pièces de forte épaisseur. Nous avons notamment suivi le trajet des photons X mais également des électrons. Quelques comparaisons expérimentales ont pu être menées à l’ESRF (European Synchrotron Radiation Facility). Puis nous avons proposé une approche de simulation hybride, qui combine l'utilisation de codes déterministe et Monte Carlo pour simuler l'imagerie d'objets de forme complexe. Cette approche prend en compte la dégradation introduite par la diffusion des rayons X et la fluorescence dans l'IP ainsi que la diffusion des photons optiques dans l'IP. Les résultats de différentes configurations de simulation ont été comparés. / For over a century, film-based radiography has been used as a nondestructive testing technique for industrial inspections. With the advent of digital techniques in the medical domain, the NDT community is also considering alternative digital techniques. Computed Radiography (CR) is a cost-efficient and easy-to-implement replacement technique because it uses equipment very similar to film radiography. This technology uses flexible and reusable imaging plates (IP) as a detector to generate a latent image during x-ray exposure. With an optical scanning system, the latent image can be readout and digitized resulting in a direct digital image. CR is widely used in the medical field since it provides good performance at low energies. For industrial inspection, CR application is limited by its poor response to high energy radiation and the presence of scattering phenomena. To completely replace film radiography by such a system, its performance still needs to be improved by either finding more appropriate IPs or by optimizing operating conditions. Guidelines have been addressed in international standards to ensure a good image quality supplied by CR system, where metallic screens are recommended for the case of using high energy sources. However, the type and thickness of such a screen are not clearly defined and a large panel of possible configurations does exist. Simulation is a very useful tool to predict experimental outcomes and determine the optimal operating conditions. The Monte Carlo (MC) methods are widely accepted as the most accurate method to simulate radiation transport problems. It can give insight about physical phenomena, but due to its random nature, a large amount of computational time is required, especially for simulations involving complex geometries. Deterministic methods, on the other hand, can handle easily complex geometry, and are quite efficient. However, the estimation of scattering effects is more difficult with deterministic methods. In this thesis work, we have started with a Monte Carlo simulation study in order to investigate the physical phenomena involved in IP and in metallic screens at high energies. In particular we have studied separately the behavior of X-ray photons and electrons. Some experimental comparisons have been carried out at the European Synchrotron Radiation Facility. Then, we have proposed a hybrid simulation approach, combining the use of deterministic and Monte Carlo code, for simulating the imaging of complex shapes objects. This approach takes into account degradation introduced by X-ray scattering and fluorescence inside IP, as well as optical photons scattering during readout process. Different simulation configurations have been compared.

La tomographie par rayons X ou CT pour "Computed Tomography" est un outil puissant pour caractériser et localiser les défauts internes et pour vérifier la conformité géométrique d’un objet. Contrairement au cas des applications médicales, l’objet inspecté en Contrôle Non Destructif (CND) peut être très grand et composé de matériaux de haute atténuation, auquel cas l’utilisation d’une trajectoire circulaire pour l’inspection est impossible à cause de contraintes dans l’espace. Pour cette raison, l’utilisation de bras robotisés est l’une des nouvelles tendances reconnues dans la CT, car elle autorise plus de flexibilité dans la trajectoire d’acquisition et permet donc la reconstruction 3D de régions difficilement accessibles dont la reconstruction ne pourrait pas être assurée par des systèmes de tomographie industriels classiques. Une cellule de tomographie X robotisée a été installée au CEA. La plateforme se compose de deux bras robotiques pour positionner et déplacer la source et le détecteur en vis-à-vis. Parmi les nouveaux défis posés par la tomographie robotisée, nous nous concentrons ici plus particulièrement sur la limitation de l’ouverture angulaire imposée par la configuration en raison des contraintes importantes sur le mouvement mécanique de la plateforme. Le deuxième défi majeur est la troncation des projections qui se produit lorsque l’objet est trop grand par rapport au détecteur. L’objectif principal de ce travail consiste à adapter et à optimiser des méthodes de reconstruction CT pour des trajectoires non standard. Nous étudions à la fois des algorithmes de reconstruction analytiques et itératifs. Avant d’effectuer des inspections robotiques réelles, nous comptons sur des simulations numériques pour évaluer les performances des algorithmes de reconstruction sur des configurations d’acquisition de données. Pour ce faire, nous utilisons CIVA, qui est un outil de simulation pour le CND développé au CEA et qui est capable de simuler des données de projections réalistes correspondant à des configurations d’acquisition définies par l’utilisateur. / X-ray computed tomography (CT) is a powerful tool to characterize or localize inner flaws and to verify the geometric conformity of an object. In contrast to medical applications, the scanned object in non-destructive testing (NDT) might be very large and composed of high-attenuation materials and consequently the use of a standard circular trajectory for data acquisition would be impossible due to constraints in space. For this reason, the use of robotic arms is one of the acknowledged new trends in NDT since it allows more flexibility in acquisition trajectories and therefore could be used for 3D reconstruction of hardly accessible regions that might be a major limitation of classical CT systems. A robotic X-ray inspection platform has been installed at CEA LIST. The considered system integrates two robots that move the X-ray generator and detector. Among the new challenges brought by robotic CT, we focus in this thesis more particularly on the limited access viewpoint imposed by the setup where important constraints control the mechanical motion of the platform. The second major challenge is the truncation of projections that occur when only a field-of-view (FOV) of the object is viewed by the detector. Before performing real robotic inspections, we highly rely on CT simulations to evaluate the capability of the reconstruction algorithm corresponding to a defined scanning trajectory and data acquisition configuration. For this purpose, we use CIVA which is an advanced NDT simulation platform developed at CEA and that can provide a realistic model for radiographic acquisitions and is capable of simulating the projection data corresponding to a specific CT scene defined by the user. Thus, the main objective of this thesis is to develop analytical and iterative reconstruction algorithms adapted to nonstandard trajectories and to integrate these algorithms in CIVA software as plugins of reconstruction.

Ce mémoire de thèse étudie diverses méthodes d'amélioration des interféromètres atomiques. Dans la première partie du manuscrit, nous analysons comment une détection non-destructive, au sens où elle préserve la cohérence entre les états internes de l'ensemble atomique, permet d'améliorer la sensibilité des interféromètres. Nous montrons tout d'abord, grâce à une étude théorique, que la projection du vecteur d'onde engendrée par la mesure permet de préparer des états comprimés de spin. Nous présentons ensuite la mise en œuvre de cette méthode à l'aide d'une détection reposant sur la spectroscopie par modulation de fréquence. Finalement, nous exposons quelques premières applications de cette détection non-destructive, plus précisément nous présentons la réalisation du rétroaction quantique qui protège l'état atomique contre la décohérence induite par un basculement du spin collectif, nous montrons aussi comment réaliser une boucle à verrouillage de phase où les atomes servent de référence de phase. Dans la seconde partie du manuscrit, nous présentons la réalisation tout-optique d'un condensat de Bose-Einstein dans une cavité de haute finesse, exploitant les technologies développées pour les télécommunications optiques. Nous commençons par une analyse du résonateur et des méthodes d'asservissement, nous introduisons notamment une méthode d'asservissement originale exploitant la modulation serrodyne. Enfin, nous montrons comment un condensat est obtenu par évaporation dans le mode optique de la cavité. / In this thesis, we study several methods to improve atom interferometers. In the first part of the manuscript, we analyze how a nondestructive detection, that preserves the coherence between the internal degrees of freedom in an atomic ensemble, can be used to increase the sensitivity of interferometers. We first theoretically show how the projection of the wave-function induced by the measurement prepares spin-squeezed states. We then present the implementation of this method with a detection based on the frequency modulation spectroscopy. Finally, some first applications are described, more explicitly we show how to implement a quantum feedback that preserve the atomic state against the decoherence induced by a random collective flip, we also introduce a phase-locked loop where the atomic sample is used as the phase reference. In the second part of the manuscript, we present the all-optical realization of a Bose-Einstein condensate in a high-finesse cavity using a laser system based on standard telecoms technologies. We first describe the resonator and the frequency lock of the laser on the resonance, in particular, we introduce a new stabilization method based of the serrodyne modulation. Finally, we show how the condensate is obtained from the evaporation in the cavity mode.

Este trabalho descreve o desenvolvimento de um algoritmo para avaliar o grau de corrosão em paredes de oleodutos, através de ultrassom. O algoritmo desenvolvido será implementado em um circuito digital, fazendo parte de um sistema que poderá ser embarcado em pigs ultrassônicos para inspeção de dutos à procura de corrosão. O modelo matemático conhecido como filtro casado é usado como base do algoritmo, além de um método de detecção de picos para localização dos sinais de eco de ultrassom, provenientes das reflexões do sinal nas paredes dos dutos. Um dos problemas é a definição do sinal (e do filtro) que será utilizado, mas por meio de sinais sintetizados (signal design) uma grande gama de alternativas foi analisada e apenas algumas foram selecionadas para testes em experimentos. Por último os melhores sinais foram utilizados na inspeção de um corpo de prova feito à partir de um duto corroído, e os resultados são analisados. A detecção da primeira parede do oleoduto apresenta resultados promissores em comparação com a medição tridimensional das mesmas paredes. O sistema foi capaz de reconstruir o perfil da parede interna do duto com elevada verossimilhança, mesmo em regiões severamente corroídas. A detecção da segunda parede apresenta resultados satisfatórios para baixos níveis de corrosão, porém alguns problemas foram encontrados para níveis mais altos de corrosão, prejudicando a medição direta da espessura da parede nestas condições. / This work describes the development of an algorithm to evaluate the degree of corrosion in pipeline walls, using ultrasound. The developed algorithm will be implemented in a dedicated digital circuit that can be embedded in ultrasonic pigs for pipeline corrosion inspection. The algorithm is based on the mathematic model known as matched filter, besides a peak detection method for the location of ultrasound echoes, originated from the signal reflections in the pipeline walls. One problem is the determination of the signal (and filter) that should be used, but using a signal design approach a large number of possibilities could be analyzed and only a few was tested in experiments. Finally, the best signals were used to inspect a test specimen made of a corroded pipeline section, and the results were analyzed. The first wall detection shows promising results when compared to the results of a three-dimensional inspection of the same walls. It was possible to reconstruct the profile of the internal wall with good likelihood, even at badly corroded areas. The second wall detection shows some satisfactory results on low corrosion levels, but show some problems with high corrosion levels, compromising the direct measurement of the wall of the pipeline.

Ce travail se situe à l’intersection des trois disciplines : méthodes numériques, techniques expérimentales et du machine learning, a pour but de proposer une famille de techniques d’identification par analyse inverse des lois de comportement en mécanique. Dans le domaine d’identification des matériaux, l’indentation instrumentée est particulièrement attractive, car elle permet de procéder à des essais non-destructifs sur l’échantillon ou sur une structure en service. L’essai d’indentation, similaire à un test de dureté, consiste à enfoncer la pointe de l’indenteur à une faible profondeur dans la matière tout en enregistrant le déplacement en fonction de la force appliquée. L’identification des propriétés élastoplastiques des matériaux est basée alors sur l’exploitation de la courbe force-déplacement (courbe P-h). Toutefois, le problème inverse est souvent mal posé et des problèmes d’unicité mènent à la notion de paires de "matériaux mystiques" produisant, dans des conditions d’essai donnés, des courbes P-h identiques, malgré des propriétés différentes. L’idée de notre travail est de compléter la procédure d’identification en faisant appel à des dispositifs expérimentaux récents, notamment à la microscopie laser, permettant de mesurer la carte 3D de l’empreinte résiduelle obtenue après le retrait de l’indenteur. Pour aborder la question de la richesse d’information de l’empreinte par rapport à la courbe P-h seule, nous proposons de construire, dans un espace affine réduit, la variété des formes d’empreinte admissibles au sens d’une loi de comportement et du modèle d’éléments finis de l’essai. La mesure de la dimension intrinsèque nous indique alors le nombre maximal de paramètres potentiellement identifiables. Cela nous permet de proposer et de valider numériquement des nouveaux procédés expérimentaux, plus représentatifs, à partir des données synthétiques, ainsi que des algorithmes d’identification associés. La prise en compte de l’erreur de modèle et de l’erreur de mesure, nous mène ensuite à proposer un ensemble d’algorithmes de projection d’empreintes expérimentales, réalisées en collaboration avec l’INSA de Rennes sur la variété synthétique. Nous abordons alors le problème d’identification des propriétés d’écrouissage de plusieurs matériaux de complexité croissante et départageons des "jumeaux mystiques" par des essais de multi-indentation, basés sur l’exploitation de l’empreinte seule ou en complément de la courbe P-h. / The thesis lies at the intersection of three disciplines : numerical methods, experimental techniques, and machine learning. The primary aim of this work is to develop a group of algorithms for characterization by inverse analysis of a material’s constitutive law. In the field of material characterization, indentation test is especially attractive since it is considered non-destructive, and may be performed even on a structure in service. The test, similar to a hardness test, consists in penetrating an indenter into the surface of the material. The force exerted on the indenter is recorded against the penetration depth over a series of time instants, leading to a force-displacement (P-h) curve, which is the most frequently used source of information for the identification of material properties. However, the inverse problem based solely on this curve tends to be ill-posed, leading to nonunique identification solution, i.e., the "mystical material pair", for whom the corresponding force-displacement curves are almost identical despite the very different material properties. The basic idea is then to complete the identification process with innovative experimental measurements, such as laser microscope, which allows measuring the 3D residual imprint after the withdrawal of the indenter. To address the advantage of this measurement over P-h curve, we propose to construct, within a reduced affine space, a manifold of shapes admissible to the postulated constitutive law, experimental and simulation setups, based on synthetic data. The intrinsic dimensionality of the manifold limits the number of identifiable parameters allowing to validate numerically experimental procedures. Considering both the model and measurement errors, we develop a series of local manifold learning algorithms to solve the inverse problem iteratively for experimental results obtained in cooperation with INSA de Rennes. This approach allows us to characterize diverse metallic materials of increasing complexity, based on actual experimental measurements. For example, for the Hollomon’s law, the mystical pair is alleviated in using a single imprint, while for the Voce law, a multi-depth experimental protocol is proposed to differentiate mystical siblings.

Cette thèse porte sur le développement d'un traducteur ultrasonore multi-élément capacitif à couplage air (MEACUT) et son utilisation dans le domaine du contrôle non destructif (CND) de matériaux. Un modèle est employé pour simuler numériquement ce traducteur, et pour optimiser sa conception. Un prototype est ensuite fabriqué, puis caractérisé expérimentalement pour quantifier ses performances. Son originalité réside dans le fait qu'il possède une large bande passante en fréquence, tout en offrant la possibilité d'une focalisation variable. Ce prototype est alors employé pour la détection d'un endommagement causé par impact, dans une plaque composite. Il est clairement constaté que la résolution spatiale du procédé d'inspection employé (C-scan) est fortement améliorée grâce aux performances techniques du MEACUT. Enfin, un modèle hybride 3D est développé pour simuler, rapidement et intégralement, ce procédé de CND. Le très bon accord obtenu entre prédictions numériques et mesures expérimentales laisse présager que cet outil de simulation pourra servir à mettre au point d'autres expérimentations de CND, qui pourront à leur tour exploiter les performances du MEACUT. / This thesis focuses on the development of a Multi-Element Air-coupled Capacitive Ultrasonic Transducer (MEACUT) and its use in the field of non-destructive testing (NDT) of materials. A numerical model is used to simulate the translator, and to optimize its design. A prototype is then built and experimentally characterized to quantify its performance. Its originality lies in the fact that it has a broad frequency bandwidth while offering the possibility of a dynamic focusing. This prototype is then used for the detection of damage caused by impact, in a composite plate. It is clearly found that the spatial resolution of the inspection process employed (C-scan) is greatly improved thanks to the technical performance MEACUT. Finally, a 3D hybrid model is developed to simulate quickly and fully, the process of NDT. The good agreement obtained between numerical predictions and experimental measurements suggests that this simulation tool can be used to develop other NDT experiments, which may in turn exploit the performance of MEACUT.

Le secteur du bâtiment représente 35% des la consommations énergétiques dans les pays membres de l’agence international de l’énergie en 2010 et 39,8% aux Etats-Unis en 2015. Plus de 50% de cette consommation a été utilisée pour la production de chaleur et de froid. Néanmoins cette consommation peut être réduite par l'amélioration la performance énergétique du bâtiment. La performance thermique de l'enveloppe du bâtiment joue un rôle primordial. Par conséquent, le diagnostic thermique de l'enveloppe du bâtiment est nécessaire pour, par exemple, la réception de nouvelles constructions, l'amélioration de la performance énergétique des anciens bâtiments, ainsi que la vente et la location des logements. Pourtant, il existe très peu de méthodes quantitatives pour la caractérisation des parois épaisses. L'objectif de cette étude est d'explorer des méthodes quantitatives innovantes de diagnostic thermique de l'enveloppe du bâtiment. Des mesures expérimentales ont été réalisées en laboratoire (à l’IFSTTAR à Nantes) et in situ (à l’IUT de Bordeaux). Différents capteurs et méthodes d'instrumentation ont été étudiés pour mesurer la densité de flux et la température de surfaces des parois, afin de procurer des recommandations pour le choix des capteurs ainsi que des protocoles de traitement de données. A partir des données mesurées (température et densité de flux des surfaces de l'enveloppe), trois approches numériques ont été proposées pour estimer des paramètres thermiques des parois multicouches épaisses : par méthode inverse, par réponse à un échelon et par réponse impulsionnelle. En outre, une méthode innovante non-destructive utilisant la rayonnement térahertz a été étudiée. Les mesures ont été effectuées au sein du laboratoire I2M. Cette méthode permet de caractériser le coefficient d'absorption des matériaux constructifs ordinaires comme isolation, plâtre, béton, bois… Elle pourrait postérieurement être combinée avec une méthode thermique pour apporter des informations complémentaires. / Buildings represent a large share in terms of energy consumption, such as 35% in the member countries of IEA (2010) and 39.8% in U.S. (2015). Climate controlling (space heating and space cooling) occupies more than half of the consumption. While this consumption can be reduced by improving the building energy efficiency, in which the thermal performance of building envelope plays a critical role. Therefore, the thermal diagnosis of building envelope is of great important, for example, in the case of new building accreditation, retrofitting energy efficiency of old building and the building resale and renting. However, very few diagnostic methods exist for the characterization of thick walls. The present measurement standards that based on steady state heat transfer regime need a long time (several days). The classical transient technologies, such as flash method, are difficult to implement on the walls because of the large thickness of walls and the complex conditions in situ. This thesis aims to explore innovative methodologies for thermal quantitative diagnosis of building envelope. Two experimental cases were carried out: one is in laboratory (IFSTTAR, Nantes) and the other is in situ (IUT, Bordeaux). Different sensors and instruments were studied to measure the wall heat flux and surface temperature, and provided some guidelines for the choice of sensors and data processing protocols as well. Using these measured data, three estimation approaches were proposed to estimate the thermal parameters of the multilayer thick wall: pulse response curve method, step response curve method and inverse method, which can be applied for different diagnostic situations. In addition, an innovative NDE (non-destructive evaluation) method using terahertz (THz) radiation was also investigated. Measurements were carried out in I2M laboratory to characterize the absorption coefficient of standard building materials (insulation, plaster, concrete, wood ...). This THz method can be combined with a previous thermal method to provide some complementary information.

Le principe du contrôle par flux de fuite magnétique (Magnetic Flux Leakage MFL) consiste à aimanter la pièce à contrôler par un champ magnétique et à détecter à l'aide d'un capteur magnétique les fuites des lignes du champ qui résultent de la présence d'un défaut dans la pièce. Dans le but d'améliorer les performances d'un dispositif de détection, le CEA et la société Vallourec collaborent pour développer un modèle numérique dédié au contrôle virtuel des défauts longitudinaux dans les tubes ferromagnétiques. Le dispositif expérimental comprend un circuit magnétique tournant à une vitesse constante autour du tube qui défile. Dans le cadre de cette thèse, on débute le problème de la modélisation sans tenir compte des effets de la vitesse de rotation, il s'agit donc de résoudre un problème d'électromagnétisme en régime magnétostatique.Pour résoudre ce problème, on propose de comparer une approche semi-analytique basée sur le formalisme des équations intégrales (EI) et une approche purement numérique utilisant les éléments finis (EF).Dans la première partie de cette thèse, après avoir établi le formalisme théorique par EI, un premier modèle considérant des matériaux ferromagnétiques à perméabilité magnétique constante (régime linéaire) a été mis en œuvre en 2D. Ce modèle a été appliqué pour un exemple de système extrait de la littérature et validé numériquement par une comparaison des résultats EI/EF. Pour une meilleure détection, il est opportun de saturer magnétiquement la pièce. Le matériau ferromagnétique est alors caractérisé par une courbe B(H) non-linéaire. Par conséquent, la deuxième partie de la thèse a été consacrée à la mise en œuvre du modèle en régime non linéaire qui tient compte de cette caractéristique.Différentes méthodes de discrétisation ont été étudiées afin de réduire le nombre d'inconnues et le temps de calcul. L'originalité de la thèse réside dans l'utilisation des fonctions d'interpolation d'ordre élevé (polynôme de Legendre) pour une discrétisation des équations intégrales par une approche de type Galerkin. Les premiers essais de validation numérique de ce modèle ont été effectués sur un système MFL simplifié. Des premiers essais de validation expérimentale pour des données obtenues par EF ont été effectués en deux phases : La première a consisté à vérifier le distribution du champ magnétique pour un tube sain et en régime magnétostatique. La deuxième phase a consisté à calculer la réponse d'un défaut dans le tube ferromagnétique en tenant en compte les effets éventuels de la rotation du circuit magnétique par rapport au tube. / The principle of the non destructive testing by magnetic flux leakage (MFL) is to magnetize the part to be inspected by a magnetic field and to detect a flaw thanks to magnetic leakage field lines due to the strong decreasing of the magnetic permeability in the flawed region. In order to improve the performance of detection, the CEA and the Vallourec society collaborate to develop a numerical model dedicated to the virtual NDT of longitudinal defects in ferromagnetic tubes. The experimental system includes a magnetic circuit rotating at a constant speed around the tube to be inspected. The modeling task is started without considering the effects of the rotational speed, so the magnetostatic regime is considered to solve the modeling problem. In the framework of this thesis, we propose to compare a semi-analytical approach based on the formalism of integral equations method (IEM) and a purely numerical approach using finite element method (FEM).In the first part of this thesis, the theoretical formalism was established. A first simple discretization scheme is been implemented in the linear regime considering a constant magnetic permeability. This first numerical model has been validated for a simplified MFL configuration extracted and modified from the literature.For better detection, it is wishable to magnetically saturate the piece under-test. The ferromagnetic material is then characterized by a B(H) curve. Therefore, the second part of the thesis was devoted to the implementation of the model in the non-linear regime that takes into account this non-linear characteristic. Different discretization schemes have been studied in order to reduce the number of unknowns and the computational time. The originality of the thesis lies in the use of basis function of high order (Legendre polynomials) associated to a Galerkin approach for the discretization of integral equations. The first numerical result has been validated on a simplified MFL system. The first results of the experimental validation based on simulated data obtain by FEM have been performed in two steps. The first one consists to verify the distribution of the magnetic field for a ferromagnetic tube without any defect and in the magnetostatic regime. The objective of the second one was to compute the response of the flaw and to evaluate the effects of the rotational speed of the magnetic circuit around the tube.

La problématique de cette thèse se place à l’interface entre le domaine scientifique du contrôle non destructif par ultrasons (CND US) et l’adéquation algorithme architecture. Le CND US comprend un ensemble de techniques utilisées pour examiner un matériau, qu’il soit en production ou maintenance. Afin de détecter d’éventuels défauts, de les positionner et les dimensionner, des méthodes d’imagerie et de reconstruction ont été développées au CEA-LIST, dans la plateforme logicielle CIVA.L’évolution du matériel d’acquisition entraine une augmentation des volumes de données et par conséquent nécessite toujours plus de puissance de calcul pour parvenir à des reconstructions en temps interactif. L’évolution multicoeurs des processeurs généralistes (GPP), ainsi que l’arrivée de nouvelles architectures comme les GPU rendent maintenant possible l’accélération de ces algorithmes.Le but de cette thèse est d’évaluer les possibilités d’accélération de deux algorithmes de reconstruction sur ces architectures. Ces deux algorithmes diffèrent dans leurs possibilités de parallélisation. Pour un premier, la parallélisation sur GPP est relativement immédiate, contrairement à celle sur GPU qui nécessite une utilisation intensive des instructions atomiques. Quant au second, le parallélisme est plus simple à exprimer, mais l’ordonnancement des nids de boucles sur GPP, ainsi que l’ordonnancement des threads et une bonne utilisation de la mémoire partagée des GPU sont nécessaires pour obtenir un fonctionnement efficace. Pour ce faire, OpenMP, CUDA et OpenCL ont été utilisés et comparés. L’intégration de ces prototypes dans la plateforme CIVA a mis en évidence un ensemble de problématiques liées à la maintenance et à la pérennisation de codes sur le long terme. / This thesis work is placed between the scientific domain of ultrasound non-destructive testing and algorithm-architecture adequation. Ultrasound non-destructive testing includes a group of analysis techniques used in science and industry to evaluate the properties of a material, component, or system without causing damage. In order to characterize possible defects, determining their position, size and shape, imaging and reconstruction tools have been developed at CEA-LIST, within the CIVA software platform.Evolution of acquisition sensors implies a continuous growth of datasets and consequently more and more computing power is needed to maintain interactive reconstructions. General purprose processors (GPP) evolving towards parallelism and emerging architectures such as GPU allow large acceleration possibilities than can be applied to these algorithms.The main goal of the thesis is to evaluate the acceleration than can be obtained for two reconstruction algorithms on these architectures. These two algorithms differ in their parallelization scheme. The first one can be properly parallelized on GPP whereas on GPU, an intensive use of atomic instructions is required. Within the second algorithm, parallelism is easier to express, but loop ordering on GPP, as well as thread scheduling and a good use of shared memory on GPU are necessary in order to obtain efficient results. Different API or libraries, such as OpenMP, CUDA and OpenCL are evaluated through chosen benchmarks. An integration of both algorithms in the CIVA software platform is proposed and different issues related to code maintenance and durability are discussed.

Les travaux de cette thèse portent sur la détection et le diagnostic des défauts naissants dans les systèmes d’ingénierie et industriels, par des approches statistiques non-paramétriques. Un défaut naissant est censé provoquer comme tout défaut un changement anormal dans les mesures des variables du système. Ce changement est imperceptible mais aussi imprévisible dû à l’important rapport signal-sur défaut, et le faible rapport défaut-sur-bruit caractérisant le défaut naissant. La détection et l’identification d’un changement général nécessite une approche globale qui prend en compte la totalité de la signature des défauts. Dans ce cadre, la divergence de Kullback-Leibler est proposée comme indicateur général de défauts, sensible aux petites variations anormales cachées dans les variations du bruit. Une approche d’analyse spectrale globale est également proposée pour le diagnostic de défauts ayant une signature fréquentielle. L’application de l’approche statistique globale est illustrée sur deux études différentes. La première concerne la détection et la caractérisation, par courants de Foucault, des fissures dans les structures conductrices. La deuxième application concerne le diagnostic des défauts de roulements dans les machines électriques tournantes. En outre, ce travail traite le problème d’estimation de l’amplitude des défauts naissants. Une analyse théorique menée dans le cadre d’une modélisation par analyse en composantes principales, conduit à un modèle analytique de la divergence ne dépendant que des paramètres du défaut. / This phD dissertation deals with the detection and diagnosis of incipient faults in engineering and industrial systems by non-parametric statistical approaches. An incipient fault is supposed to provoke an abnormal change in the measurements of the system variables. However, this change is imperceptible and also unpredictable due to the large signal-to-fault ratio and the low fault-to-noise ratio characterizing the incipient fault. The detection and identification of a global change require a ’global’ approach that takes into account the total faults signature. In this context, the Kullback-Leibler divergence is considered to be a ’global’ fault indicator, which is recommended sensitive to abnormal small variations hidden in noise. A ’global’ spectral analysis approach is also proposed for the diagnosis of faults with a frequency signature. The ’global’ statistical approach is proved on two application studies. The first one concerns the detection and characterization of minor cracks in conductive structures. The second application concerns the diagnosis of bearing faults in electrical rotating machines. In addition, the fault estimation problem is addressed in this work. A theoretical study is conducted to obtain an analytical model of the KL divergence, from which an estimate of the amplitude of the incipient fault is derived.