5,264 research outputs found
Dispersion and dissipation error in high-order Runge-Kutta discontinuous Galerkin discretisations of the Maxwell equations
Different time-stepping methods for a nodal high-order discontinuous Galerkin discretisation of the Maxwell equations are discussed. A comparison between the most popular choices of Runge-Kutta (RK) methods is made from the point of view of accuracy and computational work. By choosing the strong-stability-preserving Runge-Kutta (SSP-RK) time-integration method of order consistent with the polynomial order of the spatial discretisation, better accuracy can be attained compared with fixed-order schemes. Moreover, this comes without a significant increase in the computational work. A numerical Fourier analysis is performed for this Runge-Kutta discontinuous Galerkin (RKDG) discretisation to gain insight into the dispersion and dissipation properties of the fully discrete scheme. The analysis is carried out on both the one-dimensional and the two-dimensional fully discrete schemes and, in the latter case, on uniform as well as on non-uniform meshes. It also provides practical information on the convergence of the dissipation and dispersion error up to polynomial order 10 for the one-dimensional fully discrete scheme
Discontinuous Galerkin Methods with Trefftz Approximation
We present a novel Discontinuous Galerkin Finite Element Method for wave
propagation problems. The method employs space-time Trefftz-type basis
functions that satisfy the underlying partial differential equations and the
respective interface boundary conditions exactly in an element-wise fashion.
The basis functions can be of arbitrary high order, and we demonstrate spectral
convergence in the \Lebesgue_2-norm. In this context, spectral convergence is
obtained with respect to the approximation error in the entire space-time
domain of interest, i.e. in space and time simultaneously. Formulating the
approximation in terms of a space-time Trefftz basis makes high order time
integration an inherent property of the method and clearly sets it apart from
methods, that employ a high order approximation in space only.Comment: 14 pages, 12 figures, preprint submitted at J Comput Phy
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
- …