In this paper we analyze the explicit Runge-Kutta discontinuous Galerkin (RKDG) methods for the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology. The RKDG methods use a third order explicit total-variation-diminishing Runge-Kutta (TVDRK3) time discretization and upwinding numerical fluxes. By using the energy method, under a standard CFL condition, we obtain $L^2$ stability for general solutions and a priori error estimates when the solutions are smooth enough. The theoretical results are proved for piecewise polynomials with any degree $k \geq 1$. Finally, since the solutions to this system are non-negative, we discuss a positivity-preserving limiter to preserve positivity without compromising accuracy. Numerical results are provided to demonstrate these RKDG methods.

We use a suite of cosmological hydrodynamic simulations, run by two fixed grid codes, in the context of LCDM model to investigate the properties of solenoidal and dilatational, i.e, curl and compressive, motions of the intergalactic medium (IGM), and the impact of numerical viscosity on turbulence. The codes differ only in the spatial difference discretization, one is a second-order TVD scheme and the other is a fifth-order positivity-preserving WENO scheme. We find that (1) The vortical motion grows rapidly since $z=2$, and reaches $\sim 10 km/s -75 km/s$, depending on the gas density, at low redshifts. The small-scale compressive ratio $r_{CS}$ drops from $0.8$ to $0.5$ during the same period. (2) Power spectra of the solenoidal velocity possess two regimes, $\propto k^{-0.89}$ and $\propto k^{-2.02}$, while the total velocity and dilatational velocity follow the scaling $k^{-1.88}$ and $k^{-2.20}$ respectively all through the turbulent scale $\sim 2-0.2 h^{-1}Mpc$. The IGM turbulence may contain two different phases, the supersonic phase and the post-supersonic phase, simultaneously along the hierarchical structure formation history. (3) The non-thermal pressure support, measured by the vortical kinetic energy, is comparable with the thermal pressure for $\rho_b \sim 10$, or $T \sim 10^5 K$ at $z=0.25$. The deviation of the baryon fraction from the cosmic mean shows a positive correlation with the turbulence pressure support. (4) A relatively higher numerical viscosity would artificially dissipate both the compressive and vortical motions in the IGM more seriously, resulting in less developed density fluctuation and vorticity, and leading to remarkably shortened turbulence scale. Shocks in regions out of the clusters are outstandingly suppressed by the numerical viscosity since $z=2$, which could directly cause the different levels of turbulence in the IGM between the codes. We also discuss the impact of ingredients not included in our simulation.

Discontinuous Galerkin (DG) methods combine features in finite element methods (weak formulation, finite dimensional solution and test function spaces) and in finite volume methods (numerical fluxes, nonlinear limiters) and are particularly suitable for simulating convection dominated problems, such as linear and nonlinear waves including shock waves. In this article we will give a brief survey of DG methods, emphasizing their applications in computational fluid dynamics (CFD). We will discuss essential ingredients and properties of DG methods, and will also give a few examples of recent developments of DG methods which have facilitated their applications in CFD.

An important property for finite difference schemes designed on curvilinear meshes is the exact preservation of free-stream solutions. This property is difficult to fulfill for high order conservative essentially non-oscillatory (WENO) finite difference schemes. In this paper we explore an alternative flux formulation for such finite difference schemes which can preserve free-stream solutions, based on the numerical technique for the metric terms, which can be applied to this alternative flux formulation but is difficult to be applied to the standard finite difference formulation. Free-stream and vortex preservation properties are investigated, and comparison with standard finite difference WENO schemes is made. Theoretical derivation and numerical results show that the finite difference WENO schemes based on the alternative flux formulation can preserve free-stream and vortex solutions on both stationary and dynamically generalized coordinate systems, hence giving much better performance than the standard finite difference WENO schemes for such problems.

We are interested in front propagation problems in the presence of obstacles. We extend a previous work (Bokanowski, Cheng and Shu 09}), to propose a simple and direct discontinuous Galerkin (DG) method adapted to such front propagation problems. We follow the formulation of Bokanowski et al. (Bokanowski_Forcadel_Zidani_2010), leading to a level set formulation driven by $\min(u_t + H(x,\nabla u), u-g(x))=0$, where $g(x)$ is an obstacle function. The DG scheme is motivated by the variational formulation when the Hamiltonian $H$ is a linear function of $\nabla u$, corresponding to linear convection problems in presence of obstacles. The scheme is then generalized to nonlinear equations, written in an explicit form. Stability analysis are performed for the linear case with Euler forward, a Heun scheme and a Runge-Kutta third order time discretization using the technique proposed in Zhang and Shu 2010. Several numerical examples are provided to demonstrate the robustness of the method. Finally, a narrow band approach is considered in order to reduce the computational cost.