14,421 research outputs found
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 hybridizable discontinuous Galerkin method for electromagnetics with a view on subsurface applications
Two Hybridizable Discontinuous Galerkin (HDG) schemes for the solution of
Maxwell's equations in the time domain are presented. The first method is based
on an electromagnetic diffusion equation, while the second is based on
Faraday's and Maxwell--Amp\`ere's laws. Both formulations include the diffusive
term depending on the conductivity of the medium. The three-dimensional
formulation of the electromagnetic diffusion equation in the framework of HDG
methods, the introduction of the conduction current term and the choice of the
electric field as hybrid variable in a mixed formulation are the key points of
the current study. Numerical results are provided for validation purposes and
convergence studies of spatial and temporal discretizations are carried out.
The test cases include both simulation in dielectric and conductive media
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
- …