In this paper, we apply discontinuous Galerkin (DG) methods to solve two model equations: Krause's consensus models and pressureless Euler equations. These two models are used to describe the collisions of particles, and the distributions can be identified as density functions. If the particles are placed at a single point, then the density function turns out to be a $\delta$-function and is difficult to be well approximated numerically. In this paper, we use DG method to approximate such a singularity and demonstrate the good performance of the scheme. Since the density functions are always positive, we apply a positivity-preserving limiter to them. Moreover, for pressureless Euler equations, the velocity satisfies the maximum principle. We also construct special limiters to fulfill this requirement.

This work is devoted to the numerical computation of complex band structure $\k=\k(\omega)\in\mathbb C^3$ for positive $\omega$ of three dimensional isotropic dispersive or non-dispersive photonic crystals from the perspective of structured quadratic eigenvalue problems (QEPs). Our basic strategy is to fix two degrees of freedom in $\k\in\mathbb C^3$ and to view the remaining one as the eigenvalue of a quadratic operator pencil derived from Maxwell's equations. Then Yee's scheme is employed to discretize $\nabla\times$ and $\k\times$ operators in this quadratic operator pencil. Distinct from the others' works which either ignore or directly exploit the Hamiltonian structure of the spectrum of the resulting QEP, we reformulate this QEP into an equivalent $\top$-palindromic QEP to facilitate the use of superior structure-preserving algorithms. Ultimately we rely on the structured Arnoldi algorithm, namely the G$\top$SHIRA algorithm, to compute eigenvalues of a $\top$-skew-Hamiltonian pair which are near or in $[-2,2]$, a much narrower region than the whole positive real axis in the origin problem. Moreover, to accelerate the inner iterations of the G$\top$SHIRA algorithm, we propose the preconditioning technique, making most of the eigenmatrix, which can essentially be seen as the Kronecker product of three discrete Fourier transformation matrices, of the commutative discretized $\partial_x,\partial_y,\partial_z$ operators. The advantage of our method is discussed in detail and corroborated by several numerical results.

The main purpose of this paper is to give stability analysis and error estimates of the local discontinuous Galerkin (LDG) methods coupled with three specific implicit-explicit (IMEX) Runge-Kutta time discretization methods up to third order accuracy, for solving one-dimensional time-dependent linear fourth order partial differential equations. In the time discretization, all the lower order derivative terms are treated explicitly and the fourth order derivative term is treated implicitly. By the aid of energy analysis, we show that the IMEX-LDG schemes are unconditionally energy stable, in the sense that the time step $\dt$ is only required to be upper-bounded by a constant which is independent of the mesh size $h$. The optimal error estimate is also derived by the aid of the elliptic projection and the adjoint argument. Numerical experiments are given to verify that the corresponding IMEX-LDG schemes can achieve optimal error accuracy.

Poisson equation on point cloud with Dirichlet boundary condition plays important role in many problems. In this paper, we use the volume constraint proposed by Du et.al to handle the Dirichlet boundary condition in the point integral method for Poisson equation on point cloud. We prove that the solution given by volume constraint converges to the true solution as the point cloud converges to the underlying smooth manifold.

For numerical simulation of detonation, computational cost using uniform meshes is large due to the vast separation in both time and space scales. Adaptive mesh refinement (AMR) is advantageous for problems with vastly different scales. This paper aims to propose an AMR method with high order accuracy for numerical investigation of multi-dimensional detonation. A well-designed AMR method based on finite difference weighted essentially non-oscillatory (WENO) scheme, named as AMR\&WENO is proposed. A new cell-based data structure is used to organize the adaptive meshes. The new data structure makes it possible for cells to communicate with each other quickly and easily. In order to develop an AMR method with high order accuracy, high order prolongations in both space and time are utilized in the data prolongation procedure. Based on the message passing interface (MPI) platform, we have developed a workload balancing parallel AMR\&WENO code using the Hilbert space-filling curve algorithm. Our numerical experiments with detonation simulations indicate that the AMR\&WENO is accurate and have a high resolution. Moreover, we evaluate and compare the performance between the uniform mesh WENO scheme and the parallel AMR\&WENO method. The comparison results provide us further insight into the high performance of the parallel AMR\&WENO method.