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).

26 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
14 votes
0 answers
516 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},...
Justin Dong's user avatar
5 votes
0 answers
215 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 ...
JPL's user avatar
  • 81
3 votes
0 answers
72 views

Typo in a-priori error estimate in a Discontinuous Galerkin paper

I'm looking at this famous paper which is available in the link below: Franco Brezzi, LD Marini, Endre Süli, Discontinuous Galerkin methods for first-order hyperbolic problems, Mathematical Models ...
FEGirl's user avatar
  • 323
3 votes
0 answers
99 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}} $$ ($...
Zxcvasdf's user avatar
  • 141
3 votes
0 answers
137 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 ...
hari's user avatar
  • 95
3 votes
0 answers
164 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 ...
JPL's user avatar
  • 81
2 votes
0 answers
91 views

Modal representations of nodal tensor product Galerkin elements

Nodal discontinuous Galerkin methods on simplices, like those described in Hesthaven and Warburton, have the nice property that the number of nodes is equal to the minimum number needed to represent a ...
Aurelius's user avatar
  • 2,365
2 votes
0 answers
120 views

FEM libraries with trace spaces

To implement hybridizable discontinuous Galerkin methods, one needs finite element spaces defined on the skeleton of the mesh. deal.II has support for HDG through FE_FaceQ class which provides ...
cfdlab's user avatar
  • 3,028
2 votes
0 answers
642 views

Jacobian-Free Newton-Krylov vs explicitly forming jacobian in DG

For a given discontinuous galerkin (DG) implementation for Navier-Stokes, targeting 10,000 to 1,000,000 4th order cells in 3D, I'm using PETSc's suite of linear/non-linear solvers on the back-end. It ...
Aurelius's user avatar
  • 2,365
1 vote
0 answers
95 views

Galerkin projection in AMG

In the context of Classical AMG for elliptic problems discretised with finite elements (DG or CG), one has the (fine) matrix of the problem, say $A_0$, and the coarser operators of the hierarchy $\{...
FEGirl's user avatar
  • 323
1 vote
0 answers
136 views

Are both centered and upwind numerical fluxes correct in DGTD?

In discontinuous Galerkin time domain (DGTD) method, a critical concept is the numerical flux that is used to link neighbouring elements. The numerical flux is however not unique. The popular choices ...
Pu Zhang's user avatar
  • 302
1 vote
0 answers
101 views

Is it possible to use a fixed point iteration for solving this nonlinear system?

Consider the following differential equation \begin{align} \frac{\partial f(u)}{\partial x} &= g(x), \ \ x\in [x_{L},x_{R}] \label{Eq2.2} \\ u(x_{L}) &= g_{1} \end{align} where $f(u)$ is a ...
TheComander's user avatar
1 vote
0 answers
72 views

Solution predictors for accelerating convergence in nonlinear FEM

I am looking for the details of commonly-used predictors for accelerating the convergence of iterations using Newton-Raphson scheme for nonlinear problems in FEM. I am looking specifically for static ...
Chenna K's user avatar
  • 934
1 vote
0 answers
54 views

Cell-centered DG extension to the two-point flux approximation scheme

A current problem that I am working on requires me to compute the solution from the heat diffusion evolution on a discontinuous function. More precisely - I have a Delaunay triangulation and within ...
lightxbulb's user avatar
  • 2,112
1 vote
0 answers
77 views

Question regarding 1D implementation of the DG method

I'm pretty new to the DG method and have been writing a 1D code to help me understand the coding aspect. With respect a reference, I've been following these notes https://www3.nd.edu/~zxu2/...
Motig's user avatar
  • 11
1 vote
0 answers
110 views

Do DG methods for the Helmholtz equation always return positive quantities?

Helmholtz Diffusion equation with reaction term: $$ k\Delta u + u = f ~ \text{in} ~\Omega $$ $$ \nabla u \cdot \mathbf{n} = 0 ~ \text{in} ~\partial \Omega $$ For sufficiently small $k$ (relative to ...
balborian's user avatar
  • 601
1 vote
0 answers
53 views

How do we implement Parameter free generalised Moment limiter in 1D Case in Discontinuous Galerkin methods?

I am referring to this paper:- "A Parameter-Free Generalized Moment Limiter for High- Order Methods on Unstructured Grids " by Michael Yang and Z.J. Wang. http://dept.ku.edu/~cfdku/papers/AIAA-2009-...
Manish's user avatar
  • 163
1 vote
0 answers
107 views

Discontinuos Galerkin Method - inhomogeneous advection problem

I'm currently trying to get into this topic. I've learned that the basic scheme for the advection problem ($D_{x}u+a*D_{t}u=0$) can be solved in a scheme like $$ M^{k}\frac{d}{dt}u^{k}_{h}-(S^{k})^{T}...
NG2207's user avatar
  • 11
1 vote
0 answers
380 views

How to perform Taylor series expansion consisting of cell averaged derivatives in a computational element?

I have encountered the term cell "Cell averaged derivatives" and "Taylor series expansion using cell averaged derivatives about centroid" in the context of polynomial representation in any cell (...
Manish's user avatar
  • 163
1 vote
0 answers
154 views

How do we derive the elemental equation for Discontinuous Galerkin method using Centered Numerical Flux?

If we take a 1st order polynomial approximation in each cell, we can find the (2x2) -mass matrix, differentiation matrix, and flux matrix through the integration of Lagrange polynomials. However, I ...
Manish's user avatar
  • 163
1 vote
0 answers
87 views

Time integration for elastodynamics

I'd like to solve the elastodynamics equation for a problem with fast scales from an impact (although no shocks). $\rho\ddot{\mathbf{u}} - \nabla \cdot \sigma(\mathbf{u}) = \mathbf{f}$ In structural ...
balborian's user avatar
  • 601
1 vote
0 answers
131 views

understanding interior penalty jump of basis function

If cardinal basis functions are used (i.e. $\psi_{ij}=1$ iff $i=j$, and 0 otherwise) in interior penalty methods for elliptic equations such as SIPG & NIPG, shoudn't the jump in basis functions ...
danny's user avatar
  • 233
1 vote
0 answers
423 views

Solving the Jump function in discontinuous Galerkin method

I am trying to solve the elasticity equation using Crouzeix-Raviart element. Since this can not be solved using continuous FEM I am trying to solve this by discontinuous Galerkin method. In the ...
Jose's user avatar
  • 11
0 votes
0 answers
59 views

solve a coupled PDE system with some discontinuity by a mixed FEM

$$ \begin{cases} u_t=f(u,v,\nabla{v}),\\ \Delta{[g(u)v]}=0. \end{cases}$$ I want to solve the above PDE system by a mixed FEM, that is, $u_t=f(u,v,\nabla{v})$ by discontinuous Galerkin (DG), where $f$...
feynman's user avatar
  • 269
0 votes
0 answers
155 views

Discontinuous Galerkin failing to converge Euler equations under p-refinement

I am solving the steady state compressible Euler equations in conservation form in 2D in a rectangular domain with a discontinuous Galerkin (DG) code I have developed. As a test, the boundary ...
Wil's user avatar
  • 63
0 votes
0 answers
52 views

Can you describe the Galerkin numerical method to solve the wave equation?

How would you describe the Galerkin method to solving the 3D wave equation $$u_{tt}= c^2\Delta u$$ to someone who wants to implement it immediately? More precisely, we want to solve the Cauchy problem ...
NotaChoice's user avatar