A unified structure preserving scheme for a multi-species model with a gradient flow structure and nonlocal interactions via singular kernels

Abstract

In this paper, we consider a nonlinear and nonlocal parabolic model for multi-species ionic fluids and introduce a semi-implicit finite volume scheme, which is second order accurate in space, first order in time and satisfies the following properties: positivity preserving, mass conservation and energy dissipation. Besides, our scheme involves a fast algorithm on the convolution terms with singular but integrable kernels, which otherwise impedes the accuracy and efficiency of the whole scheme. Error estimates on the fast convolution algorithm are shown next. Numerous numerical tests are provided to demonstrate the properties, such as unconditional stability, order of convergence, energy dissipation and the complexity of the fast convolution algorithm. Furthermore, extensive numerical experiments are carried out to explore the modeling effects in specific examples, such as, the steric repulsion, the concentration of ions at the boundary and the blowup phenomenon of the Keller-Segel equations.

arXiv comment27 pages, 17 figures

References in corpus (6)

• Diffuse-Charge Dynamics in Electrochemical Systems
• A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure
• A free energy satisfying finite difference method for Poisson--Nernst--Planck equations
• A free energy satisfying discontinuous Galerkin method for one-dimensional Poisson--Nernst--Planck systems
• Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system
• Field model for complex ionic fluids: analytical properties and numerical investigation