No abstract uploaded!

This paper deals with the numerical resolution of kinetic models for systems of self-propelled particles subject to alignment interaction and attraction-repulsion. We focus on the kinetic model considered in \cite{DM, DLMP} where alignment is taken into account in addition of an attraction-repulsion interaction potential. We apply a discontinuous Galerkin method for the free transport and non-local drift velocity together with a spectral method for the velocity variable. Then, we analyse consistency and stability of the semi-discrete scheme. We propose several numerical experiments which provide a solid validation of the method and its underlying concepts.

In this paper, we study the superconvergence behavior of the semi-discrete discontinuous Galerkin (DG) method for {\color{red} scalar} nonlinear hyperbolic equations in one spatial dimension. Superconvergence results for problems with fixed and alternating wind directions are established. On the one hand, we prove that, if the wind direction is fixed (i.e., the derivative of the flux function is bounded away from zero), both the cell average error and numerical flux error at cell interfaces converge at a rate of $2k+1$ when upwind fluxes and piecewise polynomials of degree $k$ are used. Moreover, we also prove that the function value approximation of the DG solution is superconvergent at interior right Radau points, and the derivative value approximation is superconvergent at interior left Radau points, with an order of $k+2$ and $k+1$, respectively. As a byproduct, we show a $k+2$-th order superconvergence of the DG solution towards the Gauss-Radau projection of the exact solution. On the other hand, superconvergence results for problems with alternating wind directions (i.e., the derivative of the flux function either changes sign or otherwise achieves the value zero in the domain) are also established. To be more precise, we first prove that the DG flux function is superconvergent towards a particular flux function of the exact solution, with an order of $k+2$, when Godunov fluxes are used. We then prove that the highest superconvergence rate of the DG solution itself is $k+\frac32$ when sonic points (i.e., the derivative of the flux function achieves zero) appear in the computational domain. As byproducts, we obtain superconvergence properties for the DG solution and the DG flux function at special points and for cell average. Numerical experiments demonstrate that most of our results are optimal, i.e., the superconvergence rates are sharp.

We analyze the central discontinuous Galerkin (DG) method for time-dependent linear conservation laws. In one dimension, optimal a priori $L^2$ error estimates of order $k+1$ are obtained for the semidiscrete scheme when piecewise polynomials of degree at most $k$ ($k\geq0$) are used on overlapping uniform meshes. %Our analysis is valid for both periodic boundary conditions and %for inflow-outflow boundary conditions. We then extend the analysis to multidimensions on uniform Cartesian meshes when piecewise tensor product polynomials are used on overlapping meshes. Numerical experiments are given to demonstrate the theoretical results.

We present a discontinuous Galerkin (DG) scheme with suitable quadrature rules \cite{Tchen2017entropy} for ideal compressible magnetohydrodynamic (MHD) equations on structural meshes. The semi-discrete scheme is analyzed to be entropy stable by using the symmetrizable version of the equations as introduced by Godunov \cite{godunov1972symmetric}, the entropy stable DG framework with suitable quadrature rules \cite{Tchen2017entropy}, the entropy conservative flux in \cite{chandrashekar2016entropy} inside each cell and the entropy dissipative approximate Godunov type numerical flux at cell interfaces to make the scheme entropy stable. {\color{blue}{ The main difficulty in the generalization of the results in \cite{Tchen2017entropy} is the appearance of the non-conservative ``source terms'' added in the modified MHD model introduced by Godunov \cite{godunov1972symmetric}, which do not exist in the general hyperbolic system studied in \cite{Tchen2017entropy}. Special care must be taken to discretize these ``source terms'' adequately so that the resulting DG scheme satisfies entropy stability. }} Total variation diminishing / bounded (TVD/TVB) limiters and bound-preserving limiters are applied to control spurious oscillations. We demonstrate the accuracy and robustness of this new scheme on standard MHD examples.