No abstract uploaded!

In this paper, we design and analyze third order positivity-preserving discontinuous Galerkin (DG) schemes for solving the time-dependent system of Poisson–Nernst–Planck (PNP) equations, which have found much use in diverse applications. Our DG method with Euler forward time discretization is shown to preserve the positivity of cell averages at all time steps. The positivity of numerical solutions is then restored by a scaling limiter in reference to positive weighted cell averages. The method is also shown to preserve steady states. Numerical examples are presented to demonstrate the third order accuracy and illustrate the positivity-preserving property in both one and two dimensions.

In this article, we give a brief overview on high order accurate shock capturing schemes with the aim of applications in compressible turbulence simulations. The emphasis is on the basic methodology and recent algorithm developments for two classes of high order methods: the weighted essentially non-oscillatory (WENO) and discontinuous Galerkin (DG) methods.

We propose efficient Eulerian methods for approximating the finite-time Lyapunov exponent (FTLE). The idea is to compute the related flow map using the Level Set Method and the Liouville equation. There are several advantages of the proposed approach. Unlike the usual Lagrangian-type computations, the resulting method requires the velocity field defined only at discrete locations. No interpolation of the velocity field is needed. Also, the method automatically stops a particle trajectory in the case when the ray hits the boundary of the computational domain. The computational complexity of the algorithm is <i>O</i>(<i>x</i>^(<i>d</i>+1)) with <i>d</i> the dimension of the physical space. Since there are the same number of mesh points in the <i>x</i><i>t</i> space, the computational complexity of the proposed Eulerian approach is optimal in the sense that each grid point is visited for only <i>O</i>(1) time. We also extend the algorithm to compute the FTLE

In this paper, we analyze the stability of the fourth order Runge-Kutta method for integrating semi-discrete approximations of time-dependent partial differential equations. Our study focuses on linear problems and covers general semi-bounded spatial discretizations. A counter example is given to show that the classical four-stage fourth order Runge-Kutta method can not preserve the one-step strong stability, even though the ordinary differential equation system is energy-decaying. But with an energy argument, we show that the strong stability property holds in two steps under an appropriate time step constraint. Based on this fact, the stability extends to general well-posed linear systems. As an application, we utilize the results to examine the stability of the fourth order Runge-Kutta approximations of several specific method of lines schemes for hyperbolic problems, including the spectral Galerkin method and the discontinuous Galerkin method.