A Discontinuous Galerkin Time Domain method (DGTD), using a fourth order Runge-Kutta time-stepping of Maxwell's equations, was applied to the simulation of the optical properties of dielectric particles in two-dimensional (2-D) geometry. As examples of the numerical implementation of this method, the single-scattering properties of 2D circular and hexagonal particles are presented. In the case of circular particles, the scattering phase matrix was computed using the DGTD method and compared with the exact solution. For hexagonal particles, the DGTD method was used to compute single-scattering properties of randomly oriented 2-D hexagonal ice crystals, and results were compared with those calculated using a geometric optics method. Both shortwave (visible) and longwave (infrared) cases are considered, with particle size parameters 50 and 100. Ice in shortwave and longwave cases is absorptive and non-absorptive, respectively. The comparisons between DG solutions and the exact solutions in computing the optical properties of circular ice crystals reveal the applicability of the DG method to calculations of both absorptive and non-absorptive particles. In the hexagonal case scattering results are also presented as a function of both incident and scattering angles, revealing structure apparently not reported before. Using the geometric optics method we are able to interpret this structure in terms of contributions from varying numbers of internal reflections within the crystal.

The discontinuous Galerkin method (DGM) is applied to a number of problems in computational electromagnetics. This is achieved by obtaining numerical solutions to Maxwell's equations using the DGM. The aim of these simulations is to highlight the strengths of the method while showing its resilience in handling problems other schemes may not be able to accurately model. Although no method will ever be the best choice for every problem in electromagnetics, the discontinuous Galerkin method is able to accurately approximate any problem, although the computational costs can make the scheme impractical for some. Like other time domain schemes, the DGM becomes inefficient on large domains where the solution contains small wavelengths. We demonstrate that all of the different types of boundary conditions in electromagnetic wave propagation can be implemented into the DGM. Reflection and transmission boundaries fit easily into the framework, whereas perfect absorption requires a more advanced technique known as the perfectly matched layer. We begin by simulating mirrors with several different geometries, and analyze how the DGM method performs, and how it offers a more complete evaluation of the behavior in this problem than some other methods. Since Maxwell's equations describe the macroscopic features of electromagnetics, our simulations are able to capture the wave features of electromagnetics, such as interference and diffraction. We demonstrate this by accurately modelling Young's double slit experiment, a classic experiment which features well understood interference and diffraction phenomena. We also extend the basic electromagnetic wave propagation simulations to include situations where the waves travel into new media. The formulation of the DGM for Maxwell's equations allows the numerical solutions to accurately resolve the features at the interface of two media as predicted by the Fresnel coefficients. This allows the DGM to model lenses and other sources of refraction. We predict that the DGM will become an increasingly valuable method for computational electromagnetics because of its wide range of applicability as well as the lack of undesirable features in the numerical solutions. Furthermore, the only limiting factor for applying DGM, its computational cost, will become less influential as computing power continues to increase, allowing us to apply the DGM to an increasing set of applications.

Le but de cette thèse est l'établissement de formulations discontinues de Galerkin pour trois problèmes dérivant des équations de Maxwell.<br />En outre nous présentons de nouvelles inégalités de type Poincré Friedrichs et un couplage entre une méthode discontinue de Galerkin et une méthode intégrale de frontière.

The discontinuous Galerkin method (DGM) is applied to a number of problems in computational electromagnetics. This is achieved by obtaining numerical solutions to Maxwell's equations using the DGM. The aim of these simulations is to highlight the strengths of the method while showing its resilience in handling problems other schemes may not be able to accurately model. Although no method will ever be the best choice for every problem in electromagnetics, the discontinuous Galerkin method is able to accurately approximate any problem, although the computational costs can make the scheme impractical for some. Like other time domain schemes, the DGM becomes inefficient on large domains where the solution contains small wavelengths. We demonstrate that all of the different types of boundary conditions in electromagnetic wave propagation can be implemented into the DGM. Reflection and transmission boundaries fit easily into the framework, whereas perfect absorption requires a more advanced technique known as the perfectly matched layer. We begin by simulating mirrors with several different geometries, and analyze how the DGM method performs, and how it offers a more complete evaluation of the behavior in this problem than some other methods. Since Maxwell's equations describe the macroscopic features of electromagnetics, our simulations are able to capture the wave features of electromagnetics, such as interference and diffraction. We demonstrate this by accurately modelling Young's double slit experiment, a classic experiment which features well understood interference and diffraction phenomena. We also extend the basic electromagnetic wave propagation simulations to include situations where the waves travel into new media. The formulation of the DGM for Maxwell's equations allows the numerical solutions to accurately resolve the features at the interface of two media as predicted by the Fresnel coefficients. This allows the DGM to model lenses and other sources of refraction. We predict that the DGM will become an increasingly valuable method for computational electromagnetics because of its wide range of applicability as well as the lack of undesirable features in the numerical solutions. Furthermore, the only limiting factor for applying DGM, its computational cost, will become less influential as computing power continues to increase, allowing us to apply the DGM to an increasing set of applications.

Shock capturing has been a challenge for computational fluid dynamicists over the years. This article deals with discontinuous Galerkin method to solve the hyperbolic equations in which solutions may develop discontinuities in finite time. The high order discontinuous Galerkin method combining the basis of finite volume and finite element methods has shown a lot of attractive features for a wide range of applications. Various techniques proposed in the literature to deal with discontinuities basically reduce the order of interpolation in the region around these discontinuities. The accuracy of the scheme therefore may be degraded in the vicinity of the shock. The proposed method resolves the discontinuities presented in the solution by applying viscosity into the shock-containing elements. The discontinuity is spread over a distance and is well approximated in the space of interpolation functions. The technique of adding viscosity to the system and the indicator based on the expansion coefficients of the solution are presented. A number of numerical examples in one and two dimensions is carried out to show the capability of the scheme for shock capturing. / Singapore-MIT Alliance (SMA)

In this dissertation we discuss bilinear immersed finite elements (IFE) for solving interface problems. The related research works can be categorized into three aspects: (1) the construction of the bilinear immersed finite element spaces; (2) numerical methods based on these IFE spaces for solving interface problems; and (3) the corresponding error analysis. All of these together form a solid foundation for the bilinear IFEs. The research on immersed finite elements is motivated by many real world applications, in which a simulation domain is often formed by several materials separated from each other by curves or surfaces while a mesh independent of interface instead of a body-fitting mesh is preferred. The bilinear IFE spaces are nonconforming finite element spaces and the mesh can be independent of interface. The error estimates for the interpolation of a Sobolev function in a bilinear IFE space indicate that this space has the usual approximation capability expected from bilinear polynomials, which is <i>O</i>(<i>h</i>²) in <i>L</i>² norm and <i>O</i>(<i>h</i>) in <i>H</i>¹ norm. Then the immersed spaces are applied in Galerkin, finite volume element (FVE) and discontinuous Galerkin (DG) methods for solving interface problems. Numerical examples show that these methods based on the bilinear IFE spaces have the same optimal convergence rates as those based on the standard bilinear finite element for solutions with certain smoothness. For the symmetric selective immersed discontinuous Galerkin method based on bilinear IFE, we have established its optimal convergence rate. For the Galerkin method based on bilinear IFE, we have also established its convergence. One of the important advantages of the discontinuous Galerkin method is its flexibility for both <i>p</i> and <i>h</i> mesh refinement. Because IFEs can use a mesh independent of interface, such as a structured mesh, the combination of a DG method and IFEs allows a flexible adaptive mesh independent of interface to be used for solving interface problems. That is, a mesh independent of interface can be refined wherever needed, such as around the interface and the singular source. We also develop an efficient selective immersed discontinuous Galerkin method. It uses the sophisticated discontinuous Galerkin formulation only around the locations needed, but uses the simpler Galerkin formulation everywhere else. This selective formulation leads to an algebraic system with far less unknowns than the immersed DG method without scarifying the accuracy; hence it is far more efficient than the conventional discontinuous Galerkin formulations. / Ph. D.

This work deals with numerical modeling of electromagnetic induction in 3D environment with heterogeneous conductivity. We develop a program to solve Maxwell's equations in quasistatic approximation by using Continuous and Discontinuous Finite Elements. Their implementation in the numerical library deal.ii is discussed. The obtained numerical results are compared with each other and also with a quasianalytic solution for an environment with 1D heterogeneous conductivity. We discuss different numerical methods, limits of our code for practical use and possible future enhancements.

In this thesis, we analyze theoretical properties of the discontinuous Galerkin method (DGM) and propose novel approaches to implementation with the aim to increase its efficiency. First, we derive explicit expressions for the eigenvalues (spectrum) of the discontinuous Galerkin spatial discretization applied to the linear advection equation. We show that the eigenvalues are related to the subdiagonal [p/p+1] Pade approximation of exp(-z) when the p-th degree basis functions are used. Then, we extend the analysis to nonuniform meshes where both the size of elements and the composition of the mesh influence the spectrum. We show that the spectrum depends on the ratio of the size of the largest to the smallest cell as well as the number of cells of different types. We find that the spectrum grows linearly as a function of the proportion of small cells present in the mesh when the size of small cells is greater than some critical value. When the smallest cells are smaller than this critical value, the corresponding eigenvalues lie outside of the main spectral curve. Numerical examples on nonuniform meshes are presented to show the improvement on the time step restriction. In particular, this result can be used to improve the time step restriction on Cartesian grids. Finally, we present a discontinuous Galerkin method for solutions of the Euler equations on Cartesian grids with embedded geometries. Cutting an embedded geometry out of the Cartesian grid creates cut cells, which are difficult to deal with for two reasons. One is the restrictive CFL number and the other is the integration on irregularly shaped cells. We use explicit time integration employing cell merging to avoid restrictively small time steps. We provide an algorithm for splitting complex cells into triangles and use standard quadrature rules on these for numerical integration. To avoid the loss of accuracy due to straight sided grids, we employ the curvature boundary conditions. We show that the proposed method is robust and high-order accurate.

In this thesis, we analyze theoretical properties of the discontinuous Galerkin method (DGM) and propose novel approaches to implementation with the aim to increase its efficiency. First, we derive explicit expressions for the eigenvalues (spectrum) of the discontinuous Galerkin spatial discretization applied to the linear advection equation. We show that the eigenvalues are related to the subdiagonal [p/p+1] Pade approximation of exp(-z) when the p-th degree basis functions are used. Then, we extend the analysis to nonuniform meshes where both the size of elements and the composition of the mesh influence the spectrum. We show that the spectrum depends on the ratio of the size of the largest to the smallest cell as well as the number of cells of different types. We find that the spectrum grows linearly as a function of the proportion of small cells present in the mesh when the size of small cells is greater than some critical value. When the smallest cells are smaller than this critical value, the corresponding eigenvalues lie outside of the main spectral curve. Numerical examples on nonuniform meshes are presented to show the improvement on the time step restriction. In particular, this result can be used to improve the time step restriction on Cartesian grids. Finally, we present a discontinuous Galerkin method for solutions of the Euler equations on Cartesian grids with embedded geometries. Cutting an embedded geometry out of the Cartesian grid creates cut cells, which are difficult to deal with for two reasons. One is the restrictive CFL number and the other is the integration on irregularly shaped cells. We use explicit time integration employing cell merging to avoid restrictively small time steps. We provide an algorithm for splitting complex cells into triangles and use standard quadrature rules on these for numerical integration. To avoid the loss of accuracy due to straight sided grids, we employ the curvature boundary conditions. We show that the proposed method is robust and high-order accurate.

A look into an implementation of turbulence model into the ADGFEM code for viscous flow. Discretization, theory background and development of the method will be carried during this thesis. Also some numerical examples of the application of the code will be provided. 1