506 research outputs found
hp-Version Composite Discontinuous Galerkin Methods for Elliptic Problems on Complicated Domains
In this paper we introduce the hp-version discontinuous Galerkin composite finite element method for the discretization of second-order elliptic partial differential equations. This class of methods allows for the approximation of problems posed on computational domains which may contain a huge number of local geometrical features, or microstructures. While standard numerical methods can be devised for such problems, the computational effort may be extremely high, as the minimal number of elements needed to represent the underlying domain can be very large. In contrast, the minimal dimension of the underlying composite finite element space is independent of the number of geometric features. The key idea in the construction of this latter class of methods is that the computational domain Ω is no longer resolved by the mesh; instead, the finite element basis (or shape) functions are adapted to the geometric details present in Ω. In this paper, we extend these ideas to the discontinuous Galerkin setting, based on employing the hp-version of the finite element method. Numerical experiments highlighting the practical application of the proposed numerical scheme will be presented
An agglomeration-based massively parallel non-overlapping additive Schwarz preconditioner for high-order discontinuous Galerkin methods on polytopic grids
In this article we design and analyze a class of two-level non-overlapping
additive Schwarz preconditioners for the solution of the linear system of
equations stemming from discontinuous Galerkin discretizations of second-order
elliptic partial differential equations on polytopic meshes. The preconditioner
is based on a coarse space and a non-overlapping partition of the computational
domain where local solvers are applied in parallel. In particular, the coarse
space can potentially be chosen to be non-embedded with respect to the finer
space; indeed it can be obtained from the fine grid by employing agglomeration
and edge coarsening techniques. We investigate the dependence of the condition
number of the preconditioned system with respect to the diffusion coefficient
and the discretization parameters, i.e., the mesh size and the polynomial
degree of the fine and coarse spaces. Numerical examples are presented which
confirm the theoretical bounds
Discontinuous Galerkin Methods for the Biharmonic Problem on Polygonal and Polyhedral Meshes
We introduce an -version symmetric interior penalty discontinuous
Galerkin finite element method (DGFEM) for the numerical approximation of the
biharmonic equation on general computational meshes consisting of
polygonal/polyhedral (polytopic) elements. In particular, the stability and
-version a-priori error bound are derived based on the specific choice of
the interior penalty parameters which allows for edges/faces degeneration.
Furthermore, by deriving a new inverse inequality for a special class {of}
polynomial functions (harmonic polynomials), the proposed DGFEM is proven to be
stable to incorporate very general polygonal/polyhedral elements with an
\emph{arbitrary} number of faces for polynomial basis with degree . The
key feature of the proposed method is that it employs elemental polynomial
bases of total degree , defined in the physical coordinate
system, without requiring the mapping from a given reference or canonical
frame. A series of numerical experiments are presented to demonstrate the
performance of the proposed DGFEM on general polygonal/polyhedral meshes
Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains
The numerical approximation of partial differential equations (PDEs) posed on complicated geometries, which include a large number of small geometrical features or microstructures, represents a challenging computational problem. Indeed, the use of standard mesh generators, employing simplices or tensor product elements, for example, naturally leads to very fine finite element meshes, and hence the computational effort required to numerically approximate the underlying PDE problem may be prohibitively expensive. As an alternative approach, in this article we present a review of composite/agglomerated discontinuous Galerkin finite element methods (DGFEMs) which employ general polytopic elements. Here, the elements are typically constructed as the union of standard element shapes; in this way, the minimal dimension of the underlying composite finite element space is independent of the number of geometrical features. In particular, we provide an overview of hp-version inverse estimates and approximation results for general polytopic elements, which are sharp with respect to element facet degeneration. On the basis of these results, a priori error bounds for the hp-DGFEM approximation of both second-order elliptic and first-order hyperbolic PDEs will be derived. Finally, we present numerical experiments which highlight the practical application of DGFEMs on meshes consisting of general polytopic elements
hp-Adaptive composite discontinuous Galerkin methods for elliptic eigenvalue problems on complicated domains
In this paper we develop the a posteriori error estimation of hp-adaptive discontinuous Galerkin composite finite element methods (DGFEMs) for the discretization of second-order elliptic eigenvalue problems. DGFEMs allow for the approximation of problems posed on computational domains which may contain local geometric features. The dimension of the composite finite element space is independent of the number of geometric features. This is in contrast with standard finite element methods, as the minimal number of elements needed to represent the underlying domain can be very large and so the dimension of the finite element space. Computable upper bounds on the error for both eigenvalues and eigenfunctions are derived. Numerical experiments highlighting the practical application of the proposed estimators within an automatic hp-adaptive refinement procedure will be presented
Adaptive discontinuous Galerkin methods on polytopic meshes
In this article we consider the application of discontinuous Galerkin finite element methods, defined on agglomerated meshes consisting of general polytopic elements, to the numerical approximation of partial differential equation problems posed on complicated geometries. Here, we assume that the underlying computational domain may be accurately represented by a geometry-conforming fine mesh; the resulting coarse mesh is then constructed based on employing standard graph partitioning algorithms. To improve the accuracy of the computed numerical approximation, we consider the development of goal-oriented adaptation techniques within an automatic mesh refinement strategy. In this setting, elements marked for refinement are subdivided by locally constructing finer agglomerates; should further resolution of the underlying fine mesh T_f be required, then adaptive refinement of T_f will also be undertaken. As an example of the application of these techniques, we consider the numerical approximation of the linear elasticity equations for a homogeneous isotropic material. In particular, the performance of the proposed adaptive refinement algorithm is studied for the computation of the (scaled) effective Young's modulus of a section of trabecular bone
Adaptive discontinuous Galerkin methods on polytopic meshes
In this article we consider the application of discontinuous Galerkin finite element methods, defined on agglomerated meshes consisting of general polytopic elements, to the numerical approximation of partial differential equation problems posed on complicated geometries. Here, we assume that the underlying computational domain may be accurately represented by a geometry-conforming fine mesh; the resulting coarse mesh is then constructed based on employing standard graph partitioning algorithms. To improve the accuracy of the computed numerical approximation, we consider the development of goal-oriented adaptation techniques within an automatic mesh refinement strategy. In this setting, elements marked for refinement are subdivided by locally constructing finer agglomerates; should further resolution of the underlying fine mesh T_f be required, then adaptive refinement of T_f will also be undertaken. As an example of the application of these techniques, we consider the numerical approximation of the linear elasticity equations for a homogeneous isotropic material. In particular, the performance of the proposed adaptive refinement algorithm is studied for the computation of the (scaled) effective Young's modulus of a section of trabecular bone
Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains
The numerical approximation of partial differential equations (PDEs) posed on complicated geometries, which include a large number of small geometrical features or microstructures, represents a challenging computational problem. Indeed, the use of standard mesh generators, employing simplices or tensor product elements, for example, naturally leads to very fine finite element meshes, and hence the computational effort required to numerically approximate the underlying PDE problem may be prohibitively expensive. As an alternative approach, in this article we present a review of composite/agglomerated discontinuous Galerkin finite element methods (DGFEMs) which employ general polytopic elements. Here, the elements are typically constructed as the union of standard element shapes; in this way, the minimal dimension of the underlying composite finite element space is independent of the number of geometrical features. In particular, we provide an overview of hp-version inverse estimates and approximation results for general polytopic elements, which are sharp with respect to element facet degeneration. On the basis of these results, a priori error bounds for the hp-DGFEM approximation of both second-order elliptic and first-order hyperbolic PDEs will be derived. Finally, we present numerical experiments which highlight the practical application of DGFEMs on meshes consisting of general polytopic elements
- …