In this paper, we develop a discontinuous Galerkin (DG) method to solve the ideal special relativistic hydrodynamics (RHD) and design a bound-preserving (BP) limiter for this scheme by extending the idea in (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 8918-8934). For RHD, the density and pressure are positive and the velocity is bounded by the speed of light. One difficulty in numerically solving the RHD in its conservative form is that the failure of preserving these physical bounds will result in ill-posedness of the problem and blowup of the code, especially in extreme relativistic cases. The standard way in dealing with this difficulty is to add extra numerical dissipation, while in doing so there is no guarantee in maintaining the high order of accuracy. Our BP limiter has the following features. It can theoretically guarantee to preserve the physical bounds for the numerical solution and maintain its designed high order accuracy. The limiter is local to the cell and hence is very easy to implement. Moreover, it renders $L^1$-stability to the numerical scheme. Numerical experiments are performed to demonstrate the good performance of this bound-preserving DG scheme. Even though we only discuss the BP limiter for DG schemes, it can be applied to high order finite volume schemes, such as weighted essentially non-oscillatory (WENO) finite volume schemes as well.

No abstract uploaded!

No abstract uploaded!

No abstract uploaded!

We propose a semi-discrete scheme for 2D Keller-Segel equations based on a symmetrization reformation, which is equivalent to the convex splitting method and is free of any nonlinear solver. We show that, this new scheme is unconditionally stable as long as the initial condition does not exceed certain threshold, and it asymptotically preserves the quasi-static limit in the transient regime. Furthermore, we prove that the fully discrete scheme is conservative and positivity preserving, which makes it ideal for simulations. The analogical schemes for the radial symmetric cases and the subcritical degenerate cases are also presented and analyzed. With extensive numerical tests, we verify the claimed properties of the methods and demonstrate their superiority in various challenging applications.