In this paper we generalize a new type of limiters based on the weighted essentially non-oscillatory (WENO) finite volume methodology for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving nonlinear hyperbolic conservation laws, which were recently developed for structured meshes, to two-dimensional unstructured triangular meshes. The key idea of such limiters is to use the entire polynomials of the DG solutions from the troubled cell and its immediate neighboring cells, and then apply the classical WENO procedure to form a convex combination of these polynomials based on smoothness indicators and nonlinear weights, with suitable adjustments to guarantee conservation. The main advantage of this new limiter is its simplicity in implementation, especially for the unstructured meshes considered in this paper, as only information from immediate neighbors is needed and the usage of complicated geometric information of the meshes is largely avoided. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good performance of this procedure.

Two-dimensional interaction between two Taylor vortices is simulated systematically through solving the two-dimensional, unsteady compressible Navier-Stokes equations using a fifth order weighted essentially nonoscillatory finite difference scheme. The main purpose of this study is to reveal the mechanism of sound generation in two-dimensional interaction of two Taylor vortices. Based on an extensive parameter study on the evolution of the vorticity field, we classify the interaction of two Taylor vortices into four types. The first type is the interaction of two counter-rotating vortices with similar strengths. The second type is the interaction of two co-rotating vortices without merging. The third type is the merging of two co-rotating vortices. The fourth type is the interaction of two vortices with a large difference in their strengths or scales. The mechanism of sound generation is analyzed.

Homotopy continuation is an efficient tool for solving polynomial systems. Its efficiency relies on utilizing adaptive stepsize and adaptive precision path tracking, and endgames. In this article, we apply homotopy continuation to solve steady state problems of hyperbolic conservation laws. A third-order accurate finite difference weighted essentially non-oscillatory (WENO) scheme with Lax-Friedrichs flux splitting is utilized to derive the difference equation. This new approach is free of the CFL condition constraint. Extensive numerical examples in both scalar and system test problems in one and two dimensions demonstrate the efficiency and robustness of the new method.

We investigate the angular distribution of Ly$\alpha$ photons transferring in or emergent from an optically thick medium. Since the evolutions of specific intensity $I$ in the frequency space and the angular space are coupled with each other, we first develop the WENO numerical solver in order to find the time-dependent solutions of the integro-differential equation of $I$ in frequency and angular space simultaneously. We first show that the solutions with the Eddington approximation, which assume $I$ to be linearly dependent on the angular variable $\mu$, yield similar frequency profiles of the photon flux as that without the Eddington approximation. However, the solutions of the $\mu$ distribution evolution are significantly different from that given by Eddington approximation. First, the angular distribution of $I$ are found to be substantially dependent on the frequency of photons. For photons with the resonant frequency $\nu_0$, $I$ contains only a linear term of $\mu$. For photons with frequency at the double peaks of the flux, the $\mu$-distribution is highly anisotropic, in which most photons are in the direction of radial forward. Moreover, either at $\nu_0$ or at the double peaks, the $\mu$ distributions actually are independent of the initial $\mu$ distribution of photons of the source. This is because the photons with frequency either of $\nu_0$ or of the double peaks undergo the process of forgetting their initial conditions due to the resonant scattering. We also show that the optically thick medium is a collimator of photons at the double peaks. Photons from the double peaks form a forward beam with very small spread angle.

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.