A Positive and Energy Stable Numerical Scheme for the Poisson-Nernst-Planck-Cahn-Hilliard Equations with Steric Interactions

02/22/2020
by   Yiran Qian, et al.
0

We consider numerical methods for the Poisson-Nernst-Planck-Cahn-Hilliard (PNPCH) equations with steric interactions. We propose a novel energy stable numerical scheme that respects mass conservation and positivity at the discrete level. Existence and uniqueness of the solution to the proposed nonlinear scheme are established by showing that the solution is a unique minimizer of a convex functional over a closed, convex domain. The positivity of numerical solutions is further theoretically justified by the singularity of the entropy terms, which prevents the minimizer from approaching zero concentrations. A further numerical analysis proves discrete free-energy dissipation. Extensive numerical tests are performed to validate that the numerical scheme is first-order accurate in time and second-order accurate in space, and is capable of preserving the desired properties, such as mass conservation, positivity, and free energy dissipation, at the discrete level. Moreover, the PNPCH equations and the proposed scheme are applied to study charge dynamics and self-assembled nanopatterns in highly concentrated electrolytes that are widely used in electrochemical energy devices. Numerical results demonstrate that the PNPCH equations and our numerical scheme are able to capture nanostructures, such as lamellar patterns and labyrinthine patterns in electric double layers and the bulk, and multiple time relaxation with multiple time scales. In addition, we numerically characterize the interplay between cross steric interactions of short range and the concentration gradient regularization, and their impact on the development of nanostructures in the equilibrium state.

READ FULL TEXT

page 17

page 19

research
09/07/2023

Second-order, Positive, and Unconditional Energy Dissipative Scheme for Modified Poisson-Nernst-Planck Equations

First-order energy dissipative schemes in time are available in literatu...
research
05/03/2022

A dynamic mass transport method for Poisson-Nernst-Planck equations

A dynamic mass-transport method is proposed for approximately solving th...
research
02/16/2021

A positivity-preserving and convergent numerical scheme for the binary fluid-surfactant system

In this paper, we develop a first order (in time) numerical scheme for t...
research
06/30/2020

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

In this paper, we consider a nonlinear and nonlocal parabolic model for ...
research
02/26/2020

A corrected decoupled scheme for chemotaxis models

The main purpose of this paper is to present a new corrected decoupled s...
research
04/25/2022

Structure-preserving numerical method for Ampere-Nernst-Planck model

Charge dynamics play essential role in many practical applications such ...

Please sign up or login with your details

Forgot password? Click here to reset