In literature Trefftz method normally has geometric (exponential) convergence. Recently many scholars have found that spectral method in some cases can converge faster than exponential, which is called super-geometric convergence. Since Trefftz method can be regarded as a kind of spectral method, we expect it might possess super-geometric convergence too. In this thesis, we classify all types of super-geometric convergence and compare their speeds. We develop a method to decide the convergent type of given error data. Finally we can observe in many numerical experiments the super-geometric convergence of Trefftz method to solve Helmholtz boundary value problems.

In the thesis, we propose a newly-developed method called the method of Connected Local Fields (CLF) to analyze opto-electromagnetic passive devices. The method of CLF somewhat resembles a hybrid between the finite difference and pseudo-spectral methods. For opto-electromagnetic passive devices, our primary concern is their steady state behavior, or narrow-band characteristics, so we use a frequency-domain method, in which the system is governed by the Helmholtz equation. The essence of CLF is to use the intrinsic general solution of the Helmholtz equation to expand the local fields on the compact stencil. The original equation can then be transformed into the discretized form called LFE-9 (in 2-D case), and the intrinsic reconstruction formulae describing each overlapping local region can be obtained. Further, we present rigorous analysis of the numerical dispersion equation of LFE-9, by means of first-order approximation, and acquire the closed-form formula of the relative numerical dispersion error. We are thereby able to grasp the tangible influences brought both by the sampling density as well as the propagation direction of plane wave on dispersion error. In our dispersion analysis, we find that the LFE-9 formulation achieves the sixth-order accuracy: the theoretical highest order for discretizing elliptic partial differential equations on a compact nine-point stencil. Additionally, the relative dispersion error of LFE-9 is less than 1%, given that sampling density greater than 2.1 points per wavelength. At this point, the sampling density is nearing that of the Nyquist-Shannon sampling limit, and therefore computational efforts can be significantly reduced.

For Laplace's eigenvalue problems, this thesis presents new algorithms of the Trefftz method (i.e. the boundary approximation method), which solve the Helmholtz equation and then use a iteration process to yield approximate eigenvalues and eigenfunctions. The new iteration method has superlinear convergence rates and gives a better performance in numerical testing, compared with the other popular methods of rootfinding. Moreover, piecewise particular solutions are used for a basic model of eigenvalue problems on the unit square with the Dirichlet condition. Numerical experiments are also conducted for the eigenvalue problems with singularities. Our new algorithms using piecewise particular solutions are well suited to seek very accurate solutions of eigenvalue problems, in particular those with multiple singularities, interfaces and those on unbounded domains. Using piecewise particular solutions has also the advantage to solve complicated problems because uniform particular solutions may not always exist for the entire solution domain.

In the thesis, the error and stability analysis is made for the 2D Helmholtz equation by the method of fundamental solutions (MFS) using both Bessel and Neumann functions. The bounds of errors in bounded simply-connected domains are derived, while the bounds of condition number are derived only for disk domains. The MFS using Bessel functions is more efficient than the MFS using Neumann functions. Interestingly, for the MFS using Bessel functions, the radius R of the source points is not necessarily larger than the maximal radius r_max of the solution domain. This is against the traditional condition: r_max < R for MFS. Numerical experiments are carried out to support the analysis and conclusions made.

Dans cette thèse, nous présentons des méthodes asymptotiques <br />mathématiquement justifiées permettant de connaître les champs <br />électromagnétiques dans des milieux à couches minces hétérogènes. <br />La motivation de ce travail est le calcul du champ électrique dans des <br />cellules biologiques composées d'un cytoplasme conducteur entouré <br />d'une fine membrane très isolante. <br />Nous remplaçons la membrane, lorsque son épaisseur est infiniment <br />petite, par des conditions de transmission ou des conditions aux <br />limites appropriées et nous estimons l'erreur commise par ces <br />approximations.<br /> Pour les basses fréquences, nous considérons l'équation quasistatique<br />donnant le potentiel dont dérive le champ. A l'aide d'un <br />calcul en géométrie circulaire nous obtenons les expressions explicites<br /> du potentiel et nous en déduisons les asymptotiques du champ <br />électrique, en fonction de l'épaisseur de la couche mince, avec des <br />estimations de l'erreur. Nous estimons ensuite la différence entre le <br />champ réel et le champ statique. Puis nous généralisons notre <br />développement asymptotique à une géométrie quelconque. <br /> La deuxième partie de cette thèse traite des moyennes fréquences : <br />nous donnons le développement asymptotique de la solution de <br />l'équation de Helmholtz lorsque l'épaisseur de la membrane tend vers <br />0. Tous ces précédents résultats sont illustrés par des calculs par <br />éléments finis.<br /> Enfin, pour les hautes fréquences, nous construisons une condition <br />d'impédance pseudodifférentielle permettant de concentrer l'effet de <br />la couche sur son bord intérieur. Nous concluons cette thèse par un <br />problème de diffraction à haute fréquence d'une onde incidente par <br />un disque de petite taille. A l'aide d'une analyse pseudodifférentielle, <br />nous bornons la norme de la trace du champ diffracté à distance fixe <br />de l'inhomogénéité en fonction de la taille de l'objet et de l'onde <br />incidente.

Thesis (Ph. D.)--University of Washington, 2003. / Vita. Includes bibliographical references (p. 134-137).

No description available.

Thesis (MEng) -- University of Stellenbosch, 1991. / ENGLISH ABSTRACT: A two-dimensional Coupled Element Method (CEM) for solving electromagnetic scattering problems involving lossy, inhomogeneous, arbitrarily shaped cylinders, was investigated and implemented. The CEM uses the Finite Element Method (FEM) to approximate the fields in and around the scatterer and the Boundary Element Method (BEM) to approximate the far-field values. The basic CEM theory is explained using the special, static electric field problem involving the solution of Laplace's equation. This theory is expanded to incorporate scattering problems, involving the solution of the Helmholtz equation. This is done for linear as well as quadratic elements. Some of the important algorithms used to implement the CEM theory are discussed. Analytical solutions for a round, homogeneous- and one layer coated PC cylinder are discussed and obtained. The materials used in these analytical solutions can be lossy as well as chiral. The CEM is validated by comparing near- and far-field results to the analytical solution. A comparison between linear and quadratic elements is also made. The theory of the CEM is further expanded to incorporate scattering from chiral media / AFRIKAANSE OPSOMMING: 'n Gekoppelde Element Metode (GEM) wat elektromagnetiese weerkaatsingsprobleme, van verlieserige, nie-homogene, arbitrere voorwerpe kan oplos, is ondersoek en geimplimenteer. Die GEM gebruik die Eindige Element Metode (EEM) om die velde in en om die voorwerp te benader. 'n Grenselementmetode word gebruik om die vervelde te benader. Die basiese teorie van die GEM word verduidelik deur die toepassing daarvan op die spesiale geval van 'n statiese elektriese veld- probleem. Hierdie probleem verlang die oplossing van Laplace se vergelyking. Die teorie word uitgebrei om weerkaatsingsprobleme te kan hanteer. Die weerkaatsingsprobleme verlang die oplossing van 'n Helmholtz-vergelyking. Hierdie teorie word ontwikkel vir lineere sowel as kwadratiese elemente. Van die belangrike algoritmes wat gebruik is om die GEM-teorie te implimenteer, word bespreek. Analietise oplossings vir ronde, homogene en eenlaag bedekte perfek geleidende silinders word bespreek en verkry. Die material wat in die oplossings gebruik word, kan verlieserig of kiraal wees. Die GEM word bekragtig deur naby- en verveld resultate te vergelyk met ooreenkomstige aitalitiese oplossings. Die lineere en kwadratiese element- resultate word ook met mekaar vergelyk. Die GEM-teorie is verder uitgebrei sodat weerkaatsing vanaf kirale materiale ook hanteer kan word.

No description available.

Dans cette thèse, on s'intéresse à la propagation d'ondes en milieu fortement hétérogène modélisée par l'équation d'Helmholtz. Les méthodes numériques permettant de résoudre ce problème souffrent de dispersion numérique, en particulier à haute fréquence. Ce phénomène, appelé "effet de pollution", est largement analysé dans la littérature quand le milieu de propagation est homogène et l'utilisation de "méthodes d'ordre élevé" est souvent proposée pour minimiser ce problème. Dans ce travail, on s'intéresse à un milieu de propagation hétérogène, cas pour lequel on dispose de moins de connaissances. On propose d'adapter des méthodes éléments finis d'ordre élevé pour résoudre l'équation d'Helmholtz en milieu hétérogène, afin de réduire l'effet de pollution. Les méthodes d'ordre élevé étant généralement basées sur des maillages "larges", une stratégie multi-échelle originale est développée afin de prendre en compte des hétérogénéités de petite échelle. La convergence de la méthode est démontrée. En particulier, on montre que la méthode est robuste vis-a-vis de l'effet de pollution. D'autre part, on applique la méthode a plusieurs cas-tests numériques. On s'intéresse d'abord à des problèmes académiques, qui permettent de valider la théorie de convergence développée. On considère ensuite des cas-tests "industriels" appliqués à la Géophysique. Ces derniers nous permettent de conclure que la méthode multi-échelle proposée est plus performante que les éléments finis "classiques" et que des problèmes 3D réalistes peuvent être considérés. / The main objective of this work is the design of an efficient numerical strategy to solve the Helmholtz equation in highly heterogeneous media. We propose a methodology based on coarse meshes and high order polynomials together with a special quadrature scheme to take into account fine scale heterogeneities. The idea behind this choice is that high order polynomials are known to be robust with respect to the pollution effect and therefore, efficient to solve wave problems in homogeneous media. In this work, we are able to extend so-called "asymptotic error-estimate" derived for problems homogeneous media to the case of heterogeneous media. These results are of particular interest because they show that high order polynomials bring more robustness with respect to the pollution effect even if the solution is not regular, because of the fine scale heterogeneities. We propose special quadrature schemes to take int account fine scale heterogeneities. These schemes can also be seen as an approximation of the medium parameters. If we denote by h the finite-element mesh step and by e the approximation level of the medium parameters, we are able to show a convergence theorem which is explicit in terms of h, e and f, where f is the frequency. The main theoretical results are further validated through numerical experiments. 2D and 3D geophysica benchmarks have been considered. First, these experiments confirm that high-order finite-elements are more efficient to approximate the solution if they are coupled with our multiscale strategy. This is in agreement with our results about the pollution effect. Furthermore, we have carried out benchmarks in terms of computational time and memory requirements for 3D problems. We conclude that our multiscale methodology is able to greatly reduce the computational burden compared to the standard finite-element method