L’île de La Réunion a été couverte par une vaste campagne de géophysique héliportée en 2014. Le dispositif SkyTEM a permis de cartographier l’anomalie du champ magnétique et la résistivité électrique du sous-sol jusqu’à 300 m de profondeur, et ce, sur l’ensemble de l’île. L’acquisition de ce jeu de données met en évidence la problématique de l’interprétation des données électromagnétiques héliportées à l’échelle régionale. L’interprétation d’un modèle 3D de résistivité d’une superficie 2 500 km2 nécessite la mise en place de méthodologies innovantes pour le confronter aux données géologiques, climatiques et hydrogéologiques. L’objectif de la thèse est de valoriser ce jeu de données géophysique à l’échelle régionale en intégrant l’analyse hydrogéologique à différentes échelles spatiales. Pour répondre à ces problématiques, cette thèse s’organise autour d’une méthodologie hydrogéophysique multi-échelle qui permet d’étudier la répartition et le fonctionnement des aquifères à des échelles inédites jusqu’alors. Pour interpréter les 350 000 sondages électromagnétiques de la campagne, une approche statistique a été développée pour synthétiser une information sur les contrastes 3D en 2D. Les résultats, cohérents avec les connaissances géologiques et hydrogéologiques de l’île, ont permis d’analyser l’impact des paramètres géologiques et climatiques régionaux sur la répartition des aquifères et les processus d’altération. Sur le littoral, les données de résistivité sont confrontées à des logs de conductivité électrique de l’eau souterraine afin de valider l’interprétation de la géophysique. Le modèle de résistivité permet ainsi de cartographier et d’analyser les paramètres de contrôle du phénomène d’intrusion saline à l’échelle de l’île. En altitude, l’utilisation conjointe des données électromagnétiques et magnétiques héliportées permet de caractériser la géométrie des structures géologiques et d’imager les profils et paléo-profils d’altération météorique. / In 2014, Reunion Island was covered by a vast airborne geophysics campaign. The SkyTEM device mapped magnetic anomalies and bulk resistivity at an extra 300 m depth throughout the island. The resulting dataset revealed the issue of interpreting airborne electromagnetic data at the regional scale. The interpretation of a 3D resistivity model wrapping 2 500 km2 requires the implementation of innovative methodologies to compare it with geological, climatic and hydrogeological data. The aim of the thesis is to enhance this geophysical dataset at the regional scale by integrating hydrogeological analysis at different spatial scales. To answer these issues, this thesis revolves around a multi-scale hydrogeophysical methodology, which allows to analyze how aquifers are distributed and their behaviors on unprecedented scales. A statistical approach has been developed to synthesize 3D geo-electrical contrasts in 2D to interpret the 350 000 soundings. All outcomes are consistent with the geological and hydrogeological knowledge of the island, enabling us to study how regional geological and climatic parameters impact the distribution of aquifers and weathering processes. In the coastal zone, resistivity data are compared with groundwater electrical conductivity logs to validate the interpretation of geophysics. Thus, the resistivity model helped mapping saltwater intrusion and analyzing its control parameters at the island scale. Inland, the combined use of airborne electromagnetic and magnetic data made it possible to characterize the geometry of geological structures and to image weathering profiles and paleo-profiles. Combining results gathered on various scales allows the analysis of weathering impact on aquifers of shield volcanoes. On the coastal zone, un-weathered and highly permeable lithological units (up to 10-1 m/s) are subject to saltwater intrusion, even with recharge rates up to 6 to 8 m per year. With decreasing transmissivity, aquifers become less exposed to saltwater intrusion. According to geophysics outcomes, the contacts between different eruptive phase products at high altitude are characterized by specific weathering paleo-profiles which create permeability contrasts conducive to perched aquifers establishment. Examining the geophysical response in different areas of the island made it possible to analyze the dynamics of weathering processes, initially controlled by the age of the formations and secondly, by the rainfall rates and the temperature. These results showed that weathering and structure of basaltic volcanoes depend on their age and on windward/leeward slopes. Thus, our research helped improve the accuracy of the conceptual hydrogeological model including geological structuring and time evolution of basaltic volcanic aquifers in an insular environment.

The development of reliable and cost-saving methods to select pre-clinically new potential drugs with unknown and original mechanisms of action for cancer therapy has become crucial. Previous investigations demonstrated that infrared spectra of cancer cells exposed to well-known anticancer drugs used in the clinics provide a global fingerprint of all the metabolic modifications induced. Nowadays several natural products have been recognized for their medicinal values. Polyphenolic compounds constitute one of the largest groups of plant metabolites and many studies have demonstrated their anticancer properties at multiple steps of carcinogenesis. Taking into account the large diversity of polyphenolic structures in nature and their numerous targets against carcinogenesis, the step of selection becomes essential as it is virtually impossible to classify them using traditional classification techniques. While FTIR spectroscopy appears to have a definite potential to sort anticancer drugs on the basis of the metabolic modifications they induced, the present challenge in this thesis is to evaluate the drug-induced spectral changes in cancer cells on a larger scale. The coupling of FTIR spectroscopy with a high throughput screening extension could become a useful method to generate drug classifications based on the “modes of action”. In a first step, the IC50 was evaluated for each polyphenol to normalize all spectral experiments that will be carried out on breast cancer cells. The first experiments revealed dispersive artifacts (wide variations of absorbance distribution and Mie scattering effects) contributing to spectral measurements. To minimize them, the best option was covering uniformly the surface of spectral measurement with cells (detailed in chapter III). Once the protocol refined, it was applied to the study of MDA-MB-231 breast cancer cells exposed for 24 hrs to 15 polyphenols. Through unsupervised and supervised statistical analyses, a distinction between polyphenol treatments could be well established. Complex effects of polyphenols on cancer cells were revealed, suggesting that mechanisms specific to each polyphenol were evidenced by the whole infrared spectrum. Clustering of polyphenol-induced spectral signatures by Hierarchical Cluster Analysis indicates that some of the polyphenols share similar effects on MDA-MB-231 (detailed in chapter IV). This experiment was then extended to another breast cancer cell line, MCF-7. It demonstrated a cell line dependency to polyphenolic treatments. Finally a subcellular investigation of treated MCF-7 breast cancer cells in their live state was done using Raman imaging. A distinction between nucleus and cytoplasm of treated cells brought supplementary information regarding the effect of polyphenols, leading to subcellular biological assumptions on polyphenol effects. These results paves the way for a new classification based on infrared spectral signatures that reflect the overall chemical modifications experienced by the cells when exposed to drugs. / Le développement de méthodes fiables et peu cher pour sélectionner de potentiels nouveaux médicaments présentant un mécanisme d’action original et inconnu avant toute étape clinique devient crucial. De précédentes études ont pu démontrer que les spectres infrarouges de cellules cancéreuses exposées à des agents anticancéreux connus et utilisé dans le monde clinique fournissent une empreinte globale de toutes les modifications métaboliques induites. La spectroscopie infrarouge est un outil innovant qui semble prometteur pour offrir un aperçu global des processus biologiques et physiologiques qui sont menés par un médicament dans des cellules cancéreuses. De nos jours, de nombreux produits naturels ont été reconnus pour leurs propriétés médicinales. Les polyphénols constituent l’un des plus vastes groupes de métabolites végétaux et de nombreuses études ont démontré leurs propriétés anticancéreuses à de multiples étapes de la carcinogénèse. En prenant en compte la très grande diversité de structures polyphénoliques existantes dans la nature et leurs nombreuses cibles anti-tumorales, l’étape de sélection est devenue essentielle comme il est virtuellement impossible de les classifier grâce à des techniques de classification traditionnelles telles que les études –omiques. Dès lors, le défi de cette thèse est d’évaluer les variations spectrales induites par un polyphénol dans des cellules cancéreuses à une plus grande échelle. Le couplage de la spectroscopie IRTF avec une extension de criblage de haut débit pourrait devenir une méthode utile pour générer des classifications de molécules sur base de leur « modes d’action ». Dans un premier temps, la concentration qui inhibe 50% de la croissance des cellules cancéreuses fut déterminée pour chaque polyphénol et chaque lignée de cellules cancéreuses. Le traitement des cellules à cette concentration permet une normalisation interne des expériences réalisées ultérieurement en spectroscopie infrarouge. Une fois le protocole établi, la lignée MDA-MB-231 fut exposée durant 24 heures à 15 polyphénols différents. Au moyen d’analyses statistiques multivariées supervisées et non supervisées, une distinction parmi les polyphénols a pu être établie et des effets complexes des polyphénols sur les cellules cancéreuses ont pu être révélés, suggérant des mécanismes d’action spécifiques à chaque polyphénol mis en évidence par spectroscopie infrarouge. Finalement, une étude subcellulaire sur cellules vivantes fut réalisée par imagerie Raman sur une seconde lignée de cellules cancéreuses mammaires appelées MCF-7. Cela permis de compléter en partie l’information macroscopique offerte par la spectroscopie infrarouge par une information microscopique sur l’effet de certains polyphénols. Cette thèse a ouvert la voie pour de nouvelles techniques de classification d’agents anticancéreux basées sur la spectroscopie infrarouge, technique sensible à l’ensemble des modifications chimiques subies par des cellules. / Option Chimie du Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Cette thèse aborde la problématique de la modélisation multiphysique en génie électrique, avec une application à l’étude des vibrations d’origine électromagnétique des cages de développantes. Cette étude comporte quatre parties : la construction de la densité de courant, le calcul des forces locales, le transfert de solutions entre maillages et la résolution des problèmes couplés. Un premier enjeu est de correctement représenter les courants, cette opération est effectuée en deux étapes : la construction de la densité de courant et l’annulation de la divergence. Si des structures complexes sont utilisées, l’imposition du courant ne peut pas toujours être réalisée à l’aide de méthodes analytiques. Une méthode basée sur une résolution électrocinétique ainsi qu’une méthode purement géométrique sont testées. Cette dernière donne des résultats plus proches de la densité de courant réelle. Parmi les nombreuses méthodes de calcul de forces, les méthodes des travaux virtuels et des forces de Laplace, considérées par la littérature comme les plus adaptées au calcul des forces locales, ont été étudiées. Nos travaux ont montré que bien que les forces de Laplace sont particulièrement précises, elles ne sont pas valables si la perméabilité n’est plus homogène. Ainsi, la méthode des travaux virtuels, applicable de manière universelle, est préférée. Afin de modéliser des problèmes multi-physiques complexes à l’aide de plusieurs codes de calculs dédiés, des méthodes de transferts entre maillages non conformes ont été développées. Les procédures d’interpolations, les méthodes localement conservatives et les projections orthogonales sont comparées. Les méthodes d’interpolations sont réputées rapides mais très diffusives tandis que les méthodes de projections sont considérées comme les plus précises. La méthode localement conservative peut être vue comme produisant des résultats comparables aux méthodes de projections, mais évite l’assemblage et la résolution de systèmes linéaires. La modélisation des problèmes multi-physiques est abordée à l’aide des méthodes de transferts de solutions. Pour une classe de problème donnée, l’assemblage d’un schéma de couplage n’est pas unique. Des tests sur des cas analytiques sont réalisés afin de déterminer, pour plusieurs types de couplages, les stratégies les plus appropriées.Ces travaux ont permis une application à la modélisation magnéto-mécanique des cages de développantes est présentée. / The modeling of multi-phycics problems in electrical engineering is presented, with an application to the numerical computation of vibrations within the end windings of large turbo-generators. This study is divided into four parts: the impositions of current density, the computation of local forces, the transfer of data between disconnected meshes, and the computation of multi-physics problems using weak coupling, Firstly, the representation of current density within numerical models is presented. The process is decomposed into two stages: the construction of the initial current density, and the determination of a divergence-free field. The representation of complex geometries makes the use of analytical methods impossible. A method based on an electrokinetical problem is used and a fully geometrical method are tested. The geometrical method produces results closer to the real current density than the electrokinetical problem. Methods to compute forces are numerous, and this study focuses on the virtual work principle and the Laplace force considering the recommendations of the literature. Laplace force is highly accurate but is applicable only if the permeability is uniform. The virtual work principle is finally preferred as it appears as the most general way to compute local forces. Mesh-to-mesh data transfer methods are developed to compute multi-physics models using multiples meshes adapted to the subproblems and multiple computational software. The interpolation method, a locally conservative projection, and an orthogonal projection are compared. Interpolation method is said to be fast but highly diffusive, and the orthogonal projections are highly accurate. The locally conservative method produces results similar to the orthogonal projection but avoid the assembly of linear systems. The numerical computation of multi-physical problems using multiple meshes and projections is then presented. However for a given class of problems, there is not an unique coupling scheme possible. Analytical tests are used to determine, for different class of problems, the most accurate scheme. Finally, numerical computations applied to the structure of end-windings is presented.

Cette thèse est dédiée à la résolution numérique tridimensionnelle des équations de Maxwell harmoniques, par des méthodes de décomposition de domaine couplant des résolutions par équations intégrales entre elles. Pour traiter les problèmes de diffraction d'ondes, la méthode des équations intégrales est un outil précieux. Elle consiste à paramétrer le champ électromagnétique solution par une source définie sur la surface de l'objet diffractant, solution d'une nouvelle équation linéaire (l'équation intégrale). Pour des applications à haute fréquence, le grand nombre d'inconnues (de l'ordre du million) nous oblige à utiliser un solveur itératif pour résoudre l'équation intégrale. Le problème du conditionnement des systèmes linéaires est alors crucial. De récents développements ont permis de construire une équation intégrale performante (la GCSIE) et de conditionnement stable avec la montée en fréquence. Cependant, la présence d'une cavité large et résonnante dans l'objet diffractant (telle que la cavité moteur d'un avion) dégrade le conditionnement de cette équation. Nous proposons deux méthodes de décomposition de domaine (DDM) afin de découpler le problème de la cavité du problème extérieur. La première (DDM en Y) s'exprime en fonction des opérateurs Dirichlet-to-Neumann Y, qui sont synthétisés via la résolution de problèmes métalliques par équations intégrales dans chaque sous-domaine. La seconde (DDM en S) s'exprime en fonction des opérateurs de scattering S, synthétisés par résolution de problèmes de type métal-impédant, donc bien posés à toute fréquence. La DDM en S permet ainsi de se débarrasser des phénomènes de résonance dans les cavités. Nous proposons dans un premier temps un préconditionneur analytique pour la DDM en Y, basé sur l'opérateur électromagnétique de simple couche. Nous calculons ensuite les modes guidés le long d'un cylindre infini tangent à la cavité près de l'interface, et nous diagonalisons les opérateurs Dirichlet-to-Neumann et scattering dans la base des traces de modes guidés sur l'interface. On extrait de cette étude deux préconditionneurs spectraux respectivement pour la DDM en Y et la DDM en S. Les résultats numériques confirment l'efficacité des préconditionneurs proposés / This work is dedicated to the numerical solution of the tridimensional harmonic Maxwell equations, using domain decomposition methods coupling integral equations between them. To deal with scattering problems, integral equations methods are a precious tool. They allow to look for the electromagnetic field by parameterizing it with a source only defined on the boundary of the scattering object, solution of a new linear equation (the integral equation). For applications at high frequency, the great number of unknowns forces the use of iterative methods. To accelerate the solution of integral equations, one moreover has to ensure the good condition number of the linear systems, or to propose well-suited preconditioners. An efficient method, the GCSIE, was developed in Onera. It is an intrinsically well-conditioned integral equation whose condition number remains stable whith the frequency increase. However, the existence of large and resonant cavities (such as air intakes) deteriorates the condition number. In order to circumvent this problem, we propose two domain decomposition methods (DDM) allowing to decouple the exterior problem from the problem of the cavities. The first one (Y-DDM) is based on Dirichlet-to-Neumann operators Y, which are built through the solution of metallic problems using integral equations in each subdomain. The second one (S-DDM) is based on scattering operators S, built through the solution of problems of metallic-impedant type, which are well-posed at any frequency. The S-DDM allows to avoid the resonance phenomena inside the cavities. First, we propose an analytic preconditioner for the Y-DDM, based on the electromagnetic single layer operator. We then calculate the modes guided along an artificial infinite cylinder, that is tangent to the cavity near the interface. We diagonalize the Dirichlet-to-Neumann and scattering operators in the basis of the traces of the guided modes on the interface. We deduce from this study two spectral preconditioners for the Y-DDM and the S-DDM. The numerical results confirm the efficiency of the employed preconditioners.

La méthode des éléments frontières pour l’équation des ondes (BEM) est utilisée en acoustique eten électromagnétisme pour simuler la propagation d’une onde avec une discrétisation en temps(TD). Elle permet d’obtenir un résultat pour plusieurs fréquences à partir d’une seule résolution.Dans cette thèse, nous nous intéressons à l’implémentation efficace d’un simulateur TD-BEM sousdifférents angles. Nous décrivons le contexte de notre étude et la formulation utilisée qui s’exprimesous la forme d’un système linéaire composé de plusieurs matrices d’interactions/convolutions.Ce système est naturellement calculé en utilisant l’opérateur matrice/vecteur creux (SpMV). Nousavons travaillé sur la limite du SpMV en étudiant la permutation des matrices et le comportementde notre implémentation aidé par la vectorisation sur CPU et avec une approche par bloc surGPU. Nous montrons que cet opérateur n’est pas approprié pour notre problème et nous proposonsde changer l’ordre de calcul afin d’obtenir une matrice avec une structure particulière.Cette nouvelle structure est appelée une matrice tranche et se calcule à l’aide d’un opérateur spécifique.Nous décrivons des implémentations optimisées sur architectures modernes du calculhaute-performance. Le simulateur résultant est parallélisé avec une approche hybride (mémoirespartagées/distribuées) sur des noeuds hétérogènes, et se base sur une nouvelle heuristique pouréquilibrer le travail entre les processeurs. Cette approche matricielle a une complexité quadratiquesi bien que nous avons étudié son accélération par la méthode des multipoles rapides (FMM). Nousavons tout d’abord travaillé sur la parallélisation de l’algorithme de la FMM en utilisant différentsparadigmes et nous montrons comment les moteurs d’exécution sont adaptés pour relâcher le potentielde la FMM. Enfin, nous présentons des résultats préliminaires d’un simulateur TD-BEMaccéléré par FMM . / The time-domain BEM for the wave equation in acoustics and electromagnetism is used to simulatethe propagation of a wave with a discretization in time. It allows to obtain several frequencydomainresults with one solve. In this thesis, we investigate the implementation of an efficientTD-BEM solver using different approaches. We describe the context of our study and the TD-BEMformulation expressed as a sparse linear system composed of multiple interaction/convolutionmatrices. This system is naturally computed using the sparse matrix-vector product (SpMV). Wework on the limits of the SpMV kernel by looking at the matrix reordering and the behavior of ourSpMV kernels using vectorization (SIMD) on CPUs and an advanced blocking-layout on NvidiaGPUs. We show that this operator is not appropriate for our problem, and we then propose toreorder the original computation to get a special matrix structure. This new structure is called aslice matrix and is computed with a custom matrix/vector product operator. We present an optimizedimplementation of this operator on CPUs and Nvidia GPUs for which we describe advancedblocking schemes. The resulting solver is parallelized with a hybrid strategy above heterogeneousnodes and relies on a new heuristic to balance the work among the processing units. Due tothe quadratic complexity of this matrix approach, we study the use of the fast multipole method(FMM) for our time-domain BEM solver. We investigate the parallelization of the general FMMalgorithm using several paradigms in both shared and distributed memory, and we explain howmodern runtime systems are well-suited to express the FMM computation. Finally, we investigatethe implementation and the parametrization of an FMM kernel specific to our TD-BEM, and weprovide preliminary results.

Les procédés de traitement de surface sont utilisés à l'échelle industrielle pour améliorer les performances de pièces mécaniques en introduisant des contraintes résiduelles de compression. Cette mise en compression de surface permet de limiter l'amorçage et la propagation de fissures dans le matériau. Ceci permet d'augmenter de façon significative la durée de vie en fatigue des pièces mécaniques ainsi traitées. L'utilisation de ces procédés dans l'industrie a démontré leur efficacité, mais aussi leurs limitations et inconvénients. Les défauts récurrents consistent en une profondeur traitée faible, une dégradation de l'état de surface (rugosité), des difficultés de contrôle, une contamination du matériau traité, etc. Ces défauts ont conduit à l'élaboration de nouveaux procédés innovants qui permettent de meilleures performances en évitant certains des inconvénients succinctement évoqués. Parmi ces procédés innovants, le traitement de surface par impulsion électromagnétique semble particulièrement intéressant. Ce procédé met en œuvre un puissant champ magnétique transitoire pour engendrer des forces de Laplace dans une pièce métallique et induire des contraintes résiduelles. Il n'existe que peu d'informations dans la littérature et il n'existe aucun dispositif expérimental de ce procédé. Cette thèse est dédiée à la conception et la réalisation d'un prototype de mise en compression électromagnétique. Le premier chapitre de cette thèse est un état de l'art des technologies de mise en compression et du procédé de mise en compression par impulsion électromagnétique. Ainsi, les besoins de ce procédé sont identifiés et les technologies pouvant répondre à ces besoins sont explorées. Le deuxième chapitre, après une sélection de la structure globale du dispositif, va consister aux dimensionnements des éléments du prototype EMP. Cette étude commencera avec une étude sur l'inducteur qui va être utilisé avant de continuer sur le dimensionnement du stockage d'énergie et de l'interrupteur de décharge. Afin de valider le dimensionnement des composants précédents, une simulation électromagnétique 3D du système est réalisée. L'assemblage du prototype est présenté dans le troisième chapitre ainsi qu'une première campagne d'essai sur un alliage d'aluminium. Deux types d'éprouvettes sont testées : une éprouvette fine pour vérifier visuellement la mise en compression (essai Almen) et une éprouvette massive afin d'évaluer la profondeur traitée. Une modélisation multiphysique 3D du procédé est réalisée afin de corréler ces résultats avec l'expérience. Dans un dernier chapitre, une étude exploratoire est menée sur un matériau ferromagnétique, le mumétal, pour visualiser l'influence des contraintes résiduelles sur les propriétés magnétiques de ce dernier. / Penning processes are widely used in industries to apply compressive residual stresses into the most solicited part of mechanical pieces. In that way, the compressive residual stresses limit the priming and the propagation of micro-cracks in the material. This increases significantly the lifespan of the treated mechanical piece under fatigue stresses. These existing peening processes have proved their efficiency and also their limitations and weaknesses. The main recurrent defaults are a shallow depth of treatment, a degradation of the surface condition, a random control of the treatment, a material contamination, etc. These problems have led towards the development of news innovative peening processes which allow better performance avoiding some previous defaults briefly evoked. Among these news processes, the electromagnetic peening process seems especially interesting. This process uses high energy electromagnetic fields to induce Lorentz forces into a metallic piece and thus residual stresses. Actually, there is not much information about this process in the literature and no prototype was ever built. The work of this thesis is dedicated to development and realization of an electromagnetic peening prototype. The first chapter of this thesis adresses the state of the art of major peening processes actually in industrial use. Next, the electromagnetic peening process, or EMP process, is described and the electrical needs are exposed. A second state of the art is made about the technological solutions to respond to the EMP needs. The second chapter is about the conception of the EMP prototype with the electrical structure adopted in the previous chapter. The first step is about the inductor sizing to generate an electromagnetic field sufficient enough for a peening application. Next, the storage system is designed depending on the inductor parameters and finally the closing switch is created considering the electrical parameters used for the EMP process. To validate the previous results, a 3D electromagnetic simulation is done. The prototype assembly is presented in the third chapter and also the first experimental test on the EMP prototype. To begin with, an aluminium alloy with low yield strength is selected to be treated. Two different samples forms are used, a thin one, to realize a similar test to the Almen test and thick one to check the EMP depth of treatment. A 3D multiphysics simulation of these experiments is made and these numeric results are next correlated to the experimental ones. In the fourth chapter, an exploratory study is realized on the effects of the residual stresses on magnetic properties of ferromagnetic material, the mumetal

L'évolution de la technologie CMOS autorise la réduction des dimensions transverses du pixel d'imageur qui atteignent aujourd'hui 1.75μm. Son épaisseur, en revanche, ne peut être réduite sans modifier le principe de filtrage existant, assuré par des résines colorées absorbantes dont l'épaisseur fixe le profil spectral en transmission. Ces résines ont, de plus, une faible tenue en température et doivent être positionnées au sommet de l'empilement métallo-diélectrique que constitue le pixel. Ces deux facteurs contribuent à la dégradation de la résolution optique du capteur. Dans ce mémoire, nous discutons des possibilités de les remplacer par des structures diffractives métalliques bi-périodiques, moins épaisses et plus résistantes aux traitements thermiques. Pour cela, une nouvelle formulation de la Méthode des Eléments Finis a été introduite. Cette technique, très générale, de résolution d'équations aux dérivées partielles est adaptée, de par sa flexibilité, à la richesse géométrique du pixel CMOS et à la complexité des structures diffractives pressenties. Dans un premier temps, cette méthode a été validée numériquement dans le cas scalaire, puis expérimentalement grâce à des photodiodes de test munies de réseaux mono-dimensionnels. Dans un second temps, la formulation a été généralisée au cas vectoriel de la diffraction d'une onde plane d'incidence et de polarisation arbitraires par un réseau bi-dimensionnel quelconque. Enfin, ce modèle est exploité pour dégager un jeu de paramètres opto-géométriques de réseaux croisés permettant de filtrer la lumière dans ses composantes Rouge, Verte et Bleue.

Nous nous intéressons à l'imagerie d'un objet enfoui dans un milieu multi-couches inhomogène, avec des données ne contenant pas la phase. Nous résolvons un problème direct modèle de propagation des ondes dans un tel milieu, à l'aide de l'analyse asymptotique et des équations intégrales. Puis nous développons des algorithmes de reconstruction à base de dérivée topologique et des techniques de l'optimisation de forme.

I Première partie : Contribution à l'étude des fibres optiques microstructurées 1 Introduction à l'étude des fibres optiques microstructurées 2 Two numerical methods to study microstructured optical fibres 3 Main properties of microstructured optical fibres Conclusion II Deuxième partie : sélection d'articles III Troisième partie : annexes, bibliographie et index A Some details of the Multipole Method derivation B Computing band diagram using the Finite Element Method Bibliographie Index

Ce mémoire traite de la contribution des méthodes géophysiques à la détection et à la cartographie des cavités souterraines de la Médina de Béni-Mellal. Il s'inscrit dans un projet de réhabilitation de la Médina, à l'initiative de l'Agence Urbaine de Béni-Mellal et des autorités municipale, régionale et nationale. La méthodologie mise en œuvre est celle proposée dans le Guide technique élaboré par le Laboratoire Central des Ponts et Chaussées de Nantes (2004). Le premier chapitre pose le problème des cavités, dites "Khefs", en relation avec l'histoire et la géographie de la médina. Le second chapitre traite des phases d'étude préliminaires à l'application de la géophysique. Une étude critique des documents et rapports antérieurs et des visites d'affleurements géologiques et de cavités, ont modifié profondément les hypothèses sur l'origine et la nature des cavités. Contrairement à ce qui a souvent été dit et écrit, il n'y a aucune cavité d'origine naturelle liée à la karstification. Il s'agit de cavités creusées par l'homme dans des matériaux détritiques et lacustres pour en extraire les matériaux de construction de la médina elle-même, sous la forme d'exploitation classique de type "chambres et piliers". La dégradation des cavités est également liée à l'activité anthropique. Le troisième chapitre traite du choix des méthodes géophysiques mises en œuvre et présente les résultats de trois d'entre elles. La sismique réfraction n'est pas recommandée. La tomographie électrique et l'électromagnétisme en champ proche sont en revanche recommandés. Sept exemples de profils de tomographie 2D sont présentés, qui montrent clairement la présence d'anomalies résistantes attendues : leur toit est situé entre 0,50 m et 1, 50 m de profondeur ; leur extension latérale est inférieure à 5m. Cinq profils et 2 cartes électromagnétiques démontrent la capacité à cartographier les variations de la géologie superficielle et la sensibilité aux contacts verticaux entre les vides et les matériaux détritiques en place. Le quatrième chapitre présente les résultats obtenus avec le radar géologique sur une zone expérimentale de la médina. Des exemples spécifiques sont extraits des 270 profiles qui ont été enregistrés, avec une antenne de 400 MHz pour la totalité et une antenne de 200 MHz pour quelques uns. Les résultats ont été contrôlés par 17 forages. La confirmation est satisfaisante : 7 forages ont rencontré les vides tel que prévu, 2 les ont rencontrés à une profondeur supérieure à la profondeur d'investigation, 5 ont mis en évidence un réflecteur géologique comme prévu, 2 n'ont rencontré aucun réflecteur comme prévu et enfin pour 1 seul d'entre eux, le forage destructif n'a pas permis d'identifier le réflecteur attendu. Un début d'étude sur les traces et sur les spectres, en particulier les spectres glissants, semble ouvrir des voies prometteuses pour une interprétation des radargrammes, avec un nombre limité de forages ultérieurs. Le cinquième chapitre présente une synthèse des résultats et répond notamment aux questions posées par les directeurs du projet. Une carte de toutes les cavités détectées par l'ensemble des techniques est présentée. Il y en a près de 200. Aucune zone de la Médina n'est épargnée. En revanche, cette présence n'interdit pratiquement aucun projet d'aménagement ou d'urbanisme. A chaque usage d'un lieu, un traitement spécifique est proposé, dont le coût devrait rester faible par comparaison avec celui de la construction elle-même. Des propositions sont également faites concernant les réseaux, en bon accord avec des solutions de principe proposées précédemment par des bureaux d'étude compétents. En bref, les cavités passent du statut de danger inéluctable et imprévisible, à celui d'éléments apprivoisés de la vie de la médina.