The present work aims at enhancing the external quantum efficiencies of ultra-violet (UV) sensitive photodetectors (PDs) and light emitting diodes (LEDs)for any light polarization. Deep UV solid state devices are made out of AlGaN or MgZnO and their performances suffer from the high resistivity of their p-doped regions. They require transparent p-contacts; yet the most commonly used transparent contacts have low transmission in the UV: indium tin oxide (ITO) and nickel-gold (Ni/Au 5/5 nms) transmit less than 50% and 30% respectively at 300 nm. Here we investigate the use of surface plasmons (SPs) to design transparent p-contacts for AlGaN devices in the deep UV region of the spectrum. The appeal of using surface plasmon coupling arose from the local electromagnetic field enhancement near the metal surface as well as the increase in interaction time between the field and semiconductor if placed on top of a semiconductor. An in/out-coupling mechanism is achieved by using a grating consisting of two perpendicularly oriented sets of parallel aluminum lines with periods as low as 250 nm. The incident light is first coupled into SPs at the air/aluminum interface which then re-radiate at the aluminum/AlGaN interface and the photons energy is transferred to SP polaritons (SPPs) and back to photons. High transmission can be achieved not only at normal incidence but for a wider range of incident angles. iv A finite difference time domain (FDTD) package from R-Soft was used to simulate and design such aluminum gratings with transparency as high as 100% with tunable peak wavelength, bandwidth and angular acceptance. A rigorous coupled wave analysis (RCWA) was developed in Matlab to validate the FDTD results. The high UV transparency meshes were then fabricated using an e-beam assisted lithography lift-off process. Their electrical and optical properties were investigated. The electrical characterization was very encouraging; the sheet resistances of these meshes were lower than those of the conventionally used transparent contacts. The optical transmissions were lower than expected and the causes for the lower measurements have been investigated. The aluminum oxidation, the large metal grain size and the line edge roughness were identified as the main factors of inconsistency and solutions are proposed to improve these shortcomings. The effect of aluminum oxidation was calculated and the passivation of aluminum with SiO2 was evaluated as a solution. A cold deposition of aluminum reduced the aluminum grain size from 60 nm to 20 nm and the roughness from 5 nm to 0.5 nm. Furthermore, replacing the conventional lift-off process by a dry back-etch process led to much smoother metal line edges and much high optical transparency. The optical measurements were consistent with the simulations. Therefore, reduced roughness and smooth metal line edges were found to be especially critical considerations for deep UV application of the meshes.

The depths to buried bodies, characterized by anomalous gravity and magnetic properties, are determined by a combination of two numerical techniques. An upward continuation integral is solved by a method by Paul and Nagy using elemental squares and low order polynomials to describe the behavior of the gravity or magnetic data between observed data points. Downward continuation of the magnetic or gravity data is done by a finite difference technique as described by Bullard and Cooper. The applicability of the techniques are determined by comparison to depths determined by other means over the same anomalies and by comparison to various rule-of-thumb methods prevalent in the geophysical literature. The relative speed and cost of the particular computer system used is also considered in the applicability. The results show that although the initial costs of the computer program are high, the combined technique is as good as and at times better than the rule-of-thumb methods in determining the depth to the anomaly-causing body and is useful when more than just an approximate depth is of interest.

This thesis explores the links between different ways of modelling the physical world. Finite difference numerical simulations may be used to encode the behaviour of physical systems, allowing us to gain insight into their workings and even to predict their behaviour. Similarly, one can investigate the properties of gravitational black holes through the use of analogue black holes, physical systems which share at least some part of the physics of the astronomical objects. Concentrating on black hole analogues using Bose-Einstein condensates, I show how simulations of these systems may be greatly assisted through the use of a proper absorbing boundary condition, the Perfectly Matched Layer. Such a boundary condition allows the effcient truncation of the computational domain, both saving computational time and increasing accuracy. I then apply this technique to the simulation of the supersonic flow of a Bose-Einstein condensate through a Laval nozzle, a black hole analogue, showing that such a flow should be stable and observable in the laboratory. Moving to a related system, I investigate the optical analogue of the Iordanskii force - the friction resulting from interaction between excitations in a superfluid's normal component and a superfluid vortex - through the simulation of such a vortex in a Bose-Einstein condensate illuminated by slow light, which is light whose group velocity is on the order of metres per second. The interaction of the slow light with the vortex should produce a momentum transfer due to the optical Aharonov-Bohm effect, exerting a force on the vortex. The coupled system of equations describing the condensate-slow light system is simulated, giving some surprising results.

Les préoccupations environnementales actuelles visent à réduire les consommations énergétiques. Dans une démarche d’amélioration des bâtiments existants, l'étude du comportement thermique d’une paroi n'est pas aisée du fait de la méconnaissance de ses propriétés thermophysiques réelles. Ces paramètres sont pourtant prépondérants pour la phase d'optimisation économique des opérations de réhabilitation ou pour vérifier ses performances in situ. Il apparaît donc important de pouvoir caractériser les parois de bâtiment en place. Notre travail vise à développer une méthode de caractérisation thermique d’une paroi adaptée aux applications in situ basée sur une approche active. Le principe d'identification consiste à solliciter thermiquement une face d’accès en imposant un flux de chaleur sous forme d’un créneau et à étudier la réponse en température enregistrée par thermographie infrarouge sur l’autre face. A partir de signaux de flux et de températures mesurés aux limites de la paroi, les propriétés thermophysiques de la paroi seront estimées par méthode inverse. Nous nous sommes dans un premier temps intéressés aux parois homogènes. Le schéma d’inversion est construit autour d’un modèle numérique décrivant la réponse de la paroi suivant la méthode des différences finies en 1D. L’identification de la conductivité thermique et de la chaleur volumique de la paroi est réalisée en optimisant le groupement de paramètres qui permet de minimiser l’écart entre la température normalisée mesurée et la température normalisée simulée. Le coefficient d’échange surfacique global est également identifié à partir du même essai. Dans ce travail, la méthode a été appliquée à une paroi homogène en carreaux de plâtre mise en place au laboratoire. Elle a une épaisseur de 6.5 cm. Cette technique a été utilisée pour les parois multicouches de bâtiments. Les résultats issus de cette procédure d’inversion ont été comparés à des valeurs de référence obtenues à partir d’une procédure classique (NF EN 12664-méthode fluxmétrique). Une bonne concordance des résultats est obtenue. Une autre partie représente les essais in situ. / Current environmental concerns are intended to reduce energy consumption. In a process of improving existing buildings, the study of the thermal behavior of a wall is not easy because of the ignorance of its real thermophysical properties. These parameters are yet to dominate the economic optimization phase of the rehabilitation or to check its performance in situ. It therefore appears important to characterize the walls of existing building. Our work aims to develop a method of thermal characterization of a wall suitable for in situ applications based on an active approach. The principle of identification is to apply a heat-face access by imposing a heat flux in the form of a pulse and to study the temperature response recorded by infrared thermography on the other side. From signal flow and temperature measured at the limits of the wall, the thermophysical properties of the wall will be estimated by inverse method. We are at present interested in homogeneous walls. The inversion scheme is built around a digital model describing the response of the wall following the finite difference method in 1D. The identification of the thermal conductivity and heat volume of the wall is achieved by optimizing the group of parameters which minimizes the normalized difference between the temperature measured and the temperature standard simulated. The overall Global exchange coefficient is also identified from the same test. In this work, the method was applied to a homogeneous wall tile plaster introduction to the laboratory. It has a thickness of 6.5 cm. This technique was used for multilayer walls of buildings. The results of this inversion procedure were compared with reference values obtained from a standard procedure (DIN EN 12664-flow meter methods). A good agreement is obtained. Another part is the in situ tests.

Les temps de trajet, amplitudes et angles à la source des ondes sismiques sont utilisés dans de nombreuses applications telles que la migration, la tomographie, l'estimation de la sensibilité de détection et la localisation des microséismes. Dans le contexte de la microsismicité, il est nécessaire de calculer en quasi temps réel ces attributs avec précision. Nous avons développé ici un ensemble d'algorithmes rapides et précis en 3D pour des modèles à fort contraste de vitesse.Nous présentons une nouvelle méthode pour calculer les temps de trajet, les amplitudes et les angles à la source des ondes correspondant aux premières arrivées. Plus précisément, nous résolvons l'équation Eikonal, l'équation de transport et l'équation des angles en nous basant sur une approche par différences finies pour des modèles de vitesse en 3D. Nous proposons une nouvelle méthode hybride qui bénéficie des avantages respectifs de plusieurs approches existantes de résolution de l'équation Eikonal. En particulier, les approches classiques proposent généralement de résoudre directement les équations et font l'approximation localement d'une onde plane. Cette approximation n'est pas bien adaptée au voisinage de la source car la courbure du front d'onde est importante. Des erreurs de temps de trajet sont alors générées près de la position de la source, puis propagées à travers tout le modèle de vitesse. Ceci empêche de calculer correctement les amplitudes et les angles à la source puisqu'ils reposent sur les gradients des temps. Nous surmontons cette difficulté en introduisant les opérateurs sphériques ; plus précisément nous reformulons les temps de trajet, amplitudes et angles à la source par la méthode des perturbations.Nous validons nos nouvelles méthodes pour différents modèles à fort contraste de vitesse en 2D et 3D et montrons notre contribution par rapport aux approches existantes. Nos résultats sont similaires à ceux calculés en utilisant la modélisation de la forme d'onde totale alors qu'ils sont bien moins coûteux en temps de calcul. Ces résultats ouvrent donc de nouvelles perspectives pour de nombreuses applications telles que la migration, l'estimation de la sensibilité de détection et l'inversion des mécanismes au foyer. / Traveltimes, amplitudes and take-off angles of seismic waves are used in many applications such as migration, tomography, detection sensitivity estimation and microseism location. In the microseismicty context it is necessary to compute in near real time accurately these attributes. Here we developed a set of fast and accurate algorithms in 3D for highly contrasted velocity models.We present a new accurate method for computing first arrival traveltimes, amplitudes and take-off angles; more precisely we solve the Eikonal, transport and take-off angle equations based on a finite difference approach for 3D velocity models. We propose a new hybrid method that benefits from the advantages of several existing Eikonal solvers. Common approaches that solve directly these equations assume that we are locally propagating a plane wave. This approximation is not well adapted in the neighborhood of the source since the wavefront curvature is important. Travel times errors are generated near the source position and then propagated through the whole velocity model. This prevents from properly calculating the amplitudes and the take-off angles since they rely on the travel time gradients that are not accurate. We overcome this difficulty by introducing spherical operators. Indeed we reformulate the traveltimes, amplitudes and take-off angles with the perturbation method.We validate our new methods on various highly contrasted velocity models in 2D and 3D and show our contribution compared to other existing approaches. Our results are similar to those computed using full waveform modeling while they are obtained in a much shorter CPU time. These results open thus new perspectives for several applications such as migration, detection sensitivity estimation and focal mechanism inversion.

Le présent travail de thèse porte sur la modélisation et la simulation, en domaine temporel, de l'interaction entre un faisceau d'électrons et une onde hyperfréquence dans la structure à onde lente d'un TOP à hélice. Puisque le TOP est un instrument surdimensionné, les modèles non-stationnaires généraux utilisés dans les codes commerciaux nécessitent de trop grosses ressources de calcul pour pouvoir être utilisés en un temps raisonnable dans un but de conception. Il est donc nécessaire de faire appel à des modèles spécialisés. Durant cette thèse, nous nous sommes intéressés au "modèle discret non-stationnaire d'excitation d'un guide d'onde Périodique" de S. Kuznetsov. En 2007, N. Ryskin et al. avaient prouvé que ce modèle pouvait convenablement s'appliquer aux TOP à cavités couplées dans le cadre d'une application à une dimension du modèle. Lors de cette thèse, nous avons démontré, via le développement d'un code à une dimension (HelL-1D), que le modèle discret s'applique convenablement aux TOP à hélice. L'implémentation de ce modèle, dans un code à deux dimensions (HelL-2D) a, elle aussi, été effectuée. Enfin, nous avons développé une méthode permettant de contrôler de manière quantitative, dans le modèle discret, les phénomènes de réflexions aux extrémités de la ligne à retard, qui peuvent jouer un rôle important dans la stabilité de l'instrument. / This Ph.D. work deals with the time domain modeling and simulation of the electron beam/wave interaction in the slow-wave structure of a helix traveling-wave tube. Since a TWT is a device of which the geometry is oversized then the commercial software based on non-stationary general models needs so much computational resources that it cannot be used for design activities. During this Ph.D., we focused on the so called Kuznetsov's "discrete model of excitation of a periodic waveguide" which is a specialized model of beam/wave interaction in TWT. By 2007, N. Ryskin et al., showed that this model could conveniently apply to TWT with coupled cavities structure in one dimension. During this thesis, we first demonstrated that the discrete model could also apply to helix TWT with a sufficient (1%) accuracy via the development of a one-dimensional software called HelL-1D. We also implemented the discrete model for helix TWT in two dimensions (HelL-2D code). Finally, we developed a method, which is an extension of the discrete model, and which permits to take into account quantitatively the reflection phenomena at the terminations of a slow-wave structure in this model. This last study was very important since the stability of TWT strongly depends on this parameter.

Uma grande limitação dos métodos de diferenças finitas é que eles estão restritos a malhas e domínios retangulares. Para descrever escoamentos em domínios complexos, como, por exemplo, problemas com superfícies livres, faz-se necessário o uso de técnicas acessórias. O método de interfaces imersas é uma dessas técnicas. Nesse trabalho, primeiramente foi desenvolvido um método de projeção, totalmente livre de pressão, para as equações de Navier-Stokes com variáveis primitivas em malha deslocada. Esse método é baseado em diferenças finitas compactas, possuindo segunda ordem temporal e quarta ordem espacial. Esse método foi combinado com o método de interface imersa de Linnick e Fasel [2] para resolver numericamente as equações de Stokes com quarta ordem de precisão. A verificação do código foi feita por meio do método das soluções manufaturadas e da comparação com resultados de outros autores em problemas clássicos da literatura. / A great limitation of finite differences methods is that they are restricted to retangular meshes and domains. In order to describe flows in complex domains, e.g. free surface problems, it is necessary to use accessory techniques. The immersed interface method is one of such techniques. In the present work, firstly, a projection method was developed, which is completely pressure-free, for the Navier-Stokes equations with primitive variables in a staggered mesh. This method is based on compact finite differences, with temporal second-order precision and spatial foruth-order precision. This method was combined with the immersed interface method from Linnick e Fasel [2] in order to numerically solve the Stokes equations with fourth-order precision. The verification of the code was performed with the manufactured solutions method and by comparing results with other authors for some classical problems in the literature.

Made available in DSpace on 2015-03-04T18:50:39Z (GMT). No. of bitstreams: 1 thesis.pdf: 4735629 bytes, checksum: efd1eaf8c6c6b15b99b84eb0d25bd102 (MD5) Previous issue date: 2012-09-06 / Coordenacao de Aperfeicoamento de Pessoal de Nivel Superior / Multiphase flows in porous media are modeled by a system of partial differential equations and the study of the numerical approximation to the solutions of these plays a crucial role in the simulation and prediction of problems that are of great practical interest and of economic and social impact, such as secondary oil recovery, geological storage of CO2 and transport of pollutants in aquifers. The goal of this work is the development of a three-dimensional numerical simulator that precisely evaluates the transport of two immiscible fluids in a heterogeneous porous media using multithread parallel programming to shared memory multiprocessors computers. The system of partial differential equations is decomposed into a elliptic subsystem used to determine the velocity field and into a hyperbolic equation (nonlinear) to determine the transport of the fluid phases. The approximation to the solution of the latter one is calculated using a high order non-oscillatory finite-differences numerical method based on central schemes that allows a semi-discrete formulation which an extension that enables to work with variable space coefficients. Numerical experiments on three-dimensional models were performed considering linear and nonlinear flow problems in typical settings of oil reservoirs simulations. The results were satisfactory since they presented mass conservation, precise capture of shock waves and small numeric diffusion, regardless of the time step. / Escoamentos multifásicos em meios porosos são modelados por um sistema de equações diferenciais parciais e o estudo da aproximação das soluções dessas equações desempenha papel crucial na simulação e previsão de problemas de grande interesse prático e impacto econômico e social, tais como a recuperação secundária de petróleo, o armazenamento geológico de CO2 e o transporte de poluentes em aquíferos. O presente trabalho tem como objetivo o desenvolvimento de um simulador numérico tridimensional para avaliar com precisão o transporte de dois fluidos imiscíveis em um meio poroso heterogêneo e que utiliza computação paralela multithread para computadores multiprocessados de memória compartilhada. O sistema de equações diferenciais parciais é decomposto em um subsistema elíptico para a determinação do campo de velocidades dos fluidos e uma equação hiperbólica não-linear para o transporte das fases fluidas. Para esta última, foi utilizado um método numérico de volumes finitos, não-oscilatório de alta ordem baseado em esquemas centrais e que admite uma formulação semi-discreta com coeficientes variáveis no espaço. Experimentos numéricos em modelos tridimensionais foram realizados considerando problemas de escoamentos lineares e não lineares postos em configurações típicas de simulação de reservatórios de petróleo. Os resultados mostraram-se satisfatórios por apresentarem conservação da massa, boa captura das ondas de choque e pequena difusão numérica, independente do passo de tempo.

Dans la présente étude, la conductivité thermique et le module d'élasticité de revêtementsd’YPSZ élaborés par projection plasma ont été prédits par modélisations numériques 2D et3D de type différences finies et éléments finis.L'influence de la résolution d'image, de la taille et de la valeur du seuil sur les propriétésprédites du revêtement a été étudiée. En outre, les effets de la méthode numérique et du typede condition aux limites ont été étudiés. En particulier, la quantification de l'effet Knudsen(effet de raréfaction) sur le transfert de chaleur à travers une structure poreuse a été réaliséepar modélisation numérique en combinaison avec l'analyse d'image. Les conductivitéseffectives obtenues par modélisation 3D s'avèrent plus élevées que celles obtenues en 2D, etaussi en meilleur accord avec les résultats mesurés. Une corrélation 2D/3D a été trouvéepour la modélisation de la conductivité thermique : cette corrélation permet de prédire lesvaleurs 3D à partir des valeurs calculées en 2D. / In the present study, the thermal conductivity and elastic modulus of thermal spray YPSZcoatings were predicted by 2D and 3D finite differences and finite elements numericalmodeling based on cross-sectional images.The influence of the image resolution, size and threshold on the predicted properties of thecoating was studied. Moreover, the effects of the numerical method and of the boundarycondition were investigated. In particular, the quantification of the Knudsen effect(rarefaction effect) on the heat transfer through a porous structure was realized by numericalmodeling in combination with image analysis. The predicted thermal conductivities obtainedby 3D modeling were found to be higher than those obtained by 2D modeling, and in betteragreement with the measured results. A 2D/3D correlation was sucessfully found for themodeling of thermal conductivity: this correlation allows predicting 3D computed valuesfrom 2D ones.

Le réverbérateur à plaque est un dispositif électromécanique de réverbération artificielle utilisé dans les studios d'enregistrement pour traiter les signaux audio enregistrés dépourvus d'effet de salle. Ce dispositif simule la réverbération à partir des vibrations de flexion d'une plaque mince. L'objectif de cette thèse est la simulation par modèles physiques d'un réverbérateur à plaque générique permettant à l'utilisateur le choix des principaux paramètres : caractéristiques géométriques et physiques du matériau de la plaque, conditions aux limites, positions de l'excitateur et des accéléromètres, et l'influence d'une plaque poreuse à proximité. Outre la compréhension des phénomènes entrant en jeu dans de tels dispositifs, cette approche permet la simulation de réverbérateurs à plaque existants comme l'EMT-140 ou l'EMT-240, mais aussi de réverbérateurs avec d'autres jeux de paramètres réels ou fictifs. La démarche entreprise s'articule en trois points : la modélisation du dispositif, les mesures sur un réverbérateur EMT-140 et la définition des outils numériques pour la synthèse sonore dans le domaine audible. La modélisation développée relève de la dynamique et de la vibroacoustique, et porte une attention particulière sur les mécanismes d'amortissement en raison de leur grande influence dans la perception de l'effet de réverbération. La démarche expérimentale, effectuée sur un réverbérateur EMT-140, permet l'identification des paramètres d'un dispositif réel et l'identification des limites de validité de la modélisation. L'outil de simulation dans le domaine temporel repose sur la méthode des différences finies.