In this dissertation we present an analysis for the discontinuous Galerkin discretization error of multi-dimensional first-order linear symmetric and symmetrizable hyperbolic systems of conservation laws. We explicitly write the leading term of the local DG error, which is spanned by Legendre polynomials of degree p and p+1 when p-th degree polynomial spaces are used for the solution. For special hyperbolic systems, where the coefficient matrices are nonsingular, we show that the leading term of the error is spanned by (p+1)-th degree Radau polynomials. We apply these asymptotic results to observe that projections of the error are pointwise O(h^p+2)-superconvergent in some cases and establish superconvergence results for some integrals of the error. We develop an efficient implicit residual-based a posteriori error estimation scheme by solving local finite element problems to compute estimates of the leading term of the discretization error. For smooth solutions we obtain error estimates that converge to the true error under mesh refinement. We first show these results for linear symmetric systems that satisfy certain assumptions, then for general linear symmetric systems. We further generalize these results to linear symmetrizable systems by considering an equivalent symmetric formulation, which requires us to make small modifications in the error estimation procedure. We also investigate the behavior of the discretization error when the Lax-Friedrichs numerical flux is used, and we construct asymptotically exact a posteriori error estimates. While no superconvergence results can be obtained for this flux, the error estimation results can be recovered in most cases. These error estimates are used to drive h- and p-adaptive algorithms and assess the numerical accuracy of the solution. We present computational results for different fluxes and several linear and nonlinear hyperbolic systems in one, two and three dimensions to validate our theory. Examples include the wave equation, Maxwell's equations, and the acoustic equation. / Ph. D.

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.

In this thesis, we present a simple and efficient \emph{a posteriori} error estimation procedure for a discontinuous finite element method applied to scalar first-order hyperbolic problems on structured and unstructured tetrahedral meshes. We present a local error analysis to derive a discontinuous Galerkin orthogonality condition for the leading term of the discretization error and find basis functions spanning the error for several finite element spaces. We describe an implicit error estimation procedure for the leading term of the discretization error by solving a local problem on each tetrahedron. Numerical computations show that the implicit \emph{a posteriori} error estimation procedure yields accurate estimates for linear and nonlinear problems with smooth solutions. Furthermore, we show the performance of our error estimates on problems with discontinuous solutions. We investigate pointwise superconvergence properties of the discontinuous Galerkin (DG) method using enriched polynomial spaces. We study the effect of finite element spaces on the superconvergence properties of DG solutions on each class and type of tetrahedral elements. We show that, using enriched polynomial spaces, the discretization error on tetrahedral elements having one inflow face, is O(h^{p+2}) superconvergent on the three edges of the inflow face, while on elements with one inflow and one outflow faces the DG solution is O(h^{p+2}) superconvergent on the outflow face in addition to the three edges of the inflow face. Furthermore, we show that, on tetrahedral elements with two inflow faces, the DG solution is O(h^{p+2}) superconvergent on the edge shared by two of the inflow faces. On elements with two inflow and one outflow faces and on elements with three inflow faces, the DG solution is O(h^{p+2}) superconvergent on two edges of the inflow faces. We also show that using enriched polynomial spaces lead to a simpler{a posterior error estimation procedure. Finally, we extend our error analysis for the discontinuous Galerkin method applied to linear three-dimensional hyperbolic systems of conservation laws with smooth solutions. We perform a local error analysis by expanding the local error as a series and showing that its leading term is O( h^{p+1}). We further simplify the leading term and express it in terms of an optimal set of polynomials which can be used to estimate the error. / Ph. D.

We present immersed discontinuous Galerkin finite element methods for one and two dimensional acoustic wave propagation problems in inhomogeneous media where elements are allowed to be cut by the material interface. The proposed methods use the standard discontinuous Galerkin finite element formulation with polynomial approximation on elements that contain one fluid while on interface elements containing more than one fluid they use specially-built piecewise polynomial shape functions that satisfy appropriate interface jump conditions. The finite element spaces on interface elements satisfy physical interface conditions from the acoustic problem in addition to extended conditions derived from the system of partial differential equations. Additional curl-free and consistency conditions are added to generate bilinear and biquadratic piecewise shape functions for two dimensional problems. We established the existence and uniqueness of one dimensional immersed finite element shape functions and existence of two dimensional bilinear immersed finite element shape functions for the velocity. The proposed methods are tested on one dimensional problems and are extended to two dimensional problems where the problem is defined on a domain split by an interface into two different media. Our methods exhibit optimal $O(h^{p+1})$ convergence rates for one and two dimensional problems. However it is observed that one of the proposed methods is not stable for two dimensional interface problems with high contrast media such as water/air. We performed an analysis to prove that our immersed Petrov-Galerkin method is stable for interface problems with high jumps across the interface. Local time-stepping and parallel algorithms are used to speed up computation. Several realistic interface problems such as ether/glycerol, water/methyl-alcohol and water/air with a circular interface are solved to show the stability and robustness of our methods. / Ph. D.

In this work, we develop a method for solving viscoelastic fluid flows using the Navier-Stokes equations coupled with the Oldroyd-B model. We solve the Navier-Stokes equations in skew-symmetric form using the mixed finite element method, and we solve the Oldroyd-B model using the discontinuous Galerkin method. The Crank-Nicolson scheme is used for the temporal discretization of the Navier-Stokes equations in order to achieve a second-order accuracy in time, while the optimal third-order total-variation diminishing Runge-Kutta scheme is used for the temporal discretization of the Oldroyd-B equation. The overall accuracy in time is therefore limited to second-order due to the Crank-Nicolson scheme; however, a third-order Runge-Kutta scheme is implemented for greater stability over lower order Runge-Kutta schemes. We test our numerical method using the 2D cavity flow benchmark problem and compare results generated with those found in literature while discussing the influence of mesh refinement on suppressing oscillations in the polymer stress. / Master of Science / Viscoelastic fluids are a type of non-Newtonian fluid of great importance to the study of fluid flows. Such fluids exhibit both viscous and elastic behaviors. We develop a numerical method to solve the partial differential equations governing viscoelastic fluid flows using various finite element methods. Our method is then validated using previous numerical results in literature.

This research is mainly focused on the development of the adaptive Cartesian grid methods for compressibl  e flow. At first, the ghost cell method and its applications for inviscid compressible flow on adaptive tree Cartesian grid are developed. The proposed method is successfully used to evaluate various inviscid compressible flows around complex bodies. The mass conservation of the method is also studied by numerical analysis. The extension to three-dimensional flow is presented. Then, an h-adaptive Runge–Kutta discontinuous Galerkin (RKDG) method is presented in detail for the development of high accuracy numerical method under Cartesian grid. This method combined with the ghost cell immersed boundary method is also validated by well documented test problems involving both steady and unsteady compressible flows over complex bodies in a wide range of Mach numbers. In addition, in order to suppress the failure of preserving positivity of density or pressure, which may cause blow-ups of the high order numerical algorithms, a positivity-preserving limiter technique coupled with h-adaptive RKDG method is developed. Such a method has been successfully implemented to study flows with the large Mach number, strong shock/obstacle interactions and shock diffraction. The extension of the method to viscous flow under the adaptive Cartesian grid with hybrid overlapping bodyfitted grid is developed. The method is validated by benchmark problems and has been successfully implemented to study airfoil with ice accretion. Finally, based on an open source code, the detached eddy simulation (DES) is developed for massive separation flow, and it is used to perform the research on aerodynamic performance analysis over the wing with ice accretion.

In the present work we implemented parallel version of a computational fluid dynamics code. This code is based on Discontinuous Galerkin Method (DGM), which is due to its favourable properties suitable for parallelization. In the work we describe the Navier-Stokes equations and their discretization using DGM. We explain the advantages of usage of the DGM and formulate the serial algorithm. Next we focus on the parallel implementation of the algorithm and several particular issues connected to the parallelization. We present the numerical experiments showing the efficiency of the parallel code in the last chapter.

This thesis is concerned on the simulation of wet steam flow using discontinuous Galerkin method. Wet steam flow equations consist of Naviere-Stokes equations for compressible flow and Hill's equations for condensation of water vapor. The first part of this thesis describes the mathematical formulation of wet steam model and the derivation of Hill's equations. The model equations are discretized with the aid of discontinuous Galerkin method and backward difference formula which leads to implicit scheme represented by nonlinear algebraic system. This system is solved using Newton-like method. The derived scheme was implemented in program ADGFEM which is used for solving non-stationary convective-diffusive problems. The numerical results are presented in the last part of this thesis. 1

Distributed optimal control problems with transient convection dominated diffusion convection reaction equations are considered. The problem is discretized in space by using three types of discontinuous Galerkin (DG) method: symmetric interior penalty Galerkin (SIPG), nonsymmetric interior penalty Galerkin (NIPG), incomplete interior penalty Galerkin (IIPG). For time discretization, Crank-Nicolson and backward Euler methods are used. The discretize-then-optimize approach is used to obtain the finite dimensional problem. For one-dimensional unconstrained problem, Newton-Conjugate Gradient method with Armijo line-search. For two-dimensional control constrained problem, active-set method is applied. A priori error estimates are derived for full discretized optimal control problem. Numerical results for one and two-dimensional distributed optimal control problems for diffusion convection equations with boundary layers confirm the predicted orders derived by a priori error estimates.

In this thesis, we analyze an adaptive discontinuous finite element method for symmetric second order linear elliptic operators. Moreover, we obtain a fully computable convergence analysis on the broken energy seminorm in first order symmetric interior penalty discontin- uous Galerkin finite element approximations of this problem. The method is formulated on nonconforming meshes made of triangular elements with first order polynomial in two di- mension. We use an estimator which is completely free of unknown constants and provide a guaranteed numerical bound on the broken energy norm of the error. This estimator is also shown to provide a lower bound for the broken energy seminorm of the error up to a constant and higher order data oscillation terms. Consequently, the estimator yields fully reliable, quantitative error control along with efficiency. As a second problem, explicit expression for constants of the inverse inequality are given in 1D, 2D and 3D. Increasing mathematical analysis of finite element methods is motivating the inclusion of mesh dependent terms in new classes of methods for a variety of applications. Several inequalities of functional analysis are often employed in convergence proofs. Inverse estimates have been used extensively in the analysis of finite element methods. It is char- acterized as tools for the error analysis and practical design of finite element methods with terms that depend on the mesh parameter. Sharp estimates of the constants of this inequality is provided in this thesis.