2,900 research outputs found
On the stability of projection methods for the incompressible Navier-Stokes equations based on high-order discontinuous Galerkin discretizations
The present paper deals with the numerical solution of the incompressible
Navier-Stokes equations using high-order discontinuous Galerkin (DG) methods
for discretization in space. For DG methods applied to the dual splitting
projection method, instabilities have recently been reported that occur for
coarse spatial resolutions and small time step sizes. By means of numerical
investigation we give evidence that these instabilities are related to the
discontinuous Galerkin formulation of the velocity divergence term and the
pressure gradient term that couple velocity and pressure. Integration by parts
of these terms with a suitable definition of boundary conditions is required in
order to obtain a stable and robust method. Since the intermediate velocity
field does not fulfill the boundary conditions prescribed for the velocity, a
consistent boundary condition is derived from the convective step of the dual
splitting scheme to ensure high-order accuracy with respect to the temporal
discretization. This new formulation is stable in the limit of small time steps
for both equal-order and mixed-order polynomial approximations. Although the
dual splitting scheme itself includes inf-sup stabilizing contributions, we
demonstrate that spurious pressure oscillations appear for equal-order
polynomials and small time steps highlighting the necessity to consider inf-sup
stability explicitly.Comment: 31 page
A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations
A discontinuous Galerkin method for the ideal 5 moment two-fluid plasma
system is presented. The method uses a second or third order discontinuous
Galerkin spatial discretization and a third order TVD Runge-Kutta time stepping
scheme. The method is benchmarked against an analytic solution of a dispersive
electron acoustic square pulse as well as the two-fluid electromagnetic shock
and existing numerical solutions to the GEM challenge magnetic reconnection
problem. The algorithm can be generalized to arbitrary geometries and three
dimensions. An approach to maintaining small gauge errors based on error
propagation is suggested.Comment: 40 pages, 18 figures
Well-balanced -adaptive and moving mesh space-time discontinuous Galerkin method for the shallow water equations
In this article we introduce a well-balanced discontinuous Galerkin method for the shallow water equations on moving meshes. Particular emphasis will be given on -adaptation in which mesh points of an initially uniform mesh move to concentrate in regions where interesting behaviour of the solution is observed. Obtaining well-balanced numerical schemes for the shallow water equations on fixed meshes is nontrivial and has been a topic of much research. In [S. Rhebergen, O. Bokhove, J.J.W. van der Vegt, Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations, J. Comput. Phys. 227 (2008) 1887–1922] we introduced a well-balanced discontinuous Galerkin method using the theory of weak solutions for nonconservative products introduced in [G. Dal Maso, P.G. LeFloch, F. Murat, Definition and weak stability of nonconservative products, J. Math. Pures Appl. 74 (1995) 483–548]. In this article we continue this approach and prove well-balancedness of a discontinuous Galerkin method for the shallow water equations on moving meshes. Numerical simulations are then performed to verify the -adaptive method in combination with the space-time discontinuous Galerkin method against analytical solutions and showing its robustness on more complex problems
- …