Le Contrôle Non Destructif (CND) de matériaux composites multicouches pour des problèmes de qualité, viabilité, sécurité et disponibilité des systèmes qui impliquent des pièces fabriquées dans les industries aéronautiques et de l’automobile est devenu une tâche essentielle aujourd’hui. L'objectif visé par cette thèse est l’imagerie électromagnétique de structures complexes multicouches anisotropes, de plus en plus utilisées dans des applications, et encore source de sérieux défis à l'étape de leur modélisation et encore plus à l'étape souvent en enfance de leur imagerie. En utilisant une vaste gamme de fréquences, qui va des courants de Foucault jusqu’aux micro-ondes, il y a un fort besoin de rendre disponibles des procédures de modélisation et d'imagerie qui sont robustes, rapides, précises et utiles à la décision des utilisateurs finaux sur des défauts potentiels, tant donc en basse fréquence (BF) (matériaux conducteurs, type fibre de carbone) qu’en haute fréquence (HF) (matériaux diélectriques, type fibre de verre). De plus, il est important d'obtenir des résultats en des temps brefs. Cependant, cela nécessite la connaissance d’une réponse précise à des sources externes aux multicouches, en considérant les couches des composites comme non endommagées ou endommagées : on parle donc de solution du problème direct, avec le cas particulier de sources élémentaires conduisant aux dyades de Green (DGF). La modélisation et la simulation numérique du problème direct sont gérés principalement via une solution au premier ordre de la formulation intégrale de contraste de source impliquant le tenseur de dépolarisation des défauts, quand ceux-ci sont assez petits vis-à-vis de l’épaisseur de peau locale (cas BF) ou de la longueur d'onde locale (cas HF). La précision des DGF doit nécessairement être assurée alors, même si les sources se situent loin de l'origine, ce qui donne un spectre de dyades qui oscille très rapidement. La technique d'interpolation-intégration dite de Padua-Domínguez est ainsi introduite dans le but d'évaluer de façon efficace des intégrales fortement oscillantes.Néanmoins, les matériaux composites peuvent souffrir de divers défauts, lors du processus de fabrication ou pendant leurs utilisations. Vides d’air, cavités remplies de liquide, fissures, etc., peuvent affecter le fonctionnement correct des structures composites. Il est donc indispensable de pouvoir détecter la présence des défauts. Ici, l’insistance est sur la méthode bien connue d’imagerie dite MUltiple SIgnal Classification (MUSIC), qui est basée sur la décomposition en valeurs singulières (SVD) des DGF ; celle-ci est développée afin de localiser les positions de multiples petits défauts volumiques en interaction faible enfouis dans des milieux anisotropes uniaxiaux. Le principal inconvénient de la méthode MUSIC est cependant sa sensibilité par rapport au bruit. Par conséquent, des méthodes MUSIC avec une résolution améliorée et la Recursively Applied and Projected (RAP) MUSIC sont introduites afin de surmonter un tel inconvénient de l'algorithme standard et de fournir des résultats de qualité avec une meilleure résolution. De nombreuses simulations numériques illustrent ces investigations. / Non-Destructive Testing/Evaluation (NdT/E) of multi-layered composite materials for problems of quality, viability, safety and availability of systems involving manufactured parts (in aeronautics and in automotive industry, as a good example) has become an interesting and challenging task nowadays. The focus of the PhD thesis is on the electromagnetic (EM) imaging of complex anisotropic multi-slab composite panels as increasingly encountered in applications, yet source of strong challenges at modeling stage and even more at often-in-infancy imaging stage. From eddy-currents to microwaves, there is a strong need to make available modeling and imaging procedures that are robust, fast, accurate and useful to potential end-users’ decision about potential defects both at low-frequency (LF) (conductive materials, carbon-fiber like) and high-frequency (HF) (dielectric materials, glass-fiber like). Moreover, it is important to get the results in close-to-real-time. However, this requires an accurate response to external sources of the multilayers, considering the layers which these composite structures are made of as undamaged or damaged. The modeling at forward stage is managed via a first-order solution involving the dyadic Green’s functions (DGF) of the layers along with the depolarization tensor of the assumed defects when they are small enough vis-à-vis the skin depth (LF case) or the wavelength (HF case). The accuracy of the DGF has to be ensured even if the sources lie far away from the origin, which yields a fast-oscillating spectrum of the dyads. The Padua-Domínguez interpolation-integration technique is introduced herein in order to evaluate in an effective fashion fast-oscillating integrals.Damages or disorders, which these composite structures may suffer from, are of many kinds. One could mention voids, fluid-filled cavities or uniaxial defects with obvious impacts on the electromagnetic and geometric parameters of the multilayers. That is, the task to make available to end-users imaging algorithms tailored to detect the presence of defects. The well-known standard MUltiple SIgnal Classification (MUSIC) algorithm, which is based on the Singular Value Decomposition (SVD) of such DGF, is here applied to localize the positions of small multiple defects with weak interaction embedded in anisotropic uniaxial media. The main drawback of MUSIC is its sensitivity with respect to the noise. Therefore, MUSIC with enhanced resolution and Recursively Applied and Projected (RAP) MUSIC are introduced to overcome such a drawback of the standard algorithm and to provide quality results with better resolution.

Les radars millimétriques en bande W (75-110 GHz) sont en plein essor, grâce notamment aux progrès des circuits intégrés, permettant de réaliser des systèmes compacts à bas coût et haute résolution due à la courte longueur d’onde. Dans un premier temps, ces systèmes ont été utilisés à des fins de détection et de localisation, avec à terme, pour objectif l’identification. Ainsi, des systèmes d’imagerie radar ont été développés, notamment grâce à l’imagerie qualitative, basée par exemple sur l’imagerie radar par synthèse d’ouverture (SAR). Cependant, afin de reconstruire les propriétés électromagnétiques des objets pour une identification complète, il est nécessaire de développer des algorithmes de reconstruction quantitatifs. Le travail présenté dans ce manuscrit est de poser les bases d’un système d’imagerie qualitative et quantitative en gamme millimétrique pour la détection et l’identification des objets sur les pistes d’aéroport par tomographie, tenant compte de la polarisation de l’onde incidente. Au cours de cette thèse, un outil de simulation permettant de la résolution des problèmes direct et inverse, pour les deux types de polarisation à deux dimensions 2D-TE et 2D-TM, basé sur la méthode des moments (MoM) a été développé. La première étape a consisté en la validation du problème direct en effectuant des comparaisons numériques avec des solutions analytiques pour des cibles canoniques. Ensuite, des mesures expérimentales ont été effectuées et comparées aux résultats numériques. Enfin, les résultats des reconstructions obtenus ont permis de valider l’algorithme de reconstruction 2D développé pour l’imagerie quantitative. / Millimeter-wave radar systems in W-band (75-110 GHz) are booming, due to advances in integrated circuits, allowing the fabrication of low-cost and high-resolution compact systems, thanks to the short wavelength. First, these systems were used for detecting and localizing purposes, with the aim of identification. Thus, imaging radar systems have been developed, especially using qualitative imaging, based for example, on Synthetic Aperture Radar (SAR). Nevertheless, in order to reconstruct the electromagnetic properties of objects, for a complete identification, we must develop quantitative reconstruction algorithms. The work presented in this manuscript is to give the basis of a qualitative and quantitative millimeter wave imaging system for detecting and identifying foreign debris on airport runways using tomography, taking into account the polarization of the incident wave. In this thesis, a simulation tool for solving forward and inverse problems, for the two-dimensional polarization cases 2D-TM and 2D-TE, based on the method of moments (MoM) has been developed. The first step was to study the validation of the direct problem by comparing numerical results with analytical solutions for canonical targets. Then, experimental measurements have ben carried out and compared with numerical results. Finally, reconstruction results obtained have validated the reconstruction algorithm developed for quantitative imaging.

The capability to detect, track and monitor human activities behind building walls and other non-line-of-sight environments is an important component of security and surveillance operations. Over the years, both ultrawideband and Doppler based radar techniques have been researched and developed for tracking humans behind walls. In particular, Doppler radars capture some interesting features of the human radar returns called microDopplers that arise from the dynamic movements of the different body parts. All the current research efforts have focused on building hardware sensors with very specific capabilities. This dissertation focuses on developing a physics based Doppler radar simulator to generate the dynamic signatures of complex human motions in nonline-of-sight environments. The simulation model incorporates dynamic human motion, electromagnetic scattering mechanisms, channel propagation effects and radar sensor parameters. Detailed, feature-by-feature analyses of the resulting radar signatures are carried out to enhance our fundamental understanding of human sensing using radar. First, a methodology for simulating the radar returns from complex human motions in free space is presented. For this purpose, computer animation data from motion capture technologies are exploited to describe the human movements. Next, a fast, simple, primitive-based electromagnetic model is used to simulate the human body. The microDopplers of several human motions such as walking, running, crawling and jumping are generated by integrating the animation models of humans with the electromagnetic model of the human body. Next, a methodology for generating the microDoppler radar signatures of humans moving behind walls is presented. This involves combining wall propagation functions derived from the finite-difference time-domain (FDTD) simulation with the free space radar simulations of humans. The resulting hybrid simulator of the human and wall is used to investigate the effects of both homogeneous and inhomogeneous walls on human microDopplers. The results are further corroborated by basic point-scatterer analysis of different wall effects. The wall studies are followed by an analysis of the effects of flat grounds on human radar signatures. The ground effect is modeled using the method of images and a ground reflection coefficient. A suitable Doppler radar testbed is developed in the laboratory for simulation validation. Measured data of different human activities are collected in both line-of-sight and through-wall environments and the resulting microDoppler signatures are compared with the simulation results. The human microDopplers are best observed in the joint timefrequency space. Hence, suitable joint time-frequency transforms are investigated for improving the display and the readability of both simulated and measured spectrograms. Finally, two new Doppler radar paradigms are considered. First, a scenario is considered where multiple, spatially distributed Doppler radars are used to measure the microDopplers of a moving human from different viewing angles. The possibility of using these microDoppler data for estimating the positions of different point scatterers on the human body is investigated. Second, a scenario is considered where multiple Doppler radars are collocated in a two-dimensional (2-D) array configuration. The possibility of generating frontal images of human movements using joint Doppler and 2-D spatial beamforming is considered. The performance of this concept is compared with that of conventional 2-D array processing without Doppler processing. / text

No description available.

Cette thèse a consisté à élaborer une méthode qui permet de résoudre l’équation intégrale comportant comme inconnues les courants et les charges introduite récemment par Taskinen et Ylä-Oijala par une méthode d’éléments frontière sans aucune contrainte de continuité au niveau des interfaces des éléments aussi bien pour les courants que pour les charges. Nous avons d’abord montré comment on pouvait construire cette équation de façon simple et similaire à celle des formulations intégrales usuelles en imposant au problème intérieur relatif au système de Picard, qui est en fait une extension du système de Maxwell, des conditions aux limites adéquates. Pour des géométries régulières de l’objet diffractant, nous avons établi de façon théorique la stabilité et la convergence des schémas numériques ci-dessus en montrant que cette équation peut être décomposée sous la forme d’un système elliptique coercif et d’un opérateur compact dans le cadre des fonctions de carré intégrable.Toute cette étude a été confirmée par des tests numériques tridimensionnels. Comme pour les équations intégrales usuelles de seconde espèce, le cadre théorique valable pour des surfaces régulières ne l’est plus pour des surfaces avec des singularités. L’utilisation formelle de cette équation,pour des surfaces singulières, a donné des résultats entachés d’erreur. Nous avons mis en évidence l’origine des instabilités numériques à l’origine de ces erreurs lorsque les géométries sont singulières en développant une version bidimensionnelle de cette équation. Cette version nous a permis en particulier de montrer que les instabilités étaient dues à des oscillations parasites concentrées autour des singularités de la géométrie. Dans ce cadre nous avons pu mettre en oeuvre plus aisément des approches pour supprimer ou atténuer ces oscillations parasites ou leur effet sur les calculs en champ lointain. Nous avons montré qu’un procédé d’augmentation des degrés de liberté pour la charge par rapport au courant pouvait sensiblement réduire ces instabilités. A la suite de l’amélioration observée sur les résultats dans le cas 2D, nous avons transposé cette procédure au cas tridimensionnel. A travers divers tests, nous avons constaté l’amélioration de la qualité de l’approximation amenée par la procédure de stabilisation / The objective of this thesis was to develop a method that solves the integral equation whose unknowns are the currents and the charges, recently introduced by Taskinen and Ylä-Oijala, by a boundary element method without any continuity constraint at the interfaces of the elements,for both the unknowns. We first show how to construct this equation in a simple way, similar tothe usual integral formulations, through imposing to the internal problem related to the Picard system,which is an extension of the Maxwell system, appropriate boundary conditions. For regular geometries, we have established a theoretical background ensuring the stability and the convergence of numerical scheme, by proving that this equation can be decomposed in a coercive elliptic and a compact parts in the context of square integrable functions. Our study was validated by three-dimensional numerical tests. In the case of usual integral equations of the second kind, the theoretical background for smooth surfaces is no longer valid when the surfaces is singular. The formal use of this equation for singular surfaces gave erroneous results. We pointed out the origin of numerical instabilities bydeveloping a two-dimensional version of this equation. This version has allowed us to show that the instabilities were due to parasitic oscillations accumulating on the geometrical singularities. In this context, we have implemented some approaches to reduce this parasitic oscillations on the calculations in the far field.We have shown that the method of increasing the freedom degrees for the charges relatively to the current could significantly reduces these instabilities. As a result, we have implemented this procedure in three-dimensional case. Throughout various tests, we noted the improvement on the approximation brough bay to the stabilization procedure

No description available.

O espalhamento eletromagnético por uma esfera com propriedades ópticas e raio arbitrários, conhecido como espalhamento de Lorenz-Mie, ou por um cilindro circular infinito, pode ser resolvido analiticamente e é comumente tratado dentro da abordagem de espalhadores dielétricos. Na região óptica, tanto meio circundante quanto partícula espalhadora possuem o mesmo valor de permeabilidade magnética. A ausência do magnetismo nessa região do espectro torna o índice de refração relativo entre os meios interno e externo ao espalhador homogêneo equivalente ao respectivo índice de impedância óptica. Em regiões espectrais de micro-ondas ou radiofrequências, entretanto, materiais ferro- e ferrimagnéticos podem exibir valores absolutos de permeabilidade magnética extremamente elevados, reduzindo então a impedância óptica em comparação ao valor correspondente de índice de refração relativo. Uma característica marcante vinculada ao magnetismo no centro espalhador é que pequenas partículas comparadas com o comprimento de onda (parâmetros de tamanho na região de Rayleigh) podem apresentar grandes seções de choque de extinção a despeito de suas pequenas seções de choque geométricas. Isso torna possível, fisicamente, a presença de picos de ressonância morfológica na energia eletromagnética interna ao centro espalhador mesmo na região de parâmetros de tamanho inferiores à unidade. Em especial, mostramos que essa energia eletromagnética possui, no regime de fraca absorção, uma relação funcional simples com o comprimento de onda incidente e a seção de choque de absorção da partícula espalhadora, independentemente do formato geométrico da mesma. No espalhamento por uma coleção de partículas magnéticas, a velocidade de transporte de energia pode ser estimada a partir da energia eletromagnética média que é armazenada no interior de um centro espalhador isolado. Dessa maneira, a validade da relação universal que encontramos entre o fator de aumento da energia eletromagnética interna ao centro espalhador e sua correspondente seção de choque de absorção (no regime de baixa absorção óptica) implica na possibilidade de estimarmos de maneira simples a velocidade de transporte de energia em um meio desordenado. Um resultado decorrente dessa aproximação é que mesmo na região de Rayleigh a velocidade de transporte de energia em um meio magnético desordenado é drasticamente reduzida, levando, por conseguinte, à redução do coeficiente de difusividade dos fótons no meio. O estudo analítico e numérico da energia eletromagnética armazenada por uma esfera e por um cilindro magnéticos irradiados por ondas planas homogêneas são os temas em foco nesta dissertação de Mestrado. / Electromagnetic scattering by a sphere with arbitrary optical properties and radius, known as the Lorenz-Mie scattering, or by an infinite right circular cylinder can be solved analytically and is widely treated in the approach of dielectric scatterers. In the optical range, both embedding medium and scattering particle have the same magnetic permeability. The absence of magnetism in this spectral range leads to the equivalence between the relative refraction and impedance indices associated with the scatterers. However, in microwave or radio-frequency ranges, ferro- and ferrimagnetic materials can exhibit extremely huge values of magnetic permeability, which reduce the optical impedance in comparison to the corresponding value of relative refraction index. One striking feature associated with the magnetism in the scatterer is that particles smaller than the wavelength (Rayleigh size region) can present large extinction cross sections in despite of their small geometric cross sections. This becomes physically possible the presence of morphology-dependent resonances in the electromagnetic energy within the scatterer even in size parameters region smaller than unity. In particular, we show that this time-averaged electromagnetic energy has, in the weak absorption regime, a simple functional relation with the incident wavelength and the scatterer absorption cross section which does not depend on the shape of the scatterer. In the multiple scattering regime, the energy-transport velocity can be estimated from the time-averaged electromagnetic energy stored in a single scatterer. Thereby, the validity of the universal relation between the internal energy-enhancement factor and the absorption cross section respective to an arbitrary scattering center (in the weak absorption regime) implies that the energy-transport velocity in disordered media can be evaluated in a simple way. From this approximation, we obtain that even in the Rayleigh size region the energy-transport velocity in disordered magnetic media is dramatically reduced, which consequently leads to a reduction of the diffusion coefficient of the photons. The analytical and numerical studies of the time-averaged electromagnetic energy within magnetic isotropic spheres and cylinders irradiated by plane waves are our aim in this Master\'s degree dissertation.

O espalhamento eletromagnético por uma esfera com propriedades ópticas e raio arbitrários, conhecido como espalhamento de Lorenz-Mie, ou por um cilindro circular infinito, pode ser resolvido analiticamente e é comumente tratado dentro da abordagem de espalhadores dielétricos. Na região óptica, tanto meio circundante quanto partícula espalhadora possuem o mesmo valor de permeabilidade magnética. A ausência do magnetismo nessa região do espectro torna o índice de refração relativo entre os meios interno e externo ao espalhador homogêneo equivalente ao respectivo índice de impedância óptica. Em regiões espectrais de micro-ondas ou radiofrequências, entretanto, materiais ferro- e ferrimagnéticos podem exibir valores absolutos de permeabilidade magnética extremamente elevados, reduzindo então a impedância óptica em comparação ao valor correspondente de índice de refração relativo. Uma característica marcante vinculada ao magnetismo no centro espalhador é que pequenas partículas comparadas com o comprimento de onda (parâmetros de tamanho na região de Rayleigh) podem apresentar grandes seções de choque de extinção a despeito de suas pequenas seções de choque geométricas. Isso torna possível, fisicamente, a presença de picos de ressonância morfológica na energia eletromagnética interna ao centro espalhador mesmo na região de parâmetros de tamanho inferiores à unidade. Em especial, mostramos que essa energia eletromagnética possui, no regime de fraca absorção, uma relação funcional simples com o comprimento de onda incidente e a seção de choque de absorção da partícula espalhadora, independentemente do formato geométrico da mesma. No espalhamento por uma coleção de partículas magnéticas, a velocidade de transporte de energia pode ser estimada a partir da energia eletromagnética média que é armazenada no interior de um centro espalhador isolado. Dessa maneira, a validade da relação universal que encontramos entre o fator de aumento da energia eletromagnética interna ao centro espalhador e sua correspondente seção de choque de absorção (no regime de baixa absorção óptica) implica na possibilidade de estimarmos de maneira simples a velocidade de transporte de energia em um meio desordenado. Um resultado decorrente dessa aproximação é que mesmo na região de Rayleigh a velocidade de transporte de energia em um meio magnético desordenado é drasticamente reduzida, levando, por conseguinte, à redução do coeficiente de difusividade dos fótons no meio. O estudo analítico e numérico da energia eletromagnética armazenada por uma esfera e por um cilindro magnéticos irradiados por ondas planas homogêneas são os temas em foco nesta dissertação de Mestrado. / Electromagnetic scattering by a sphere with arbitrary optical properties and radius, known as the Lorenz-Mie scattering, or by an infinite right circular cylinder can be solved analytically and is widely treated in the approach of dielectric scatterers. In the optical range, both embedding medium and scattering particle have the same magnetic permeability. The absence of magnetism in this spectral range leads to the equivalence between the relative refraction and impedance indices associated with the scatterers. However, in microwave or radio-frequency ranges, ferro- and ferrimagnetic materials can exhibit extremely huge values of magnetic permeability, which reduce the optical impedance in comparison to the corresponding value of relative refraction index. One striking feature associated with the magnetism in the scatterer is that particles smaller than the wavelength (Rayleigh size region) can present large extinction cross sections in despite of their small geometric cross sections. This becomes physically possible the presence of morphology-dependent resonances in the electromagnetic energy within the scatterer even in size parameters region smaller than unity. In particular, we show that this time-averaged electromagnetic energy has, in the weak absorption regime, a simple functional relation with the incident wavelength and the scatterer absorption cross section which does not depend on the shape of the scatterer. In the multiple scattering regime, the energy-transport velocity can be estimated from the time-averaged electromagnetic energy stored in a single scatterer. Thereby, the validity of the universal relation between the internal energy-enhancement factor and the absorption cross section respective to an arbitrary scattering center (in the weak absorption regime) implies that the energy-transport velocity in disordered media can be evaluated in a simple way. From this approximation, we obtain that even in the Rayleigh size region the energy-transport velocity in disordered magnetic media is dramatically reduced, which consequently leads to a reduction of the diffusion coefficient of the photons. The analytical and numerical studies of the time-averaged electromagnetic energy within magnetic isotropic spheres and cylinders irradiated by plane waves are our aim in this Master\'s degree dissertation.

Les instruments, satellites ou systèmes aéroportés, actuellement utilisés pour la détection et la caractérisation d'hydrocarbure sur la mer sont basés sur des moyens optiques ou radars. Ces moyens présentent une performance dégradée due à une fréquence encore trop importante de fausses alarmes ou à un temps de traitement des données trop conséquent. Les méthodes de détection, d'identification et de quantification des fuites d'hydrocarbures offshores peuvent donc être améliorées en associant robustesse et réactivité. Cette amélioration suppose une compréhension approfondie des phénomènes océanographiques et électromagnétiques à l'œuvre dans cette scène particulière. La thèse s'appuie sur des données regroupant des images optiques et SAR aéroportées ou satellites ainsi que des mesures réalisées en laboratoire. Ce jeu de données permet de vérifier la cohérence des résultats obtenus par modélisation. L'objectif de la thèse est de distinguer une surface de mer polluée d'une surface de mer propre à l'aide de la signature électromagnétique de la surface totale puis de détailler le type et la quantité d'hydrocarbure présent. La thèse se divise en deux domaines, à savoir modélisation océanographique et modélisation électromagnétique. La modélisation océanographique intègre la simulation de la surface rugueuse imitant une surface de mer propre, et polluée. Cette surface de mer doit être générée sur une superficie importante et doit conserver une résolution restituant les petites vagues avec un temps de génération minimal. La partie électromagnétique est centrée sur les modèles asymptotiques de diffusion des ondes électromagnétiques par une interface rugueuse. Ces modèles sont adaptés au contexte de la thèse, complexité de la scène et rapidité du traitement, mais nécessitent plusieurs hypothèses pour être appliqués. / Satellites or airborne systems currently used for the detection and characterization of oil slicks on sea surface are based on optical or radar means. These means have a lack of performance due to a too high frequency of false alarms or to an excessively long data processing time. The methods for detecting, identifying and quantifying offshore pollutant can therefore be improved by combining robustness and reactivity. This improvement implies an in-depth understanding of the oceanographic and electromagnetic phenomena at work in this particular scene. The thesis is based on data gathering aerial and satellite images and SAR as well as measurements carried out in laboratory. This dataset makes it possible to check the consistency of the results obtained by modeling. The objective of the thesis is to distinguish a polluted sea surface from a clean sea surface using the electromagnetic signature of the total surface and then to detail the type and quantity of pollutant. The thesis is divided into two domains, namely oceanographic modeling and electromagnetic modeling. Oceanographic modeling integrates the simulation of the rough surface imitating a clean or polluted sea surface. This sea surface must be generated over a large area with a thin resolution. The electromagnetic part is centered on the asymptotic models for the electromagnetic waves diffraction by a rough interface. These models are adapted to the context of the thesis, the complexity of the scene and the speed of processing, but require several hypotheses to be applied.

No description available.