2,261 research outputs found
On the convergence of a shock capturing discontinuous Galerkin method for nonlinear hyperbolic systems of conservation laws
In this paper, we present a shock capturing discontinuous Galerkin (SC-DG)
method for nonlinear systems of conservation laws in several space dimensions
and analyze its stability and convergence. The scheme is realized as a
space-time formulation in terms of entropy variables using an entropy stable
numerical flux. While being similar to the method proposed in [14], our
approach is new in that we do not use streamline diffusion (SD) stabilization.
It is proved that an artificial-viscosity-based nonlinear shock capturing
mechanism is sufficient to ensure both entropy stability and entropy
consistency, and consequently we establish convergence to an entropy
measure-valued (emv) solution. The result is valid for general systems and
arbitrary order discontinuous Galerkin method.Comment: Comments: Affiliations added Comments: Numerical results added,
shortened proo
A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes
We are interested in the approximation of a steady hyperbolic problem. In
some cases, the solution can satisfy an additional conservation relation, at
least when it is smooth. This is the case of an entropy. In this paper, we
show, starting from the discretisation of the original PDE, how to construct a
scheme that is consistent with the original PDE and the additional conservation
relation. Since one interesting example is given by the systems endowed by an
entropy, we provide one explicit solution, and show that the accuracy of the
new scheme is at most degraded by one order. In the case of a discontinuous
Galerkin scheme and a Residual distribution scheme, we show how not to degrade
the accuracy. This improves the recent results obtained in [1, 2, 3, 4] in the
sense that no particular constraints are set on quadrature formula and that a
priori maximum accuracy can still be achieved. We study the behavior of the
method on a non linear scalar problem. However, the method is not restricted to
scalar problems
Affordable, Entropy Conserving and Entropy Stable Flux Functions for the Ideal MHD Equations
In this work, we design an entropy stable, finite volume approximation for
the ideal magnetohydrodynamics (MHD) equations. The method is novel as we
design an affordable analytical expression of the numerical interface flux
function that discretely preserves the entropy of the system. To guarantee the
discrete conservation of entropy requires the addition of a particular source
term to the ideal MHD system. Exact entropy conserving schemes cannot dissipate
energy at shocks, thus to compute accurate solutions to problems that may
develop shocks, we determine a dissipation term to guarantee entropy stability
for the numerical scheme. Numerical tests are performed to demonstrate the
theoretical findings of entropy conservation and robustness.Comment: arXiv admin note: substantial text overlap with arXiv:1509.06902;
text overlap with arXiv:1007.2606 by other author
On discretely entropy conservative and entropy stable discontinuous Galerkin methods
High order methods based on diagonal-norm summation by parts operators can be
shown to satisfy a discrete conservation or dissipation of entropy for
nonlinear systems of hyperbolic PDEs. These methods can also be interpreted as
nodal discontinuous Galerkin methods with diagonal mass matrices. In this work,
we describe how use flux differencing, quadrature-based projections, and
SBP-like operators to construct discretely entropy conservative schemes for DG
methods under more arbitrary choices of volume and surface quadrature rules.
The resulting methods are semi-discretely entropy conservative or entropy
stable with respect to the volume quadrature rule used. Numerical experiments
confirm the stability and high order accuracy of the proposed methods for the
compressible Euler equations in one and two dimensions
A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting
An explicit moving boundary method for the numerical solution of
time-dependent hyperbolic conservation laws on grids produced by the
intersection of complex geometries with a regular Cartesian grid is presented.
As it employs directional operator splitting, implementation of the scheme is
rather straightforward. Extending the method for static walls from Klein et
al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme
calculates fluxes needed for a conservative update of the near-wall cut-cells
as linear combinations of standard fluxes from a one-dimensional extended
stencil. Here the standard fluxes are those obtained without regard to the
small sub-cell problem, and the linear combination weights involve detailed
information regarding the cut-cell geometry. This linear combination of
standard fluxes stabilizes the updates such that the time-step yielding
marginal stability for arbitrarily small cut-cells is of the same order as that
for regular cells. Moreover, it renders the approach compatible with a wide
range of existing numerical flux-approximation methods. The scheme is extended
here to time dependent rigid boundaries by reformulating the linear combination
weights of the stabilizing flux stencil to account for the time dependence of
cut-cell volume and interface area fractions. The two-dimensional tests
discussed include advection in a channel oriented at an oblique angle to the
Cartesian computational mesh, cylinders with circular and triangular
cross-section passing through a stationary shock wave, a piston moving through
an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil
profile.Comment: 30 pages, 27 figures, 3 table
- …