A Discontinuous Galerkin Method for the Subjective Surfaces Problem

The work formulates and evaluates the local discontinuous Galerkin method for the subjective surfaces problem based on the curvature driven level set equation. A new mixed formulation simplifying the treatment of nonlinearities is proposed. The numerical algorithm is evaluated using several artificial and realistic test cases.     

Discontinuous Galerkin Method with Staggered Hybridization for a Class of Nonlinear Stokes Equations

In this paper, we present a discontinuous Galerkin method with staggered hybridization to discretize a class of nonlinear Stokes equations in two dimensions. The utilization of staggered hybridization is new and this approach combines the features of traditional hybridization method and staggered discontinuous Galerkin method. The main idea of our method is to use hybrid variables to impose the staggered continuity conditions instead of enforcing them in the approximation space. Therefore, our method enjoys some distinctive advantages, including mass conservation, optimal convergence and preservation of symmetry of the stress tensor. We will also show that, one can obtain superconvergent and strongly divergence-free velocity by applying a local postprocessing technique on the approximate solution. We will analyze the stability and derive a priori error estimates of the proposed scheme. The resulting nonlinear system is solved by using the Newton’s method, and some numerical results will be demonstrated to confirm the theoretical rates of convergence and superconvergence.     

Energy Norm shape A Posteriori Error Estimation for Mixed Discontinuous Galerkin Approximations of the Stokes Problem

In this paper, we develop the a posteriori error estimation of mixed discontinuous Galerkin finite element approximations of the Stokes problem. In particular, we derive computable upper bounds on the error, measured in terms of a natural (mesh-dependent) energy norm. This is done by rewriting the underlying method in a non-consistent form using appropriate lifting operators, and by employing a decomposition result for the discontinuous spaces. A series of numerical experiments highlighting the performance of the proposed a posteriori error estimator on adaptively refined meshes are presented.Paul Houston - Funded by the EPSRC (Grant GR/R76615). Thomas P. Wihler - Funded by the Swiss National Science Foundation (Grant PBEZ2-102321). This revised version was published online in July 2005 with corrected volume and issue numbers.     

A Uniform Additive Schwarz Preconditioner for High-Order Discontinuous Galerkin Approximations of Elliptic Problems

In this paper we design and analyze a uniform preconditioner for a class of high-order Discontinuous Galerkin schemes. The preconditioner is based on a space splitting involving the high-order conforming subspace and results from the interpretation of the problem as a nearly-singular problem. We show that the proposed preconditioner exhibits spectral bounds that are uniform with respect to the discretization parameters, i.e., the mesh size, the polynomial degree and the penalization coefficient. The theoretical estimates obtained are supported by numerical tests.     

Local Discontinuous Galerkin Method for Incompressible Miscible Displacement Problem in Porous Media

In this paper, we develop local discontinuous Galerkin method for the two-dimensional coupled system of incompressible miscible displacement problem. Optimal error estimates in \(L^{\infty }(0, T; L^{2})\) for concentration c, \(L^{2}(0, T; L^{2})\) for \(\nabla c\) and \(L^{\infty }(0, T; L^{2})\) for velocity \(\mathbf{u}\) are derived. The main techniques in the analysis include the treatment of the inter-element jump terms which arise from the discontinuous nature of the numerical method, the nonlinearity, and the coupling of the models. The main difficulty is how to treat the inter-element discontinuities of two independent solution variables (one from the flow equation and the other from the transport equation) at cell interfaces. Numerical experiments are shown to demonstrate the theoretical results.     

A Discontinuous Galerkin Method for Three-Dimensional Shallow Water Equations

We describe the application of a local discontinuous Galerkin method to the numerical solution of the three-dimensional shallow water equations. The shallow water equations are used to model surface water flows where the hydrostatic pressure assumption is valid. The authors recently developed a DG\linebreak method for the depth-integrated shallow water equations. The method described here is an extension of these ideas to non-depth-integrated models. The method and its implementation are discussed, followed by numerical examples on several test problems.This revised version was published online in July 2005 with corrected volume and issue numbers.     

A Hybrid Staggered Discontinuous Galerkin Method for KdV Equations

A hybrid staggered discontinuous Galerkin method is developed for the Korteweg–de Vries equation. The equation is written into a system of first order equations by introducing auxiliary variables. Two sets of finite element functions are introduced to approximate the solution and the auxiliary variables. The staggered continuity of the two finite element function spaces gives a natural flux condition and trace value on the element boundaries in the derivation of Galerkin approximation. On the other hand, to deal with the third order derivative term an hybridization idea is used and additional flux unknowns are introduced. The auxiliary variables can be eliminated in each element and the resulting algebraic system on the solution and the additional flux unknowns is solved. Stability of the semi discrete form is proven for various boundary conditions. Numerical results present the optimal order of \(L^2\)-errors of the proposed method for a given polynomial order.     

Discontinuous Galerkin Method for Investigating Ice Strength

This paper discusses the numerical modeling of various ice-strength measurement experiments, including uniaxial compression and bending, and it also compares the data obtained by field and numerical experiments. Numerical simulation is based on a dynamic system of continuum mechanics equations with ice considered as an elasto-plastic medium with brittle and crushing fracture dynamic criteria. The simulation software developed by the authors is based on the discontinuous Galerkin method and runs on high-performance systems with a distributed memory. Estimating the explicit values used by the mathematical models poses a major problem because some of them cannot be directly measured in field experiments due to the multiple interferences of physical processes. In practice, it is only possible to directly measure their total influence. However, this problem can be solved by comparing the numerical experiment with the field data. As a result of this work, the elasto-plastic ice model is verified and some missing physical properties are obtained by the numerical experiments.     

An A Posteriori Error Estimate for the Local Discontinuous Galerkin Method

In this work an a posteriori global error estimate for the Local Discontinuous Galerkin (LDG) applied to a linear second order elliptic problem is analyzed. Using a mixed formulation, an upper bound of the error in the primal variable is derived from explicit computations. Finally, a local adaptive scheme based on explicit error estimators is studied numerically using one dimensional problems.This revised version was published online in July 2005 with corrected volume and issue numbers.     

A Hybridized Discontinuous Galerkin Method with Reduced Stabilization

In this paper, we propose a hybridized discontinuous Galerkin (HDG) method with reduced stabilization for the Poisson equation. The reduce stabilization proposed here enables us to use piecewise polynomials of degree \(k\) and \(k-1\) for the approximations of element and inter-element unknowns, respectively, unlike the standard HDG methods. We provide the error estimates in the energy and \(L^2\) norms under the chunkiness condition. In the case of \(k=1\), it can be shown that the proposed method is closely related to the Crouzeix–Raviart nonconforming finite element method. Numerical results are presented to verify the validity of the proposed method.     

A high order Discontinuous Galerkin Finite Element solver for the incompressible Navier–Stokes equations

The paper presents an unsteady high order Discontinuous Galerkin (DG) solver that has been developed, verified and validated for the solution of the two-dimensional incompressible Navier–Stokes equations. A second order stiffly stable method is used to discretise the equations in time. Spatial discretisation is accomplished using a modal DG approach, in which the inter-element fluxes are approximated using the Symmetric Interior Penalty Galerkin formulation. The non-linear terms in the Navier–Stokes equations are expressed in the convective form and approximated through the Lesaint–Raviart fluxes modified for DG methods.Verification of the solver is performed for a series of test problems; purely elliptic, unsteady Stokes and full Navier–Stokes. The resulting method leads to a stable scheme for the unsteady Stokes and Navier–Stokes equations when equal order approximation is used for velocity and pressure. For the validation of the full Navier–Stokes solver, we consider unsteady laminar flow past a square cylinder at a Reynolds number of 100 (unsteady wake). The DG solver shows favourably comparisons to experimental data and a continuous Spectral code.     

Error Analysis for a Hybridizable Discontinuous Galerkin Method for the Helmholtz Equation

Finite element methods for acoustic wave propagation problems at higher frequency result in very large matrices due to the need to resolve the wave. This problem is made worse by discontinuous Galerkin methods that typically have more degrees of freedom than similar conforming methods. However hybridizable discontinuous Galerkin methods offer an attractive alternative because degrees of freedom in each triangle can be cheaply removed from the global computation and the method reduces to solving only for degrees of freedom on the skeleton of the mesh. In this paper we derive new error estimates for a hybridizable discontinuous Galerkin scheme applied to the Helmholtz equation. We also provide extensive numerical results that probe the optimality of these results. An interesting observation is that, after eliminating the internal element degrees of freedom, the condition number of the condensed hybridized system is seen to be almost independent of the wave number.     

A Fully-Discrete Local Discontinuous Galerkin Method for Convection-Dominated Sobolev Equation

In this paper we shall present, for the convection-dominated Sobolev equations, the fully-discrete numerical scheme based on the local discontinuous Galerkin (LDG) finite element method and the third-order explicitly total variation diminishing Runge-Kutta (TVDRK3) time marching. A priori error estimate is obtained for any piecewise polynomials of degree at most k≥1, under the general spatial-temporal restriction. The bounded constant in error estimate is independent of the reciprocal of the diffusion and dispersion coefficients, after removing the effect of smoothness of the exact solution. Finally some numerical results are given to verify the presented conclusion.     

A Sequential Discontinuous Galerkin Method for the Coupling of Flow and Geomechanics

A decoupled and sequential numerical method is proposed and analyzed for solving the linear poroelasticity equations. Unlike other splitting approaches, this method is not iterative, which results in a speed-up of the computational time. The interior penalty discontinuous Galerkin method is employed for the spatial discretization and is combined with the backward Euler method for the time discretization. We provide a convergence analysis of the scheme along with numerical results that confirm the theoretical results.