Return to search

Basis Functions With Divergence Constraints for the Finite Element Method

Maxwell's equations are a system of partial differential equations of vector fields. Imposing the constitutive relations for material properties yields equations for the curl and divergence of the electric and magnetic fields. The curl and divergence equations must be solved simultaneously, which is not the same as solving three separate scalar problems in each component of the vector field.

This thesis describes a new method for solving partial differential equations of vector fields using the finite element method. New basis functions are used to solve the curl equation while allowing the divergence to be set as a constraint. The basis functions are defined on a mesh of bricks and the method is applicable for geometries that conform to a Cartesian coordinate system. The basis functions are a combination of cubic Hermite splines and second order Lagrange interpolation polynomials. The method yields a linearly independent set of constraints for the divergence, which is modelled to second order accuracy within each brick.

Mesh refinement is accomplished by dividing selected bricks into $2\times 2\times 2$ smaller bricks of equal size. The change in the node pattern at an interface where mesh refinement occurs necessitates a modified implementation of the divergence constraints as well as additional constraints for hanging nodes. The mesh can be refined to an arbitrary number of levels.

The basis functions can exactly model the discontinuity in the normal component of the field at a planar interface. The method is modified to solve problems with singularities at material boundaries that form $90^{\circ}$ edges and corners.

The primary test problem of the new basis functions is to obtain the resonant frequencies and fields of three-dimensional cavities. The new basis functions can resolve physical solutions and non-physical, spurious modes. The eigenvalues obtained with the new method are in good agreement with exact solutions and experimental values in cases where they exist. There is also good agreement with results from second-order edge elements that are obtained with the software package HFSS.

Finally, the method is modified to solve problems in cylindrical coordinates provided the domain does not contain the coordinate axis.

Identiferoai:union.ndltd.org:LACETR/oai:collectionscanada.gc.ca:OTU.1807/34847
Date19 December 2012
CreatorsPinciuc, Christopher
ContributorsDawson, Francis P.
Source SetsLibrary and Archives Canada ETDs Repository / Centre d'archives des thèses électroniques de Bibliothèque et Archives Canada
Languageen_ca
Detected LanguageEnglish
TypeThesis

Page generated in 0.0018 seconds