Made available in DSpace on 2014-12-17T14:55:39Z (GMT). No. of bitstreams: 1 LincolnMA.pdf: 1667449 bytes, checksum: b8113389f31903ba22cf94dbc22192c1 (MD5) Previous issue date: 2009-08-13 / Conselho Nacional de Desenvolvimento Cient?fico e Tecnol?gico / This work presents a theoretical and numerical analysis of structures using frequency selective surfaces applied on patch antennas. The FDTD method is used to determine the time domain reflected fields. Applications of frequency selective surfaces and patch antennas cover a wide area of telecommunications, especially mobile communications, filters and WB antennas. scattering parameters are obteained from Fourier Transformer of transmited and reflected fields in time domain. The PML are used as absorbing boundary condition, allowing the determination of the fields with a small interference of reflections from discretized limit space. Rectangular patches are considered on dielectric layer and fed by microstrip line. Frequency selective surfaces with periodic and quasi-periodic structures are analyzed on both sides of antenna. A literature review of the use of frequency selective surfaces in patch antennas are also performed. Numerical results are also compared with measured results for return loss of analyzed structures. It is also presented suggestions of continuity to this work / Este trabalho apresenta uma an?lise te?rica e num?rica de estruturas que utilizam superf?cies seletivas de frequ?ncia aplicadas a antenas do tipo patch. Para isso, ? utilizado o m?todo das diferen?as finitas no dom?nio do tempo (FDTD) visando determinar os campos refletidos a partir de uma onda plana incidente no dom?nio do tempo. As aplica??es das superf?cies seletivas de freq??ncia e antenas patch abrangem uma grande ?rea das Telecomunica??es, principalmente em comunica??es m?veis e v?o desde filtros at? as antenas banda larga. Especificamente, a an?lise usa os campos transmitidos e refletidos obtidos no dom?nio do tempo, em conjunto com transformada de Fourier permitindo a obten??o dos par?metros de transmiss?o da antena. A condi??o de contorno absorvedora utilizada foi a de camada perfeitamente casada (PML), permitindo a determina??o num?rica dos campos com uma quantidade menor de interfer?ncias provenientes de reflex?es nos limites do espa?o discretizado. S?o considerados patches retangulares condutores sobre uma camada diel?trica e alimentados por linha de microfita. Foram analisadas superf?cies seletivas de frequ?ncia peri?dicas e quase peri?dicas tanto no plano de terra quanto no plano do pr?prio patch. ? realizada uma revis?o bibliogr?fica a respeito da utiliza??o de superf?cies seletivas de frequ?ncia em antenas patch. Tamb?m s?o comparados resultados num?ricos e medidos para a perda de retorno das estruturas analisadas. S?o apresentadas, ainda, sugest?es de continuidade para este trabalho

Neste trabalho foi realizado um levantamento GPR 4D em laboratório no qual foi simulado o vazamento de óleo de uma ETD - Estação de Transmissão e Distribuição de Energia Elétrica. Os estudos consistiram em derramar 15 litros de óleo em três experimentos: tanque contendo areia seca, tanque contendo areia úmida mais gradiente hidráulico e tanque contendo areia úmida. Em todos os experimentos o objetivo foi detectar a pluma de contaminação e avaliar a migração desse óleo com o tempo. Os dados GPR 4D foram adquiridos utilizando-se uma antena blindada de 400 MHz modelo SIR-3000 equipamento da GSSI. O monitoramento no tempo de aquisição dos dados variou de 2 minutos até 12 dias. No experimento com o tanque preenchido com areia seca não foi possível determinar a pluma de contaminação devido ao baixo contraste entre as propriedades físicas do meio e a pluma contaminante. No experimento em que o tanque foi preenchido com areia úmida e havia a presença de um gradiente hidráulico, foi possível determinar a migração da pluma de contaminação a partir do 5º dia. Para finalizar, o 3º experimento onde o tanque estava cheio de areia úmida também não foi possível identificar a pluma, devido ao baixo contraste entre as propriedades físicas. Nesta pesquisa também foram realizadas modelagens numéricas 2D utilizando o método FDTD, onde foram simulados os três experimentos, considerando posições intermediárias da pluma de contaminação em função do tempo. Os resultados das modelagens foram concordantes com os resultados obtidos com dados reais e ajudaram a definir o padrão de reflexão da pluma contaminante. Os promissores resultados indicam que a determinação efetiva da pluma contaminante de óleo numa ETD é possível desde o meio que esteja saturado com água. Portanto, para garantir o sucesso das pesquisas GPR para fins de mapeamento e delineamento de plumas de contaminação em ETDs recomenda-se que a aquisição dos dados seja feita após um período de chuvas. / In this work a 4D GPR survey was accomplished in a laboratory in which a oil leak from a ETD Estação de Transmissão e Distribuição of electric energy was simulated. The studies consisted in an oil spill (15 liters) in three different experiments: i) tank containing dry sand, ii) tank containing wet sand with a hydraulic gradient and iii) tank containing wet sand. In all experiments the objective was detect the contamination plume and evaluate the migration of the oil with the time. The 4D GPR data was acquired using a shielded antenna of 400 MHz model SIR-3000 of the GSSI equipment. The time monitoring of the data acquisition ranged from 2 minutes up to 12 days. In the experiment-1 with the tank filled with dry sand wasnt possible to determine the plume contamination due to the low contrast of the physical properties between the environment and the plume. In the experiment-2 where the tank filled with the dry sand in the presence of a hydraulic gradient, was possible to determine the migration of the plume from the fifth day. Finally, the experiment-3 where the tank was filled with wet sand wasnt possible to identify the plume, due to the low contrast between the physical properties. This research also conducted 2D numerical modeling using the FDTD method, where the three experiments were simulated, considering intermediary positions of the contamination plume with respect with time. The results of the modeling were consistent with the results obtained with real data and helped to define a pattern of reflection of the contamination plume. The promising results indicate that the effective determination of the contaminating plume of oil in an ETD is possible provided that the environment be saturated with water. Therefore, to guarantee the success of researches with GPR for mapping and delineation of contamination plumes in ETDs its recommended that the data acquisition be realized in the rainy season.

Notre recherche porte sur la réponse optique des nanostructures et nous sommes certainement intéressés à la modélisation de ces structures afin d'améliorer le confinement de la lumière. Ce confinement est un des paramètres qui conduisent à l'exaltation des effets optique linéaires et non linéaires, simultanément. Notre travail est divisé en deux sections, qui présentent deux effets optiques diffèrents basées sur le confinement des champs à l’intérieur du structure. Dans la première section, la réponse optique d’une nano-antenne à ouverture en forme de nœud papillon (BNA : Bowtie nano-aperture) sera étudiée en fonction de l’indice de réfraction du milieu. L’étude discute la variation de la longueur d’onde de résonance ainsi que l’intensité du champ confiné au milieu du gap de la BNA en fonction de la distance qui sépare l’antenne d’un substrat placé en face. L’étude prend le cas d’une BNA gravée au bout d’une fibre optique métallisée. Une étude numérique complète a été réalisée par des simulations numériques basées sur la méthode des différences finies FDTD-3D (Finite Difference Time Domaine – Three dimension). Le code utilisé est développé au sein du département d’Optique de l’Institut FEMTO-ST. Notre modèle numérique décrits bien la géométrie de l’antenne ainsi que la pointe. De même les propriétés optiques de la couche métallique de la sonde sont bien décrites à travers un modèle de dispersion (model de Drude). Une étude expérimentale a été réalisée en plus pour une validation des résultats collecte théoriquement, l’étude a considéré le cas de deux BNA avec des paramètres géométrique différents, néanmoins les résultats obtenues se correspond fortement aux résultats obtenues numériquement. Dans la deuxième partie, l'effet électro-optique des nanostructures sera étudié. L’effet électro-optique ou effet Pockels consiste d’une variation linéaire de l'indice de réfraction d'un milieu non-linéaire en fonction d'un champ électrique extérieur. Cependant, la variation est reliée au tenseur de susceptibilité non linéaire du deuxième ordre, ainsi, cette effet ce produit seulement dans les matériaux non-centrosymétrique. Dans notre étude les nanostructures sont fabrique avec le Niobate de Lithium (LN) qui est considéré comme le plus avantageux diélectrique pour l’exaltation des effets non linéaire grâce à ces propriétés (acousto-optique, électro-optique, piézoélectrique …). L’étude est complétée numériquement garce a des simulations basées sur la FDTD en tenant compte de la polarisation du champ. Au début, une validation de l’utilisation de la FDTD pour estimer l’effet électro-optique intrinsèque du LN a été achevée, l’étude considère un milieu homogène (pas de confinement des champs dans le milieu) qui a pour indice de réfraction celle du LN. Pour des structure qui présentes un confinement du champ plusieurs approximation ont été suggère dans cette partie. En outre, nous présentons un nouveau modèle auto-cohérent dans lequel la variation de l'indice de réfraction est modifiée au cours de la simulation. Plusieurs structure ont été discutés (réflecteur de Bragg, structure à cavité et 2D cristal photonique) qui présentent des facteurs de confinement différents. Une étude comparative entre les différentes méthodes, montre que la différence entre les résultats de chaque hypothèse devient plus important proportionnellement au facteur de confinement. / Our research is concerned with the optical response of nano-structures by modeling them in order to enhance the confinement of light in these structures, which leads to the exaltation of linear and nonlinear optical effects.Our work is divided into two sections, which are based on the enhancement of the electric field inside the structure. In the first section, we study the optical properties of a Bowtie Nano-aperture, BNA, as a function of the refractive index of the surrounding medium. The study discusses the variation of the resonance wavelength and the intensity of the enhanced field in the gap of the BNA as a function of the distance from a sample placed in front of our BNA. The BNA is engraved at the apex of a metallic coated fiber tip. In this section a theoretical study was achieved using the Finite Difference Time Domain method FDTD in which we implement a Drude dispersion model to faithfully describe the optical properties of metals. In addition, a validating experimental study was achieved and a high accordance between both results is recorded.In the second section, the electro optical effect of nano-structures is studied. Electro-optical effect or Pockels effect is the variation of the refractive index of a nonlinear media as a function of an applied external electric field. The electro-optical effect is a linear variation of the media refractive index. However it is also related to the second order nonlinear susceptibility tensor, thus it becomes a nonlinear effect that only occurs in non-centrosymetric material. In our study we chose the case of a nano-structure fabricated with Lithium Niobate. Lithium Niobate is widely used in photonic applications due to its electro-optical, acousto-optical and nonlinear optical properties. We present a theoretical study of the electro-optical effects using the FDTD simulation method. We started by approving the ability to use the FDTD to calculate the refractive index variation in bulk Lithium Niobate then we suggest different approximations to estimate the refractive index variation when the light is confined inside the structure. In addition we suggest a new self-consistent method in which the variation of the refractive index is modified during the simulation. The study shows a comparison between different assumptions (used in previous research) and the self-consistent method for various structures, like Bragg reflectors, cavity structures and 2D photonic crystals. The study shows that the difference between the results of each assumption becomes greater when the optical confinement in the structure becomes more important.

Dans le domaine de la Compatibilité ElectroMagnétique (CEM), la connaissance du niveau d’immunité ou de susceptibilité d’un équipement embarqué est un paramètre important à contrôler et à gérer tout au long de sa durée de vie. En particulier,lors des opérations de maintenance de cet équipement, il est important de veiller au maintien de son niveau de susceptibilité pour garantir son fonctionnement. Pour répondre à ce besoin, les travaux de cette thèse proposent une méthode expérimentale originale de contrôle rapide et facile à mettre en œuvre dans le cas d’essais de susceptibilité rayonnée en haute fréquence. Cette méthode est basée sur l’observation de fréquences harmoniques rayonnées par un équipement lorsque celui-ci est soumis une illumination parfaitement monochromatique en haute fréquence. Pour cela, on fait l’hypothèse que l’amplitude de ces fréquences est d’autant plus grande que l’équipement électronique est en dysfonctionnement. Pour vérifier cette hypothèse, un dispositif expérimental adapté à une illumination en CRBM a tout d’abord été développé et étudié afin de mettre en œuvre cette méthode. Ensuite, après avoir mis en évidence sur plusieurs configurations de cartes électronique l’observation d’harmoniques, nous nous sommes intéressés à une modélisation du phénomène physique pour vérifier celui-ci par la simulation. Ace titre, une modélisation FDTD des expériences a été réalisée en 3D et une comparaison mesures/calcul a été effectuée. Ce travail a permis de montrer qu’à partir du des équations de Maxwell, des harmoniques identiques à celles mesurées étaient rayonnées et que la méthode de contrôle de susceptibilité rayonnée que nous proposions était donc tout à fait viable. Concernant la phase de modélisation, nous nous sommes intéressés essentiellement à représenter le fonctionnement du circuit logique et non à entrer finement dans le détail de chaque composant électronique le constituant. Pour les applications retenues, nous avons donc porté notre effort sur la prise en compte d’un signal numérique représenté sous forme d’une suite de bits, ainsi que sur l’introduction dans le modèle FDTD de Maxwell de quelques composants non linéaires comme la diode et un inverseur CMOS. La généralisation du modèle d’inverseur peut s’appliquer sans difficultés à d’autres composants CMOS, comme les mémoires de type SRAM que nous avons étudiées expérimentalement. / In the frame of ElectroMagnetic Compatibility (EMC), the knowledge of EM immunity or susceptibility levels of on-board equipment is an important parameter to control and manage throughout the whole life of equipment. In fact, during maintenance phases of this equipment, it is important to ensure that these levels are maintained in order to guarantee its operation. To meet this need, this thesis proposes to focus on a novel and easy-to-implement experimental method for a quick diagnosis of equipment EM radiated susceptibility at high frequency.. This method is based on the observation of the harmonic frequencies radiated by an equipment device when this one is summited to a perfectly monochromatic illumination at high frequency. We make the hypothesis that the amplitude of those frequencies increases when the electronic equipment is in a dysfunction state. For this, a experimental method suitable for MSRCs has first been developed and studied to highlight the principle of this method. Then, having observed experimentally the radiation phenomenon of harmonic frequencies on several configurations of electronic cards, we have developed a 3D FDTD model to simulate the phenomenon. With this model, a modeling of the experiments was carried out and comparisons between measurements and calculations were performed. Calculation results coming from FDTD Maxwell’s equation models showed trends and functional behaviors identical to those observed in measurements which demonstrated that the proposed control method was fully viable. Regarding the modeling phase, we concentrated on models describing the operation of the logic circuit but we did not to enter into the physical description of each electronic component. Our selected applications have therefore focused on the inclusion in the Maxwell FDTD model of a digital signal represented as a series of bits, as well as the introduction of some components such as a diode and a CMOS inverter. The generalization of the proposed inverter model can be easily extended to other CMOS components such as the SRAM memories we studied in experimental tests.

The efficiency of the conventional, explicit finite difference time domain (FDTD)method is constrained by the upper limit on the temporal discretization, imposed by the Courant–Friedrich–Lewy (CFL) stability condition. Therefore, there is a growing interest in overcoming this limitation by employing unconditionally stable FDTD methods for which time-step and space-step can be independently chosen. Unconditionally stable Crank Nicolson method has not been widely used in time domain electromagnetics despite its high accuracy and low anisotropy. There has been no work on the Crank Nicolson FDTD (CN–FDTD) method for frequency dependent medium. In this thesis a new three-dimensional frequency dependent CN–FDTD (FD–CN–FDTD) method is proposed. Frequency dependency of single–pole Debye materials is incorporated into the CN–FDTD method by means of an auxiliary differential formulation. In order to provide a convenient and straightforward algorithm, Mur’s first-order absorbing boundary conditions are used in the FD–CN–FDTD method. Numerical tests validate and confirm that the FD–CN–FDTD method is unconditionally stable beyond the CFL limit. The proposed method yields a sparse system of linear equations which can be solved by direct or iterative methods, but numerical experiments demonstrate that for large problems of practical importance iterative solvers are to be used. The FD–CN–FDTD sparse matrix is diagonally dominant when the time-stepis near the CFL limit but the diagonal dominance of the matrix deteriorates with the increase of the time-step, making the solution time longer. Selection of the matrix solver to handle the FD–CN–FDTD sparse system is crucial to fully harness the advantages of using larger time-step, because the computational costs associated with the solver must be kept as low as possible. Two best–known iterative solvers, Bi-Conjugate Gradient Stabilised (BiCGStab) and Generalised Minimal Residual (GMRES), are extensively studied in terms of the number of iteration requirements for convergence, CPU time and memory requirements. BiCGStab outperforms GMRES in every aspect. Many of these findings do not match with the existing literature on frequency–independent CN–FDTD method and the possible reasons for this are pointed out. The proposed method is coded in Fortran and major implementation techniques of the serial code as well as its parallel implementation in Open Multi-Processing (OpenMP) are presented. As an application, a simulation model of the human body is developed in the FD–CN–FDTD method and numerical simulation of the electromagnetic wave propagation inside the human head is shown. Finally, this thesis presents a new method modifying the frequency dependent alternating direction implicit FDTD (FD–ADI–FDTD) method. Although the ADI–FDTD method provides a computationally affordable approximation of the CN–FDTD method, it exhibits a loss of accuracy with respect to the CN-FDTD method which may become severe for some practical applications. The modified FD–ADI–FDTD method can improve the accuracy of the normal FD–ADI–FDTD method without significantly increasing the computational costs.

Tato diplomová práce se zabývá vlastnostmi plazmonických antén v oblasti vysokých vlnových délek. Důraz je kladen na popis rezonančních vlastností jednotlivých antén i antén uspořádaných do periodických polí. Těžiště práce spočívá v počítačovém modelování navýšení magnetického pole v blízkosti antén, které lze využít ve vysokofrekvenční elektronové paramagnetické rezonanci (HFEPR) k zesílení měřeného signálu. Autor se zabývá kvantifikací zesílení v anténách s odlišnou geometrií a navrhuje i geometrii vlastní. Značná část práce se také věnuje snaze rozlišit příspěvky k navýšení magnetického pole od různých zdrojů při měření HFEPR v uspořádání s dvojitou transmisí záření.

Selected methods of Raman spectroscopy, like surface-enhanced Raman spectroscopy and single molecule Raman spectroscopy, are described in this diploma work. The basis of two methods for numerical modelling of optical properties of micro- and nanoparticles are prefaced. The methods are Discrete Dipole Approximation and Finite Difference Time Domain. Micro- and nanoparticles are used in surface enhanced Raman spectroscopy and other nanospectroscopic methods. Further, the main instrumentation needed for Raman spectroscopy is described. The first part of experimental section of this work is numerical modelling of photonic nanojet that occures behind dielectric microparticles. This phenomenon leads to a new technique of confocal microscopy with Raman spectra measuring. The second experimental section contains results of Raman spectra measurement with beta-carotene and surface-enhanced Raman spectra of beta-carotene in silver-sol solution.

This thesis is dealing with problem and solution of simulation in room acoustics. At the beginning is introduction with sound as waves and his behavior in closed space. Following part of text is dealing with computional methods in room acoustics, statistics methods , ray-based methods, wave-based methods. Following are focused FDTD wave method, which serves as the basis for creating simulation algorithm.Last part of this work is practical sample of MATLAB aplication enviroment for simulation of sound waves in closed room by explicit sub-methods of method FDTD. The last section contains an example and discussion of the results of simulations.

This thesis aims to (i) design and optimize the geometry of magneto-photonic crystal based on ferromagnetic garnet in order to resonantly enhance the magneto-optical re- sponse, (ii) to explore the possibility of using magnetic shape memory alloy to build an optically active photonic element, using advanced FDTD modeling. A Faraday rotation of 180◦ was reached but with low values of transmissivity. An investigation of the origin of such high values of Faraday rotation led to a conclusion that such structure has to be highly sensitive towards a change of a refractive index of its surroundings. This was confirmed, and so further development of this structure can lead to an efficient concentra- tion detector. Three designs of optically active element utilizing deformation of magnetic shape memory material in the external magnetic field were numerically simulated. Two designs (photonic crystal with cylindric holes in a hexagonal lattice and self-standing foil with cylindric holes in a square lattice) proved to be efficient and worth of further development. 1

The aim of this thesis was to make a practical tool for low frequency analysis in room acoustics.The need arises from Acad’s experience that their results from simulations using raytracing software deviate in the lower frequencies when compared to field measurements inrooms. The tool was programmed in Matlab and utilizes the Finite Difference Time Domain (FDTD) method, which is a form of rapid finite element analysis in the time domain.A number of tests have been made to investigate the practical limitations of the FDTD method, such as numerical errors caused by sound sources, discretization and simulation time. Boundary conditions, with and without frequency dependence, have been analysed bycomparing results from simulations of a virtual impedance tube and reverberation room to analytical solutions. These tests show that the use of the FDTD method appears well suited for the purpose of the tool.A field test was made to verify that the tool enables easy and relatively quick simulations of real rooms, with results well in line with measured acoustic parameters. Comparisons of the results from using the FDTD method, ray-tracing and finite elements (FEM) showed goodcorrelation. This indicates that the deviations Acad experience between simulated results and field measurements are most likely caused by uncertainties in the sound absorption data used for low frequencies rather than by limitations in the ray-tracing software. The FDTDtool might still come in handy for more complex models, where edge diffraction is a more important factor, or simply as a means for a “second opinion” to ray-tracing - in general FEM is too time consuming a method to be used on a daily basis.Auxiliary tools made for importing models, providing output data in the of room acoustic parameters, graphs and audio files are not covered in detail here, as these lay outside the scope of this thesis. / Målet för detta examensarbete var att undersöka möjligheten att programmera ett praktisktanvändbart verktyg för lågfrekvensanalys inom rumsakustik. Behovet uppstår från Acadserfarenhet att resultat från simuleringar med hjälp av strålgångsmjukvara avviker i lågfrekvensområdeti jämförelse med fältmätningar i färdigställda rum. Verktyget är programmerati Matlab och använder Finite Difference Time Domain (FDTD) metoden, vilket är en typav snabb finita elementanalys i tidsdomänen.En rad tester har genomförts för att se metodens praktiska begräsningar orsakade av numeriskafel vid val av ljudkälla, diskretisering och simuleringstid. Randvillkor, med och utanfrekvensberoende, har analyserats genom jämförelser av simulerade resultat i virtuella impedansröroch efterklangsrum mot analytiska beräkningar. Testerna visar att FDTD-metodentycks fungerar väl för verktygets tilltänkta användningsområde.Ett fälttest genomfördes för att verifiera att det med verktyget är möjligt att enkelt och relativtsnabbt simulera resultat som väl matcher uppmätta rumsakustiska parametrar. Jämförelsermellan FDTD-metoden och resultat beräknade med strålgångsanalys och finita elementmetoden(FEM) visade även på god korrelation. Detta indikerar att de avvikelser Acaderfar mellan simulerade resultat och fältmätningar troligen orsakas av osäkerheter i den ingåendeljudabsorptionsdata som används för låga frekvenser, snarare än av begränsningar istrålgångsmjukvaran. Verktyget kan fortfarande komma till användning för mer komplexamodeller, där kantdiffraktion är en viktigare faktor, eller helt enkelt som ett sätt att få ett”andra utlåtande” till resultaten från strålgångsmjukvaran då FEM-analys generellt är en förtidskrävande metod för att användas på daglig basis.Kringverktyg skapade för t.ex. import av modeller, utdata i form av rumsakustiska parametrar,grafer och ljudfiler redovisas inte i detalj i denna rapport eftersom dessa ligger utanförexamensarbetet.