Alexander Jaust

Alexander Jaust, Balthasar Reuter, Vadym Aizinger et al.

The third paper in our series on open source MATLAB / GNU Octave implementation of the discontinuous Galerkin (DG) method(s) focuses on a hybridized formulation. The main aim of this ongoing work is to develop rapid prototyping techniques covering a range of standard DG methodologies and suitable for small to medium sized applications. Our FESTUNG package relies on fully vectorized matrix / vector operations throughout, and all details of the implementation are fully documented. Once again, great care was taken to maintain a direct mapping between discretization terms and code routines as well as to ensure full compatibility to GNU Octave. The current work formulates a hybridized DG scheme for linear advection problem, describes hybrid approximation spaces on the mesh skeleton, and compares the performance of this discretization to the standard (element-based) DG method for different polynomial orders.

Jochen Schütz, David C. Seal, Alexander Jaust

In this work, we construct novel discretizations for the unsteady convection-diffusion equation. Our discretization relies on multiderivative time integrators together with a novel discretization that reduces the total number of unknowns for the solver. These type of temporal discretizations come from an umbrella class of methods that include Lax-Wendroff (Taylor) as well as Runge-Kutta methods as special cases. We include two-point collocation methods with multiple time derivatives as well as a sixth-order fully implicit collocation method that only requires a total of three stages. Numerical results for a number of sample linear problems indicate the expected order of accuracy and indicate we can take arbitrarily large time steps.

Alexander Jaust, Jochen Schütz, David C. Seal

In this paper we apply implicit two-derivative multistage time integrators to viscous conservation laws in one and two dimensions. The one dimensional solver discretizes space with the classical discontinuous Galerkin (DG) method, and the two dimensional solver uses a hybridized discontinuous Galerkin (HDG) spatial discretization for efficiency. We propose methods that permit us to construct implicit solvers using each of these spatial discretizations, wherein a chief difficulty is how to handle the higher derivatives in time. The end result is that the multiderivative time integrator allows us to obtain high-order accuracy in time while keeping the number of implicit stages at a minimum. We show numerical results validating and comparing methods.