We propose second order accurate discontinuous Galerkin (DG) schemes which satisfy a strict maximum principle for general nonlinear convection-diffusion equations on unstructured triangular meshes. Motivated by genuinely high order maximum-principle-satisfying DG schemes for hyperbolic conservation laws, we prove that under suitable time step restriction for forward Euler time stepping, for general nonlinear convection-diffusion equations, the same scaling limiter coupled with second order DG methods preserves the physical bounds indicated by the initial condition while maintaining uniform second order accuracy. Similar to the purely convection cases, the limiters are mass conservative and easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. Following the idea in earlier work, we extend the schemes to two-dimensional convection-diffusion equations on triangular meshes. There are no geometric constraints on the mesh such as angle acuteness. Numerical results including incompressible Navier-Stokes equations are presented to validate and demonstrate the effectiveness of the numerical methods.

In this paper, we apply the local discontinuous Galerkin (LDG) method to 2D Keller--Segel (KS) chemotaxis model. We improve the results upon (Y. Epshteyn and A. Kurganov, SIAM Journal on Numerical Analysis, 47 (2008), 368-408) and give optimal rate of convergence under special finite element spaces before the blow-up occurs (the exact solutions are smooth). Moreover, to construct physically relevant numerical approximations, we consider $P^1$ LDG scheme and develop a positivity-preserving limiter to the scheme, extending the idea in (Y. Zhang, X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 8918-8934). With this limiter, we can prove the $L^1$-stability of the numerical scheme. Numerical experiments are performed to demonstrate the good performance of the positivity-preserving LDG scheme. Moreover, it is known that the chemotaxis model will yield blow-up solutions under certain initial conditions. We numerically demonstrate how to find the approximate blow-up time by using the $L^2$-norm of the $L^1$-stable numerical solution.

In applications such as astrophysics and inertial confinement fusion, there are many three-dimensional cylindrical-symmetric multi-material problems which are usually simulated by Lagrangian schemes in the two-dimensional cylindrical coordinates. For this type of simulation, a critical issue for the schemes is to keep spherical symmetry in the cylindrical coordinate system if the original physical problem has this symmetry. In the past decades, several Lagrangian schemes with such symmetry property have been developed, but all of them are only first order accurate. In this paper, we develop a second order cell-centered Lagrangian scheme for solving compressible Euler equations in cylindrical coordinates, based on the control volume discretizations, which is designed to have uniformly second order accuracy and capability to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. The scheme maintains several good properties such as conservation for mass, momentum and total energy, and the geometric conservation law. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of accuracy, symmetry, non-oscillation and robustness. The advantage of higher order accuracy is demonstrated in these examples.

The sequential regularization method is a reformulation of the unsteady Navier–Stokes equations from the viewpoint of constrained dynamical systems or the approximate Helmholtz–Hodge projection. In this paper we study the long time behavior of the sequential regularization formulation. We give a uniform-in-time estimate between the solution of the reformulated system and that of the Navier–Stokes equations. We also conduct an error analysis for the temporal discrete system and show that the error bound is independent of time. A couple of long time flow examples are computed to demonstrate this method.

The convergence to steady state solutions of the Euler equations for weighted compact nonlinear schemes (WCNS) [Deng X. and Zhang H. (2000), J. Comput. Phys. 165, 22-44 and Zhang S., Jiang S and Shu C.-W. (2008), J. Comput. Phys. 227, 7294-7321] is studied through numerical tests. Like most other shock capturing schemes, WCNS also suffers from the problem that the residue can not settle down to machine zero for the computation of the steady state solution which contains shock waves but hangs at the truncation error level. In this paper, the techniques studied in [Zhang S. and Shu. C.-W. (2007), J. Sci. Comput. 31, 273-305 and Zhang S., Jiang S and Shu. C.-W. (2011), J. Sci. Comput. 47, 216-238], to improve the convergence to steady state solutions for WENO schemes, are generalized to the WCNS. Detailed numerical studies in one and two dimensional cases are performed. Numerical tests demonstrate the effectiveness of these techniques when applied to WCNS. The residue of various order WCNS can settle down to machine zero for typical cases while the small post-shock oscillations can be removed.