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.

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)

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

This thesis studies the numerical solution of non-linear convection-diffusion problems using the space- time discontinuous Galerkin method, which perfectly suits the space as well as time local adaptation. We aim to develop a posteriori error estimates reflecting the spatial, temporal, and algebraic errors. These estimates are based on the measurement of the residuals in dual norms. We derive these estimates and numerically verify their properties. Finally, we derive an adaptive algorithm and apply it to the numerical simulation of non-stationary viscous compressible flows. Powered by TCPDF (www.tcpdf.org)

No description available.