Questions tagged [discontinuous-galerkin]

Questions about analysis, implementation or application of Galerkin methods for partial differential equations using piecewise functions that are not globally continuous (and hence require surface terms on element boundaries in addition to the usual volume terms occurring in finite element methods).

Filter by
Sorted by
Tagged with
19
votes
2answers
4k views

Discontinuous Galerkin: Nodal vs Modal advantages and disadvantages

There are two general approaches to representing solutions in the discontinuous galerkin method: nodal and modal. Modal: Solutions are represented by sums of modal coefficients multiplied by a set of ...
15
votes
1answer
2k views

Visualizing discontinuous Galerkin/finite element data

I would like to visualize simulation results, obtained using the discontinuous Galerkin (DG) approach, within ParaView. Similarly to finite volume methods, the problem domain is divided into cube-...
14
votes
3answers
1k views

Role of the numerical flux in DG-FEM

I am learning the theory behind DG-FEM methods using the Hesthaven/Warburton book and I am a bit confused about the role of the 'numerical flux.' I apologize if this is a basic question, but I have ...
13
votes
0answers
425 views

Sequential approach to solving coupled PDEs

I'm dealing with a coupled system of three transient, non-linear convection-diffusion equations. Let's just say to simplify the problem that they take the following form: $$ -\nabla\cdot(D_{1}(u_{2},...
11
votes
1answer
461 views

DG local equation, how to interpret mean-averaged test function

In the paper http://www.sciencedirect.com/science/article/pii/S0045782509003521, an HDG element-local equation is described on page 584 equation (4), with one of the equations taking the following ...
9
votes
5answers
1k views

Is discontinuous Galerkin really any more parallelizable than continuous Galerkin?

I've always heard that easy parallelization was one of the advantages of DG methods, but I don't really see why any of those reasons don't also apply to continuous Galerkin.
9
votes
1answer
665 views

CFL condition in Discontinuous Galerkin schemes

I have implemented an ADER-Discontinuous Galerkin scheme for the resolution of linear systems of conservation laws of the type of $\partial_t U + A \partial_x U + B \partial_y U=0 $ and observed that ...
8
votes
1answer
148 views

Testing and visualizing large index arrays

I will be implementing nodal discontinuous Galerkin method soon, and having done this before I know the basic indexing arrays I will need to compute, given a mesh and polynomial data. The problem I ...
7
votes
3answers
225 views

Has a uniform estimate in k of the inf-sup constant for hp-DG methods for the Stokes problem been established?

In Theorem 6.2 of their 2003 paper on "Mixed hp-DGFEM for incompressible flows. SIAM J. Numer. Anal., 40(6), 2171–2194", D. Schötzau, Ch. Schwab, and A. Toselli prove a bound of the $\inf$-$\sup$ ...
7
votes
1answer
248 views

Size of jump for piecewise discontinuous approximations

If one has a sufficiently smooth function $u$ that is approximated by a piecewise constant function $u_h=\Pi^0_h u$ on a mesh of cell size $h$ (where $\Pi^0_h$ is the $L_2$ projection onto the ...
7
votes
2answers
134 views

Convergence of adaptive finite elements with inexact solves

I'm working on some adaptive discontinuous Galerkin codes for time harmonic wave propagation, currently just Helmholtz, but will be branching out once I have a working prototype in this case. There ...
7
votes
1answer
302 views

Solving the Advection Equation with Forcing using the Discontinuous Galerkin Method

I've been learning about the Discontinous Galkerin Method by reading the book by Hesthaven and Warburton and have ran into a problem with the advection equation with forcing $u_t + u_x = g(x,t)$ ...
6
votes
2answers
807 views

Local inversion of small matrices on GPUs?

I don't know much about GPU computing at the moment, so please pardon the simple question. Can one invert local matrices in parallel on the GPU? CUBLAS doesn't seem to support factorization, and most ...
6
votes
1answer
273 views

Velocity-Stress formulation of Elastodynamics/Wave Equation for beginner

I'm used to displacement forumlation of elastic wave equation: $$ \nabla \cdot \sigma (u) + F = \rho \ddot{ u } $$ where $u$ is the primary variable. Recenty I started experimenting with DG and in ...
6
votes
1answer
132 views

$L^\infty$ stability property of an ODE

Suppose we have the initial-value problem on $(0,L)$: $$ \frac{d u(x)}{d x} = f(x) u(x),\, \qquad x\in\Omega,\,~~ u(0) = u_0, $$ I am reading a claim that says if we multiply the ODE by $u$ and ...
6
votes
1answer
587 views

How to project a vector into the H(div) space (in the context of finite elements)?

Say I have a simple elliptic PDE: $$ -\nabla\cdot(K\nabla p) = f \;\;\;\text{in}\;\Omega $$ with the appropriate boundary conditions. I solve for $p$ using a FEM (a discontinuous Galerkin method to ...
5
votes
1answer
274 views

Stability of hyperbolic PDE and DG-FEM

In the book of Hesthaven and Warburton on discontinuous Galerlkin methods in example 2.3 (regarding solutions of the wave equation), the authors regard the following PDE: $$\frac{\partial u }{\...
5
votes
2answers
678 views

Absorbing boundary conditions for acoustics in Discontinuous Galerkin

Note: I'm trying to implement a Discontinuous Galerkin method, as kind of a way to learn about these things. As of now, I've taken the acoustic wave equation $c^2 \nabla \cdot \nabla u(x,t) - \frac{\...
5
votes
1answer
460 views

Choosing the penalty for Discontinuous Galerkin

So I am attempting to solve a 3D poisson problem with Discontinuous Galerkin (interior penalty method). The weak form (written in FEniCS) is as following: ...
5
votes
3answers
265 views

How to efficiently structure simulation data in memory for cells with varying degrees of freedom?

For a discontinuous Galerkin-based simulation I need to store cell-based simulation data in memory. Since the order of the polynomial approximation $N_p$ may vary between cells, I wonder what the most ...
5
votes
1answer
336 views

Gauss-Lobatto quadrature and nodal points for FEM

By using the Legendre-Gauss-Lobatto (LGL) quadrature formula (QF) and LGL nodal points one achives a diagonal mass-matrix for finite element problems. (More specifically, the spectral element method.) ...
5
votes
1answer
295 views

Proof of CFL condition for RKDG scheme

The cfl condition for linear advection equation $$ u_t + a u_x = 0 $$ using a DG method of degree $k$ polynomials, upwind flux and an RK scheme of $k+1$ stage/accuracy is stated to be $\frac{1}{2k+1}$...
5
votes
1answer
325 views

Are additional penalty terms necessary to solve elliptic PDE's with DG-FEM?

After a cursory glance at several references to discontinuous galerkin finite element methods for the elliptic poisson PDE, i notice that all of them emphasize using penalty methods where an ...
5
votes
1answer
298 views

Numerical quadrature in Discontinuous Galerkin

I would like to know which is the best way to integrate numerically Legendre polynomials. I am building up a Discontinuous Galerkin CFD code for which Legendre polynomials are used as basis functions ...
5
votes
1answer
995 views

Which finite-element framework allows the easy implementation of HDG methods?

I have tinkered around with the implementation of hybridizable Discontinuous Galerkin (HDG) methods in DUNE-PDELab and DUNE-FEM for my PhD but since then I have not had the time to work on these codes ...
5
votes
0answers
173 views

How to compute numerical fluxes in the local discontinuous galerkin method for poisson equation 1D

Some days ago I began to study the local discontinuous galerkin (LDG) method, this is my first time working with a discontinuous method, so I decided to solve the poisson equation in 1D to learn the ...
5
votes
0answers
202 views

Discontinuous Galerkin for flow through porous media

I am struggling with DG methods for 2 phase flow through porous media. I managed to get the global pressure, total flux equations to work with an unconditionally stable mixed FE DG formulation as ...
4
votes
2answers
2k views

The meaning of conservative discretization in Galerkin FEM and Discontinuous Galerkin

I do understand the meanning of "conservative discretization" within the FVM/FDM framework, indeed it is well explained in this post. Now, according to the table in this slide (pp.8), it concludes: ...
4
votes
1answer
272 views

DG-FEM integration by parts

I am going through the book of Hesthaven and Warburton on discontinuous Galerkin methods. I have difficulties understanding some basic steps in the calculations. Consider the PDE: $$\frac{\partial u}...
4
votes
1answer
613 views

How to handle inflow and outflow boundaries for a non-linear convection-diffusion equation (DGFEM)

Following "A conservative DGM for Convection-Diffusion and Navier-Stokes Problems" (Oden and Baumann), if we have a linear convection-diffusion equation of the following form: $$ \nabla\cdot(\mathbf{b}...
4
votes
1answer
311 views

Quadrature and quadrature-free discontinuous galerkin method for non-linear PDE

Quadrature-free DG method using nodal Lagrangian basis are computationally very efficient. I have seen many papers using this method for linear PDE but almost no literature for non-linear PDE like ...
3
votes
1answer
463 views

Intro to DG Finite Element methods

I wrote a number of 1D/2D FE and FD programs as a bachelor student, but the main problem I continually came into contact with was gradient shocks related to convection/diffusion problems in convection-...
3
votes
2answers
618 views

What's the difference between C0 penalty methods and Discontinuous Galerkin methods?

I am trying to understand the DG FEM methods, but I got lost in their definitions. In some papers I read that the "C0 penalty method" is one example of the DG method, but sometimes they are separated ...
3
votes
2answers
1k views

Slow convergence of Newton's method for finite elements

The application is a simple non-linear advection diffusion problem (steady state) using DGFEM. My error at each iteration is given by $$ e_{n+1} = ||\mathbf{J}^{-1}(\mathbf{u}_{n})\mathbf{F}(\mathbf{u}...
3
votes
1answer
4k views

what is the difference between non-conformal and conformal?

So far what I understand is that two neighbouring elements are conformal if their edges and faces match exactly, whereas with non-conformal elements this is not the case. For instance, h-refinement ...
3
votes
2answers
492 views

Nodal DG method and limiters for hyperbolic conservation laws

All the papers I have seen on DG methods for hyperbolic conservation laws together with limiters to compute discontinuous solutions make use of Taylor polynomial basis (Pk basis) or Legendre ...
3
votes
1answer
349 views

Convergence of interior penalty DG methods

I’m currently having some issues with my routine for the linear advection-diffusion problem. The model problem is as follows: $$ \nabla\cdot(\mathbf{s} u) - \nabla\cdot(\kappa\nabla u) = f, \;\;\;\...
3
votes
1answer
187 views

calculation of the right hand side of DG FEM (with code)

I got stuck with Hestaven/Warburton's dG-FEM Matlab code. Starting with the file AdvecRHS1D.m, we see in line 11 ...
3
votes
1answer
144 views

Discontinuous Galerkin FEM : Control points are mid-points of edges instead of nodes

I am thinking to use discontinuous galerkin FEM (DGFEM) method to estimate discontinuous displacement field $u: \Omega \rightarrow \mathbb{R}^2$ at the crack surface of a material. The domain is ...
3
votes
1answer
494 views

Recommendations on FEM software for implementing Nitsche's method on interfaces between matching meshes?

Suppose: I have two domains, $\Omega_{1} = [0, 1/2] \times [0, 1]$ and $\Omega_{2} = [1/2, 1] \times [0, 1]$. The domains share an interface $\Gamma = \{1/2\} \times [0, 1] = \partial\Omega_{1} \cap \...
3
votes
2answers
1k views

How to properly use polynomial projection to get values at visualization nodes?

I am trying to implementing a nodal discontinuous Galerkin spectral element method for linear and non-linear systems of equations. The solution at each time step is given at ...
3
votes
1answer
98 views

In mixed elliptic formulation, what are the weakest requirements to ensure the flux is in $H^1$?

In the book Mixed Finite Element Methods and Applications by Boffi, Brezzi, and Fortin there is a pretty long discussion about why the Raviart-Thomas (RT) projection is only defined for functions in $...
3
votes
1answer
210 views

Is it possible to solve Euler equation numerically without using any flux limiter (in DG scheme)?

I have recently learned about Discontinuous Galerkin method to solve differential equations and I was trying to implement it to solve Euler equation. For now, consider the standard Sod Shock Tube Case....
3
votes
0answers
62 views

Explicit DG time step restriction for compressible Navier-Stokes equations

Hesthaven's book 1 mentions the following time step restriction for Navier-Stokes equations (see (7.32) in 2008 edition) $$ \Delta t \approx \frac{h}{N^2} \frac{C}{|u| + |c| + \frac {N^2 \mu}{h}} $$ ($...
3
votes
0answers
111 views

Piecewise Constant Enrichment of the continuous galerkin method

I am interested to study crack propagation in a hyperelastic material in a variational setting. The crack surface exhibits a jump discontinuity. The function space for displacement field should ...
3
votes
0answers
691 views

1D Discontinuous Galerkin - Lagrange vs Legendre Basis

Consider the 1D advection equation in its strong and weak forms $$\ u_t + a u_x = 0 $$ $$\ \int_{x_{j-0.5}}^{x_{j+0.5}} w u_t \ dx - a \int_{x_{j-0.5}}^{x_{j+0.5}} w_x u \ dx + a [w(x)\hat{u}]_{...
3
votes
0answers
136 views

Mixed DG for Poisson with mixed BC's

I am trying to find a good reference on a proper weak formulation for mixed DG (Raviart Thomas and DG) formulation for a Poisson equation with mixed boundary conditions. Can anyone suggest a good ...
2
votes
1answer
779 views

Discontinuous Galerkin method VS Continuous Galerkin method Degrees of freedom

I was looking into the Book of Riviere " Discontinuous Galerkin Methods for solving Elliptic and Parabolic Equations". In the comparaison of section 2.12 (copied below), the example of rectangular ...
2
votes
1answer
630 views

Spectral methods, Spectral Volume methods, Spectral Difference methods

Could someone explain the link (if any) between the spectral methods (SM), as presented for example here and the so called spectral volume methods (SV) and spectral difference methods (SD) for CFD ? ...
2
votes
1answer
125 views

Differentiation Matrix In DG-FEM - Hesthaven/Warburton

In the book of Hesthaven and Warburton on discontinual Galerkin methods the authors give motivation to the differentiation matrix (page 52), referred to as $D_r(i,j)=\frac{dl_j}{dr}|_{r_i}$ where $l_i(...