Numerical Methods for Partial Differential Equations | 2019
Implicit–explicit schemes for nonlinear nonlocal equations with a gradient flow structure in one space dimension
Abstract
Nonlinear convection-diffusion equations with nonlocal flux and possibly degenerate diffusion arise in various contexts including interacting gases, porous media flows, and collective behavior in biology. Their numerical solution by an explicit finite difference method is costly due to the necessity of discretizing a local spatial convolution for each evaluation of the convective numerical flux, and due to the disadvantageous Courant-Friedrichs-Lewy (CFL) condition incurred by the diffusion term. Based on explicit schemes for such models devised in [J.A. Carrillo, A. Chertock, Y. Huang, Commun. Comput. Phys. vol. 17 (2015) pp. 233– 258] a second-order implicit-explicit Runge-Kutta (IMEX-RK) method can be formulated. This method avoids the restrictive time step limitation of explicit schemes since the diffusion term is handled implicitly, but entails the necessity to solve nonlinear algebraic systems in every time step. It is proven that this method is well defined. Numerical experiments illustrate that for fine discretizations it is more efficient in terms of reduction of error versus CPU time than the original explicit method. One of the test cases is given by a strongly degenerate parabolic, nonlocal equation modelling aggregation [F. Betancourt, R. Bürger, K.H. Karlsen, Commun. Math. Sci. vol. 9 (2011) pp. 711–742]. This model can be transformed to a local partial differential equation that can be solved numerically easily to generate a reference solution for the IMEX-RK method, but is limited to one space dimension.