1,409 research outputs found
Runge-Kutta discontinuous Galerkin method using WENO limiters
In [J. Qiu, C.-W. Shu, Runge-Kutta discontinuous Galerkin method using WENO limiters, SIAM Journal on Scientific Computing 26 (2005) 907-929], Qiu and Shu investigated using weighted essentially non-oscillatory (WENO) finite volume methodology as limiters for the Runge-Kutta discontinuous Galerkin (RKDG) methods for solving nonlinear hyperbolic conservation law systems on structured meshes. In this continuation paper, we extend the method to solve two-dimensional problems on unstructured meshes, with the goal of obtaining a robust and high order limiting procedure to simultaneously obtain uniform high order accuracy and sharp, nonoscillatory shock transition for RKDG methods. Numerical results are provided to illustrate the behavior of this procedure. (C) 2008 Elsevier Inc. All rights reserved
Embedded discontinuous Galerkin transport schemes with localised limiters
Motivated by finite element spaces used for representation of temperature in
the compatible finite element approach for numerical weather prediction, we
introduce locally bounded transport schemes for (partially-)continuous finite
element spaces. The underlying high-order transport scheme is constructed by
injecting the partially-continuous field into an embedding discontinuous finite
element space, applying a stable upwind discontinuous Galerkin (DG) scheme, and
projecting back into the partially-continuous space; we call this an embedded
DG scheme. We prove that this scheme is stable in L2 provided that the
underlying upwind DG scheme is. We then provide a framework for applying
limiters for embedded DG transport schemes. Standard DG limiters are applied
during the underlying DG scheme. We introduce a new localised form of
element-based flux-correction which we apply to limiting the projection back
into the partially-continuous space, so that the whole transport scheme is
bounded. We provide details in the specific case of tensor-product finite
element spaces on wedge elements that are discontinuous P1/Q1 in the horizontal
and continuous P2 in the vertical. The framework is illustrated with numerical
tests
SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
We introduce a new relativistic astrophysics code, SpECTRE, that combines a
discontinuous Galerkin method with a task-based parallelism model. SpECTRE's
goal is to achieve more accurate solutions for challenging relativistic
astrophysics problems such as core-collapse supernovae and binary neutron star
mergers. The robustness of the discontinuous Galerkin method allows for the use
of high-resolution shock capturing methods in regions where (relativistic)
shocks are found, while exploiting high-order accuracy in smooth regions. A
task-based parallelism model allows efficient use of the largest supercomputers
for problems with a heterogeneous workload over disparate spatial and temporal
scales. We argue that the locality and algorithmic structure of discontinuous
Galerkin methods will exhibit good scalability within a task-based parallelism
framework. We demonstrate the code on a wide variety of challenging benchmark
problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the
code's scalability including its strong scaling on the NCSA Blue Waters
supercomputer up to the machine's full capacity of 22,380 nodes using 671,400
threads.Comment: 41 pages, 13 figures, and 7 tables. Ancillary data contains
simulation input file
High Order Maximum Principle Preserving Semi-Lagrangian Finite Difference WENO schemes for the Vlasov Equation
In this paper, we propose the parametrized maximum principle preserving (MPP)
flux limiter, originally developed in [Z. Xu, Math. Comp., (2013), in press],
to the semi- Lagrangian finite difference weighted essentially non-oscillatory
scheme for solving the Vlasov equation. The MPP flux limiter is proved to
maintain up to fourth order accuracy for the semi-Lagrangian finite difference
scheme without any time step restriction. Numerical studies on the
Vlasov-Poisson system demonstrate the performance of the proposed method and
its ability in preserving the positivity of the probability distribution
function while maintaining the high order accuracy
- …