In this paper, we develop a class of nonlinear compact schemes based on our previous linear central compact schemes with spectral-like resolution [X. Liu, S. Zhang, H. Zhang and C-W. Shu, A new class of central compact schemes with spectral-like resolution I: Linear schemes, Journal of Computational Physics 248 (2013) 235-256]. In our approach, we compute the flux derivatives on the cell-nodes by the physical fluxes on the cell nodes and numerical fluxes on the cell centers. To acquire the numerical fluxes on the cell centers, we perform a weighted hybrid interpolation of an upwind interpolation and a central interpolation. Through systematic analysis and numerical tests, we show that our nonlinear compact scheme has high order, high resolution and low dissipation, and has the same ability to capture strong discontinuities as regular weighted essentially non-oscillatory (WENO) schemes. It is a good choice for the simulation of multiscale problems with shock waves.

No abstract uploaded!

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 this paper, we develop the high order weighted essentially non-oscillatory (WENO) schemes for solving the Degasperis-Procesi (DP) equation, including finite volume (FV) and finite difference (FD) methods. The DP equation contains nonlinear high order derivatives, and possibly discontinuous or sharp transition solutions. The finite volume method is de- signed based on the total variation bounded property of the DP equation. And the finite difference method is constructed based on the L2 stability of the DP equation. Due to the adoption of the WENO reconstruction, both schemes are arbitrary high order accuracy and shock capturing. The numerical simulation results for different types of solutions of the nonlinear Degasperis-Procesi equation are provided to illustrate the accuracy and capability of the methods.

No abstract uploaded!