This thesis presents a very efficient hybrid method for analysis and design of optical and passive photonic devices. The main focus is on unbounded wave structures. This class of photonic systems are in general very large in terms of the wavelength of the driving optical sources. The size of the problem space makes the electromagnetic modelling of these structure a very challenging problem. Our approach and main contribution has been to combine or hybridize three methods that together can handle this class of photonic structures as a whole. <br /><br /> The basis of the hybrid method is a novel Gaussian Beam Tracing method GBT. Gaussian Beams (GB) are very suitable elementary functions for tracing and tracking purposes due to their finite extent and the fact that they are good approximations for actual laser beams. The GBT presented in this thesis is based on the principle of phase matching. This method can be used to model the reflection and refraction of Gaussian beams from general curved surfaces as long as the curvature of the surface is relatively small. It can also model wave propagation in free space. The developed GBT is extremely fast as it essentially uses simple algebraic equations to find the parameters of the reflected and refracted beams once the parameters of the incident beam is known. Therefore sections of the systems whose dimensions are large relative to the optical wavelength are simulated by the GBT method. <br /><br /> Fields entering a photonic system may not possess an exact Gaussian profile. For example if an aperture limits the input laser to the system, the field is no longer a GB. In these and other similar cases the field at some aperture plane needs to be expanded into a sum of GBs. Gabor expansion has been used for this purpose. This method allows any form of field distribution on a flat or curved surface to be expanded into a sum of GBs. The resultant GBs are then launched inside the system and tracked by GBT. Calculation of the coefficients of the Gabor series is very fast (1-2 minutes on a typical computer for most applications). <br /><br /> In some cases the dimensions or physical properties of structures do not allow the application of the GBT method. For example if the curvature of a surface is very large (or its radius of curvature is very small) or if the surface contains sharp edges or sub-wavelength dimensions GBT is no longer valid. In these cases we have utilized the Finite Difference Time Domain method (FDTD). FDTD is a rigorous and very accurate full wave electromagnetic solver. The time domain form of Maxwell's equations are discretized and solved. No matrix inversion is needed for this method. If the size of the structure that needs to be analyzed is large relative to the wavelength FDTD can become increasingly time consuming. Nevertheless once a structure is simulated using FDTD for a given input, the output is expanded using Gabor expansion and the resultant beams can then be efficiently propagated through any desired system using GBT. For example if a diffraction grating is illuminated by some source, once the reflection is found using FDTD, it can be propagated very efficiently through any kind of lens or prism (or other optical structures) using GBT. Therefore the overall computational efficiency of the hybrid method is very high compared to other methods.

This thesis presents a very efficient hybrid method for analysis and design of optical and passive photonic devices. The main focus is on unbounded wave structures. This class of photonic systems are in general very large in terms of the wavelength of the driving optical sources. The size of the problem space makes the electromagnetic modelling of these structure a very challenging problem. Our approach and main contribution has been to combine or hybridize three methods that together can handle this class of photonic structures as a whole. <br /><br /> The basis of the hybrid method is a novel Gaussian Beam Tracing method GBT. Gaussian Beams (GB) are very suitable elementary functions for tracing and tracking purposes due to their finite extent and the fact that they are good approximations for actual laser beams. The GBT presented in this thesis is based on the principle of phase matching. This method can be used to model the reflection and refraction of Gaussian beams from general curved surfaces as long as the curvature of the surface is relatively small. It can also model wave propagation in free space. The developed GBT is extremely fast as it essentially uses simple algebraic equations to find the parameters of the reflected and refracted beams once the parameters of the incident beam is known. Therefore sections of the systems whose dimensions are large relative to the optical wavelength are simulated by the GBT method. <br /><br /> Fields entering a photonic system may not possess an exact Gaussian profile. For example if an aperture limits the input laser to the system, the field is no longer a GB. In these and other similar cases the field at some aperture plane needs to be expanded into a sum of GBs. Gabor expansion has been used for this purpose. This method allows any form of field distribution on a flat or curved surface to be expanded into a sum of GBs. The resultant GBs are then launched inside the system and tracked by GBT. Calculation of the coefficients of the Gabor series is very fast (1-2 minutes on a typical computer for most applications). <br /><br /> In some cases the dimensions or physical properties of structures do not allow the application of the GBT method. For example if the curvature of a surface is very large (or its radius of curvature is very small) or if the surface contains sharp edges or sub-wavelength dimensions GBT is no longer valid. In these cases we have utilized the Finite Difference Time Domain method (FDTD). FDTD is a rigorous and very accurate full wave electromagnetic solver. The time domain form of Maxwell's equations are discretized and solved. No matrix inversion is needed for this method. If the size of the structure that needs to be analyzed is large relative to the wavelength FDTD can become increasingly time consuming. Nevertheless once a structure is simulated using FDTD for a given input, the output is expanded using Gabor expansion and the resultant beams can then be efficiently propagated through any desired system using GBT. For example if a diffraction grating is illuminated by some source, once the reflection is found using FDTD, it can be propagated very efficiently through any kind of lens or prism (or other optical structures) using GBT. Therefore the overall computational efficiency of the hybrid method is very high compared to other methods.

While exact methods, such as DDA or T-matrix, can be applied to particles withsizes comparable to the wavelength, computational demands mean that they are size limited. For particles much larger than the wavelength, the Geometric Optics approximation can be employed, but in doing so wave effects, such as interference and diffraction, are ignored. In between these two size extremes there exists a need for computational techniques which are capable of handling the wide array of ice crystal shapes and sizes that are observed in cirrus clouds. The Beam Tracing model developed within this project meets these criteria. It combines aspects of geometric optics and physical optics. Beam propagation is handled by Snell's law and the law of reflection. A beam is divided into reflected and transmitted components each time a crystal facet is illuminated. If the incident beam illuminates multiple facets it is split, with a new beam being formed for each illuminated facet. The phase-dependent electric field amplitude of the beams is known from their ampli- tude (Jones) matrices. These are modified by transmission and reflection matrices, whose elements are Fresnel amplitude coefficients, each time a beam intersects a crystal facet. Phase tracing is carried out for each beam by considering the path that its 'centre ray' would have taken. The local near-field is then mapped, via a surface integral formulation of a vector Kirchhoff diffraction approximation, to the far-field. Once in the far-field the four elements of the amplitude matrix are trans- formed into the sixteen elements of the scattering matrix via known relations. The model is discussed in depth, with details given on its implementation. The physical basis of the model is given through a discussion of Ray Tracing and how this leads to the notion of Beam Tracing. The beam splitting algorithm is described for convex particles followed by the necessary adaptations for concave and/or ab- sorbing particles. Once geometric aspects have been established details are given as to how physical properties of beams are traced including: amplitude, phase and power. How diffraction is implemented in the model is given along with a review of existing diffraction implementations. Comparisons are given, first against a modified Ray Tracing code to validate the geometric optics aspects of the model. Then, specific examples are given for the cases of transparent, pristine, smooth hexagonal columns of four different sizes and orientations; a highly absorbing, pristine, smooth hexagonal column and a highly absorbing, indented, smooth hexagonal column. Analysis of two-dimensional and one-dimensional intensity distributions and degree of linear polarisation results are given for each case and compared with results acquired through use of the Amster- dam Discrete-Dipole Approximation (ADDA) code; with good agreement observed. To the author's best knowledge, the Beam Tracer developed here is unique in its ability to handle concave particles; particles with complex structures and the man- ner in which beams are divided into sub-beams of quasi-constant intensity when propagating in an absorbing medium. One of the model's potential applications is to create a database of known particle scattering patterns, for use in aiding particle classification from images taken by the Small Ice Detector (SID) in-situ probe. An example of creating such a database for hexagonal columns is given.

No description available.

[pt] Durante muito tempo, a simulação computacional de fenômenos acústicos tem sido utilizada principalmente no projeto e estudo da acústica de ambientes. Recentemente, no entanto, podemos ver um maior interesse na utilização dessas simulações como forma de aumentar a sensação de imersão em ambientes virtuais. De forma geral, podemos dizer que um ambiente acústico virtual deve ser capaz de realizar duas tarefas: simular a propagação do som em um ambiente e ser capaz de reproduzi-lo com seu conteúdo espacial, isto é, reproduzi-lo de forma a permitir o reconhecimento da direção de propagação do som. Esta dissertação trata desses dois assuntos. São revistos os algoritmos mais comuns para o cálculo da propagação do som e, brevemente, as formas utilizadas para reproduzir áudio com conteúdo espacial. Também é apresentada a implementação de um ambiente acústico virtual, baseado nos algoritmos de beam tracing, que simula a propagação do som em ambientes bidimensionais. Como grande parte do cálculo de propagação é realizada em uma etapa de pré-processamento, o ambiente acústico virtual implementado trata apenas de fontes fixas no espaço. Os caminhos de propagação calculados são compostos de reflexões especulares e difrações do som. / [en] For a long time, computational simulation of acoustic phenomena has been used mainly in the design and study of the acoustic properties of concert and lecture halls. Recently, however, there is a growing interest in the use of such simulations in virtual environments in order to enhance users` immersion experience. Generally, we can say that a virtual acoustic environment must be able to accomplish two tasks: simulating the propagation of sound in an environment and reproducing audio with spatial content, that is, in a way that it allows the recognition of the direction of sound propagation. These tasks are the topic of the present dissertation. We begin with a revision of the most common algorithms for the simulation of sound propagation and, briefly, of the reproduction of audio with spatial content. We then present the implementation of a virtual acoustic environment, based on beam tracing algorithms, which simulates the propagation of sound waves in two-dimensional environments. As most of the computation is made in a pre-processing stage, the virtual acoustic environment implemented is appropriate only for spatially fixed sound sources. The propagation paths computed are made of specular reflections and of diffractions.

Pour commercialiser un véhicule, les constructeurs automobiles doivent se soumettre à la réglementation sur le bruit extérieur. Le règlement de la commission économique pour l'Europe, ECE R51.03, spécifie les niveaux admissibles que peut rayonner un véhicule automobile en roulage. Ce règlement est entré en vigueur depuis le 1er juillet 2016 pour remplacer l'ancien règlement ECE R51.02 (changement de méthode d’essai et sévérisation des niveaux de bruit admissibles). La diminution drastique des niveaux sonores tolérés se fait en trois étapes : passage de 74 dB (A) sous l'ancien règlement, à 68 dB (A) en 2024. Par conséquent, les constructeurs ainsi que les fournisseurs automobiles seront confrontés à un grand défi pour atteindre cet objectif. Ainsi, l'objectif de ces travaux de thèse consiste à développer une aide à la modélisation totale du bruit de passage d’un véhicule, comme le préconisent les essais réglementaires. Le but est de construire des modèles optimaux pour prévoir et évaluer avec précision le bruit que peut rayonner un véhicule en roulage plus tôt dans son cycle de développement, i.e. avant l'étape d'industrialisation. Il faut alors se placer dans la recherche d'un compromis entre précision des estimations, sensibilité aux paramètres, robustesse de la méthode et efficacité numérique. / Currently, to put a vehicle on market, car manufacturers must comply to a certification test of exterior noise. The regulation of the United Nations Economic Commission for Europe, ECE R51-03, specifies permissible levels a rolling motor vehicle can emit. This regulation is applied since July 1st, 2016, to replace the old regulation ECE R51-02 (test method change and tightening of permissible levels). The drastic reduction in noise levels will be done in 3 steps: from 74 dB (A) under the old regulation to 68 dB (A) in 2024. Therefore, manufacturers as well as their suppliers will face a great challenge to achieve this goal. The objective of this thesis is to develop an aid to the modeling of the pass-by noise of a vehicle, as called for in regulatory testing. The goal is to predict and evaluate accurately the noise emissions earlier in the vehicle development cycle, i.e. before the industrialization stage. We must then seek a trade-off between accuracy of estimates, sensitivity to parameters, robustness of the method and numerical efficiency.