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.

We establish uniform error bounds of time-splitting Fourier pseudospectral (TSFP) methods for the nonlinear Klein–Gordon equation (NKGE) with weak power-type nonlinearity and O(1) initial data, while the nonlinearity strength is characterized by εp with a constant p ∈ N+ and a dimensionless parameter ε ∈ (0, 1], for the long-time dynamics up to the time at O(ε^{−β}) with 0 ≤ β ≤ p. In fact, when 0 < ε ≪ 1, the problem is equivalent to the long-time dynamics of NKGE with small initial data and O(1) nonlinearity strength, while the amplitude of the initial data (and the solution) is at O(ε). By reformulating the NKGE into a relativistic nonlinear Schr ̈odinger equation, we adapt the TSFP method to discretize it numerically. By using the method of mathematical induction to bound the numerical solution, we prove uniform error bounds at O(h^m + ε^{p−β} τ^2) of the TSFP method with h mesh size, τ time step and m ≥ 2 depending on the regularity of the solution. The error bounds are uniformly accurate for the long-time simulation up to the time at O(ε^{−β}) and uniformly valid for ε ∈ (0, 1]. Especially, the error bounds are uniformly at the second order rate for the large time step τ = O(ε^{−(p−β)/2}) in the parameter regime 0 ≤ β < p. Numerical results are reported to confirm our error bounds in the long-time regime. Finally, the TSFP method and its error bounds are extended to a highly oscillatory complex NKGE which propagates waves with wavelength at O(1) in space and O(ε^β) in time and wave velocity at O(ε^{−β}).

In this paper, we present a constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) for parabolic equations with multiscale coefficients, arising from applications in porous media. We will present the construction of CEM-GMsFEM and rigorously analyze its convergence for the parabolic equations. The convergence rate is characterized by the coarse grid size and the eigenvalue decay of local spectral problems, but is independent of the scale length and contrast of the media. The analysis shows that the method has a first order convergence rate with respect to coarse grid size in the energy norm and second order convergence rate with respect to coarse grid size in L^2 norm under some appropriate assumptions. For the temporal discretization, finite difference techniques are used and the convergence analysis of full discrete scheme is given. Moreover, a posteriori error estimator is

We prove the convergence of a discontinuous Galerkin method approximating the 2-D incompressible Euler equations with discontinuous initial vorticity:$omega_0 in L^2(Omega)$. Furthermore, when $omega_0 in L^infty(Omega)$, the whole sequence is shown to be strongly convergent. This is the first convergence result in numerical approximations of this general class of discontinuous flows. Some important flows such as vortex patches belong to this class.

In this paper, we apply local discontinuous Galerkin methods to the compressible wormhole propagation. Optimal error estimates for the pressure, velocity, porosity and concentration in different norms are established on non-uniform grids. Numerical experiments are presented to verify the theoretical analysis and show the good performance of the LDG scheme for compressible wormhole propagation.