803 research outputs found
Stabilization in relation to wavenumber in HDG methods
Simulation of wave propagation through complex media relies on proper
understanding of the properties of numerical methods when the wavenumber is
real and complex. Numerical methods of the Hybrid Discontinuous Galerkin (HDG)
type are considered for simulating waves that satisfy the Helmholtz and Maxwell
equations. It is shown that these methods, when wrongly used, give rise to
singular systems for complex wavenumbers. A sufficient condition on the HDG
stabilization parameter for guaranteeing unique solvability of the numerical
HDG system, both for Helmholtz and Maxwell systems, is obtained for complex
wavenumbers. For real wavenumbers, results from a dispersion analysis are
presented. An asymptotic expansion of the dispersion relation, as the number of
mesh elements per wave increase, reveal that some choices of the stabilization
parameter are better than others. To summarize the findings, there are values
of the HDG stabilization parameter that will cause the HDG method to fail for
complex wavenumbers. However, this failure is remedied if the real part of the
stabilization parameter has the opposite sign of the imaginary part of the
wavenumber. When the wavenumber is real, values of the stabilization parameter
that asymptotically minimize the HDG wavenumber errors are found on the
imaginary axis. Finally, a dispersion analysis of the mixed hybrid
Raviart-Thomas method showed that its wavenumber errors are an order smaller
than those of the HDG method
Finite volume methods for unidirectional dispersive wave models
We extend the framework of the finite volume method to dispersive
unidirectional water wave propagation in one space dimension. In particular we
consider a KdV-BBM type equation. Explicit and IMEX Runge-Kutta type methods
are used for time discretizations. The fully discrete schemes are validated by
direct comparisons to analytic solutions. Invariants conservation properties
are also studied. Main applications include important nonlinear phenomena such
as dispersive shock wave formation, solitary waves and their various
interactions.Comment: 25 pages, 12 figures, 51 references. Other authors papers can be
downloaded at http://www.lama.univ-savoie.fr/~dutykh
Finite volume methods for unidirectional dispersive wave model
We extend the framework of the finite volume method to dispersive unidirectional water wave propagation in one space dimension. In particular, we consider a KdV–BBM-type equation. Explicit and implicit–explicit Runge–Kutta-type methods are used for time discretizations. The fully discrete schemes are validated by direct comparisons to analytic solutions. Invariants’ conservation properties are also studied. Main applications include important nonlinear phenomena such as dispersive shock wave formation, solitary waves, and their various interaction
A Space-Time Discontinuous Galerkin Trefftz Method for time dependent Maxwell's equations
We consider the discretization of electromagnetic wave propagation problems
by a discontinuous Galerkin Method based on Trefftz polynomials. This method
fits into an abstract framework for space-time discontinuous Galerkin methods
for which we can prove consistency, stability, and energy dissipation without
the need to completely specify the approximation spaces in detail. Any method
of such a general form results in an implicit time-stepping scheme with some
basic stability properties. For the local approximation on each space-time
element, we then consider Trefftz polynomials, i.e., the subspace of
polynomials that satisfy Maxwell's equations exactly on the respective element.
We present an explicit construction of a basis for the local Trefftz spaces in
two and three dimensions and summarize some of their basic properties. Using
local properties of the Trefftz polynomials, we can establish the
well-posedness of the resulting discontinuous Galerkin Trefftz method.
Consistency, stability, and energy dissipation then follow immediately from the
results about the abstract framework. The method proposed in this paper
therefore shares many of the advantages of more standard discontinuous Galerkin
methods, while at the same time, it yields a substantial reduction in the
number of degrees of freedom and the cost for assembling. These benefits and
the spectral convergence of the scheme are demonstrated in numerical tests
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
Localized solutions and filtering mechanisms for the discontinuous Galerkin semi-discretizations of the 1-d wave equation
We perform a complete Fourier analysis of the semi-discrete 1-d wave equation
obtained through a P1 discontinuous Galerkin (DG) approximation of the
continuous wave equation on an uniform grid. The resulting system exhibits the
interaction of two types of components: a physical one and a spurious one,
related to the possible discontinuities that the numerical solution allows.
Each dispersion relation contains critical points where the corresponding group
velocity vanishes. Following previous constructions, we rigorously build wave
packets with arbitrarily small velocity of propagation concentrated either on
the physical or on the spurious component. We also develop filtering mechanisms
aimed at recovering the uniform velocity of propagation of the continuous
solutions. Finally, some applications to numerical approximation issues of
control problems are also presented.Comment: 6 pages, 2 figure
- …