75 research outputs found
Stabilization in relation to wavenumber in HDG methods
Simulation of wave propagation through complex media relies on proper
understanding of the properties of numerical methods when the wavenumber is
real and complex. Numerical methods of the Hybrid Discontinuous Galerkin (HDG)
type are considered for simulating waves that satisfy the Helmholtz and Maxwell
equations. It is shown that these methods, when wrongly used, give rise to
singular systems for complex wavenumbers. A sufficient condition on the HDG
stabilization parameter for guaranteeing unique solvability of the numerical
HDG system, both for Helmholtz and Maxwell systems, is obtained for complex
wavenumbers. For real wavenumbers, results from a dispersion analysis are
presented. An asymptotic expansion of the dispersion relation, as the number of
mesh elements per wave increase, reveal that some choices of the stabilization
parameter are better than others. To summarize the findings, there are values
of the HDG stabilization parameter that will cause the HDG method to fail for
complex wavenumbers. However, this failure is remedied if the real part of the
stabilization parameter has the opposite sign of the imaginary part of the
wavenumber. When the wavenumber is real, values of the stabilization parameter
that asymptotically minimize the HDG wavenumber errors are found on the
imaginary axis. Finally, a dispersion analysis of the mixed hybrid
Raviart-Thomas method showed that its wavenumber errors are an order smaller
than those of the HDG method
A hybridizable discontinuous Galerkin method for electromagnetics with a view on subsurface applications
Two Hybridizable Discontinuous Galerkin (HDG) schemes for the solution of
Maxwell's equations in the time domain are presented. The first method is based
on an electromagnetic diffusion equation, while the second is based on
Faraday's and Maxwell--Amp\`ere's laws. Both formulations include the diffusive
term depending on the conductivity of the medium. The three-dimensional
formulation of the electromagnetic diffusion equation in the framework of HDG
methods, the introduction of the conduction current term and the choice of the
electric field as hybrid variable in a mixed formulation are the key points of
the current study. Numerical results are provided for validation purposes and
convergence studies of spatial and temporal discretizations are carried out.
The test cases include both simulation in dielectric and conductive media
Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft
HDGlab: An Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB
This paper presents HDGlab, an open source MATLAB implementation of the hybridisable discontinuous Galerkin (HDG) method. The main goal is to provide a detailed description of both the HDG method for elliptic problems and its implementation available in HDGlab. Ultimately, this is expected to make this relatively new advanced discretisation method more accessible to the computational engineering community. HDGlab presents some features not available in other implementations of the HDG method that can be found in the free domain. First, it implements high-order polynomial shape functions up to degree nine, with both equally-spaced and Fekete nodal distributions. Second, it supports curved isoparametric simplicial elements in two and three dimensions. Third, it supports non-uniform degree polynomial approximations and it provides a flexible structure to devise degree adaptivity strategies. Finally, an interface with the open-source high-order mesh generator Gmsh is provided to facilitate its application to practical engineering problems
- …