Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA robust discontinuous Galerkin scheme on anisotropic meshes
114
0
0.0
(
0
)
Added by
Takuya Tsuchiya
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Discontinuous Galerkin (DG) methods are extensions of the usual Galerkin finite element methods. Although there are vast amount of studies on DG methods, most of them have assumed shape-regularity conditions on meshes for both theoretical error analysis and practical computations. In this paper, we present a new symmetric interior penalty DG scheme with a modified penalty term. We show that, without imposing the shape-regularity condition on the meshes, the new DG scheme inherits all of the good properties of standard DG methods, and is thus robust on anisotropic meshes. Numerical experiments confirm the theoretical error estimates obtained.
rate research
Read More
In this paper we consider a level set reinitialization technique based on a high-order, local discontinuous Galerkin method on unstructured triangular meshes. A finite volume based subcell stabilization is used to improve the nonlinear stability of the method. Instead of the standard hyperbolic level set reinitialization, the flow of time Eikonal equation is discretized to construct an approximate signed distance function. Using the Eikonal equation removes the regularization parameter in the standard approach which allows more predictable behavior and faster convergence speeds around the interface. This makes our approach very efficient especially for banded level set formulations. A set of numerical experiments including both smooth and non-smooth interfaces indicate that the method experimentally achieves design order accuracy.
We present a discontinuous Galerkin internal-penalty scheme that is applicable to a large class of linear and non-linear elliptic partial differential equations. The scheme constitutes the foundation of the elliptic solver for the SpECTRE numerical relativity code. As such it can accommodate (but is not limited to) elliptic problems in linear elasticity, general relativity and hydrodynamics, including problems formulated on a curved manifold. We provide practical instructions that make the scheme functional in a production code, such as instructions for imposing a range of boundary conditions, for implementing the scheme on curved and non-conforming meshes and for ensuring the scheme is compact and symmetric so it may be solved more efficiently. We report on the accuracy of the scheme for a suite of numerical test problems.
In this article, using the weighted discrete least-squares, we propose a patch reconstruction finite element space with only one degree of freedom per element. As the approximation space, it is applied to the discontinuous Galerkin methods with the upwind scheme for the steady-state convection-diffusion-reaction problems over polytopic meshes. The optimal error estimates are provided in both diffusion-dominated and convection-dominated regimes. Furthermore, several numerical experiments are presented to verify the theoretical error estimates, and to well approximate boundary layers and/or internal layers.
We introduce a new stabilization for discontinuous Galerkin methods for the Poisson problem on polygonal meshes, which induces optimal convergence rates in the polynomial approximation degree $p$. In the setting of [S. Bertoluzza and D. Prada, A polygonal discontinuous Galerkin method with minus one stabilization, ESAIM Math. Mod. Numer. Anal. (DOI: 10.1051/m2an/2020059)], the stabilization is obtained by penalizing, in each mesh element $K$, a residual in the norm of the dual of $H^1(K)$. This negative norm is algebraically realized via the introduction of new auxiliary spaces. We carry out a $p$-explicit stability and error analysis, proving $p$-robustness of the overall method. The theoretical findings are demonstrated in a series of numerical experiments.
A new scheme for communication between overset grids using subcells and Weighted Essentially Non Oscillatory (WENO) reconstruction for two-dimensional problems has been proposed. The effectiveness of this procedure is demonstrated using the discontinuous Galerkin method (DGM). This scheme uses WENO reconstruction using cell averages by dividing the immediate neighbors into subcells to find the degrees of freedom in cells near the overset interface. This also has the added advantage that it also works as a limiter if a discontinuity passes through the overset interface. Accuracy tests to demonstrate the maintenance of higher order are provided. Results containing shocks are also provided to demonstrate the limiter aspect of the data communication procedure.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
